Phylogenetic position of Myriapoda revealed by 454 transcriptome sequencing

The structural and functional modularity of ovarian follicle epithelium in the pill-millipede Hyleoglomeris japonica Verhoeff, 1936 (Diplopoda: Glomerida: Glomeridae)

Ovarian somatic tissues typically surround developing oocytes and play a crucial role in oogenesis across various metazoans, often displaying structural properties specific to their functions. However, there is an absence of evident structural modularity in the follicle epithelium of Myriapoda. We report here two structurally and developmentally distinct domains within the follicle epithelium of the Japanese pill millipede, Hyleoglomeris japonica. The follicle epithelium of H. japonica exhibits a thick cell mass at the apex of the follicle. These cells harbor abundant rough endoplasmic reticulum, mitochondria, Golgi complexes, and numerous microvilli, indicative of synthetic/secretory activities. Moreover, their height increases as oogenesis progresses. In contrast, another region of the epithelium lacks these features. Our findings highlight the presence of structural and functional modularity in the follicle epithelium of H. japonica. We suggest classifying the follicle epithelium of Myriapoda into three types: homogenous epithelia with enhanced synthetic activities, homogenous epithelia with diminished such activities, and heterogeneous epithelia with varying synthetic activities. These findings prompt a reevaluation of the nature of ovarian somatic tissues in Myriapoda as well as in Arthropoda.

Re-evaluating and dating myriapod diversification with phylotranscriptomics under a regime of dense taxon sampling

Recent transcriptomic studies of myriapod phylogeny have been based on relatively small datasets with <40 myriapod terminals and variably supported or contradicted the traditional morphological groupings of Progoneata and Dignatha. Here we amassed a large dataset of 104 myriapod terminals, including multiple species for each of the four myriapod classes. Across the tree, most nodes are stable and well supported. Most analyses across a range of gene occupancy levels provide moderate to strong support for a deep split of Myriapoda into Symphyla + Pauropoda (=Edafopoda) and an uncontradicted grouping of Chilopoda + Diplopoda (=Pectinopoda nov.), as in other recent transcriptome-based analyses; no analysis recovers Progoneata or Dignatha as clades. As in all recent multi-locus and phylogenomic studies, chilopod interrelationships resolve with Craterostigmus excluded from Amalpighiata rather than uniting with other centipedes with maternal brood care in Phylactometria. Diplopod ordinal interrelationships are largely congruent with morphology-based classifications. Chilognathan clades that are not invariably advocated by morphologists include Glomerida + Glomeridesmida, such that the volvation-related characters of pill millipedes may be convergent, and Stemmiulida + Polydesmida more closely allied to Juliformia than to Callipodida + Chordeumatida. The latter relationship implies homoplasy in spinnerets and contradicts Nematophora. A time-tree with nodes calibrated by 25 myriapod and six outgroup fossil terminals recovers Cambrian-Ordovician divergences for the deepest splits in Myriapoda, Edafopoda and Pectinopoda, predating the terrestrial fossil record of myriapods as in other published chronograms, whereas age estimates within Chilopoda and Diplopoda overlap with or do not appreciably predate the calibration fossils. The grouping of Chilopoda and Diplopoda is recovered in all our analyses and is formalized as Pectinopoda nov., named for the shared presence of mandibular comb lamellae. New taxonomic proposals for Chilopoda based on uncontradicted clades are Tykhepoda nov. for the three blind families of Scolopendromorpha that share a "sieve-type" gizzard, and Taktikospina nov. for Scolopendromorpha to the exclusion of Mimopidae.

Was There a Cambrian Explosion on Land? The Case of Arthropod Terrestrialization

Arthropods, the most diverse form of macroscopic life in the history of the Earth, originated in the sea. Since the early Cambrian, at least ~518 million years ago, these animals have dominated the oceans of the world. By the Silurian-Devonian, the fossil record attests to arthropods becoming the first animals to colonize land, However, a growing body of molecular dating and palaeontological evidence suggests that the three major terrestrial arthropod groups (myriapods, hexapods, and arachnids), as well as vascular plants, may have invaded land as early as the Cambrian-Ordovician. These dates precede the oldest fossil evidence of those groups and suggest an unrecorded continental "Cambrian explosion" a hundred million years prior to the formation of early complex terrestrial ecosystems in the Silurian-Devonian. We review the palaeontological, phylogenomic, and molecular clock evidence pertaining to the proposed Cambrian terrestrialization of the arthropods. We argue that despite the challenges posed by incomplete preservation and the scarcity of early Palaeozoic terrestrial deposits, the discrepancy between molecular clock estimates and the fossil record is narrower than is often claimed. We discuss strategies for closing the gap between molecular clock estimates and fossil data in the evolution of early ecosystems on land.

De novo metatranscriptomic exploration of gene function in the millipede holobiont

Invertebrate-microbial associations are widespread in the biosphere and are often related to the function of novel genes, fitness advantages, and even speciation events. Despite ~ 13,000 species of millipedes identified across the world, millipedes and their gut microbiota are markedly understudied compared to other arthropods. Exploring the contribution of individual host-associated microbes is often challenging as many are uncultivable. In this study, we conducted metatranscriptomic profiling of different body segments of a millipede at the holobiont level. This is the first reported transcriptome assembly of a tropical millipede Telodeinopus aoutii (Demange, 1971), as well as the first study on any Myriapoda holobiont. High-throughput RNA sequencing revealed that Telodeinopus aoutii contained > 90% of the core Arthropoda genes. Proteobacteria, Bacteroidetes, Firmicutes, and Euryarchaeota represented dominant and functionally active phyla in the millipede gut, among which 97% of Bacteroidetes and 98% of Firmicutes were present exclusively in the hindgut. A total of 37,831 predicted protein-coding genes of millipede holobiont belonged to six enzyme classes. Around 35% of these proteins were produced by microbiota in the hindgut and 21% by the host in the midgut. Our results indicated that although major metabolic pathways operate at the holobiont level, the involvement of some host and microbial genes are mutually exclusive and microbes predominantly contribute to essential amino acid biosynthesis, short-chain fatty acid metabolism, and fermentation.

Novel mitochondrial gene rearrangements pattern in the millipede Polydesmus sp. GZCS-2019 and phylogenetic analysis of the Myriapoda

The subphylum Myriapoda included four extant classes (Chilopoda, Symphyla, Diplopoda, and Pauropoda). Due to the limitation of taxon sampling, the phylogenetic relationships within Myriapoda remained contentious, especially for Diplopoda. Herein, we determined the complete mitochondrial genome of Polydesmus sp. GZCS-2019 (Myriapoda: Polydesmida) and the mitochondrial genomes are circular molecules of 15,036 bp, with all genes encoded on + strand. The A+T content is 66.1%, making the chain asymmetric, and exhibits negative AT-skew (-0.236). Several genes rearrangements were detected and we propose a new rearrangement model: "TD (N\R) L + C" based on the genome-scale duplication + (non-random/random) loss + recombination. Phylogenetic analyses demonstrated that Chilopoda and Symphyla both were monophyletic group, whereas Pauropoda was embedded in Diplopoda to form the Dignatha. Divergence time showed the first split of Myriapoda occurred between the Chilopoda and other classes (Wenlock period of Silurian). We combine phylogenetic analysis, divergence time, and gene arrangement to yield valuable insights into the evolutionary history and classification relationship of Myriapoda and these results support a monophyletic Progoneata and the relationship (Chilopoda + (Symphyla + (Diplopoda + Pauropoda))) within myriapod. Our results help to better explain the gene rearrangement events of the invertebrate mitogenome and lay the foundation for further phylogenetic study of Myriapoda.

Phylogenetic analyses suggest centipede venom arsenals were repeatedly stocked by horizontal gene transfer

Venoms have evolved over a hundred times in animals. Venom toxins are thought to evolve mostly by recruitment of endogenous proteins with physiological functions. Here we report phylogenetic analyses of venom proteome-annotated venom gland transcriptome data, assisted by genomic analyses, to show that centipede venoms have recruited at least five gene families from bacterial and fungal donors, involving at least eight horizontal gene transfer events. These results establish centipedes as currently the only known animals with venoms used in predation and defence that contain multiple gene families derived from horizontal gene transfer. The results also provide the first evidence for the implication of horizontal gene transfer in the evolutionary origin of venom in an animal lineage. Three of the bacterial gene families encode virulence factors, suggesting that horizontal gene transfer can provide a fast track channel for the evolution of novelty by the exaptation of bacterial weapons into animal venoms.

The tracheal system of scutigeromorph centipedes and the evolution of respiratory systems of myriapods

The tracheal system of scutigeromorph centipedes (Chilopoda) is special, as it consists of dorsally arranged unpaired spiracles. In this study, we investigate the tracheal systems of five different scutigeromorph species. They are strikingly similar to each other but depict unique characters compared to the tracheal systems of pleurostigmophoran centipedes, which has engendered an ongoing debate over a single versus independent origin of tracheal systems in Chilopoda. Up to now, only the respiratory system of Scutigera coleoptrata was investigated intensively using LM-, TEM-, and SEM-techniques. We supplement this with data for species from all three families of Scutigeromorpha. These reveal interspecific differences in atrial width and the shape and branching pattern of the tracheal tubules. Further, we investigated the tracheal system of Scutigera coleoptrata with three additional techniques: light sheet microscopy, microCT and synchrotron radiation based microCT analysis. This set of techniques allows a comparison between fresh versus fixed and dried material. The question of a unique vs. multiple origin of tracheal systems in centipedes and in Myriapoda as a whole is discussed with regard to their structural similarities and differences and the presence of hemocyanin as an oxygen carrier. We used morphological and molecular data and the fossil record to evaluate the alternative hypotheses.

Four myriapod relatives - but who are sisters? No end to debates on relationships among the four major myriapod subgroups

Background

Phylogenetic relationships among the myriapod subgroups Chilopoda, Diplopoda, Symphyla and Pauropoda are still not robustly resolved. The first phylogenomic study covering all subgroups resolved phylogenetic relationships congruently to morphological evidence but is in conflict with most previously published phylogenetic trees based on diverse molecular data. Outgroup choice and long-branch attraction effects were stated as possible explanations for these incongruencies. In this study, we addressed these issues by extending the myriapod and outgroup taxon sampling using transcriptome data.

Results

We generated new transcriptome data of 42 panarthropod species, including all four myriapod subgroups and additional outgroup taxa. Our taxon sampling was complemented by published transcriptome and genome data resulting in a supermatrix covering 59 species. We compiled two data sets, the first with a full coverage of genes per species (292 single-copy protein-coding genes), the second with a less stringent coverage (988 genes). We inferred phylogenetic relationships among myriapods using different data types, tree inference, and quartet computation approaches. Our results unambiguously support monophyletic Mandibulata and Myriapoda. Our analyses clearly showed that there is strong signal for a single unrooted topology, but a sensitivity of the position of the internal root on the choice of outgroups. However, we observe strong evidence for a clade Pauropoda+Symphyla, as well as for a clade Chilopoda+Diplopoda.

Conclusions

Our best quartet topology is incongruent with current morphological phylogenies which were supported in another phylogenomic study. AU tests and quartet mapping reject the quartet topology congruent to trees inferred with morphological characters. Moreover, quartet mapping shows that confounding signal present in the data set is sufficient to explain the weak signal for the quartet topology derived from morphological characters. Although outgroup choice affects results, our study could narrow possible trees to derivatives of a single quartet topology. For highly disputed relationships, we propose to apply a series of tests (AU and quartet mapping), since results of such tests allow to narrow down possible relationships and to rule out confounding signal.

Aquatic stem group myriapods close a gap between molecular divergence dates and the terrestrial fossil record

Identifying marine or freshwater fossils that belong to the stem groups of the major terrestrial arthropod radiations is a longstanding challenge. Molecular dating and fossils of their pancrustacean sister group predict that myriapods originated in the Cambrian, much earlier than their oldest known fossils, but uncertainty about stem group Myriapoda confounds efforts to resolve the timing of the group's terrestrialization. Among a small set of candidates for membership in the stem group of Myriapoda, the Cambrian to Triassic euthycarcinoids have repeatedly been singled out. The only known Devonian euthycarcinoid, Heterocrania rhyniensis from the Rhynie and Windyfield cherts hot spring complex in Scotland, reveals details of head structures that constrain the evolutionary position of euthycarcinoids. The head capsule houses an anterior cuticular tentorium, a feature uniquely shared by myriapods and hexapods. Confocal microscopy recovers myriapod-like characters of the preoral chamber, such as a prominent hypopharynx supported by tentorial bars and superlinguae between the mandibles and hypopharynx, reinforcing an alliance between euthycarcinoids and myriapods recovered in recent phylogenetic analysis. The Cambrian occurrence of the earliest euthycarcinoids supplies the oldest compelling evidence for an aquatic stem group for either Myriapoda or Hexapoda, previously a lacuna in the body fossil record of these otherwise terrestrial lineages until the Silurian and Devonian, respectively. The trace fossil record of euthycarcinoids in the Cambrian and Ordovician reveals amphibious locomotion in tidal environments and fills a gap between molecular estimates for myriapod origins in the Cambrian and a post-Ordovician crown group fossil record.

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.

A new lineage of segmented RNA viruses infecting animals

Metagenomic sequencing has revolutionised our knowledge of virus diversity, with new virus sequences being reported faster than ever before. However, virus discovery from metagenomic sequencing usually depends on detectable homology: without a sufficiently close relative, so-called 'dark' virus sequences remain unrecognisable. An alternative approach is to use virus-identification methods that do not depend on detecting homology, such as virus recognition by host antiviral immunity. For example, virus-derived small RNAs have previously been used to propose 'dark' virus sequences associated with the Drosophilidae (Diptera). Here, we combine published Drosophila data with a comprehensive search of transcriptomic sequences and selected meta-transcriptomic datasets to identify a completely new lineage of segmented positive-sense single-stranded RNA viruses that we provisionally refer to as the Quenyaviruses. Each of the five segments contains a single open reading frame, with most encoding proteins showing no detectable similarity to characterised viruses, and one sharing a small number of residues with the RNA-dependent RNA polymerases of single- and double-stranded RNA viruses. Using these sequences, we identify close relatives in approximately 20 arthropods, including insects, crustaceans, spiders, and a myriapod. Using a more conserved sequence from the putative polymerase, we further identify relatives in meta-transcriptomic datasets from gut, gill, and lung tissues of vertebrates, reflecting infections of vertebrates or of their associated parasites. Our data illustrate the utility of small RNAs to detect viruses with limited sequence conservation, and provide robust evidence for a new deeply divergent and phylogenetically distinct RNA virus lineage.

Enzymatic activities in the digestive tract of spirostreptid and spirobolid millipedes (Diplopoda: Spirostreptida and Spirobolida)

Millipedes represent a model for the study of organic matter transformation, animal-microbial interactions, and compartmentalisation of digestion. The activity of saccharidases (amylase, laminarinase, cellulase, xylanase, chitinase, maltase, cellobiase, and trehalase) and protease were measured in the midgut and hindgut contents and walls of the millipedes Archispirostreptus gigas and Epibolus pulchripes. Assays done at pH 4 and 7 confirmed activities of all enzymes except xylanase. Hydrolysing of starch and laminarin prevailed. The hindgut of E. pulchripes was shorter, less differentiated. Micro-apocrine secretion was observed only in the midgut of A. gigas. Merocrine secretion was present in midgut and hindgut of E. pulchripes, and in the pyloric valve and anterior hindgut of A. gigas. Alpha-polysaccharidases were mostly active in the midgut content and walls, with higher activity at pH 4. The low activity of amylase (A. gigas) and laminarinase (E. pulchripes) in midgut tissue may indicate their synthesis in salivary glands. Cellulases were found in midgut. Chitinases, found in midgut content and tissue (E. pulchripes) or concentrated in the midgut wall (A. gigas), were more active at an acidic pH. Polysaccharidases were low in hindguts. Protease shows midgut origin and alkaline activity extending to the hindgut in E. pulchripes, whereas in A. gigas it is of salivary gland origin and acid activity restricted to the midgut. Some disaccharidases, with more alkaline activity, showed less apparent midgut-hindgut differences. It may indicate an axial separating of the primary and secondary digestion along the intestinal pH gradient or the presence of enzymes of hindgut parasites.

A millipede compound eye mediating low-resolution vision

Millipedes are a species-rich and ancient arthropod clade which typically bear a pair of lateral compound eyes with a small number of large facets. To understand the visual tasks that underlie the evolution of millipede eyes, their spatial resolving performance is of key importance. We here investigate the spatial resolution of the millipede Cylindroiulus punctatus using behavioural assays. Individual animals were placed in the centre of a cylindrical arena under bright downwelling light, with dark stimuli of varying angular dimensions placed on the arena wall. We used continuous isoluminant stimuli based on a difference of Gaussians signal to test for orientation to the dark target via object taxis. Headings of individual animals were tracked in relation to the stimuli to determine whether the animals oriented towards the stimulus. We implemented a multilevel logistic regression model to identify the arc width of the stimulus that animals could resolve. We then modelled the angular sensitivity needed to identify this. We also related the visual performance to the 3D anatomy of the eye. We found that C. punctatus can resolve a stimulus of 56° period (sufficient to detect a 20° dark target). Assuming a contrast threshold of 10%, this requires a receptor acceptance angle of 72° or narrower. Spatial resolving power this low would only suffice for the simplest visual tasks, such as shelter-seeking.

Convergent recruitment of adamalysin-like metalloproteases in the venom of the red bark centipede (Scolopocryptops sexspinosus)

Many venom proteins have presumably been convergently recruited by taxa from diverse venomous lineages. These toxic proteins have characteristics that allow them to remain stable in solution and have a high propensity for toxic effects on prey and/or potential predators. Despite this well-established convergent toxin recruitment, some toxins seem to be lineage specific. To further investigate the toxic proteins found throughout venomous lineages, venom proteomics and venom-gland transcriptomics were performed on two individual red bark centipedes (Scolopocryptops sexspinosus). Combining the protein phenotype with the transcript genotype resulted in the first in-depth venom characterization of S. sexspinosus, including 72 venom components that were identified in both the transcriptome and proteome and 1468 nontoxin transcripts identified in the transcriptome. Ten different toxin families were represented in the venom and venom gland with the majority of the toxins belonging to metalloproteases, CAPS (cysteine-rich secretory protein, antigen 5, and pathogenesis-related 1 proteins), and β-pore-forming toxins. Nine of these toxin families shared a similar proteomic structure to venom proteins previously identified from other centipedes. However, the most highly expressed toxin family, the adamalysin-like metalloproteases, has until now only been observed in the venom of snakes. We confirmed adamalysin-like metalloprotease activity by means of in vivo functional assays. The recruitment of an adamalysin-like metalloprotease into centipede venom represents a striking case of convergent evolution.

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.

Evolutionary dynamics of origin and loss in the deep history of phospholipase D toxin genes

Background

Venom-expressed sphingomyelinase D/phospholipase D (SMase D/PLD) enzymes evolved from the ubiquitous glycerophosphoryl diester phosphodiesterases (GDPD). Expression of GDPD-like SMaseD/PLD toxins in both arachnids and bacteria has inspired consideration of the relative contributions of lateral gene transfer and convergent recruitment in the evolutionary history of this lineage. Previous work recognized two distinct lineages, a SicTox-like (ST-like) clade including the arachnid toxins, and an Actinobacterial-toxin like (AT-like) clade including the bacterial toxins and numerous fungal homologs.

Results

Here we expand taxon sampling by homology detection to discover new GDPD-like SMase D/PLD homologs. The ST-like clade now includes homologs in a wider variety of arthropods along with a sister group in Cnidaria; the AT-like clade now includes additional fungal phyla and proteobacterial homologs; and we report a third clade expressed in diverse aquatic metazoan taxa, a few single-celled eukaryotes, and a few aquatic proteobacteria. GDPD-like SMaseD/PLDs have an ancient presence in chelicerates within the ST-like family and ctenophores within the Aquatic family. A rooted phylogenetic tree shows that the three clades derived from a basal paraphyletic group of proteobacterial GDPD-like SMase D/PLDs, some of which are on mobile genetic elements. GDPD-like SMase D/PLDs share a signature C-terminal motif and a shortened βα1 loop, features that distinguish them from GDPDs. The three major clades also have active site loop signatures that distinguish them from GDPDs and from each other. Analysis of molecular phylogenies with respect to organismal relationships reveals a dynamic evolutionary history including both lateral gene transfer and gene duplication/loss.

Conclusions

The GDPD-like SMaseD/PLD enzymes derive from a single ancient ancestor, likely proteobacterial, and radiated into diverse organismal lineages at least in part through lateral gene transfer.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1302-2) contains supplementary material, which is available to authorized users.

Morphological study of the ovary in Hanseniella caldaria (Myriapoda; Symphyla): The position of oocyte-growth and evolution of ovarian structure in Arthropoda

In Arthropoda, the ovary is classified into Chelicerata-type and Mandibulata-type, based on the oocyte-growth position within the ovary. By contrast, oocytes of Diplopoda and Chilopoda grow within the hemocoelic space. However, as the position of oocyte-growth in Symphyla and Pauropoda has not been confirmed, whether the hemocoelic nature of oocyte-growth is common among myriapods remains ambiguous. This study described the ovarian structure of Hanseniella caldaria to reveal the oocyte-growth position in Symphyla. The oocyte is surrounded by the follicle epithelium, and the inner surface of the follicle epithelium, i.e., the space between follicle cells and oocytes, is lined with a basement membrane. The follicle epithelial layer continues to the ovarian epithelium via the follicle extension with a continuous layer of basement membrane. Data on the architecture of the follicle suggest that the follicle pouch opens to the hemocoel. Hence, the oocyte of H. caldaria grows within the hemocoelic space. Based on our findings in H. caldaria and previous studies in a millipede and in centipedes, the hemocoelic nature of oocyte-growth is considered as a common feature among myriapods and a synapomorphy of the Myriapoda for which morphological synapomorphies have been ambiguous.

Venom-gland transcriptomics and venom proteomics of the giant Florida blue centipede, Scolopendra viridis

The limited number of centipede venom characterizations have revealed a rich diversity of toxins, and recent work has suggested centipede toxins may be more rapidly diversifying than previously considered. Additionally, many identified challenges in venomics research, including assembly and annotation methods, toxin quantification, and the ability to provide biological or technical replicates, have yet to be addressed in centipede venom characterizations. We performed high-throughput, quantifiable transcriptomic and proteomic methods on two individual Scolopendra viridis centipedes from North Florida. We identified 39 toxins that were proteomically confirmed, and 481 nontoxins that were expressed in the venom gland of S. viridis. The most abundant toxins expressed in the venom of S. viridis belonged to calcium and potassium ion-channel toxins, venom allergens, metalloproteases, and β-pore forming toxins. We compared our results to the previously characterized S. viridis from Morelos, Mexico, and found only five proteomically confirmed toxins in common to both localities, suggesting either extreme toxin divergence within S. viridis, or that these populations may represent entirely different species. By using multiple assembly and annotation methods, we generated a comprehensive and quantitative reference transcriptome and proteome of a Scolopendromorpha centipede species, while overcoming some of the challenges present in venomics research.

Diversity, evolution, and function of myriapod hemocyanins

Background

Hemocyanin transports O₂ in the hemolymph of many arthropod species. Such respiratory proteins have long been considered unnecessary in Myriapoda. As a result, the presence of hemocyanin in Myriapoda has long been overlooked. We analyzed transcriptome and genome sequences from all major myriapod taxa – Chilopoda, Diplopoda, Symphyla, and Pauropoda – with the aim of identifying hemocyanin-like proteins.

Results

We investigated the genomes and transcriptomes of 56 myriapod species and identified 46 novel full-length hemocyanin subunit sequences in 20 species of Chilopoda, Diplopoda, and Symphyla, but not Pauropoda. We found in Cleidogona sp. (Diplopoda, Chordeumatida) a hemocyanin-like sequence with mutated copper-binding centers, which cannot bind O₂. An RNA-seq approach showed markedly different hemocyanin mRNA levels from ~ 6 to 25,000 reads per kilobase per million reads. To evaluate the contribution of hemocyanin to O₂ transport, we specifically studied the hemocyanin of the centipede Scolopendra dehaani. This species harbors two distinct hemocyanin subunits with low expression levels. We showed cooperative O₂ binding in the S. dehaani hemolymph, indicating that hemocyanin supports O₂ transport even at low concentration. Further, we demonstrated that hemocyanin is > 1500-fold more highly expressed in the fertilized egg than in the adult.

Conclusion

Hemocyanin was most likely the respiratory protein in the myriapod stem-lineage, but multiple taxa may have independently lost hemocyanin and thus the ability of efficient O₂ transport. In myriapods, hemocyanin is much more widespread than initially appreciated. Some myriapods express hemocyanin only at low levels, which are, nevertheless, sufficient for O₂ supply. Notably, also in myriapods, a non-respiratory protein similar to insect storage hexamerins evolved from the hemocyanin.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1221-2) contains supplementary material, which is available to authorized users.

Comparison of ventral organ development across Pycnogonida (Arthropoda, Chelicerata) provides evidence for a plesiomorphic mode of late neurogenesis in sea spiders and myriapods

BACKGROUND

Comparative studies of neuroanatomy and neurodevelopment provide valuable information for phylogenetic inference. Beyond that, they reveal transformations of neuroanatomical structures during animal evolution and modifications in the developmental processes that have shaped these structures. In the extremely diverse Arthropoda, such comparative studies contribute with ever-increasing structural resolution and taxon coverage to our understanding of nervous system evolution. However, at the neurodevelopmental level, in-depth data remain still largely confined to comparably few laboratory model organisms. Therefore, we studied postembryonic neurogenesis in six species of the bizarre Pycnogonida (sea spiders), which - as the likely sister group of all remaining chelicerates - promise to illuminate neurodevelopmental changes in the chelicerate lineage.

RESULTS

We performed in vivo cell proliferation experiments with the thymidine analogs 5-bromo-2'-deoxyuridine and 5-ethynl-2'-deoxyuridine coupled to fluorescent histochemical staining and immunolabeling, in order to compare ventral nerve cord anatomy and to localize and characterize centers of postembryonic neurogenesis. We report interspecific differences in the architecture of the subesophageal ganglion (SEG) and show the presence of segmental "ventral organs" (VOs) that act as centers of neural cell production during gangliogenesis. These VOs are either incorporated into the ganglionic soma cortex or found on the external ganglion surface. Despite this difference, several shared features support homology of the two VO types, including (1) a specific arrangement of the cells around a small central cavity, (2) the presence of asymmetrically dividing neural stem cell-like precursors, (3) the migration of newborn cells along corresponding pathways into the cortex, and (4) the same VO origin and formation earlier in development.

CONCLUSIONS

Evaluation of our findings relative to current hypotheses on pycnogonid phylogeny resolves a bipartite SEG and internal VOs as plesiomorphic conditions in pycnogonids. Although chelicerate taxa other than Pycnogonida lack comparable VOs, they are a characteristic feature of myriapod gangliogenesis. Accordingly, we propose internal VOs with neurogenic function to be part of the ground pattern of Arthropoda. Further, our findings illustrate the importance of dense sampling in old arthropod lineages - even if as gross-anatomically uniform as Pycnogonida - in order to reliably differentiate plesiomorphic from apomorphic neurodevelopmental characteristics prior to outgroup comparison.

Step-wise evolution of complex chemical defenses in millipedes: a phylogenomic approach

With fossil representatives from the Silurian capable of respiring atmospheric oxygen, millipedes are among the oldest terrestrial animals, and likely the first to acquire diverse and complex chemical defenses against predators. Exploring the origin of complex adaptive traits is critical for understanding the evolution of Earth's biological complexity, and chemical defense evolution serves as an ideal study system. The classic explanation for the evolution of complexity is by gradual increase from simple to complex, passing through intermediate "stepping stone" states. Here we present the first phylogenetic-based study of the evolution of complex chemical defenses in millipedes by generating the largest genomic-based phylogenetic dataset ever assembled for the group. Our phylogenomic results demonstrate that chemical complexity shows a clear pattern of escalation through time. New pathways are added in a stepwise pattern, leading to greater chemical complexity, independently in a number of derived lineages. This complexity gradually increased through time, leading to the advent of three distantly related chemically complex evolutionary lineages, each uniquely characteristic of each of the respective millipede groups.

Phylogenomics illuminates the backbone of the Myriapoda Tree of Life and reconciles morphological and molecular phylogenies

The interrelationships of the four classes of Myriapoda have been an unresolved question in arthropod phylogenetics and an example of conflict between morphology and molecules. Morphology and development provide compelling support for Diplopoda (millipedes) and Pauropoda being closest relatives, and moderate support for Symphyla being more closely related to the diplopod-pauropod group than any of them are to Chilopoda (centipedes). In contrast, several molecular datasets have contradicted the Diplopoda-Pauropoda grouping (named Dignatha), often recovering a Symphyla-Pauropoda group (named Edafopoda). Here we present the first transcriptomic data including a pauropod and both families of symphylans, allowing myriapod interrelationships to be inferred from phylogenomic data from representatives of all main lineages. Phylogenomic analyses consistently recovered Dignatha with strong support. Taxon removal experiments identified outgroup choice as a critical factor affecting myriapod interrelationships. Diversification of millipedes in the Ordovician and centipedes in the Silurian closely approximates fossil evidence whereas the deeper nodes of the myriapod tree date to various depths in the Cambrian-Early Ordovician, roughly coinciding with recent estimates of terrestrialisation in other arthropod lineages, including hexapods and arachnids.

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.

Serotonergic neurons in the ventral nerve cord of Chilopoda - a mandibulate pattern of individually identifiable neurons

BACKGROUND

Given the numerous hypotheses concerning arthropod phylogeny, independent data are needed to supplement knowledge based on traditional external morphology and modern molecular sequence information. One promising approach involves comparisons of the structure and development of the nervous system. Along these lines, the morphology of serotonin-immunoreactive neurons in the ventral nerve cord has been investigated in numerous tetraconate taxa (Crustacea and Hexapoda). It has been shown that these neurons can be identified individually due to their comparably low number, characteristic soma position, and neurite morphology, thus making it possible to establish homologies at the single cell level. Within Chilopoda (centipedes), detailed analyses of major branching patterns of serotonin-immunoreactive neurons are missing, but are crucial for developing meaningful conclusions on the homology of single cells.

RESULTS

In the present study, we re-investigated the distribution and projection patterns of serotonin-immunoreactive neurons in the ventral nerve cord of three centipede species: Scutigera coleoptrata, Lithobius forficatus, and Scolopendra oraniensis. The centipede serotonergic system in the ventral nerve cord contains defined groups of individually identifiable neurons. An anterior and two posterior immunoreactive neurons per hemiganglion with contralateral projections, a pair of ipsilateral projecting lateral neurons (an autapomorphic character for Chilopoda), as well as a postero-lateral group of an unclear number of cells are present in the ground pattern of Chilopoda.

CONCLUSIONS

Comparisons to the patterns of serotonin-immunoreactive neurons of tetraconate taxa support the homology of anterior and posterior neurons. Our results thus support a sister group relationship of Myriapoda and Tetraconata and, further, a mandibulate ground pattern of individually identifiable serotonin-immunoreactive neurons in the ventral nerve cord. Medial neurons are not considered to be part of the tetraconate ground pattern, but could favor the 'Miracrustacea hypothesis', uniting Remipedia, Cephalocarida, and Hexapoda.

Complete mitochondrial genomes of two flat-backed millipedes by next-generation sequencing (Diplopoda, Polydesmida)

A lack of mitochondrial genome data from myriapods is hampering progress across genetic, systematic, phylogenetic and evolutionary studies. Here, the complete mitochondrial genomes of two millipedes, Asiomorphacoarctata Saussure, 1860 (Diplopoda: Polydesmida: Paradoxosomatidae) and Xystodesmus sp. (Diplopoda: Polydesmida: Xystodesmidae) were assembled with high coverage using Illumina sequencing data. The mitochondrial genomes of the two newly sequenced species are circular molecules of 15,644 bp and 15,791 bp, within which the typical mitochondrial genome complement of 13 protein-coding genes, 22 tRNAs and two ribosomal RNA genes could be identified. The mitochondrial genome of Asiomorphacoarctata is the first complete sequence in the family Paradoxosomatidae (Diplopoda: Polydesmida) and the gene order of the two flat-backed millipedes is novel among known myriapod mitochondrial genomes. Unique translocations have occurred, including inversion of one half of the two genomes with respect to other millipede genomes. Inversion of the entire side of a genome (trnF-nad5-trnH-nad4-nad4L, trnP, nad1-trnL2-trnL1-rrnL-trnV-rrnS, trnQ, trnC and trnY) could constitute a common event in the order Polydesmida. Last, our phylogenetic analyses recovered the monophyletic Progoneata, subphylum Myriapoda and four internal classes.

A flexible genetic toolkit for arthropod neurogenesis

Arthropods show considerable variations in early neurogenesis. This includes the pattern of specification, division and movement of neural precursors and progenitors. In all metazoans with nervous systems, including arthropods, conserved genes regulate neurogenesis, which raises the question of how the various morphological mechanisms have emerged and how the same genetic toolkit might generate different morphological outcomes. Here I address this question by comparing neurogenesis across arthropods and show how variations in the regulation and function of the neural genes might explain this phenomenon and how they might have facilitated the evolution of the diverse morphological mechanisms of neurogenesis.

Exploring Phylogenetic Relationships within Myriapoda and the Effects of Matrix Composition and Occupancy on Phylogenomic Reconstruction

Myriapods, including the diverse and familiar centipedes and millipedes, are one of the dominant terrestrial arthropod groups. Although molecular evidence has shown that Myriapoda is monophyletic, its internal phylogeny remains contentious and understudied, especially when compared to those of Chelicerata and Hexapoda. Until now, efforts have focused on taxon sampling (e.g., by including a handful of genes from many species) or on maximizing matrix size (e.g., by including hundreds or thousands of genes in just a few species), but a phylogeny maximizing sampling at both levels remains elusive. In this study, we analyzed 40 Illumina transcriptomes representing 3 of the 4 myriapod classes (Diplopoda, Chilopoda, and Symphyla); 25 transcriptomes were newly sequenced to maximize representation at the ordinal level in Diplopoda and at the family level in Chilopoda. Ten supermatrices were constructed to explore the effect of several potential phylogenetic biases (e.g., rate of evolution, heterotachy) at 3 levels of gene occupancy per taxon (50%, 75%, and 90%). Analyses based on maximum likelihood and Bayesian mixture models retrieved monophyly of each myriapod class, and resulted in 2 alternative phylogenetic positions for Symphyla, as sister group to Diplopoda + Chilopoda, or closer to Diplopoda, the latter hypothesis having been traditionally supported by morphology. Within centipedes, all orders were well supported, but 2 deep nodes remained in conflict in the different analyses despite dense taxon sampling at the family level. Relationships among centipede orders in all analyses conducted with the most complete matrix (90% occupancy) are at odds not only with the sparser but more gene-rich supermatrices (75% and 50% supermatrices) and with the matrices optimizing phylogenetic informativeness or most conserved genes, but also with previous hypotheses based on morphology, development, or other molecular data sets. Our results indicate that a high percentage of ribosomal proteins in the most complete matrices, in conjunction with distance from the root, can act in concert to compromise the estimated relationships within the ingroup. We discuss the implications of these findings in the context of the ever more prevalent quest for completeness in phylogenomic studies.

A molecular palaeobiological exploration of arthropod terrestrialization

Understanding animal terrestrialization, the process through which animals colonized the land, is crucial to clarify extant biodiversity and biological adaptation. Arthropoda (insects, spiders, centipedes and their allies) represent the largest majority of terrestrial biodiversity. Here we implemented a molecular palaeobiological approach, merging molecular and fossil evidence, to elucidate the deepest history of the terrestrial arthropods. We focused on the three independent, Palaeozoic arthropod terrestrialization events (those of Myriapoda, Hexapoda and Arachnida) and showed that a marine route to the colonization of land is the most likely scenario. Molecular clock analyses confirmed an origin for the three terrestrial lineages bracketed between the Cambrian and the Silurian. While molecular divergence times for Arachnida are consistent with the fossil record, Myriapoda are inferred to have colonized land earlier, substantially predating trace or body fossil evidence. An estimated origin of myriapods by the Early Cambrian precedes the appearance of embryophytes and perhaps even terrestrial fungi, raising the possibility that terrestrialization had independent origins in crown-group myriapod lineages, consistent with morphological arguments for convergence in tracheal systems.This article is part of the themed issue 'Dating species divergences using rocks and clocks'.

XX/XY System of Sex Determination in the Geophilomorph Centipede Strigamia maritima

We show that the geophilomorph centipede Strigamia maritima possesses an XX/XY system of sex chromosomes, with males being the heterogametic sex. This is, to our knowledge, the first report of sex chromosomes in any geophilomorph centipede. Using the recently assembled Strigamia genome sequence, we identified a set of scaffolds differentially represented in male and female DNA sequence. Using quantitative real-time PCR, we confirmed that three candidate X chromosome-derived scaffolds are present at approximately twice the copy number in females as in males. Furthermore, we confirmed that six candidate Y chromosome-derived scaffolds contain male-specific sequences. Finally, using this molecular information, we designed an X chromosome-specific DNA probe and performed fluorescent in situ hybridization against mitotic and meiotic chromosome spreads to identify the Strigamia XY sex-chromosome pair cytologically. We found that the X and Y chromosomes are recognizably different in size during the early pachytene stage of meiosis, and exhibit incomplete and delayed pairing.

Genomic insights into the Ixodes scapularis tick vector of Lyme disease

Ticks transmit more pathogens to humans and animals than any other arthropod. We describe the 2.1 Gbp nuclear genome of the tick, Ixodes scapularis (Say), which vectors pathogens that cause Lyme disease, human granulocytic anaplasmosis, babesiosis and other diseases. The large genome reflects accumulation of repetitive DNA, new lineages of retro-transposons, and gene architecture patterns resembling ancient metazoans rather than pancrustaceans. Annotation of scaffolds representing ∼57% of the genome, reveals 20,486 protein-coding genes and expansions of gene families associated with tick-host interactions. We report insights from genome analyses into parasitic processes unique to ticks, including host 'questing', prolonged feeding, cuticle synthesis, blood meal concentration, novel methods of haemoglobin digestion, haem detoxification, vitellogenesis and prolonged off-host survival. We identify proteins associated with the agent of human granulocytic anaplasmosis, an emerging disease, and the encephalitis-causing Langat virus, and a population structure correlated to life-history traits and transmission of the Lyme disease agent.

Phylotranscriptomic Analysis Based on Coalescence was Less Influenced by the Evolving Rates and the Number of Genes: A Case Study in Ericales

Advances in high-throughput sequencing have generated a vast amount of transcriptomic data that are being increasingly used in phylogenetic reconstruction. However, processing the vast datasets for a huge number of genes and even identifying optimal analytical methodology are challenging. Through de novo sequenced and retrieved data from public databases, we identified 221 orthologous protein-coding genes to reconstruct the phylogeny of Ericales, an order characterized by rapid ancient radiation. Seven species representing different families in Ericales were used as in-groups. Both concatenation and coalescence methods yielded the same well-supported topology as previous studies, with only two nodes conflicting with previously reported relationships. The results revealed that a partitioning strategy could improve the traditional concatenation methodology. Rapidly evolving genes negatively affected the concatenation analysis, while slowly evolving genes slightly affected the coalescence analysis. The coalescence methods usually accommodated rate heterogeneity better and required fewer genes to yield well-supported topologies than the concatenation methods with both real and simulated data.

Unlocking the early fossil record of the arthropod central nervous system

Extant panarthropods (euarthropods, onychophorans and tardigrades) are hallmarked by stunning morphological and taxonomic diversity, but their central nervous systems (CNS) are relatively conserved. The timing of divergences of the ground pattern CNS organization of the major panarthropod clades has been poorly constrained because of a scarcity of data from their early fossil record. Although the CNS has been documented in three-dimensional detail in insects from Cenozoic ambers, it is widely assumed that these tissues are too prone to decay to withstand other styles of fossilization or geologically older preservation. However, Cambrian Burgess Shale-type compressions have emerged as sources of fossilized brains and nerve cords. CNS in these Cambrian fossils are preserved as carbon films or as iron oxides/hydroxides after pyrite in association with carbon. Experiments with carcasses compacted in fine-grained sediment depict preservation of neural tissue for a more prolonged temporal window than anticipated by decay experiments in other media. CNS and compound eye characters in exceptionally preserved Cambrian fossils predict divergences of the mandibulate and chelicerate ground patterns by Cambrian Stage 3 (ca 518 Ma), a dating that is compatible with molecular estimates for these splits.

Foreseeing fates: a commentary on Manton (1928) 'On the embryology of a mysid crustacean, Hemimysis lamornae'

Sidnie Manton became best known for her work on arthropod locomotion, and for proposing radical views on the evolution of arthropods that were accepted for a generation. However, her early training was as an embryologist, and the work that she carried out at the beginning of her career still stands as one of the major twentieth century contributions to the study of crustacean embryology. Here, I review her first major paper, largely completed while she was a graduate student, describing embryonic development in Hemimysis lamornae, a small shrimp-like animal found in the seas around the UK. The clarity of her writing and the quality of her figures set a standard that laid the basis for subsequent work, and although not all of her conclusions have stood the test of time, they remain a standard reference for work today. This commentary was written to celebrate the 350th anniversary of the journal Philosophical Transactions of the Royal Society.

Neuropeptide discovery in Symphylella vulgaris (Myriapoda, Symphyla): In silico prediction of the first myriapod peptidome

Arthropods have contributed greatly to our understanding of peptidergic control of physiology and behavior, and being the largest and most diverse animal phylum, represent a model for investigating peptide hormone evolution. Surprisingly, one arthropod subphylum, the Myriapoda, is uninvestigated in terms of its peptide hormones. The public deposition of a transcriptome for Symphylella vulgaris, a pseudocentipede, provides a means for peptide discovery in myriapods. Here, in silico transcriptome mining was used to identify 47 S. vulgaris neuropeptide-encoding transcripts within this dataset. The identified transcripts allowed for the deduction of 31 unique pre/preprohormone sequences, with 97 distinct mature peptides predicted from the deduced proteins. The predicted S. vulgaris peptidome includes members of the adipokinetic hormone/red pigment concentrating hormone, adipokinetic hormone-corazonin-like peptide, allatostatin A, allatostatin C (AST-C), allatotropin, CCHamide, crustacean cardioactive peptide, GSEFLamide, insulin-like peptide, intocin, proctolin, pyrokinin, short neuropeptide F, SIFamide and sulfakinin families. This is the first, and thus far only, peptidome predicted for a myriapod. Of particular note were a modified AST-C, TYWKQCAFNAVSRFamide, that lacks one of two cysteine residues (i.e. one at position 13) stereotypically present in members of this peptide family (and hence is missing the disulfide bridge that spans these residues) and a SIFamide, PPFNGSIFamide, that is truncated due to a lysine for arginine substitution in the dibasic residue pair commonly located at positions 3 and 4 of stereotypical full-length isoforms (e.g. the crustacean peptide GYRKPPFNGSIFamide). The peptides predicted here represent the only extant resource for initiating investigations of native peptidergic signaling in the Myriapoda.

Serotonin-immunoreactivity in the ventral nerve cord of Pycnogonida--support for individually identifiable neurons as ancestral feature of the arthropod nervous system

BACKGROUND

The arthropod ventral nerve cord features a comparably low number of serotonin-immunoreactive neurons, occurring in segmentally repeated arrays. In different crustaceans and hexapods, these neurons have been individually identified and even inter-specifically homologized, based on their soma positions and neurite morphologies. Stereotypic sets of serotonin-immunoreactive neurons are also present in myriapods, whereas in the investigated chelicerates segmental neuron clusters with higher and variable cell numbers have been reported. This led to the suggestion that individually identifiable serotonin-immunoreactive neurons are an apomorphic feature of the Mandibulata. To test the validity of this neurophylogenetic hypothesis, we studied serotonin-immunoreactivity in three species of Pycnogonida (sea spiders). This group of marine arthropods is nowadays most plausibly resolved as sister group to all other extant chelicerates, rendering its investigation crucial for a reliable reconstruction of arthropod nervous system evolution.

RESULTS

In all three investigated pycnogonids, the ventral walking leg ganglia contain different types of serotonin-immunoreactive neurons, the somata of which occurring mostly singly or in pairs within the ganglionic cortex. Several of these neurons are readily and consistently identifiable due to their stereotypic soma position and characteristic neurite morphology. They can be clearly homologized across different ganglia and different specimens as well as across the three species. Based on these homologous neurons, we reconstruct for their last common ancestor (presumably the pycnogonid stem species) a minimal repertoire of at least seven identified serotonin-immunoreactive neurons per hemiganglion. Beyond that, each studied species features specific pattern variations, which include also some neurons that were not reliably labeled in all specimens.

CONCLUSIONS

Our results unequivocally demonstrate the presence of individually identifiable serotonin-immunoreactive neurons in the pycnogonid ventral nerve cord. Accordingly, the validity of this neuroanatomical feature as apomorphy of Mandibulata is questioned and we suggest it to be ancestral for arthropods instead. The pronounced disparities between the segmental pattern in pycnogonids and the one of studied euchelicerates call for denser sampling within the latter taxon. By contrast, overall similarities between the pycnogonid and myriapod patterns may be indicative of single cell homologies in these two taxa. This notion awaits further substantiation from future studies.

Genome of the Rusty Millipede, Trigoniulus corallinus, Illuminates Diplopod, Myriapod, and Arthropod Evolution

The increasing availability of genomic information from the Arthropoda continues to revolutionize our understanding of the biology of this most diverse animal phylum. However, our sampling of arthropod diversity remains uneven, and key clade such as the Myriapoda are severely underrepresented. Here we present the genome of the cosmopolitanly distributed Rusty Millipede Trigoniulus corallinus, which represents the first diplopod genome to be published, and the second example from the Myriapoda as a whole. This genomic resource contains the majority of core eukaryotic genes (94.3%), and key transcription factor classes that were thought to be lost in the Ecdysozoa. Mitochondrial genome and gene family (transcription factor, Dscam, circadian clock-driving protein, odorant receptor cassette, bioactive compound, and cuticular protein) analyses were also carried out to shed light on their states in the Diplopoda and Myriapoda. The ready availability of T. corallinus recommends it as a new model for evolutionary developmental biology, and the data set described here will be of widespread utility in investigating myriapod and arthropod genomics and evolution.

The complete mitochondrial genome of the geophilomorph centipede Strigamia maritima

Strigamia maritima (Myriapoda; Chilopoda) is a species from the soil-living order of geophilomorph centipedes. The Geophilomorpha is the most speciose order of centipedes with over a 1000 species described. They are notable for their large number of appendage bearing segments and are being used as a laboratory model to study the embryological process of segmentation within the myriapods. Using a scaffold derived from the recently published genome of Strigamia maritima that contained multiple mitochondrial protein-coding genes, here we report the complete mitochondrial genome of Strigamia, the first from any geophilomorph centipede. The mitochondrial genome of S. maritima is a circular molecule of 14,938 base pairs, within which we could identify the typical mitochondrial genome complement of 13 protein-coding genes and 2 ribosomal RNA genes. Sequences resembling 16 of the 22 transfer RNA genes typical of metazoan mitochondrial genomes could be identified, many of which have clear deviations from the standard ‘cloverleaf’ secondary structures of tRNA. Phylogenetic trees derived from the concatenated alignment of protein-coding genes of S. maritima and >50 other metazoans were unable to resolve the Myriapoda as monophyletic, but did support a monophyletic group of chilopods: Strigamia was resolved as the sister group of the scolopendromorph Scolopocryptos sp. and these two (Geophilomorpha and Scolopendromorpha), along with the Lithobiomorpha, formed a monophyletic group the Pleurostigmomorpha. Gene order within the S. maritima mitochondrial genome is unique compared to any other arthropod or metazoan mitochondrial genome to which it has been compared. The highly unusual organisation of the mitochondrial genome of Strigamia maritima is in striking contrast with the conservatively evolving nuclear genome: sampling of more members of this order of centipedes will be required to see whether this unusual organization is typical of the Geophilomorpha or results from a more recent reorganisation in the lineage leading to Strigamia.

The first myriapod genome sequence reveals conservative arthropod gene content and genome organisation in the centipede Strigamia maritima

Myriapods (e.g., centipedes and millipedes) display a simple homonomous body plan relative to other arthropods. All members of the class are terrestrial, but they attained terrestriality independently of insects. Myriapoda is the only arthropod class not represented by a sequenced genome. We present an analysis of the genome of the centipede Strigamia maritima. It retains a compact genome that has undergone less gene loss and shuffling than previously sequenced arthropods, and many orthologues of genes conserved from the bilaterian ancestor that have been lost in insects. Our analysis locates many genes in conserved macro-synteny contexts, and many small-scale examples of gene clustering. We describe several examples where S. maritima shows different solutions from insects to similar problems. The insect olfactory receptor gene family is absent from S. maritima, and olfaction in air is likely effected by expansion of other receptor gene families. For some genes S. maritima has evolved paralogues to generate coding sequence diversity, where insects use alternate splicing. This is most striking for the Dscam gene, which in Drosophila generates more than 100,000 alternate splice forms, but in S. maritima is encoded by over 100 paralogues. We see an intriguing linkage between the absence of any known photosensory proteins in a blind organism and the additional absence of canonical circadian clock genes. The phylogenetic position of myriapods allows us to identify where in arthropod phylogeny several particular molecular mechanisms and traits emerged. For example, we conclude that juvenile hormone signalling evolved with the emergence of the exoskeleton in the arthropods and that RR-1 containing cuticle proteins evolved in the lineage leading to Mandibulata. We also identify when various gene expansions and losses occurred. The genome of S. maritima offers us a unique glimpse into the ancestral arthropod genome, while also displaying many adaptations to its specific life history.

Arthropods are the most abundant animals on earth. Among them, insects clearly dominate on land, whereas crustaceans hold the title for the most diverse invertebrates in the oceans. Much is known about the biology of these groups, not least because of genomic studies of the fruit fly Drosophila, the water flea Daphnia, and other species used in research. Here we report the first genome sequence from a species belonging to a lineage that has previously received very little attention—the myriapods. Myriapods were among the first arthropods to invade the land over 400 million years ago, and survive today as the herbivorous millipedes and venomous centipedes, one of which—Strigamia maritima—we have sequenced here. We find that the genome of this centipede retains more characteristics of the presumed arthropod ancestor than other sequenced insect genomes. The genome provides access to many aspects of myriapod biology that have not been studied before, suggesting, for example, that they have diversified receptors for smell that are quite different from those used by insects. In addition, it shows specific consequences of the largely subterranean life of this particular species, which seems to have lost the genes for all known light-sensing molecules, even though it still avoids light.

A transcriptome approach to ecdysozoan phylogeny

The monophyly of Ecdysozoa, which comprise molting phyla, has received strong support from several lines of evidence. However, the internal relationships of Ecdysozoa are still contended. We generated expressed sequence tags from a priapulid (penis worm), a kinorhynch (mud dragon), a tardigrade (water bear) and five chelicerate taxa by 454 transcriptome sequencing. A multigene alignment was assembled from 63 taxa, which comprised after matrix optimization 24,249 amino acid positions with high data density (2.6% gaps, 19.1% missing data). Phylogenetic analyses employing various models support the monophyly of Ecdysozoa. A clade combining Priapulida and Kinorhyncha (i.e. Scalidophora) was recovered as the earliest branch among Ecdysozoa. We conclude that Cycloneuralia, a taxon erected to combine Priapulida, Kinorhyncha and Nematoda (and others), are paraphyletic. Rather Arthropoda (including Onychophora) are allied with Nematoda and Tardigrada. Within Arthropoda, we found strong support for most clades, including monophyletic Mandibulata and Pancrustacea. The phylogeny within the Euchelicerata remained largely unresolved. There is conflicting evidence on the position of tardigrades: While Bayesian and maximum likelihood analyses of only slowly evolving genes recovered Tardigrada as a sister group to Arthropoda, analyses of the full data set, and of subsets containing genes evolving at fast and intermediate rates identified a clade of Tardigrada and Nematoda. Notably, the latter topology is also supported by the analyses of indel patterns.