BACKGROUND
general hypotheses on arthropod phylogeny are rapidly being altered by dna sequence data  <cit> .
for instance, the atelocerata concept held that hexapods and myriapods are united in one clade, but under the influence of molecular data  this concept was replaced by the view that crustaceans and hexapods constitute a monophyletic group, which is known as pancrustacea .
another recently proposed, but still highly debated viewpoint is the diphyletic origin of hexapoda, which was initially raised by nardi and co-workers in  <dig>  <cit> .
based on four mitochondrial genes, they  <cit>  observed that two species of collembola  branched off before the other pancrustacean groups that were included in their study , suggesting paraphyly of hexapoda.
their thesis was that the six-legged body plan of collembola and other hexapods evolved at least twice: once in the group of wingless hexapods and another time in the true insects.
the conclusions of nardi et al.
<cit>  resulted in a vivid scientific debate, and many studies have addressed the phylogenetic placement of collembola since then.
some authors focused on mitochondrial sequences, others analyzed nuclear genes.
additional mitochondrial sequences confirmed that, due to the placement of collembola separate from the other hexapods, hexapoda are indeed not monophyletic  <cit> .
however, after thorough analyses exploring the effects of outgroup and gene choice, sequence handling and optimality criteria on inferred trees, cameron and co-workers  <cit>  concluded that the mitochondrial data as available at the time were inadequate to fully resolve hexapod relationships  <cit> .
hassinin  <cit>  arrived at a similar conclusion in a more recent study focusing on the effects of reverse strand-bias.
most recently, carapelli and co-workers  <cit>  reported new analyses on a very large dataset, consisting of no fewer than a hundred almost-complete mitochondrial genomes.
these new analyses, which were based on a novel model of amino acid sequence evolution , supported the non-monophyly of hexapod groups.
it has gradually become clear in pancrustacean phylogeny that nuclear and mitochondrial datasets tell different stories, and often result in different conclusions  <cit> .
remarkably, studies that addressed the question using nuclear genomic data  indicate that the collembola group between crustaceans and insects and that hexapoda is monophyletic  <cit> .
however, most of those studies included a relatively small number of loci  <cit> , most likely because obtaining data on protein-encoding dna sequences is not always straightforward for groups for which little genomic information is available.
here we try to fill this gap by re-evaluating the position of collembola using a relatively large number of nuclear protein-encoding sequences that are, although all for ribosomal proteins, assumed to be distributed throughout the genome .
several authors have shown that publicly available data can be useful when conducting a large-scale phylogenetic study , and that expressed sequence tags  can be extremely valuable for phylogenetic purposes  <cit> .
here, we combine data from a recently finished est sequencing project on the collembolan folsomia candida  <cit> , with data on  <dig> ecdysozoan species  available in the public genbank repository  <cit> , and with data from a smaller est dataset of the collembolan orchesella cincta.
we focus on ribosomal proteins to prevent the problem of analyzing paralogous genes .
RESULTS
in total, gene-sequences for  <dig> ribosomal proteins were obtained from the folsomia candida est dataset.
this is almost two-thirds of the total set of  <dig> ribosomal proteins  <cit>  found in the genome of drosophila melanogaster.
four d. melanogaster ribosomal protein sequences  showed high similarity with two, instead of one f. candida transcript cluster in the est dataset.
comparison of the f. candida transcripts with those of d. melanogaster revealed insertions/deletions resulting in frame shifts in one of the two f. candida est clusters for rpl <dig>  rpl <dig> and rpl <dig>  transcripts with a frame-shift were discarded.
two highly diverse f. candida est clusters  showed homology with d. melanogaster rpl <dig>  the f. candida rpl <dig> singleton sequence was excluded from further analysis.
the discarded rpl <dig>  rpl <dig>  rpl <dig> and rpl <dig> transcripts may stem from duplications in the f. candida genome , or from constitutively expressed pseudogenes.
this situation may be analogous to the apparent amplification of many mammalian ribosomal proteins; for instance, the human genome contains over  <dig> ribosomal protein pseudogenes  <cit> .
still, it seems that only one copy of each ribosomal protein is actually functional  <cit> .
as described in the methods section, the remaining  <dig> ribosomal protein sequences were used to retrieve ribosomal protein sequence information on  <dig> additional ecdysozoan species.
in addition, ribosomal protein sequences of d. melanogaster, apis mellifera and caenorhabditis elegans were retrieved from the ribosomal protein gene-database .
the number of usable  ribosomal protein gene sequences that were obtained per species ranged from two  to  <dig>  .
redundancy for a given ribosomal protein gene in a given species was often low, and many gene sequences were represented by one or a few est sequences only.
it should be mentioned that due to this rather low sequence coverage the dataset is most probably not free from sequencing errors.
furthermore, none of the  <dig> ribosomal proteins that were included in the dataset were observed in all of the  <dig> species investigated .
in summary, for each ribosomal protein gene information was available for 25% to 86% of the taxa.
left: the inferred ribosomal protein sequences  that were included in the concatenated dataset.
only ribosomal proteins present in the folsomia candida est dataset were included in the analysis.
occurrence*: the number of species for which was data available for a certain ribosomal protein.
# of variable sites: positions that constitute more than one amino acid in the different alignments.
numbers were calculated in mega <dig>  <cit> .
right: all species that were included in the analyses.
species marked with ** were excluded from the final analyses because they contained too few sequences.
occurrence***: the number of genes  that were available for a specific taxon.
****: rps27a and rpl <dig> are fused to ubiquitin  <cit> .
we calculated "similarity" values between the amino acid sequences of f. candida and three well-represented species , daphnia magna and a. mellifera).
these values were mapped onto a ternary graph .
almost all points cluster in the lower region of the ternary graph, showing that for almost all genes the distance between f. candida and c. elegans is greater than the distance between f. candida and a. mellifera.
the graph also shows that most genes of f. candida are more "similar" to a. mellifera, while some have more in common with d. magna.
the individual alignments were concatenated and phylogenetic analyses were conducted to investigate the position of collembola.
two species were excluded from the analyses .
the final alignment had a length of  <dig> inferred amino acids, representing in total  <dig>  nucleotides.
information was available for 66% of the amino acid positions.
likelihood mapping was applied to obtain estimates of phylogenetic signal.
the concatenated dataset contained more phylogenetic signal  than each of the independent ribosomal protein alignments .
the trees obtained by the different tree-reconstruction algorithms were highly comparable .
in all reconstructions , chelicerata and pancrustacea each formed a monophyletic group, with relatively high support .
the two branchiopods included in this study  grouped together, and remained separate from the other crustaceans .
the relationships within the insecta were weakly resolved; however, diptera was recovered as a monophyletic clade, as were lepidoptera and coleoptera .
however, the hemiptera were resolved as a paraphyletic group.
homalodisca coagulata grouped with locusta migratoria , rather than with the other hemipterans acyrthosypon pisum and toxoptera citricida.
the highly supported, but obviously incorrect, positioning of homalodisca coagulata does not seem to be an artefact of the method that allowed for missing data, since all three hemiptera, as well as locusta migratoria, were represented by a large number of ribosomal protein gene sequences .
the incorrect placement of h. coagulata could be a consequence of the inability of ribosomal protein genes to resolve more recent evolutionary splits, which may be a trade-off of their suitability for deeper phylogenies.
hexapoda was clearly monophyletic: both collembolans  grouped together and formed the sister-group to the insecta in all analyses conducted .
several c. elegans and d. melanogaster ribosomal proteins are duplicated .
mp analyses of a second concatenated dataset that contained d. melanogaster homologs for ribosomal protein rps <dig>  rps15a, rps <dig>  rps <dig>  rpl <dig> and rpl10a resulted in a similar topology .
discussion
in this study we reassessed the position of collembola, using  genes for  <dig> nuclear encoded proteins.
the main result of our study is clear evidence of monophyly of hexapoda.
all phylogenetic reconstruction methods employed in this study support this hypothesis .
based on our nuclear dataset we conclude that the six-legged body plan, as found among insects and collembola, evolved only once in the course of evolution.
this is in contrast to results obtained using large mitochondrial datasets  <cit>  that by and large suggest that the characteristic hexapod body plan was acquired in parallel by collembola and insects due to convergent evolution, rather than by descent.
discrepancies between pancrustacean relationships as revealed by either nuclear or mitochondrial datasets seem almost universal.
it is of major importance to focus on the causes of these discrepancies, and whether or not one of the two types of markers is superior.
elaborate discussions on the 'pros and cons' of one or both of the two different markers, and the possible approaches on how to correct for ambiguous signals are given in several recent papers  <cit> .
comparative studies that contrast nuclear and mitochondrial datasets suggest that nuclear markers are preferred in deep arthropodan molecular phylogenetics, as mitochondrial genes tend to be more substitutionally biased and evolve  in a much faster way  <cit> .
already in  <dig> curole and kocher  <cit>  stated in a review paper that the value of mitochondrial genes in deep-level phylogeny is debatable and that "controversial" mitochondrial dna  results should be verified with nuclear encoded genes.
this was also the final conclusion of springer and co-workers  <cit> .
these authors compared the usability of nuclear and mitochondrial encoded genes in inferring deep-level mammalian phylogenies.
the authors report that nuclear encoded genes  outperform mitochondrial markers in resolving deep splits.
springer and co-workers suggest that the reason for this dissimilarity in resolving-power might be found, among others, in the rate of nucleotide substitution  <cit> .
still, although the nuclear protein-encoding sequences in the study of springer et al.
<cit>  outperformed the mitochondrial genes, mtdna-based studies are not necessarily useless for deep phylogeny.
they are only problematic if mitochondrial genomes evolve at such a rate that saturation of substitutions makes actual phylogenetic signals from deeper nodes hazy  <cit> .
otherwise, analyses using appropriate models should still be able to retrieve a plausible tree  <cit> .
in a recent study, kjer and honeycutt  <cit>  used an approach that included all data found in mitochondrial genomes .
after applying a site-specific rate model, these authors retrieved a phylogenetic tree of mammals that was in accordance with recent nuclear dna based phylogenies  <cit> .
when investigating cheliceratan relationships jones et al.
<cit>  arrived at a comparable conclusion.
these authors state that mtdna can be applied in molecular phylogenetics, but only when an appropriate substitution model  is used.
these authors state as a final remark that earlier mtdna studies that focused on deep-phylogenetic questions should be thoroughly re-evaluated  <cit> .
however, such models of mitochondrial sequence evolution might first need to be developed before collembola can be placed with certainty in the arthropod phylogenetic tree.
as mentioned before, carapelli and co-workers  <cit>  investigated an innovative pancrustacean-model of mitochondrial protein change.
this model significantly aided the tree building, but did not yield a monophyletic hexapoda  <cit> .
an advantage of ribosomal protein genes is that the sequences of different species can be relatively easily homologized due to their conserved nature.
however, there are also disadvantages.
although ribosomal protein genes are distributed all over the genome, they definitely do not evolve independently.
coevolving sites are known to exist in ribosomal proteins  <cit> .
for example, amino acid residues that are near trna binding sites in the ribosome appear to evolve in a related manner  <cit> .
it has to be mentioned that we included only two collembola in our analyses.
preferably, more springtail species, and maybe even more importantly, proturans and diplurans, should be included.
those latter basal hexapod groups were excluded from the current analysis as they lack available  data.
while earlier work suggests that proturan and dipluran genes might be fairly divergent from other arthropods  <cit> , this and other papers  suggest that it should be relatively easy to obtain phylogenetically relevant sequence information on those groups by est sequencing.
another intriguing result of this study is the non-monophyly of the crustaceans.
the branchiopods d. magna and a. franciscana clustered with the hexapods rather than with the other crustaceans in the malacostracan group.
this is in accordance with studies by regier and co-workers  <cit>  and mallat and giribet  <cit> , which suggests that the hexapod lineage evolved from within the crustaceans  <cit> .
the observed close relationship between hexapods and branchiopods, in combination with some other characteristics, made glenner and co-workers  <cit>  suggest that branchiopod groups colonized terrestrial ecosystems as insects.
as a final remark we would like to point out that this study shows that collembola occupy a crucial position.
obtaining additional  sequence information on collembola, as well as other basal hexapods  will definitely result in a better understanding of the phylogenetic origin of insects.
CONCLUSIONS
the phylogenetic efforts presented here clearly show that collembola is a sister group of insecta .
our results reinforce the discrepancy between results obtained using mitochondrial and nuclear datasets.
it seems of major importance to unravel the underlying causes of the disagreements observed, or otherwise focus on nuclear encoded genes.

 general hypotheses on arthropod phylogeny are rapidly being altered by dna sequence data. for instance, the atelocerata concept held that hexapods and myriapods are united in one clade, but under the influence of molecular data this concept was replaced by the view that crustaceans and hexapods constitute a monophyletic group, which is known as pancrustacea. another recently proposed, but still highly debated viewpoint is the diphyletic origin of hexapoda, which was initially raised by nardi and co-workers in. based on four mitochondrial genes, they observed that two species of collembola branched off before the other pancrustacean groups that were included in their study, suggesting paraphyly of hexapoda. their thesis was that the six-legged body plan of collembola and other hexapods evolved at least twice: once in the group of wingless hexapods and another time in the true insects. the conclusions of nardi et al. resulted in a vivid scientific debate, and many studies have addressed the phylogenetic placement of collembola since then. some authors focused on mitochondrial sequences, others analyzed nuclear genes. additional mitochondrial sequences confirmed that, due to the placement of collembola separate from the other hexapods, hexapoda are indeed not monophyletic. however, after thorough analyses exploring the effects of outgroup and gene choice, sequence handling and optimality criteria on inferred trees, cameron and co-workers concluded that the mitochondrial data as available at the time were inadequate to fully resolve hexapod relationships. hassinin arrived at a similar conclusion in a more recent study focusing on the effects of reverse strand-bias. most recently, carapelli and co-workers reported new analyses on a very large dataset, consisting of no fewer than a hundred almost-complete mitochondrial genomes. these new analyses, which were based on a novel model of amino acid sequence evolution, supported the non-monophyly of hexapod groups. it has gradually become clear in pancrustacean phylogeny that nuclear and mitochondrial datasets tell different stories, and often result in different conclusions. remarkably, studies that addressed the question using nuclear genomic data indicate that the collembola group between crustaceans and insects and that hexapoda is monophyletic. however, most of those studies included a relatively small number of loci, most likely because obtaining data on protein-encoding dna sequences is not always straightforward for groups for which little genomic information is available. here we try to fill this gap by re-evaluating the position of collembola using a relatively large number of nuclear protein-encoding sequences that are, although all for ribosomal proteins, assumed to be distributed throughout the genome. several authors have shown that publicly available data can be useful when conducting a large-scale phylogenetic study, and that expressed sequence tags can be extremely valuable for phylogenetic purposes. here, we combine data from a recently finished est sequencing project on the collembolan folsomia candida, with data on ecdysozoan species available in the public genbank repository, and with data from a smaller est dataset of the collembolan orchesella cincta. we focus on ribosomal proteins to prevent the problem of analyzing paralogous genes.  in total, gene-sequences for ribosomal proteins were obtained from the folsomia candida est dataset. this is almost two-thirds of the total set of ribosomal proteins found in the genome of drosophila melanogaster. four d. melanogaster ribosomal protein sequences showed high similarity with two, instead of one f. candida transcript cluster in the est dataset. comparison of the f. candida transcripts with those of d. melanogaster revealed insertions/deletions resulting in frame shifts in one of the two f. candida est clusters for rpl rpl and rpl transcripts with a frame-shift were discarded. two highly diverse f. candida est clusters showed homology with d. melanogaster rpl the f. candida rpl singleton sequence was excluded from further analysis. the discarded rpl rpl rpl and rpl transcripts may stem from duplications in the f. candida genome, or from constitutively expressed pseudogenes. this situation may be analogous to the apparent amplification of many mammalian ribosomal proteins; for instance, the human genome contains over ribosomal protein pseudogenes. still, it seems that only one copy of each ribosomal protein is actually functional. as described in the methods section, the remaining ribosomal protein sequences were used to retrieve ribosomal protein sequence information on additional ecdysozoan species. in addition, ribosomal protein sequences of d. melanogaster, apis mellifera and caenorhabditis elegans were retrieved from the ribosomal protein gene-database. the number of usable ribosomal protein gene sequences that were obtained per species ranged from two to. redundancy for a given ribosomal protein gene in a given species was often low, and many gene sequences were represented by one or a few est sequences only. it should be mentioned that due to this rather low sequence coverage the dataset is most probably not free from sequencing errors. furthermore, none of the ribosomal proteins that were included in the dataset were observed in all of the species investigated. in summary, for each ribosomal protein gene information was available for 25% to 86% of the taxa. left: the inferred ribosomal protein sequences that were included in the concatenated dataset. only ribosomal proteins present in the folsomia candida est dataset were included in the analysis. occurrence*: the number of species for which was data available for a certain ribosomal protein. # of variable sites: positions that constitute more than one amino acid in the different alignments. numbers were calculated in mega. right: all species that were included in the analyses. species marked with ** were excluded from the final analyses because they contained too few sequences. occurrence***: the number of genes that were available for a specific taxon. ****: rps27a and rpl are fused to ubiquitin. we calculated "similarity" values between the amino acid sequences of f. candida and three well-represented species, daphnia magna and a. mellifera). these values were mapped onto a ternary graph. almost all points cluster in the lower region of the ternary graph, showing that for almost all genes the distance between f. candida and c. elegans is greater than the distance between f. candida and a. mellifera. the graph also shows that most genes of f. candida are more "similar" to a. mellifera, while some have more in common with d. magna. the individual alignments were concatenated and phylogenetic analyses were conducted to investigate the position of collembola. two species were excluded from the analyses. the final alignment had a length of inferred amino acids, representing in total nucleotides. information was available for 66% of the amino acid positions. likelihood mapping was applied to obtain estimates of phylogenetic signal. the concatenated dataset contained more phylogenetic signal than each of the independent ribosomal protein alignments. the trees obtained by the different tree-reconstruction algorithms were highly comparable. in all reconstructions, chelicerata and pancrustacea each formed a monophyletic group, with relatively high support. the two branchiopods included in this study grouped together, and remained separate from the other crustaceans. the relationships within the insecta were weakly resolved; however, diptera was recovered as a monophyletic clade, as were lepidoptera and coleoptera. however, the hemiptera were resolved as a paraphyletic group. homalodisca coagulata grouped with locusta migratoria, rather than with the other hemipterans acyrthosypon pisum and toxoptera citricida. the highly supported, but obviously incorrect, positioning of homalodisca coagulata does not seem to be an artefact of the method that allowed for missing data, since all three hemiptera, as well as locusta migratoria, were represented by a large number of ribosomal protein gene sequences. the incorrect placement of h. coagulata could be a consequence of the inability of ribosomal protein genes to resolve more recent evolutionary splits, which may be a trade-off of their suitability for deeper phylogenies. hexapoda was clearly monophyletic: both collembolans grouped together and formed the sister-group to the insecta in all analyses conducted. several c. elegans and d. melanogaster ribosomal proteins are duplicated. mp analyses of a second concatenated dataset that contained d. melanogaster homologs for ribosomal protein rps rps15a, rps rps rpl and rpl10a resulted in a similar topology. discussion in this study we reassessed the position of collembola, using genes for nuclear encoded proteins. the main result of our study is clear evidence of monophyly of hexapoda. all phylogenetic reconstruction methods employed in this study support this hypothesis. based on our nuclear dataset we conclude that the six-legged body plan, as found among insects and collembola, evolved only once in the course of evolution. this is in contrast to results obtained using large mitochondrial datasets that by and large suggest that the characteristic hexapod body plan was acquired in parallel by collembola and insects due to convergent evolution, rather than by descent. discrepancies between pancrustacean relationships as revealed by either nuclear or mitochondrial datasets seem almost universal. it is of major importance to focus on the causes of these discrepancies, and whether or not one of the two types of markers is superior. elaborate discussions on the'pros and cons' of one or both of the two different markers, and the possible approaches on how to correct for ambiguous signals are given in several recent papers. comparative studies that contrast nuclear and mitochondrial datasets suggest that nuclear markers are preferred in deep arthropodan molecular phylogenetics, as mitochondrial genes tend to be more substitutionally biased and evolve in a much faster way. already in curole and kocher stated in a review paper that the value of mitochondrial genes in deep-level phylogeny is debatable and that "controversial" mitochondrial dna results should be verified with nuclear encoded genes. this was also the final conclusion of springer and co-workers. these authors compared the usability of nuclear and mitochondrial encoded genes in inferring deep-level mammalian phylogenies. the authors report that nuclear encoded genes outperform mitochondrial markers in resolving deep splits. springer and co-workers suggest that the reason for this dissimilarity in resolving-power might be found, among others, in the rate of nucleotide substitution. still, although the nuclear protein-encoding sequences in the study of springer et al. outperformed the mitochondrial genes, mtdna-based studies are not necessarily useless for deep phylogeny. they are only problematic if mitochondrial genomes evolve at such a rate that saturation of substitutions makes actual phylogenetic signals from deeper nodes hazy. otherwise, analyses using appropriate models should still be able to retrieve a plausible tree. in a recent study, kjer and honeycutt used an approach that included all data found in mitochondrial genomes. after applying a site-specific rate model, these authors retrieved a phylogenetic tree of mammals that was in accordance with recent nuclear dna based phylogenies. when investigating cheliceratan relationships jones et al. arrived at a comparable conclusion. these authors state that mtdna can be applied in molecular phylogenetics, but only when an appropriate substitution model is used. these authors state as a final remark that earlier mtdna studies that focused on deep-phylogenetic questions should be thoroughly re-evaluated. however, such models of mitochondrial sequence evolution might first need to be developed before collembola can be placed with certainty in the arthropod phylogenetic tree. as mentioned before, carapelli and co-workers investigated an innovative pancrustacean-model of mitochondrial protein change. this model significantly aided the tree building, but did not yield a monophyletic hexapoda. an advantage of ribosomal protein genes is that the sequences of different species can be relatively easily homologized due to their conserved nature. however, there are also disadvantages. although ribosomal protein genes are distributed all over the genome, they definitely do not evolve independently. coevolving sites are known to exist in ribosomal proteins. for example, amino acid residues that are near trna binding sites in the ribosome appear to evolve in a related manner. it has to be mentioned that we included only two collembola in our analyses. preferably, more springtail species, and maybe even more importantly, proturans and diplurans, should be included. those latter basal hexapod groups were excluded from the current analysis as they lack available data. while earlier work suggests that proturan and dipluran genes might be fairly divergent from other arthropods, this and other papers suggest that it should be relatively easy to obtain phylogenetically relevant sequence information on those groups by est sequencing. another intriguing result of this study is the non-monophyly of the crustaceans. the branchiopods d. magna and a. franciscana clustered with the hexapods rather than with the other crustaceans in the malacostracan group. this is in accordance with studies by regier and co-workers and mallat and giribet, which suggests that the hexapod lineage evolved from within the crustaceans. the observed close relationship between hexapods and branchiopods, in combination with some other characteristics, made glenner and co-workers suggest that branchiopod groups colonized terrestrial ecosystems as insects. as a final remark we would like to point out that this study shows that collembola occupy a crucial position. obtaining additional sequence information on collembola, as well as other basal hexapods will definitely result in a better understanding of the phylogenetic origin of insects.  the phylogenetic efforts presented here clearly show that collembola is a sister group of insecta. our results reinforce the discrepancy between results obtained using mitochondrial and nuclear datasets. it seems of major importance to unravel the underlying causes of the disagreements observed, or otherwise focus on nuclear encoded genes. 
in recent years, several new hypotheses on phylogenetic relations among arthropods have been proposed on the basis of dna sequences. in total ribosomal proteins were obtained for the collembolan folsomia candida. although ribosomal proteins in general may not evolve independently, they once more appear highly valuable for phylogenetic reconstruction. here, we re-evaluate the position of collembola using ribosomal protein gene sequences. each ribosomal protein gene was available for 25% to 86% of the taxa. this underpins the inconsistency between nuclear and mitochondrial datasets when analyzing pancrustacean relationships. caution is needed when applying mitochondrial markers in deep phylogeny.
 in recent years, several new hypotheses on phylogenetic relations among arthropods have been proposed on the basis of dna sequences. one of the challenged hypotheses is the monophyly of hexapods. this discussion originated from analyses based on mitochondrial dna datasets that, due to an unusual positioning of collembola, suggested that the hexapod body plan evolved at least twice. here, we re-evaluate the position of collembola using ribosomal protein gene sequences.  in total ribosomal proteins were obtained for the collembolan folsomia candida. these sequences were aligned with sequence data on other ecdysozoans. each ribosomal protein gene was available for 25% to 86% of the taxa. however, the total sequence information was unequally distributed over the taxa and ranged between 4% and 100%. a concatenated dataset was constructed, of which ~66% of the positions were filled. phylogenetic tree reconstructions, using maximum likelihood, maximum parsimony, and bayesian methods, resulted in a topology that supports monophyly of hexapoda.  although ribosomal proteins in general may not evolve independently, they once more appear highly valuable for phylogenetic reconstruction. our analyses clearly suggest that hexapoda is monophyletic. this underpins the inconsistency between nuclear and mitochondrial datasets when analyzing pancrustacean relationships. caution is needed when applying mitochondrial markers in deep phylogeny. 