Xenoturbellida and Acoelomorpha are marine worms with contentious ancestry. Both were originally associated with the flatworms (Platyhelminthes), but molecular data haverevised their phylogenetic positions, generally linking Xenoturbellida to the deuterostomes1,2 and positioning the Acoelomorpha as the most basally branching bilaterian group(s)3–6. Recent phylogenomic data suggested that Xenoturbellida and Acoelomorpha are sister taxa and together constitute an early branch of Bilateria7. Here we assemble three independent data sets—mitochondrial genes, a phylogenomic data set of 38,330 amino-acid positions and new microRNA (miRNA) complements—and show that the position of Acoelomorpha is strongly affected by a long-branch attraction (LBA) artefact. When we minimize LBA we find consistent support for a position of both acoelomorphs and Xenoturbella within the deuterostomes. The most likely phylogeny links Xenoturbella and Acoelomorpha in a clade we call Xenacoelomorpha. The Xenacoelomorpha is the sister group of the Ambulacraria (hemichordates and echinoderms). We show that analyses of miRNA complements8 have been affected by character loss in the acoels and that both groups possess one miRNA and the gene Rsb66 otherwise specific to deuterostomes. In addition, Xenoturbella shares one miRNA with the ambulacrarians, and two with the acoels. This phylogeny makes sense of the shared characteristics of Xenoturbellida and Acoelomorpha, such as ciliary ultrastructure and diffuse nervous system, and implies the loss of various deuterostome characters in the Xenacoelomorpha including coelomic cavities, through gut and gill slits.

No abstract

No abstract

Meiotic recombination is a fundamental cellular process, with important consequences for evolution and genome integrity. However, we know little about how recombination rates vary across the genomes of most species and the molecular and evolutionary determinants of this variation. The honeybee, Apis mellifera, has extremely high rates of meiotic recombination, although the evolutionary causes and consequences of this are unclear. Here we use patterns of linkage disequilibrium in whole genome resequencing data from 30 diploid honeybees to construct a fine-scale map of rates of crossing over in the genome. We find that, in contrast to vertebrate genomes, the recombination landscape is not strongly punctate. Crossover rates strongly correlate with levels of genetic variation, but not divergence, which indicates a pervasive impact of selection on the genome. Germ-line methylated genes have reduced crossover rate, which could indicate a role of methylation in suppressing recombination. Controlling for the effects of methylation, we do not infer a strong association between gene expression patterns and recombination. The site frequency spectrum is strongly skewed from neutral expectations in honeybees: rare variants are dominated by AT-biased mutations, whereas GC-biased mutations are found at higher frequencies, indicative of a major influence of GC-biased gene conversion (gBGC), which we infer to generate an allele fixation bias 5 – 50 times the genomic average estimated in humans. We uncover further evidence that this repair bias specifically affects transitions and favours fixation of CpG sites. Recombination, via gBGC, therefore appears to have profound consequences on genome evolution in honeybees and interferes with the process of natural selection. These findings have important implications for our understanding of the forces driving molecular evolution.

Background The ability to generate long sequencing reads and access long-range linkage information is revolutionizing the quality and completeness of genome assemblies. Here we use a hybrid approach that combines data from four genome sequencing and mapping technologies to generate a new genome assembly of the honeybee Apis mellifera . We first generated contigs based on PacBio sequencing libraries, which were then merged with linked-read 10x Chromium data followed by scaffolding using a BioNano optical genome map and a Hi-C chromatin interaction map, complemented by a genetic linkage map. Results Each of the assembly steps reduced the number of gaps and incorporated a substantial amount of additional sequence into scaffolds. The new assembly (Amel_HAv3) is significantly more contiguous and complete than the previous one (Amel_4.5), based mainly on Sanger sequencing reads. N50 of contigs is 120-fold higher (5.381 Mbp compared to 0.053 Mbp) and we anchor > 98% of the sequence to chromosomes. All of the 16 chromosomes are represented as single scaffolds with an average of three sequence gaps per chromosome. The improvements are largely due to the inclusion of repetitive sequence that was unplaced in previous assemblies. In particular, our assembly is highly contiguous across centromeres and telomeres and includes hundreds of AvaI and AluI repeats associated with these features. Conclusions The improved assembly will be of utility for refining gene models, studying genome function, mapping functional genetic variation, identification of structural variants, and comparative genomics. Electronic supplementary material The online version of this article (10.1186/s12864-019-5642-0) contains supplementary material, which is available to authorized users.

The native range of the honeybee Apis mellifera encompasses Europe, Africa, and the Middle East, whereas the nine other species of Apis are found exclusively in Asia. It is therefore commonly assumed that A. mellifera arose in Asia and expanded into Europe and Africa. However, other hypotheses for the origin of A. mellifera have also been proposed based on phylogenetic trees constructed from genetic markers. In particular, an analysis based on >1000 single‐nucleotide polymorphism markers placed the root of the tree of A. mellifera subspecies among samples from Africa, suggestive of an out‐of‐Africa expansion. Here, we re‐evaluate the evidence for this and other hypotheses by testing the robustness of the tree topology to different tree‐building methods and by removing specimens with a potentially hybrid background. These analyses do not unequivocally place the root of the tree of A. mellifera subspecies within Africa, and are potentially consistent with a variety of hypotheses for honeybee evolution, including an expansion out of Asia. Our analyses also support high divergence between western and eastern European populations of A. mellifera, suggesting they are likely derived from two distinct colonization routes, although the sources of these expansions are still unclear.

Acoela are marine microscopic worms currently thought to be the sister taxon of all other bilaterians. Acoels have long been used as models in evolutionary scenarios, and generalized conclusions about acoel and bilaterian ancestral features are frequently drawn from studies of single acoel species. There is no extensive phylogenetic study of Acoela and the taxonomy of the 380 species is chaotic. Here we use two nuclear ribosomal genes and one mitochondrial gene in combination with 37 morphological characters in an analysis of 126 acoel terminals (about one-third of the described species) to estimate the phylogeny and character evolution of Acoela. We present an estimate of posterior probabilities for ancestral character states at 31 control nodes in the phylogeny. The overall reconstruction signal based on the shape of the posterior distribution of character states was computed for all morphological characters and control nodes to assess how well these were reconstructed. The body-wall musculature appears more clearly reconstructed than the reproductive organs. Posterior similarity to the root was calculated by averaging the divergence between the posterior distributions at the nodes and the root over all morphological characters. Diopisthoporidae is the sister group to all other acoels and has the highest posterior similarity to the root. Convolutidae, including several "model" acoels, is most divergent. Finally, we present a phylogenetic classification of Acoela down to the family level where six previous family level taxa are synonymized.

The transition to a vermiform body shape is one of the most important events in animal evolution, having led to the impressive radiation of Bilateria. However, the sister group of Bilateria has remained obscure. Cladistic analyses of morphology indicate that Ctenophora is the sister group of Bilateria. Previous analyses of SSU rRNA sequences have yielded conflicting results; in many studies Ctenophora forms the sister group of Cnidaria + Bilateria, but in others the ctenophores group with poriferans. Here we re-examine the SSU sequence by analyzing a dataset with 528 metazoan + outgroup sequences, including almost 120 poriferan and diploblast sequences. We use parsimony ratchet and jackknife methods, as well as Bayesian methods, to analyze the data. The results indicate strong phylogenetic signals for a cnidarian + bilaterian group and for the comb jellies to have branched off early within a group uniting all epithelial animals [(Ct,(Cn,Bi))]. We demonstrate the importance of inclusive taxonomic coverage of ribosomal sequences for resolving this problematic part of the metazoan tree: topological stability increases dramatically with the addition of taxa, and the jackknife frequencies of the internal nodes uniting the lineages [(Cn,Bi) and ((Ct,(Cn,Bi))] also increase. We consider the reconstructed topology to represent the current best hypothesis of the interrelationships of these old lineages. Some morphological features supporting alternative hypotheses are discussed in the light of this result.