Mind the Outgroup and Bare Branches in Total-Evidence Dating: a Case Study of Pimpliform Darwin Wasps (Hymenoptera, Ichneumonidae).

Taxon sampling is a central aspect of phylogenetic study design, but it has received limited attention in the context of total-evidence dating, a widely used dating approach that directly integrates molecular and morphological information from extant and fossil taxa. We here assess the impact of commonly employed outgroup sampling schemes and missing morphological data in extant taxa on age estimates in a total-evidence dating analysis under the uniform tree prior. Our study group is Pimpliformes, a highly diverse, rapidly radiating group of parasitoid wasps of the family Ichneumonidae. We analyze a data set comprising 201 extant and 79 fossil taxa, including the oldest fossils of the family from the Early Cretaceous and the first unequivocal representatives of extant subfamilies from the mid-Paleogene. Based on newly compiled molecular data from ten nuclear genes and a morphological matrix that includes 222 characters, we show that age estimates become both older and less precise with the inclusion of more distant and more poorly sampled outgroups. These outgroups not only lack morphological and temporal information but also sit on long terminal branches and considerably increase the evolutionary rate heterogeneity. In addition, we discover an artifact that might be detrimental for total-evidence dating: "bare-branch attraction," namely high attachment probabilities of certain fossils to terminal branches for which morphological data are missing. Using computer simulations, we confirm the generality of this phenomenon and show that a large phylogenetic distance to any of the extant taxa, rather than just older age, increases the risk of a fossil being misplaced due to bare-branch attraction. After restricting outgroup sampling and adding morphological data for the previously attracting, bare branches, we recover a Jurassic origin for Pimpliformes and Ichneumonidae. This first age estimate for the group not only suggests an older origin than previously thought but also that diversification of the crown group happened well before the Cretaceous-Paleogene boundary. Our case study demonstrates that in order to obtain robust age estimates, total-evidence dating studies need to be based on a thorough and balanced sampling of both extant and fossil taxa, with the aim of minimizing evolutionary rate heterogeneity and missing morphological information. [Bare-branch attraction; ichneumonids; fossils; morphological matrix; phylogeny; RoguePlots.].

Plant interactions shape pollination networks via nonadditive effects.

Plants grow in communities where they interact with other plants and with other living organisms such as pollinators. On the one hand, studies of plant-plant interactions rarely consider how plants interact with other trophic levels such as pollinators. On the other, studies of plant-animal interactions rarely deal with interactions within trophic levels such as plant-plant competition and facilitation. Thus, to what degree plant interactions affect biodiversity and ecological networks across trophic levels is poorly understood. We manipulated plant communities driven by foundation species facilitation and sampled plant-pollinator networks at fine spatial scale in a field experiment in Sierra Nevada, Spain. We found that plant-plant facilitation shaped pollinator diversity and structured pollination networks. Nonadditive effects of plant interactions on pollinator diversity and interaction diversity were synergistic in one foundation species networks while they were additive in another foundation species. Nonadditive effects of plant interactions were due to rewiring of pollination interactions. In addition, plant facilitation had negative effects on the structure of pollination networks likely due to increase in plant competition for pollination. Our results empirically demonstrate how different network types are coupled, revealing pervasive consequences of interaction chains in diverse communities.

More on the Best Evolutionary Rate for Phylogenetic Analysis.

The accumulation of genome-scale molecular data sets for nonmodel taxa brings us ever closer to resolving the tree of life of all living organisms. However, despite the depth of data available, a number of studies that each used thousands of genes have reported conflicting results. The focus of phylogenomic projects must thus shift to more careful experimental design. Even though we still have a limited understanding of what are the best predictors of the phylogenetic informativeness of a gene, there is wide agreement that one key factor is its evolutionary rate; but there is no consensus as to whether the rates derived as optimal in various analytical, empirical, and simulation approaches have any general applicability. We here use simulations to infer optimal rates in a set of realistic phylogenetic scenarios with varying tree sizes, numbers of terminals, and tree shapes. Furthermore, we study the relationship between the optimal rate and rate variation among sites and among lineages. Finally, we examine how well the predictions made by a range of experimental design methods correlate with the observed performance in our simulations.We find that the optimal level of divergence is surprisingly robust to differences in taxon sampling and even to among-site and among-lineage rate variation as often encountered in empirical data sets. This finding encourages the use of methods that rely on a single optimal rate to predict a gene's utility. Focusing on correct recovery either of the most basal node in the phylogeny or of the entire topology, the optimal rate is about 0.45 substitutions from root to tip in average Yule trees and about 0.2 in difficult trees with short basal and long-apical branches, but all rates leading to divergence levels between about 0.1 and 0.5 perform reasonably well.Testing the performance of six methods that can be used to predict a gene's utility against our simulation results, we find that the probability of resolution, signal-noise analysis, and Fisher information are good predictors of phylogenetic informativeness, but they require specification of at least part of a model tree. Likelihood quartet mapping also shows very good performance but only requires sequence alignments and is thus applicable without making assumptions about the phylogeny. Despite them being the most commonly used methods for experimental design, geometric quartet mapping and the integration of phylogenetic informativeness curves perform rather poorly in our comparison. Instead of derived predictors of phylogenetic informativeness, we suggest that the number of sites in a gene that evolve at near-optimal rates (as inferred here) could be used directly to prioritize genes for phylogenetic inference. In combination with measures of model fit, especially with respect to compositional biases and among-site and among-lineage rate variation, such an approach has the potential to greatly improve marker choice and should be tested on empirical data.

Post-Eocene climate change across continental Australia and the diversification of Australasian spiny trapdoor spiders (Idiopidae: Arbanitinae).

The formation and spread of the Australian arid zone during the Neogene was a profoundly transformative event in the biogeographic history of Australia, resulting in extinction or range contraction in lineages adapted to mesic habitats, as well as diversification and range expansion in arid-adapted taxa (most of which evolved from mesic ancestors). However, the geographic origins of the arid zone biota are still relatively poorly understood, especially among highly diverse invertebrate lineages, many of which are themselves poorly documented at the species level. Spiny trapdoor spiders (Idiopidae: Arbanitinae) are one such lineage, having mesic 'on-the-continent' Gondwanan origins, while also having experienced major arid zone radiations in select clades. In this study, we present new orthologous nuclear markers for the phylogenetic inference of mygalomorph spiders, and use them to infer the phylogeny of Australasian Idiopidae with a 12-gene parallel tagged amplicon next-generation sequencing approach. We use these data to test the mode and timing of diversification of arid-adapted idiopid lineages across mainland Australia, and employ a continent-wide sampling of the fauna's phylogenetic and geographic diversity to facilitate ancestral area inference. We further explore the evolution of phenotypic and behavioural characters associated with both arid and mesic environments, and test an 'out of south-western Australia' hypothesis for the origin of arid zone clades. Three lineages of Idiopidae are shown to have diversified in the arid zone during the Miocene, one (genus Euoplos) exclusively in Western Australia. Arid zone Blakistonia likely had their origins in South Australia, whereas in the most widespread genus Aganippe, a more complex scenario is evident, with likely range expansion from southern Western Australia to southern South Australia, from where the bulk of the arid zone fauna then originated. In Aganippe, remarkable adaptations to phragmotic burrow-plugging in transitional arid zone taxa have evolved twice independently in Western Australia, while in Misgolas and Cataxia, burrow door-building behaviours have likely been independently lost at least three times in the eastern Australian mesic zone. We also show that the presence of idiopids in New Zealand (Cantuaria) is likely to be the result of recent dispersal from Australia, rather than ancient continental vicariance. By providing the first comprehensive, continental synopsis of arid zone biogeography in an Australian arachnid lineage, we show that the diversification of arbanitine Idiopidae was intimately associated with climate shifts during the Neogene, resulting in multiple Mio-Pliocene radiations.

Total-Evidence Dating under the Fossilized Birth-Death Process.

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.

Molecular and morphological phylogenetics of chelonine parasitoid wasps (Hymenoptera: Braconidae), with a critical assessment of divergence time estimations.

Parasitoid wasps of the subfamily Cheloninae are both species rich and poorly known. Although the taxonomy of Cheloninae appears to be relatively stable, there is no clear understanding of relationships among higher-level taxa. We here applied molecular phylogenetic analyses using three markers (COI, EF1α, 28S) and 37 morphological characters to elucidate the evolution and systematics of these wasps. Analyses were based on 83 specimens representing 13 genera. All genera except Ascogaster, Phanerotoma, and Pseudophanerotoma formed monophyletic groups; Furcidentia (stat. rev.) is raised to generic rank. Neither Chelonus (Chelonus) nor Chelonus (Microchelonus) were recovered as monophyletic, but together formed a monophyletic lineage. The tribes Chelonini and Odontosphaeropygini formed monophyletic groups, but the Phanerotomini sensu Zettel and Pseudophanerotomini were retrieved as either para- or polyphyletic. The genera comprising the former subfamily Adeliinae were confirmed as being nested within the Cheloninae. To estimate the age of the subfamily, we used 16 fossil taxa. Three approaches were compared: fixed-rate dating, node dating, and total-evidence dating, with age estimates differing greatly between the three methods. Shortcomings of each approach in relation to our dataset are discussed, and none of the age estimates is deemed sufficiently reliable. Given that most dating studies use a single method only, in most cases without presenting analyses on the sensitivity to priors, it is likely that numerous age estimates in the literature suffer from a similar lack of robustness. We argue for a more rigorous approach to dating analyses and for a faithful presentation of uncertainties in divergence time estimates. Given the results of the phylogenetic analysis the following taxonomic changes are proposed: Furcidentia Zettel (stat. rev.), previously treated as a subgenus of Pseudophanerotoma Zettel is raised to generic rank; Microchelonus Szépligeti (syn. nov.), variously treated by previous authors, is proposed as a junior synonym of Chelonus Jurine; the following subgenera of Microchelonus - Baculonus Braet & van Achterberg (syn. nov.), Carinichelonus Tobias (syn. nov.) and Scabrichelonus He, Chen & van Achterberg (syn. nov.), are proposed as junior synonyms of Chelonus; a number of new species names are proposed due to homonyms resulting from the above changes and these are listed in the paper.

A Nonstationary Markov Model Detects Directional Evolution in Hymenopteran Morphology.

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.

Does early learning drive ecological divergence during speciation processes in parasitoid wasps?

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.

Dating the diversification of the major lineages of Passeriformes (Aves).

BACKGROUND: The avian Order Passeriformes is an enormously species-rich group, which comprises almost 60% of all living bird species. This diverse order is believed to have originated before the break-up of Gondwana in the late Cretaceous. However, previous molecular dating studies have relied heavily on the geological split between New Zealand and Antarctica, assumed to have occurred 85-82 Mya, for calibrating the molecular clock and might thus be circular in their argument. RESULTS: This study provides a time-scale for the evolution of the major clades of passerines using seven nuclear markers, five taxonomically well-determined passerine fossils, and an updated interpretation of the New Zealand split from Antarctica 85-52 Mya in a Bayesian relaxed-clock approach. We also assess how different interpretations of the New Zealand-Antarctica vicariance event influence our age estimates. Our results suggest that the diversification of Passeriformes began in the late Cretaceous or early Cenozoic. Removing the root calibration for the New Zealand-Antarctica vicariance event (85-52 Mya) dramatically increases the 95% credibility intervals and leads to unrealistically old age estimates. We assess the individual characteristics of the seven nuclear genes analyzed in our study. Our analyses provide estimates of divergence times for the major groups of passerines, which can be used as secondary calibration points in future molecular studies. CONCLUSIONS: Our analysis takes recent paleontological and geological findings into account and provides the best estimate of the passerine evolutionary time-scale currently available. This time-scale provides a temporal framework for further biogeographical, ecological, and co-evolutionary studies of the largest bird radiation, and adds to the growing support for a Cretaceous origin of Passeriformes.

Revision of the Western Palaearctic Diplazontinae (Hymenoptera, Ichneumonidae).

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 exemplifies an integrative approach to systematics and taxonomy.

Review of the Diplazontinae (Hymenoptera, Ichneumonidae) of the Kuril islands, with descriptions of two new species.

A sample of 346 specimens of Diplazontinae parasitoid wasps from the Kuril islands was studied. Twenty-six species are reported, Tymmophorus gelidus Dasch for the first time for the Eastern Palaearctic. Two new species are described, Diplazon kurilensis sp. n. and Homotropus formosus sp. n. Diplazon urupensis Uchida is removed from synonymy (stat. rev.), and Promethes persulcatus Nakanishi is suggested as a synonym of Promethes bridgmani Fitton. Reasons are discussed for the large proportion of species with a Holarctic or even multi-regional distribution in the sample, which amounts to17 of the species or 65%.

A macroevolutionary role for chromosomal fusion and fission in Erebia butterflies.

The impact of large-scale chromosomal rearrangements, such as fusions and fissions, on speciation is a long-standing conundrum. We assessed whether bursts of change in chromosome numbers resulting from chromosomal fusion or fission are related to increased speciation rates in Erebia, one of the most species-rich and karyotypically variable butterfly groups. We established a genome-based phylogeny and used state-dependent birth-death models to infer trajectories of karyotype evolution. We demonstrated that rates of anagenetic chromosomal changes (i.e., along phylogenetic branches) exceed cladogenetic changes (i.e., at speciation events), but, when cladogenetic changes occur, they are mostly associated with chromosomal fissions rather than fusions. We found that the relative importance of fusion and fission differs among Erebia clades of different ages and that especially in younger, more karyotypically diverse clades, speciation is more frequently associated with cladogenetic chromosomal changes. Overall, our results imply that chromosomal fusions and fissions have contrasting macroevolutionary roles and that large-scale chromosomal rearrangements are associated with bursts of species diversification.

The evolution of host associations in the parasitic wasp genus Ichneumon (Hymenoptera: Ichneumonidae): convergent adaptations to host pupation sites.

BACKGROUND: The diversification of organisms with a parasitic lifestyle is often tightly linked to the evolution of their host associations. If a tight host association exists, closely related species tend to attack closely related hosts; host associations are less stable if associations are determined by more plastic traits like parasitoid searching and oviposition behaviour. The pupal-parasitoids of the genus Ichneumon attack a variety of macrolepidopteran hosts. They are either monophagous or polyphagous, and therefore offer a promissing system to investigate the evolution of host associations. Ichneumon was previously divided into two groups based on general body shape; however, a stout shape has been suggested as an adaptation to buried host pupation sites, and might thus not represent a reliable phylogenetic character. RESULTS: We here reconstruct the first molecular phylogeny of the genus Ichneumon using two mitochondrial (CO1 and NADH1) and one nuclear marker (28S). The resulting phylogeny only supports monophyly of Ichneumon when Ichneumon lugens Gravenhorst, 1829 (formerly in Chasmias, stat. rev.) and Ichneumon deliratorius Linnaeus, 1758 (formerly Coelichneumon) are included. Neither parasitoid species that attack hosts belonging to one family nor those attacking butterflies (Rhopalocera) form monophyletic clades. Ancestral state reconstructions suggest multiple transitions between searching for hosts above versus below ground and between a stout versus elongated body shape. A model assuming correlated evolution between the two characters was preferred over independent evolution of host-searching niche and body shape. CONCLUSIONS: Host relations, both in terms of phylogeny and ecology, evolved at a high pace in the genus Ichneumon. Numerous switches between hosts of different lepidopteran families have occurred, a pattern that seems to be the rule among idiobiont parasitoids. A stout body and antennal shape in the parasitoid female is confirmed as an ecological adaptation to host pupation sites below ground and has evolved convergently several times. Morphological characters that might be involved in adaptation to hosts should be avoided as diagnostic characters for phylogeny and classification, as they can be expected to show high levels of homoplasy.

Convergent intron gains in hymenopteran elongation factor-1α.

The eukaryotic translation elongation factor-1α gene (eEF1A) has been used extensively in higher level phylogenetics of insects and other groups, despite being present in two or more copies in several taxa. Orthology assessment has relied heavily on the position of introns, but the basic assumption of low rates of intron loss and absence of convergent intron gains has not been tested thoroughly. Here, we study the evolution of eEF1A based on a broad sample of taxa in the insect order Hymenoptera. The gene is universally present in two copies - F1 and F2 - both of which apparently originated before the emergence of the order. An elevated ratio of non-synonymous versus synonymous substitutions and differences in rates of amino acid replacements between the copies suggest that they evolve independently, and phylogenetic methods clearly cluster the copies separately. The F2 copy appears to be ancient; it is orthologous with the copy known as F1 in Diptera, and is likely present in most insect orders. The hymenopteran F1 copy, which may or may not be unique to this order, apparently originated through retroposition and was originally intron free. During the evolution of the Hymenoptera, it has successively accumulated introns, at least three of which have appeared at the same position as introns in the F2 copy or in eEF1A copies in other insects. The sites of convergent intron gain are characterized by highly conserved nucleotides that strongly resemble specific intron-associated sequence motifs, so-called proto-splice sites. The significant rate of convergent intron gain renders intron-exon structure unreliable as an indicator of orthology in eEF1A, and probably also in other protein-coding genes.

A total-evidence approach to dating with fossils, applied to the early radiation of the hymenoptera.

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

Impact of increasing morphological information by micro-CT scanning on the phylogenetic placement of Darwin wasps (Hymenoptera, Ichneumonidae) in amber.

The correct interpretation of fossils and their reliable taxonomic placements are fundamental for understanding the evolutionary history of biodiversity. Amber inclusions often preserve more morphological information than compression fossils, but are often partially hidden or distorted, which can impede taxonomic identification. Here, we studied four new fossil species of Darwin wasps from Baltic and Dominican amber, using micro computed tomography (micro-CT) scans and 3D reconstructions to accurately interpret and increase the availability of morphological information. We then infer their taxonomic placement in a Bayesian phylogenetic analysis by combining morphological and molecular data of extant and fossil Darwin wasps and evaluate the impact and usefulness of the additional information from micro-CT scanning. The results show that although we gained significant morphological information from micro-CT scanning, especially concerning measurements and hidden dorsal and ventral structures, this did not impact subfamily-level placement for any of the four fossils. However, micro-CT scanning improved the precision of fossil placements at the genus level, which might be key in future dating and diversification analyses. Finally, we describe the four new fossil species as Rhyssa gulliveri sp. nov. in Rhyssinae, Triclistus levii sp. nov. in Metopiinae, Firkantus freddykruegeri gen. et. sp. nov. in Pimplinae and Magnocula sarcophaga gen. et sp. nov. in Phygadeuontinae. The first two species are the first known representatives of the subfamilies Rhyssinae and Metopiinae in amber. Supplementary Information: The online version contains supplementary material available at 10.1186/s13358-023-00294-2.

Classifying fossil Darwin wasps (Hymenoptera: Ichneumonidae) with geometric morphometrics of fore wings.

Linking fossil species to the extant diversity is often a difficult task, and the correct interpretation of character evidence is crucial for assessing their taxonomic placement. Here, we make use of geometric morphometrics of fore wings to help classify five fossil Darwin wasps from the Early Eocene Fur Formation in Denmark into subfamilies and often tribes. We compile a reference dataset with 342 fore wings of nine extant subfamilies and nine relevant fossil species. Since geometric morphometrics was mostly ignored in the past in Darwin wasp classification, the dataset is first used to examine differences and similarities in wing venation among subfamilies. In a next step, we used the reference dataset to inform the classification of the fossil species, which resulted in the description of one new genus and five new species, Crusopimpla weltii sp. nov., Ebriosa flava gen. et sp. nov., Entypoma? duergari sp. nov., Lathrolestes? zlatorog sp. nov., and Triclistus bibori sp. nov., in four different subfamilies. Carefully assessing data quality, we show that the fore wing venation of fossil Darwin wasps is surprisingly suitable to assign them to a subfamily or even lower taxonomic level, especially when used in conjunction with characters from other parts of the body to narrow down a candidate set of potential subfamilies and tribes. Our results not only demonstrate a fast and useful approach to inform fossil classification but provide a basis for future investigations into evolutionary changes in fore wings of ichneumonids. The high informativeness of wing venation for classification furthermore could be harvested for phylogenetic analyses, which are otherwise often hampered by homoplasy in this parasitoid wasp family.

Many evolutionary roads led to virus domestication in ichneumonoid parasitoid wasps.

The investigation of endogenous viral elements (EVEs) has historically focused on only a few lineages of parasitoid wasps, with negative results consistently underreported. Recent studies show that multiple viral lineages were integrated in at least seven instances in Ichneumonoidea and may be much more widespread than previously thought. Increasingly affordable genomic and bioinformatic approaches have made it feasible to search for viral sequences within wasp genomes, opening an extremely promising research avenue. Advances in wasp phylogenetics have shed light on the evolutionary history of EVE integration, although many questions remain. Phylogenetic proximity can be used as a guide to facilitate targeted screening, to estimate the number and age of integration events and to identify taxa involved in major host switches.

An evaluation of phylogenetic informativeness profiles and the molecular phylogeny of diplazontinae (Hymenoptera, Ichneumonidae).

How to quantify the phylogenetic information content of a data set is a longstanding question in phylogenetics, influencing both the assessment of data quality in completed studies and the planning of future phylogenetic projects. Recently, a method has been developed that profiles the phylogenetic informativeness (PI) of a data set through time by linking its site-specific rates of change to its power to resolve relationships at different timescales. Here, we evaluate the performance of this method in the case of 2 standard genetic markers for phylogenetic reconstruction, 28S ribosomal RNA and cytochrome oxidase subunit 1 (CO1) mitochondrial DNA, with maximum parsimony, maximum likelihood, and Bayesian analyses of relationships within a group of parasitoid wasps (Hymenoptera: Ichneumonidae, Diplazontinae). Retrieving PI profiles of the 2 genes from our own and from 3 additional data sets, we find that the method repeatedly overestimates the performance of the more quickly evolving CO1 compared with 28S. We explore possible reasons for this bias, including phylogenetic uncertainty, violation of the molecular clock assumption, model misspecification, and nonstationary nucleotide composition. As none of these provides a sufficient explanation of the observed discrepancy, we use simulated data sets, based on an idealized setting, to show that the optimum evolutionary rate decreases with increasing number of taxa. We suggest that this relationship could explain why the formula derived from the 4-taxon case overrates the performance of higher versus lower rates of evolution in our case and that caution should be taken when the method is applied to data sets including more than 4 taxa.

Incongruent molecular and morphological variation in the crab spider Synemaglobosum (Araneae, Thomisidae) in Europe.

Establishing species boundaries is one of the challenges taxonomists around the world have been tackling for centuries. The relation between intraspecific and interspecific variability is still under discussion and in many taxa it remains understudied. Here the hypothesis of single versus multiple species of the crab spider Synemaglobosum (Fabricius) is tested. The wide distribution range as well as its high morphological variability makes this species an interesting candidate for re-evaluation using an integrative approach. This study combines information from barcoding, phylogenetic reconstruction based on mitochondrial CO1 and ITS2 of more than 60 specimens collected over a wide range of European localities, and morphology. The findings show deep clades with up to 6% mean pairwise distance in the CO1 barcode without any biogeographical pattern. The nuclear ITS2 gene did not support the CO1 clades. Morphological assessment of somatic and genital characters in males and females and a morphometric analysis of the male palp uncovered high intraspecific variation that does not match the CO1 or ITS2 phylogenies or biogeography either. Screening for endosymbiotic Wolbachia bacteria was conducted and only a single infected specimen was found. Several scenarios might explain these inconsistent patterns. While the deep divergences in the barcoding marker might suggest cryptic or ongoing speciation or geographical isolation in the past, the lack of congruent variation in the nuclear ITS2 gene or the studied morphological character systems, especially the male palp, indicates that S.globosum might simply be highly polymorphic both in terms of its mtDNA and morphology. Therefore, more data on ecology and behaviour and full genome sequences are necessary to ultimately resolve this taxonomically intriguing case.