Revealing pancrustacean relationships: phylogenetic analysis of ribosomal protein genes places Collembola (springtails) in a monophyletic Hexapoda and reinforces the discrepancy between mitochondrial and nuclear DNA markers

The functional evolution of collembolan Ubx on the regulation of abdominal appendage formation

Folsomia candida is a tiny soil-living arthropod belonging to the Collembola, which is an outgroup to Insecta. It resembles insects as having a pair of antennae and three pairs of thorax legs, while it also possesses three abdominal appendages: a ventral tube located in the first abdominal segment (A1), a retinaculum in A3, and a furca in A4. Collembolan Ubx and AbdA specify abdominal appendages, but they are unable to repress appendage marker gene Dll. The genetic basis of collembolan appendage formation and the mechanisms by which Ubx and AbdA regulate Dll transcription and appendage development remains unknown. In this study, we analysed the developmental transcriptomes of F. candida and identified candidate appendage formation genes, including Ubx (FcUbx). The expression data revealed the dominance of Dll over Ubx during the embryonic 3.5 and 4.5 days, suggesting that Ubx is deficient in suppressing Dll at early appendage formation stages. Furthermore, via electrophoretic mobility shift assays and dual luciferase assays, we found that the binding and repression capacity of FcUbx on Drosophila Dll resembles those of the longest isoform of Drosophila Ubx (DmUbx_Ib), while the regulatory mechanism of the C-terminus of FcUbx on Dll repression is similar to that of the crustacean Artemia franciscana Ubx (AfUbx), demonstrating that the function of collembolan Ubx is intermediate between that of Insecta and Crustacea. In summary, our study provides novel insights into collembolan appendage formation and sheds light on the functional evolution of Ubx. Additionally, we propose a model that collembolan Ubx regulates abdominal segments in a context-specific manner.

Mitogenome selection in the evolution of key ecological strategies in the ancient hexapod class Collembola

A longstanding question in evolutionary biology is how natural selection and environmental pressures shape the mitochondrial genomic architectures of organisms. Mitochondria play a pivotal role in cellular respiration and aerobic metabolism, making their genomes functionally highly constrained. Evaluating selective pressures on mitochondrial genes can provide functional and ecological insights into the evolution of organisms. Collembola (springtails) are an ancient hexapod group that includes the oldest terrestrial arthropods in the fossil record, and that are closely associated with soil environments. Of interest is the diversity of habitat stratification preferences (life forms) exhibited by different species within the group. To understand whether signals of positive selection are linked to the evolution of life forms, we analysed 32 published Collembola mitogenomes in a phylomitogenomic framework. We found no evidence that signatures of selection are correlated with the evolution of novel life forms, but rather that mutations have accumulated as a function of time. Our results highlight the importance of nuclear-mitochondrial interactions in the evolution of collembolan life forms and that mitochondrial genomic data should be interpreted with caution, as complex selection signals may complicate evolutionary inferences.

Species-specific transcriptomic responses in Daphnia magna exposed to a bio-plastic production intermediate

Hydroxymethylfurfural (HMF) is a plant-based chemical building block that could potentially substitute petroleum-based equivalents, yet ecotoxicological data of this compound is currently limited. In this study, the effects of HMF on the reproduction and survival of Daphnia magna were assessed through validated ecotoxicological tests. The mechanism of toxicity was determined by analysis of transcriptomic responses induced by exposure to different concentrations of HMF using RNA sequencing. HMF exerted toxicity to D. magna with an EC₅₀ for effects on reproduction of 17.2 mg/l. HMF exposure affected molecular pathways including sugar and polysaccharide metabolism, lipid metabolism, general stress metabolism and red blood cell metabolism, although most molecular pathways affected by HMF exposure were dose specific. Hemoglobin genes, however, responded in a sensitive and dose-related manner. No induction of genes involved in the xenobiotic metabolism or oxidative stress metabolism pathway could be observed, which contrasted earlier observations on transcriptional responses of the terrestrial model Folsomia candida exposed to the same compound in a similar dose. We found 4189 orthologue genes between D. magna and F. candida, yet only twenty-one genes of those orthologues were co-regulated in both species. The contrasting transcriptional responses to the same compound exposed at a similar dose between D. magna and F. candida indicates limited overlap in stress responses among soil and aquatic invertebrates. The dose-related expression of hemoglobin provides further support for using hemoglobin expression as a biomarker for general stress responses in daphnids.

Multi-gene phylogeny and divergence estimations for Evaniidae (Hymenoptera)

Ensign wasps (Hymenoptera: Evaniidae) develop as predators of cockroach eggs (Blattodea), have a wide distribution and exhibit numerous interesting biological phenomena. The taxonomy of this lineage has been the subject of several recent, intensive efforts, but the lineage lacked a robust phylogeny. In this paper we present a new phylogeny, based on increased taxonomic sampling and data from six molecular markers (mitochondrial 16S and COI, and nuclear markers 28S, RPS23, CAD, and AM2), the latter used for the first time in phylogenetic reconstruction. Our intent is to provide a robust phylogeny that will stabilize and facilitate revision of the higher-level classification. We also show the continued utility of molecular motifs, especially the presence of an intron in the RPS23 fragments of certain taxa, to diagnose evaniid clades and assist with taxonomic classification. Furthermore, we estimate divergence times among evaniid lineages for the first time, using multiple fossil calibrations. Evaniidae radiated primarily in the Early Cretaceous (134.1-141.1 Mya), with and most extant genera diverging near the K-T boundary. The estimated phylogeny reveals a more robust topology than previous efforts, with the recovery of more monophyletic taxa and better higher-level resolution. The results facilitate a change in ensign wasp taxonomy, with Parevania, and Papatuka, syn. nov. becoming junior synonyms of Zeuxevania, and Acanthinevania, syn. nov. being designated as junior synonym of Szepligetella. We transfer 30 species to Zeuxevania, either reestablishing past combinations or as new combinations. We also transfer 20 species from Acanthinevania to Szepligetella as new combinations.

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.

Collembolan Transcriptomes Highlight Molecular Evolution of Hexapods and Provide Clues on the Adaptation to Terrestrial Life

Background

Collembola (springtails) represent a soil-living lineage of hexapods in between insects and crustaceans. Consequently, their genomes may hold key information on the early processes leading to evolution of Hexapoda from a crustacean ancestor.

Method

We assembled and annotated transcriptomes of the Collembola Folsomia candida and Orchesella cincta, and performed comparative analysis with protein-coding gene sequences of three crustaceans and three insects to identify adaptive signatures associated with the evolution of hexapods within the pancrustacean clade.

Results

Assembly of the springtail transcriptomes resulted in 37,730 transcripts with predicted open reading frames for F. candida and 32,154 for O. cincta, of which 34.2% were functionally annotated for F. candida and 38.4% for O. cincta. Subsequently, we predicted orthologous clusters among eight species and applied the branch-site test to detect episodic positive selection in the Hexapoda and Collembola lineages. A subset of 250 genes showed significant positive selection along the Hexapoda branch and 57 in the Collembola lineage. Gene Ontology categories enriched in these genes include metabolism, stress response (i.e. DNA repair, immune response), ion transport, ATP metabolism, regulation and development-related processes (i.e. eye development, neurological development).

Conclusions

We suggest that the identified gene families represent processes that have played a key role in the divergence of hexapods within the pancrustacean clade that eventually evolved into the most species-rich group of all animals, the hexapods. Furthermore, some adaptive signatures in collembolans may provide valuable clues to understand evolution of hexapods on land.

The ovary structure and oogenesis in the basal crustaceans and hexapods. Possible phylogenetic significance

Recent large-scale phylogenetic analyses of exclusively molecular or combined molecular and morphological characters support a close relationship between Crustacea and Hexapoda. The growing consensus on this phylogenetic link is reflected in uniting both taxa under the name Pancrustacea or Tetraconata. Several recent molecular phylogenies have also indicated that the monophyletic hexapods should be nested within paraphyletic crustaceans. However, it is still contentious exactly which crustacean taxon is the sister group to Hexapoda. Among the favored candidates are Branchiopoda, Malacostraca, Remipedia and Xenocarida (Remipedia + Cephalocarida). In this context, we review morphological and ultrastructural features of the ovary architecture and oogenesis in these crustacean groups in search of traits potentially suitable for phylogenetic considerations. We have identified a suite of morphological characters which may prove useful in further comparative studies.

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 complete mitochondrial genome of the black mud crab, Scylla serrata (Crustacea: Brachyura: Portunidae) and its phylogenetic position among (pan)crustaceans

The black mud crab, Scylla serrata (Forskål 1775), is the most economically important edible crab in South-East Asia. In the present study, the complete mitochondrial genome of black mud crab, S. serrata, was determined with the sequential polymerase chain reaction and primer walking sequencing. The complete mitochondrial genome was 15,721 bp in length with an A+T content of 69.2 % and contained 37 mitochondrial genes (13 protein coding genes (PCGs), 2 ribosomal RNA genes and 22 transfer RNA genes) and a control region (CR). The analysis of the CR sequence shows that it contains a multitude of repetitive fragments which can fold into hairpin-like or secondary structures and conserved elements as in other arthropods. The gene order of S. serrata mainly retains as the pancrustacean ground pattern, except for a single translocation of trnH. The gene arrangement of S. serrata appears to be a typical feature of portunid crabs. Phylogenetic analyses with concatenated amino acid sequences of 12 PCGs establishes that S. serrata in a well-supported monophyletic Portunidae and is consistent with previous morphological classification. Moreover, the phylogenomic results strongly support monophyletic Pancrustacea (Hexapoda plus "Crustaceans"). Within Pancrustacea, this study identifies Malacostraca + Entomostraca and Branchiopoda as the sister group to Hexapoda, which confirms that "Crustacea" is not monophyletic. Cirripedia + Remipedia appear to be a basal lineage of Pancrustacea. The present study also provides considerable data for the application of both population and phylogenetic studies of other crab species.

Advances in insect phylogeny at the dawn of the postgenomic era

Most species on Earth are insects and thus, understanding their evolutionary relationships is key to understanding the evolution of life. Insect relationships are increasingly well supported, due largely to technological advances in molecular sequencing and phylogenetic computational analysis. In this postgenomic era, insect systematics will be furthered best by integrative methods aimed at hypothesis corroboration from molecular, morphological, and paleontological evidence. This review of the current consensus of insect relationships provides a foundation for comparative study and offers a framework to evaluate incoming genomic evidence. Notable recent phylogenetic successes include the resolution of Holometabola, including the identification of the enigmatic Strepsiptera as a beetle relative and the early divergence of Hymenoptera; the recognition of hexapods as a crustacean lineage within Pancrustacea; and the elucidation of Dictyoptera orders, with termites placed as social cockroaches. Regions of the tree that require further investigation include the earliest winged insects (Palaeoptera) and Polyneoptera (orthopteroid lineages).

What is the phylogenetic signal limit from mitogenomes? The reconciliation between mitochondrial and nuclear data in the Insecta class phylogeny

BACKGROUND

Efforts to solve higher-level evolutionary relationships within the class Insecta by using mitochondrial genomic data are hindered due to fast sequence evolution of several groups, most notably Hymenoptera, Strepsiptera, Phthiraptera, Hemiptera and Thysanoptera. Accelerated rates of substitution on their sequences have been shown to have negative consequences in phylogenetic inference. In this study, we tested several methodological approaches to recover phylogenetic signal from whole mitochondrial genomes. As a model, we used two classical problems in insect phylogenetics: The relationships within Paraneoptera and within Holometabola. Moreover, we assessed the mitochondrial phylogenetic signal limits in the deeper Eumetabola dataset, and we studied the contribution of individual genes.

RESULTS

Long-branch attraction (LBA) artefacts were detected in all the datasets. Methods using Bayesian inference outperformed maximum likelihood approaches, and LBA was avoided in Paraneoptera and Holometabola when using protein sequences and the site-heterogeneous mixture model CAT. The better performance of this method was evidenced by resulting topologies matching generally accepted hypotheses based on nuclear and/or morphological data, and was confirmed by cross-validation and simulation analyses. Using the CAT model, the order Strepsiptera was recovered as sister to Coleoptera for the first time using mitochondrial sequences, in agreement with recent results based on large nuclear and morphological datasets. Also the Hymenoptera-Mecopterida association was obtained, leaving Coleoptera and Strepsiptera as the basal groups of the holometabolan insects, which coincides with one of the two main competing hypotheses. For the Paraneroptera, the currently accepted non-monophyly of Homoptera was documented as a phylogenetic novelty for mitochondrial data. However, results were not satisfactory when exploring the entire Eumetabola, revealing the limits of the phylogenetic signal that can be extracted from Insecta mitogenomes. Based on the combined use of the five best topology-performing genes we obtained comparable results to whole mitogenomes, highlighting the important role of data quality.

CONCLUSION

We show for the first time that mitogenomic data agrees with nuclear and morphological data for several of the most controversial insect evolutionary relationships, adding a new independent source of evidence to study relationships among insect orders. We propose that deeper divergences cannot be inferred with the current available methods due to sequence saturation and compositional bias inconsistencies. Our exploratory analysis indicates that the CAT model is the best dealing with LBA and it could be useful for other groups and datasets with similar phylogenetic difficulties.

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.

Mitogenome of the small abalone Haliotis diversicolor Reeve and phylogenetic analysis within Gastropoda

The complete mitochondrial coding regions of three small abalones Haliotis diversicolor Reeve, two collected from Vietnam and one from southern China, were successfully sequenced. The molecular feature of the mitochondrial genome is identical with the general description of the family Haliotidae mtDNAs that have been sequenced so far. The sequenced nucleotides are 16,186-16,266bp in length. The mitogenome encodes 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. Between adjacent genes trnH and nad5 there is an overlapping region. Comparison of the ratios of Ka and Ks among four species of Haliotis (H. diversicolor, H. discus hannai, H. rubra and H. tuberculata tuberculata) reveals that values of Ka/Ks in some NADH dehydrogenase and ATPase genes such as nad2, nad6 and atp8 are higher than those in other mitochondrial genes. Genome-wide gene arrangement among four species of Haliotis has been compared. Although all gene arrangement is the same in H. diversicolor, H. discus hannai and H. rubra, the location of trnS₂ and trnF in H. tuberculata tuberculata are inter-exchanged. Both gene arrangement and phylogenetic analysis support that the family Haliotidae is at a relatively primordial phylogenetic position in Gastropoda. Through alignment between Vietnam and southern China individuals, 111 SNPs are detected, most SNPs are synonymous mutations, and 7, 94 and 10 SNPs are observed at NCR, protein-coding region and RNA region, respectively. The result of SNP analysis also demonstrates that the difference is mainly in some NADH dehydrogenase genes between the Vietnam and southern China individuals.

A phylogenomic approach to resolve the arthropod tree of life

Arthropods were the first animals to conquer land and air. They encompass more than three quarters of all described living species. This extraordinary evolutionary success is based on an astoundingly wide array of highly adaptive body organizations. A lack of robustly resolved phylogenetic relationships, however, currently impedes the reliable reconstruction of the underlying evolutionary processes. Here, we show that phylogenomic data can substantially advance our understanding of arthropod evolution and resolve several conflicts among existing hypotheses. We assembled a data set of 233 taxa and 775 genes from which an optimally informative data set of 117 taxa and 129 genes was finally selected using new heuristics and compared with the unreduced data set. We included novel expressed sequence tag (EST) data for 11 species and all published phylogenomic data augmented by recently published EST data on taxonomically important arthropod taxa. This thorough sampling reduces the chance of obtaining spurious results due to stochastic effects of undersampling taxa and genes. Orthology prediction of genes, alignment masking tools, and selection of most informative genes due to a balanced taxa-gene ratio using new heuristics were established. Our optimized data set robustly resolves major arthropod relationships. We received strong support for a sister group relationship of onychophorans and euarthropods and strong support for a close association of tardigrades and cycloneuralia. Within pancrustaceans, our analyses yielded paraphyletic crustaceans and monophyletic hexapods and robustly resolved monophyletic endopterygote insects. However, our analyses also showed for few deep splits that were recently thought to be resolved, for example, the position of myriapods, a remarkable sensitivity to methods of analyses.

Loss of lipid synthesis as an evolutionary consequence of a parasitic lifestyle

Evolutionary loss of traits can result from negative selection on a specific phenotype, or if the trait is selectively neutral, because the phenotype associated with the trait has become redundant. Even essential traits may be lost, however, if the resulting phenotypic deficiencies can be compensated for by the environment or a symbiotic partner. Here we demonstrate that loss of an essential me-tabolic trait in parasitic wasps has evolved through environmental compensation. We tested 24 species for the ability to synthesize lipids de novo and collected additional data from the literature. We found the majority of adult parasitoid species to be incapable of synthesizing lipids, and phylogenetic analyses showed that the evolution of lack of lipogenesis is concurrent with that of parasitism in insects. Exploitive host manipulation, in which the host is forced to synthesize lipids to the benefit of the parasitoid, presumably facilitates loss of lipogenesis through environmental compensation. Lipogenesis re-evolved in a small number of parasitoid species, particularly host generalists. The wide range of host species in which generalists are able to develop may impede effective host manipulation and could have resulted in regaining of lipogenic ability in generalist parasitoids. As trait loss through environmental compensation is unnoticed at the phenotypic level, it may be more common than currently anticipated, especially in species involved in intricate symbiotic relationships with other species.

Parametric and non-parametric masking of randomness in sequence alignments can be improved and leads to better resolved trees

Background

Methods of alignment masking, which refers to the technique of excluding alignment blocks prior to tree reconstructions, have been successful in improving the signal-to-noise ratio in sequence alignments. However, the lack of formally well defined methods to identify randomness in sequence alignments has prevented a routine application of alignment masking. In this study, we compared the effects on tree reconstructions of the most commonly used profiling method (GBLOCKS) which uses a predefined set of rules in combination with alignment masking, with a new profiling approach (ALISCORE) based on Monte Carlo resampling within a sliding window, using different data sets and alignment methods. While the GBLOCKS approach excludes variable sections above a certain threshold which choice is left arbitrary, the ALISCORE algorithm is free of a priori rating of parameter space and therefore more objective.

Results

ALISCORE was successfully extended to amino acids using a proportional model and empirical substitution matrices to score randomness in multiple sequence alignments. A complex bootstrap resampling leads to an even distribution of scores of randomly similar sequences to assess randomness of the observed sequence similarity. Testing performance on real data, both masking methods, GBLOCKS and ALISCORE, helped to improve tree resolution. The sliding window approach was less sensitive to different alignments of identical data sets and performed equally well on all data sets. Concurrently, ALISCORE is capable of dealing with different substitution patterns and heterogeneous base composition. ALISCORE and the most relaxed GBLOCKS gap parameter setting performed best on all data sets. Correspondingly, Neighbor-Net analyses showed the most decrease in conflict.

Conclusions

Alignment masking improves signal-to-noise ratio in multiple sequence alignments prior to phylogenetic reconstruction. Given the robust performance of alignment profiling, alignment masking should routinely be used to improve tree reconstructions. Parametric methods of alignment profiling can be easily extended to more complex likelihood based models of sequence evolution which opens the possibility of further improvements.

400 million years on six legs: on the origin and early evolution of Hexapoda

Identifying the unambiguous sister group to the hexapods has been elusive. Traditional concepts include the Myriapoda (the Tracheata/Atelocerata hypothesis), but recent molecular studies consistently indicate it is the Crustacea, either in part or entirety (the Pancrustacea/Tetraconata hypothesis). The morphological evidence in support of Tracheata is reviewed, and most features are found to be ambiguous (i.e., losses, poorly known and surveyed structures, and probable convergences), though some appear to be synapomorphic, such as tentorial structure and the presence of styli and eversible vesicles. Other morphological features, particularly the structure of the eyes and nervous system, support Pancrustacea, as does consistent molecular evidence (which is reviewed and critiqued). Suggestions are made regarding hexapod-crustacean limb homologies. Relationships among basal (apterygote) hexapods are reviewed, and critical Paleozoic fossils are discussed. Despite the scarceness of Devonian hexapods, major lineages like Collembola and even dicondylic Insecta appeared in the Early Devonian; stem-group and putative Archaeognatha are known from the Carboniferous through Permian and the Late Devonian, respectively. Thus, the earliest divergences of hexapods were perhaps Late Silurian, which is considerably younger than several estimates made using molecular data.

Arthropod phylogeny: an overview from the perspectives of morphology, molecular data and the fossil record

Monophyly of Arthropoda is emphatically supported from both morphological and molecular perspectives. Recent work finds Onychophora rather than Tardigrada to be the closest relatives of arthropods. The status of tardigrades as panarthropods (rather than cycloneuralians) is contentious from the perspective of phylogenomic data. A grade of Cambrian taxa in the arthropod stem group includes gilled lobopodians, dinocaridids (e.g., anomalocaridids), fuxianhuiids and canadaspidids that inform on character acquisition between Onychophora and the arthropod crown group. A sister group relationship between Crustacea (itself likely paraphyletic) and Hexapoda is retrieved by diverse kinds of molecular data and is well supported by neuroanatomy. This clade, Tetraconata, can be dated to the early Cambrian by crown group-type mandibles. The rival Atelocerata hypothesis (Myriapoda+Hexapoda) has no molecular support. The basal node in the arthropod crown group is embroiled in a controversy over whether myriapods unite with chelicerates (Paradoxopoda or Myriochelata) or with crustaceans and hexapods (Mandibulata). Both groups find some molecular and morphological support, though Mandibulata is presently the stronger morphological hypothesis. Either hypothesis forces an unsampled ghost lineage for Myriapoda from the Cambrian to the mid Silurian.

Arthropod phylogeny revisited, with a focus on crustacean relationships

Higher-level arthropod phylogenetics is an intensely active field of research, not least as a result of the hegemony of molecular data. However, not all areas of arthropod phylogenetics have so far received equal attention. The application of molecular data to infer a comprehensive phylogeny of Crustacea is still in its infancy, and several emerging results are conspicuously at odds with morphology-based studies. In this study, we present a series of molecular phylogenetic analyses of 88 arthropods, including 57 crustaceans, representing all the major lineages, with Onychophora and Tardigrada as outgroups. Our analyses are based on published and new sequences for two mitochondrial markers, 16S rDNA and cytochrome c oxidase subunit I (COI), and the nuclear ribosomal gene 18S rDNA. We designed our phylogenetic analyses to assess the effects of different strategies of sequence alignment, alignment masking, nucleotide coding, and model settings. Our comparisons show that alignment optimization of ribosomal markers based on secondary structure information can have a radical impact on phylogenetic reconstruction. Trees based on optimized alignments recover monophyletic Arthropoda (excluding Onychophora), Pancrustacea, Malacostraca, Insecta, Myriapoda and Chelicerata, while Maxillopoda and Hexapoda emerge as paraphyletic groups. Our results are unable to resolve the highest-level relationships within Arthropoda, and none of our trees supports the monophyly of Myriochelata or Mandibulata. We discuss our results in the context of both the methodological variations between different analyses, and of recently proposed phylogenetic hypotheses. This article offers a preliminary attempt to incorporate the large diversity of crustaceans into a single molecular phylogenetic analysis, assessing the robustness of phylogenetic relationships under varying analysis parameters. It throws into sharp relief the relative strengths and shortcomings of the combined molecular data for assessing this challenging phylogenetic problem, and thereby provides useful pointers for future studies.

Higher-level crustacean phylogeny: consensus and conflicting hypotheses

This paper presents an overview of current hypotheses of higher-level crustacean phylogeny in order to assist and help focus further research. It concentrates on hypotheses proposed or debated in the recent literature based on morphological, molecular and combined evidence phylogenetic analyses. It can be concluded that crustacean phylogeny remains essentially unresolved. Conflict is rife, irrespective of whether one compares different morphological studies, molecular studies, or both. Using the number of recently proposed alternative sister group hypotheses for each of the major tetraconatan taxa as a rough estimate of phylogenetic uncertainty, it can be concluded that the phylogenetic position of Malacostraca remains the most problematic, closely followed by Branchiopoda, Cephalocarida, Remipedia, Ostracoda, Branchiura, Copepoda and Hexapoda. Future progress will depend upon a broader taxon sampling in molecular analyses, and the further exploration of new molecular phylogenetic markers. However, the need for continued revision and expansion of morphological datasets remains undiminished given the conspicuous lack of agreement between molecules and morphology for positioning several taxa. In view of the unparalleled morphological diversity of Crustacea, and the likely nesting of Hexapoda somewhere within Crustacea, working out a detailed phylogeny of Tetraconata is a crucial step towards understanding arthropod body plan evolution.

Molecular evolution of the crustacean hyperglycemic hormone family in ecdysozoans

BACKGROUND

Crustacean Hyperglycemic Hormone (CHH) family peptides are neurohormones known to regulate several important functions in decapod crustaceans such as ionic and energetic metabolism, molting and reproduction. The structural conservation of these peptides, together with the variety of functions they display, led us to investigate their evolutionary history. CHH family peptides exist in insects (Ion Transport Peptides) and may be present in all ecdysozoans as well. In order to extend the evolutionary study to the entire family, CHH family peptides were thus searched in taxa outside decapods, where they have been, to date, poorly investigated.

RESULTS

CHH family peptides were characterized by molecular cloning in a branchiopod crustacean, Daphnia magna, and in a collembolan, Folsomia candida. Genes encoding such peptides were also rebuilt in silico from genomic sequences of another branchiopod, a chelicerate and two nematodes. These sequences were included in updated datasets to build phylogenies of the CHH family in pancrustaceans. These phylogenies suggest that peptides found in Branchiopoda and Collembola are more closely related to insect ITPs than to crustacean CHHs. Datasets were also used to support a phylogenetic hypothesis about pancrustacean relationships, which, in addition to gene structures, allowed us to propose two evolutionary scenarios of this multigenic family in ecdysozoans.

CONCLUSIONS

Evolutionary scenarios suggest that CHH family genes of ecdysozoans originate from an ancestral two-exon gene, and genes of arthropods from a three-exon one. In malacostracans, the evolution of the CHH family has involved several duplication, insertion or deletion events, leading to neuropeptides with a wide variety of functions, as observed in decapods. This family could thus constitute a promising model to investigate the links between gene duplications and functional divergence.

Arthropod relationships revealed by phylogenomic analysis of nuclear protein-coding sequences

The remarkable antiquity, diversity and ecological significance of arthropods have inspired numerous attempts to resolve their deep phylogenetic history, but the results of two decades of intensive molecular phylogenetics have been mixed. The discovery that terrestrial insects (Hexapoda) are more closely related to aquatic Crustacea than to the terrestrial centipedes and millipedes (Myriapoda) was an early, if exceptional, success. More typically, analyses based on limited samples of taxa and genes have generated results that are inconsistent, weakly supported and highly sensitive to analytical conditions. Here we present strongly supported results from likelihood, Bayesian and parsimony analyses of over 41 kilobases of aligned DNA sequence from 62 single-copy nuclear protein-coding genes from 75 arthropod species. These species represent every major arthropod lineage, plus five species of tardigrades and onychophorans as outgroups. Our results strongly support Pancrustacea (Hexapoda plus Crustacea) but also strongly favour the traditional morphology-based Mandibulata (Myriapoda plus Pancrustacea) over the molecule-based Paradoxopoda (Myriapoda plus Chelicerata). In addition to Hexapoda, Pancrustacea includes three major extant lineages of 'crustaceans', each spanning a significant range of morphological disparity. These are Oligostraca (ostracods, mystacocarids, branchiurans and pentastomids), Vericrustacea (malacostracans, thecostracans, copepods and branchiopods) and Xenocarida (cephalocarids and remipedes). Finally, within Pancrustacea we identify Xenocarida as the long-sought sister group to the Hexapoda, a result confirming that 'crustaceans' are not monophyletic. These results provide a statistically well-supported phylogenetic framework for the largest animal phylum and represent a step towards ending the often-heated, century-long debate on arthropod relationships.

The complete mitochondrial genome of Symphylella sp. (Myriapoda: Symphyla): Extensive gene order rearrangement and evidence in favor of Progoneata

We determined the complete 14,667bp mitochondrial DNA sequence of Symphylella sp., the first representative of the Scolopendrellidae (Arthropoda: Myriapoda: Symphyla). With respect to the ancestral arthropod mitochondrial gene order, two protein-coding genes, the rRNAs and 10 of the tRNAs appear to be rearranged. This rearrangement is novel in the arthropods and genes with identical transcriptional polarity are clustered except for trnE, trnN and putative control region (CR), resembling two previously reported diplopod genomes. A duplication/loss (random and non-random)-recombination model was proposed to account for the generation of the gene order in Symphylella sp. All phylogenetic analysis yielded strong support for a clade of Symphyla plus Diplopoda, i.e., Progoneata. However, the phylogenetic position of Myriapoda within Arthropoda remains unclear. The amino acid dataset gives strong support for an affinity to Pancrustacea, while the nucleotide dataset weakly supports Myriapoda grouped with Chelicerata.