Phylogenies are usually dated by calibrating interior nodes against the fossil record. This relies on indirect methods that, in the worst case, misrepresent the fossil information. Here, we contrast such node dating with an approach that includes fossils along with the extant taxa in a Bayesian total-evidence analysis. As a test case, we focus on the early radiation of the Hymenoptera, mostly documented by poorly preserved impression fossils that are difficult to place phylogenetically. Specifically, we compare node dating using nine calibration points derived from the fossil record with total-evidence dating based on 343 morphological characters scored for 45 fossil (4--20 complete) and 68 extant taxa. In both cases we use molecular data from seven markers (∼5 kb) for the extant taxa. Because it is difficult to model speciation, extinction, sampling, and fossil preservation realistically, we develop a simple uniform prior for clock trees with fossils, and we use relaxed clock models to accommodate rate variation across the tree. Despite considerable uncertainty in the placement of most fossils, we find that they contribute significantly to the estimation of divergence times in the total-evidence analysis. In particular, the posterior distributions on divergence times are less sensitive to prior assumptions and tend to be more precise than in node dating. The total-evidence analysis also shows that four of the seven Hymenoptera calibration points used in node dating are likely to be based on erroneous or doubtful assumptions about the fossil placement. With respect to the early radiation of Hymenoptera, our results suggest that the crown group dates back to the Carboniferous, ∼309 Ma (95% interval: 291--347 Ma), and diversified into major extant lineages much earlier than previously thought, well before the Triassic. [Bayesian inference; fossil dating; morphological evolution; relaxed clock; statistical phylogenetics.]

Many species, including humans, have dramatically expanded their range in the past, and such range expansions had certainly an impact on their genetic diversity. For example, mutations arising in populations at the edge of a range expansion can sometimes surf on the wave of advance and thus reach a larger spatial distribution and a much higher frequency than would be expected in stationary populations. We study here this surfing phenomenon in more detail, by performing extensive computer simulations under a two-dimensional stepping-stone model. We find that the probability of survival of a new mutation depends to a large degree on its proximity to the edge of the wave. Demographic factors such as deme size, migration rate, and local growth rate also influence the fate of these new mutations. We also find that the final spatial and frequency distributions depend on the local deme size of a subdivided population. This latter result is discussed in the light of human expansions in Europe as it should allow one to distinguish between mutations having spread with Paleolithic or Neolithic expansions. By favoring the spread of new mutations, a consequence of the surfing phenomenon is to increase the rate of evolution of spatially expanding populations.

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

Bayesian total-evidence dating involves the simultaneous analysis of morphological data from the fossil record and morphological and sequence data from recent organisms, and it accommodates the uncertainty in the placement of fossils while dating the phylogenetic tree. Due to the flexibility of the Bayesian approach, total-evidence dating can also incorporate additional sources of information. Here, we take advantage of this and expand the analysis to include information about fossilization and sampling processes. Our work is based on the recently described fossilized birth–death (FBD) process, which has been used to model speciation, extinction, and fossilization rates that can vary over time in a piecewise manner. So far, sampling of extant and fossil taxa has been assumed to be either complete or uniformly at random, an assumption which is only valid for a minority of data sets. We therefore extend the FBD process to accommodate diversified sampling of extant taxa, which is standard practice in studies of higher-level taxa. We verify the implementation using simulations and apply it to the early radiation of Hymenoptera (wasps, ants, and bees). Previous total-evidence dating analyses of this data set were based on a simple uniform tree prior and dated the initial radiation of extant Hymenoptera to the late Carboniferous (309 Ma). The analyses using the FBD prior under diversified sampling, however, date the radiation to the Triassic and Permian (252 Ma), slightly older than the age of the oldest hymenopteran fossils. By exploring a variety of FBD model assumptions, we show that it is mainly the accommodation of diversified sampling that causes the push toward more recent divergence times. Accounting for diversified sampling thus has the potential to close the long-discussed gap between rocks and clocks. We conclude that the explicit modeling of fossilization and sampling processes can improve divergence time estimates, but only if all important model aspects, including sampling biases, are adequately addressed.

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.

I revise the subfamily Diplazontinae to include 99 Western Palaearctic species, review morphological characters useful for species delimitation and identification, and clarify the status of some morphologically similar taxa using molecular approaches. Illustrated, dichotomous keys to the Western Palaearctic genera and species of the subfamily are presented, and the utility of the molecular markers CO1 and ITS2 for species delimitation in Diplazontinae is discussed. Seven new species are described, Diplazon flixi sp. nov., Diplazon nordicus sp. nov., Diplazon parvus sp. nov., Diplazon zetteli sp. nov., Eurytyloides umbrinus sp. nov., Sussaba roberti sp. nov., and Woldstedtius bauri sp. nov. The Nearctic Sussaba cultriformis (Ashmead), formerly a subspecies of Sussaba dorsalis (Holmgren), is raised to species rank. The following taxa are valid species and hereby removed from synonymy: Episemura ensata (Bauer), stat. rev.; Homotropus frontorius (Thunberg), stat. rev.; Syrphoctonus desvignesii (Marshall), stat. rev.; Syrphophilus scabriculus (Holmgren), stat. rev.; and Tymmophorus suspiciosus (Brischke), stat. rev. Nineteen new synonyms are established: Bioblapsis mallochi Rotheray of Bioblapsis cultiformis (Davis), syn. nov.; Bioblapsis tricincta Ashmead of Syrphophilus scabriculus (Holmgren), syn. nov.; Diplazon bachmaieri Diller of Diplazon angustus Dasch, syn. nov.; Diplazon fechteri Diller of Diplazon cascadensis Dasch, syn. nov.; Homocidus brevis Hedwig of Homotropus pictus (Gravenhorst); Homocidus rubiginosum Schmiedeknecht of Enizemum scutellare (Lange), syn. nov.; Homocidus simulans Stelfox of Homotropus collinus (Stelfox), syn. nov.; Homotropus crassicrus Thomson and Homotropus nudus Dasch of Homotropus dimidiatus (Schrank), syn. nov.; Homocidus asyntactus Schmiedeknecht of Homotropus crassicornis Thomson, syn. nov.; Homocidus subopacus Stelfox and Homotropus quadrangularis Dasch of Homotropus frontorius (Thunberg), syn. nov.; Homocidus impolitus Stelfox of Homotropus pallipes (Gravenhorst), syn. nov.; Homotropus incisus Thomson and Homotropus reflexus Morley of Homotropus pectoralis (Provancher), syn. nov.; Tryphon nigricornis Zetterstedt, a former synonym of H. dimidiatus Schrank, is a synonym of Homotropus pictus (Gravenhorst), syn. nov.; Homotropus fraudulentus Dasch and Homotropus neopulcher Horstmann of Syrphoctonus desvignesii (Marshall), syn. nov.; Homotropus eximius Habermehl of Syrphoctonus tarsatorius (Panzer), syn. nov. The following new combinations are established: Bioblapsis cultiformis (Davis), comb. nov.; Homotropus collinus (Stelfox), comb. nov.; Homotropus dimidiatus (Schrank), comb. nov.; Homotropus frontorius (Thunberg), comb. nov.; Homotropus pectoralis (Provancher), comb. nov.; Homotropus strigator (Fabricius), comb. nov.; Homotropus sundevalli (Holmgren), comb. nov. The present revision is the first comprehensive treatment of the Western Palaearctic Diplazontinae, provides the basis for taxonomic, faunistic, ecological and evolutionary studies in these hoverfly parasitoids, and exemp...

Directional evolution has played an important role in shaping the morphological, ecological, and molecular diversity of life. However, standard substitution models assume stationarity of the evolutionary process over the time scale examined, thus impeding the study of directionality. Here we explore a simple, nonstationary model of evolution for discrete data, which assumes that the state frequencies at the root differ from the equilibrium frequencies of the homogeneous evolutionary process along the rest of the tree (i.e., the process is nonstationary, nonreversible, but homogeneous). Within this framework, we develop a Bayesian approach for testing directional versus stationary evolution using a reversible-jump algorithm. Simulations show that when only data from extant taxa are available, the success in inferring directionality is strongly dependent on the evolutionary rate, the shape of the tree, the relative branch lengths, and the number of taxa. Given suitable evolutionary rates (0.1–0.5 expected substitutions between root and tips), accounting for directionality improves tree inference and often allows correct rooting of the tree without the use of an outgroup. As an empirical test, we apply our method to study directional evolution in hymenopteran morphology. We focus on three character systems: wing veins, muscles, and sclerites. We find strong support for a trend toward loss of wing veins and muscles, while stationarity cannot be ruled out for sclerites. Adding fossil and time information in a total-evidence dating approach, we show that accounting for directionality results in more precise estimates not only of the ancestral state at the root of the tree, but also of the divergence times. Our model relaxes the assumption of stationarity and reversibility by adding a minimum of additional parameters, and is thus well suited to studying the nature of the evolutionary process in data sets of limited size, such as morphology and ecology.

Central to the concept of ecological speciation is the evolution of ecotypes, i.e. groups of individuals occupying different ecological niches. However, the mechanisms behind the first step of separation, the switch of individuals into new niches, are unclear. One long-standing hypothesis, which was proposed for insects but never tested, is that early learning causes new ecological preferences, leading to a switch into a new niche within one generation. Here, we show that a host switch occurred within a parasitoid wasp, which is associated with the ability for early learning and the splitting into separate lineages during speciation. Lariophagus distinguendus consists of two genetically distinct lineages, most likely representing different species. One attacks drugstore beetle larvae (Stegobium paniceum (L.)), which were probably the ancestral host of both lineages. The drugstore beetle lineage has an innate host preference that cannot be altered by experience. In contrast, the second lineage is found on Sitophilus weevils as hosts and changes its preference by early learning. We conclude that a host switch has occurred in the ancestor of the second lineage, which must have been enabled by early learning. Because early learning is widespread in insects, it might have facilitated ecological divergence and associated speciation in this hyperdiverse group.