Mitigating Anticipated Effects of Systematic Errors Supports Sister-Group Relationship between Xenacoelomorpha and Ambulacraria

A Phylogenomic Backbone for Acoelomorpha Inferred From Transcriptomic Data

Xenacoelomorpha are mostly microscopic, morphologically simple worms, lacking many structures typical of other bilaterians. Xenacoelomorphs-which include three main groups, namely Acoela, Nemertodermatida, and Xenoturbella-have been proposed to be an early diverging Bilateria, sister to protostomes and deuterostomes, but other phylogenomic analyses have recovered this clade nested within the deuterostomes, as sister to Ambulacraria. The position of Xenacoelomorpha within the metazoan tree has understandably attracted a lot of attention, overshadowing the study of phylogenetic relationships within this group. Given that Xenoturbella includes only six species whose relationships are well understood, we decided to focus on the most speciose Acoelomorpha (Acoela + Nemertodermatida). Here, we have sequenced 29 transcriptomes, doubling the number of sequenced species, to infer a backbone tree for Acoelomorpha based on genomic data. The recovered topology is mostly congruent with previous studies. The most important difference is the recovery of Paratomella as the first off-shoot within Acoela, dramatically changing the reconstruction of the ancestral acoel. Besides, we have detected incongruence between the gene trees and the species tree, likely linked to incomplete lineage sorting, and some signal of introgression between the families Dakuidae and Mecynostomidae, which hampers inferring the correct placement of this family and, particularly, of the genus Notocelis. We have also used this dataset to infer for the first time diversification times within Acoelomorpha, which coincide with known bilaterian diversification and extinction events. Given the importance of morphological data in acoelomorph phylogenetics, we tested several partitions and models. Although morphological data failed to recover a robust phylogeny, phylogenetic placement has proven to be a suitable alternative when a reference phylogeny is available.

Ediacaran origin and Ediacaran-Cambrian diversification of Metazoa

The timescale of animal diversification has been a focus of debate over how evolutionary history should be calibrated to geologic time. Molecular clock analyses have invariably estimated a Cryogenian or Tonian origin of animals while unequivocal animal fossils first occur in the Ediacaran. However, redating of key Ediacaran biotas and the discovery of several Ediacaran crown-Metazoa prompt recalibration of molecular clock analyses. We present revised fossil calibrations and use them in molecular clock analyses estimating the timescale of metazoan evolutionary history. Integrating across uncertainties including phylogenetic relationships, clock model, and calibration strategy, we estimate Metazoa to have originated in the early Ediacaran, Eumetazoa in the middle Ediacaran, and Bilateria in the upper Ediacaran, with many crown-phyla originating across the Ediacaran-Cambrian interval or elsewise fully within the Cambrian. These results are in much closer accord with the fossil record, coinciding with marine oxygenation, but they reject a literal reading of the fossil record.

Evolution of the Cdk4/6-Cdkn2 system in invertebrates

The cell cycle is driven by cyclin-dependent kinases (Cdks). The decision whether the cell cycle proceeds is made during G1 phase, when Cdk4/6 functions. Cyclin-dependent kinase inhibitor 2 (Cdkn2) is a specific inhibitor of Cdk4/6, and their interaction depends on D84 in Cdkn2 and R24/31 in Cdk4/6. This knowledge is based mainly on studies in mammalian cells. Here, we comprehensively analyzed Cdk4/6 and Cdkn2 in invertebrates and found that Cdk4/6 was present in most of the investigated phyla, but the distribution of Cdkn2 was rather uneven among and within the phyla. The positive charge of R24/R31 in Cdk4/6 was conserved in all analyzed species in phyla with Cdkn2. The presence of Cdkn2 and the conservation of the positive charge were statistically correlated. We also found that Cdkn2 has been tightly linked to Fas associated factor 1 (Faf1) during evolution. We discuss potential interactions between Cdkn2 and Cdk4/6 in evolution and the possible cause of the strong conservation of the microsynteny.

Acoelomorph flatworm monophyly is a long-branch attraction artefact obscuring a clade of Acoela and Xenoturbellida

Acoelomorpha is a broadly accepted clade of bilaterian animals made up of the fast-evolving, morphologically simple, mainly marine flatworm lineages Acoela and Nemertodermatida. Phylogenomic studies support Acoelomorpha's close relationship with the slowly evolving and similarly simplistic Xenoturbella, together forming the phylum Xenacoelomorpha. The phylogenetic placement of Xenacoelomorpha amongst bilaterians is controversial, with some studies supporting Xenacoelomorpha as the sister group to all other bilaterians, implying that their simplicity may be representative of early bilaterians. Others propose that this placement is an error resulting from the fast-evolving Acoelomorpha, and instead suggest that they are the degenerate sister group to Ambulacraria. Perhaps as a result of this debate, internal xenacoelomorph relationships have been somewhat overlooked at a phylogenomic scale. Here, I employ a highly targeted approach to detect and overcome possible phylogenomic error in the relationship between Xenoturbella and the fast-evolving acoelomorph flatworms. The results indicate that the subphylum Acoelomorpha is a long-branch attraction artefact obscuring a previously undiscovered clade comprising Xenoturbella and Acoela, which I name Xenacoela. The findings also suggest that Xenacoelomorpha is not the sister group to all other bilaterians. This study provides a template for future efforts aimed at discovering and correcting unrecognized long-branch attraction artefacts throughout the tree of life.

A taxon-rich and genome-scale phylogeny of Opisthokonta

Ancient divergences within Opisthokonta-a major lineage that includes organisms in the kingdoms Animalia, Fungi, and their unicellular relatives-remain contentious. To assess progress toward a genome-scale Opisthokonta phylogeny, we conducted the most taxon rich phylogenomic analysis using sets of genes inferred with different orthology inference methods and established the geological timeline of Opisthokonta diversification. We also conducted sensitivity analysis by subsampling genes or taxa from the full data matrix based on filtering criteria previously shown to improve phylogenomic inference. We found that approximately 85% of internal branches were congruent across data matrices and the approaches used. Notably, the use of different orthology inference methods was a substantial contributor to the observed incongruence: analyses using the same set of orthologs showed high congruence of 97% to 98%, whereas different sets of orthologs resulted in somewhat lower congruence (87% to 91%). Examination of unicellular Holozoa relationships suggests that the instability observed across varying gene sets may stem from weak phylogenetic signals. Our results provide a comprehensive Opisthokonta phylogenomic framework that will be useful for illuminating ancient evolutionary episodes concerning the origin and diversification of the 2 major eukaryotic kingdoms and emphasize the importance of investigating effects of orthology inference on phylogenetic analyses to resolve ancient divergences.

Insights into early animal evolution from the genome of the xenacoelomorph worm Xenoturbella bocki

The evolutionary origins of Bilateria remain enigmatic. One of the more enduring proposals highlights similarities between a cnidarian-like planula larva and simple acoel-like flatworms. This idea is based in part on the view of the Xenacoelomorpha as an outgroup to all other bilaterians which are themselves designated the Nephrozoa (protostomes and deuterostomes). Genome data can provide important comparative data and help understand the evolution and biology of enigmatic species better. Here, we assemble and analyze the genome of the simple, marine xenacoelomorph Xenoturbella bocki, a key species for our understanding of early bilaterian evolution. Our highly contiguous genome assembly of X. bocki has a size of ~111 Mbp in 18 chromosome-like scaffolds, with repeat content and intron, exon, and intergenic space comparable to other bilaterian invertebrates. We find X. bocki to have a similar number of genes to other bilaterians and to have retained ancestral metazoan synteny. Key bilaterian signaling pathways are also largely complete and most bilaterian miRNAs are present. Overall, we conclude that X. bocki has a complex genome typical of bilaterians, which does not reflect the apparent simplicity of its body plan that has been so important to proposals that the Xenacoelomorpha are the simple sister group of the rest of the Bilateria.

The brittle star genome illuminates the genetic basis of animal appendage regeneration

Species within nearly all extant animal lineages are capable of regenerating body parts. However, it remains unclear whether the gene expression programme controlling regeneration is evolutionarily conserved. Brittle stars are a species-rich class of echinoderms with outstanding regenerative abilities, but investigations into the genetic bases of regeneration in this group have been hindered by the limited genomic resources. Here we report a chromosome-scale genome assembly for the brittle star Amphiura filiformis. We show that the brittle star genome is the most rearranged among echinoderms sequenced so far, featuring a reorganized Hox cluster reminiscent of the rearrangements observed in sea urchins. In addition, we performed an extensive profiling of gene expression during brittle star adult arm regeneration and identified sequential waves of gene expression governing wound healing, proliferation and differentiation. We conducted comparative transcriptomic analyses with other invertebrate and vertebrate models for appendage regeneration and uncovered hundreds of genes with conserved expression dynamics, particularly during the proliferative phase of regeneration. Our findings emphasize the crucial importance of echinoderms to detect long-range expression conservation between vertebrates and classical invertebrate regeneration model systems.

A wound-induced differentiation trajectory for neurons

Animals capable of whole-body regeneration can replace any missing cell type and regenerate fully functional new organs, including new brains, de novo. The regeneration of a new brain requires the formation of diverse neural cell types and their assembly into an organized structure with correctly wired circuits. Recent work in various regenerative animals has revealed transcriptional programs required for the differentiation of distinct neural subpopulations, however, how these transcriptional programs are initiated in response to injury remains unknown. Here, we focused on the highly regenerative acoel worm, Hofstenia miamia, to study wound-induced transcriptional regulatory events that lead to the production of neurons and subsequently a functional brain. Footprinting analysis using chromatin accessibility data on a chromosome-scale genome assembly revealed that binding sites for the Nuclear Factor Y (NFY) transcription factor complex were significantly bound during regeneration, showing a dynamic increase in binding within one hour upon amputation specifically in tail fragments, which will regenerate a new brain. Strikingly, NFY targets were highly enriched for genes with neuronal function. Single-cell transcriptome analysis combined with functional studies identified soxC+ stem cells as a putative progenitor population for multiple neural subtypes. Further, we found that wound-induced soxC expression is likely under direct transcriptional control by NFY, uncovering a mechanism for the initiation of a neural differentiation pathway by early wound-induced binding of a transcriptional regulator.

Chromosome-level genome assemblies of 2 hemichordates provide new insights into deuterostome origin and chromosome evolution

Deuterostomes are a monophyletic group of animals that includes Hemichordata, Echinodermata (together called Ambulacraria), and Chordata. The diversity of deuterostome body plans has made it challenging to reconstruct their ancestral condition and to decipher the genetic changes that drove the diversification of deuterostome lineages. Here, we generate chromosome-level genome assemblies of 2 hemichordate species, Ptychodera flava and Schizocardium californicum, and use comparative genomic approaches to infer the chromosomal architecture of the deuterostome common ancestor and delineate lineage-specific chromosomal modifications. We show that hemichordate chromosomes (1N = 23) exhibit remarkable chromosome-scale macrosynteny when compared to other deuterostomes and can be derived from 24 deuterostome ancestral linkage groups (ALGs). These deuterostome ALGs in turn match previously inferred bilaterian ALGs, consistent with a relatively short transition from the last common bilaterian ancestor to the origin of deuterostomes. Based on this deuterostome ALG complement, we deduced chromosomal rearrangement events that occurred in different lineages. For example, a fusion-with-mixing event produced an Ambulacraria-specific ALG that subsequently split into 2 chromosomes in extant hemichordates, while this homologous ALG further fused with another chromosome in sea urchins. Orthologous genes distributed in these rearranged chromosomes are enriched for functions in various developmental processes. We found that the deeply conserved Hox clusters are located in highly rearranged chromosomes and that maintenance of the clusters are likely due to lower densities of transposable elements within the clusters. We also provide evidence that the deuterostome-specific pharyngeal gene cluster was established via the combination of 3 pre-assembled microsyntenic blocks. We suggest that since chromosomal rearrangement events and formation of new gene clusters may change the regulatory controls of developmental genes, these events may have contributed to the evolution of diverse body plans among deuterostomes.

Single cell atlas of Xenoturbella bocki highlights limited cell-type complexity

Phylogenetic analyses over the last two decades have united a few small, and previously orphan clades, the nematodermatids, acoels and xenoturbelids, into the phylum Xenacoelomorpha. Some phylogenetic analyses support a sister relationship between Xenacoelomorpha and Ambulacraria (Xenambulacraria), while others suggest that Xenacoelomorpha may be sister to the rest of the Bilateria (Nephrozoa). An understanding of the cell type complements of Xenacoelomorphs is essential to assessing these alternatives as well as to our broader understanding of bilaterian cell type evolution. Employing whole organism single-cell RNA-seq in the marine xenacoelomorph worm Xenoturbella bocki, we show that Xenambulacrarian nerve nets share regulatory features and a peptidergic identity with those found in cnidarians and protostomes and more broadly share muscle and gland cell similarities with other metazoans. Taken together, these data are consistent with broad homologies of animal gland, muscle, and neurons as well as more specific affinities between Xenoturbella and acoel gut and epidermal tissues, consistent with the monophyly of Xenacoelomorpha.

Structure of putative epidermal sensory receptors in an acoel flatworm, Praesagittifera naikaiensis

Acoel flatworms possess epidermal sensory-receptor cells on their body surfaces and exhibit behavioral repertoires such as geotaxis and phototaxis. Acoel epidermal sensory receptors should be mechanical and/or chemical receptors; however, the mechanisms of their sensory reception have not been elucidated. We examined the three-dimensional relationship between epidermal sensory receptors and their innervation in an acoel flatworm, Praesagittifera naikaiensis. The distribution of the sensory receptors was different between the ventral and dorsal sides of worms. The nervous system was mainly composed of a peripheral nerve net, an anterior brain, and three pairs of longitudinal nerve cords. The nerve net was located closer to the body surface than the brain and the nerve cords. The sensory receptors have neural connections with the nerve net in the entire body of worms. We identified five homologs of polycystic kidney disease (PKD): PKD1-1, PKD1-2, PKD1-3, PKD1-4, and, PKD2, from the P. naikaiensis genome. All of these PKD genes were implied to be expressed in the epidermal sensory receptors of P. naikaiensis. PKD1-1 and PKD2 were dispersed across the entire body of worms. PKD1-2, PKD1-3, and PKD1-4 were expressed in the anterior region of worms. PKD1-4 was also expressed around the mouth opening. Our results indicated that P. naikaiensis possessed several types of epidermal sensory receptors to convert various environmental stimuli into electrical signals via the PKD channels and transmit the signals to afferent nerve and/or effector cells.

Evolution of the Ikaros family transcription factors: From a deuterostome ancestor to humans

Ikaros family proteins (Ikaros, Helios, Aiolos, Eos) are zinc finger transcription factors essential for the development and function of the adaptive immune system. They also control developmental events in neurons and other cell types, suggesting that they possess crucial functions across disparate cell types. These functions are likely shared among the organisms in which these factors exist, and it is thus important to obtain a view of their distribution and conservation across organisms. How this family evolved remains poorly understood. Here we mined protein, mRNA and DNA databases to identify proteins with DNA-binding domains homologous to that of Ikaros. We show that Ikaros-related proteins exist in organisms from all four deuterostome phyla (chordates, echinoderms, hemichordates, xenacoelomorpha), but not in more distant groups. While most non-vertebrates have a single family member, this family grew to six members in the acoel worm Hofstenia miamia, three in jawless and four in jawed vertebrates. Most residues involved in DNA contact from zinc fingers 2 to 4 were identical across the Ikaros family, suggesting conserved mechanisms for target sequence recognition. Further, we identified a novel KRKxxxPxK/R motif that inhibits DNA binding in vitro which was conserved across the deuterostome phyla. We also identified a EψψxxxψM(D/E)QAIxxAIxYLGA(D/E)xL motif conserved among human Ikaros, Aiolos, Helios and subsets of chordate proteins, and motifs that are specific to subsets of vertebrate family members. Some of these motifs are targets of mutations in human patients. Finally we show that the atypical family member Pegasus emerged only in vertebrates, which is consistent with its function in bone. Our data provide a novel evolutionary perspective for Ikaros family proteins and suggest that they have conserved regulatory functions across deuterostomes.

DNA Sequences Are as Useful as Protein Sequences for Inferring Deep Phylogenies

Inference of deep phylogenies has almost exclusively used protein rather than DNA sequences based on the perception that protein sequences are less prone to homoplasy and saturation or to issues of compositional heterogeneity than DNA sequences. Here, we analyze a model of codon evolution under an idealized genetic code and demonstrate that those perceptions may be misconceptions. We conduct a simulation study to assess the utility of protein versus DNA sequences for inferring deep phylogenies, with protein-coding data generated under models of heterogeneous substitution processes across sites in the sequence and among lineages on the tree, and then analyzed using nucleotide, amino acid, and codon models. Analysis of DNA sequences under nucleotide-substitution models (possibly with the third codon positions excluded) recovered the correct tree at least as often as analysis of the corresponding protein sequences under modern amino acid models. We also applied the different data-analysis strategies to an empirical dataset to infer the metazoan phylogeny. Our results from both simulated and real data suggest that DNA sequences may be as useful as proteins for inferring deep phylogenies and should not be excluded from such analyses. Analysis of DNA data under nucleotide models has a major computational advantage over protein-data analysis, potentially making it feasible to use advanced models that account for among-site and among-lineage heterogeneity in the nucleotide-substitution process in inference of deep phylogenies.

The draft genome of the microscopic Nemertoderma westbladi sheds light on the evolution of Acoelomorpha genomes

Background: Xenacoelomorpha is a marine clade of microscopic worms that is an important model system for understanding the evolution of key bilaterian novelties, such as the excretory system. Nevertheless, Xenacoelomorpha genomics has been restricted to a few species that either can be cultured in the lab or are centimetres long. Thus far, no genomes are available for Nemertodermatida, one of the group's main clades and whose origin has been dated more than 400 million years ago. Methods: DNA was extracted from a single specimen and sequenced with HiFi following the PacBio Ultra-Low DNA Input protocol. After genome assembly, decontamination, and annotation, the genome quality was benchmarked using two acoel genomes and one Illumina genome as reference. The gene content of three cnidarians, three acoelomorphs, four deuterostomes, and eight protostomes was clustered in orthogroups to make inferences of gene content evolution. Finally, we focused on the genes related to the ultrafiltration excretory system to compare patterns of presence/absence and gene architecture among these clades. Results: We present the first nemertodermatid genome sequenced from a single specimen of Nemertoderma westbladi. Although genome contiguity remains challenging (N50: 60 kb), it is very complete (BUSCO: 80.2%, Metazoa; 88.6%, Eukaryota) and the quality of the annotation allows fine-detail analyses of genome evolution. Acoelomorph genomes seem to be relatively conserved in terms of the percentage of repeats, number of genes, number of exons per gene and intron size. In addition, a high fraction of genes present in both protostomes and deuterostomes are absent in Acoelomorpha. Interestingly, we show that all genes related to the excretory system are present in Xenacoelomorpha except Osr, a key element in the development of these organs and whose acquisition seems to be interconnected with the origin of the specialised excretory system. Conclusion: Overall, these analyses highlight the potential of the Ultra-Low Input DNA protocol and HiFi to generate high-quality genomes from single animals, even for relatively large genomes, making it a feasible option for sequencing challenging taxa, which will be an exciting resource for comparative genomics analyses.

Confusion will be my epitaph: genome-scale discordance stifles phylogenetic resolution of Holothuroidea

Sea cucumbers (Holothuroidea) are a diverse clade of echinoderms found from intertidal waters to the bottom of the deepest oceanic trenches. Their reduced skeletons and limited number of phylogenetically informative traits have long obfuscated morphological classifications. Sanger-sequenced molecular datasets have also failed to constrain the position of major lineages. Noteworthy, topological uncertainty has hindered a resolution for Neoholothuriida, a highly diverse clade of Permo-Triassic age. We perform the first phylogenomic analysis of Holothuroidea, combining existing datasets with 13 novel transcriptomes. Using a highly curated dataset of 1100 orthologues, our efforts recapitulate previous results, struggling to resolve interrelationships among neoholothuriid clades. Three approaches to phylogenetic reconstruction (concatenation under both site-homogeneous and site-heterogeneous models, and coalescent-aware inference) result in alternative resolutions, all of which are recovered with strong support and across a range of datasets filtered for phylogenetic usefulness. We explore this intriguing result using gene-wise log-likelihood scores and attempt to correlate these with a large set of gene properties. While presenting novel ways of exploring and visualizing support for alternative trees, we are unable to discover significant predictors of topological preference, and our efforts fail to favour one topology. Neoholothuriid genomes seem to retain an amalgam of signals derived from multiple phylogenetic histories.

Independent rediploidization masks shared whole genome duplication in the sturgeon-paddlefish ancestor

Whole genome duplication (WGD) is a dramatic evolutionary event generating many new genes and which may play a role in survival through mass extinctions. Paddlefish and sturgeon are sister lineages that both show genomic evidence for ancient WGD. Until now this has been interpreted as two independent WGD events due to a preponderance of duplicate genes with independent histories. Here we show that although there is indeed a plurality of apparently independent gene duplications, these derive from a shared genome duplication event occurring well over 200 million years ago, likely close to the Permian-Triassic mass extinction period. This was followed by a prolonged process of reversion to stable diploid inheritance (rediploidization), that may have promoted survival during the Triassic-Jurassic mass extinction. We show that the sharing of this WGD is masked by the fact that paddlefish and sturgeon lineage divergence occurred before rediploidization had proceeded even half-way. Thus, for most genes the resolution to diploidy was lineage-specific. Because genes are only truly duplicated once diploid inheritance is established, the paddlefish and sturgeon genomes are thus a mosaic of shared and non-shared gene duplications resulting from a shared genome duplication event.

A wound-induced differentiation trajectory for neurons

Animals capable of whole-body regeneration can replace any missing cell type and regenerate fully-functional new organs, de novo . The regeneration of a new brain requires the formation of diverse neuronal cell types and their assembly into an organized structure and correctly-wired circuits. Recent work in various regenerative animals has revealed transcriptional programs required for the differentiation of distinct neuronal subpopulations, however how these transcriptional programs are initiated upon amputation remains unknown. Here, we focused on the highly regenerative acoel worm, Hofstenia miamia , to study wound-induced transcriptional regulatory events that lead to the production of neurons. Footprinting analysis using chromatin accessibility data on an improved genome assembly revealed that binding sites for the NFY transcription factor complex were significantly bound during regeneration, showing a dynamic increase in binding within one hour upon amputation specifically in tail fragments, which will regenerate a new brain. Strikingly, NFY targets were highly enriched for genes with neuronal functional. Single-cell transcriptome analysis combined with functional studies identified sox4 + stem cells as the likely progenitor population for multiple neuronal subtypes. Further, we found that wound-induced sox4 expression is likely under direct transcriptional control by NFY, uncovering a mechanism for how early wound-induced binding of a transcriptional regulator results in the initiation of a neuronal differentiation pathway.

Highlights

A new chromosome-scale assembly for Hofstenia enables comprehensive analysis of transcription factor binding during regeneration NFY motifs become dynamically bound by 1hpa in regenerating tail fragments, particularly in the loci of neural genes A sox4 + neural-specialized stem cell is identified using scRNA-seq sox4 is wound-induced and required for differentiation of multiple neural cell types NFY regulates wound-induced expression of sox4 during regeneration.

Acoel single-cell atlas reveals expression dynamics and heterogeneity of adult pluripotent stem cells

Adult pluripotent stem cell (aPSC) populations underlie whole-body regeneration in many distantly-related animal lineages, but how the underlying cellular and molecular mechanisms compare across species is unknown. Here, we apply single-cell RNA sequencing to profile transcriptional cell states of the acoel worm Hofstenia miamia during postembryonic development and regeneration. We identify cell types shared across stages and their associated gene expression dynamics during regeneration. Functional studies confirm that the aPSCs, also known as neoblasts, are the source of differentiated cells and reveal transcription factors needed for differentiation. Subclustering of neoblasts recovers transcriptionally distinct subpopulations, the majority of which are likely specialized to differentiated lineages. One neoblast subset, showing enriched expression of the histone variant H3.3, appears to lack specialization. Altogether, the cell states identified in this study facilitate comparisons to other species and enable future studies of stem cell fate potentials.

Cambrian stem-group ambulacrarians and the nature of the ancestral deuterostome

Deuterostomes are characterized by some of the most widely divergent body plans in the animal kingdom. These striking morphological differences have hindered efforts to predict ancestral characters, with the origin and earliest evolution of the group remaining ambiguous. Several iconic Cambrian fossils have been suggested to be early deuterostomes and hence could help elucidate ancestral character states. However, their phylogenetic relationships are controversial. Here, we describe new, exceptionally preserved specimens of the discoidal metazoan Rotadiscus grandis from the early Cambrian Chengjiang biota of China. These reveal a previously unknown double spiral structure, which we interpret as a chordate-like covering to a coelomopore, located adjacent to a horseshoe-shaped tentacle complex. The tentacles differ in key aspects from those seen in lophophorates and are instead more similar to the tentacular systems of extant pterobranchs and echinoderms. Thus, Rotadiscus exhibits a chimeric combination of ambulacrarian and chordate characters. Phylogenetic analyses recover Rotadiscus and closely related fossil taxa as stem ambulacrarians, filling a significant morphological gap in the deuterostome tree of life. These results allow us to reconstruct the ancestral body plans of major clades of deuterostomes, revealing that key traits of extant forms, such as a post-anal region, gill bars, and a U-shaped gut, evolved through convergence.

Peripheral and central employment of acid-sensing ion channels during early bilaterian evolution

Nervous systems are endowed with rapid chemosensation and intercellular signaling by ligand-gated ion channels (LGICs). While a complex, bilaterally symmetrical nervous system is a major innovation of bilaterian animals, the employment of specific LGICs during early bilaterian evolution is poorly understood. We therefore questioned bilaterian animals' employment of acid-sensing ion channels (ASICs), LGICs that mediate fast excitatory responses to decreases in extracellular pH in vertebrate neurons. Our phylogenetic analysis identified an earlier emergence of ASICs from the overarching DEG/ENaC (degenerin/epithelial sodium channel) superfamily than previously thought and suggests that ASICs were a bilaterian innovation. Our broad examination of ASIC gene expression and biophysical function in each major bilaterian lineage of Xenacoelomorpha, Protostomia, and Deuterostomia suggests that the earliest bilaterian ASICs were probably expressed in the periphery, before being incorporated into the brain as it emerged independently in certain deuterostomes and xenacoelomorphs. The loss of certain peripheral cells from Ecdysozoa after they separated from other protostomes likely explains their loss of ASICs, and thus the absence of ASICs from model organisms Drosophila and Caenorhabditis elegans. Thus, our use of diverse bilaterians in the investigation of LGIC expression and function offers a unique hypothesis on the employment of LGICs in early bilaterian evolution.

Induced spawning with gamete release from body ruptures during reproduction of Xenoturbella bocki

Xenoturbella is a marine invertebrate with a simple body plan, with recent phylogenomic studies suggesting that it forms the phylum Xenacoelomorpha together with the acoelomorphs. The phylogenetic position of the phylum is still under debate, whether it is an early branching bilaterian or a sister group to the Ambulacraria. Phylogenetic traits often appear during development, and larva resembling the cnidarian planula has been reported for Xenoturbella. However, subsequent developmental studies on Xenoturbella have been scarce. This is mainly due to the difficulties in collecting and keeping adult animals, resulting in the lack of data on the reproduction of the animal, such as the breeding season and the spawning pattern. Here we report on the reproduction of X. bocki and confirm that its breeding season is winter. Spawning induction resulted in gametes being released from body ruptures and not the mouth. No evidence supported the animal as a simultaneous hermaphrodite.

Genome assembly of the acoel flatworm Symsagittifera roscoffensis, a model for research on body plan evolution and photosymbiosis

Symsagittifera roscoffensis is a well-known member of the order Acoela that lives in symbiosis with the algae Tetraselmis convolutae during its adult stage. Its natural habitat is the eastern coast of the Atlantic, where at specific locations thousands of individuals can be found, mostly, lying in large pools on the surface of sand at low tide. As a member of the Acoela it has been thought as a proxy for ancestral bilaterian animals; however, its phylogenetic position remains still debated. In order to understand the basic structural characteristics of the acoel genome, we sequenced and assembled the genome of aposymbiotic species S. roscoffensis. The size of this genome was measured to be in the range of 910-940 Mb. Sequencing of the genome was performed using PacBio Hi-Fi technology. Hi-C and RNA-seq data were also generated to scaffold and annotate it. The resulting assembly is 1.1 Gb large (covering 118% of the estimated genome size) and highly continuous, with N50 scaffold size of 1.04 Mb. The repetitive fraction of the genome is 61%, of which 85% (half of the genome) are LTR retrotransposons. Genome-guided transcriptome assembly identified 34,493 genes, of which 29,351 are protein coding (BUSCO score 97.6%), and 30.2% of genes are spliced leader trans-spliced. The completeness of this genome suggests that it can be used extensively to characterize gene families and conduct accurate phylogenomic reconstructions.

Phylogenomics and the flowering plant tree of life

The advances accelerated by next-generation sequencing and long-read sequencing technologies continue to provide an impetus for plant phylogenetic study. In the past decade, a large number of phylogenetic studies adopting hundreds to thousands of genes across a wealth of clades have emerged and ushered plant phylogenetics and evolution into a new era. In the meantime, a roadmap for researchers when making decisions across different approaches for their phylogenomic research design is imminent. This review focuses on the utility of genomic data (from organelle genomes, to both reduced representation sequencing and whole-genome sequencing) in phylogenetic and evolutionary investigations, describes the baseline methodology of experimental and analytical procedures, and summarizes recent progress in flowering plant phylogenomics at the ordinal, familial, tribal, and lower levels. We also discuss the challenges, such as the adverse impact on orthology inference and phylogenetic reconstruction raised from systematic errors, and underlying biological factors, such as whole-genome duplication, hybridization/introgression, and incomplete lineage sorting, together suggesting that a bifurcating tree may not be the best model for the tree of life. Finally, we discuss promising avenues for future plant phylogenomic studies.

Worms and gills, plates and spines: the evolutionary origins and incredible disparity of deuterostomes revealed by fossils, genes, and development

Deuterostomes are the major division of animal life which includes sea stars, acorn worms, and humans, among a wide variety of ecologically and morphologically disparate taxa. However, their early evolution is poorly understood, due in part to their disparity, which makes identifying commonalities difficult, as well as their relatively poor early fossil record. Here, we review the available morphological, palaeontological, developmental, and molecular data to establish a framework for exploring the origins of this important and enigmatic group. Recent fossil discoveries strongly support a vermiform ancestor to the group Hemichordata, and a fusiform active swimmer as ancestor to Chordata. The diverse and anatomically bewildering variety of forms among the early echinoderms show evidence of both bilateral and radial symmetry. We consider four characteristics most critical for understanding the form and function of the last common ancestor to Deuterostomia: Hox gene expression patterns, larval morphology, the capacity for biomineralization, and the morphology of the pharyngeal region. We posit a deuterostome last common ancestor with a similar antero-posterior gene regulatory system to that found in modern acorn worms and cephalochordates, a simple planktonic larval form, which was later elaborated in the ambulacrarian lineage, the ability to secrete calcium minerals in a limited fashion, and a pharyngeal respiratory region composed of simple pores. This animal was likely to be motile in adult form, as opposed to the sessile origins that have been historically suggested. Recent debates regarding deuterostome monophyly as well as the wide array of deuterostome-affiliated problematica further suggest the possibility that those features were not only present in the last common ancestor of Deuterostomia, but potentially in the ur-bilaterian. The morphology and development of the early deuterostomes, therefore, underpin some of the most significant questions in the study of metazoan evolution.

Exploring genome gene content and morphological analysis to test recalcitrant nodes in the animal phylogeny

An accurate phylogeny of animals is needed to clarify their evolution, ecology, and impact on shaping the biosphere. Although datasets of several hundred thousand amino acids are nowadays routinely used to test phylogenetic hypotheses, key deep nodes in the metazoan tree remain unresolved: the root of animals, the root of Bilateria, and the monophyly of Deuterostomia. Instead of using the standard approach of amino acid datasets, we performed analyses of newly assembled genome gene content and morphological datasets to investigate these recalcitrant nodes in the phylogeny of animals. We explored extensively the choices for assembling the genome gene content dataset and model choices of morphological analyses. Our results are robust to these choices and provide additional insights into the early evolution of animals, they are consistent with sponges as the sister group of all the other animals, the worm-like bilaterian lineage Xenacoelomorpha as the sister group of the other Bilateria, and tentatively support monophyletic Deuterostomia.

Resolving tricky nodes in the tree of life through amino acid recoding

Genomic data allowed a detailed resolution of the Tree of Life, but "tricky nodes" such as the root of the animals remain unresolved. Genome-scale datasets are heterogeneous as genes and species are exposed to different pressures, and this can negatively impacts phylogenetic accuracy. We use simulated genomic-scale datasets and show that recoding amino acid data improves accuracy when the model does not account for the compositional heterogeneity of the amino acid alignment. We apply our findings to three datasets addressing the root of the animal tree, where the debate centers on whether sponges (Porifera) or comb jellies (Ctenophora) represent the sister of all other animals. We show that results from empirical data follow predictions from simulations and suggest that, at the least in phylogenies inferred from amino acid sequences, a placement of the ctenophores as sister to all the other animals is best explained as a tree reconstruction artifact.

Phylogenomics: Is less more when using large-scale datasets?

Phylogenetic studies have traditionally placed the simple Xenoacoelomorph worms as the sister group of all other animals with bilateral body symmetry. A new study shows that misidentification of orthologous genes might have been the source of at least some support for this placement.

Embryonic origins of adult pluripotent stem cells

Although adult pluripotent stem cells (aPSCs) are found in many animal lineages, mechanisms for their formation during embryogenesis are unknown. Here, we leveraged Hofstenia miamia, a regenerative worm that possesses collectively pluripotent aPSCs called neoblasts and produces manipulable embryos. Lineage tracing and functional experiments revealed that one pair of blastomeres gives rise to cells that resemble neoblasts in distribution, behavior, and gene expression. In Hofstenia, aPSCs include transcriptionally distinct subpopulations that express markers associated with differentiated tissues; our data suggest that despite their heterogeneity, aPSCs are derived from one lineage, not from multiple tissue-specific lineages during development. Next, we combined single-cell transcriptome profiling across development with neoblast cell-lineage tracing and identified a molecular trajectory for neoblast formation that includes transcription factors Hes, FoxO, and Tbx. This identification of a cellular mechanism and molecular trajectory for aPSC formation opens the door for in vivo studies of aPSC regulation and evolution.

Filtering artifactual signal increases support for Xenacoelomorpha and Ambulacraria sister relationship in the animal tree of life

Conflicting studies place a group of bilaterian invertebrates containing xenoturbellids and acoelomorphs, the Xenacoelomorpha, as either the primary emerging bilaterian phylum^1,2,3,4,5,6 or within Deuterostomia, sister to Ambulacraria.^7,8,9,10,11 Although their placement as sister to the rest of Bilateria supports relatively simple morphology in the ancestral bilaterian, their alternative placement within Deuterostomia suggests a morphologically complex ancestral bilaterian along with extensive loss of major phenotypic traits in the Xenacoelomorpha. Recent studies have questioned whether Deuterostomia should be considered monophyletic at all.^10,12,13 Hidden paralogy and poor phylogenetic signal present a major challenge for reconstructing species phylogenies.^{14,15,16,17,18} Here, we assess whether these issues have contributed to the conflict over the placement of Xenacoelomorpha. We reanalyzed published datasets, enriching for orthogroups whose gene trees support well-resolved clans elsewhere in the animal tree.¹⁶ We find that most genes in previously published datasets violate incontestable clans, suggesting that hidden paralogy and low phylogenetic signal affect the ability to reconstruct branching patterns at deep nodes in the animal tree. We demonstrate that removing orthogroups that cannot recapitulate incontestable relationships alters the final topology that is inferred, while simultaneously improving the fit of the model to the data. We discover increased, but ultimately not conclusive, support for the existence of Xenambulacraria in our set of filtered orthogroups. At a time when we are progressing toward sequencing all life on the planet, we argue that long-standing contentious issues in the tree of life will be resolved using smaller amounts of better quality data that can be modeled adequately.¹⁹.

Seeking Sense in the Hox Gene Cluster

The Hox gene cluster, responsible for patterning of the head-tail axis, is an ancestral feature of all bilaterally symmetrical animals (the Bilateria) that remains intact in a wide range of species. We can say that the Hox cluster evolved successfully only once since it is commonly the same in all groups, with labial-like genes at one end of the cluster expressed in the anterior embryo, and Abd-B-like genes at the other end of the cluster expressed posteriorly. This review attempts to make sense of the Hox gene cluster and to address the following questions. How did the Hox cluster form in the protostome-deuterostome last common ancestor, and why was this with a particular head-tail polarity? Why is gene clustering usually maintained? Why is there collinearity between the order of genes along the cluster and the positions of their expressions along the embryo? Why do the Hox gene expression domains overlap along the embryo? Why have vertebrates duplicated the Hox cluster? Why do Hox gene knockouts typically result in anterior homeotic transformations? How do animals adapt their Hox clusters to evolve new structural patterns along the head-tail axis?

A Practical Guide to Design and Assess a Phylogenomic Study

Over the last decade, molecular systematics has undergone a change of paradigm as high-throughput sequencing now makes it possible to reconstruct evolutionary relationships using genome-scale datasets. The advent of "big data" molecular phylogenetics provided a battery of new tools for biologists but simultaneously brought new methodological challenges. The increase in analytical complexity comes at the price of highly specific training in computational biology and molecular phylogenetics, resulting very often in a polarized accumulation of knowledge (technical on one side and biological on the other). Interpreting the robustness of genome-scale phylogenetic studies is not straightforward, particularly as new methodological developments have consistently shown that the general belief of "more genes, more robustness" often does not apply, and because there is a range of systematic errors that plague phylogenomic investigations. This is particularly problematic because phylogenomic studies are highly heterogeneous in their methodology, and best practices are often not clearly defined. The main aim of this article is to present what I consider as the ten most important points to take into consideration when planning a well-thought-out phylogenomic study and while evaluating the quality of published papers. The goal is to provide a practical step-by-step guide that can be easily followed by nonexperts and phylogenomic novices in order to assess the technical robustness of phylogenomic studies or improve the experimental design of a project.

Renewed perspectives on the sedentary-pelagic last common bilaterian ancestor

Various evaluations of the last common bilaterian ancestor (lcba) currently suggest that it resembled either a microscopic, non-segmented motile adult; or, on the contrary, a complex segmented adult motile urbilaterian. These fundamental inconsistencies remain largely unexplained. A majority of multidisciplinary data regarding sedentary adult ancestral bilaterian organization is overlooked. The sedentary-pelagic model is supported now by a number of novel developmental, paleontological and molecular phylogenetic data: (1) data in support of sedentary sponges, in the adult stage, as sister to all other Metazoa; (2) a similarity of molecular developmental pathways in both adults and larvae across sedentary sponges, cnidarians, and bilaterians; (3) a cnidarian-bilaterian relationship, including a unique sharing of a bona fide Hox-gene cluster, of which the evolutionary appearance does not connect directly to a bilaterian motile organization; (4) the presence of sedentary and tube-dwelling representatives of the main bilaterian clades in the early Cambrian; (5) an absence of definite taxonomic attribution of Ediacaran taxa reconstructed as motile to any true bilaterian phyla; (6) a similarity of tube morphology (and the clear presence of a protoconch-like apical structure of the Ediacaran sedentary Cloudinidae) among shells of the early Cambrian, and later true bilaterians, such as semi-sedentary hyoliths and motile molluscs; (7) recent data that provide growing evidence for a complex urbilaterian, despite a continuous molecular phylogenetic controversy. The present review compares the main existing models and reconciles the sedentary model of an urbilaterian and the model of a larva-like lcba with a unified sedentary(adult)-pelagic(larva) model of the lcba.

Polyzoa is back: The effect of complete gene sets on the placement of Ectoprocta and Entoprocta

The phylogenomic approach has largely resolved metazoan phylogeny and improved our knowledge of animal evolution based on morphology, paleontology, and embryology. Nevertheless, the placement of two major lophotrochozoan phyla, Entoprocta (Kamptozoa) and Ectoprocta (Bryozoa), remains highly controversial: Originally considered as a single group named Polyzoa (Bryozoa), they were separated on the basis of morphology. So far, each new study of lophotrochozoan evolution has still consistently proposed different phylogenetic positions for these groups. Here, we reinvestigated the placement of Entoprocta and Ectoprocta using highly complete datasets with rigorous contamination removal. Our results from maximum likelihood, Bayesian, and coalescent analyses strongly support the topology in which Entoprocta and Bryozoa form a distinct clade, placed as a sister group to all other lophotrochozoan clades: Annelida, Mollusca, Brachiopoda, Phoronida, and Nemertea. Our study favors the evolutionary scenario where Entoprocta, Cycliophora, and Bryozoa constitute one of the earliest branches among Lophotrochozoa and thus supports the Polyzoa hypothesis.

A phylogenomic backbone for gastropod molluscs

Gastropods have survived several mass extinctions during their evolutionary history resulting in extraordinary diversity in morphology, ecology, and developmental modes, which complicate the reconstruction of a robust phylogeny. Currently, gastropods are divided into six subclasses: Caenogastropoda, Heterobranchia, Neomphaliones, Neritimorpha, Patellogastropoda, and Vetigastropoda. Phylogenetic relationships among these taxa historically lack consensus, despite numerous efforts using morphological and molecular information. We generated sequence data for transcriptomes derived from twelve taxa belonging to clades with little or no prior representation in previous studies in order to infer the deeper cladogenetic events within Gastropoda and, for the first time, infer the position of the deep-sea Neomphaliones using a phylogenomic approach. We explored the impact of missing data, homoplasy, and compositional heterogeneity on the inferred phylogenetic hypotheses. We recovered a highly supported backbone for gastropod relationships that is congruent with morphological and mitogenomic evidence, in which Patellogastropoda, true limpets, are the sister lineage to all other gastropods (Orthogastropoda) which are divided into two main clades (i) Vetigastropoda s.l. (including Pleurotomariida + Neomphaliones) and (ii) Neritimorpha + (Caenogastropoda + Heterobranchia). As such, our results support the recognition of five subclasses (or infraclasses) in Gastropoda: Patellogastropoda, Vetigastropoda, Neritimorpha, Caenogastropoda and Heterobranchia.

Conflict over the eukaryote root resides in strong outliers, mosaics and missing data sensitivity of site-specific (CAT) mixture models

Phylogenetic reconstruction using concatenated loci ("phylogenomics" or "supermatrix phylogeny") is a powerful tool for solving evolutionary splits that are poorly resolved in single gene/protein trees (SGTs). However, recent phylogenomic attempts to resolve the eukaryote root have yielded conflicting results, along with claims of various artefacts hidden in the data. We have investigated these conflicts using two new methods for assessing phylogenetic conflict. ConJak uses whole marker (gene or protein) jackknifing to assess deviation from a central mean for each individual sequence, while ConWin uses a sliding window to screen for incongruent protein fragments (mosaics). Both methods allow selective masking of individual sequences or sequence fragments in order to minimize missing data, an important consideration for resolving deep splits with limited data. Analyses focused on a set of 76 eukaryotic proteins of bacterial-ancestry previously used in various combinations to assess the branching order among the three major divisions of eukaryotes: Amorphea (mainly animals, fungi and Amoebozoa), Diaphoretickes (most other well-known eukaryotes and nearly all algae) and Excavata, represented here by Discoba (Jakobida, Heterolobosea, and Euglenozoa). ConJak analyses found strong outliers to be concentrated in under-sampled lineages, while ConWin analyses of Discoba, the most under-sampled of the major lineages, detected potentially incongruent fragments scattered throughout. Phylogenetic analyses of the full data using an LG-gamma model support a Discoba sister scenario (neozoan-excavate root), which rises to 99-100% bootstrap support with data masked according to either protocol. However, analyses with two site-specific (CAT) mixture models yielded widely inconsistent results and a striking sensitivity to missing data. The neozoan-excavate root places Amorphea and Diaphoretickes as more closely related to each other than either is to Discoba, a fundamental relationship that should remain unaffected by additional taxa.

Evolution and loss of ß-catenin and TCF-dependent axis specification in insects

Mechanisms and evolution of primary axis specification in insects are discussed in the context of the roles of ß-catenin and TCF in polarizing metazoan embryos. Three hypotheses are presented. First, insects with sequential segmentation and posterior growth use cell-autonomous mechanisms for establishing embryo polarity via the nuclear ratio of ß-catenin and TCF. Second, TCF homologs establish competence for anterior specification. Third, the evolution of simultaneous segmentation mechanisms, also known as long-germ development, resulted in primary axis specification mechanisms that are independent of ß-catenin but reliant on TCF, a condition that preceded the frequent replacement of anterior determinants in long germ insects.

Studying development, regeneration, stem cells, and more in the acoel Hofstenia miamia

Acoel worms represent an enigmatic lineage of animals (Acoelomorpha) that has danced around the tree of animal life. Morphology-based classification placed them as flatworms (Phylum Platyhelminthes), with much of their biology being interpreted as a variation on what is observed in better-studied members of that phylum. However, molecular phylogenies suggest that acoels belong to a clade (Xenacoelomorpha) that could be a sister group to other animals with bilateral symmetry (Bilateria) or could belong within deuterostomes, closely related to a group that includes sea stars (Ambulacraria). This change in phylogenetic position has led to renewed interest in the biology of acoels, which can now offer insights into the evolution of many bilaterian traits. The acoel Hofstenia miamia has emerged as a powerful new research organism that enables mechanistic studies of xenacoelomorph biology, especially of developmental and regenerative processes. This article explains the motivation for developing Hofstenia as a new model system, describes Hofstenia biology, highlights the tools and resources that make Hofstenia a good research organism, and considers the questions that Hofstenia is well-positioned to answer. Finally, looking to the future, this article serves as an invitation to new and established scientists to join the growing community of researchers studying this exciting model system.

How to build phylogenetic species trees with OMA

Knowledge of species phylogeny is critical to many fields of biology. In an era of genome data availability, the most common way to make a phylogenetic species tree is by using multiple protein-coding genes, conserved in multiple species. This methodology is composed of several steps: orthology inference, multiple sequence alignment and inference of the phylogeny with dedicated tools. This can be a difficult task, and orthology inference, in particular, is usually computationally intensive and error prone if done ad hoc. This tutorial provides protocols to make use of OMA Orthologous Groups, a set of genes all orthologous to each other, to infer a phylogenetic species tree. It is designed to be user-friendly and computationally inexpensive, by providing two options: (1) Using only precomputed groups with species available on the OMA Browser, or (2) Computing orthologs using OMA Standalone for additional species, with the option of using precomputed orthology relations for those present in OMA. A protocol for downstream analyses is provided as well, including creating a supermatrix, tree inference, and visualization. All protocols use publicly available software, and we provide scripts and code snippets to facilitate data handling. The protocols are accompanied with practical examples.

Emergence of distinct syntenic density regimes is associated with early metazoan genomic transitions

Background

Animal genomes are strikingly conserved in terms of local gene order (microsynteny). While some of these microsyntenies have been shown to be coregulated or to form gene regulatory blocks, the diversity of their genomic and regulatory properties across the metazoan tree of life remains largely unknown.

Results

Our comparative analyses of 49 animal genomes reveal that the largest gains of synteny occurred in the last common ancestor of bilaterians and cnidarians and in that of bilaterians. Depending on their node of emergence, we further show that novel syntenic blocks are characterized by distinct functional compositions (Gene Ontology terms enrichment) and gene density properties, such as high, average and low gene density regimes. This is particularly pronounced among bilaterian novel microsyntenies, most of which fall into high gene density regime associated with higher gene coexpression levels. Conversely, a majority of vertebrate novel microsyntenies display a low gene density regime associated with lower gene coexpression levels.

Conclusions

Our study provides first evidence for evolutionary transitions between different modes of microsyntenic block regulation that coincide with key events of metazoan evolution. Moreover, the microsyntenic profiling strategy and interactive online application (Syntenic Density Browser, available at: http://synteny.csb.univie.ac.at/) we present here can be used to explore regulatory properties of microsyntenic blocks and predict their coexpression in a wide-range of animal genomes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08304-2.

Somatostatin-type and allatostatin-C-type neuropeptides are paralogous and have opposing myoregulatory roles in an echinoderm

Somatostatin (SS) and allatostatin-C (ASTC) are related neuropeptide hormones that act as inhibitory regulators of physiological processes in chordates (e.g., humans) and protostome invertebrates (e.g., insects), respectively. We have discovered that echinoderms (e.g., starfish) uniquely have both SS-type and ASTC-type neuropeptides, which act as inhibitory and excitatory regulators of muscle activity, respectively. Our findings suggest that SS-type and ASTC-type neuropeptides evolved by duplication of a common ancestral encoding gene. Then, one of the neuropeptides was lost in protostomes and chordates, probably because of their functional redundancy as inhibitory regulators. Conversely, the unique retention of both neuropeptide types in echinoderms may be explained by evolution of an excitatory role for ASTC-type neuropeptides mediated by yet-to-be-determined signaling mechanisms.

Somatostatin (SS) and allatostatin-C (ASTC) are inhibitory neuropeptides in chordates and protostomes, respectively, which hitherto were identified as orthologs. However, echinoderms have two SS/ASTC-type neuropeptides (SS1 and SS2), and here, our analysis of sequence data indicates that SS1 is an ortholog of ASTC and SS2 is an ortholog of SS. The occurrence of both SS-type and ASTC-type neuropeptides in echinoderms provides a unique context to compare their physiological roles. Investigation of the expression and actions of the ASTC-type neuropeptide ArSS1 in the starfish Asterias rubens revealed that it causes muscle contraction (myoexcitation), contrasting with myoinhibitory effects of the SS-type neuropeptide ArSS2. Our findings suggest that SS-type and ASTC-type neuropeptides are paralogous and originated by gene duplication in a common ancestor of the Bilateria, with only one type being retained in chordates (SS) and protostomes (ASTC) but with both types being retained in echinoderms. Loss of ASTC-type and SS-type neuropeptides in chordates and protostomes, respectively, may have been due to their functional redundancy as inhibitory regulators of physiological processes. Conversely, the retention of both neuropeptide types in echinoderms may be a consequence of the evolution of a myoexcitatory role for ASTC-type neuropeptides mediated by as yet unknown signaling mechanisms.

Studying Xenacoelomorpha WBR Using Isodiametra pulchra

Xenacoelomorpha are a phylogenetically and biologically interesting, but severely understudied group of worm-like animals. Among them, the acoel Isodiametra pulchra has been shown to be amenable to experimental work, including the study of stem cells and regeneration. The animal is capable of regenerating the posterior part of the body, but not its head. Here, methods such as nucleic acid extractions, in situ hybridisation, RNA interference, antibody and cytochemical stainings, and the general handling of the animals are presented.

Functional Conservation and Genetic Divergence of Chordate Glycinergic Neurotransmission: Insights from Amphioxus Glycine Transporters

Glycine is an important neurotransmitter in vertebrates, performing both excitatory and inhibitory actions. Synaptic levels of glycine are tightly controlled by the action of two glycine transporters, GlyT1 and GlyT2, located on the surface of glial cells and neurons, respectively. Only limited information is available on glycinergic neurotransmission in invertebrates, and the evolution of glycinergic neurotransmission is poorly understood. Here, by combining phylogenetic and gene expression analyses, we characterized the glycine transporter complement of amphioxus, an important invertebrate model for studying the evolution of chordates. We show that amphioxus possess three glycine transporter genes. Two of these (GlyT2.1 and GlyT2.2) are closely related to GlyT2 of vertebrates, whereas the third (GlyT) is a member of an ancestral clade of deuterostome glycine transporters. GlyT2.2 expression is predominantly non-neural, whereas GlyT and GlyT2.1 are widely expressed in the amphioxus nervous system and are differentially expressed, respectively, in neurons and glia. Vertebrate glycinergic neurons express GlyT2 and glia GlyT1, suggesting that the evolution of the chordate glycinergic system was accompanied by a paralog-specific inversion of gene expression. Despite this genetic divergence between amphioxus and vertebrates, we found strong evidence for conservation in the role glycinergic neurotransmission plays during larval swimming, the implication being that the neural networks controlling the rhythmic movement of chordate bodies may be homologous.

Transgenesis in the acoel worm Hofstenia miamia

The acoel worm Hofstenia miamia, which can replace tissue lost to injury via differentiation of a population of stem cells, has emerged as a new research organism for studying regeneration. To enhance the depth of mechanistic studies in this system, we devised a protocol for microinjection into embryonic cells that resulted in stable transgene integration into the genome and generated animals with tissue-specific fluorescent transgene expression in epidermis, gut, and muscle. We demonstrate that transgenic Hofstenia are amenable to the isolation of specific cell types, investigations of regeneration, tracking of photoconverted molecules, and live imaging. Further, our stable transgenic lines revealed insights into the biology of Hofstenia, including a high-resolution three-dimensional view of cell morphology and the organization of muscle as a cellular scaffold for other tissues. Our work positions Hofstenia as a powerful system with multiple toolkits for mechanistic investigations of development, whole-body regeneration, and stem cell biology.

Citrullination Was Introduced into Animals by Horizontal Gene Transfer from Cyanobacteria

Protein posttranslational modifications add great sophistication to biological systems. Citrullination, a key regulatory mechanism in human physiology and pathophysiology, is enigmatic from an evolutionary perspective. Although the citrullinating enzymes peptidylarginine deiminases (PADIs) are ubiquitous across vertebrates, they are absent from yeast, worms, and flies. Based on this distribution PADIs were proposed to have been horizontally transferred, but this has been contested. Here, we map the evolutionary trajectory of PADIs into the animal lineage. We present strong phylogenetic support for a clade encompassing animal and cyanobacterial PADIs that excludes fungal and other bacterial homologs. The animal and cyanobacterial PADI proteins share functionally relevant primary and tertiary synapomorphic sequences that are distinct from a second PADI type present in fungi and actinobacteria. Molecular clock calculations and sequence divergence analyses using the fossil record estimate the last common ancestor of the cyanobacterial and animal PADIs to be less than 1 billion years old. Additionally, under an assumption of vertical descent, PADI sequence change during this evolutionary time frame is anachronistically low, even when compared with products of likely endosymbiont gene transfer, mitochondrial proteins, and some of the most highly conserved sequences in life. The consilience of evidence indicates that PADIs were introduced from cyanobacteria into animals by horizontal gene transfer (HGT). The ancestral cyanobacterial PADI is enzymatically active and can citrullinate eukaryotic proteins, suggesting that the PADI HGT event introduced a new catalytic capability into the regulatory repertoire of animals. This study reveals the unusual evolution of a pleiotropic protein modification.

Lower statistical support with larger datasets: insights from the Ochrophyta radiation

It is commonly assumed that increasing the number of characters has the potential to resolve evolutionary radiations. Here, we studied photosynthetic stramenopiles (Ochrophyta) using alignments of heterogeneous origin mitochondrion, plastid, and nucleus. Surprisingly while statistical support for the relationships between the six major Ochrophyta lineages increases when comparing the mitochondrion (6,762 sites) and plastid (21,692 sites) trees, it decreases in the nuclear (209,105 sites) tree. Statistical support is not simply related to the data set size but also to the quantity of phylogenetic signal available at each position and our ability to extract it. Here, we show that this ability for current phylogenetic methods is limited, because conflicting results were obtained when varying taxon sampling. Even though the use of a better fitting model improved signal extraction and reduced the observed conflicts, the plastid data set provided higher statistical support for the ochrophyte radiation than the larger nucleus data set. We propose that the higher support observed in the plastid tree is due to an acceleration of the evolutionary rate in one short deep internal branch, implying that more phylogenetic signal per position is available to resolve the Ochrophyta radiation in the plastid than in the nuclear data set. Our work therefore suggests that, in order to resolve radiations, beyond the obvious use of data sets with more positions, we need to continue developing models of sequence evolution that better extract the phylogenetic signal and design methods to search for genes/characters that contain more signal specifically for short internal branches.

Phylogenomic Subsampling and the Search for Phylogenetically Reliable Loci

Phylogenomic subsampling is a procedure by which small sets of loci are selected from large genome-scale data sets and used for phylogenetic inference. This step is often motivated by either computational limitations associated with the use of complex inference methods or as a means of testing the robustness of phylogenetic results by discarding loci that are deemed potentially misleading. Although many alternative methods of phylogenomic subsampling have been proposed, little effort has gone into comparing their behavior across different data sets. Here, I calculate multiple gene properties for a range of phylogenomic data sets spanning animal, fungal, and plant clades, uncovering a remarkable predictability in their patterns of covariance. I also show how these patterns provide a means for ordering loci by both their rate of evolution and their relative phylogenetic usefulness. This method of retrieving phylogenetically useful loci is found to be among the top performing when compared with alternative subsampling protocols. Relatively common approaches such as minimizing potential sources of systematic bias or increasing the clock-likeness of the data are found to fare worse than selecting loci at random. Likewise, the general utility of rate-based subsampling is found to be limited: loci evolving at both low and high rates are among the least effective, and even those evolving at optimal rates can still widely differ in usefulness. This study shows that many common subsampling approaches introduce unintended effects in off-target gene properties and proposes an alternative multivariate method that simultaneously optimizes phylogenetic signal while controlling for known sources of bias.

Phylogenomic Subsampling and the Search for Phylogenetically Reliable Loci

Phylogenomic subsampling is a procedure by which small sets of loci are selected from large genome-scale data sets and used for phylogenetic inference. This step is often motivated by either computational limitations associated with the use of complex inference methods or as a means of testing the robustness of phylogenetic results by discarding loci that are deemed potentially misleading. Although many alternative methods of phylogenomic subsampling have been proposed, little effort has gone into comparing their behavior across different data sets. Here, I calculate multiple gene properties for a range of phylogenomic data sets spanning animal, fungal, and plant clades, uncovering a remarkable predictability in their patterns of covariance. I also show how these patterns provide a means for ordering loci by both their rate of evolution and their relative phylogenetic usefulness. This method of retrieving phylogenetically useful loci is found to be among the top performing when compared with alternative subsampling protocols. Relatively common approaches such as minimizing potential sources of systematic bias or increasing the clock-likeness of the data are found to fare worse than selecting loci at random. Likewise, the general utility of rate-based subsampling is found to be limited: loci evolving at both low and high rates are among the least effective, and even those evolving at optimal rates can still widely differ in usefulness. This study shows that many common subsampling approaches introduce unintended effects in off-target gene properties and proposes an alternative multivariate method that simultaneously optimizes phylogenetic signal while controlling for known sources of bias.

Genomics and transcriptomics of epizoic Seisonidea (Rotifera, syn. Syndermata) reveal strain formation and gradual gene loss with growing ties to the host

Background

Seisonidea (also Seisonacea or Seisonidae) is a group of small animals living on marine crustaceans (Nebalia spec.) with only four species described so far. Its monophyletic origin with mostly free-living wheel animals (Monogononta, Bdelloidea) and endoparasitic thorny-headed worms (Acanthocephala) is widely accepted. However, the phylogenetic relationships inside the Rotifera-Acanthocephala clade (Rotifera sensulato or Syndermata) are subject to ongoing debate, with consequences for our understanding of how genomes and lifestyles might have evolved. To gain new insights, we analyzed first drafts of the genome and transcriptome of the key taxon Seisonidea.

Results

Analyses of gDNA-Seq and mRNA-Seq data uncovered two genetically distinct lineages in Seison nebaliae Grube, 1861 off the French Channel coast. Their mitochondrial haplotypes shared only 82% sequence identity despite identical gene order. In the nuclear genome, distinct linages were reflected in different gene compactness, GC content and codon usage. The haploid nuclear genome spans ca. 46 Mb, of which 96% were reconstructed. According to ~ 23,000 SuperTranscripts, gene number in S. nebaliae should be within the range published for other members of Rotifera-Acanthocephala. Consistent with this, numbers of metazoan core orthologues and ANTP-type transcriptional regulatory genes in the S. nebaliae genome assembly were between the corresponding numbers in the other assemblies analyzed. We additionally provide evidence that a basal branching of Seisonidea within Rotifera-Acanthocephala could reflect attraction to the outgroup. Accordingly, rooting via a reconstructed ancestral sequence led to monophyletic Pararotatoria (Seisonidea+Acanthocephala) within Hemirotifera (Bdelloidea+Pararotatoria).

Conclusion

Matching genome/transcriptome metrics with the above phylogenetic hypothesis suggests that a haploid nuclear genome of about 50 Mb represents the plesiomorphic state for Rotifera-Acanthocephala. Smaller genome size in S. nebaliae probably results from subsequent reduction. In contrast, genome size should have increased independently in monogononts as well as bdelloid and acanthocephalan stem lines. The present data additionally indicate a decrease in gene repertoire from free-living to epizoic and endoparasitic lifestyles. Potentially, this reflects corresponding steps from the root of Rotifera-Acanthocephala via the last common ancestors of Hemirotifera and Pararotatoria to the one of Acanthocephala. Lastly, rooting via a reconstructed ancestral sequence may prove useful in phylogenetic analyses of other deep splits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07857-y.