A congruent solution to arthropod phylogeny: phylogenomics, microRNAs and morphology support monophyletic Mandibulata

Mitochondrial data are not suitable for resolving placental mammal phylogeny

Mitochondrial data have traditionally been used in reconstructing a variety of species phylogenies. The low rates of recombination and thorough characterization of mitochondrial data across vertebrate species make it a particularly attractive phylogenetic marker. The relatively low number of fully sequenced mammal genomes and the lack of extensive sampling within Superorders have posed a serious problem for reaching agreement on the placement mammal species. The use of mitochondrial data sequences from large numbers of mammals could serve to circumvent the taxon-sampling deficit. Here we assess the suitability of mitochondrial data as a phylogenetic marker in mammal phylogenetics. MtDNA datasets of mammal origin have been filtered as follows: (i) we have sampled sparsely across the phylogenetic tree, (ii) we have constrained our sampling to genes with high taxon coverage, (iii) we have categorised rates across sites in a phylogeny independent manner and have removed fast evolving sites, and (iv), we have sampled from very shallow divergence times to reduce phylogenetic conflict. However, topologies obtained using these filters are not consistent with previous studies and are discordant across different genes. Individual mitochondrial genes, and indeed all mitochondrial genes analysed as a supermatrix, resulted in poor resolution of the species phylogeny. Overall, our study highlights the limitations of mitochondrial data, not only for resolving deep divergences and but also for shallow divergences in the mammal phylogeny.

The HR97 (NR1L) group of nuclear receptors: a new group of nuclear receptors discovered in Daphnia species

The recently sequenced Daphnia pulex genome revealed the NR1L nuclear receptor group consisting of three novel receptors. Phylogenetic studies show that this group is related to the NR1I group (CAR/PXR/VDR) and the NR1J group (HR96), and were subsequently named HR97a/b/g. Each of the HR97 paralogs from Daphnia magna, a commonly used crustacean in toxicity testing, was cloned, sequenced, and partially characterized. Phylogenetic analysis indicates that the HR97 receptors are present in primitive arthropods such as the chelicerates but lost in insects. qPCR and immunohistochemistry demonstrate that each of the receptors is expressed near or at reproductive maturity, and that HR97g, the most ancient of the HR97 receptors, is primarily expressed in the gastrointestinal tract, mandibular region, and ovaries, consistent with a role in reproduction. Transactivation assays using an HR97a/b/g-GAL4 chimera indicate that unlike Daphnia HR96 that is promiscuous with respect to ligand recognition, the HR97 receptors do not respond to many of the ligands that activate CAR/PXR/HR96 nuclear receptors. Only three putative ligands of HR97 receptors were identified in this study: pyriproxyfen, methyl farnesoate, and arachidonic acid. Only arachidonic acid, which acts as an inverse agonist, alters HR97g activity at concentrations that would be considered within physiologically relevant ranges. Overall, this study demonstrates that, although closely related to the promiscuous receptors in the NR1I and NR1J groups, the HR97 receptors are mostly likely not multi-xenobiotic sensors, but rather may perform physiological functions, potentially in reproduction, unique to crustaceans and other non-insect arthropod groups.

A critical appraisal of the use of microRNA data in phylogenetics

Recent progress in resolving the tree of life continues to expose relationships that resist resolution, which drives the search for novel sources of information to solve these difficult phylogenetic problems. A recent example, the presence and absence of microRNA families, has been vigorously promoted as an ideal source of phylogenetic data and has been applied to several perennial phylogenetic problems. The utility of such data for phylogenetic inference hinges critically both on developing stochastic models that provide a reasonable description of the process that give rise to these data, and also on the careful validation of those models in real inference scenarios. Remarkably, however, the statistical behavior and phylogenetic utility of microRNA data have not yet been rigorously characterized. Here we explore the behavior and performance of microRNA presence/absence data under a variety of evolutionary models and reexamine datasets from several previous studies. We find that highly heterogeneous rates of microRNA gain and loss, pervasive secondary loss, and sampling error collectively render microRNA-based inference of phylogeny difficult. Moreover, our reanalyses fundamentally alter the conclusions for four of the five studies that we reexamined. Our results indicate that the capacity of miRNA data to resolve the tree of life has been overstated, and we urge caution in their application and interpretation.

Conserved temporal patterns of microRNA expression in Drosophila support a developmental hourglass model

The spatiotemporal control of gene expression is crucial for the successful completion of animal development. The evolutionary constraints on development are particularly strong for the mid-embryonic stage when body segments are specified, as evidenced by a high degree of morphological and protein-coding gene conservation during this period—a phenomenon known as the developmental hourglass. The discovery of microRNA-mediated gene control revealed an entirely new layer of complexity of the molecular networks that orchestrate development. However, the constraints on microRNA developmental expression and evolution, and the implications for animal evolution are less well understood. To systematically explore the conservation of microRNAs during development, we carried out a genome-wide comparative study of microRNA expression levels throughout the ontogenesis of two divergent fruit flies, Drosophila melanogaster and D. virilis. We show that orthologous microRNAs display highly similar temporal profiles regardless of their mutation rates, suggesting that the timely expression of microRNA genes can be more constrained than their sequence. Furthermore, transitions between key developmental events in the different species are accompanied by conserved shifts in microRNA expression profiles, with the mid-embryonic period between gastrulation and segmentation characterized by the highest similarity of microRNA expression. The conservation of microRNA expression therefore displays an hourglass pattern similar to that observed for protein-coding genes.

Germ cells of the centipede Strigamia maritima are specified early in embryonic development

We provide the first systematic description of germ cell development with molecular markers in a myriapod, the centipede Strigamia maritima. By examining the expression of Strigamia vasa and nanos orthologues, we find that the primordial germ cells are specified from at least the blastoderm stage. This is a much earlier embryonic stage than previously described for centipedes, or any other member of the Myriapoda. Using these genes as markers, and taking advantage of the developmental synchrony of Strigamia embryos within single clutches, we are able to track the development of the germ cells throughout embryogenesis. We find that the germ cells accumulate at the blastopore; that the cells do not internalize through the hindgut, but rather through the closing blastopore; and that the cells undergo a long-range migration to the embryonic gonad. This is the first evidence for primordial germ cells displaying these behaviours in any myriapod. The myriapods are a phylogenetically important group in the arthropod radiation for which relatively little developmental data is currently available. Our study provides valuable comparative data that complements the growing number of studies in insects, crustaceans and chelicerates, and is important for the correct reconstruction of ancestral states and a fuller understanding of how germ cell development has evolved in different arthropod lineages.

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.

L.U.St: a tool for approximated maximum likelihood supertree reconstruction

Background

Supertrees combine disparate, partially overlapping trees to generate a synthesis that provides a high level perspective that cannot be attained from the inspection of individual phylogenies. Supertrees can be seen as meta-analytical tools that can be used to make inferences based on results of previous scientific studies. Their meta-analytical application has increased in popularity since it was realised that the power of statistical tests for the study of evolutionary trends critically depends on the use of taxon-dense phylogenies. Further to that, supertrees have found applications in phylogenomics where they are used to combine gene trees and recover species phylogenies based on genome-scale data sets.

Results

Here, we present the L.U.St package, a python tool for approximate maximum likelihood supertree inference and illustrate its application using a genomic data set for the placental mammals. L.U.St allows the calculation of the approximate likelihood of a supertree, given a set of input trees, performs heuristic searches to look for the supertree of highest likelihood, and performs statistical tests of two or more supertrees. To this end, L.U.St implements a winning sites test allowing ranking of a collection of a-priori selected hypotheses, given as a collection of input supertree topologies. It also outputs a file of input-tree-wise likelihood scores that can be used as input to CONSEL for calculation of standard tests of two trees (e.g. Kishino-Hasegawa, Shimidoara-Hasegawa and Approximately Unbiased tests).

Conclusion

This is the first fully parametric implementation of a supertree method, it has clearly understood properties, and provides several advantages over currently available supertree approaches. It is easy to implement and works on any platform that has python installed.

Availability: bitBucket page - https://afro-juju@bitbucket.org/afro-juju/l.u.st.git.

Contact: Davide.Pisani@bristol.ac.uk.

Heterochrony in mandible development of larval shrimp (Decapoda: Caridea)--a comparative morphological SEM study of two carideans

Mandible development in the larval stages I-V of two palaemonid shrimp species, Palaemon elegans and Macrobrachium amazonicum, was analyzed using scanning electron microscopy, light microscopy, and confocal laser scanning microscopy. In contrast to the zoea I of P. elegans, first-stage larvae of M. amazonicum are nonfeeding. At hatching, the morphology of the mandibles is fully expressed in P. elegans, while it appears underdeveloped in M. amazonicum, presenting only small precursors of typical caridean features. In successive zoeal stages, both species show similar developmental changes, but the mandibular characters of the larvae in M. amazonicum were delayed compared to the equivalent stages in P. elegans, especially in the development of submarginal setae and mandible size. In conclusion, our results indicate heterochrony (postdisplacement) of mandible development in M. amazonicum compared to that in P. elegans, which is related to initial lack of mandible functionality or planktivorous feeding at hatching, respectively. This conclusion is supported by comparison with other palaemonid zoeae exhibiting different feeding modes. Our data suggest that an evolutionary ground pattern of mandible morphology is present even in species with nonfeeding first-stage larvae.

Joined at the hip: linked characters and the problem of missing data in studies of disparity

Paleontological investigations into morphological diversity, or disparity, are often confronted with large amounts of missing data. We illustrate how missing discrete data affect disparity using a novel simulation for removing data based on parameters from published datasets that contain both extinct and extant taxa. We develop an algorithm that assesses the distribution of missing characters in extinct taxa, and simulates data loss by applying that distribution to extant taxa. We term this technique "linkage." We compare differences in disparity metrics and ordination spaces produced by linkage and random character removal. When we incorporated linkage among characters, disparity metrics declined and ordination spaces shrank at a slower rate with increasing missing data, indicating that correlations among characters govern the sensitivity of disparity analysis. We also present and test a new disparity method that uses the linkage algorithm to correct for the bias caused by missing data. We equalized proportions of missing data among time bins before calculating disparity, and found that estimates of disparity changed when missing data were taken into account. By removing the bias of missing data, we can gain new insights into the morphological evolution of organisms and highlight the detrimental effects of missing data on disparity analysis.

Connecting ears to eye muscles: evolution of a 'simple' reflex arc

Developmental and evolutionary data from vertebrates are beginning to elucidate the origin of the sensorimotor pathway that links gravity and motion detection to image-stabilizing eye movements--the vestibulo-ocular reflex (VOR). Conserved transcription factors coordinate the development of the vertebrate ear into three functional sensory compartments (graviception/translational linear acceleration, angular acceleration and sound perception). These sensory components connect to specific populations of vestibular and auditory projection neurons in the dorsal hindbrain through undetermined molecular mechanisms. In contrast, a molecular basis for the patterning of the vestibular projection neurons is beginning to emerge. These are organized through the actions of rostrocaudally and dorsoventrally restricted transcription factors into a 'hodological mosaic' within which coherent and largely segregated subgroups are specified to project to different targets in the spinal cord and brain stem. A specific set of these regionally diverse vestibular projection neurons functions as the central element that transforms vestibular sensory signals generated by active and passive head and body movements into motor output through the extraocular muscles. The large dynamic range of motion-related sensory signals requires an organization of VOR pathways as parallel, frequency-tuned, hierarchical connections from the sensory periphery to the motor output. We suggest that eyes, ears and functional connections subserving the VOR are vertebrate novelties that evolved into a functionally coherent motor control system in an almost stereotypic organization across vertebrate taxa.

The oldest known priapulid-like scalidophoran animal and its implications for the early evolution of cycloneuralians and ecdysozoans

Morphological phylogenetic analyses suggest that scalidophorans (priapulids, loriciferans, and kinorhynchs) and nematoids (nematodes and nematomorphs) form the ecdysozoan clade Cycloneuralia, which is a sister group to panarthropods. It has been proposed that extant priapulids and Cambrian priapulid-like scalidophorans, because of their conserved evolution, have the potential to illuminate the ancestral morphology, ecology, and developmental biology of highly derived ecdysozoans such as nematods and arthropods. As such, Cambrian fossils, particularly Markuelia and possibly olivooids, can inform the early evolution of scalidophorans, cycloneuralians, and ecdysozoans. However, the scalidophoran Markuelia is known exclusively as embryo fossils, and the olivooids have been alternatively interpreted as cnidarians or cycloneuralians. Here, we describe a post-embryonic scalidophoran fossil Eopriapulites sphinx new genus and species, which represents the oldest known scalidophoran, from the early Cambrian Period (∼535 Ma) in South China. E. sphinx is similar to modern scalidophorans in having an introvert armed with hollow scalids, a collar with coronal scalids, and a pharynx with pharyngeal teeth, but its scalids and pharyngeal teeth are arranged in a hexaradial pattern. Phylogenetically resolved as a stem-group scalidophoran, E. sphinx shares a hexaradial pattern with the hexaradial arrangement of certain anatomical structures in kinorhynchs, loriciferans, nematoids, and Cambrian fossils such as Eolympia pediculata, which could also be a scalidophoran. Thus, the bodyplan of ancestral cycloneuralians may have had a component of hexaradial symmetry (i.e., some but not necessarily all anatomical parts are hexaradially arranged). If panarthropods are nested within paraphyletic cycloneuralians, as several molecular phylogenetic analyses suggest, the ancestral ecdysozoans may have been a legless worm possibly with a component of hexaradial symmetry.

Arthropod fossil data increase congruence of morphological and molecular phylogenies

The relationships of major arthropod clades have long been contentious, but refinements in molecular phylogenetics underpin an emerging consensus. Nevertheless, molecular phylogenies have recovered topologies that morphological phylogenies have not, including the placement of hexapods within a paraphyletic Crustacea, and an alliance between myriapods and chelicerates. Here we show enhanced congruence between molecular and morphological phylogenies based on 753 morphological characters for 309 fossil and Recent panarthropods. We resolve hexapods within Crustacea, with remipedes as their closest extant relatives, and show that the traditionally close relationship between myriapods and hexapods is an artefact of convergent character acquisition during terrestrialisation. The inclusion of fossil morphology mitigates long-branch artefacts as exemplified by pycnogonids: when fossils are included, they resolve with euchelicerates rather than as a sister taxon to all other euarthropods.

The 'ventral organs' of Pycnogonida (Arthropoda) are neurogenic niches of late embryonic and post-embryonic nervous system development

Early neurogenesis in arthropods has been in the focus of numerous studies, its cellular basis, spatio-temporal dynamics and underlying genetic network being by now comparably well characterized for representatives of chelicerates, myriapods, hexapods and crustaceans. By contrast, neurogenesis during late embryonic and/or post-embryonic development has received less attention, especially in myriapods and chelicerates. Here, we apply (i) immunolabeling, (ii) histology and (iii) scanning electron microscopy to study post-embryonic ventral nerve cord development in Pseudopallene sp., a representative of the sea spiders (Pycnogonida), the presumable sister group of the remaining chelicerates. During early post-embryonic development, large neural stem cells give rise to additional ganglion cell material in segmentally paired invaginations in the ventral ectoderm. These ectodermal cell regions - traditionally designated as 'ventral organs' - detach from the surface into the interior and persist as apical cell clusters on the ventral ganglion side. Each cluster is a post-embryonic neurogenic niche that features a tiny central cavity and initially still houses larger neural stem cells. The cluster stays connected to the underlying ganglionic somata cortex via an anterior and a posterior cell stream. Cell proliferation remains restricted to the cluster and streams, and migration of newly produced cells along the streams seems to account for increasing ganglion cell numbers in the cortex. The pycnogonid cluster-stream-systems show striking similarities to the life-long neurogenic system of decapod crustaceans, and due to their close vicinity to glomerulus-like neuropils, we consider their possible involvement in post-embryonic (perhaps even adult) replenishment of olfactory neurons - as in decapods. An instance of a potentially similar post-embryonic/adult neurogenic system in the arthropod outgroup Onychophora is discussed. Additionally, we document two transient posterior ganglia in the ventral nerve cord of Pseudopallene sp. and evaluate this finding in light of the often discussed reduction of a segmented 'opisthosoma' during pycnogonid evolution.

Phylogenetic position of Myriapoda revealed by 454 transcriptome sequencing

Myriapods had been considered closely allied to hexapods (insects and relatives). However, analyses of molecular sequence data have consistently placed Myriapoda either as a sister group of Pancrustacea, comprising crustaceans and hexapods, and thereby supporting the monophyly of Mandibulata, or retrieved Myriapoda as a sister group of Chelicerata (spiders, ticks, mites and allies). In addition, the relationships among the four myriapod groups (Pauropoda, Symphyla, Diplopoda, Chilopoda) are unclear. To resolve the phylogeny of myriapods and their relationship to other main arthropod groups, we collected transcriptome data from the symphylan Symphylella vulgaris, the centipedes Lithobius forficatus and Scolopendra dehaani, and the millipedes Polyxenus lagurus, Glomeris pustulata and Polydesmus angustus by 454 sequencing. We concatenated a multiple sequence alignment that contained 1550 orthologous single copy genes (1,109,847 amino acid positions) from 55 euarthropod and 14 outgroup taxa. The final selected alignment included 181 genes and 37,425 amino acid positions from 55 taxa, with eight myriapods and 33 other euarthropods. Bayesian analyses robustly recovered monophyletic Mandibulata, Pancrustacea and Myriapoda. Most analyses support a sister group relationship of Symphyla in respect to a clade comprising Chilopoda and Diplopoda. Inclusion of additional sequence data from nine myriapod species resulted in an alignment with poor data density, but broader taxon average. With this dataset we inferred Diplopoda+Pauropoda as closest relatives (i.e., Dignatha) and recovered monophyletic Helminthomorpha. Molecular clock calculations suggest an early Cambrian emergence of Myriapoda ∼513 million years ago and a late Cambrian divergence of myriapod classes. This implies a marine origin of the myriapods and independent terrestrialization events during myriapod evolution.

Genomic sequence and experimental tractability of a new decapod shrimp model, Neocaridina denticulata

The speciose Crustacea is the largest subphylum of arthropods on the planet after the Insecta. To date, however, the only publically available sequenced crustacean genome is that of the water flea, Daphnia pulex, a member of the Branchiopoda. While Daphnia is a well-established ecotoxicological model, previous study showed that one-third of genes contained in its genome are lineage-specific and could not be identified in any other metazoan genomes. To better understand the genomic evolution of crustaceans and arthropods, we have sequenced the genome of a novel shrimp model, Neocaridina denticulata, and tested its experimental malleability. A library of 170-bp nominal fragment size was constructed from DNA of a starved single adult and sequenced using the Illumina HiSeq2000 platform. Core eukaryotic genes, the mitochondrial genome, developmental patterning genes (such as Hox) and microRNA processing pathway genes are all present in this animal, suggesting it has not undergone massive genomic loss. Comparison with the published genome of Daphnia pulex has allowed us to reveal 3750 genes that are indeed specific to the lineage containing malacostracans and branchiopods, rather than Daphnia-specific (E-value: 10⁻⁶). We also show the experimental tractability of N. denticulata, which, together with the genomic resources presented here, make it an ideal model for a wide range of further aquacultural, developmental, ecotoxicological, food safety, genetic, hormonal, physiological and reproductive research, allowing better understanding of the evolution of crustaceans and other arthropods.

Evolutionary relationships of the critically endangered frog Ericabatrachus baleensis Largen, 1991 with notes on incorporating previously unsampled taxa into large-scale phylogenetic analyses

BACKGROUND

The phylogenetic relationships of many taxa remain poorly known because of a lack of appropriate data and/or analyses. Despite substantial recent advances, amphibian phylogeny remains poorly resolved in many instances. The phylogenetic relationships of the Ethiopian endemic monotypic genus Ericabatrachus has been addressed thus far only with phenotypic data and remains contentious.

RESULTS

We obtained fresh samples of the now rare and Critically Endangered Ericabatrachus baleensis and generated DNA sequences for two mitochondrial and four nuclear genes. Analyses of these new data using de novo and constrained-tree phylogenetic reconstructions strongly support a close relationship between Ericabatrachus and Petropedetes, and allow us to reject previously proposed alternative hypotheses of a close relationship with cacosternines or Phrynobatrachus.

CONCLUSIONS

We discuss the implications of our results for the taxonomy, biogeography and conservation of E. baleensis, and suggest a two-tiered approach to the inclusion and analyses of new data in order to assess the phylogenetic relationships of previously unsampled taxa. Such approaches will be important in the future given the increasing availability of relevant mega-alignments and potential framework phylogenies.

Onychophoran Hox genes and the evolution of arthropod Hox gene expression

Introduction

Onychophora is a relatively small phylum within Ecdysozoa, and is considered to be the sister group to Arthropoda. Compared to the arthropods, that have radiated into countless divergent forms, the onychophoran body plan is overall comparably simple and does not display much in-phylum variation. An important component of arthropod morphological diversity consists of variation of tagmosis, i.e. the grouping of segments into functional units (tagmata), and this in turn is correlated with differences in expression patterns of the Hox genes. How these genes are expressed in the simpler onychophorans, the subject of this paper, would therefore be of interest in understanding their subsequent evolution in the arthropods, especially if an argument can be made for the onychophoran system broadly reflecting the ancestral state in the arthropods.

Results

The sequences and embryonic expression patterns of the complete set of ten Hox genes of an onychophoran (Euperipatoides kanangrensis) are described for the first time. We find that they are all expressed in characteristic patterns that suggest a function as classical Hox genes. The onychophoran Hox genes obey spatial colinearity, and with the exception of Ultrabithorax (Ubx), they all have different and distinct anterior expression borders. Notably, Ubx transcripts form a posterior to anterior gradient in the onychophoran trunk. Expression of all onychophoran Hox genes extends continuously from their anterior border to the rear end of the embryo.

Conclusions

The spatial expression pattern of the onychophoran Hox genes may contribute to a combinatorial Hox code that is involved in giving each segment its identity. This patterning of segments in the uniform trunk, however, apparently predates the evolution of distinct segmental differences in external morphology seen in arthropods. The gradient-like expression of Ubx may give posterior segments their specific identity, even though they otherwise express the same set of Hox genes. We suggest that the confined domains of Hox gene expression seen in arthropods evolved from an ancestral onychophoran-like Hox gene pattern. Reconstruction of the ancestral arthropod Hox pattern and comparison with the patterns in the different arthropod classes reveals phylogenetic support for Mandibulata and Tetraconata, but not Myriochelata and Atelocerata.

Production of a reference transcriptome and transcriptomic database (EdwardsiellaBase) for the lined sea anemone, Edwardsiella lineata, a parasitic cnidarian

Background

The lined sea anemone Edwardsiella lineata is an informative model system for evolutionary-developmental studies of parasitism. In this species, it is possible to compare alternate developmental pathways leading from a larva to either a free-living polyp or a vermiform parasite that inhabits the mesoglea of a ctenophore host. Additionally, E. lineata is confamilial with the model cnidarian Nematostella vectensis, providing an opportunity for comparative genomic, molecular and organismal studies.

Description

We generated a reference transcriptome for E. lineata via high-throughput sequencing of RNA isolated from five developmental stages (parasite; parasite-to-larva transition; larva; larva-to-adult transition; adult). The transcriptome comprises 90,440 contigs assembled from >15 billion nucleotides of DNA sequence. Using a molecular clock approach, we estimated the divergence between E. lineata and N. vectensis at 215–364 million years ago. Based on gene ontology and metabolic pathway analyses and gene family surveys (bHLH-PAS, deiodinases, Fox genes, LIM homeodomains, minicollagens, nuclear receptors, Sox genes, and Wnts), the transcriptome of E. lineata is comparable in depth and completeness to N. vectensis. Analyses of protein motifs and revealed extensive conservation between the proteins of these two edwardsiid anemones, although we show the NF-κB protein of E. lineata reflects the ancestral structure, while the NF-κB protein of N. vectensis has undergone a split that separates the DNA-binding domain from the inhibitory domain. All contigs have been deposited in a public database (EdwardsiellaBase), where they may be searched according to contig ID, gene ontology, protein family motif (Pfam), enzyme commission number, and BLAST. The alignment of the raw reads to the contigs can also be visualized via JBrowse.

Conclusions

The transcriptomic data and database described here provide a platform for studying the evolutionary developmental genomics of a derived parasitic life cycle. In addition, these data from E. lineata will aid in the interpretation of evolutionary novelties in gene sequence or structure that have been reported for the model cnidarian N. vectensis (e.g., the split NF-κB locus). Finally, we include custom computational tools to facilitate the annotation of a transcriptome based on high-throughput sequencing data obtained from a “non-model system.”

House dust mite (Der p 10) and crustacean allergic patients may react to food containing Yellow mealworm proteins

SCOPE

Due to the imminent growth of the world population, shortage of protein sources for human consumption will arise in the near future. Alternative and sustainable protein sources (e.g. insects) are being explored for the production of food and feed. In this project, the safety of Yellow mealworms (Tenebrio molitor L.) for human consumption was tested using approaches as advised by the European Food Safety Authority for allergenicity risk assessment.

METHODS AND RESULTS

Different Yellow mealworm protein fractions were prepared, characterised, and tested for cross-reactivity using sera from patients with an inhalation or food allergy to biologically related species (House dust mite (HDM) and crustaceans) by immunoblotting and basophil activation. Furthermore, the stability was investigated using an in vitro pepsin digestion test. IgE from HDM- and crustacean allergic patients cross-reacted with Yellow mealworm proteins. This cross-reactivity was functional, as shown by the induction of basophil activation. The major cross-reactive proteins were identified as tropomyosin and arginine kinase, which are well known allergens in arthropods. These proteins were moderately stable in the pepsin stability test.

CONCLUSION

Based on these cross-reactivity studies, there is a realistic possibility that HDM- and crustacean allergic patients may react to food containing Yellow mealworm proteins.

Subdivision of arthropod cap-n-collar expression domains is restricted to Mandibulata

BACKGROUND

The monophyly of Mandibulata - the division of arthropods uniting pancrustaceans and myriapods - is consistent with several morphological characters, such as the presence of sensory appendages called antennae and the eponymous biting appendage, the mandible. Functional studies have demonstrated that the patterning of the mandible requires the activity of the Hox gene Deformed and the transcription factor cap-n-collar (cnc) in at least two holometabolous insects: the fruit fly Drosophila melanogaster and the beetle Tribolium castaneum. Expression patterns of cnc from two non-holometabolous insects and a millipede have suggested conservation of the labral and mandibular domains within Mandibulata. However, the activity of cnc is unknown in crustaceans and chelicerates, precluding understanding of a complete scenario for the evolution of patterning of this appendage within arthropods. To redress these lacunae, here we investigate the gene expression of the ortholog of cnc in Parhyale hawaiensis, a malacostracan crustacean, and two chelicerates: the harvestman Phalangium opilio, and the scorpion Centruroides sculpturatus.

RESULTS

In the crustacean P. hawaiensis, the segmental expression of Ph-cnc is the same as that reported previously in hexapods and myriapods, with two distinct head domains in the labrum and the mandibular segment. In contrast, Po-cnc and Cs-cnc expression is not enriched in the labrum of either chelicerate, but instead is expressed at comparable levels in all appendages. In further contrast to mandibulate orthologs, the expression domain of Po-cnc posterior to the labrum is not confined within the expression domain of Po-Dfd.

CONCLUSIONS

Expression data from two chelicerate outgroup taxa suggest that the signature two-domain head expression pattern of cnc evolved at the base of Mandibulata. The observation of the archetypal labral and mandibular segment domains in a crustacean exemplar supports the synapomorphic nature of mandibulate cnc expression. The broader expression of Po-cnc with respect to Po-Dfd in chelicerates further suggests that the regulation of cnc by Dfd was also acquired at the base of Mandibulata. To test this hypothesis, future studies examining panarthropod cnc evolution should investigate expression of the cnc ortholog in arthropod outgroups, such as Onychophora and Tardigrada.

Sternites and spiracles - the unclear homology of ventral sclerites in the basal millipede order Glomeridesmida (Myriapoda, Diplopoda)

We report the discovery of a ventral plate in the basal and little-known chilognath millipede order Glomeridesmida. This ventral plate, interpreted here as a 'true sternite', is clearly separate from both the coxa and the more lateral stigma-carrying plates commonly referred to as 'diplopod sternites'. Therefore, the lateral, stigma-carrying plates of the Diplopoda, previously referred to as sternites, are not sternal elements, but subcoxal elements associated with the limb base. This discovery changes the nomenclature used for the ventral plates in Diplopoda, with the formerly named 'sternite' better referred to as 'stigma-carrying plate'. In helminthomorph Diplopoda, the stigma-carrying plates are apparently secondarily fused with the sternite. The main argument for the independent evolution of tracheae in insects and myriapods, the different location of their respective spiracles, no longer holds true. In all Myriapoda and Hexapoda the spiracles associated with subcoxal elements are located lateral to the limb base. This discovery shows that the arguments for an independent origin of tracheae in insects and myriapods are not uncontestable.

Embryonic neurogenesis in Pseudopallene sp. (Arthropoda, Pycnogonida) includes two subsequent phases with similarities to different arthropod groups

BACKGROUND

Studies on early neurogenesis have had considerable impact on the discussion of the phylogenetic relationships of arthropods, having revealed striking similarities and differences between the major lineages. In Hexapoda and crustaceans, neurogenesis involves the neuroblast, a type of neural stem cell. In each hemi-segment, a set of neuroblasts produces neural cells by repeated asymmetrical and interiorly directed divisions. In Euchelicerata and Myriapoda, neurogenesis lacks neural stem cells, featuring instead direct immigration of neural cell groups from fixed sites in the neuroectoderm. Accordingly, neural stem cells were hitherto assumed to be an evolutionary novelty of the Tetraconata (Hexapoda + crustaceans). To further test this hypothesis, we investigated neurogenesis in Pycnogonida, or sea spiders, a group of marine arthropods with close affinities to euchelicerates.

RESULTS

We studied neurogenesis during embryonic development of Pseudopallene sp. (Callipallenidae), using fluorescent histochemical staining and immunolabelling. Embryonic neurogenesis has two phases. The first phase shows notable similarities to euchelicerates and myriapods. These include i) the lack of morphologically different cell types in the neuroectoderm; ii) the formation of transiently identifiable, stereotypically arranged cell internalization sites; iii) immigration of predominantly post-mitotic ganglion cells; and iv) restriction of tangentially oriented cell proliferation to the apical cell layer. However, in the second phase, the formation of a central invagination in each hemi-neuromere is accompanied by the differentiation of apical neural stem cells. The latter grow in size, show high mitotic activity and an asymmetrical division mode. A marked increase of ganglion cell numbers follows their differentiation. Directly basal to the neural stem cells, an additional type of intermediate neural precursor is found.

CONCLUSIONS

Embryonic neurogenesis of Pseudopallene sp. combines features of central nervous system development that have been hitherto described separately in different arthropod taxa. The two-phase character of pycnogonid neurogenesis calls for a thorough reinvestigation of other non-model arthropods over the entire course of neurogenesis. With the currently available data, a common origin of pycnogonid neural stem cells and tetraconate neuroblasts remains unresolved. To acknowledge this, we present two possible scenarios on the evolution of arthropod neurogenesis, whereby Myriapoda play a key role in the resolution of this issue.

A developmental study of serotonin-immunoreactive neurons in the embryonic brain of the marbled crayfish and the migratory locust: evidence for a homologous protocerebral group of neurons

It is well established that the brains of adult malacostracan crustaceans and winged insects display distinct homologies down to the level of single neuropils such as the central complex and the optic neuropils. We wanted to know if developing insect and crustacean brains also share similarities and therefore have explored how neurotransmitter systems arise during arthropod embryogenesis. Previously, Sintoni et al. (2007) had already reported a homology of an individually identified cluster of neurons in the embryonic crayfish and insect brain, the secondary head spot cells that express the Engrailed protein. In the present study, we have documented the ontogeny of the serotonergic system in embryonic brains of the Marbled Crayfish in comparison to Migratory Locust embryos using immunohistochemical methods combined with confocal laser-scan microscopy. In both species, we found a cluster of early emerging serotonin-immunoreactive neurons in the protocerebrum with neurites that cross to the contralateral brain hemisphere in a characteristic commissure suggesting a homology of this cell cluster. Our study is a first step towards a phylogenetic analysis of neurotransmitter system development and shows that, as for the ventral nerve cord, traits related to neurogenesis in the brain can provide valuable hints for resolving the much debated question of arthropod phylogeny.

Molecular phylogenetic analyses support the monophyly of Hexapoda and suggest the paraphyly of Entognatha

BACKGROUND

Molecular phylogenetic analyses have revealed that Hexapoda and Crustacea form a common clade (the Pancrustacea), which is now widely accepted among zoologists; however, the origin of Hexapoda remains unresolved. The main problems are the unclear relationships among the basal hexapod lineages, Protura (proturans), Collembola (springtails), Diplura (diplurans), and Ectognatha (bristletails, silverfishes, and all winged insects). Mitogenomic analyses have challenged hexapod monophyly and suggested the reciprocal paraphyly of Hexapoda and Crustacea, whereas studies based on nuclear molecular data support the monophyletic origin of hexapods. Additionally, there are significant discrepancies with respect to these issues between the results of morphological and molecular studies. To investigate these problems, we performed phylogenetic analyses of Pancrustacea based on the protein sequences of three orthologous nuclear genes encoding the catalytic subunit of DNA polymerase delta and the largest and second largest subunits of RNA polymerase II from 64 species of arthropods, including representatives of all hexapod orders.

RESULTS

Phylogenetic analyses were conducted based on the inferred amino acid (aa) sequences (~3400 aa in total) of the three genes using the maximum likelihood (ML) method and Bayesian inference. Analyses were also performed with additional datasets generated by excluding long-branch taxa or by using different outgroups. These analyses all yielded essentially the same results. All hexapods were clustered into a common clade, with Branchiopoda as its sister lineage, whereas Crustacea was paraphyletic. Within Hexapoda, the lineages Ectognatha, Palaeoptera, Neoptera, Polyneoptera, and Holometabola were each confirmed to be monophyletic with robust support, but monophyly was not supported for Entognatha (Protura + Collembola + Diplura), Ellipura (Protura + Collembola), or Nonoculata (Protura + Diplura). Instead, our results showed that Protura is the sister lineage to all other hexapods and that Diplura or Diplura + Collembola is closely related to Ectognatha.

CONCLUSION

This is the first study to include all hexapod orders in a phylogenetic analysis using multiple nuclear protein-coding genes to investigate the phylogeny of Hexapoda, with an emphasis on Entognatha. The results strongly support the monophyletic origin of hexapods but reject the monophyly of Entognatha, Ellipura, and Nonoculata. Our results provided the first molecular evidence in support of Protura as the sister group to other hexapods. These findings are expected to provide additional insights into the origin of hexapods and the processes involved in the adaptation of insects to life on land.

The expression of limb gap genes in the mite Archegozetes longisetosus reveals differential patterning mechanisms in chelicerates

The modular organization of arthropod limbs has lead to the evolution of a diversity of appendages within this phylum. A conserved trait within the arthropods is the utilization of a conserved set of regulatory genes that specify the appendage podomeres along the proximo-distal axis, termed the limb gap genes. These include extradenticle, homothorax, dachshund, and Distal-less. The deployment of these genes in the most basally branching arthropod group, the chelicerates, has only been studied in detail in two chelicerate groups, the harvestmen and spiders. Given the broad range of appendage diversity within the chelicerates, comparative studies of gap gene deployment in other chelicerates groups is needed. We therefore followed limb gap gene expression in a member of the largest chelicerate group, Acari, the oribatid mite Archegozetes longisetosus. We show that in contrast to many arthropod species, A. longisetosus expresses homothorax and extradenticle exclusively in the proximal portion of the appendages, which refutes the hypothesis of a sister-group relationship between chelicerates and myriapods. We also provide evidence that mites posses the ancestral chelicerate condition of possessing three-segmented chelicerae, which also express the gene dachshund. This adds support to the hypothesis that a cheliceral dachshund domain is ancestral to arachnids. Lastly, we provide evidence that the suppression of the fourth pair of walking legs, a putative synapomorphy for Acari, is accomplished by repressing the development of the medial and distal regions of the limb.

Rates of phenotypic and genomic evolution during the Cambrian explosion

The near-simultaneous appearance of most modern animal body plans (phyla) ~530 million years ago during the Cambrian explosion is strong evidence for a brief interval of rapid phenotypic and genetic innovation, yet the exact speed and nature of this grand adaptive radiation remain debated. Crucially, rates of morphological evolution in the past (i.e., in ancestral lineages) can be inferred from phenotypic differences among living organisms-just as molecular evolutionary rates in ancestral lineages can be inferred from genetic divergences. We here employed Bayesian and maximum likelihood phylogenetic clock methods on an extensive anatomical and genomic data set for arthropods, the most diverse phylum in the Cambrian and today. Assuming an Ediacaran origin for arthropods, phenotypic evolution was ~4 times faster, and molecular evolution ~5.5 times faster, during the Cambrian explosion compared to all subsequent parts of the Phanerozoic. These rapid evolutionary rates are robust to assumptions about the precise age of arthropods. Surprisingly, these fast early rates do not change substantially even if the radiation of arthropods is compressed entirely into the Cambrian (~542 mega-annum [Ma]) or telescoped into the Cryogenian (~650 Ma). The fastest inferred rates are still consistent with evolution by natural selection and with data from living organisms, potentially resolving "Darwin's dilemma." However, evolution during the Cambrian explosion was unusual (compared to the subsequent Phanerozoic) in that fast rates were present across many lineages.

Substantial loss of conserved and gain of novel MicroRNA families in flatworms

Recent studies on microRNA (miRNA) evolution focused mainly on the comparison of miRNA complements between animal clades. However, evolution of miRNAs within such groups is poorly explored despite the availability of comparable data that in some cases lack only a few key taxa. For flatworms (Platyhelminthes), miRNA complements are available for some free-living flatworms and all major parasitic lineages, except for the Monogenea. We present the miRNA complement of the monogenean flatworm Gyrodactylus salaris that facilitates a comprehensive analysis of miRNA evolution in Platyhelminthes. Using the newly designed bioinformatics pipeline miRCandRef, the miRNA complement was disentangled from next-generation sequencing of small RNAs and genomic DNA without a priori genome assembly. It consists of 39 miRNA hairpin loci of conserved miRNA families, and 22 novel miRNAs. A comparison with the miRNA complements of Schmidtea mediterranea (Turbellaria), Schistosoma japonicum (Trematoda), and Echinococcus granulosus (Cestoda) reveals a substantial loss of conserved bilaterian, protostomian, and lophotrochozoan miRNAs. Eight of the 46 expected conserved miRNAs were lost in all flatworms, 16 in Neodermata and 24 conserved miRNAs could not be detected in the cestode and the trematode. Such a gradual loss of miRNAs has not been reported before for other animal phyla. Currently, little is known about miRNAs in Platyhelminthes, and for the majority of the lost miRNAs there is no prediction of function. As suggested earlier they might be related to morphological simplifications. The presence and absence of 153 conserved miRNAs was compared for platyhelminths and 32 other metazoan taxa. Phylogenetic analyses support the monophyly of Platyhelminthes (Turbellaria + Neodermata [Monogenea {Trematoda + Cestoda}]).

miRNAs: small genes with big potential in metazoan phylogenetics

microRNAs (miRNAs) are a key component of gene regulatory networks and have been implicated in the regulation of virtually every biological process found in multicellular eukaryotes. What makes them interesting from a phylogenetic perspective is the high conservation of primary sequence between taxa, their accrual in metazoan genomes through evolutionary time, and the rarity of secondary loss in most metazoan taxa. Despite these properties, the use of miRNAs as phylogenetic markers has not yet been discussed within a clear conceptual framework. Here we highlight five properties of miRNAs that underlie their utility in phylogenetics: 1) The processes of miRNA biogenesis enable the identification of novel miRNAs without prior knowledge of sequence; 2) The continuous addition of miRNA families to metazoan genomes through evolutionary time; 3) The low level of secondary gene loss in most metazoan taxa; 4) The low substitution rate in the mature miRNA sequence; and 5) The small probability of convergent evolution of two miRNAs. Phylogenetic analyses using both Bayesian and parsimony methods on a eumetazoan miRNA data set highlight the potential of miRNAs to become an invaluable new tool, especially when used as an additional line of evidence, to resolve previously intractable nodes within the tree of life.

Arthropod phylogenetics in light of three novel millipede (myriapoda: diplopoda) mitochondrial genomes with comments on the appropriateness of mitochondrial genome sequence data for inferring deep level relationships

BACKGROUND

Arthropods are the most diverse group of eukaryotic organisms, but their phylogenetic relationships are poorly understood. Herein, we describe three mitochondrial genomes representing orders of millipedes for which complete genomes had not been characterized. Newly sequenced genomes are combined with existing data to characterize the protein coding regions of myriapods and to attempt to reconstruct the evolutionary relationships within the Myriapoda and Arthropoda.

RESULTS

The newly sequenced genomes are similar to previously characterized millipede sequences in terms of synteny and length. Unique translocations occurred within the newly sequenced taxa, including one half of the Appalachioria falcifera genome, which is inverted with respect to other millipede genomes. Across myriapods, amino acid conservation levels are highly dependent on the gene region. Additionally, individual loci varied in the level of amino acid conservation. Overall, most gene regions showed low levels of conservation at many sites. Attempts to reconstruct the evolutionary relationships suffered from questionable relationships and low support values. Analyses of phylogenetic informativeness show the lack of signal deep in the trees (i.e., genes evolve too quickly). As a result, the myriapod tree resembles previously published results but lacks convincing support, and, within the arthropod tree, well established groups were recovered as polyphyletic.

CONCLUSIONS

The novel genome sequences described herein provide useful genomic information concerning millipede groups that had not been investigated. Taken together with existing sequences, the variety of compositions and evolution of myriapod mitochondrial genomes are shown to be more complex than previously thought. Unfortunately, the use of mitochondrial protein-coding regions in deep arthropod phylogenetics appears problematic, a result consistent with previously published studies. Lack of phylogenetic signal renders the resulting tree topologies as suspect. As such, these data are likely inappropriate for investigating such ancient relationships.

How are comparative genomics and the study of microRNAs changing our views on arthropod endocrinology and adaptations to the environment?

As the last few decades of work has shown, precise regulation of biosynthesis and release of arthropod hormones is essential to cope with environmental stresses and challenges. In crustaceans and insects, the sesquiterpenoids methyl farnesoate (MF), farnesoic acid (FA) and juvenile hormone (JH) regulate many developmental, physiological, and reproductive processes. In this review, we discuss how comparative genomics has and will impact our views on arthropod endocrinology. We will also highlight the current knowledge of regulation of genes involved in arthropod hormone biosynthesis by microRNAs, and describe the potential insights into arthropod endocrinology, evolution, and adaptation that are likely to come from the study of microRNAs.

MicroRNAs support the monophyly of enteropneust hemichordates

Understanding the evolutionary history of deuterostomes requires elucidating the phylogenetic interrelationships amongst the constituent taxa. Although the monophyly and interrelationships among the three principal groups-the chordates, the echinoderms, and the hemichordates-are well established, as are the internal relationships among the echinoderm and chordate taxa, the interrelationships among the principal groups of hemichordates-the harrimaniid enteropneusts, the ptychoderid enteropneusts, and the pterobranchs-remain unresolved. Depending on the study some find enteropneusts paraphyletic with pterobranchs (e.g., Cephalodiscus) more closely related to the harrimaniid enteropneusts (e.g., Saccoglossus) than either are to the ptychoderid enteropneusts (e.g., Ptychodera), whereas other studies support a monophyletic Enteropneusta. To try and resolve between these two competing hypotheses, we turned to microRNAs, small ∼22 nt non-coding RNA genes that have been shown to shed insight into particularly difficult phylogenetic questions. Using deep sequencing we characterized the small RNA repertoires of two hemichordate species, Cephalodiscus hodgsoni and Ptychodera flava, and the crinoid echinoderm Antedon mediterranea, and combined our results with the described complements of the hemichordate Saccoglossus kowalevskii, the sea urchin Strongylocentrotus purpuratus, and the starfish Patiria miniata. Our data unambiguously support the monophyly of Enteropneusts as S. kowalevskii shares 12 miRNA sequences with P. flava that are not present in the C. hodgsoni or A. mediterranea libraries, and have never been reported from another metazoan taxon. Thus, these data resolve the phylogenetic position of pterobranchs, ultimately allowing for a better understanding of body plan evolution throughout the deuterostomes.

aLeaves facilitates on-demand exploration of metazoan gene family trees on MAFFT sequence alignment server with enhanced interactivity

We report a new web server, aLeaves (http://aleaves.cdb.riken.jp/), for homologue collection from diverse animal genomes. In molecular comparative studies involving multiple species, orthology identification is the basis on which most subsequent biological analyses rely. It can be achieved most accurately by explicit phylogenetic inference. More and more species are subjected to large-scale sequencing, but the resultant resources are scattered in independent project-based, and multi-species, but separate, web sites. This complicates data access and is becoming a serious barrier to the comprehensiveness of molecular phylogenetic analysis. aLeaves, launched to overcome this difficulty, collects sequences similar to an input query sequence from various data sources. The collected sequences can be passed on to the MAFFT sequence alignment server (http://mafft.cbrc.jp/alignment/server/), which has been significantly improved in interactivity. This update enables to switch between (i) sequence selection using the Archaeopteryx tree viewer, (ii) multiple sequence alignment and (iii) tree inference. This can be performed as a loop until one reaches a sensible data set, which minimizes redundancy for better visibility and handling in phylogenetic inference while covering relevant taxa. The work flow achieved by the seamless link between aLeaves and MAFFT provides a convenient online platform to address various questions in zoology and evolutionary biology.

The impact of paralogy on phylogenomic studies - a case study on annelid relationships

Phylogenomic studies based on hundreds of genes derived from expressed sequence tags libraries are increasingly used to reveal the phylogeny of taxa. A prerequisite for these studies is the assignment of genes into clusters of orthologous sequences. Sophisticated methods of orthology prediction are used in such analyses, but it is rarely assessed whether paralogous sequences have been erroneously grouped together as orthologous sequences after the prediction, and whether this had an impact on the phylogenetic reconstruction using a super-matrix approach. Herein, I tested the impact of paralogous sequences on the reconstruction of annelid relationships based on phylogenomic datasets. Using single-partition analyses, screening for bootstrap support, blast searches and pruning of sequences in the supermatrix, wrongly assigned paralogous sequences were found in eight partitions and the placement of five taxa (the annelids Owenia, Scoloplos, Sthenelais and Eurythoe and the nemertean Cerebratulus) including the robust bootstrap support could be attributed to the presence of paralogous sequences in two partitions. Excluding these sequences resulted in a different, weaker supported placement for these taxa. Moreover, the analyses revealed that paralogous sequences impacted the reconstruction when only a single taxon represented a previously supported higher taxon such as a polychaete family. One possibility of a priori detection of wrongly assigned paralogous sequences could combine 1) a screening of single-partition analyses based on criteria such as nodal support or internal branch length with 2) blast searches of suspicious cases as presented herein. Also possible are a posteriori approaches in which support for specific clades is investigated by comparing alternative hypotheses based on differences in per-site likelihoods. Increasing the sizes of EST libraries will also decrease the likelihood of wrongly assigned paralogous sequences, and in the case of orthology prediction methods like HaMStR it is likewise decreased by using more than one reference taxon.

Molecular timetrees reveal a Cambrian colonization of land and a new scenario for ecdysozoan evolution

Ecdysozoans have been key components of ecosystems since the early Cambrian, when trilobites and soft-bodied Burgess Shale-type ecdysozoans dominated marine animal communities. Even today, the most abundant animals on Earth are either nematode worms or plankton-forming crustaceans, whereas the most diverse are the insects. Throughout geological time, several ecdysozoan lineages independently colonized land, shaping both marine and terrestrial ecosystems and providing an adequate environment for successive animal terrestrialization. The timing of these events is largely uncertain and has been investigated only partially using molecular data. Here we present a timescale of ecdysozoan evolution based on multiple molecular data sets, the most complete set of fossil calibrations to date, and a thorough series of validation analyses. Results converge on an Ediacaran origin of all major ecdysozoan lineages (∼587-543 million years ago [mya]), followed by a fast Cambrian radiation of the pancrustaceans (∼539-511 mya), a Cambro-Ordovician colonization of land of different arthropod lineages (∼510-471 mya), and a relatively recent radiation of extant nematodes, onychophorans, and tardigrades (∼442 mya). Arthropods colonized land nearly synchronously with land plants. Further diversification within flying insects, nematodes and onychophorans might be related to the evolution of vascular plants and forests.

Cambrian bivalved arthropod reveals origin of arthrodization

Extant arthropods are diverse and ubiquitous, forming a major constituent of most modern ecosystems. Evidence from early Palaeozoic Konservat Lagerstätten indicates that this has been the case since the Cambrian. Despite this, the details of arthropod origins remain obscure, although most hypotheses regard the first arthropods as benthic predators or scavengers such as the fuxianhuiids or megacheirans ('great-appendage' arthropods). Here, we describe a new arthropod from the Tulip Beds locality of the Burgess Shale Formation (Cambrian, series 3, stage 5) that possesses a weakly sclerotized thorax with filamentous appendages, encased in a bivalved carapace, and a strongly sclerotized, elongate abdomen and telson. A cladistic analysis resolved this taxon as the basal-most member of a paraphyletic grade of nekto-benthic forms with bivalved carapaces. This grade occurs at the base of Arthropoda (panarthropods with arthropodized trunk limbs) and suggests that arthrodization (sclerotization and jointing of the exoskeleton) evolved to facilitate swimming. Predatory and fully benthic habits evolved later in the euarthropod stem-lineage and are plesiomorphically retained in pycnogonids (sea spiders) and euchelicerates (horseshoe crabs and arachnids).

Resolving discrepancy between nucleotides and amino acids in deep-level arthropod phylogenomics: differentiating serine codons in 21-amino-acid models

BACKGROUND

In a previous study of higher-level arthropod phylogeny, analyses of nucleotide sequences from 62 protein-coding nuclear genes for 80 panarthopod species yielded significantly higher bootstrap support for selected nodes than did amino acids. This study investigates the cause of that discrepancy.

METHODOLOGY/PRINCIPAL FINDINGS

The hypothesis is tested that failure to distinguish the serine residues encoded by two disjunct clusters of codons (TCN, AGY) in amino acid analyses leads to this discrepancy. In one test, the two clusters of serine codons (Ser1, Ser2) are conceptually translated as separate amino acids. Analysis of the resulting 21-amino-acid data matrix shows striking increases in bootstrap support, in some cases matching that in nucleotide analyses. In a second approach, nucleotide and 20-amino-acid data sets are artificially altered through targeted deletions, modifications, and replacements, revealing the pivotal contributions of distinct Ser1 and Ser2 codons. We confirm that previous methods of coding nonsynonymous nucleotide change are robust and computationally efficient by introducing two new degeneracy coding methods. We demonstrate for degeneracy coding that neither compositional heterogeneity at the level of nucleotides nor codon usage bias between Ser1 and Ser2 clusters of codons (or their separately coded amino acids) is a major source of non-phylogenetic signal.

CONCLUSIONS

The incongruity in support between amino-acid and nucleotide analyses of the forementioned arthropod data set is resolved by showing that "standard" 20-amino-acid analyses yield lower node support specifically when serine provides crucial signal. Separate coding of Ser1 and Ser2 residues yields support commensurate with that found by degenerated nucleotides, without introducing phylogenetic artifacts. While exclusion of all serine data leads to reduced support for serine-sensitive nodes, these nodes are still recovered in the ML topology, indicating that the enhanced signal from Ser1 and Ser2 is not qualitatively different from that of the other amino acids.

The Cambrian conundrum: early divergence and later ecological success in the early history of animals

Diverse bilaterian clades emerged apparently within a few million years during the early Cambrian, and various environmental, developmental, and ecological causes have been proposed to explain this abrupt appearance. A compilation of the patterns of fossil and molecular diversification, comparative developmental data, and information on ecological feeding strategies indicate that the major animal clades diverged many tens of millions of years before their first appearance in the fossil record, demonstrating a macroevolutionary lag between the establishment of their developmental toolkits during the Cryogenian [(850 to 635 million years ago (Ma)], and the later ecological success of metazoans during the Ediacaran (635 to 541 Ma) and Cambrian (541 to 488 Ma) periods. We argue that this diversification involved new forms of developmental regulation, as well as innovations in networks of ecological interaction within the context of permissive environmental circumstances.

Cap'n'collar differentiates the mandible from the maxilla in the beetle Tribolium castaneum

BACKGROUND

The biting mandible of the arthropods is thought to have evolved in the ancestor of the insects, crustaceans and myriapods: the Mandibulata. A unique origin suggests a common set of developmental genes will be required to pattern the mandible in different arthropods. To date we have functional studies on patterning of the mandibular segment of Drosophila melanogaster showing in particular the effects of the gene cap'n'collar (cnc), however, the dipteran head is far from representative of insects or of more distantly related mandibulates; Drosophila does not even possess a mandibular appendage. To study the development of a more representative insect mandible, we chose the red flour beetle Tribolium castaneum and investigated the function of the Tribolium orthologs of cap'n'collar (Tc-cnc) and the Hox gene Deformed (Tc-Dfd). In order to determine the function of Tc-cnc and Tc-Dfd, transcripts were knocked down by maternal RNA interference (RNAi). The effects of gene knockdown were examined in the developing embryos and larvae. The effect of Tc-cnc and Tc-Dfd knockdown on the expression of other genes was determined by using in situ hybridization on Tribolium embryos.

RESULTS

Our analyses show that Tc-cnc is required for specification of the identity of the mandibular segment of Tribolium and differentiates the mandible from maxillary identity. Loss of Tc-cnc function results in a transformation of the mandible to maxillary identity as well as deletion of the labrum. Tc-Dfd and the Tribolium homolog of proboscipedia (Tc-mxp = maxillopedia), Hox genes that are required to pattern the maxillary appendage, are expressed in a maxilla-like manner in the transformed mandible. Tribolium homologs of paired (Tc-prd) and Distal-less (Tc-Dll) that are expressed in the endites and telopodites of embryonic appendages are also expressed in a maxilla-like manner in the transformed mandible.We also show that Tc-Dfd is required to activate the collar of Tc-cnc expression in the mandibular segment but not the cap expression in the labrum. Tc-Dfd is also required for the activation of Tc-prd in the endites of the mandible and maxillary appendages.

CONCLUSIONS

Tc-cnc is necessary for patterning the mandibular segment of Tribolium. Together, Tc-cnc and Tc-Dfd cooperate to specify mandibular identity, as in Drosophila. Expression patterns of the homologs of cnc and Dfd are conserved in mandibulate arthropods suggesting that the mandible specifying function of cnc is likely to be conserved across the mandibulate arthropods.