- Letter
- Published:
Xenacoelomorpha is the sister group to Nephrozoa
- Johanna Taylor Cannon1,
- Bruno Cossermelli Vellutini2,
- Julian Smith3,
- Fredrik Ronquist1,
- Ulf Jondelius1 &
- …
- Andreas Hejnol2
Naturevolume 530, pages89–93 (2016)Cite this article
16kAccesses
373Citations
137Altmetric
Subjects
Abstract
The position of Xenacoelomorpha in the tree of life remains a major unresolved question in the study of deep animal relationships1. Xenacoelomorpha, comprising Acoela, Nemertodermatida, andXenoturbella, are bilaterally symmetrical marine worms that lack several features common to most other bilaterians, for example an anus, nephridia, and a circulatory system. Two conflicting hypotheses are under debate: Xenacoelomorpha is the sister group to all remaining Bilateria (= Nephrozoa, namely protostomes and deuterostomes)2,3 or is a clade inside Deuterostomia4. Thus, determining the phylogenetic position of this clade is pivotal for understanding the early evolution of bilaterian features, or as a case of drastic secondary loss of complexity. Here we show robust phylogenomic support for Xenacoelomorpha as the sister taxon of Nephrozoa. Our phylogenetic analyses, based on 11 novel xenacoelomorph transcriptomes and using different models of evolution under maximum likelihood and Bayesian inference analyses, strongly corroborate this result. Rigorous testing of 25 experimental data sets designed to exclude data partitions and taxa potentially prone to reconstruction biases indicates that long-branch attraction, saturation, and missing data do not influence these results. The sister group relationship between Nephrozoa and Xenacoelomorpha supported by our phylogenomic analyses implies that the last common ancestor of bilaterians was probably a benthic, ciliated acoelomate worm with a single opening into an epithelial gut, and that excretory organs, coelomic cavities, and nerve cords evolved after xenacoelomorphs separated from the stem lineage of Nephrozoa.
This is a preview of subscription content,access via your institution
Access options
Subscription info for Japanese customers
We have a dedicated website for our Japanese customers. Please go tonatureasia.com to subscribe to this journal.
Buy this article
- Purchase on SpringerLink
- Instant access to the full article PDF.
¥ 4,980
Prices may be subject to local taxes which are calculated during checkout




Similar content being viewed by others
Accession codes
References
Dunn, C. W., Giribet, G., Edgecombe, G. D. & Hejnol, A. Animal phylogeny and its evolutionary implications.Annu. Rev. Ecol. Evol. Syst.45, 371–395 (2014)
Hejnol, A. et al. Assessing the root of bilaterian animals with scalable phylogenomic methods.Proc. R. Soc. B276, 4261–4270 (2009)
Srivastava, M., Mazza-Curll, K. L., van Wolfswinkel, J. C. & Reddien, P. W. Whole-body acoel regeneration is controlled by Wnt and Bmp-Admp signaling.Curr. Biol.24, 1107–1113 (2014)
Philippe, H. et al. Acoelomorph flatworms are deuterostomes related toXenoturbella .Nature470, 255–258 (2011)
Nielsen, C.Animal Evolution: Interrelationships of the Living Phyla (Oxford Univ. Press, 2012)
Ehlers, U.Das phylogenetische System der Plathelminthes (G. Fischer, 1985)
Smith, J. P. S., III, Tyler, S. & Rieger, R. M. Is the Turbellaria polyphyletic?Hydrobiologia132, 13–21 (1986)
Haszprunar, G. Plathelminthes and Plathelminthomorpha — paraphyletic taxa.J. Zool. Syst. Evol. Res.34, 41–48 (1996)
Ruiz-Trillo, I., Riutort, M., Littlewood, D. T. J., Herniou, E. A. & Baguña, J. Acoel flatworms: earliest extant bilaterian metazoans, not members of Platyhelminthes.Science283, 1919–1923 (1999)
Steinböck, O. Ergebnisse einer von E. Reisinger & O. Steinböck mit Hilfe des Rask-Örsted fonds durchgefuhrten Reise in Grönland 1926. 2.Nemertoderma bathycola nov. gen. nov. spec., eine eigenartige Turbellarie aus der Tiefe der Diskobay: nebst einem Beitrag zur Kenntnis des Nemertinenepithels.Vidensk. Medd. Dan. Naturhist. Foren.90, 47–84 (1930)
Paps, J., Baguñà, J. & Riutort, M. Bilaterian phylogeny: a broad sampling of 13 nuclear genes provides a new Lophotrochozoa phylogeny and supports a paraphyletic basal acoelomorpha.Mol. Biol. Evol.26, 2397–2406 (2009)
Jondelius, U., Ruiz-Trillo, I., Baguñà, J. & Riutort, M. The Nemertodermatida are basal bilaterians and not members of the Platyhelminthes.Zool. Scr.31, 201–215 (2002)
Westblad, E.Xenoturbella bocki n.g, n.sp, a peculiar, primitive turbellarian type.Ark. Zool.1, 3–29 (1949)
Franzén, Å. & Afzelius, B. A. The ciliated epidermis ofXenoturbella bocki (Platyhelminthes, Xenoturbellida) with some phylogenetic considerations.Zool. Scr.16, 9–17 (1987)
Ehlers, U. & Sopott-Ehlers, B. Ultrastructure of the subepidermal musculature ofXenoturbella bocki, the adelphotaxon of the Bilateria.Zoomorphology117, 71–79 (1997)
Norén, M. & Jondelius, U.Xenoturbella’s molluscan relatives….Nature390, 31–32 (1997)
Bourlat, S. J. et al. Deuterostome phylogeny reveals monophyletic chordates and the new phylum Xenoturbellida.Nature444, 85–88 (2006)
Bourlat, S. J., Rota-Stabelli, O., Lanfear, R. & Telford, M. J. The mitochondrial genome structure ofXenoturbella bocki (phylum Xenoturbellida) is ancestral within the deuterostomes. BMCEvol. Biol.9, 107 (2009)
Mwinyi, A. et al. The phylogenetic position of Acoela as revealed by the complete mitochondrial genome ofSymsagittifera roscoffensis. BMCEvol. Biol.10, 309 (2010)
Ruiz-Trillo, I., Riutort, M., Fourcade, H. M., Baguñà, J. & Boore, J. L. Mitochondrial genome data support the basal position of Acoelomorpha and the polyphyly of the Platyhelminthes.Mol. Phylogenet. Evol.33, 321–332 (2004)
Rouse, G., Wilson, N. G., Carvajal, J. I. & Vrijenhoek, R. C. New deep-sea species ofXenoturbella and the position of Xenacoelomorpha.Naturehttp://dx.doi.org/10.1038/nature16545 (this issue)
Thomson, R. C., Plachetzki, D. C., Mahler, D. L. & Moore, B. R. A critical appraisal of the use of microRNA data in phylogenetics.Proc. Natl Acad. Sci. USA111, E3659–E3668 (2014)
Jondelius, U., Wallberg, A., Hooge, M. & Raikova, O. I. How the worm got its pharynx: phylogeny, classification and Bayesian assessment of character evolution in Acoela.Syst. Biol.60, 845–871 (2011)
Le, S. Q., Dang, C. C. & Gascuel, O. Modeling protein evolution with several amino acid replacement matrices depending on site rates.Mol. Biol. Evol.29, 2921–2936 (2012)
Dunn, C. W. et al. Broad phylogenomic sampling improves resolution of the animal tree of life.Nature452, 745–749 (2008)
Mirarab, S. et al. ASTRAL: genome-scale coalescent-based species tree estimation.Bioinformatics30, i541–i548 (2014)
Bleidorn, C. et al. On the phylogenetic position of Myzostomida: can 77 genes get it wrong? BMCEvol. Biol.9, 150 (2009)
Whelan, N. V., Kocot, K. M., Moroz, L. L. & Halanych, K. M. Error, signal, and the placement of Ctenophora sister to all other animals.Proc. Natl Acad. Sci. USA112, 5773–5778 (2015)
Hejnol, A. & Martindale, M. Q. Acoel development supports a simple planula-like urbilaterian.Phil. Trans. R. Soc. B363, 1493–1501 (2008)
Laumer, C. E. et al. Spiralian phylogeny informs the evolution of microscopic lineages.Curr. Biol.25, 2000–2006 (2015)
Grabherr, M. G. et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome.Nature Biotechnol.29, 644–652 (2011)
Ebersberger, I., Strauss, S. & von Haeseler, A. HaMStR: profile hidden markov model based search for orthologs in ESTs.BMC Evol. Biol.9, 157 (2009)
Lechner, M. et al. Proteinortho: detection of (co-)orthologs in large-scale analysis.BMC Bioinformatics12, 124 (2011)
Kocot, K. M. et al. Phylogenomics reveals deep molluscan relationships.Nature477, 452–456 (2011)
Cannon, J. T. et al. Phylogenomic resolution of the hemichordate and echinoderm clade.Curr. Biol.24, 2827–2832 (2014)
Katoh, K., Kuma, K., Toh, H. & Miyata, T. MAFFT version 5: improvement in accuracy of multiple sequence alignment.Nucleic Acids Res.33, 511–518 (2005)
Misof, B. & Misof, K. A Monte Carlo approach successfully identifies randomness in multiple sequence alignments: a more objective means of data exclusion.Syst. Biol.58, 21–34 (2009)
Price, M. N., Dehal, P. S. & Arkin, A. P . FastTree 2—approximately maximum-likelihood trees for large alignments.PLoS One 5, e9490 (2010)
Kocot, K. M., Citarella, M. R., Moroz, L. L. & Halanych, K. M. PhyloTreePruner: a phylogenetic tree-based approach for selection of orthologous sequences for phylogenomics.Evol. Bioinform. Online9, 429–435 (2013)
Rice, P., Longden, I. & Bleasby, A. EMBOSS: the European Molecular Biology Open Software Suite.Trends Genet.16, 276–277 (2000)
Struck, T. H. TreSpEx-Detection of misleading signal in phylogenetic reconstructions based on tree information.Evol. Bioinform. Online10, 51–67 (2014)
Criscuolo, A. & Gribaldo, S. BMGE (Block Mapping and Gathering with Entropy): a new software for selection of phylogenetic informative regions from multiple sequence alignments.BMC Evol. Biol.10, 210 (2010)
Stamatakis, A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies.Bioinformatics30, 1312–1313 (2014)
Darriba, D., Taboada, G. L., Doallo, R. & Posada, D. ProtTest 3: fast selection of best-fit models of protein evolution.Bioinformatics27, 1164–1165 (2011)
Lartillot, N., Rodrigue, N., Stubbs, D. & Richer, J. PhyloBayes MPI: phylogenetic reconstruction with infinite mixtures of profiles in a parallel environment.Syst. Biol.62, 611–615 (2013)
Ronquist, F. et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space.Syst. Biol.61, 539–542 (2012)
Acknowledgements
The Swedish Research Council provided funding for U.J. and J.T.C. (grant 2012-3913) and F.R. (grant 2014-5901). A.H. received support from the Sars Core budget and Marie Curie Innovative Training Networks ‘NEPTUNE’ (FP7-PEOPLE-2012-ITN 317172) and FP7-PEOPLE-2009-RG 256450. We thank N. Lartillot and K. Kocot for discussions. Hejnol laboratory members K. Pang and A. Børve assisted with RNA extraction; A. Boddington, J. Bengtsen and A. Elde assisted with culture forIsodiametra pulchra andConvolutriloba macropyga. Thanks to W. Sterrer for collection ofSterreria sp. andAscoparia sp., and to R. Janssen for findingX. bocki. The Sven Lovén Centre of Marine Sciences Kristineberg, University of Gothenburg, and the Interuniversity Institute of Marine Sciences in Eilat provided logistical support for field collection. S. Baldauf assisted with laboratory space and resources for complementary DNA synthesis. We thank K. Larsson for the original illustrations. Computations were performed on resources provided by the Swedish National Infrastructure for Computing (SNIC). Transcriptome assembly, data set construction, RAxML and PhyloBayes analyses were performed using resources provided through Uppsala Multidisciplinary Center for Advanced Computational Science (UPPMAX) under project b2013077, and MrBayes analyses were run under project snic2014-1-323.
Author information
Authors and Affiliations
Naturhistoriska Riksmuseet, PO Box 50007, Stockholm, SE-104 05, Sweden
Johanna Taylor Cannon, Fredrik Ronquist & Ulf Jondelius
Sars International Centre for Marine Molecular Biology, University of Bergen, Thormøhlensgate 55, Bergen, 5008, Norway
Bruno Cossermelli Vellutini & Andreas Hejnol
Department of Biology, Winthrop University, 701 Oakland Avenue, Rock Hill, 29733, South Carolina, USA
Julian Smith
- Johanna Taylor Cannon
Search author on:PubMed Google Scholar
- Bruno Cossermelli Vellutini
Search author on:PubMed Google Scholar
- Julian Smith
Search author on:PubMed Google Scholar
- Fredrik Ronquist
Search author on:PubMed Google Scholar
- Ulf Jondelius
Search author on:PubMed Google Scholar
- Andreas Hejnol
Search author on:PubMed Google Scholar
Contributions
J.T.C., U.J., B.C.V., and A.H. conceived and designed the study. U.J. and A.H. collected several specimens and J.S. III collectedDiopisthoporus gymnopharyngeus specimens. J.T.C. and B.C.V. performed molecular work and RNA sequencing assembly. J.T.C. assembled the datasets and performed phylogenetic analyses. F.R. conducted Bayesian phylogenetic analyses using MrBayes. All authors contributed to writing the manuscript.
Corresponding authors
Correspondence toJohanna Taylor Cannon orAndreas Hejnol.
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Additional information
Sequence data have been deposited in the NCBI Sequence Read Archive under BioProjectPRJNA295688. Data matrices and trees from this study are available from the Dryad Digital Repository (http://datadryad.org) under DOI10.5061/dryad.493b7.
Extended data figures and tables
Extended Data Figure 1 Maximum likelihood topology of metazoan relationships inferred from 336 genes from the best-sampled 56 taxa.
Maximum likelihood tree is shown as inferred using the LG + I + Γ model for each gene partition, and 100 bootstrap replicates. Filled blue circles represent 100% bootstrap support. The length of the matrix is 81,451 amino acids and overall matrix completeness is 89%.
Extended Data Figure 2 Maximum likelihood topology of metazoan relationships inferred from 881 genes and 77 taxa.
Maximum likelihood tree is shown as inferred using the LG + I + Γ model for each gene partition, and 100 bootstrap replicates. Filled blue circles represent 100% bootstrap support. The length of the matrix is 337,954 amino acids and overall matrix completeness is 62%.
Extended Data Figure 3 Maximum likelihood topology of metazoan relationships inferred from 212 genes with Acoelomorpha removed.
Maximum likelihood tree is shown as inferred using the LG + I + Γ model for each gene partition, and 100 bootstrap replicates. Filled blue circles represent 100% bootstrap support. The length of the matrix is 43,942 amino acids and overall matrix completeness is 70%.
Extended Data Figure 4 ASTRAL species tree, constructed from 212 input partial gene trees inferred in RAxML version 8.0.20.
Nodal support values reflect the frequency of splits in trees constructed by ASTRAL from 100 bootstrap replicate gene trees.
Extended Data Figure 5 Bayesian inference topology of metazoan relationships inferred on the basis of 212 genes and 78 taxa.
Results are shown from MrBayes analyses of four independent Metropolis-coupled chains run for 4,000,000 generations, with sampling every 500 generations. Amino-acid data were back-translated to nucleotides and analysed under an independent substitution model.
Supplementary information
Supplementary Information
This file contains a Supplementary Discussion, Supplementary Tables 1-3 and Supplementary Figures 1-22. (PDF 5072 kb)
Rights and permissions
About this article
Cite this article
Cannon, J., Vellutini, B., Smith, J.et al. Xenacoelomorpha is the sister group to Nephrozoa.Nature530, 89–93 (2016). https://doi.org/10.1038/nature16520
Received:
Accepted:
Published:
Issue date:
Share this article
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
This article is cited by
Structure of putative epidermal sensory receptors in an acoel flatworm, Praesagittifera naikaiensis
- Tosuke Sakagami
- Kaho Watanabe
- Motonori Ando
Cell and Tissue Research (2024)
Highly conserved and extremely evolvable: BMP signalling in secondary axis patterning of Cnidaria and Bilateria
- David Mörsdorf
- Paul Knabl
- Grigory Genikhovich
Development Genes and Evolution (2024)
Induced spawning with gamete release from body ruptures during reproduction of Xenoturbella bocki
- Hiroaki Nakano
- Ako Nakano
- Michael C. Thorndyke
Communications Biology (2023)
Emergence of distinct syntenic density regimes is associated with early metazoan genomic transitions
- Nicolas S. M. Robert
- Fatih Sarigol
- Oleg Simakov
BMC Genomics (2022)


