A transcriptome approach to ecdysozoan phylogeny

Developmental and genomic insight into the origin of the tardigrade body plan

Tardigrada is an ancient lineage of miniaturized animals. As an outgroup of the well-studied Arthropoda and Onychophora, studies of tardigrades hold the potential to reveal important insights into body plan evolution in Panarthropoda. Previous studies have revealed interesting facets of tardigrade development and genomics that suggest that a highly compact body plan is a derived condition of this lineage, rather than it representing an ancestral state of Panarthropoda. This conclusion was based on studies of several species from Eutardigrada. We review these studies and expand on them by analyzing the publicly available genome and transcriptome assemblies of Echiniscus testudo, a representative of Heterotardigrada. These new analyses allow us to phylogenetically reconstruct important features of genome evolution in Tardigrada. We use available data from tardigrades to interrogate several recent models of body plan evolution in Panarthropoda. Although anterior segments of panarthropods are highly diverse in terms of anatomy and development, both within individuals and between species, we conclude that a simple one-to-one alignment of anterior segments across Panarthropoda is the best available model of segmental homology. In addition to providing important insight into body plan diversification within Panarthropoda, we speculate that studies of tardigrades may reveal generalizable pathways to miniaturization.

De novo assembly of transcriptomes and differential gene expression analysis using short-read data from emerging model organisms - a brief guide

Many questions in biology benefit greatly from the use of a variety of model systems. High-throughput sequencing methods have been a triumph in the democratization of diverse model systems. They allow for the economical sequencing of an entire genome or transcriptome of interest, and with technical variations can even provide insight into genome organization and the expression and regulation of genes. The analysis and biological interpretation of such large datasets can present significant challenges that depend on the 'scientific status' of the model system. While high-quality genome and transcriptome references are readily available for well-established model systems, the establishment of such references for an emerging model system often requires extensive resources such as finances, expertise and computation capabilities. The de novo assembly of a transcriptome represents an excellent entry point for genetic and molecular studies in emerging model systems as it can efficiently assess gene content while also serving as a reference for differential gene expression studies. However, the process of de novo transcriptome assembly is non-trivial, and as a rule must be empirically optimized for every dataset. For the researcher working with an emerging model system, and with little to no experience with assembling and quantifying short-read data from the Illumina platform, these processes can be daunting. In this guide we outline the major challenges faced when establishing a reference transcriptome de novo and we provide advice on how to approach such an endeavor. We describe the major experimental and bioinformatic steps, provide some broad recommendations and cautions for the newcomer to de novo transcriptome assembly and differential gene expression analyses. Moreover, we provide an initial selection of tools that can assist in the journey from raw short-read data to assembled transcriptome and lists of differentially expressed genes.

New insights into mesoderm and endoderm development, and the nature of the onychophoran blastopore

BACKGROUND

Early during onychophoran development and prior to the formation of the germ band, a posterior tissue thickening forms the posterior pit. Anterior to this thickening forms a groove, the embryonic slit, that marks the anterior-posterior orientation of the developing embryo. This slit is by some authors considered the blastopore, and thus the origin of the endoderm, while others argue that the posterior pit represents the blastopore. This controversy is of evolutionary significance because if the slit represents the blastopore, then this would support the amphistomy hypothesis that suggests that a slit-like blastopore in the bilaterian ancestor evolved into protostomy and deuterostomy.

RESULTS

In this paper, we summarize our current knowledge about endoderm and mesoderm development in onychophorans and provide additional data on early endoderm- and mesoderm-determining marker genes such as Blimp, Mox, and the T-box genes.

CONCLUSION

We come to the conclusion that the endoderm of onychophorans forms prior to the development of the embryonic slit, and thus that the slit is not the primary origin of the endoderm. It is thus unlikely that the embryonic slit represents the blastopore. We suggest instead that the posterior pit indeed represents the lips of the blastopore, and that the embryonic slit (and surrounding tissue) represents a morphologically superficial archenteron-like structure. We conclude further that both endoderm and mesoderm development are under control of conserved gene regulatory networks, and that many of the features found in arthropods including the model Drosophila melanogaster are likely derived.

Musculature of an Early Cambrian cycloneuralian animal

Cycloneuralians are ecdysozoans with a fossil record extending to the Early Cambrian Fortunian Age and represented mostly by cuticular integuments. However, internal anatomies of Fortunian cycloneuralians are virtually unknown, hampering our understanding of their functional morphology and phylogenetic relationships. Here we report the exceptional preservation of cycloneuralian introvert musculature in Fortunian rocks of South China. The musculature consists of an introvert body-wall muscular grid of four circular and 36 radially arranged longitudinal muscle bundles, as well as an introvert circular muscle associated with 19 roughly radially arranged, short retractors. Collectively, these features support at least a scalidophoran affinity, and the absence of muscles associated with a mouth cone and scalids further indicates a priapulan affinity. As in modern scalidophorans, the fossil musculature, and particularly the introvert circular muscle retractors, may have controlled introvert inversion and facilitated locomotion and feeding. This work supports the evolution of scalidophoran-like or priapulan-like introvert musculature in cycloneuralians at the beginning of the Cambrian Period.

Neglected no longer: Phylogenomic resolution of higher-level relationships in Solifugae

Advanced sequencing technologies have expedited resolution of higher-level arthropod relationships. Yet, dark branches persist, principally among groups occurring in cryptic habitats. Among chelicerates, Solifugae ("camel spiders") is the last order lacking a higher-level phylogeny and have thus been historically characterized as "neglected [arachnid] cousins". Though renowned for aggression, remarkable running speed, and xeric adaptation, inferring solifuge relationships has been hindered by inaccessibility of diagnostic morphological characters, whereas molecular investigations have been limited to one of 12 recognized families. Our phylogenomic dataset via capture of ultraconserved elements sampling all extant families recovered a well-resolved phylogeny, with two distinct groups of New World taxa nested within a broader Paleotropical radiation. Divergence times using fossil calibrations inferred that Solifugae radiated by the Permian, and most families diverged prior to the Paleogene-Cretaceous extinction, likely driven by continental breakup. We establish Boreosolifugae new suborder uniting five Laurasian families, and Australosolifugae new suborder uniting seven Gondwanan families using morphological and biogeographic signal.

A novel, conserved and possibly functional motif "WHWGHTW" in mitochondrial transcription across Bilateria

The animal mitogenomes which undergone a reductive evolution has an obvious loss of coding capacity compared to their known closest relatives, but it has not yet been fully investigated why and how the intergenic regions do not encode protein and have no known functions, are stably maintained, replicated, and transmitted by the genome. These relatively small intergenic regions may not be under neutral evolution and they may have functional and/or regulatory roles that have yet to be identified. Here, the distribution pattern, sequence content and location of a novel sequence motif of 'WWWGHTW' were bioinformatically investigated and characterised by constructing a sampling mitogenome dataset of 1889 species from 14 phyla representing the clade of Bilateria. This motif is reverse complementary of the previously described DmTTF binding sequence and found in the nd4L- (X) -trnT gene cluster. This cluster commonly exhibits a strand displacement region and an intergenic region among the bilaterian superphylums, particularly in Ecdysozoa. This motif may be accepted as a substrate providing binding sites for the specific interaction with transcription factors because of (i) its reverse complementarity of previously described DmTTF binding sequence, and (ii) the possession of G and T nucleotides in the fourth and sixth positions, (iii) the bias on T and G nucleotides instead of C and A in the degenerated positions. This suggestion is also supported by the presence of a strand displacement region in the nd4L- (X) -trnT gene cluster, particularly in Ecdysozoa consisting of the most rearranged mitogenomes among the bilaterian superphylums.

The unbearable uncertainty of panarthropod relationships

Panarthropoda, the clade comprising the phyla Onychophora, Tardigrada and Euarthropoda, encompasses the largest majority of animal biodiversity. The relationships among the phyla are contested and resolution is key to understanding the evolutionary assembly of panarthropod bodyplans. Molecular phylogenetic analyses generally support monophyly of Onychophora and Euarthropoda to the exclusion of Tardigrada (Lobopodia hypothesis), which is also supported by some analyses of morphological data. However, analyses of morphological data have also been interpreted to support monophyly of Tardigrada and Euarthropoda to the exclusion of Onychophora (Tactopoda hypothesis). Support has also been found for a clade of Onychophora and Tardigrada that excludes Euarthropoda (Protarthropoda hypothesis). Here we show, using a diversity of phylogenetic inference methods, that morphological datasets cannot discriminate statistically between the Lobopodia, Tactopoda and Protarthropoda hypotheses. Since the relationships among the living clades of panarthropod phyla cannot be discriminated based on morphological data, we call into question the accuracy of morphology-based phylogenies of Panarthropoda that include fossil species and the evolutionary hypotheses based upon them.

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.

How many long branch orders occur in Chelicerata? Opposing effects of Palpigradi and Opilioacariformes on phylogenetic stability

Excepting a handful of nodes, phylogenetic relationships between chelicerate orders remains poorly resolved, due to both the incidence of long branch attraction artifacts and the limited sampling of key lineages. It has recently been shown that increasing representation of basal nodes plays an outsized role in resolving the higher-level placement of long branch chelicerate orders. Two lineages have been consistently undersampled in chelicerate phylogeny. First, sampling of the miniaturized order Palpigradi has been restricted to a fragmentary transcriptome of a single species. Second, sampling of Opilioacariformes, a rarely encountered and key group of Parasitiformes, has been restricted to a single exemplar. These two lineages exhibit dissimilar properties with respect to branch length; Opilioacariformes shows relatively low evolutionary rate compared to other Parasitiformes, whereas Palpigradi possibly acts as another long branch order (an effect that may be conflated with the degree of missing data). To assess these properties and their effects on tree stability, we constructed a phylogenomic dataset of Chelicerata wherein both lineages were sampled with three terminals, increasing the representation of these lineages per locus. We examined the effect of subsampling phylogenomic matrices using (1) taxon occupancy, (2) evolutionary rate, and (3) a principal components-based approach. We further explored the impact of taxon deletion experiments that mitigate the effect of long branches. Here, we show that Palpigradi constitutes a fourth long branch chelicerate order (together with Acariformes, Parasitiformes, and Pseudoscorpiones), which further destabilizes the chelicerate backbone topology. By contrast, the slow-evolving Opilioacariformes were consistently recovered within Parasitiformes, with certain subsampling practices recovering their placement as the sister group to the remaining Parasitiformes. Whereas the inclusion of Opilioacariformes always resulted in the non-monophyly of Acari with support, deletion of Opilioacariformes from datasets consistently incurred the monophyly of Acari except in matrices constructed on the basis of evolutionary rate. Our results strongly suggest that Acari is an artifact of long-branch attraction.

Phylogenomic Analysis of Velvet Worms (Onychophora) Uncovers an Evolutionary Radiation in the Neotropics

Onychophora ("velvet worms") are charismatic soil invertebrates known for their status as a "living fossil," their phylogenetic affiliation to arthropods, and their distinctive biogeographic patterns. However, several aspects of their internal phylogenetic relationships remain unresolved, limiting our understanding of the group's evolutionary history, particularly with regard to changes in reproductive mode and dispersal ability. To address these gaps, we used RNA sequencing and phylogenomic analysis of transcriptomes to reconstruct the evolutionary relationships and infer divergence times within the phylum. We recovered a fully resolved and well-supported phylogeny for the circum-Antarctic family Peripatopsidae, which retains signals of Gondwanan vicariance and showcases the evolutionary lability of reproductive mode in the family. Within the Neotropical clade of Peripatidae, though, we found that amino acid-translated sequence data masked nearly all phylogenetic signal, resulting in highly unstable and poorly supported relationships. Analyses using nucleotide sequence data were able to resolve many more relationships, though we still saw discordant phylogenetic signal between genes, probably indicative of a rapid, mid-Cretaceous radiation in the group. Finally, we hypothesize that the unique reproductive mode of placentotrophic viviparity found in all Neotropical peripatids may have facilitated the multiple inferred instances of over-water dispersal and establishment on oceanic islands.

Taxonomic Sampling and Rare Genomic Changes Overcome Long-Branch Attraction in the Phylogenetic Placement of Pseudoscorpions

Long-branch attraction is a systematic artifact that results in erroneous groupings of fast-evolving taxa. The combination of short, deep internodes in tandem with long-branch attraction artifacts has produced empirically intractable parts of the Tree of Life. One such group is the arthropod subphylum Chelicerata, whose backbone phylogeny has remained unstable despite improvements in phylogenetic methods and genome-scale data sets. Pseudoscorpion placement is particularly variable across data sets and analytical frameworks, with this group either clustering with other long-branch orders or with Arachnopulmonata (scorpions and tetrapulmonates). To surmount long-branch attraction, we investigated the effect of taxonomic sampling via sequential deletion of basally branching pseudoscorpion superfamilies, as well as varying gene occupancy thresholds in supermatrices. We show that concatenated supermatrices and coalescent-based summary species tree approaches support a sister group relationship of pseudoscorpions and scorpions, when more of the basally branching taxa are sampled. Matrix completeness had demonstrably less influence on tree topology. As an external arbiter of phylogenetic placement, we leveraged the recent discovery of an ancient genome duplication in the common ancestor of Arachnopulmonata as a litmus test for competing hypotheses of pseudoscorpion relationships. We generated a high-quality developmental transcriptome and the first genome for pseudoscorpions to assess the incidence of arachnopulmonate-specific duplications (e.g., homeobox genes and miRNAs). Our results support the inclusion of pseudoscorpions in Arachnopulmonata (new definition), as the sister group of scorpions. Panscorpiones (new name) is proposed for the clade uniting Scorpiones and Pseudoscorpiones.

Ancestral morphology of Ecdysozoa constrained by an early Cambrian stem group ecdysozoan

BACKGROUND

Ecdysozoa are the moulting protostomes, including arthropods, tardigrades, and nematodes. Both the molecular and fossil records indicate that Ecdysozoa is an ancient group originating in the terminal Proterozoic, and exceptional fossil biotas show their dominance and diversity at the beginning of the Phanerozoic. However, the nature of the ecdysozoan common ancestor has been difficult to ascertain due to the extreme morphological diversity of extant Ecdysozoa, and the lack of early diverging taxa in ancient fossil biotas.

RESULTS

Here we re-describe Acosmia maotiania from the early Cambrian Chengjiang Biota of Yunnan Province, China and assign it to stem group Ecdysozoa. Acosmia features a two-part body, with an anterior proboscis bearing a terminal mouth and muscular pharynx, and a posterior annulated trunk with a through gut. Morphological phylogenetic analyses of the protostomes using parsimony, maximum likelihood and Bayesian inference, with coding informed by published experimental decay studies, each placed Acosmia as sister taxon to Cycloneuralia + Panarthropoda-i.e. stem group Ecdysozoa. Ancestral state probabilities were calculated for key ecdysozoan nodes, in order to test characters inferred from fossils to be ancestral for Ecdysozoa. Results support an ancestor of crown group ecdysozoans sharing an annulated vermiform body with a terminal mouth like Acosmia, but also possessing the pharyngeal armature and circumoral structures characteristic of Cambrian cycloneuralians and lobopodians.

CONCLUSIONS

Acosmia is the first taxon placed in the ecdysozoan stem group and provides a constraint to test hypotheses on the early evolution of Ecdysozoa. Our study suggests acquisition of pharyngeal armature, and therefore a change in feeding strategy (e.g. predation), may have characterised the origin and radiation of crown group ecdysozoans from Acosmia-like ancestors.

Arthropod Origins: Integrating Paleontological and Molecular Evidence

Phylogenomics underpins a stable and mostly well-resolved hypothesis for the interrelationships of extant arthropods. Exceptionally preserved fossils are integrated into this framework by coding their morphological characters, as exemplified by total-evidence dating approaches that treat fossils as dated tips in analyses numerically dominated by molecular data. Cambrian fossils inform on the sequence of character acquisition in the arthropod stem group and in the stems of its main extant clades. The arthropod head problem incorporates unique appendage combinations and remains of the nervous system in fossils into a scheme mostly based on neuroanatomy and Hox expression domains for extant forms. Molecular estimates of arthropod origins in the Cryogenian or Ediacaran predate a coherent picture from the arthropod fossil record, which commences as trace fossils in the earliest Cambrian. Probabilistic morphological clock analysis of trilobites, which exemplify the earliest arthropod body fossils, supports a Cambrian origin, without the need to posit an unfossilized Ediacaran history.

Naturally occurring fluorescence protects the eutardigrade Paramacrobiotus sp. from ultraviolet radiation

Naturally occurring fluorescence has been observed in multiple species ranging from bacteria to birds. In macroscopic animals such as birds, fluorescence provides a visual communication signal. However, the functional significance of this phenomenon is unknown in most cases. Though photoprotection is attributed to fluorescence under ultraviolet (UV) light in some organisms, it lacks direct experimental evidence. Here, we demonstrate naturally occurring fluorescence under UV light in a eutardigrade belonging to the genus Paramacrobiotus. Using a natural variant that lacks fluorescence, we show that the fluorescence confers tolerance to lethal UV radiation. Remarkably, the fluorescent extract from Paramacrobiotus sp. could protect the UV-sensitive tardigrade Hypsibius exemplaris and nematode Caenorhabditis elegans from germicidal UV radiation. We propose that Paramacrobiotus sp. possess a protective fluorescent shield that absorbs harmful UV radiation and emits harmless blue light.

Function and Distribution of the Wamide Neuropeptide Superfamily in Metazoans

The Wamide neuropeptide superfamily is of interest due to its distinctive functions in regulating life cycle transitions, metamorphic hormone signaling, and several aspects of digestive system function, from gut muscle contraction to satiety and fat storage. Due to variation among researchers in naming conventions, a global view of Wamide signaling in animals in terms of conservation or diversification of function is currently lacking. Here, I summarize the phylogenetic distribution of Wamide neuropeptides based on current data and describe recent findings in the areas of Wamide receptors and biological functions. Common trends that emerge across Cnidarians and protostomes are the presence of multiple Wamide receptors within a single organism, and the fact that Wamide signaling likely functions across an extensive variety of biological systems, including visual, circadian, and reproductive systems. Important areas of focus for future research are the further identification of Wamide-receptor pairs, confirmation of the phylogenetic distribution of Wamides through largescale sequencing and mass spectrometry, and assignment of different functions to specific subsets of Wamide-expressing neurons. More extensive study of Wamide signaling throughout larval development in a greater number of phyla is also important in order to understand the role of Wamides in hormonal regulation. Defining the evolution and function of neuropeptide signaling in animal nervous systems will benefit from an increased understanding of Wamide function and signaling mechanisms in a wider variety of organisms, beyond the traditional model systems.

Cuticular reticulation replicates the pattern of epidermal cells in lowermost Cambrian scalidophoran worms

The cuticle of ecdysozoans (Panarthropoda, Scalidophora, Nematoida) is secreted by underlying epidermal cells and renewed via ecdysis. We explore here the relationship between epidermis and external cuticular ornament in stem-group scalidophorans from the early Cambrian of China (Kuanchuanpu Formation; ca 535 Ma) that had two types of microscopic polygonal cuticular networks with either straight or microfolded boundaries. Detailed comparisons with modern scalidophorans (priapulids) indicate that these networks faithfully replicate the cell boundaries of the epidermis. This suggests that the cuticle of early scalidophorans formed through the fusion between patches of extracellular material secreted by epidermal cells, as observed in various groups of present-day ecdysozoans, including arthropods. Key genetic, biochemical and mechanical processes associated with ecdysis and cuticle formation seem to have appeared very early (at least not later than 535 Ma) in the evolution of ecdysozoans. Microfolded reticulation is likely to be a mechanical response to absorbing contraction exerted by underlying muscles. The polygonal reticulation in early and extant ecdysozoans is clearly a by-product of the epidermal cell pavement and interacted with the sedimentary environment.

The Unique Antimicrobial Recognition and Signaling Pathways in Tardigrades with a Comparison Across Ecdysozoa

Tardigrades are microscopic animals known to withstand unfavorable abiotic conditions. These animals are also constantly exposed to biotic stresses, including parasites and internal microbiomes. However, the tardigrade immune mechanisms against these biotic stresses are largely uncharacterized. Due to the contentious phylogenetic position of tardigrades, it is not intuitive whether they possess an immune system more similar to that of arthropods (e.g., Toll, Imd, and JNK pathways of the Drosophila melanogaster antimicrobial response) or to that of nematodes (e.g., the Tir-1/Nsy-1/Sek-1/Pmk-1/Atf-7 signaling cassette [called Tir-1 pathway here]) in Caenorhabditis elegans). In this study, comparative genomic analyses were conducted to mine homologs of canonical D. melanogaster and C. elegans immune pathway genes from eight tardigrades (Echiniscoides cf. sigismundi, Echiniscus testudo, Hypsibius exemplaris, Mesobiotus philippinicus, Milnesium tardigradum, Paramacrobiotus richtersi, Richtersius cf. coronifer, and Ramazzottius varieornatus) and four non-arthropod ecdysozoans (two onychophorans: Epiperipatus sp. and Opisthopatus kwazululandi; one nematomorph: Paragordius varius; and one priapulan: Priapulus caudatus) in order to provide insights into the tardigrade antimicrobial system. No homologs of the intracellular components of the Toll pathway were detected in any of the tardigrades examined. Likewise, no homologs of most of the Imd pathway genes were detected in any of the tardigrades or any of the other non-arthropod ecdysozoans. Both the JNK and Tir-1 pathways, on the other hand, were found to be conserved across ecdysozoans. Interestingly, tardigrades had no detectable homologs of NF-κB, the major activator of antimicrobial response gene expression. Instead, tardigrades appear to possess NF-κB distantly related NFAT homologs. Overall, our results show that tardigrades have a unique gene pathway repertoire that differs from that of other ecdysozoans. Our study also provides a framework for future studies on tardigrade immune responses.

Major Revisions in Arthropod Phylogeny Through Improved Supermatrix, With Support for Two Possible Waves of Land Invasion by Chelicerates

Deep phylogeny involving arthropod lineages is difficult to recover because the erosion of phylogenetic signals over time leads to unreliable multiple sequence alignment (MSA) and subsequent phylogenetic reconstruction. One way to alleviate the problem is to assemble a large number of gene sequences to compensate for the weakness in each individual gene. Such an approach has led to many robustly supported but contradictory phylogenies. A close examination shows that the supermatrix approach often suffers from two shortcomings. The first is that MSA is rarely checked for reliability and, as will be illustrated, can be poor. The second is that, to alleviate the problem of homoplasy at the third codon position of protein-coding genes due to convergent evolution of nucleotide frequencies, phylogeneticists may remove or degenerate the third codon position but may do it improperly and introduce new biases. We performed extensive reanalysis of one of such "big data" sets to highlight these two problems, and demonstrated the power and benefits of correcting or alleviating these problems. Our results support a new group with Xiphosura and Arachnopulmonata (Tetrapulmonata + Scorpiones) as sister taxa. This favors a new hypothesis in which the ancestor of Xiphosura and the extinct Eurypterida (sea scorpions, of which many later forms lived in brackish or freshwater) returned to the sea after the initial chelicerate invasion of land. Our phylogeny is supported even with the original data but processed with a new "principled" codon degeneration. We also show that removing the 1673 codon sites with both AGN and UCN codons (encoding serine) in our alignment can partially reconcile discrepancies between nucleotide-based and AA-based tree, partly because two sequences, one with AGN and the other with UCN, would be identical at the amino acid level but quite different at the nucleotide level.

Ordered phylogenomic subsampling enables diagnosis of systematic errors in the placement of the enigmatic arachnid order Palpigradi

The miniaturized arachnid order Palpigradi has ambiguous phylogenetic affinities owing to its odd combination of plesiomorphic and derived morphological traits. This lineage has never been sampled in phylogenomic datasets because of the small body size and fragility of most species, a sampling gap of immediate concern to recent disputes over arachnid monophyly. To redress this gap, we sampled a population of the cave-inhabiting species Eukoenenia spelaea from Slovakia and inferred its placement in the phylogeny of Chelicerata using dense phylogenomic matrices of up to 1450 loci, drawn from high-quality transcriptomic libraries and complete genomes. The complete matrix included exemplars of all extant orders of Chelicerata. Analyses of the complete matrix recovered palpigrades as the sister group of the long-branch order Parasitiformes (ticks) with high support. However, sequential deletion of long-branch taxa revealed that the position of palpigrades is prone to topological instability. Phylogenomic subsampling approaches that maximized taxon or dataset completeness recovered palpigrades as the sister group of camel spiders (Solifugae), with modest support. While this relationship is congruent with the location and architecture of the coxal glands, a long-forgotten character system that opens in the pedipalpal segments only in palpigrades and solifuges, we show that nodal support values in concatenated supermatrices can mask high levels of underlying topological conflict in the placement of the enigmatic Palpigradi.

Comparative myoanatomy of Tardigrada: new insights from the heterotardigrades Actinarctus doryphorus (Tanarctidae) and Echiniscoides sigismundi (Echiniscoididae)

BACKGROUND

Tardigrada is a group of microscopic invertebrates distributed worldwide in permanent and temporal aquatic habitats. Famous for their extreme stress tolerance, tardigrades are also of interest due to their close relationship with Arthropoda and Cycloneuralia. Despite recent efforts in analyzing the musculature of a number of tardigrade species, data on the class Heterotardigrada remain scarce. Aiming to expand the current morphological framework, and to promote the use of muscular body plans in elucidating tardigrade phylogeny, the myoanatomy of two heterotardigrades, Actinarctus doryphorus and Echiniscoides sigismundi, was analyzed by cytochemistry, scanning electron and confocal laser scanning microscopy and 3D imaging. We discuss our findings with reference to other tardigrades and internal phylogenetic relationships of the phylum.

RESULTS

We focus our analyses on the somatic musculature, which in tardigrades includes muscle groups spanning dorsal, ventral, and lateral body regions, with the legs being musculated by fibers belonging to all three groups. A pronounced reduction of the trunk musculature is seen in the dorsoventrally compressed A. doryphorus, a species that generally has fewer cuticle attachment sites as compared to E. sigismundi and members of the class Eutardigrada. Interestingly, F-actin positive signals were found in the head appendages of A. doryphorus. Our analyses further indicate that cross-striation is a feature common to the somatic muscles of heterotardigrades and that E. sigismundi-as previously proposed for other echiniscoidean heterotardigrades-has relatively thick somatic muscle fibers.

CONCLUSIONS

We provide new insights into the myoanatomical differences that characterize distinct evolutionary lineages within Tardigrada, highlighting characters that potentially can be informative in future phylogenetic analyses. We focus our current analyses on the ventral trunk musculature. Our observations suggest that seven paired ventromedian attachment sites anchoring a large number of muscles can be regarded as part of the ground pattern of Tardigrada and that fusion and reduction of cuticular attachment sites is a derived condition. Specifically, the pattern of these sites differs in particular details between tardigrade taxa. In the future, a deeper understanding of the tardigrade myoanatomical ground pattern will require more investigations in order to include all major tardigrade lineages.

Origin of ecdysis: fossil evidence from 535-million-year-old scalidophoran worms

With millions of extant species, ecdysozoans (Scalidophora, Nematoida and Panarthropoda) constitute a major portion of present-day biodiversity. All ecdysozoans secrete an exoskeletal cuticle which must be moulted periodically and replaced by a larger one. Although moulting (ecdysis) has been recognized in early Palaeozoic panarthropods such as trilobites and basal groups such as anomalocaridids and lobopodians, the fossil record lacks clear evidence of ecdysis in early scalidophorans, largely because of difficulties in recognizing true exuviae. Here, we describe two types of exuviae in microscopic scalidophoran worms from the lowermost Cambrian Kuanchuanpu Formation ( ca 535 Ma) of China and reconstruct their moulting process. These basal scalidophorans moulted in a manner similar to that of extant priapulid worms, extricating themselves smoothly from their old tubular cuticle or turning their exuviae inside out like the finger of a glove. This is the oldest record of moulting in ecdysozoans. We also discuss the origin of ecdysis in the light of recent molecular analyses and the significance of moulting in the early evolution of animals.

Revisiting metazoan phylogeny with genomic sampling of all phyla

Proper biological interpretation of a phylogeny can sometimes hinge on the placement of key taxa-or fail when such key taxa are not sampled. In this light, we here present the first attempt to investigate (though not conclusively resolve) animal relationships using genome-scale data from all phyla. Results from the site-heterogeneous CAT + GTR model recapitulate many established major clades, and strongly confirm some recent discoveries, such as a monophyletic Lophophorata, and a sister group relationship between Gnathifera and Chaetognatha, raising continued questions on the nature of the spiralian ancestor. We also explore matrix construction with an eye towards testing specific relationships; this approach uniquely recovers support for Panarthropoda, and shows that Lophotrochozoa (a subclade of Spiralia) can be constructed in strongly conflicting ways using different taxon- and/or orthologue sets. Dayhoff-6 recoding sacrifices information, but can also reveal surprising outcomes, e.g. full support for a clade of Lophophorata and Entoprocta + Cycliophora, a clade of Placozoa + Cnidaria, and raising support for Ctenophora as sister group to the remaining Metazoa, in a manner dependent on the gene and/or taxon sampling of the matrix in question. Future work should test the hypothesis that the few remaining uncertainties in animal phylogeny might reflect violations of the various stationarity assumptions used in contemporary inference methods.

Embryonic expression of priapulid Wnt genes

Posterior elongation of the developing embryo is a common feature of animal development. One group of genes that is involved in posterior elongation is represented by the Wnt genes, secreted glycoprotein ligands that signal to specific receptors on neighbouring cells and thereby establish cell-to-cell communication. In segmented animals such as annelids and arthropods, Wnt signalling is also likely involved in segment border formation and regionalisation of the segments. Priapulids represent unsegmented worms that are distantly related to arthropods. Despite their interesting phylogenetic position and their importance for the understanding of ecdysozoan evolution, priapulids still represent a highly underinvestigated group of animals. Here, we study the embryonic expression patterns of the complete sets of Wnt genes in the priapulids Priapulus caudatus and Halicryptus spinulosus. We find that both priapulids possess a complete set of 12 Wnt genes. At least in Priapulus, most of these genes are expressed in and around the posterior-located blastopore and thus likely play a role in posterior elongation. Together with previous work on the expression of other genetic factors such as caudal and even-skipped, this suggests that posterior elongation in priapulids is under control of the same (or very similar) conserved gene regulatory network as in arthropods.

Improved phylogenomic sampling of free-living nematodes enhances resolution of higher-level nematode phylogeny

Background

Nematodes are among the most diverse and abundant metazoans on Earth, but research on them has been biased toward parasitic taxa and model organisms. Free-living nematodes, particularly from the clades Enoplia and Dorylaimia, have been underrepresented in genome-scale phylogenetic analyses to date, leading to poor resolution of deep relationships within the phylum.

Results

We supplemented publicly available data by sequencing transcriptomes of nine free-living nematodes and two important outgroups and conducted a phylum-wide phylogenomic analysis including a total of 108 nematodes. Analysis of a dataset generated using a conservative orthology inference strategy resulted in a matrix with a high proportion of missing data and moderate to weak support for branching within and placement of Enoplia. A less conservative orthology inference approach recovered more genes and resulted in higher support for the deepest splits within Nematoda, recovering Enoplia as the sister taxon to the rest of Nematoda. Relationships within major clades were similar to those found in previously published studies based on 18S rDNA.

Conclusions

Expanded transcriptome sequencing of free-living nematodes has contributed to better resolution among deep nematode lineages, though the dataset is still strongly biased toward parasites. Inclusion of more free-living nematodes in future phylogenomic analyses will allow a clearer understanding of many interesting aspects of nematode evolution, such as morphological and molecular adaptations to parasitism and whether nematodes originated in a marine or terrestrial environment.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1444-x) contains supplementary material, which is available to authorized users.

Phylogenomics of the longitarsal Colossendeidae: The evolutionary history of an Antarctic sea spider radiation

Sea spiders (Pycnogonida) constitute a group of marine benthic arthropods that has a particularly high species diversity in the Southern Ocean. The "longitarsal" group of the sea spider family Colossendeidae is especially abundant in this region. However, this group also includes some representatives from other oceans, which raises the question where the group originates from. Therefore, we here investigated the phylogeny of the group with a hybrid enrichment approach that yielded a dataset of 1607 genes and over one million base pairs. We obtained a well-resolved phylogeny of the group, which is mostly consistent with morphological data. The data support an Antarctic origin of the longitarsal Colossendeidae and multiple dispersal events to other regions, which occurred at different timescales. This scenario is consistent with evidence found in other groups of marine invertebrates and highlights the role of the Southern Ocean as a source for non-Antarctic biota, especially of the deep sea. Our results suggest an initially slow rate of diversification followed by a more rapid radiation possibly correlated with the mid-Miocene cooling of Antarctica, similar to what is found in other taxa.

Neuroanatomy of mud dragons: a comprehensive view of the nervous system in Echinoderes (Kinorhyncha) by confocal laser scanning microscopy

Background

The Scalidophora (Kinorhyncha, Loricifera and Priapulida) have an important phylogenetic position as early branching ecdysozoans, yet the architecture of their nervous organ systems is notably underinvestigated. Without such information, and in the absence of a stable phylogenetic context, we are inhibited from producing adequate hypotheses about the evolution and diversification of ecdysozoan nervous systems. Here, we utilize confocal laser scanning microscopy to characterize serotonergic, tubulinergic and FMRFamidergic immunoreactivity patterns in a comparative neuroanatomical study with three species of Echinoderes, the most speciose, abundant and diverse genus within Kinorhyncha.

Results

Neuroanatomy in Echinoderes as revealed by acetylated α-tubulin immunoreactivity includes a circumpharyngeal brain and ten neurite bundles in the head region that converge into five longitudinal nerves within the trunk. The ventral nerve cord is ganglionated, emerging from the brain with two connectives that converge in trunk segments 2–3, and diverge again within segment 8. The longitudinal nerves and ventral nerve cord are connected by two transverse neurites in segments 2–9. Differences among species correlate with the number, position and innervation of cuticular structures along the body. Patterns of serotoninergic and FMRFamidergic immunoreactivity correlate with the position of the brain neuropil and the ventral nerve cord. Distinct serotonergic and FMRFamidergic somata are associated with the brain neuropil and specific trunk segments along the ventral nerve cord.

Conclusions

Neural architecture is highly conserved across all three species, suggesting that our results reveal a pattern that is common to more than 40% of the species within Kinorhyncha. The nervous system of Echinoderes is segmented along most of the trunk; however, posterior trunk segments exhibit modifications that are likely associated with sensorial, motor or reproductive functions. Although all kinorhynchs show some evidence of an externally segmented trunk, it is unclear whether external segmentation matches internal segmentation of nervous and muscular organ systems across Kinorhyncha, as we observed in Echinoderes. The neuroanatomical data provided in this study not only expand the limited knowledge on kinorhynch nervous systems but also establish a comparative morphological framework within Scalidophora that will support broader inferences about the evolution of neural architecture among the deepest branching lineages of the Ecdysozoa.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1405-4) contains supplementary material, which is available to authorized users.

A Critical Appraisal of the Placement of Xiphosura (Chelicerata) with Account of Known Sources of Phylogenetic Error

Horseshoe crabs (Xiphosura) are traditionally regarded as sister group to the clade of terrestrial chelicerates (Arachnida). This hypothesis has been challenged by recent phylogenomic analyses, but the non-monophyly of Arachnida has consistently been disregarded as artifactual. We re-evaluated the placement of Xiphosura among chelicerates using the most complete phylogenetic data set to date, expanding outgroup sampling, and including data from whole genome sequencing projects. In spite of uncertainty in the placement of some arachnid clades, all analyses show Xiphosura consistently nested within Arachnida as the sister group to Ricinulei (hooded tick spiders). It is apparent that the radiation of arachnids is an old one and occurred over a brief period of time, resulting in several consecutive short internodes, and thus is a potential case for the confounding effects of incomplete lineage sorting (ILS). We simulated coalescent gene trees to explore the effects of increasing levels of ILS on the placement of horseshoe crabs. In addition, common sources of systematic error were evaluated, as well as the effects of fast-evolving partitions and the dynamics of problematic long branch orders. Our results indicated that the placement of horseshoe crabs cannot be explained by missing data, compositional biases, saturation, or ILS. Interrogation of the phylogenetic signal showed that the majority of loci favor the derived placement of Xiphosura over a monophyletic Arachnida. Our analyses support the inference that horseshoe crabs represent a group of aquatic arachnids, comparable to aquatic mites, breaking a long-standing paradigm in chelicerate evolution and altering previous interpretations of the ancestral transition to the terrestrial habitat. Future studies testing chelicerate relationships should approach the task with a sampling strategy where the monophyly of Arachnida is not held as the premise.

The Natural History of Teneurins: A Billion Years of Evolution in Three Key Steps

The entire evolutionary history of the animal gene family, Teneurin, can be summed up in three key steps, plus three salient footnotes. In a shared ancestor of all bilaterians, the first step began with gene fusions that created a protein with an amino-terminal intracellular domain bridged via a single transmembrane helix to extracellular EGF-like domains. This first step was completed with a further gene fusion: an additional carboxy-terminal stretch of about 2000 amino acids (aa) was adopted, as-a-whole, from bacteria. The 2000 aa structure in Teneurin was recently solved in three dimensions. The 2000 aa region appears in a number of bacteria, yet was co-opted solely into Teneurin, and into no other eukaryotic proteins. Outside of bilaterian animals, no Teneurins exist, with a “Monosiga brevicollis caveat” brought below, as ‘the third footnote.” Subsequent to the “urTeneurin’s” genesis-by-fusions, all bilaterians bore a single Teneurin gene, always encoding an extraordinarily conserved Type II transmembrane protein with invariant domain content and order. The second key step was a duplication that led to an exception to singleton Teneurin genomes. A pair of Teneurin paralogs, Ten-a and Ten-m, are found in representatives of all four Arthropod sub-phyla, in: insects, crustaceans, myriapods, and chelicerates. In contrast, in every other protostome species’ genome, including those of all non-Arthropod ecdysozoan phyla, only a single Teneurin gene occurs. The closest, sister, phylum of arthropods, the Onychophorans (velvet worms), bear a singleton Teneurin. Ten-a and Ten-m therefore arose from a duplication in an urArthropod only after Arthropods split from Onychophorans, but before the splits that led to the four Arthropod sub-phyla. The third key step was a quadruplication of Teneurins at the root of vertebrate radiation. Four Teneurin paralogs (Teneurins 1 through 4) arose first by a duplication of a single chordate gene likely leading to one 1/4–type gene, and one 2/3-type gene: the two copies found in extant jawless vertebrates. Relatively soon thereafter, a second duplication round yielded the -1, -2, -3, and -4 paralog types now found in all jawed vertebrates, from sharks to humans. It is possible to assert that these duplication events correlate well to the Ohno hypothesized 2R (two round) vertebrate whole genome duplication (WGD), as refined in more recent treatments. The quadruplication can therefore be placed at approximately 400 Myr ago. Echinoderms, hemichordates, cephalochordates, and urochordates have only a single copy of Teneurin in their genomes. These deuterostomes and non-vertebrate chordates provide the anchor showing that the quadruplication happened at the root of vertebrates. A first footnote must be brought concerning some of the ‘invertebrate’ relatives of vertebrates, among Deuterostomes. A family of genes that encode 7000 aa proteins was derived from, but is distinct from, the Teneurin family. This distinct family arose early in deuterostomes, yet persists today only in hemichordate and cephalochordate genomes. They are named here TRIPs (Teneurin-related immense proteins). As a second of three ‘footnotes’: a limited number of species exist with additional Teneurin gene copies. However, these further duplications of Teneurins occur for paralog types (a, m, or 1–4) only in specific lineages within Arthropods or Vertebrates. All examples are paralog duplications that evidently arose in association with lineage specific WGDs. The increased Teneurin paralog numbers correlate with WGDs known and published in bony fish, Xenopus, plus select Chelicerates lineages and Crustaceans. The third footnote, alluded to above, is that a Teneurin occurs in one unicellular species: Monosiga brevicollis. Teneurins are solely a metazoan, bilaterian-specific family, to the exclusion of the Kingdoms of prokaryotes, plants, fungi, and protists. The single exception occurs among the unicellular, opisthokont, closest relatives of metazoans, the choanoflagellates. There is a Teneurin in Monosiga brevicollis, one species of the two fully sequenced choanoflagellate species. In contrast, outside of triploblast-bilaterians, there are no Teneurins in any diploblast genomes, including even sponges – those metazoans closest to choanoflagellates. Perhaps the ‘birth’ of the original Teneurin occurred in a shared ancestor of M. brevicollis and metazoans, then was lost in M. brevicollis’ sister species, and was serially and repeatedly lost in all diploblast metazoans. Alternatively, and as favored above, it first arose in the ‘urBilaterian,’ then was subsequently acquired from some bilaterian via horizontal transfer by a single choanoflagellate clade. The functional partnership of Teneurins and Latrophilins was discovered in rodents through the LPH1-TENM2 interaction. Recent work extends this to further members of each family. Surveying when the interacting domains of Teneurins and Latrophilins co-exist within different organisms can give an indication of how widespread their functional cooperation might be, across bilaterians. Paralog number for the two families is relatively correlated among bilaterians, and paralog numbers underwent co-increase in the WGDs mentioned above. With co-increasing paralog numbers, the possible combinatorial pairs grow factorially. This should have a significant impact for increasing nervous system complexity. The 3 key events in the ‘natural history’ of the Teneurins and their Latrophilin partners coincide with the ascendance of particularly successful metazoan clades: bilaterians; arthropods; and vertebrates. Perhaps we can attribute some of this success to the unique Teneurin family, and to its partnership with Latrophilins.

Was the ancestral panarthropod mouth ventral or terminal?

It has recently been suggested that the panarthropod mouth was ancestrally terminal, based on the assumption that the ancestral tardigrade had a terminal mouth and on the observations of a terminal mouth in adults of some stem-group fossils. This is shown to be unlikely, and it is concluded that the ancestral panarthropod had a ventral mouth.

Taxon-specific expansion and loss of tektins inform metazoan ciliary diversity

Background

Cilia and flagella are complex cellular structures thought to have first evolved in a last ciliated eukaryotic ancestor due to the conserved 9 + 2 microtubule doublet structure of the axoneme and associated proteins. The Tektin family of coiled-coil domain containing proteins was previously identified in cilia of organisms as diverse as green algae and sea urchin. While studies have shown that some Tektins are necessary for ciliary function, there has been no comprehensive phylogenetic survey of tektin genes. To fill this gap, we sampled tektin sequences broadly among metazoan and unicellular lineages in order to determine how the tektin gene complements evolved in over 100 different extant species.

Results

Using Bayesian and Maximum Likelihood analyses, we have ascertained with high confidence that all metazoan tektins arose from a single ancestral tektin gene in the last common ancestor of metazoans and choanoflagellates. Gene duplications gave rise to two tektin genes in the metazoan ancestor, and a subsequent expansion to three and four tektin genes in early bilaterian ancestors. While all four tektin genes remained highly conserved in most deuterostome and spiralian species surveyed, most tektin genes in ecdysozoans are highly derived with extensive gene loss in several lineages including nematodes and some crustaceans. In addition, while tektin-1, − 2, and − 4 have remained as single copy genes in most lineages, tektin-3/5 has been duplicated independently several times, notably at the base of the spiralian, vertebrate and hymenopteran (Ecdysozoa) clades.

Conclusions

We provide a solid description of tektin evolution supporting one, two, three, and four ancestral tektin genes in a holozoan, metazoan, bilaterian, and nephrozoan ancestor, respectively. The isolated presence of tektin in a cryptophyte and a chlorophyte branch invokes events of horizontal gene transfer, and that the last common ciliated eukaryotic ancestor lacked a tektin gene. Reconstructing the evolutionary history of the tektin complement in each extant metazoan species enabled us to pinpoint lineage specific expansions and losses. Our analysis will help to direct future studies on Tektin function, and how gain and loss of tektin genes might have contributed to the evolution of various types of cilia and flagella.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1360-0) contains supplementary material, which is available to authorized users.

Tetraconatan phylogeny with special focus on Malacostraca and Branchiopoda: highlighting the strength of taxon-specific matrices in phylogenomics

Understanding the evolution of Tetraconata or Pancrustacea-the clade that includes crustaceans and insects-requires a well-resolved hypothesis regarding the relationships within and among its constituent taxa. Here, we assembled a taxon-rich phylogenomic dataset focusing on crustacean lineages based solely on genomes and new-generation Illumina-generated transcriptomes, including 89 representatives of Tetraconata. This constitutes, to our knowledge, the first phylogenomic study specifically addressing internal relationships of Malacostraca (with 26 species included) and Branchiopoda (36 species). Seven matrices comprising 81-684 orthogroups and 17 690-242 530 amino acid positions were assembled and analysed under five different analytical approaches. To maximize gene occupancy and to improve resolution, taxon-specific matrices were designed for Malacostraca and Branchiopoda. Key tetraconatan taxa (i.e. Oligostraca, Multicrustacea, Branchiopoda, Malacostraca, Thecostraca, Copepoda and Hexapoda) were monophyletic and well supported. Within Branchiopoda, Phyllopoda, Diplostraca, Cladoceromorpha and Cladocera were monophyletic. Within Malacostraca, the clades Eumalacostraca, Decapoda and Reptantia were well supported. Recovery of Caridoida or Peracarida was highly dependent on the analysis for the complete matrix, but it was consistently monophyletic in the malacostracan-specific matrices. From such examples, we demonstrate that taxon-specific matrices and particular evolutionary models and analytical methods, namely CAT-GTR and Dayhoff recoding, outperform other approaches in resolving certain recalcitrant nodes in phylogenomic analyses.

Embryonic expression patterns and phylogenetic analysis of panarthropod sox genes: insight into nervous system development, segmentation and gonadogenesis

BACKGROUND

Sox (Sry-related high-mobility-group box) genes represent important factors in animal development. Relatively little, however, is known about the embryonic expression patterns and thus possible function(s) of Sox genes during ontogenesis in panarthropods (Arthropoda+Tardigrada+Onychophora). To date, studies have been restricted exclusively to higher insects, including the model system Drosophila melanogaster, with no comprehensive data available for any other arthropod group, or any tardigrade or onychophoran.

RESULTS

This study provides a phylogenetic analysis of panarthropod Sox genes and presents the first comprehensive analysis of embryonic expression patterns in the flour beetle Tribolium castaneum (Hexapoda), the pill millipede Glomeris marginata (Myriapoda), and the velvet worm, Euperipatoides kanangrensis (Onychophora). 24 Sox genes were identified and investigated: 7 in Euperipatoides, 8 in Glomeris, and 9 in Tribolium. Each species possesses at least one ortholog of each of the five expected Sox gene families, B, C, D, E, and F, many of which are differentially expressed during ontogenesis.

CONCLUSION

Sox gene expression (and potentially function) is highly conserved in arthropods and their closest relatives, the onychophorans. Sox B, C and D class genes appear to be crucial for nervous system development, while the Sox B genes Dichaete (D) and Sox21b likely play an additional conserved role in panarthropod segmentation. The Sox B gene Sox21a likely has a conserved function in foregut and Malpighian tubule development, at least in Hexapoda. The data further suggest that Sox D and E genes are involved in mesoderm differentiation, and that Sox E genes are involved in gonadal development. The new data expand our knowledge about the expression and implied function of Sox genes to Mandibulata (Myriapoda+Pancrustacea) and Panarthropoda (Arthropoda+Onychophora).

Transcriptomic Analysis of Pseudoscorpion Venom Reveals a Unique Cocktail Dominated by Enzymes and Protease Inhibitors

Transcriptomic and genomic analyses have illuminated the diversity of venoms in three of the four venomous arachnid orders (scorpions, spiders, and ticks). To date, no venom gland transcriptome analysis has been available for pseudoscorpions, the fourth venomous arachnid lineage. To redress this gap, we sequenced an mRNA library generated from the venom glands of the species Synsphyronus apimelus (Garypidae). High-throughput sequencing by the Illumina protocol, followed by de novo assembly, resulted in a total of 238,331 transcripts. From those, we annotated 131 transcripts, which code for putative peptides/proteins with similar sequences to previously reported venom components available from different arachnid species in protein databases. Transcripts putatively coding for enzymes showed the richest diversity, followed by other venom components such as peptidase inhibitors, cysteine-rich peptides, and thyroglobulin 1-like peptides. Only 11 transcripts were found that code for putatively low molecular mass spider toxins. This study constitutes the first report of the diversity of components within pseudoscorpion venom.

Resolving Relationships among the Megadiverse Butterflies and Moths with a Novel Pipeline for Anchored Phylogenomics

The advent of next-generation sequencing technology has allowed for thecollection of large portions of the genome for phylogenetic analysis. Hybrid enrichment and transcriptomics are two techniques that leverage next-generation sequencing and have shown much promise. However, methods for processing hybrid enrichment data are still limited. We developed a pipeline for anchored hybrid enrichment (AHE) read assembly, orthology determination, contamination screening, and data processing for sequences flanking the target "probe" region. We apply this approach to study the phylogeny of butterflies and moths (Lepidoptera), a megadiverse group of more than 157,000 described species with poorly understood deep-level phylogenetic relationships. We introduce a new, 855 locus AHE kit for Lepidoptera phylogenetics and compare resulting trees to those from transcriptomes. The enrichment kit was designed from existing genomes, transcriptomes, and expressed sequence tags and was used to capture sequence data from 54 species from 23 lepidopteran families. Phylogenies estimated from AHE data were largely congruent with trees generated from transcriptomes, with strong support for relationships at all but the deepest taxonomic levels. We combine AHE and transcriptomic data to generate a new Lepidoptera phylogeny, representing 76 exemplar species in 42 families. The tree provides robust support for many relationships, including those among the seven butterfly families. The addition of AHE data to an existing transcriptomic dataset lowers node support along the Lepidoptera backbone, but firmly places taxa with AHE data on the phylogeny. Combining taxa sequenced for AHE with existing transcriptomes and genomes resulted in a tree with strong support for (Calliduloidea $+$ Gelechioidea $+$ Thyridoidea) $+$ (Papilionoidea $+$ Pyraloidea $+$ Macroheterocera). To examine the efficacy of AHE at a shallow taxonomic level, phylogenetic analyses were also conducted on a sister group representing a more recent divergence, the Saturniidae and Sphingidae. These analyses utilized sequences from the probe region and data flanking it, nearly doubled the size of the dataset; resulting trees supported new phylogenetics relationships, especially within the Saturniidae and Sphingidae (e.g., Hemarina derived in the latter). We hope that our data processing pipeline, hybrid enrichment gene set, and approach of combining AHE data with transcriptomes will be useful for the broader systematics community.

A revised dated phylogeny of scorpions: Phylogenomic support for ancient divergence of the temperate Gondwanan family Bothriuridae

The scorpion family Bothriuridae occupies a subset of landmasses formerly constituting East and West temperate Gondwana, but its relationship to other scorpion families is in question. Whereas morphological data have strongly supported a sister group relationship of Bothriuridae and the superfamily Scorpionoidea, a recent phylogenomic analysis recovered a basal placement of bothriurids within Iurida, albeit sampling only a single exemplar. Here we reexamined the phylogenetic placement of the family Bothriuridae, sampling six bothriurid exemplars representing both East and West Gondwana, using transcriptomic data. Our results demonstrate that the sister group relationship of Bothriuridae to the clade ("Chactoidea" + Scorpionoidea) is supported by the inclusion of additional bothriurid taxa, and that this placement is insensitive to matrix completeness or partitioning by evolutionary rate. We also estimated divergence times within the order Scorpiones using multiple fossil calibrations, to infer whether the family Bothriuridae is sufficiently old to be characterized as a true Gondwanan lineage. We show that scorpions underwent ancient diversification between the Devonian and early Carboniferous. The age interval of the bothriurids sampled (a derived group that excludes exemplars from South Africa) spans the timing of breakup of temperate Gondwana.

Fossils and the Evolution of the Arthropod Brain

The discovery of fossilized brains and ventral nerve cords in lower and mid-Cambrian arthropods has led to crucial insights about the evolution of their central nervous system, the segmental identity of head appendages and the early evolution of eyes and their underlying visual systems. Fundamental ground patterns of lower Cambrian arthropod brains and nervous systems correspond to the ground patterns of brains and nervous systems belonging to three of four major extant panarthropod lineages. These findings demonstrate the evolutionary stability of early neural arrangements over an immense time span. Here, we put these fossil discoveries in the context of evidence from cladistics, as well as developmental and comparative neuroanatomy, which together suggest that despite many evolved modifications of neuropil centers within arthropod brains and ganglia, highly conserved arrangements have been retained. Recent phylogenies of the arthropods, based on fossil and molecular evidence, and estimates of divergence dates, suggest that neural ground patterns characterizing onychophorans, chelicerates and mandibulates are likely to have diverged between the terminal Ediacaran and earliest Cambrian, heralding the exuberant diversification of body forms that account for the Cambrian Explosion.

Mandibulate convergence in an armoured Cambrian stem chelicerate

BACKGROUND

Chelicerata represents a vast clade of mostly predatory arthropods united by a distinctive body plan throughout the Phanerozoic. Their origins, however, with respect to both their ancestral morphological features and their related ecologies, are still poorly understood. In particular, it remains unclear whether their major diagnostic characters were acquired early on, and their anatomical organization rapidly constrained, or if they emerged from a stem lineage encompassing an array of structural variations, based on a more labile "panchelicerate" body plan.

RESULTS

In this study, we reinvestigated the problematic middle Cambrian arthropod Habelia optata Walcott from the Burgess Shale, and found that it was a close relative of Sanctacaris uncata Briggs and Collins (in Habeliida, ord. nov.), both retrieved in our Bayesian phylogeny as stem chelicerates. Habelia possesses an exoskeleton covered in numerous spines and a bipartite telson as long as the rest of the body. Segments are arranged into three tagmata. The prosoma includes a reduced appendage possibly precursor to the chelicera, raptorial endopods connected to five pairs of outstandingly large and overlapping gnathobasic basipods, antennule-like exopods seemingly dissociated from the main limb axis, and, posteriorly, a pair of appendages morphologically similar to thoracic ones. While the head configuration of habeliidans anchors a seven-segmented prosoma as the chelicerate ground pattern, the peculiar size and arrangement of gnathobases and the presence of sensory/tactile appendages also point to an early convergence with the masticatory head of mandibulates.

CONCLUSIONS

Although habeliidans illustrate the early appearance of some diagnostic chelicerate features in the evolution of euarthropods, the unique convergence of their cephalons with mandibulate anatomies suggests that these traits retained an unusual variability in these taxa. The common involvement of strong gnathal appendages across non-megacheirans Cambrian taxa also illustrates that the specialization of the head as the dedicated food-processing tagma was critical to the emergence of both lineages of extant euarthropods-Chelicerata and Mandibulata-and implies that this diversification was facilitated by the expansion of durophagous niches.

The Origin of Animal Multicellularity and Cell Differentiation

Over 600 million years ago, animals evolved from a unicellular or colonial organism whose cell(s) captured bacteria with a collar complex, a flagellum surrounded by a microvillar collar. Using principles from evolutionary cell biology, we reason that the transition to multicellularity required modification of pre-existing mechanisms for extracellular matrix synthesis and cytokinesis. We discuss two hypotheses for the origin of animal cell types: division of labor from ancient plurifunctional cells and conversion of temporally alternating phenotypes into spatially juxtaposed cell types. Mechanistic studies in diverse animals and their relatives promise to deepen our understanding of animal origins and cell biology.

Expression and function of spineless orthologs correlate with distal deutocerebral appendage morphology across Arthropoda

The deutocerebral (second) head segment is putatively homologous across Arthropoda, in spite of remarkable disparity of form and function of deutocerebral appendages. In Mandibulata this segment bears a pair of sensory antennae, whereas in Chelicerata the same segment bears a pair of feeding appendages called chelicerae. Part of the evidence for the homology of deutocerebral appendages is the conserved function of homothorax (hth), which has been shown to specify antennal or cheliceral fate in the absence of Hox signaling, in both mandibulate and chelicerate exemplars. However, the genetic basis for the morphological disparity of antenna and chelicera is not understood. To test whether downstream targets of hth have diverged in a lineage-specific manner, we examined the evolution of the function and expression of spineless (ss), which in two holometabolous insects is known to act as a hth target and distal antennal determinant. Toward expanding phylogenetic representation of gene expression data, here we show that strong expression of ss is observed in developing antennae of a hemimetabolous insect, a centipede, and an amphipod crustacean. By contrast, ss orthologs are not expressed throughout the cheliceral limb buds of spiders or harvestmen during developmental stages when appendage fate is specified. RNA interference-mediated knockdown of ss in Oncopeltus fasciatus, which bears a simple plesiomorphic antenna, resulted in homeotic distal antenna-to-leg transformation, comparable to data from holometabolous insect counterparts. Knockdown of hth in Oncopeltus fasciatus abrogated ss expression, suggesting conservation of upstream regulation. These data suggest that ss may be a flagellar (distal antennal) determinant more broadly, and that this function was acquired at the base of Mandibulata.

Comparative genomics of the tardigrades Hypsibius dujardini and Ramazzottius varieornatus

Tardigrada, a phylum of meiofaunal organisms, have been at the center of discussions of the evolution of Metazoa, the biology of survival in extreme environments, and the role of horizontal gene transfer in animal evolution. Tardigrada are placed as sisters to Arthropoda and Onychophora (velvet worms) in the superphylum Panarthropoda by morphological analyses, but many molecular phylogenies fail to recover this relationship. This tension between molecular and morphological understanding may be very revealing of the mode and patterns of evolution of major groups. Limnoterrestrial tardigrades display extreme cryptobiotic abilities, including anhydrobiosis and cryobiosis, as do bdelloid rotifers, nematodes, and other animals of the water film. These extremophile behaviors challenge understanding of normal, aqueous physiology: how does a multicellular organism avoid lethal cellular collapse in the absence of liquid water? Meiofaunal species have been reported to have elevated levels of horizontal gene transfer (HGT) events, but how important this is in evolution, and particularly in the evolution of extremophile physiology, is unclear. To address these questions, we resequenced and reassembled the genome of H. dujardini, a limnoterrestrial tardigrade that can undergo anhydrobiosis only after extensive pre-exposure to drying conditions, and compared it to the genome of R. varieornatus, a related species with tolerance to rapid desiccation. The 2 species had contrasting gene expression responses to anhydrobiosis, with major transcriptional change in H. dujardini but limited regulation in R. varieornatus. We identified few horizontally transferred genes, but some of these were shown to be involved in entry into anhydrobiosis. Whole-genome molecular phylogenies supported a Tardigrada+Nematoda relationship over Tardigrada+Arthropoda, but rare genomic changes tended to support Tardigrada+Arthropoda.

Segmentation in Tardigrada and diversification of segmental patterns in Panarthropoda

The origin and diversification of segmented metazoan body plans has fascinated biologists for over a century. The superphylum Panarthropoda includes three phyla of segmented animals-Euarthropoda, Onychophora, and Tardigrada. This superphylum includes representatives with relatively simple and representatives with relatively complex segmented body plans. At one extreme of this continuum, euarthropods exhibit an incredible diversity of serially homologous segments. Furthermore, distinct tagmosis patterns are exhibited by different classes of euarthropods. At the other extreme, all tardigrades share a simple segmented body plan that consists of a head and four leg-bearing segments. The modular body plans of panarthropods make them a tractable model for understanding diversification of animal body plans more generally. Here we review results of recent morphological and developmental studies of tardigrade segmentation. These results complement investigations of segmentation processes in other panarthropods and paleontological studies to illuminate the earliest steps in the evolution of panarthropod body plans.