Johannes Bergsten

A review of long-branch attraction

The history of long‐branch attraction, and in particular methods suggested to detect and avoid the artifact to date, is reviewed. Methods suggested to avoid LBA‐artifacts include excluding long‐branch taxa, excluding faster evolving third codon positions, using inference methods less sensitive to LBA such as likelihood, the Aguinaldo et al. approach, sampling more taxa to break up long branches and sampling more characters especially of another kind, and the pros and cons of these are discussed. Methods suggested to detect LBA are numerous and include methodological disconcordance, RASA, separate partition analyses, parametric simulation, random outgroup sequences, long‐branch extraction, split decomposition and spectral analysis. Less than 10 years ago it was doubted if LBA occurred in real datasets. Today, examples are numerous in the literature and it is argued that the development of methods to deal with the problem is warranted. A 16 kbp dataset of placental mammals and a morphological and molecular combined dataset of gall waSPS are used to illustrate the particularly common problem of LBA of problematic ingroup taxa to outgroups. The preferred methods of separate partition analysis, methodological disconcordance, and long branch extraction are used to demonstrate detection methods. It is argued that since outgroup taxa almost always represent long branches and are as such a hazard towards misplacing long branched ingroup taxa, phylogenetic analyses should always be run with and without the outgroups included. This will detect whether only the outgroup roots the ingroup or if it simultaneously alters the ingroup topology, in which case previous studies have shown that the latter is most often the worse. Apart from that LBA to outgroups is the major and most common problem; scanning the literature also detected the ill advised comfort of high support values from thousands of characters, but very few taxa, in the age of genomics. Taxon sampling is crucial for an accurate phylogenetic estimate and trust cannot be put on whole mitochondrial or chloroplast genome studies with only a few taxa, despite their high support values. The placental mammal example demonstrates that parsimony analysis will be prone to LBA by the attraction of the tenrec to the distant marsupial outgroups. In addition, the murid rodents, creating the classic “the guinea‐pig is not a rodent” hypothesis in 1996, are also shown to be attracted to the outgroup by nuclear genes, although including the morphological evidence for rodents and Glires overcomes the artifact. The gall wasp example illustrates that Bayesian analyses with a partition‐specific GTR + Γ + I model give a conflicting resolution of clades, with a posterior probability of 1.0 when comparing ingroup alone versus outgroup rooted topologies, and this is due to long‐branch attraction to the outgroup.

A Comprehensive Phylogeny of Beetles Reveals the Evolutionary Origins of a Superradiation

Beetles represent almost one-fourth of all described species, and knowledge about their relationships and evolution adds to our understanding of biodiversity. We performed a comprehensive phylogenetic analysis of Coleoptera inferred from three genes and nearly 1900 species, representing more than 80% of the worlds recognized beetle families. We defined basal relationships in the Polyphaga supergroup, which contains over 300,000 species, and established five families as the earliest branching lineages. By dating the phylogeny, we found that the success of beetles is explained neither by exceptional net diversification rates nor by a predominant role of herbivory and the Cretaceous rise of angiosperms. Instead, the pre-Cretaceous origin of more than 100 present-day lineages suggests that beetle species richness is due to high survival of lineages and sustained diversification in a variety of niches.

The Effect of Geographical Scale of Sampling on DNA Barcoding

Abstract Eight years after DNA barcoding was formally proposed on a large scale, CO1 sequences are rapidly accumulating from around the world. While studies to date have mostly targeted local or regional species assemblages, the recent launch of the global iBOL project (International Barcode of Life), highlights the need to understand the effects of geographical scale on Barcodings goals. Sampling has been central in the debate on DNA Barcoding, but the effect of the geographical scale of sampling has not yet been thoroughly and explicitly tested with empirical data. Here, we present a CO1 data set of aquatic predaceous diving beetles of the tribe Agabini, sampled throughout Europe, and use it to investigate how the geographic scale of sampling affects 1) the estimated intraspecific variation of species, 2) the genetic distance to the most closely related heterospecific, 3) the ratio of intraspecific and interspecific variation, 4) the frequency of taxonomically recognized species found to be monophyletic, and 5) query identification performance based on 6 different species assignment methods. Intraspecific variation was significantly correlated with the geographical scale of sampling (R-square = 0.7), and more than half of the species with 10 or more sampled individuals (N = 29) showed higher intraspecific variation than 1% sequence divergence. In contrast, the distance to the closest heterospecific showed a significant decrease with increasing geographical scale of sampling. The average genetic distance dropped from > 7% for samples within 1 km, to < 3.5% for samples up to > 6000 km apart. Over a third of the species were not monophyletic, and the proportion increased through locally, nationally, regionally, and continentally restricted subsets of the data. The success of identifying queries decreased with increasing spatial scale of sampling; liberal methods declined from 100% to around 90%, whereas strict methods dropped to below 50% at continental scales. The proportion of query identifications considered uncertain (more than one species < 1% distance from query) escalated from zero at local, to 50% at continental scale. Finally, by resampling the most widely sampled species we show that even if samples are collected to maximize the geographical coverage, up to 70 individuals are required to sample 95% of intraspecific variation. The results show that the geographical scale of sampling has a critical impact on the global application of DNA barcoding. Scale-effects result from the relative importance of different processes determining the composition of regional species assemblages (dispersal and ecological assembly) and global clades (demography, speciation, and extinction). The incorporation of geographical information, where available, will be required to obtain identification rates at global scales equivalent to those in regional barcoding studies. Our result hence provides an impetus for both smarter barcoding tools and sprouting national barcoding initiatives—smaller geographical scales deliver higher accuracy.

Speciation and DNA barcodes: testing the effects of dispersal on the formation of discrete sequence clusters

Large-scale sequencing of short mtDNA fragments for biodiversity inventories (‘DNA barcoding’) indicates that sequence variation in animal mtDNA is highly structured and partitioned into discrete genetic clusters that correspond broadly to species-level entities. Here we explore how the migration rate, an important demographic parameter that is directly related to population isolation, might affect variation in the strength of mtDNA clustering among taxa. Patterns of mtDNA variation were investigated in two groups of beetles that both contain lineages occupying habitats predicted to select for different dispersal abilities: predacious diving beetles (Dytiscidae) in the genus Bidessus from lotic and lentic habitats across Europe and darkling beetles (Tenebrionidae) in the genus Eutagenia from sand and other soil types in the Aegean Islands. The degree of genetic clustering was determined using the recently developed ‘mixed Yule coalescent’ (MYC) model that detects the transition from between-species to within-population branching patterns. Lineages from presumed stable habitats, and therefore displaying lower dispersal ability and migration rates, showed greater levels of mtDNA clustering and geographical subdivision than their close relatives inhabiting ephemeral habitats. Simulations of expected patterns of mtDNA variation under island models showed that MYC clusters are only detected when the migration rates are much lower than the value of Nm=1 typically used to define the threshold for neutral genetic divergence. Therefore, discrete mtDNA clusters provide strong evidence for independently evolving populations or species, but their formation is suppressed even under very low levels of dispersal.

Bayesian Tests of Topology Hypotheses with an Example from Diving Beetles

We review Bayesian approaches to model testing in general and to the assessment of topological hypotheses in particular. We show that the standard way of setting up Bayes factor tests of the monophyly of a group, or the placement of a sample sequence in a known reference tree, can be misleading. The reason for this is related to the well-known dependency of Bayes factors on model-specific priors. Specifically, when testing tree hypotheses it is important that each hypothesis is associated with an appropriate tree space in the prior. This can be achieved by using appropriately constrained searches or by filtering trees in the posterior sample, but in a more elaborate way than typically implemented. If it is difficult to find the appropriate tree sets to be contrasted, then the posterior model odds may be more informative than the Bayes factor. We illustrate the recommended techniques using an empirical test case addressing the issue of whether two genera of diving beetles (Coleoptera: Dytiscidae), Suphrodytes and Hydroporus, should be synonymized. Our refined Bayes factor tests, in contrast to standard analyses, show that there is strong support for Suphrodytes nesting inside Hydroporus, and the genera are therefore synonymized. [Bayes factor; Coleoptera; Dytiscidae; marginal likelihood; model testing; posterior odds; reversible-jump MCMC; stepping-stone sampling.]

Nonrandom mating preserves intrasexual polymorphism and stops population differentiation in sexual conflict

Evolutionary conflict between the sexes is predicted to lead to sexual arms races in which male adaptations for acquiring mates (“offense” traits) are met by female counteradaptations—for example, to reduce mating rate (“defense” traits). Such coevolutionary chases may be perpetual. However, we show here that the coevolutionary process may also lead to a stable state in which multiple offense‐defense trait pairs are maintained. This type of polymorphism below the species level is a result of sexual conflict in combination with nonrandom mating. Our results show that if nonrandom mating occurs with respect to male and female conflict traits, genetic correlations will act to stabilize the trait frequencies so that all morphs are maintained. We discuss the results in special relation to the evolution of female polymorphism in diving beetles and argue that the process we describe may be a general force that maintains polymorphism in other taxa as well.

Whole-community DNA barcoding reveals a spatio-temporal continuum of biodiversity at species and genetic levels

A correlation of species and genetic diversity has been proposed but not uniformly supported. Large-scale DNA barcoding provides qualitatively novel data to test for correlations across hierarchical levels (genes, genealogies and species), and may help to unveil the underlying processes. Here we analyse sequence variation in communities of aquatic beetles across Europe (>5,000 individuals) to test for self-similarity of beta diversity patterns at multiple hierarchical levels. We show that community similarity at all levels decreases exponentially with geographic distance, and initial similarity is correlated with the lineage age, consistent with a molecular clock. Log-log correlations between lineage age, number of lineages, and range sizes, reveal a fractal geometry in time and space, indicating a spatio-temporal continuum of biodiversity across scales. Simulations show that these findings mirror dispersal-constrained models of haplotype distributions. These novel macroecological patterns may be explained by neutral evolutionary processes, acting continuously over time to produce multi-scale regularities of biodiversity.

Phylogeny of diving beetles reveals a coevolutionary arms race between the sexes.

Background Darwin illustrated his sexual selection theory with male and female morphology of diving beetles, but maintained a cooperative view of their interaction. Present theory suggests that instead sexual conflict should be a widespread evolutionary force driving both intersexual coevolutionary arms races and speciation. Methodology/Principal Findings We combined Bayesian phylogenetics, complete taxon sampling and a multi-gene approach to test the arms race scenario on a robust diving beetle phylogeny. As predicted, suction cups in males and modified dorsal surfaces in females showed a pronounced coevolutionary pattern. The female dorsal modifications impair the attachment ability of male suction cups, but each antagonistic novelty in females corresponds to counter-differentiation of suction cups in males. Conclusions A recently diverged sibling species pair in Japan is possibly one consequence of this arms race and we suggest that future studies on hypoxia might reveal the key to the extraordinary selection for female counter-adaptations in diving beetles.

Phylogeny and classification of the tribe Hydaticini (Coleoptera: Dytiscidae): partition choice for Bayesian analysis with multiple nuclear and mitochondrial protein‐coding genes

A phylogenetic analysis of the diving beetle tribe Hydaticini Sharp (Coleoptera: Dytiscidae: Dytiscinae) is presented based on data from adult morphology, two nuclear (histone III and wingless) and two mitochondrial (cytochrome c oxidase I and II) protein‐coding genes. We explore how to best partition a data set of multiple nuclear and mitochondrial protein‐coding genes by using Bayes factor and a penalized modification of Bayes Factor. Ten biologically relevant partitioning strategies were identified ranging from all DNA analysed under a single model to each codon position of each gene treated with a separate model. Model selection criteria AIC, AICc, BIC and four ways of traversing parameter space in a hierarchical likelihood ratio test were applied to each partition. All unique partitioning and model combinations were analysed with Bayesian methods. Results show that partitioning by codon position and genome source (nuclear vs. mitochondrial) is strongly favoured over partitioning by gene. We also find evidence that Bayes Factor can penalize overparameterization even when comparing nested models. Species groups showing a strong geographical pattern were generally highly supported, however, the sister group relationship of an isolated Madagascan and Australian species were shown to be artefactual with a long‐branch extraction test. The following conclusions were supported in both the selected method of partitioning the Bayesian analysis and combined parsimony analyses: (i) the tribe Hydaticini is monophyletic (ii) the genus Hydaticus Leach is paraphyletic with respect to Prodaticus Sharp (iii) the subgenus Hydaticus (Hydaticus) is monophyletic, and (iv) the subgenus H. (Guignotites) Brinck is paraphyletic with respect to Prodaticus and the subgenera H. (Pleurodytes) Régimbart and H. (Hydaticinus) Guignot. Based on these results, Hydaticus and Prodaticus are each recognized as valid genera and Guignotites, Hydaticinus and Pleurodytes are each placed as junior synonyms of Prodaticus (new synonymies).

Bayesian species delimitation reveals generalist and specialist parasitic wasps on Galerucella beetles (Chrysomelidae): sorting by herbivore or plant host

BackgroundTo understand the ecological and evolutionary consequences of species interactions in food webs necessitates that interactions are properly identified. Genetic analyses suggest that many supposedly generalist parasitoid species should rather be defined as multiple species with a more narrow diet, reducing the probability that such species may mediate indirect interactions such as apparent competition among hosts. Recent studies showed that the parasitoid Asecodes lucens mediate apparent competition between two hosts, Galerucella tenella and G. calmariensis, affecting both interaction strengths and evolutionary feedbacks. The same parasitoid was also recorded from other species in the genus Galerucella, suggesting that similar indirect effects may also occur for other species pairs.MethodsTo explore the possibility of such interactions, we sequenced mitochondrial and nuclear genetic markers to resolve the phylogeny of both host and parasitoid and to test the number of parasitoid species involved. We thus collected 139 Galerucella larvae from 8 host plant species and sequenced 31 adult beetle and 108 parasitoid individuals.ResultsThe analysis of the Galerucella data, that also included sequences from previous studies, verified the five species previously documented as reciprocally monophyletic, but the Bayesian species delimitation for A. lucens suggested 3–4 cryptic taxa with a more specialised host use than previously suggested. The gene data analyzed under the multispecies coalescent model allowed us to reconstruct the species tree phylogeny for both host and parasitoid and we found a fully congruent coevolutionary pattern suggesting that parasitoid speciation followed upon host speciation.ConclusionUsing multilocus sequence data in a Bayesian species delimitation analysis we propose that hymenopteran parasitoids of the genus Asecodes that infest Galerucella larvae constitute at least three species with narrow diet breath. The evolution of parasitoid Asecodes and host Galerucella show a fully congruent coevolutionary pattern. This finding strengthens the hypothesis that the parasitoid in host search uses cues of the host rather than more general cues of both host and plant.