A 454 sequencing approach for large scale phylogenomic analysis of the common emperor scorpion (Pandinus imperator)

The separation between mRNA-ends is more variable than expected

Effective circularization of mRNA molecules is a key step for the efficient initiation of translation. Research has shown that the intrinsic separation of the ends of mRNA molecules is rather small, suggesting that intramolecular arrangements could provide this effective circularization. Considering that the innate proximity of RNA ends might have important unknown biological implications, we aimed to determine whether the close proximity of the ends of mRNA molecules is a conserved feature across organisms and gain further insights into the functional effects of the proximity of RNA ends. To do so, we studied the secondary structure of 274 full native mRNA molecules from 17 different organisms to calculate the contour length (CL) of the external loop as an index of their end-to-end separation. Our computational predictions show bigger variations (from 0.59 to 31.8 nm) than previously reported and also than those observed in random sequences. Our results suggest that separations larger than 18.5 nm are not favored, whereas short separations could be related to phenotypical stability. Overall, our work implies the existence of a biological mechanism responsible for the increase in the observed variability, suggesting that the CL features of the exterior loop could be relevant for the initiation of translation and that a short CL could contribute to the stability of phenotypes.

Function and evolution of the aquaporin IsAQP1 in the Lyme disease vector Ixodes scapularis

Ticks are important vectors of pathogenic viruses, bacteria, and protozoans to humans, wildlife, and domestic animals. Due to their life cycles, ticks face significant challenges related to water homeostasis. When blood-feeding, they must excrete water and ions, but when off-host (for stretches lasting several months), they must conserve water to avoid desiccation. Aquaporins (AQPs), a family of membrane-bound water channels, are key players in osmoregulation in many animals but remain poorly characterized in ticks. Here, we bioinformatically identified AQP-like genes from the deer tick Ixodes scapularis and used phylogenetic approaches to map the evolution of the aquaporin gene family in arthropods. Most arachnid AQP-like sequences (including those of I. scapularis) formed a monophyletic group clustered within aquaglycerolporins (GLPs) from bacteria to vertebrates. This gene family is absent from insects, revealing divergent evolutionary paths for AQPs in different hematophagous arthropods. Next, we sequenced the full-length cDNA of I. scapularis aquaporin 1 (IsAQP1) and expressed it heterologously in Xenopus oocytes to functionally characterize its permeability to water and solutes. Additionally, we examined IsAQP1 expression across different life stages and adult female organs. We found IsAQP1 is an efficient water channel with high expression in salivary glands prior to feeding, suggesting it plays a role in osmoregulation before or during blood feeding. Its functional properties are unique: unlike most GLPs, IsAQP1 has low glycerol permeability, and unlike most AQPs, it is insensitive to mercury. Together, our results suggest IsAQP1 plays an important role in tick water balance physiology and that it may hold promise as a target of novel vector control efforts.

A Cambrian-Ordovician Terrestrialization of Arachnids

Understanding the temporal context of terrestrialization in chelicerates depends on whether terrestrial groups, the traditional Arachnida, have a single origin and whether or not horseshoe crabs are primitively or secondarily marine. Molecular dating on a phylogenomic tree that recovers arachnid monophyly, constrained by 27 rigorously vetted fossil calibrations, estimates that Arachnida originated during the Cambrian or Ordovician. After the common ancestor colonized the land, the main lineages appear to have rapidly radiated in the Cambrian-Ordovician boundary interval, coinciding with high rates of molecular evolution. The highest rates of arachnid diversification are detected between the Permian and Early Cretaceous. A pattern of ancient divergence estimates for terrestrial arthropod groups in the Cambrian while the oldest fossils are Silurian (seen in both myriapods and arachnids) is mirrored in the molecular and fossil records of land plants. We suggest the discrepancy between molecular and fossil evidence for terrestrialization is likely driven by the extreme sparseness of terrestrial sediments in the rock record before the late Silurian.

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

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

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

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

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.

Comparative transcriptome reveal the potential adaptive evolutionary genes in Andrias davidianus

To search the evidence of molecular evolution mechanism for aquatic and cave habitat in Andrias davidianus, the evolution analysis was carried out among several species transcriptome data. The transcriptome data of Notophthalmus viridescens, Xenopus tropicalis, Cynops pyrrhogaster, Hynobius chinensis and A. davidianus were obtained from the Genbank and reassembled except Xenopus tropicalis. The BLAST search of transcriptome data obtained 1244 single-copy orthologous genes among five species. A phylogenetic tree showed A. davidianus to have the closest relationship to H. chinensis. Fourteen positively selected genes were detected in A. davidianus and N. vridescens group and fifteen in A. davidianus and H. chinensis group. Five genes were shared in the both groups which involved in the immune system, suggesting that A. davidianus adaptation to an aquatic and cave environment required rapid evolution of the immune system compared to N. viridescens and H. chinensis.

Evolutionary origin of type IV classical cadherins in arthropods

Background

Classical cadherins are a metazoan-specific family of homophilic cell-cell adhesion molecules that regulate morphogenesis. Type I and type IV cadherins in this family function at adherens junctions in the major epithelial tissues of vertebrates and insects, respectively, but they have distinct, relatively simple domain organizations that are thought to have evolved by independent reductive changes from an ancestral type III cadherin, which is larger than derived paralogs and has a complicated domain organization. Although both type III and type IV cadherins have been identified in hexapods and branchiopods, the process by which the type IV cadherin evolved is still largely unclear.

Results

Through an analysis of arthropod genome sequences, we found that the only classical cadherin encoded in chelicerate genomes was the type III cadherin and that the two type III cadherin genes found in the spider Parasteatoda tepidariorum genome exhibited a complex yet ancestral exon-intron organization in arthropods. Genomic and transcriptomic data from branchiopod, copepod, isopod, amphipod, and decapod crustaceans led us to redefine the type IV cadherin category, which we separated into type IVa and type IVb, which displayed a similar domain organization, except type IVb cadherins have a larger number of extracellular cadherin (EC) domains than do type IVa cadherins (nine versus seven). We also showed that type IVa cadherin genes occurred in the hexapod, branchiopod, and copepod genomes whereas only type IVb cadherin genes were present in malacostracans. Furthermore, comparative characterization of the type IVb cadherins suggested that the presence of two extra EC domains in their N-terminal regions represented primitive characteristics. In addition, we identified an evolutionary loss of two highly conserved cysteine residues among the type IVa cadherins of insects.

Conclusions

We provide a genomic perspective of the evolution of classical cadherins among bilaterians, with a focus on the Arthropoda, and suggest that following the divergence of early arthropods, the precursor of the insect type IV cadherin evolved through stepwise reductive changes from the ancestral type III state. In addition, the complementary distributions of polarized genomic characters related to type IVa/IVb cadherins may have implications for our interpretations of pancrustacean phylogeny.

Electronic supplementary material

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

Segmentation in Tardigrada and diversification of segmental patterns in Panarthropoda

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

Construction and Characterization of Two Novel Transcriptome Assemblies in the Congeneric Porcelain Crabs Petrolisthes cinctipes and P. manimaculis

Crustaceans have commonly been used as non-model systems in basic biological research, especially physiological regulation. With the recent and rapid adoption of functional genomic tools, crustaceans are increasingly becoming model systems for ecological investigations of development and evolution and for mechanistic examinations of genotype-phenotype interactions and molecular pathways of response to environmental stressors. Comparative transcriptomic approaches, however, remain constrained by a lack of sequence data in closely related crustacean taxa. We identify challenges in the use of functional genomics tools in comparative analysis among decapod crustacean in light of recent advances. We present RNA-seq data from two congeneric species of porcelain crabs (Petrolisthes cinctipes and P. manimaculis) used to construct two de novo transcriptome assemblies with ∼194K and ∼278K contigs, respectively. We characterize and contrast these assemblies and compare them to a previously generated EST sequence library for P. cinctipes We also discuss the potential use of these data as a case-study system in the broader context of crustacean comparative transcriptomics.

Transcriptomics and venomics: implications for medicinal chemistry

Over the last three decades, transcriptomic studies of venom gland cells have continuously evolved, opening up new possibilities for exploring the molecular diversity of animal venoms, a prerequisite for the discovery of new drug candidates and molecular phylogenetics. The molecular complexity of animal venoms is much greater than initially thought. In this review, we describe the different technologies available for transcriptomic studies of venom, from the original individual cloning approaches to the more recent global Next Generation Sequencing strategies. Our understanding of animal venoms is evolving, with the discovery of complex and diverse bio-optimized cocktails of compounds, including mostly peptides and proteins, which are now beginning to be studied by academic and industrial researchers.

Whole Transcriptome of the Venom Gland from Urodacus yaschenkoi Scorpion

Australian scorpion venoms have been poorly studied, probably because they do not pose an evident threat to humans. In addition, the continent has other medically important venomous animals capable of causing serious health problems. Urodacus yaschenkoi belongs to the most widely distributed family of Australian scorpions (Urodacidae) and it is found all over the continent, making it a useful model system for studying venom composition and evolution. This communication reports the whole set of mRNA transcripts produced by the venom gland. U. yaschenkoi venom is as complex as its overseas counterparts. These transcripts certainly code for several components similar to known scorpion venom components, such as: alpha-KTxs, beta-KTxs, calcins, protease inhibitors, antimicrobial peptides, sodium-channel toxins, toxin-like peptides, allergens, La1-like, hyaluronidases, ribosomal proteins, proteasome components and proteins related to cellular processes. A comparison with the venom gland transcriptome of Centruroides noxius (Buthidae) showed that these two scorpions have similar components related to biological processes, although important differences occur among the venom toxins. In contrast, a comparison with sequences reported for Urodacus manicatus revealed that these two Urodacidae species possess the same subfamily of scorpion toxins. A comparison with sequences of an U. yaschenkoi cDNA library previously reported by our group showed that both techniques are reliable for the description of the venom components, but the whole transcriptome generated with Next Generation Sequencing platform provides sequences of all transcripts expressed. Several of which were identified in the proteome, but many more transcripts were identified including uncommon transcripts. The information reported here constitutes a reference for non-Buthidae scorpion venoms, providing a comprehensive view of genes that are involved in venom production. Further, this work identifies new putative bioactive compounds that could be used to seed research into new pharmacological compounds and increase our understanding of the function of different ion channels.

Transcriptome analysis of scorpion species belonging to the Vaejovis genus

Scorpions belonging to the Buthidae family have traditionally drawn much of the biochemist's attention due to the strong toxicity of their venoms. Scorpions not toxic to mammals, however, also have complex venoms. They have been shown to be an important source of bioactive peptides, some of them identified as potential drug candidates for the treatment of several emerging diseases and conditions. It is therefore important to characterize the large diversity of components found in the non-Buthidae venoms. As a contribution to this goal, this manuscript reports the construction and characterization of cDNA libraries from four scorpion species belonging to the Vaejovis genus of the Vaejovidae family: Vaejovis mexicanus, V. intrepidus, V. subcristatus and V. punctatus. Some sequences coding for channel-acting toxins were found, as expected, but the main transcribed genes in the glands actively producing venom were those coding for non disulfide-bridged peptides. The ESTs coding for putative channel-acting toxins, corresponded to sodium channel β toxins, to members of the potassium channel-acting α or κ families, and to calcium channel-acting toxins of the calcin family. Transcripts for scorpine-like peptides of two different lengths were found, with some of the species coding for the two kinds. One sequence coding for La1-like peptides, of yet unknown function, was found for each species. Finally, the most abundant transcripts corresponded to peptides belonging to the long chain multifunctional NDBP-2 family and to the short antimicrobials of the NDBP-4 family. This apparent venom composition is in correspondence with the data obtained to date for other non-Buthidae species. Our study constitutes the first approach to the characterization of the venom gland transcriptome for scorpion species belonging to the Vaejovidae family.

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.

Prediction of the first neuropeptides from a member of the Remipedia (Arthropoda, Crustacea)

The Remipedia is a small, recently described crustacean class that inhabits submerged marine/anchialine cave systems. Phylogenetic and morphological investigations support a sister group relationship between these animals and the hexapods. The recent deposition of numerous (>100,000) transcriptome shotgun assembly sequences for Speleonectes cf. tulumensis provides a unique resource to identify proteins of interest from a member of the Remipedia. Here, this dataset was mined for sequences encoding putative neuropeptide pre/preprohormones, with the mature peptides predicted from the deduced precursors using an established workflow. The structures of 40 mature peptides were obtained via this strategy, including members of 11 well-known arthropod peptide families (adipokinetic hormone/corazonin-like peptide [ACP], allatostatin A, allatostatin C, diuretic hormone 31, eclosion hormone, ion transport peptide/crustacean hyperglycemic hormone, neuropeptide F, proctolin, SIFamide, sulfakinin and tachykinin-related peptide); these are the only peptides thus far described from any member of the Remipedia. Comparison of the Speleonectes isoforms with those from other crustaceans and hexapods revealed the peptidome of this species to have characteristics of both subphyla (e.g. it possesses the stereotypical decapod crustacean SIFamide and tachykinin-related peptide isoforms, while simultaneously being the only crustacean with an insect AKC). Moreover, BLAST searches in which the deduced Speleonectes precursors were compared to the pancrustacean protein database most frequently returned insect homologs as the closest matches. The peptidomic analyses presented here are consistent with the hypothesized phylogenetic position of the Remipedia within the Pancrustacea, and serve as a foundation from which to launch future investigations of peptidergic signaling in remipedes.

Scorpion venom components that affect ion-channels function

The number and types of venom components that affect ion-channel function are reviewed. These are the most important venom components responsible for human intoxication, deserving medical attention, often requiring the use of specific anti-venoms. Special emphasis is given to peptides that recognize Na(+)-, K(+)- and Ca(++)-channels of excitable cells. Knowledge generated by direct isolation of peptides from venom and components deduced from cloned genes, whose amino acid sequences are deposited into databanks are nowadays in the order of 1.5 thousands, out of an estimate biodiversity closed to 300,000. Here the diversity of components is briefly reviewed with mention to specific references. Structural characteristic are discussed with examples taken from published work. The principal mechanisms of action of the three different types of peptides are also reviewed. Na(+)-channel specific venom components usually are modifier of the open and closing kinetic mechanisms of the ion-channels, whereas peptides affecting K(+)-channels are normally pore blocking agents. The Ryanodine Ca(++)-channel specific peptides are known for causing sub-conducting stages of the channels conductance and some were shown to be able to internalize penetrating inside the muscle cells.

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.

Ordinal-level phylogenomics of the arthropod class Diplopoda (millipedes) based on an analysis of 221 nuclear protein-coding loci generated using next-generation sequence analyses

BACKGROUND

The ancient and diverse, yet understudied arthropod class Diplopoda, the millipedes, has a muddled taxonomic history. Despite having a cosmopolitan distribution and a number of unique and interesting characteristics, the group has received relatively little attention; interest in millipede systematics is low compared to taxa of comparable diversity. The existing classification of the group comprises 16 orders. Past attempts to reconstruct millipede phylogenies have suffered from a paucity of characters and included too few taxa to confidently resolve relationships and make formal nomenclatural changes. Herein, we reconstruct an ordinal-level phylogeny for the class Diplopoda using the largest character set ever assembled for the group.

METHODS

Transcriptomic sequences were obtained from exemplar taxa representing much of the diversity of millipede orders using second-generation (i.e., next-generation or high-throughput) sequencing. These data were subject to rigorous orthology selection and phylogenetic dataset optimization and then used to reconstruct phylogenies employing Bayesian inference and maximum likelihood optimality criteria. Ancestral reconstructions of sperm transfer appendage development (gonopods), presence of lateral defense secretion pores (ozopores), and presence of spinnerets were considered. The timings of major millipede lineage divergence points were estimated.

RESULTS

The resulting phylogeny differed from the existing classifications in a number of fundamental ways. Our phylogeny includes a grouping that has never been described (Juliformia+Merocheta+Stemmiulida), and the ancestral reconstructions suggest caution with respect to using spinnerets as a unifying characteristic for the Nematophora. Our results are shown to have significantly stronger support than previous hypotheses given our data. Our efforts represent the first step toward obtaining a well-supported and robust phylogeny of the Diplopoda that can be used to answer many questions concerning the evolution of this ancient and diverse animal group.

Selective neuronal staining in tardigrades and onychophorans provides insights into the evolution of segmental ganglia in panarthropods

BACKGROUND

Although molecular analyses have contributed to a better resolution of the animal tree of life, the phylogenetic position of tardigrades (water bears) is still controversial, as they have been united alternatively with nematodes, arthropods, onychophorans (velvet worms), or onychophorans plus arthropods. Depending on the hypothesis favoured, segmental ganglia in tardigrades and arthropods might either have evolved independently, or they might well be homologous, suggesting that they were either lost in onychophorans or are a synapomorphy of tardigrades and arthropods. To evaluate these alternatives, we analysed the organisation of the nervous system in three tardigrade species using antisera directed against tyrosinated and acetylated tubulin, the amine transmitter serotonin, and the invertebrate neuropeptides FMRFamide, allatostatin and perisulfakinin. In addition, we performed retrograde staining of nerves in the onychophoran Euperipatoides rowelli in order to compare the serial locations of motor neurons within the nervous system relative to the appendages they serve in arthropods, tardigrades and onychophorans.

RESULTS

Contrary to a previous report from a Macrobiotus species, our immunocytochemical and electron microscopic data revealed contralateral fibres and bundles of neurites in each trunk ganglion of three tardigrade species, including Macrobiotus cf. harmsworthi, Paramacrobiotus richtersi and Hypsibius dujardini. Moreover, we identified additional, extra-ganglionic commissures in the interpedal regions bridging the paired longitudinal connectives. Within the ganglia we found serially repeated sets of serotonin- and RFamid-like immunoreactive neurons. Furthermore, our data show that the trunk ganglia of tardigrades, which include the somata of motor neurons, are shifted anteriorly with respect to each corresponding leg pair, whereas no such shift is evident in the arrangement of motor neurons in the onychophoran nerve cords.

CONCLUSIONS

Taken together, these data reveal three major correspondences between the segmental ganglia of tardigrades and arthropods, including (i) contralateral projections and commissures in each ganglion, (ii) segmentally repeated sets of immunoreactive neurons, and (iii) an anteriorly shifted (parasegmental) position of ganglia. These correspondences support the homology of segmental ganglia in tardigrades and arthropods, suggesting that these structures were either lost in Onychophora or, alternatively, evolved in the tardigrade/arthropod lineage.

Mass fingerprinting of the venom and transcriptome of venom gland of scorpion Centruroides tecomanus

Centruroides tecomanus is a Mexican scorpion endemic of the State of Colima, that causes human fatalities. This communication describes a proteome analysis obtained from milked venom and a transcriptome analysis from a cDNA library constructed from two pairs of venom glands of this scorpion. High perfomance liquid chromatography separation of soluble venom produced 80 fractions, from which at least 104 individual components were identified by mass spectrometry analysis, showing to contain molecular masses from 259 to 44,392 Da. Most of these components are within the expected molecular masses for Na⁺- and K⁺-channel specific toxic peptides, supporting the clinical findings of intoxication, when humans are stung by this scorpion. From the cDNA library 162 clones were randomly chosen, from which 130 sequences of good quality were identified and were clustered in 28 contigs containing, each, two or more expressed sequence tags (EST) and 49 singlets with only one EST. Deduced amino acid sequence analysis from 53% of the total ESTs showed that 81% (24 sequences) are similar to known toxic peptides that affect Na⁺-channel activity, and 19% (7 unique sequences) are similar to K⁺-channel especific toxins. Out of the 31 sequences, at least 8 peptides were confirmed by direct Edman degradation, using components isolated directly from the venom. The remaining 19%, 4%, 4%, 15% and 5% of the ESTs correspond respectively to proteins involved in cellular processes, antimicrobial peptides, venom components, proteins without defined function and sequences without similarity in databases. Among the cloned genes are those similar to metalloproteinases.

Transcriptome characterisation of Pinus tabuliformis and evolution of genes in the Pinus phylogeny

BACKGROUND

The Chinese pine (Pinus tabuliformis) is an indigenous conifer species in northern China but is relatively underdeveloped as a genomic resource; thus, limiting gene discovery and breeding. Large-scale transcriptome data were obtained using a next-generation sequencing platform to compensate for the lack of P. tabuliformis genomic information.

RESULTS

The increasing amount of transcriptome data on Pinus provides an excellent resource for multi-gene phylogenetic analysis and studies on how conserved genes and functions are maintained in the face of species divergence. The first P. tabuliformis transcriptome from a normalised cDNA library of multiple tissues and individuals was sequenced in a full 454 GS-FLX run, producing 911,302 sequencing reads. The high quality overlapping expressed sequence tags (ESTs) were assembled into 46,584 putative transcripts, and more than 700 SSRs and 92,000 SNPs/InDels were characterised. Comparative analysis of the transcriptome of six conifer species yielded 191 orthologues, from which we inferred a phylogenetic tree, evolutionary patterns and calculated rates of gene diversion. We also identified 938 fast evolving sequences that may be useful for identifying genes that perhaps evolved in response to positive selection and might be responsible for speciation in the Pinus lineage.

CONCLUSIONS

A large collection of high-quality ESTs was obtained, de novo assembled and characterised, which represents a dramatic expansion of the current transcript catalogues of P. tabuliformis and which will gradually be applied in breeding programs of P. tabuliformis. Furthermore, these data will facilitate future studies of the comparative genomics of P. tabuliformis and other related species.

Neural markers reveal a one-segmented head in tardigrades (water bears)

BACKGROUND

While recent neuroanatomical and gene expression studies have clarified the alignment of cephalic segments in arthropods and onychophorans, the identity of head segments in tardigrades remains controversial. In particular, it is unclear whether the tardigrade head and its enclosed brain comprises one, or several segments, or a non-segmental structure. To clarify this, we applied a variety of histochemical and immunocytochemical markers to specimens of the tardigrade Macrobiotus cf. harmsworthi and the onychophoran Euperipatoides rowelli.

METHODOLOGY/PRINCIPAL FINDINGS

Our immunolabelling against serotonin, FMRFamide and α-tubulin reveals that the tardigrade brain is a dorsal, bilaterally symmetric structure that resembles the brain of onychophorans and arthropods rather than a circumoesophageal ring typical of cycloneuralians (nematodes and allies). A suboesophageal ganglion is clearly lacking. Our data further reveal a hitherto unknown, unpaired stomatogastric ganglion in Macrobiotus cf. harmsworthi, which innervates the ectodermal oesophagus and the endodermal midgut and is associated with the second leg-bearing segment. In contrast, the oesophagus of the onychophoran E. rowelli possesses no immunoreactive neurons, whereas scattered bipolar, serotonin-like immunoreactive cell bodies are found in the midgut wall. Furthermore, our results show that the onychophoran pharynx is innervated by a medullary loop nerve accompanied by monopolar, serotonin-like immunoreactive cell bodies.

CONCLUSIONS/SIGNIFICANCE

A comparison of the nervous system innervating the foregut and midgut structures in tardigrades and onychophorans to that of arthropods indicates that the stomatogastric ganglion is a potential synapomorphy of Tardigrada and Arthropoda. Its association with the second leg-bearing segment in tardigrades suggests that the second trunk ganglion is a homologue of the arthropod tritocerebrum, whereas the first ganglion corresponds to the deutocerebrum. We therefore conclude that the tardigrade brain consists of a single segmental region corresponding to the arthropod protocerebrum and, accordingly, that the tardigrade head is a non-composite, one-segmented structure.

A catalogue of putative unique transcripts from Douglas-fir (Pseudotsuga menziesii) based on 454 transcriptome sequencing of genetically diverse, drought stressed seedlings

Background

Douglas-fir (Pseudotsuga menziesii) extends over a wide range of contrasting environmental conditions, reflecting substantial local adaptation. For this reason, it is an interesting model species to study plant adaptation and the effects of global climate change such as increased temperatures and significant periods of drought on individual trees and the forest landscape in general. However, genomic data and tools for studying genetic variation in natural populations to understand the genetic and physiological mechanisms of adaptation are currently missing for Douglas-fir. This study represents a first step towards characterizing the Douglas-fir transcriptome based on 454 sequencing of twelve cDNA libraries. The libraries were constructed from needle and wood tissue of coastal and interior provenances subjected to drought stress experiments.

Results

The 454 sequencing of twelve normalized cDNA libraries resulted in 3.6 million reads from which a set of 170,859 putative unique transcripts (PUTs) was assembled. Functional annotation by BLAST searches and Gene Ontology mapping showed that the composition of functional classes is very similar to other plant transcriptomes and demonstrated that a large fraction of the Douglas-fir transcriptome is tagged by the PUTs. Based on evolutionary conservation, we identified about 1,000 candidate genes related to drought stress. A total number of 187,653 single nucleotide polymorphisms (SNPs) were detected by three SNP detection tools. However, only 27,688 SNPs were identified by all three methods, indicating that SNP detection depends on the particular method used. The two alleles of about 60% of the 27,688 SNPs are segregating simultaneously in both coastal and interior provenances, which indicates a high proportion of ancestral shared polymorphisms or a high level of gene flow between these two ecologically and phenotypically different varieties.

Conclusions

We established a catalogue of PUTs and large SNP database for Douglas-fir. Both will serve as a useful resource for the further characterization of the genome and transcriptome of Douglas-fir and for the analysis of genetic variation using genotyping or resequencing methods.

African swine fever virus controls the host transcription and cellular machinery of protein synthesis

Throughout a viral infection, the infected cell reprograms the gene expression pattern in order to establish a satisfactory antiviral response. African swine fever virus (ASFV), like other complex DNA viruses, sets up a number of strategies to evade the host's defense systems, such as apoptosis, inflammation and immune responses. The capability of the virus to persist in its natural hosts and in domestic pigs, which recover from infection with less virulent isolates, suggests that the virus displays effective mechanisms to escape host defense systems. ASFV has been described to regulate the activation of several transcription factors, thus regulating the activation of specific target genes during ASFV infection. Whereas some reports have concerned about anti-apoptotic ASFV genes and the molecular mechanisms by which ASFV interferes with inducible gene transcription and immune evasion, less is yet known regarding how ASFV regulates the translational machinery in infected cells, although a recent report has shown a mechanism for favored expression of viral genes based on compartmentalization of viral mRNA and ribosomes with cellular translation factors within the virus factory. The viral mechanisms involved both in the regulation of host genes transcription and in the control of cellular protein synthesis are summarized in this review.

Global transcriptome analysis of the scorpion Centruroides noxius: new toxin families and evolutionary insights from an ancestral scorpion species

Scorpion venoms have been studied for decades, leading to the identification of hundreds of different toxins with medical and pharmacological implications. However, little emphasis has been given to the description of these arthropods from cellular and evolutionary perspectives. In this report, we describe a transcriptomic analysis of the Mexican scorpion Centruroides noxius Hoffmann, performed with a pyrosequencing platform. Three independent sequencing experiments were carried out, each including three different cDNA libraries constructed from RNA extracted from the whole body of the scorpion after telson removal, and from the venom gland before and after venom extraction. Over three million reads were obtained and assembled in almost 19000 isogroups. Within the telson-specific sequences, 72 isogroups (0.4% of total unique transcripts) were found to be similar to toxins previously reported in other scorpion species, spiders and sea anemones. The annotation pipeline also revealed the presence of important elements of the small non-coding RNA processing machinery, as well as microRNA candidates. A phylogenomic analysis of concatenated essential genes evidenced differential evolution rates in this species, particularly in ribosomal proteins and proteasome components. Additionally, statistical comparison of transcript abundance before and after venom extraction showed that 3% and 2% of the assembled isogroups had higher expression levels in the active and replenishing gland, respectively. Thus, our sequencing and annotation strategies provide a general view of the cellular and molecular processes that take place in these arthropods, allowed the discovery of new pharmacological and biotechnological targets and uncovered several regulatory and metabolic responses behind the assembly of the scorpion venom. The results obtained in this report represent the first high-throughput study that thoroughly describes the universe of genes that are expressed in the scorpion Centruroides noxius Hoffmann, a highly relevant organism from medical and evolutionary perspectives.

Transcriptome analysis and SSR/SNP markers information of the blunt snout bream (Megalobrama amblycephala)

Background

Blunt snout bream (Megalobrama amblycephala) is an herbivorous freshwater fish species native to China and has been recognized as a main aquaculture species in the Chinese freshwater polyculture system with high economic value. Right now, only limited EST resources were available for M. amblycephala. Recent advances in large-scale RNA sequencing provide a fast, cost-effective, and reliable approach to generate large expression datasets for functional genomic analysis, which is especially suitable for non-model species with un-sequenced genomes.

Methodology and Principal Findings

Using 454 pyrosequencing, a total of 1,409,706 high quality reads (total length 577 Mbp) were generated from the normalized cDNA of pooled M. amblycephala individuals. These sequences were assembled into 26,802 contigs and 73,675 singletons. After BLAST searches against the NCBI non-redundant (NR) and UniProt databases with an arbitrary expectation value of E⁻¹⁰, over 40,000 unigenes were functionally annotated and classified using the FunCat functional annotation scheme. A comparative genomics approach revealed a substantial proportion of genes expressed in M. amblycephala tanscriptome to be shared across the genomes of zebrafish, medaka, tetraodon, fugu, stickleback, human, mouse, and chicken, and identified a substantial number of potentially novel M. amblycephala genes. A total number of 4,952 SSRs were found and 116 polymorphic loci have been characterized. A significant number of SNPs (25,697) and indels (23,287) were identified based on specific filter criteria in the M. amblycephala.

Conclusions

This study is the first comprehensive transcriptome analysis for a fish species belonging to the genus Megalobrama. These large EST resources are expected to be valuable for the development of molecular markers, construction of gene-based linkage map, and large-scale expression analysis of M. amblycephala, as well as comparative genome analysis for the genus Megalobrama fish species. The identified SSR and SNP markers will greatly benefit its breeding program and whole genome association studies.

The diversity and evolution of chelicerate hemocyanins

BACKGROUND

Oxygen transport in the hemolymph of many arthropod species is facilitated by large copper-proteins referred to as hemocyanins. Arthropod hemocyanins are hexamers or oligomers of hexamers, which are characterized by a high O2 transport capacity and a high cooperativity, thereby enhancing O2 supply. Hemocyanin subunit sequences had been available from horseshoe crabs (Xiphosura) and various spiders (Araneae), but not from any other chelicerate taxon. To trace the evolution of hemocyanins and the emergence of the large hemocyanin oligomers, hemocyanin cDNA sequences were obtained from representatives of selected chelicerate classes.

RESULTS

Hemocyanin subunits from a sea spider, a scorpion, a whip scorpion and a whip spider were sequenced. Hemocyanin has been lost in Opiliones, Pseudoscorpiones, Solifugae and Acari, which may be explained by the evolution of trachea (i.e., taxon Apulmonata). Bayesian phylogenetic analysis was used to reconstruct the evolution of hemocyanin subunits and a relaxed molecular clock approach was applied to date the major events. While the sea spider has a simple hexameric hemocyanin, four distinct subunit types evolved before Xiphosura and Arachnida diverged around 470 Ma ago, suggesting the existence of a 4 × 6mer at that time. Subsequently, independent gene duplication events gave rise to the other distinct subunits in each of the 8 × 6mer hemocyanin of Xiphosura and the 4 × 6mer of Arachnida. The hemocyanin sequences were used to infer the evolutionary history of chelicerates. The phylogenetic trees support a basal position of Pycnogonida, a sister group relationship of Xiphosura and Arachnida, and a sister group relationship of the whip scorpions and the whip spiders.

CONCLUSION

Formation of a complex hemocyanin oligomer commenced early in the evolution of euchelicerates. A 4 × 6mer hemocyanin consisting of seven subunit types is conserved in most arachnids since more than 400 Ma, although some entelegyne spiders display selective subunit loss and independent oligomerization. Hemocyanins also turned out to be a good marker to trace chelicerate evolution, which is, however, limited by the loss of hemocyanin in some taxa. The molecular clock calculations were in excellent agreement with the fossil record, also demonstrating the applicability of hemocyanins for such approach.

De novo assembly and characterization of a maternal and developmental transcriptome for the emerging model crustacean Parhyale hawaiensis

BACKGROUND

Arthropods are the most diverse animal phylum, but their genomic resources are relatively few. While the genome of the branchiopod Daphnia pulex is now available, no other large-scale crustacean genomic resources are available for comparison. In particular, genomic resources are lacking for the most tractable laboratory model of crustacean development, the amphipod Parhyale hawaiensis. Insight into shared and divergent characters of crustacean genomes will facilitate interpretation of future developmental, biomedical, and ecological research using crustacean models.

RESULTS

To generate a transcriptome enriched for maternally provided and zygotically transcribed developmental genes, we created cDNA from ovaries and embryos of P. hawaiensis. Using 454 pyrosequencing, we sequenced over 1.1 billion bases of this cDNA, and assembled them de novo to create, to our knowledge, the second largest crustacean genomic resource to date. We found an unusually high proportion of C2H2 zinc finger-containing transcripts, as has also been reported for the genome of the pea aphid Acyrthosiphon pisum. Consistent with previous reports, we detected trans-spliced transcripts, but found that they did not noticeably impact transcriptome assembly. Our assembly products yielded 19,067 unique BLAST hits against nr (E-value cutoff e-10). These included over 400 predicted transcripts with significant similarity to D. pulex sequences but not to sequences of any other animal. Annotation of several hundred genes revealed P. hawaiensis homologues of genes involved in development, gametogenesis, and a majority of the members of six major conserved metazoan signaling pathways.

CONCLUSIONS

The amphipod P. hawaiensis has higher transcript complexity than known insect transcriptomes, and trans-splicing does not appear to be a major contributor to this complexity. We discuss the importance of a reliable comparative genomic framework within which to consider findings from new crustacean models such as D. pulex and P. hawaiensis, as well as the need for development of further substantial crustacean genomic resources.

Scorpion and spider venom peptides: gene cloning and peptide expression

This communication reviews most of the important findings related to venom components isolated from scorpions and spiders, mainly by means of gene cloning and expression. Rather than revising results obtained by classical biochemical studies that report structure and function of venom components, here the emphasis is placed on cloning and identification of genes present in the venomous glands of these arachnids. Aspects related to cDNA library construction, specific or random ESTs cloning, transcriptome analysis, high-throughput screening, heterologous expression and folding are briefly discussed, showing some numbers of species and components already identified, but also shortly mentioning limitations and perspectives of research for the future in this field.

MicroRNAs and phylogenomics resolve the relationships of Tardigrada and suggest that velvet worms are the sister group of Arthropoda

Morphological data traditionally group Tardigrada (water bears), Onychophora (velvet worms), and Arthropoda (e.g., spiders, insects, and their allies) into a monophyletic group of invertebrates with walking appendages known as the Panarthropoda. However, molecular data generally do not support the inclusion of tardigrades within the Panarthropoda, but instead place them closer to Nematoda (roundworms). Here we present results from the analyses of two independent genomic datasets, expressed sequence tags (ESTs) and microRNAs (miRNAs), which congruently resolve the phylogenetic relationships of Tardigrada. Our EST analyses, based on 49,023 amino acid sites from 255 proteins, significantly support a monophyletic Panarthropoda including Tardigrada and suggest a sister group relationship between Arthropoda and Onychophora. Using careful experimental manipulations--comparisons of model fit, signal dissection, and taxonomic pruning--we show that support for a Tardigrada + Nematoda group derives from the phylogenetic artifact of long-branch attraction. Our small RNA libraries fully support our EST results; no miRNAs were found to link Tardigrada and Nematoda, whereas all panarthropods were found to share one unique miRNA (miR-276). In addition, Onychophora and Arthropoda were found to share a second miRNA (miR-305). Our study confirms the monophyly of the legged ecdysozoans, shows that past support for a Tardigrada + Nematoda group was due to long-branch attraction, and suggests that the velvet worms are the sister group to the arthropods.

Reevaluating the arthropod tree of life

Arthropods are the most diverse group of animals and have been so since the Cambrian radiation. They belong to the protostome clade Ecdysozoa, with Onychophora (velvet worms) as their most likely sister group and tardigrades (water bears) the next closest relative. The arthropod tree of life can be interpreted as a five-taxon network, containing Pycnogonida, Euchelicerata, Myriapoda, Crustacea, and Hexapoda, the last two forming the clade Tetraconata or Pancrustacea. The unrooted relationship of Tetraconata to the three other lineages is well established, but of three possible rooting positions the Mandibulata hypothesis receives the most support. Novel approaches to studying anatomy with noninvasive three-dimensional reconstruction techniques, the application of these techniques to new and old fossils, and the so-called next-generation sequencing techniques are at the forefront of understanding arthropod relationships. Cambrian fossils assigned to the arthropod stem group inform on the origin of arthropod characters from a lobopodian ancestry. Monophyly of Pycnogonida, Euchelicerata, Myriapoda, Tetraconata, and Hexapoda is well supported, but the interrelationships of arachnid orders and the details of crustacean paraphyly with respect to Hexapoda remain the major unsolved phylogenetic problems.

Sources of signal in 62 protein-coding nuclear genes for higher-level phylogenetics of arthropods

BACKGROUND

This study aims to investigate the strength of various sources of phylogenetic information that led to recent seemingly robust conclusions about higher-level arthropod phylogeny and to assess the role of excluding or downweighting synonymous change for arriving at those conclusions.

METHODOLOGY/PRINCIPAL FINDINGS

The current study analyzes DNA sequences from 68 gene segments of 62 distinct protein-coding nuclear genes for 80 species. Gene segments analyzed individually support numerous nodes recovered in combined-gene analyses, but few of the higher-level nodes of greatest current interest. However, neither is there support for conflicting alternatives to these higher-level nodes. Gene segments with higher rates of nonsynonymous change tend to be more informative overall, but those with lower rates tend to provide stronger support for deeper nodes. Higher-level nodes with bootstrap values in the 80% - 99% range for the complete data matrix are markedly more sensitive to substantial drops in their bootstrap percentages after character subsampling than those with 100% bootstrap, suggesting that these nodes are likely not to have been strongly supported with many fewer data than in the full matrix. Data set partitioning of total data by (mostly) synonymous and (mostly) nonsynonymous change improves overall node support, but the result remains much inferior to analysis of (unpartitioned) nonsynonymous change alone. Clusters of genes with similar nonsynonymous rate properties (e.g., faster vs. slower) show some distinct patterns of node support but few conflicts. Synonymous change is shown to contribute little, if any, phylogenetic signal to the support of higher-level nodes, but it does contribute nonphylogenetic signal, probably through its underlying heterogeneous nucleotide composition. Analysis of seemingly conservative indels does not prove useful.

CONCLUSIONS

Generating a robust molecular higher-level phylogeny of Arthropoda is currently possible with large amounts of data and an exclusive reliance on nonsynonymous change.

Applications of next generation sequencing in molecular ecology of non-model organisms

As most biologists are probably aware, technological advances in molecular biology during the last few years have opened up possibilities to rapidly generate large-scale sequencing data from non-model organisms at a reasonable cost. In an era when virtually any study organism can 'go genomic', it is worthwhile to review how this may impact molecular ecology. The first studies to put the next generation sequencing (NGS) to the test in ecologically well-characterized species without previous genome information were published in 2007 and the beginning of 2008. Since then several studies have followed in their footsteps, and a large number are undoubtedly under way. This review focuses on how NGS has been, and can be, applied to ecological, population genetic and conservation genetic studies of non-model species, in which there is no (or very limited) genomic resources. Our aim is to draw attention to the various possibilities that are opening up using the new technologies, but we also highlight some of the pitfalls and drawbacks with these methods. We will try to provide a snapshot of the current state of the art for this rapidly advancing and expanding field of research and give some likely directions for future developments.

A new view of insect-crustacean relationships II. Inferences from expressed sequence tags and comparisons with neural cladistics

The enormous diversity of Arthropoda has complicated attempts by systematists to deduce the history of this group in terms of phylogenetic relationships and phenotypic change. Traditional hypotheses regarding the relationships of the major arthropod groups (Chelicerata, Myriapoda, Crustacea, and Hexapoda) focus on suites of morphological characters, whereas phylogenomics relies on large amounts of molecular sequence data to infer evolutionary relationships. The present discussion is based on expressed sequence tags (ESTs) that provide large numbers of short molecular sequences and so provide an abundant source of sequence data for phylogenetic inference. This study presents well-supported phylogenies of diverse arthropod and metazoan outgroup taxa obtained from publicly-available databases. An in-house bioinformatics pipeline has been used to compile and align conserved orthologs from each taxon for maximum likelihood inferences. This approach resolves many currently accepted hypotheses regarding internal relationships between the major groups of Arthropoda, including monophyletic Hexapoda, Tetraconata (Crustacea + Hexapoda), Myriapoda, and Chelicerata sensu lato (Pycnogonida + Euchelicerata). "Crustacea" is a paraphyletic group with some taxa more closely related to the monophyletic Hexapoda. These results support studies that have utilized more restricted EST data for phylogenetic inference, yet they differ in important regards from recently published phylogenies employing nuclear protein-coding sequences. The present results do not, however, depart from other phylogenies that resolve Branchiopoda as the crustacean sister group of Hexapoda. Like other molecular phylogenies, EST-derived phylogenies alone are unable to resolve morphological convergences or evolved reversals and thus omit what may be crucial events in the history of life. For example, molecular data are unable to resolve whether a Hexapod-Branchiopod sister relationship infers a branchiopod-like ancestry of the Hexapoda, or whether this assemblage originates from a malacostracan-like ancestor, with the morphologically simpler Branchiopoda being highly derived. Whereas this study supports many internal arthropod relationships obtained by other sources of molecular data, other approaches are required to resolve such evolutionary scenarios. The approach presented here turns out to be essential: integrating results of molecular phylogenetics and neural cladistics to infer that Branchiopoda evolved simplification from a more elaborate ancestor. Whereas the phenomenon of evolved simplification may be widespread, it is largely invisible to molecular techniques unless these are performed in conjunction with morphology-based strategies.

The origins of the arthropod nervous system: insights from the Onychophora

A revision of evolutionary relationships of the Arthropoda has provided fresh impetus to tracing the origins of the nervous system of this group of animals: other members of the Ecdysozoa possess a markedly different type of nervous system from both the arthropods and the annelid worms, with which they were previously grouped. Given their status as favoured sister taxon of the arthropods, Onychophora (velvet worms) are a key group for understanding the evolutionary changes that have taken place in the panarthropod (Arthropoda + Onychophora + Tardigrada) lineage. This article reviews our current knowledge of the structure and development of the onychophoran nervous system. The picture that emerges from these studies is that the nervous system of the panarthropod ancestor was substantially different from that of modern arthropods: this animal probably possessed a bipartite, rather than a tripartite brain; its nerve cord displayed only a limited degree of segmentation; and neurons were more numerous but more uniform in morphology than in living arthropods. These observations suggest an evolutionary scenario, by which the arthropod nervous system evolved from a system of orthogonally crossing nerve tracts present in both a presumed protostome ancestor and many extant worm-like invertebrates, including the onychophorans.

A new view of insect-crustacean relationships I. Inferences from neural cladistics and comparative neuroanatomy

Traditional hypotheses regarding the relationships of the major arthropod lineages focus on suites of comparable characters, often those that address features of the exoskeleton. However, because of the enormous morphological variety among arthropods, external characters may lead to ambiguities of interpretation and definition, particularly when species have undergone evolutionary simplification and reversal. Here we present the results of a cladistic analysis using morphological characters associated with brains and central nervous systems, based on the evidence that cerebral organization is generally robust over geological time. Well-resolved, strongly supported phylogenies were obtained from a neuromorphological character set representing a variety of discrete neuroanatomical traits. Phylogenetic hypotheses from this analysis support many accepted relationships, including monophyletic Chelicerata, Myriapoda, and Hexapoda, paraphyletic Crustacea and the union of Hexapoda and Crustacea (Tetraconata). They also support Mandibulata (Myriapoda + Tetraconata). One problematic result, which can be explained by symplesiomorphies that are likely to have evolved in deep time, is the inability to resolve Onychophora as a taxon distinct from Arthropoda. Crucially, neuronal cladistics supports the heterodox conclusion that both Hexapoda and Malacostraca are derived from a common ancestor that possessed a suite of discrete neural centers comprising an elaborate brain. Remipedes and copepods, both resolved as basal to Branchiopoda share a neural ground pattern with Malacostraca. These findings distinguish Hexapoda (Insecta) from Branchiopoda, which is the sister group of the clade Malacostraca + Hexapoda. The present study resolves branchiopod crustaceans as descendents of an ancestor with a complex brain, which means that they have evolved secondary simplification and the loss or reduction of numerous neural systems.

De novo characterization of the gametophyte transcriptome in bracken fern, Pteridium aquilinum

Background

Because of their phylogenetic position and unique characteristics of their biology and life cycle, ferns represent an important lineage for studying the evolution of land plants. Large and complex genomes in ferns combined with the absence of economically important species have been a barrier to the development of genomic resources. However, high throughput sequencing technologies are now being widely applied to non-model species. We leveraged the Roche 454 GS-FLX Titanium pyrosequencing platform in sequencing the gametophyte transcriptome of bracken fern (Pteridium aquilinum) to develop genomic resources for evolutionary studies.

Results

681,722 quality and adapter trimmed reads totaling 254 Mbp were assembled de novo into 56,256 unique sequences (i.e. unigenes) with a mean length of 547.2 bp and a total assembly size of 30.8 Mbp with an average read-depth coverage of 7.0×. We estimate that 87% of the complete transcriptome has been sequenced and that all transcripts have been tagged. 61.8% of the unigenes had blastx hits in the NCBI nr protein database, representing 22,596 unique best hits. The longest open reading frame in 52.2% of the unigenes had positive domain matches in InterProScan searches. We assigned 46.2% of the unigenes with a GO functional annotation and 16.0% with an enzyme code annotation. Enzyme codes were used to retrieve and color KEGG pathway maps. A comparative genomics approach revealed a substantial proportion of genes expressed in bracken gametophytes to be shared across the genomes of Arabidopsis, Selaginella and Physcomitrella, and identified a substantial number of potentially novel fern genes. By comparing the list of Arabidopsis genes identified by blast with a list of gametophyte-specific Arabidopsis genes taken from the literature, we identified a set of potentially conserved gametophyte specific genes. We screened unigenes for repetitive sequences to identify 548 potentially-amplifiable simple sequence repeat loci and 689 expressed transposable elements.

Conclusions

This study is the first comprehensive transcriptome analysis for a fern and represents an important scientific resource for comparative evolutionary and functional genomics studies in land plants. We demonstrate the utility of high-throughput sequencing of a normalized cDNA library for de novo transcriptome characterization and gene discovery in a non-model plant.

The maternal and early embryonic transcriptome of the milkweed bug Oncopeltus fasciatus

BACKGROUND

Most evolutionary developmental biology ("evo-devo") studies of emerging model organisms focus on small numbers of candidate genes cloned individually using degenerate PCR. However, newly available sequencing technologies such as 454 pyrosequencing have recently begun to allow for massive gene discovery in animals without sequenced genomes. Within insects, although large volumes of sequence data are available for holometabolous insects, developmental studies of basally branching hemimetabolous insects typically suffer from low rates of gene discovery.

RESULTS

We used 454 pyrosequencing to sequence over 500 million bases of cDNA from the ovaries and embryos of the milkweed bug Oncopeltus fasciatus, which lacks a sequenced genome. This indirectly developing insect occupies an important phylogenetic position, branching basal to Diptera (including fruit flies) and Hymenoptera (including honeybees), and is an experimentally tractable model for short-germ development. 2,087,410 reads from both normalized and non-normalized cDNA assembled into 21,097 sequences (isotigs) and 112,531 singletons. The assembled sequences fell into 16,617 unique gene models, and included predictions of splicing isoforms, which we examined experimentally. Discovery of new genes plateaued after assembly of ~1.5 million reads, suggesting that we have sequenced nearly all transcripts present in the cDNA sampled. Many transcripts have been assembled at close to full length, and there is a net gain of sequence data for over half of the pre-existing O. fasciatus accessions for developmental genes in GenBank. We identified 10,775 unique genes, including members of all major conserved metazoan signaling pathways and genes involved in several major categories of early developmental processes. We also specifically address the effects of cDNA normalization on gene discovery in de novo transcriptome analyses.

CONCLUSIONS

Our sequencing, assembly and annotation framework provide a simple and effective way to achieve high-throughput gene discovery for organisms lacking a sequenced genome. These data will have applications to the study of the evolution of arthropod genes and genetic pathways, and to the wider evolution, development and genomics communities working with emerging model organisms.[The sequence data from this study have been submitted to GenBank under study accession number SRP002610 (http://www.ncbi.nlm.nih.gov/sra?term=SRP002610). Custom scripts generated are available at http://www.extavourlab.com/protocols/index.html. Seven Additional files are available.].

De novo sequencing and characterization of floral transcriptome in two species of buckwheat (Fagopyrum)

BACKGROUND

Transcriptome sequencing data has become an integral component of modern genetics, genomics and evolutionary biology. However, despite advances in the technologies of DNA sequencing, such data are lacking for many groups of living organisms, in particular, many plant taxa. We present here the results of transcriptome sequencing for two closely related plant species. These species, Fagopyrum esculentum and F. tataricum, belong to the order Caryophyllales--a large group of flowering plants with uncertain evolutionary relationships. F. esculentum (common buckwheat) is also an important food crop. Despite these practical and evolutionary considerations Fagopyrum species have not been the subject of large-scale sequencing projects.

RESULTS

Normalized cDNA corresponding to genes expressed in flowers and inflorescences of F. esculentum and F. tataricum was sequenced using the 454 pyrosequencing technology. This resulted in 267 (for F. esculentum) and 229 (F. tataricum) thousands of reads with average length of 341-349 nucleotides. De novo assembly of the reads produced about 25 thousands of contigs for each species, with 7.5-8.2× coverage. Comparative analysis of two transcriptomes demonstrated their overall similarity but also revealed genes that are presumably differentially expressed. Among them are retrotransposon genes and genes involved in sugar biosynthesis and metabolism. Thirteen single-copy genes were used for phylogenetic analysis; the resulting trees are largely consistent with those inferred from multigenic plastid datasets. The sister relationships of the Caryophyllales and asterids now gained high support from nuclear gene sequences.

CONCLUSIONS

454 transcriptome sequencing and de novo assembly was performed for two congeneric flowering plant species, F. esculentum and F. tataricum. As a result, a large set of cDNA sequences that represent orthologs of known plant genes as well as potential new genes was generated.

Interpopulation patterns of divergence and selection across the transcriptome of the copepod Tigriopus californicus

The accumulation of genetic incompatibilities between isolated populations is thought to lead to the evolution of intrinsic postzygotic isolation. The molecular basis for these mechanisms, however, remains poorly understood. The intertidal copepod Tigriopus californicus provides unique opportunities for addressing mechanistic questions regarding the early stages of speciation; hybrids between highly divergent populations are fertile and viable, but exhibit reduced fitness at the F(2) or later generations. Given the current scarcity of genomic information in taxa at incipient stages of reproductive isolation, we utilize high-throughout 454 pyrosequencing to characterize a substantial fraction of protein-coding regions (the transcriptome) of T. californicus. Our sequencing effort was divided equally between two divergent populations in order to estimate levels of divergence and to reveal patterns of selection across the transcriptome. Assembly of sequences generated over 40,000 putatively unique transcripts (unigenes) for each population, 19,622 of which were orthologous between populations. BLAST searches of public databases determined protein identity and functional features for 15,402 and 12,670 unigenes, respectively. Based on rates of nonsynonymous and synonymous substitutions in 5897 interpopulation orthologs (those >150 bp and with at least 2X coverage), we identified 229 potential targets of positive selection. Many of these genes are predicted to be involved in several metabolic processes, and to function in hydrolase, peptidase and binding activities. The library of T. californicus coding regions, annotated with their predicted functions and level of divergence, will serve as an invaluable resource for elucidating molecular mechanisms underlying the early stages of speciation.