A radiation of Psylliodes flea beetles on Brassicaceae is associated with the evolution of specific detoxification enzymes.

Flea beetles of the genus Psylliodes have evolved specialized interactions with plant species belonging to several distantly related families, mainly Brassicaceae, Solanaceae, and Fagaceae. This diverse host use indicates that Psylliodes flea beetles are able to cope with different chemical defense metabolites, including glucosinolates, the characteristic defense metabolites of Brassicaceae. Here we investigated the evolution of host use and the emergence of a glucosinolate-specific detoxification mechanism in Psylliodes flea beetles. In phylogenetic analyses, Psylliodes species clustered into four major clades, three of which contained mainly species specialized on either Brassicaceae, Solanaceae, or Fagaceae. Most members of the fourth clade have broader host use, including Brassicaceae and Poaceae as major host plant families. Ancestral state reconstructions suggest that Psylliodes flea beetles were initially associated with Brassicaceae and then either shifted to Solanaceae or Fagaceae, or expanded their host repertoire to Poaceae. Despite a putative ancestral association with Brassicaceae, we found evidence that the evolution of glucosinolate-specific detoxification enzymes coincides with the radiation of Psylliodes on Brassicaceae, suggesting that these are not required for using Brassicaceae as hosts but could improve the efficiency of host use by specialized Psylliodes species.

Evolutionary history and divergence times of Odonata (dragonflies and damselflies) revealed through transcriptomics.

Dragonflies and damselflies are among the earliest flying insects with extant representatives. However, unraveling details of their long evolutionary history, such as egg laying (oviposition) strategies, is impeded by unresolved phylogenetic relationships, particularly in damselflies. Here we present a transcriptome-based phylogenetic reconstruction of Odonata, analyzing 2,980 protein-coding genes in 105 species representing nearly all the order's families. All damselfly and most dragonfly families are recovered as monophyletic. Our data suggest a sister relationship between dragonfly families of Gomphidae and Petaluridae. According to our divergence time estimates, both crown-Zygoptera and -Anisoptera arose during the late Triassic. Egg-laying with a reduced ovipositor apparently evolved in dragonflies during the late Jurassic/early Cretaceous. Lastly, we also test the impact of fossil choice and placement, particularly, of the extinct fossil species, †Triassolestodes asiaticus, and †Proterogomphus renateae on divergence time estimates. We find placement of †Proterogomphus renateae to be much more impactful than †Triassolestodes asiaticus.

Author Correction: Phylogenomic analysis sheds light on the evolutionary pathways towards acoustic communication in Orthoptera.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Phylogenomic analysis sheds light on the evolutionary pathways towards acoustic communication in Orthoptera.

Acoustic communication is enabled by the evolution of specialised hearing and sound producing organs. In this study, we performed a large-scale macroevolutionary study to understand how both hearing and sound production evolved and affected diversification in the insect order Orthoptera, which includes many familiar singing insects, such as crickets, katydids, and grasshoppers. Using phylogenomic data, we firmly establish phylogenetic relationships among the major lineages and divergence time estimates within Orthoptera, as well as the lineage-specific and dynamic patterns of evolution for hearing and sound producing organs. In the suborder Ensifera, we infer that forewing-based stridulation and tibial tympanal ears co-evolved, but in the suborder Caelifera, abdominal tympanal ears first evolved in a non-sexual context, and later co-opted for sexual signalling when sound producing organs evolved. However, we find little evidence that the evolution of hearing and sound producing organs increased diversification rates in those lineages with known acoustic communication.

The evolution and genomic basis of beetle diversity.

The order Coleoptera (beetles) is arguably the most speciose group of animals, but the evolutionary history of beetles, including the impacts of plant feeding (herbivory) on beetle diversification, remain poorly understood. We inferred the phylogeny of beetles using 4,818 genes for 146 species, estimated timing and rates of beetle diversification using 89 genes for 521 species representing all major lineages and traced the evolution of beetle genes enabling symbiont-independent digestion of lignocellulose using 154 genomes or transcriptomes. Phylogenomic analyses of these uniquely comprehensive datasets resolved previously controversial beetle relationships, dated the origin of Coleoptera to the Carboniferous, and supported the codiversification of beetles and angiosperms. Moreover, plant cell wall-degrading enzymes (PCWDEs) obtained from bacteria and fungi via horizontal gene transfers may have been key to the Mesozoic diversification of herbivorous beetles-remarkably, both major independent origins of specialized herbivory in beetles coincide with the first appearances of an arsenal of PCWDEs encoded in their genomes. Furthermore, corresponding (Jurassic) diversification rate increases suggest that these novel genes triggered adaptive radiations that resulted in nearly half of all living beetle species. We propose that PCWDEs enabled efficient digestion of plant tissues, including lignocellulose in cell walls, facilitating the evolution of uniquely specialized plant-feeding habits, such as leaf mining and stem and wood boring. Beetle diversity thus appears to have resulted from multiple factors, including low extinction rates over a long evolutionary history, codiversification with angiosperms, and adaptive radiations of specialized herbivorous beetles following convergent horizontal transfers of microbial genes encoding PCWDEs.

Evolutionary history of Polyneoptera and its implications for our understanding of early winged insects.

Polyneoptera represents one of the major lineages of winged insects, comprising around 40,000 extant species in 10 traditional orders, including grasshoppers, roaches, and stoneflies. Many important aspects of polyneopteran evolution, such as their phylogenetic relationships, changes in their external appearance, their habitat preferences, and social behavior, are unresolved and are a major enigma in entomology. These ambiguities also have direct consequences for our understanding of the evolution of winged insects in general; for example, with respect to the ancestral habitats of adults and juveniles. We addressed these issues with a large-scale phylogenomic analysis and used the reconstructed phylogenetic relationships to trace the evolution of 112 characters associated with the external appearance and the lifestyle of winged insects. Our inferences suggest that the last common ancestors of Polyneoptera and of the winged insects were terrestrial throughout their lives, implying that wings did not evolve in an aquatic environment. The appearance of the first polyneopteran insect was mainly characterized by ancestral traits such as long segmented abdominal appendages and biting mouthparts held below the head capsule. This ancestor lived in association with the ground, which led to various specializations including hardened forewings and unique tarsal attachment structures. However, within Polyneoptera, several groups switched separately to a life on plants. In contrast to a previous hypothesis, we found that social behavior was not part of the polyneopteran ground plan. In other traits, such as the biting mouthparts, Polyneoptera shows a high degree of evolutionary conservatism unique among the major lineages of winged insects.

Climate and host-plant associations shaped the evolution of ceutorhynch weevils throughout the Cenozoic.

Using molecular phylogenetic data and methods we inferred divergence times and diversification patterns for the weevil subfamily Ceutorhynchinae in the context of host-plant associations and global climate over evolutionary time. We detected four major diversification shifts that correlate with both host shifts and major climate events. Ceutorhynchinae experienced an increase in diversification rate at ∼53 Ma, during the Early Eocene Climate Optimum, coincident with a host shift to Lamiaceae. A second major diversification phase occurred at the end of the Eocene (∼34 Ma). This contrasts with the overall deterioration in climate equability at the Eocene-Oligocene boundary, but tracks the diversification of important host plant clades in temperate (higher) latitudes, leading to increased diversification rates in the weevil clades infesting temperate hosts. A third major phase of diversification is correlated with the rising temperatures of the Late Oligocene Warming Event (∼26.5 Ma); diversification rates then declined shortly after the Middle Miocene Climate Transition (∼14.9 Ma). Our results indicate that biotic and abiotic factors together explain the evolution of Ceutorhynchinae better than each of these drivers viewed in isolation.

Timing and host plant associations in the evolution of the weevil tribe Apionini (Apioninae, Brentidae, Curculionoidea, Coleoptera) indicate an ancient co-diversification pattern of beetles and flowering plants.

Host plant shifts of insects can lead to a burst of diversification driven by their arrival in a new adaptive zone. In this context, our study aims to explore timing and patterns in the evolution of the weevil tribe Apionini (Brentidae, Curculionoidea, Coleoptera), particularly in relation to affiliations with their host plants. The classification of Apionini is difficult because of their relatively uniform appearance. Most taxa live mono- or oligophagously on members of Asteraceae or Fabaceae, but many are associated with other plant families, like Lamiaceae, Malvaceae and Polygonaceae. However, a comprehensive hypothesis of the phylogenetic relationships within the tribe Apionini is still missing. In the present study, we reconstructed trees and estimated divergence times among tribes. These results were further used to reconstruct the ancestral host plant use in Apionini weevils and to infer if the divergence timing of putative subtribes corresponds with the occurrence and radiation of their specific host plant groups. Phylogenetic analyses confirm the monophyly of most subtribes, with the exceptions of Oxystomatina, Kalcapiina and Aspidapiina. The subribe Aplemonina is inferred to be sister to all remaining Apionini. Divergence time estimates indicate the first occurrence of Apionini in the Upper Cretaceous and a simultaneous occurrence of several families of flowering plants and the occupation by Apionini weevil herbivores. These conspicuous coincidences support either an ancient co-diversification scenario or an escalating diversification in weevils induced by the radiation of flowering plants.

Not going with the flow: a comprehensive time-calibrated phylogeny of dragonflies (Anisoptera: Odonata: Insecta) provides evidence for the role of lentic habitats on diversification.

Ecological diversification of aquatic insects has long been suspected to have been driven by differences in freshwater habitats, which can be classified into flowing (lotic) waters and standing (lentic) waters. The contrasting characteristics of lotic and lentic freshwater systems imply different ecological constraints on their inhabitants. The ephemeral and discontinuous character of most lentic water bodies may encourage dispersal by lentic species in turn reducing geographical isolation among populations. Hence, speciation probability would be lower in lentic species. Here, we assess the impact of habitat use on diversification patterns in dragonflies (Anisoptera: Odonata). Based on the eight nuclear and mitochondrial genes, we inferred species diversification with a model-based evolutionary framework, to account for rate variation through time and among lineages and to estimate the impact of larval habitat on the potentially nonrandom diversification among anisopteran groups. Ancestral state reconstruction revealed lotic fresh water systems as their original primary habitat, while lentic waters have been colonized independently in Aeshnidae, Corduliidae and Libellulidae. Furthermore, our results indicate a positive correlation of speciation and lentic habitat colonization by dragonflies: speciation rates increased in lentic Aeshnidae and Libellulidae, whereas they remain mostly uniform among lotic groups. This contradicts the hypothesis of inherently lower speciation in lentic groups and suggests species with larger ranges are more likely to diversify, perhaps due to higher probability of larger areas being dissected by geographical barriers. Furthermore, larger range sizes may comprise more habitat types, which could also promote speciation by providing additional niches, allowing the coexistence of emerging species.

Tracing the radiation of Maniola (Nymphalidae) butterflies: new insights from phylogeography hint at one single incompletely differentiated species complex.

The use of DNA sequence data often leads to the recognition of cryptic species within putatively well-known taxa. The opposite case, detecting less diversity than originally described, has, however, far more rarely been documented. Maniola jurtina, the Meadow Brown butterfly, occurs all over Europe, whereas all other six species in the genus Maniola are restricted to the Mediterranean area. Among them, three are island endemics on Sardinia, Cyprus, and Chios, respectively. Maniola species are almost indistinguishable morphologically, and hybridization seems to occur occasionally. To clarify species boundaries and diversification history of the genus, we reconstructed the phylogeography and phylogeny of all seven species within Maniola analyzing 138 individuals from across its range using mitochondrial and nuclear genetic markers. Examination of variation in mitochondrial and nuclear DNA surprisingly revealed a case of taxonomic "oversplitting". The topology of the recovered phylogenetic tree is not consistent with accepted taxonomy, but rather reveals haplotype clades that are incongruent with nominal species boundaries: instead of seven species, we recognized only two major, yet incompletely segregated, lineages. Our results are consistent with the hypothesis that Maniola originated in Africa. We suggest that one lineage dispersed over the Strait of Gibraltar and the Iberian Peninsula to the west of Europe, while the other lineage spreads eastward through Asia Minor and over the Bosporus to Eastern Europe.

Phylogenomics resolves the timing and pattern of insect evolution.

Insects are the most speciose group of animals, but the phylogenetic relationships of many major lineages remain unresolved. We inferred the phylogeny of insects from 1478 protein-coding genes. Phylogenomic analyses of nucleotide and amino acid sequences, with site-specific nucleotide or domain-specific amino acid substitution models, produced statistically robust and congruent results resolving previously controversial phylogenetic relations hips. We dated the origin of insects to the Early Ordovician [~479 million years ago (Ma)], of insect flight to the Early Devonian (~406 Ma), of major extant lineages to the Mississippian (~345 Ma), and the major diversification of holometabolous insects to the Early Cretaceous. Our phylogenomic study provides a comprehensive reliable scaffold for future comparative analyses of evolutionary innovations among insects.

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.

Revival of Palaeoptera-head characters support a monophyletic origin of Odonata and Ephemeroptera (Insecta).

The earliest branching event in winged insects, one of the core problems regarding early insect evolution, was addressed using characters of the head. The head is arguably one of the most complex body regions in insects and the phylogenetic information content of its features has been demonstrated. In contrast, the wings and other body parts related to the flight apparatus and sperm transmission are not useful in the context of this problem, as the outgroups (silverfish and bristletails) are wingless and transmit spermatophores externally. Therefore, they show profound differences in the organization of the postcephalic body, and assessment of homology and subsequent phylogenetic analysis of features of these body regions is extremely difficult. The core of this study is the investigation of head structures of representatives of the major clades of dragonflies. A detailed description of the head of Lestes virens is presented and was used as a starting point for the compilation of a character set and a character state matrix for the entire Dicondylia (winged insects + silverfish), with a main focus on the placement of dragonflies and consequently the basal branching event within winged insects. Our results indicate a sister-group relationship between a clade Palaeoptera (dragonflies + mayflies) and the megadiverse monophyletic lineage Neoptera. We show that despite considerable structural similarity between the odonate and neopteran mandible, the muscle equipment in dragonflies is more plesiomorphic with respect to Dicondylia than previously known. Odonata and Ephemeroptera also share presumably derived features of the antenna, maxilla, and labial musculature. Parsimony analyses of the head data unambiguously support a clade Palaeoptera. © The Willi Hennig Society 2012.

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

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

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.

Accurate and efficient reconstruction of deep phylogenies from structured RNAs.

Ribosomal RNA (rRNA) genes are probably the most frequently used data source in phylogenetic reconstruction. Individual columns of rRNA alignments are not independent as a consequence of their highly conserved secondary structures. Unless explicitly taken into account, these correlation can distort the phylogenetic signal and/or lead to gross overestimates of tree stability. Maximum likelihood and Bayesian approaches are of course amenable to using RNA-specific substitution models that treat conserved base pairs appropriately, but require accurate secondary structure models as input. So far, however, no accurate and easy-to-use tool has been available for computing structure-aware alignments and consensus structures that can deal with the large rRNAs. The RNAsalsa approach is designed to fill this gap. Capitalizing on the improved accuracy of pairwise consensus structures and informed by a priori knowledge of group-specific structural constraints, the tool provides both alignments and consensus structures that are of sufficient accuracy for routine phylogenetic analysis based on RNA-specific substitution models. The power of the approach is demonstrated using two rRNA data sets: a mitochondrial rRNA set of 26 Mammalia, and a collection of 28S nuclear rRNAs representative of the five major echinoderm groups.

Simultaneous alignment and folding of 28S rRNA sequences uncovers phylogenetic signal in structure variation.

Secondary structure models of mitochondrial and nuclear (r)RNA sequences are frequently applied to aid the alignment of these molecules in phylogenetic analyses. Additionally, it is often speculated that structure variation of (r)RNA sequences might profitably be used as phylogenetic markers. The benefit of these approaches depends on the reliability of structure models. We used a recently developed approach to show that reliable inference of large (r)RNA secondary structures as a prerequisite of simultaneous sequence and structure alignment is feasible. The approach iteratively establishes local structure constraints of each sequence and infers fully folded individual structures by constrained MFE optimization. A comparison of structure edit distances of individual constraints and fully folded structures showed pronounced phylogenetic signal in fully folded structures. As model sequences we characterized secondary structures of 28S rRNA sequences of selected insects and examined their phylogenetic signal according to established phylogenetic hypotheses.

Phylogenetic support values are not necessarily informative: the case of the Serialia hypothesis (a mollusk phylogeny).

BACKGROUND: Molecular phylogenies are being published increasingly and many biologists rely on the most recent topologies. However, different phylogenetic trees often contain conflicting results and contradict significant background data. Not knowing how reliable traditional knowledge is, a crucial question concerns the quality of newly produced molecular data. The information content of DNA alignments is rarely discussed, as quality statements are mostly restricted to the statistical support of clades. Here we present a case study of a recently published mollusk phylogeny that contains surprising groupings, based on five genes and 108 species, and we apply new or rarely used tools for the analysis of the information content of alignments and for the filtering of noise (masking of random-like alignment regions, split decomposition, phylogenetic networks, quartet mapping). RESULTS: The data are very fragmentary and contain contaminations. We show that that signal-like patterns in the data set are conflicting and partly not distinct and that the reported strong support for a "rather surprising result" (monoplacophorans and chitons form a monophylum Serialia) does not exist at the level of primary homologies. Split-decomposition, quartet mapping and neighbornet analyses reveal conflicting nucleotide patterns and lack of distinct phylogenetic signal for the deeper phylogeny of mollusks. CONCLUSION: Even though currently a majority of molecular phylogenies are being justified with reference to the 'statistical' support of clades in tree topologies, this confidence seems to be unfounded. Contradictions between phylogenies based on different analyses are already a strong indication of unnoticed pitfalls. The use of tree-independent tools for exploratory analyses of data quality is highly recommended. Concerning the new mollusk phylogeny more convincing evidence is needed.

Can comprehensive background knowledge be incorporated into substitution models to improve phylogenetic analyses? A case study on major arthropod relationships.

BACKGROUND: Whenever different data sets arrive at conflicting phylogenetic hypotheses, only testable causal explanations of sources of errors in at least one of the data sets allow us to critically choose among the conflicting hypotheses of relationships. The large (28S) and small (18S) subunit rRNAs are among the most popular markers for studies of deep phylogenies. However, some nodes supported by this data are suspected of being artifacts caused by peculiarities of the evolution of these molecules. Arthropod phylogeny is an especially controversial subject dotted with conflicting hypotheses which are dependent on data set and method of reconstruction. We assume that phylogenetic analyses based on these genes can be improved further i) by enlarging the taxon sample and ii) employing more realistic models of sequence evolution incorporating non-stationary substitution processes and iii) considering covariation and pairing of sites in rRNA-genes. RESULTS: We analyzed a large set of arthropod sequences, applied new tools for quality control of data prior to tree reconstruction, and increased the biological realism of substitution models. Although the split-decomposition network indicated a high noise content in the data set, our measures were able to both improve the analyses and give causal explanations for some incongruities mentioned from analyses of rRNA sequences. However, misleading effects did not completely disappear. CONCLUSION: Analyses of data sets that result in ambiguous phylogenetic hypotheses demand for methods, which do not only filter stochastic noise, but likewise allow to differentiate phylogenetic signal from systematic biases. Such methods can only rely on our findings regarding the evolution of the analyzed data. Analyses on independent data sets then are crucial to test the plausibility of the results. Our approach can easily be extended to genomic data, as well, whereby layers of quality assessment are set up applicable to phylogenetic reconstructions in general.