Ancient rapid radiations of insects: challenges for phylogenetic analysis

Phylotranscriptomics: saturated third codon positions radically influence the estimation of trees based on next-gen data

Recent advancements in molecular sequencing techniques have led to a surge in the number of phylogenetic studies that incorporate large amounts of genetic data. We test the assumption that analyzing large number of genes will lead to improvements in tree resolution and branch support using moths in the superfamily Bombycoidea, a group with some interfamilial relationships that have been difficult to resolve. Specifically, we use a next-gen data set that included 19 taxa and 938 genes (∼1.2M bp) to examine how codon position and saturation might influence resolution and node support among three key families. Maximum likelihood, parsimony, and species tree analysis using gene tree parsimony, on different nucleotide and amino acid data sets, resulted in largely congruent topologies with high bootstrap support compared with prior studies that included fewer loci. However, for a few shallow nodes, nucleotide and amino acid data provided high support for conflicting relationships. The third codon position was saturated and phylogenetic analysis of this position alone supported a completely different, potentially misleading sister group relationship. We used the program RADICAL to assess the number of genes needed to fix some of these difficult nodes. One such node originally needed a total of 850 genes but only required 250 when synonymous signal was removed. Our study shows that, in order to effectively use next-gen data to correctly resolve difficult phylogenetic relationships, it is necessary to assess the effects of synonymous substitutions and third codon positions.

Molecular traces of alternative social organization in a termite genome

Although eusociality evolved independently within several orders of insects, research into the molecular underpinnings of the transition towards social complexity has been confined primarily to Hymenoptera (for example, ants and bees). Here we sequence the genome and stage-specific transcriptomes of the dampwood termite Zootermopsis nevadensis (Blattodea) and compare them with similar data for eusocial Hymenoptera, to better identify commonalities and differences in achieving this significant transition. We show an expansion of genes related to male fertility, with upregulated gene expression in male reproductive individuals reflecting the profound differences in mating biology relative to the Hymenoptera. For several chemoreceptor families, we show divergent numbers of genes, which may correspond to the more claustral lifestyle of these termites. We also show similarities in the number and expression of genes related to caste determination mechanisms. Finally, patterns of DNA methylation and alternative splicing support a hypothesized epigenetic regulation of caste differentiation.

Phylogeny and size differentially influence dissolved Cd and Zn bioaccumulation parameters among closely related aquatic insects

Evolutionarily distinct lineages can vary markedly in their accumulation of, and sensitivity to, contaminants. However, less is known about variability among closely related species. Here, we compared dissolved Cd and Zn bioaccumulation in 19 species spanning two species-rich aquatic insect families: Ephemerellidae (order Ephemeroptera (mayflies)), generalized to be metal sensitive, and Hydropsychidae (order Trichoptera (caddisflies)), generalized to be metal tolerant. Across all species, Zn and Cd uptake rate constants (k(u)s), efflux rate constants (k(e)s) and bioconcentration factors (BCFs) strongly covaried, suggesting that these metals share transport pathways in these distinct lineages. K(u)s and BCFs were substantially larger in Ephemerellidae than in Hydropsychidae, whereas k(e)s did not dramatically differ between the two families. Body size played an important role in driving ku differences among species, but had no influence on k(e)s. While familial differences in metal bioconcentration were striking, each family exhibited tremendous variability in all bioaccumulation parameters. At finer levels of taxonomic resolution (within families), phylogeny did not account for differences in metal bioaccumulation. These findings suggest that intrafamily variability can be profound and have important practical implications in that we need to better understand how well "surrogate species" represent their fellow congeners and family members.

Historical biogeography of Eastern Asian-Eastern North American disjunct Melaphidina aphids (Hemiptera: Aphididae: Eriosomatinae) on Rhus hosts (Anacardiaceae)

Intercontinental biotic disjunctions have been documented and analyzed in numerous Holarctic taxa. Patterns previously synthesized for animals compared to plants suggest that the timing of animal disjunctions are mostly Early Tertiary and were generated by migration and vicariance events occurring in the North Atlantic, while plant disjunctions are mostly Mid-Late Tertiary and imply migration and vicariance over Beringia. Melaphidina aphids (Hemiptera: Aphididae: Fordini) exhibit host-alternating life cycles comprising an obligate seasonal shift between Rhus subgenus Rhus species (Anacardiaceae) and mosses (Bryophyta). Similar to their Rhus hosts, melaphidines are distributed disjunctly between Eastern Asia and Eastern North America. We examined evolutionary relationships within Melaphidina to determine the position of the North American lineage, date its divergence from Asian relatives, and compare these results to a previous historical biogeographic study of Rhus. We sampled nine species and three subspecies representing all six genera of Melaphidina. Data included sequences of mitochondrial cytochrome c oxidase subunits I and II+leucine tRNA, cytochrome b, and nuclear elongation factor 1α genes. Phylogenetic analyses (Bayesian, maximum-likelihood, parsimony) of the combined data (3282 bp) supported the monophyly of all genera except Nurudea and Schlechtendalia, due to the position of N. ibofushi. While the exact position of the North American Melaphis was not well resolved, there was high support for a derived position within Asian taxa. The divergence of Melaphis from Asian relatives centered on the Eocene-Oligocene boundary (~33-35Ma), which coincides with closure of Beringian Land Bridge I. This also corresponded to the Asian-North American disjunction previously estimated for subgenus Rhus spp. We suggest the late-Eocene Bering Land Bridge as the most likely migration route for Melaphis ancestors, as was also hypothesized for North American Rhus ancestors. Results for the Melaphidina disjunction depart from the modal pattern in animal lineages, and present a case where insect and host-plant taxa apparently responded similarly to Tertiary climate change.

Spatiotemporal reconstruction of the Aquilegia rapid radiation through next-generation sequencing of rapidly evolving cpDNA regions

Aquilegia is a well-known model system in the field of evolutionary biology, but obtaining a resolved and well-supported phylogenetic reconstruction for the genus has been hindered by its recent and rapid diversification. Here, we applied 454 next-generation sequencing to PCR amplicons of 21 of the most rapidly evolving regions of the plastome to generate c. 24 kb of sequences from each of 84 individuals from throughout the genus. The resulting phylogeny has well-supported resolution of the main lineages of the genus, although recent diversification such as in the European taxa remains unresolved. By producing a chronogram of the whole Ranunculaceae family based on published data, we inferred calibration points for dating the Aquilegia radiation. The genus originated in the upper Miocene c. 6.9 million yr ago (Ma) in Eastern Asia, and diversification occurred c. 4.8 Ma with the split of two main clades, one colonizing North America, and the other Western Eurasia through the mountains of Central Asia. This was followed by a back-to-Asia migration, originating from the European stock using a North Asian route. These results provide the first backbone phylogeny and spatiotemporal reconstruction of the Aquilegia radiation, and constitute a robust framework to address the adaptative nature of speciation within the group.

Characterization of an ancient lepidopteran lateral gene transfer

Bacteria to eukaryote lateral gene transfers (LGT) are an important potential source of material for the evolution of novel genetic traits. The explosion in the number of newly sequenced genomes provides opportunities to identify and characterize examples of these lateral gene transfer events, and to assess their role in the evolution of new genes. In this paper, we describe an ancient lepidopteran LGT of a glycosyl hydrolase family 31 gene (GH31) from an Enterococcus bacteria. PCR amplification between the LGT and a flanking insect gene confirmed that the GH31 was integrated into the Bombyx mori genome and was not a result of an assembly error. Database searches in combination with degenerate PCR on a panel of 7 lepidopteran families confirmed that the GH31 LGT event occurred deep within the Order approximately 65–145 million years ago. The most basal species in which the LGT was found is Plutella xylostella (superfamily: Yponomeutoidea). Array data from Bombyx mori shows that GH31 is expressed, and low dN/dS ratios indicates the LGT coding sequence is under strong stabilizing selection. These findings provide further support for the proposition that bacterial LGTs are relatively common in insects and likely to be an underappreciated source of adaptive genetic material.

A unique box in 28S rRNA is shared by the enigmatic insect order Zoraptera and Dictyoptera

The position of the Zoraptera remains one of the most challenging and uncertain concerns in ordinal-level phylogenies of the insects. Zoraptera have been viewed as having a close relationship with five different groups of Polyneoptera, or as being allied to the Paraneoptera or even Holometabola. Although rDNAs have been widely used in phylogenetic studies of insects, the application of the complete 28S rDNA are still scattered in only a few orders. In this study, a secondary structure model of the complete 28S rRNAs of insects was reconstructed based on all orders of Insecta. It was found that one length-variable region, D3-4, is particularly distinctive. The length and/or sequence of D3-4 is conservative within each order of Polyneoptera, but it can be divided into two types between the different orders of the supercohort, of which the enigmatic order Zoraptera and Dictyoptera share one type, while the remaining orders of Polyneoptera share the other. Additionally, independent evidence from phylogenetic results support the clade (Zoraptera+Dictyoptera) as well. Thus, the similarity of D3-4 between Zoraptera and Dictyoptera can serve as potentially valuable autapomorphy or synapomorphy in phylogeny reconstruction. The clades of (Plecoptera+Dermaptera) and ((Grylloblattodea+Mantophasmatodea)+(Embiodea+Phasmatodea)) were also recovered in the phylogenetic study. In addition, considering the other studies based on rDNAs, this study reached the highest congruence with previous phylogenetic studies of Holometabola based on nuclear protein coding genes or morphology characters. Future comparative studies of secondary structures across deep divergences and additional taxa are likely to reveal conserved patterns, structures and motifs that can provide support for major phylogenetic lineages.

Tomographic reconstruction of neopterous carboniferous insect nymphs

Two new polyneopteran insect nymphs from the Montceau-les-Mines Lagerstätte of France are presented. Both are preserved in three dimensions, and are imaged with the aid of X-ray micro-tomography, allowing their morphology to be recovered in unprecedented detail. One-Anebos phrixos gen. et sp. nov.-is of uncertain affinities, and preserves portions of the antennae and eyes, coupled with a heavily spined habitus. The other is a roachoid with long antennae and chewing mouthparts very similar in form to the most generalized mandibulate mouthparts of extant orthopteroid insects. Computer reconstructions reveal limbs in both specimens, allowing identification of the segments and annulation in the tarsus, while poorly developed thoracic wing pads suggest both are young instars. This work describes the morphologically best-known Palaeozoic insect nymphs, allowing a better understanding of the juveniles' palaeobiology and palaeoecology. We also consider the validity of evidence from Palaeozoic juvenile insects in wing origin theories. The study of juvenile Palaeozoic insects is currently a neglected field, yet these fossils provide direct evidence on the evolution of insect development. It is hoped this study will stimulate a renewed interest in such work.

A remarkable new family of Jurassic insects (Neuroptera) with primitive wing venation and its phylogenetic position in Neuropterida

BACKGROUND

Lacewings (insect order Neuroptera), known in the fossil record since the Early Permian, were most diverse in the Mesozoic. A dramatic variety of forms ranged in that time from large butterfly-like Kalligrammatidae to minute two-winged Dipteromantispidae.

PRINCIPAL FINDINGS

We describe the intriguing new neuropteran family Parakseneuridae fam. nov. with three new genera and 15 new species from the Middle Jurassic of Daohugou (Inner Mongolia, China) and the Early/Middle Jurassic of Sai-Sagul (Kyrgyzstan): Parakseneura undula gen. et sp. nov., P. albomacula gen. et sp. nov., P. curvivenis gen. et sp. nov., P. nigromacula gen. et sp. nov., P. nigrolinea gen. et sp. nov., P. albadelta gen. et sp. nov., P. cavomaculata gen. et sp. nov., P. inflata gen. et sp. nov., P. metallica gen. et sp. nov., P. emarginata gen. et sp. nov., P. directa gen. et sp. nov., Pseudorapisma jurassicum gen. et sp. nov., P. angustipenne gen. et sp. nov., P. maculatum gen. et sp. nov. (Daohugou); Shuraboneura ovata gen. et sp. nov. (Sai-Sagul). The family comprises large neuropterans with most primitive wing venation in the order indicated by the presence of ScA and AA1+2, and the dichotomous branching of MP, CuA, CuP, AA3+4, AP1+2. The phylogenetic position of Parakseneuridae was investigated using a phylogenetic analysis of morphological scoring for 33 families of extinct and extant Neuropterida combined with DNA sequence data for representatives of all extant families. Parakseneuridae were recovered in a clade with Osmylopsychopidae, Prohemerobiidae, and Ithonidae.

CONCLUSIONS/SIGNIFICANCE

The presence of the presumed AA1+2 in wings of Parakseneuridae is a unique plesiomorphic condition hitherto unknown in Neuropterida, the clade comprising Neuroptera, Megaloptera, Raphidioptera. The relative uncertainty of phylogenetic position of Parakseneuridae and the majority of other families of Neuroptera reflects deficient paleontological data, especially from critical important periods for the order, earliest Triassic and latest Triassic/earliest Jurassic.

Insect phylogenomics: results, problems and the impact of matrix composition

In this study, we investigated the relationships among insect orders with a main focus on Polyneoptera (lower Neoptera: roaches, mantids, earwigs, grasshoppers, etc.), and Paraneoptera (thrips, lice, bugs in the wide sense). The relationships between and within these groups of insects are difficult to resolve because only few informative molecular and morphological characters are available. Here, we provide the first phylogenomic expressed sequence tags data ('EST': short sub-sequences from a c(opy) DNA sequence encoding for proteins) for stick insects (Phasmatodea) and webspinners (Embioptera) to complete published EST data. As recent EST datasets are characterized by a heterogeneous distribution of available genes across taxa, we use different rationales to optimize the data matrix composition. Our results suggest a monophyletic origin of Polyneoptera and Eumetabola (Paraneoptera + Holometabola). However, we identified artefacts of tree reconstruction (human louse Pediculus humanus assigned to Odonata (damselflies and dragonflies) or Holometabola (insects with a complete metamorphosis); mayfly genus Baetis nested within Neoptera), which were most probably rooted in a data matrix composition bias due to the inclusion of sequence data of entire proteomes. Until entire proteomes are available for each species in phylogenomic analyses, this potential pitfall should be carefully considered.

The Distribution of eIF4E-Family Members across Insecta

Insects are part of the earliest faunas that invaded terrestrial environments and are the first organisms that evolved controlled flight. Nowadays, insects are the most diverse animal group on the planet and comprise the majority of extant animal species described. Moreover, they have a huge impact in the biosphere as well as in all aspects of human life and economy; therefore understanding all aspects of insect biology is of great importance. In insects, as in all cells, translation is a fundamental process for gene expression. However, translation in insects has been mostly studied only in the model organism Drosophila melanogaster. We used all publicly available genomic sequences to investigate in insects the distribution of the genes encoding the cap-binding protein eIF4E, a protein that plays a crucial role in eukaryotic translation. We found that there is a diversity of multiple ortholog genes encoding eIF4E isoforms within the genus Drosophila. In striking contrast, insects outside this genus contain only a single eIF4E gene, related to D. melanogaster eIF4E-1. We also found that all insect species here analyzed contain only one Class II gene, termed 4E-HP. We discuss the possible evolutionary causes originating the multiplicity of eIF4E genes within the genus Drosophila.

Cretaceous origin and repeated tertiary diversification of the redefined butterflies

Although the taxonomy of the ca 18 000 species of butterflies and skippers is well known, the family-level relationships are still debated. Here, we present, to our knowledge, the most comprehensive phylogenetic analysis of the superfamilies Papilionoidea, Hesperioidea and Hedyloidea to date based on morphological and molecular data. We reconstructed their phylogenetic relationships using parsimony and Bayesian approaches. We estimated times and rates of diversification along lineages in order to reconstruct their evolutionary history. Our results suggest that the butterflies, as traditionally understood, are paraphyletic, with Papilionidae being the sister-group to Hesperioidea, Hedyloidea and all other butterflies. Hence, the families in the current three superfamilies should be placed in a single superfamily Papilionoidea. In addition, we find that Hedylidae is sister to Hesperiidae, and this novel relationship is supported by two morphological characters. The families diverged in the Early Cretaceous but diversified after the Cretaceous-Palaeogene event. The diversification of butterflies is characterized by a slow speciation rate in the lineage leading to Baronia brevicornis, a period of stasis by the skippers after divergence and a burst of diversification in the lineages leading to Nymphalidae, Riodinidae and Lycaenidae.

Evolution of a large, conserved, and syntenic gene family in insects

The Osiris gene family, first described in Drosophila melanogaster, is clustered in the genomes of all Drosophila species sequenced to date. In D. melanogaster, it explains the enigmatic phenomenon of the triplo-lethal and haploinsufficient locus Tpl. The synteny of Osiris genes in flies is well conserved, and it is one of the largest syntenic blocks in the Drosophila group. By examining the genome sequences of other insects in a wide range of taxonomic orders, we show here that the gene family is well-conserved and syntenic not only in the diptera but across the holometabolous and hemimetabolous insects. Osiris gene homologs have also been found in the expressed sequence tag sequences of various other insects but are absent from all groups that are not insects, including crustacea and arachnids. It is clear that the gene family evolved by gene duplication and neofunctionalization very soon after the divergence of the insects from other arthropods but before the divergence of the insects from one another and that the sequences and synteny have been maintained by selection ever since.

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).

Phylogeny of the mega-diverse Gelechioidea (Lepidoptera): adaptations and determinants of success

The Gelechioidea, with 18,000 described and many more unnamed species ranks among the most diverse lepidopteran superfamilies. Nevertheless, their taxonomy has remained largely unresolved, and phylogenetic affinities among gelechioid families and lower taxa have been insufficiently understood. We constructed, for the first time, a comprehensive molecular phylogeny for the Gelechioidea. We sampled seven genes, in total 5466 base pairs, of 109 gelechioid taxa representing 32 of 37 recognized subfamilies, and two outgroup taxa. We used maximum likelihood methods and Bayesian inference to construct phylogenetic trees. We found that the families Autostichidae, Lecithoceridae, Xyloryctidae, and Oecophoridae s. str., in this order, are the most basally arising clades. Elachistidae s. l. was found to be paraphyletic, with families such as Gelechiidae and Cosmopterigidae nested within it, and Parametriotinae associated with several families previously considered unrelated to them. Using the phylogenetic trees, we examined patterns of life history evolution and determinants of the success of different lineages. Gelechioids express unusually wide variability in life-history strategies, including herbivorous, saprophagous, fungivorous, and carnivorous lineages. Most species are highly specialized in diet and other life history traits. The results suggest that either saprophagy was the ancestral feeding strategy from which herbivory evolved independently on multiple occasions, or that the ancestor was herbivorous with repeated origins of saprophagy. External feeding is an ancestral trait from which internal feeding evolved independently several times. In terms of species number, saprophages are dominant in Australia, while elsewhere several phytophagous lineages have extensively specialized and diversified. Internal feeding has remained a somewhat less generally adopted feeding mode, although in a few lineages significant radiations of leaf mining species have occurred. We conclude that diverse feeding modes, specialization among saprophages, repeated shifts to phytophagy, and a generally high specialization rate on single plant species (monophagy) are the major factors behind the success of the Gelechioidea.

Increased gene sampling strengthens support for higher-level groups within leaf-mining moths and relatives (Lepidoptera: Gracillariidae)

BACKGROUND

Researchers conducting molecular phylogenetic studies are frequently faced with the decision of what to do when weak branch support is obtained for key nodes of importance. As one solution, the researcher may choose to sequence additional orthologous genes of appropriate evolutionary rate for the taxa in the study. However, generating large, complete data matrices can become increasingly difficult as the number of characters increases. A few empirical studies have shown that augmenting genes even for a subset of taxa can improve branch support. However, because each study differs in the number of characters and taxa, there is still a need for additional studies that examine whether incomplete sampling designs are likely to aid at increasing deep node resolution. We target Gracillariidae, a Cretaceous-age (~100 Ma) group of leaf-mining moths to test whether the strategy of adding genes for a subset of taxa can improve branch support for deep nodes. We initially sequenced ten genes (8,418 bp) for 57 taxa that represent the major lineages of Gracillariidae plus outgroups. After finding that many deep divergences remained weakly supported, we sequenced eleven additional genes (6,375 bp) for a 27-taxon subset. We then compared results from different data sets to assess whether one sampling design can be favored over another. The concatenated data set comprising all genes and all taxa and three other data sets of different taxon and gene sub-sampling design were analyzed with maximum likelihood. Each data set was subject to five different models and partitioning schemes of non-synonymous and synonymous changes. Statistical significance of non-monophyly was examined with the Approximately Unbiased (AU) test.

RESULTS

Partial augmentation of genes led to high support for deep divergences, especially when non-synonymous changes were analyzed alone. Increasing the number of taxa without an increase in number of characters led to lower bootstrap support; increasing the number of characters without increasing the number of taxa generally increased bootstrap support. More than three-quarters of nodes were supported with bootstrap values greater than 80% when all taxa and genes were combined. Gracillariidae, Lithocolletinae + Leucanthiza, and Acrocercops and Parectopa groups were strongly supported in nearly every analysis. Gracillaria group was well supported in some analyses, but less so in others. We find strong evidence for the exclusion of Douglasiidae from Gracillarioidea sensu Davis and Robinson (1998). Our results strongly support the monophyly of a G.B.R.Y. clade, a group comprised of Gracillariidae + Bucculatricidae + Roeslerstammiidae + Yponomeutidae, when analyzed with non-synonymous changes only, but this group was frequently split when synonymous and non-synonymous substitutions were analyzed together.

CONCLUSIONS

1) Partially or fully augmenting a data set with more characters increased bootstrap support for particular deep nodes, and this increase was dramatic when non-synonymous changes were analyzed alone. Thus, the addition of sites that have low levels of saturation and compositional heterogeneity can greatly improve results. 2) Gracillarioidea, as defined by Davis and Robinson (1998), clearly do not include Douglasiidae, and changes to current classification will be required. 3) Gracillariidae were monophyletic in all analyses conducted, and nearly all species can be placed into one of six strongly supported clades though relationships among these remain unclear. 4) The difficulty in determining the phylogenetic placement of Bucculatricidae is probably attributable to compositional heterogeneity at the third codon position. From our tests for compositional heterogeneity and strong bootstrap values obtained when synonymous changes are excluded, we tentatively conclude that Bucculatricidae is closely related to Gracillariidae + Roeslerstammiidae + Yponomeutidae.

Overcoming the effects of rogue taxa: Evolutionary relationships of the bee flies

Bombyliidae (5000 sp.), or bee flies, are a lower brachyceran family of flower-visiting flies that, as larvae, act as parasitoids of other insects. The evolutionary relationships are known from a morphological analysis that yielded minimal support for higher-level groupings. We use the protein-coding gene CAD and 28S rDNA to determine phylogeny and to test the monophyly of existing subfamilies, the divisions Tomophtalmae, and 'the sand chamber subfamilies'. Additionally, we demonstrate that consensus networks can be used to identify rogue taxa in a Bayesian framework. Pruning rogue taxa post-analysis from the final tree distribution results in increased posterior probabilities. We find 8 subfamilies to be monophyletic and the subfamilies Heterotropinae and Mythicomyiinae to be the earliest diverging lineages. The large subfamily Bombyliinae is found to be polyphyletic and our data does not provide evidence for the monophyly of Tomophthalmae or the 'sand chamber subfamilies'.

A second generation genetic map of the bumblebee Bombus terrestris (Linnaeus, 1758) reveals slow genome and chromosome evolution in the Apidae

Background

The bumblebee Bombus terrestris is an ecologically and economically important pollinator and has become an important biological model system. To study fundamental evolutionary questions at the genomic level, a high resolution genetic linkage map is an essential tool for analyses ranging from quantitative trait loci (QTL) mapping to genome assembly and comparative genomics. We here present a saturated linkage map and match it with the Apis mellifera genome using homologous markers. This genome-wide comparison allows insights into structural conservations and rearrangements and thus the evolution on a chromosomal level.

Results

The high density linkage map covers ~ 93% of the B. terrestris genome on 18 linkage groups (LGs) and has a length of 2'047 cM with an average marker distance of 4.02 cM. Based on a genome size of ~ 430 Mb, the recombination rate estimate is 4.76 cM/Mb. Sequence homologies of 242 homologous markers allowed to match 15 B. terrestris with A. mellifera LGs, five of them as composites. Comparing marker orders between both genomes we detect over 14% of the genome to be organized in synteny and 21% in rearranged blocks on the same homologous LG.

Conclusions

This study demonstrates that, despite the very high recombination rates of both A. mellifera and B. terrestris and a long divergence time of about 100 million years, the genomes' genetic architecture is highly conserved. This reflects a slow genome evolution in these bees. We show that data on genome organization and conserved molecular markers can be used as a powerful tool for comparative genomics and evolutionary studies, opening up new avenues of research in the Apidae.

The mitochondrial genome of Bacillus stick insects (Phasmatodea) and the phylogeny of orthopteroid insects

The Order Phasmatodea (stick and leaf insects) includes many well-known species of cryptic phytophagous insects. In this work, we sequenced the almost complete mitochondrial genomes of two stick insect species of the genus Bacillus. Phasmatodea pertain to the Polyneoptera, and represent one of the major clades of heterometabolous insects. Orthopteroid insect lineages arose through rapid evolutionary radiation events, which likely blurred the phylogenetic reconstructions obtained so far; we therefore performed a phylogenetic analysis to resolve and date all major splits of orthopteroid phylogeny, including the relationships between Phasmatodea and other polyneopterans. We explored several molecular models, with special reference to data partitioning, to correctly detect any phylogenetic signal lying in rough data. Phylogenetic Informativeness analysis showed that the maximum resolving power on the orthopteroid mtDNA dataset is expected for the Upper Cretaceous, about 80millionyears ago (Mya), but at least 70% of the maximum informativeness is also expected for the 150-200 Mya timespan, which makes mtDNA a suitable marker to study orthopteroid splits. A complete chronological calibration has also been computed following a Penalized Likelihood method. In summary, our analysis confirmed the monophyly of Phasmatodea, Dictyoptera and Orthoptera, and retrieved Mantophasmatodea as sister group of Phasmatodea. The origin of orthopteroid insects was also estimated to be in the Middle Triassic, while the order Phasmatodea seems to appear in the Upper Jurassic. The obtained results evidenced that mtDNA is a suitable marker to unravel the ancient splits leading to the orthopteroid orders, given a proper methodological approach.

Phylogenetic relationships among insect orders based on three nuclear protein-coding gene sequences

Many attempts to resolve the phylogenetic relationships of higher groups of insects have been made based on both morphological and molecular evidence; nonetheless, most of the interordinal relationships of insects remain unclear or are controversial. As a new approach, in this study we sequenced three nuclear genes encoding the catalytic subunit of DNA polymerase delta and the two largest subunits of RNA polymerase II from all insect orders. The predicted amino acid sequences (In total, approx. 3500 amino acid sites) of these proteins were subjected to phylogenetic analyses based on the maximum likelihood and Bayesian analysis methods with various models. The resulting trees strongly support the monophyly of Palaeoptera, Neoptera, Polyneoptera, and Holometabola, while within Polyneoptera, the groupings of Isoptera/"Blattaria"/Mantodea (Superorder Dictyoptera), Dictyoptera/Zoraptera, Dermaptera/Plecoptera, Mantophasmatodea/Grylloblattodea, and Embioptera/Phasmatodea are supported. Although Paraneoptera is not supported as a monophyletic group, the grouping of Phthiraptera/Psocoptera is robustly supported. The interordinal relationships within Holometabola are well resolved and strongly supported that the order Hymenoptera is the sister lineage to all other holometabolous insects. The other orders of Holometabola are separated into two large groups, and the interordinal relationships of each group are (((Siphonaptera, Mecoptera), Diptera), (Trichoptera, Lepidoptera)) and ((Coleoptera, Strepsiptera), (Neuroptera, Raphidioptera, Megaloptera)). The sister relationship between Strepsiptera and Diptera are significantly rejected by all the statistical tests (AU, KH and wSH), while the affinity between Hymenoptera and Mecopterida are significantly rejected only by AU and KH tests. Our results show that the use of amino acid sequences of these three nuclear genes is an effective approach for resolving the relationships of higher groups of insects.

Mitochondrial genome evolution in fire ants (Hymenoptera: Formicidae)

Background

Complete mitochondrial genome sequences have become important tools for the study of genome architecture, phylogeny, and molecular evolution. Despite the rapid increase in available mitogenomes, the taxonomic sampling often poorly reflects phylogenetic diversity and is often also biased to represent deeper (family-level) evolutionary relationships.

Results

We present the first fully sequenced ant (Hymenoptera: Formicidae) mitochondrial genomes. We sampled four mitogenomes from three species of fire ants, genus Solenopsis, which represent various evolutionary depths. Overall, ant mitogenomes appear to be typical of hymenopteran mitogenomes, displaying a general A+T-bias. The Solenopsis mitogenomes are slightly more compact than other hymentoperan mitogenomes (~15.5 kb), retaining all protein coding genes, ribosomal, and transfer RNAs. We also present evidence of recombination between the mitogenomes of the two conspecific Solenopsis mitogenomes. Finally, we discuss potential ways to improve the estimation of phylogenies using complete mitochondrial genome sequences.

Conclusions

The ant mitogenome presents an important addition to the continued efforts in studying hymenopteran mitogenome architecture, evolution, and phylogenetics. We provide further evidence that the sampling across many taxonomic levels (including conspecifics and congeners) is useful and important to gain detailed insights into mitogenome evolution. We also discuss ways that may help improve the use of mitogenomes in phylogenetic analyses by accounting for non-stationary and non-homogeneous evolution among branches.

Comprehensive gene and taxon coverage elucidates radiation patterns in moths and butterflies

Lepidoptera (butterflies and moths) represent one of the most diverse animals groups. Yet, the phylogeny of advanced ditrysian Lepidoptera, accounting for about 99 per cent of lepidopteran species, has remained largely unresolved. We report a rigorous and comprehensive analysis of lepidopteran affinities. We performed phylogenetic analyses of 350 taxa representing nearly 90 per cent of lepidopteran families. We found Ditrysia to be a monophyletic taxon with the clade Tischerioidea + Palaephatoidea being the sister group of it. No support for the monophyly of the proposed major internested ditrysian clades, Apoditrysia, Obtectomera and Macrolepidoptera, was found as currently defined, but each of these is supported with some modification. The monophyly or near-monophyly of most previously identified lepidopteran superfamilies is reinforced, but several species-rich superfamilies were found to be para- or polyphyletic. Butterflies were found to be more closely related to 'microlepidopteran' groups of moths rather than the clade Macrolepidoptera, where they have traditionally been placed. There is support for the monophyly of Macrolepidoptera when butterflies and Calliduloidea are excluded. The data suggest that the generally short diverging nodes between major groupings in basal non-tineoid Ditrysia are owing to their rapid radiation, presumably in correlation with the radiation of flowering plants.

Phylogenetic diversification patterns and divergence times in ground beetles (Coleoptera: Carabidae: Harpalinae)

BACKGROUND

Harpalinae is a species rich clade of carabid beetles with many unusual morphological forms and ecological interactions. How this diversity evolved has been difficult to reconstruct, perhaps because harpalines underwent a rapid burst of diversification early in their evolutionary history. Here we investigate the tempo of evolution in harpalines using molecular divergence dating techniques and explore the rates of lineage accumulation in harpalines and their sister group.

RESULTS

According to molecular divergence date estimates, harpalines originated in the mid Cretaceous but did not diversify extensively until the late Cretaceous or early Paleogene about 32 million years after their origin. In a relatively small window of time, harpalines underwent rapid speciation. Harpalines have a relative high net diversification rate and increased cladogenesis in some regions of the clade. We did not see a significant decrease in diversification rate through time in the MCCR test, but a model of diversification with two shift points to lower diversification rates fit the harpaline lineage accumulation through time the best.

CONCLUSIONS

Our results indicate harpalines are significantly more diverse and have higher diversification than their sistergroup. Instead of an immediate burst of explosive diversification, harpalines may have had a long "fuse" before major lineages diversified during the early Paleogene when other taxa such as mammals, birds, and some flowering plants were also rapidly diversifying.

The complete mitochondrial genome and phylogenomics of a damselfly, Euphaea formosa support a basal Odonata within the Pterygota

This study determined the first complete mitochondrial genome of a damselfly, Euphaea formosa (Insecta: Odonata: Zygoptera), and reconstructed a phylogeny based on thirteen protein-coding genes of mitochondrial genomes in twenty-five representative hexapods to examine the relationships among the basal Pterygota. The damselfly's mitochondrial genome is a circular molecule of 15,700bp long, and contains the entire set of thirty-seven genes typically found in insects. The gene arrangement, nucleotide composition, and codon usage pattern of the mitochondrial genome are similar across the three odonate species, suggesting a conserved genome evolution within the Odonata. The presence of the intergenic spacer s5 likely represents a synapomorphy for the dragonflies (Anisoptera). Maximum parsimony, maximum likelihood, and Bayesian analyses of both nucleotide and amino acid sequences cannot support the three existing phylogenetic hypotheses of the basal Pterygota (Palaeoptera, Metapterygota, and Chiastomyaria). In contrast, the phylogenetic results indicate an alternative hypothesis of a strongly supported basal Odonata and a sister relationship of the Ephemeroptera and Plecoptera. The unexpected sister Ephemeroptera+Plecoptera clade, which contradicts with the widely accepted hypothesis of a monophyletic Neoptera, requires further analyses with additional mitochondrial genome sampling at the base of the Neoptera.

Phylogenetics and biogeography of a spectacular Old World radiation of butterflies: the subtribe Mycalesina (Lepidoptera: Nymphalidae: Satyrini)

Background

Butterflies of the subtribe Mycalesina (Nymphalidae: Satyrinae) are important model organisms in ecology and evolution. This group has radiated spectacularly in the Old World tropics and presents an exciting opportunity to better understand processes of invertebrate rapid radiations. However, the generic-level taxonomy of the subtribe has been in a constant state of flux, and relationships among genera are unknown. There are six currently recognized genera in the group. Mycalesis, Lohora and Nirvanopsis are found in the Oriental region, the first of which is the most speciose genus among mycalesines, and extends into the Australasian region. Hallelesis and Bicyclus are found in mainland Africa, while Heteropsis is primarily Madagascan, with a few species in Africa. We infer the phylogeny of the group with data from three genes (total of 3139 bp) and use these data to reconstruct events in the biogeographic history of the group.

Results

The results indicate that the group Mycalesina radiated rapidly around the Oligocene-Miocene boundary. Basal relationships are unresolved, but we recover six well-supported clades. Some species of Mycalesis are nested within a primarily Madagascan clade of Heteropsis, while Nirvanopsis is nested within Lohora. The phylogeny suggests that the group had its origin either in Asia or Africa, and diversified through dispersals between the two regions, during the late Oligocene and early Miocene. The current dataset tentatively suggests that the Madagascan fauna comprises two independent radiations. The Australasian radiation shares a common ancestor derived from Asia. We discuss factors that are likely to have played a key role in the diversification of the group.

Conclusions

We propose a significantly revised classification scheme for Mycalesina. We conclude that the group originated and radiated from an ancestor that was found either in Asia or Africa, with dispersals between the two regions and to Australasia. Our phylogeny paves the way for further comparative studies on this group that will help us understand the processes underlying diversification in rapid radiations of invertebrates.

A phylogenomic approach to resolve the arthropod tree of life

Arthropods were the first animals to conquer land and air. They encompass more than three quarters of all described living species. This extraordinary evolutionary success is based on an astoundingly wide array of highly adaptive body organizations. A lack of robustly resolved phylogenetic relationships, however, currently impedes the reliable reconstruction of the underlying evolutionary processes. Here, we show that phylogenomic data can substantially advance our understanding of arthropod evolution and resolve several conflicts among existing hypotheses. We assembled a data set of 233 taxa and 775 genes from which an optimally informative data set of 117 taxa and 129 genes was finally selected using new heuristics and compared with the unreduced data set. We included novel expressed sequence tag (EST) data for 11 species and all published phylogenomic data augmented by recently published EST data on taxonomically important arthropod taxa. This thorough sampling reduces the chance of obtaining spurious results due to stochastic effects of undersampling taxa and genes. Orthology prediction of genes, alignment masking tools, and selection of most informative genes due to a balanced taxa-gene ratio using new heuristics were established. Our optimized data set robustly resolves major arthropod relationships. We received strong support for a sister group relationship of onychophorans and euarthropods and strong support for a close association of tardigrades and cycloneuralia. Within pancrustaceans, our analyses yielded paraphyletic crustaceans and monophyletic hexapods and robustly resolved monophyletic endopterygote insects. However, our analyses also showed for few deep splits that were recently thought to be resolved, for example, the position of myriapods, a remarkable sensitivity to methods of analyses.

The impact of rRNA secondary structure consideration in alignment and tree reconstruction: simulated data and a case study on the phylogeny of hexapods

The use of secondary structures has been advocated to improve both the alignment and the tree reconstruction processes of ribosomal RNA (rRNA) data sets. We used simulated and empirical rRNA data to test the impact of secondary structure consideration in both steps of molecular phylogenetic analyses. A simulation approach was used to generate realistic rRNA data sets based on real 16S, 18S, and 28S sequences and structures in combination with different branch length and topologies. Alignment and tree reconstruction performance of four recent structural alignment methods was compared with exclusively sequence-based approaches. As empirical data, we used a hexapod rRNA data set to study the influence of nucleotide interdependencies in sequence alignment and tree reconstruction. Structural alignment methods delivered significantly better sequence alignments compared with pure sequence-based methods. Also, structural alignment methods delivered better trees judged by topological congruence to simulation base trees. However, the advantage of structural alignments was less pronounced and even vanished in several instances. For simulated data, application of mixed RNA/DNA models to stems and loops, respectively, led to significantly shorter branches. The application of mixed RNA/DNA models in the hexapod analyses delivered partly implausible relationships. This can be interpreted as a stronger sensitivity of mixed model setups to nonphylogenetic signal. Secondary structure consideration clearly influenced sequence alignment and tree reconstruction of ribosomal genes. Although sequence alignment quality can considerably be improved by the use of secondary structure information, the application of mixed models in tree reconstructions needs further studies to understand the observed effects.

A multigene phylogeny of the fly superfamily Asiloidea (Insecta): Taxon sampling and additional genes reveal the sister-group to all higher flies (Cyclorrhapha)

Asiloidea are a group of 9 lower brachyceran fly families, considered to be the closest relative to the large Metazoan radiation Eremoneura (Cyclorrhapha+Empidoidea). The evidence for asiloid monophyly is limited, and few characters define the relationships between the families of Asiloidea and Eremoneura. Additionally, enigmatic genera, Hilarimorpha and Apystomyia, retain morphological characters of both asiloids and higher flies. We use the nuclear protein-coding gene CAD and 28S rDNA to test the monophyly of Asiloidea and to resolve its relationship to Eremoneura. We explore the effects of taxon sampling on support values and topological stability, the resolving power of additional genes, and hypothesis testing using four-cluster likelihood mapping. We find that: (1) the 'asiloid' genus Apystomyia is sister to Cyclorrhapha, (2) the remaining asiloids are monophyletic at the exclusion of the family Bombyliidae, and (3) our best estimate of relationships places the asiloid flies excluding Bombyliidae as the sister-group to Eremoneura, though high support is lacking.

Ribosomal protein genes of holometabolan insects reject the Halteria, instead revealing a close affinity of Strepsiptera with Coleoptera

The phylogenetic relationships among holometabolan insect orders remain poorly known, despite a wealth of previous studies. In particular, past attempts to clarify the sister-group of the enigmatic order Strepsiptera with rRNA genes have led to intense debate about long-branch attraction (the 'Strepsiptera problem'), without resolving the taxonomic question at hand. Here, we appealed to alternative nuclear sequences of 27 ribosomal proteins (RPs) to generate a data matrix of 10,731 nucleotides for 22 holometabolan taxa, including two strepsipteran species. Phylogenetic relationships among holometabolan insects were analyzed under several nucleotide-coding schemes to explore differences in signal and systematic biases. Saturation and compositional bias particularly affected third positions, which greatly differed in AT content (18-72%). Such confounding factors were best reduced by R-Y coding and removal of third codon positions, resulting in more strongly supported topologies, whereas amino acid coding gave poor resolution. The placement of Strepsiptera with Coleoptera (the Coleopterida) was recovered under most coding schemes and analytical methods, if often with modest support and ambiguity. In contrast, an alternative sister-group with Diptera (the Halteria) was only found in one analysis using parsimony, and weakly supported. The topologies here generally support a Coleoptera+Strepsiptera as sister-group to Mecopterida (Siphonaptera+Mecoptera+Diptera+Lepidoptera+Trichoptera), while Hymenoptera were always recovered as sister-group to the remaining Holometabola.

The pea aphid phylome: a complete catalogue of evolutionary histories and arthropod orthology and paralogy relationships for Acyrthosiphon pisum genes

Phylogenetic analyses serve many purposes, including the establishment of orthology relationships, the prediction of protein function and the detection of important evolutionary events. Within the context of the sequencing of the genome of the pea aphid, Acyrthosiphon pisum, we undertook a phylogenetic analysis for every protein of this species. The resulting phylome includes the evolutionary relationships of all predicted aphid proteins and their homologues among 13 other fully-sequenced arthropods and three out-group species. Subsequent analyses have revealed multiple gene expansions that are specific to aphids and have served to transfer functional annotations to 4058 pea aphid genes that display one-to-one orthology relationships with Drosophila melanogaster annotated genes. All phylogenies and alignments are accessible through the PhylomeDB database. Here we provide a description of this dataset and provide some examples on how can it be exploited.

Molecular phylogeny of the small carpenter bees (Hymenoptera: Apidae: Ceratinini) indicates early and rapid global dispersal

The small carpenter bees (tribe Ceratinini, family Apidae) are recorded from all continents except Antarctica. The Ceratinini have a near-global distribution which contrasts strongly with their sister tribe, the Allodapini which has a largely southern Old World distribution. The Ceratinini therefore provides an excellent group to understand the factors that help determine the biogeography and radiation of the bees. This is the first molecular study of ceratinine bees covering representatives from both northern and southern hemisphere Old and New World regions. We use two mitochondrial and one nuclear marker (totalling 2807 nucleotides) to examine the age, cladogenesis and historical biogeography of this tribe. Tree topology and molecular dating support an African origin at about 47 Mya with subsequent dispersal into Eurasia 44 Mya, and followed by an American invasion 32 Mya. Concentrated African and Malagasy sampling revealed there were two or three dispersals events into Madagascar ranging from 25 to 9 Mya. Lineage through time analyses suggest higher rates of cladogenesis close to the origin of the tribe, and this corresponds to both major dispersal events and divergences of lineages leading to extant subgenera. Ceratinini have potentially great importance for future studies to understand the relative roles of dispersal ability and time of origin in determining bee biogeography.