BACKGROUND
hemocyanins are large copper-proteins that transport o <dig> in the hemolymph of many arthropods and mollusks  <cit> .
however, the hemocyanins of these two phyla are structurally different and emerged independently  <cit> .
hemocyanins evolved early in the arthropod stem lineage from the phenoloxidases, which are o2-consuming enzymes involved in the melanin pathway  <cit> .
other members of the arthropod hemocyanin superfamily have lost the ability to bind copper and thus o <dig>  and gave rise to the non-respiratory pseudo-hemocyanins  in decapod crustaceans and the hexamerins in hexapods, which serve as storage proteins  <cit> .
arthropod hemocyanins form hexamers or oligo-hexamers of identical or related subunits with a molecular mass of about  <dig> kda  <cit> .
in each subunit, o2-binding is mediated by two cu+ ions, which are coordinated by six histidine residues .
based on biochemical, immunochemical and molecular phylogenetic analyses, distinct hemocyanin subunit types have been identified in chelicerata, myriapoda, crustacea and hexapoda  <cit> .
these subunits experienced an independent evolution within each of these taxa, with the exception of a more complex pattern within the pancrustacea   <cit> .
within the chelicerates, biochemical analyses have demonstrated the presence of hemocyanins in xiphosura , scorpiones, uropygi , amblypygi , and araneae , but failed to identify these respiratory proteins in pycnogonida , solifugae  and acari   <cit> .
complete sets of hemocyanin subunit sequences are available from the tarantula eurypelma californicum   <cit> , the hunting spider cupiennius salei  <cit> , the golden orb web spider nephila inaurata  <cit>  and the horse shoe crab carcinoscorpius rotundicauda .
chelicerate hemocyanins are composed of up to eight distinct subunit types and form either  <dig> ×  <dig>   <dig> ×  <dig>   <dig> ×  <dig> or  <dig> × 6mers  <cit> .
each subunit type occupies a distinct position within the native oligomer  <cit> .
the subunits have similar oxygen binding properties, but different physico-chemical characteristics  <cit>  and evolutionary origins  <cit> .
xiphosurans  have the largest hemocyanin molecules known, consisting of  <dig>  subunits and up to eight distinct subunits types in limulus polyphemus   <cit> .
scorpiones, amblypygi, uropygi, and some araneae have  <dig> × 6mer hemocyanins.
the  <dig> × 6mer hemocyanin of the tarantula e. californicum is the best studied example and comprises seven distinct subunit types   <cit> .
a similar subunit composition was found in many other araneae, the amblypygi and the uropygi  <cit> .
among the araneae, variations from this "standard" scheme have been found in the entelegyne spiders of the rta-clade .
in this large taxon,  <dig> × 6mer and  <dig> × 6mer hemocyanins occur that consist only of g-type subunits, but have lost the other six subunit types  present in other araneae  <cit> .
scorpion hemocyanins are composed of eight subunit types named  <dig>  3a, 3b, 3c,  <dig>  5a and 5b  <cit> .
immunological and structural studies have suggested the correspondence between the distinct araneaen, scorpion and limulus hemocyanin subunits  <cit> .
this indicates that the evolution of distinct subunit types preceded the separation of xiphosura and arachnida.
to understand the evolution of the complex oligomeric structures of chelicerate hemocyanins, we have obtained  <dig> novel cdna sequences of hemocyanin subunits from a sea spider, a horseshoe crab, a whip scorpion and a whip spider.
together with the previously sequenced hemocyanins from horseshoe crabs, scorpions and spiders, and those assembled from expressed sequence tags ,  <dig> full length chelicerate subunit sequences are available.
these data allow us i. to trace the hemocyanin subunit evolution, ii.
to reconstruct the emergence of the hemocyanin oligomers, and iii.
to infer chelicerate phylogeny and divergence times.
methods
sequencing of chelicerate hemocyanin cdna
full length coding sequences were available for the hemocyanins of the tarantula e. californicum  <cit> , the mangrove horseshoe crab c. rotundicauda, the golden orb web spider n. inaurata  <cit>  and the hunting spider c. salei  <cit>  .
in addition, hemocyanin primary structures had been obtained by conventional amino acid sequencing from the horseshoe crab tachypleus tridentatus   <cit>  and the scorpion androctonus australis   <cit> .
<dig>  ests were generated from total rna of the sea spider endeis spinosa  as described before  <cit> .
gene ontology assessment and blast searches identified six ests with significant similarities to arthropod hemocyanins.
assembly of the ests resulted in a single cdna sequence.
two cdna clones from the original library were selected for sequencing by a commercial service , which both yielded identical sequences .
a cdna library form the horseshoe crab l. polyphemus  total rna was prepared and screened with specific anti-limulus-hemocyanin antibodies  <cit> .
four complete hemocyanin cdnas sequences were obtained by primer walking.
the cdnas were assigned after translation to distinct subunits on the basis of known n-termini  <cit> , identifying subunits ii, iiia, iv and vi .
an additional hemocyanin cdna  was identified from a set of ests  <cit>  and the complete coding sequence was obtained by primer walking .
a cloneminer  cdna library was constructed from total rna of a female emperor scorpion pandinus imperator  and submitted to  <dig> pyrosequencing  <cit> , resulting in  <dig>  high-quality reads.
hemocyanin subunit sequences were deduced from the assembled contigs .
a single whip spider euphrynichus bacillifer  was purchased from a commercial pet supplier.
a cdna library was constructed employing the mint universal kit .
<dig>  reads were obtained from the cdna by  <dig> pyrosequencing and the hemocyanin sequences were deduced from the assembled contigs .
a single whip scorpion mastigoproctus giganteus  was purchased from a commercial pet supplier.
total rna was extracted and converted into a mint universal cdna library, from which  <dig>  reads were obtained by  <dig> pyrosequencing.
full length coding sequences of hemocyanin subunits a, d, e, f, and g, as well as partial sequences from subunits b and c were deduced from the assembled contigs.
the missing fragments of subunits b and c were obtained by rt-rcr employing gene specific primers.
the cdna fragments were cloned into pgem and sequenced.
the final subunit sequences have been submitted to the databases under the accession numbers fr865920-fr <dig>
sequence assembly and analyses
the web-based tools provided by the expasy molecular biology server of the swiss institute of bioinformatics  were used for cdna translation and the analyses of amino acid sequences.
a multiple sequence alignment of the amino acid sequences of all available arthropod hemocyanins and selected arthropod phenoloxidases was constructed employing mafft  <dig>  <cit>  with the g-ins-i routine and the blosum  <dig> matrix at http://mafft.cbrc.jp/alignment/server/.
a complete list of sequences used in this study is provided in additional file  <dig>  for the phylogenetic inferences, signal peptides as well as the n- and c-terminal extensions of some phenoloxidases and hemocyanins were excluded from the multiple sequence alignment.
the final alignment comprised  <dig> sequences and  <dig> positions .
phylogenetic analyses
the most appropriate model of amino acid sequence evolution  was selected with prottest  <cit>  using the akaike information criterion.
bayesian phylogenetic analysis was performed using mrbayes  <dig> .
<dig>  <cit> .
we assumed the wag model with a gamma distribution of substitution rates.
metropolis-coupled markov chain monte carlo sampling was performed with one cold and three heated chains.
two independent runs were performed in parallel for  <dig>  million generations until the average standard deviation of split frequencies was <  <dig> .
starting trees were random and the trees were sampled every 100th generation.
the program tracer  <dig>   was used to examine log-likelihood plots and markov chain monte carlo summaries for all parameters.
posterior probabilities were estimated on the final  <dig>  trees .
trees were displayed using the arthropod phenoloxidases as an outgroup  <cit> .
molecular clock calculations
the program phylobayes  <dig>  was used for molecular clock estimates  <cit> , employing the mrbayes consensus tree as input.
first three relaxed clock models, the lognormal autocorrelated clock model   <cit> , the cox-ingersoll-ross process   <cit>  and uncorrelated gamma multipliers   <cit>  were compared by tenfold cross-validation with eight replicates, as specified in phylobayes  <dig> .
rates across sites were modeled assuming a discrete gamma distribution with four categories.
divergence time priors were either uniform or modeled with a birth death process.
node ages were calculated using either hard constrains, which do not allow calibrated nodes to fall outside the calibration dates or soft bounds, which allows for divergence times outside the calibration interval  <cit> .
for each of these settings rates across sites were modeled assuming a discrete gamma distribution with four categories and a dirichlet process.
all calculations were run for  <dig>   cycles.
the tree was calibrated with fossil constraints  <cit> .
the maximum age for the separation of arthropod subphyla and thus the maximum age of the origin of the chelicerata was set to the base of the cambrian period  <dig> ma ago .
stratigraphic information was obtained from http://www.fossilrecord.net <cit> .
numerical ages derive from the "international stratigraphic chart"  <dig>  .
RESULTS
hemocyanin sequences
a putative hemocyanin of the sea spider e. spinosa was identified in the ests of this species  <cit> .
the esphc <dig> cdna measures  <dig>  bp and translates into a protein of  <dig> amino acids with a predicted molecular mass of  <dig>  kda .
full length cdnas of five hemocyanin subunits of the atlantic horseshoe crab l. polyphemus were obtained .
the predicted l. polyphemus hemocyanin subunits measure between  <dig> and  <dig> amino acids, with molecular masses of  <dig> - <dig>  kda .
subunits i, iia and v could not identified in the cdna library  <cit>  or in the ests  <cit> .
seven hemocyanin subunit are available from the mangrove horseshoe crab c. rotundicauda, which corresponds to the subunits i, ii, iiia, iiib, iv, v, and vi .
thus seven of the eight hemocyanin subunit types of xiphosura could be included in our analyses.
the nature of an eighth type identified exclusively in l. polyphemus  and its exact topological position within the  <dig> × 6mer remains unclear.
its n-terminal sequence resembles that of subunit iiia and it may occupy a homologous position in some hexamers  <cit> .
<dig> pyrosequencing was employed to obtain ests from the scorpion p. imperator  <cit> , the whip spider e. bacillifer , the whip scorpion m. giganteus , the pseudoscorpion chelifer cancroides , the sun spider gluvia dorsalis  and the harvestman phalangium opilio .
hemocyanin sequences were identified in e. bacillifer and m. giganteus, but not in p. opilio , g. dorsalis  or c. cancroides .
in both, e. bacillifer and m. giganteus seven hemocyanin sequences were found, which are orthologous to e. californicum hemocyanin subunits a-g .
the predicted proteins measure between  <dig> and  <dig> amino acids, with molecular masses of  <dig> - <dig>  kda .
the publicly available chelicerate ests were searched for hemocyanin sequences using the tblastn algorithm .
we identified hemocyanin sequences in the ests from the tarantula acanthoscurria gomesiana , the hunting spider c. salei , the orb web spider nephila antipodiana  and the tarantula aphonopelma sp.
.
a total of  <dig> ests from the house spider parasteatoda tepidariorum display significant similarities to hemocyanin and represent all seven subunits  found in other araneae.
however, none of the ests from c. salei, n. antipodiana, aphonopelma sp.
or p. tepidariorum could be assembled to a complete coding sequence and therefore these sequences were not included in our analyses.
lorenzini and colleagues obtained  <dig>  ests from a hemocyte library of a. gomesiana  <cit> .
a total of  <dig> ests display significant similarities with hemocyanin, as identified by blast searches.
these sequences were assembled into eight contigs.
on the basis of similarity searches, seven subunits were assigned to arachnid hemocyanin subunits a-g .
the eighth putative hemocyanin sequence  has no obvious ortholog among the chelicerate hemocyanin subunits.
phylogeny of chelicerate hemocyanin subunits
bayesian phylogenetic analyses show that the hemocyanins tree basically follows the accepted arthropod relationships on the level of the subphyla .
the euarthropods split in the two sister groups chelicerata and mandibulata; though, the clade comprising the hemocyanins of the mandibulata  is poorly supported .
myriapod and pancrustacean hemocyanins each are monophyletic, whereas the crustacean hemocyanins are not because the remipede hemocyanins were found more closely related to those of the insects than to the malacostracan hemocyanins  <cit> .
chelicerate hemocyanins are monophyletic .
the lineage leading to the hemocyanin subunit from the sea spider e. spinosa  diverges first.
within the euchelicerate hemocyanins, four well-supported clades of distinct subunit types were identified .
clade  <dig> is the first branch within the euchelicerate hemocyanins and is formed by the arachnid b/c-type and xiphosuran v/vi-type subunits.
these subunits facilitate the inter-hexamer contacts within the  <dig> × 6mer, as well as between the  <dig> × 6mer half-structures  <cit> .
while the xiphosuran subunits v and vi form a common monophyletic clade, the arachnid hemocyanin subunits c do not.
here, the sequence of the putative c-subunit of the whip spider e. bacillifer  was found more diverged, mirrored by its basal position within the b/c clade.
clade  <dig> comprises the arachnid a-type and the xiphosuran type ii subunits.
in agreement with immunological studies with the scorpion a. australis  <cit> , the scorpion p. imperator has two a-type subunits, of which pimhc3a groups with a basal position to the other arachnid a-type subunits, while the position of pimhc3b is not well resolved.
note that subunit 3b represents a unique feature of scorpions that does not occur in the other arachnid hemocyanins.
the common clade that includes the remaining euchelicerate subunits received  <dig>  bayesian support.
this clade splits into two sub-clades, leading to arachnid subunits d and f, and xiphosuran subunits iiib and iv, on the one hand , and arachnid subunits e and g, and xiphosuran subunits i and iiia, on the other .
within the arachnid d/f and e/g subtrees, respectively, the observed subunit relationships essentially mirror the expected phylogeny of the species, although the positions of the scorpion subunits are somewhat ambiguous .
the hcx sequence, which had only been identified in a. gomesiana, is included in the d/f-clade, a position that is confirmed by pairwise comparisons, revealing that agohcx displays the highest sequence similarity to the arachnid f-subunits .
exclusion of agohcx from the phylogenetic analyses results in the same topology as described here, but bayesian support slightly increased throughout the tree .
entelegyne spiders of the rta-clade, such as c. salei, are known to diverge from the arachnid standard scheme of  <dig> × 6mer hemocyanins and to possess a mixture of  <dig> × 6mers and  <dig> × 6mers  <cit> .
in agreement with previous studies  <cit> , we found that the c. salei hemocyanin subunits all belong to the arachnid g-type, whereas the other six types  appear to have been lost during the evolution of this taxon.
the c. salei hemocyanin subunits consistently group with the subunit g of n. inaurata, an entelegyne spider with a  <dig> × 6mer hemocyanin.
molecular clock analyses of chelicerate hemocyanins
a timescale of chelicerate hemocyanin evolution was inferred on the basis of the chelicerate hemocyanins described above.
due to the divergent evolutionary rates , we excluded the hemocyanin sequences of onychophora, myriapoda, crustacea and hexapoda, as well as the phenoloxidases.
cross-validation shows that that the uncorrelated gamma clock model  fits better than the log-normal model or the cir process.
the cross-validation score of ugm vs. log-normal was  <dig>  +/-  <dig>  and of ugm vs. cir  <dig>  +/-  <dig> .
further support comes from a comparison of the calculated divergence times of different orthologous subunit pairs , which are most similar under ugm.
thus ugm was applied in our calculations.
the time estimates with soft bounds and hard bounds were similar, with estimates being on average ~5% older when using hard bounds .
no notable difference was observed between the results obtained with rates modeled with the dirichlet process or with a discrete gamma distribution .
the divergence times resulting from the estimates obtained with soft bounds, a birth death process and rates modeled with the gamma distribution are displayed in figure  <dig> .
we applied five calibration points, which correspond to four upper  and two lower  bounds derived from the fossil record .
subunits with uncertain orthology were ignored.
we first calculated the divergence times of the distinct subunit types in arachnida and xiphosura .
the earliest split of euchelicerate hemocyanins is formed by the clade of arachnid b/c and xiphosuran v/vi-type subunits, which occurred ~ <dig> ma.
within this clade, the exact arrangement of the three clades consisting of i. arachnid b-type subunits, ii.
arachnid c-type subunits and xiphosuran v/vi-type subunits is not well resolved , which is reflected by a rapid diversification 437- <dig> ma.
xiphosuran subunits v and vi separated ~ <dig> ma.
the a-type subunits, which include the xiphosuran subunit ii, split ~ <dig> ma, followed by the four clades defined above, consisting of arachnid subunits d/f, arachnid e/g, xiphosuran i/iiia, and xiphosuran iiib/iv subunits, which commenced to diversify ~ <dig> ma.
arachnid d- and f-type subunits split ~ <dig> ma, arachnid e- and g-type subunits ~ <dig> ma.
separation events within the xiphosuran subunits  occurred  <dig> and  <dig> ma, respectively.
in subsequent analyses, the hemocyanin sequences were employed to estimate chelicerate divergence times .
hemocyanin-derived phylogeny of chelicerata
we specifically studied the relationships among chelicerate taxa by concatenating the seven orthologous hemocyanin subunit sequences, representing spider subunits a-g, into a single sequence alignment.
clear orthologs were available for e. californicum, a. gomesiana, n. inaurata, e. bacillifer, and m. giganteus, respectively.
there is no c-type subunit in n. inaurata, which was coded as missing data.
p. imperator has two a-type subunits , of which we selected pimhciiia on the basis of its closer relationship to araneaen a-type subunits.
in case of c. rotundicauda, clear orthology assessment was only possible for subunit a .
we assigned crohcv and vi to b and c, crohciiib and iv to d and f, and crohci and iiia to e and g-subunits.
in a second approach, we exchanged these sequences.
because the single e. spinosa hemocyanin subunit constitutes a conclusive outgroup for all seven subunit types , seven copies of this sequence were concatenated.
bayesian phylogenetic analyses were applied and the sea spider e. spinosa was used as outgroup .
all nodes displayed a bayesian support of  <dig>  and there was no effect of an exchange of the ambiguous subunits from c. rotundicauda.
the xiphosura  are the sister taxon of the arachnida.
within the arachnida, the scorpion p. imperator diverged first.
whip spiders  and whip scorpions  form a common clade .
the pedipalpi form the sister taxon of the araneae, represented by the mygalomorph spiders e. californicum and a. gomesiana on the one hand, and the entelegyne spider n. inaurata on the other.
according to the molecular clock calculations , the hemocyanins of pycnogonida and euchelicerata diverged ~ <dig> ma .
considering the different orthologous subunits, xiphosuran and arachnid hemocyanins separated between  <dig> to  <dig> ma .
we calculated that the orthologous hemocyanin subunits of scorpiones and tetrapulmonata  split ~ <dig> ma .
within the scorpions, p. imperator  and a. australis  separated ~ <dig> ma.
the hemocyanins of pedipalpi and araneae diverged ~ <dig> ma , those of amblypygi and uropygi  <dig> ma .
within the araneae, hemocyanins of n. inaurata  and the mygalomorphae  diverged ~ <dig> ma .
the hemocyanins of e. californicum and a. gomesiana split  <dig> ma .
n. inaurata subunit g and the c. salei hemocyanins separated ~ <dig> ma; the xiphosurans l. polyphemus and c. rotundicauda diverged ~ <dig> ma .
discussion
hemocyanin subunits assemble into hexamers, which may form quaternary structures comprising up to  <dig> × 6mers  <cit> .
the evolutionary advantage of large oligomers presumably lies in a higher o2-carrying capacity per mol and higher cooperativity, which also enhances the o <dig> transport, combined with a low viscosity and a low colloid-osmotic pressure of the hemolymph.
the phylogenetic tree permits inferring the origins and modifications of these complex protein structures in the chelicerates.
the presence of only a single subunit in e. spinosa along with its basal position in the tree suggests that early chelicerate hemocyanins had a simple, homo-hexameric structure .
this hypothesis is supported by the independent emergence of hemocyanin oligo-hexamers in the other arthropod subphyla, which hints to more simple hemocyanins in the last common arthropod ancestor  <cit> .
early emergence of euchelicerate hemocyanin oligomers
the phylogenetic analyses demonstrate that the early euchelicerate hemocyanin, which was already used for o <dig> supply in the last common ancestor of the arachnids and the xiphosurans more than  <dig> ma, was composed of at least four distinct subunit types.
these subunits were the ancestors of the subunits represented in clades 1- <dig>  respectively .
according to our phylogenetic reconstruction, clade  <dig>  which includes arachnid b/c and xiphosuran v/vi subunits, diverged first.
both b/c and v/vi-subunits form heterodimers, which are responsible for the contacts between the hexamers  <cit> .
thus the emergence of clade 1-subunits was most likely associated with the organization of the first oligo-hexameric hemocyanin .
this event must have taken place very early in the evolution of euchelicerates, and may be associated with significant morphological and physiological changes.
we calculated that the corresponding gene duplication occurred ~ <dig> ma .
it may be speculated that first a  <dig> × 6mer consisting of two distinct subunit types evolved.
such dodecamers occur today in various crustaceans  <cit>  and in the spiders of the rta-clade  <cit> .
alternatively, this gene duplication has already resulted in a typical chelicerate  <dig> × 6mer, which is stabilized exclusively by a central tetrameric ring of clade  <dig> subunits.
notably, stable  <dig> × 6mers were obtained in hybrid reassembly experiments from mixtures of a clade  <dig> heterodimer as "linker" and another subunit type as "hexamer former"  <cit> .
these experiments convincingly worked with scorpion heterodimer 5b-3c plus l. polyphemus subunit ii, tarantula heterodimer b-c plus scorpion subunit  <dig>  and limulus heterodimer v-vi plus scorpion subunit  <dig>
the next step in evolution was the separation of clade  <dig> from the remaining subunits ~ <dig> ma.
clade  <dig> includes arachnid subunit a and xiphosuran subunit ii.
with the exception of pimhc3b , the a-subunits follow the expected phylogeny of the euchelicerates.
in the hemocyanin quaternary structure, subunit a/ii is located at the inter-hexamer interface of the basic  <dig> × 6mer  <cit> .
it may therefore be speculated that this step was required not only for stabilization of the  <dig> × 6mer hemocyanin, but also for improving cooperativity .
an additional gene duplication event, which may have taken place ~ <dig> ma, long before the arachnid-xiphosuran split, gave rise to both, clade  <dig>  consisting of arachnid subunits d and f, and xiphosuran subunits iiib and iv, and clade  <dig>  comprising arachnid subunits e and g, and xiphosuran subunits i and iiia.
even though this largely is in accordance with previous comparative immunochemical studies  <cit> , a common origin of subunits e and i, d and iv, f and iiib, and iiia and g, respectively, as proposed before  <cit> , can now be excluded.
the separation of the arachnid subunits d and f, and e and g, respectively, occurred before scorpions and spiders diverged.
thus the last common ancestor of all arachnids had a  <dig> × 6mer of seven subunits that were similar to subunit types a-g, demonstrating the evolutionary success of this conserved structure, which has remained essentially unchanged for more than  <dig> ma.
independent but parallel evolution of hemocyanin oligomers in xiphosura and arachnida
in horseshoe crabs, a large  <dig> × 6mer evolved.
although cooperativity is not further enhanced by this step, it might have been required for reducing the osmotic pressure and the viscosity of the hemolymph  <cit> .
notably, the duplication of the subunits corresponding to arachnid b and c , d and f  and e and g  occurred independently and the separation of xiphosuran subunits occurred more than  <dig> ma later.
nevertheless, a comparable subunit diversity evolved: six distinct subunit types for each topological position within the hexamer, plus a heterogeneity of the central linker unit to form an asymmetric  <dig> × 6mer.
consequently, today both arachnid and xiphosuran hemocyanins consist of seven distinct subunits , but only subunits a and ii  are one-to-one-orthologs.
thus, there was an evolutionary pressure to maximize the distinctiveness of subunits within the hemocyanin hexamer, which may be explained by better regulatory properties.
evolution of the arachnid hemocyanins
as outlined above, the principal structure of the early arachnid hemocyanin was most likely a  <dig> × 6mer .
previous immunological and structural investigations identified orthologs between scorpions and spiders   <cit> .
the phylogenetic tree  shows that these studies were essentially correct.
we should also note that in our previous phylogenetic analyses, the subunit aauhc <dig> was assigned to the araneaen g-type subunits  <cit> , while on the basis of the structural and immunological similarities aauhc <dig> was homologized with e-type subunits  <cit> .
the new tree, which includes more sequences, suggests that the protein-based studies were correct and aauhc <dig> is indeed an e-type subunit.
the  <dig> × 6mer structure is present in the mygalomorph spiders , and is found also in many entelegynae   <cit> .
the entelegyne spiders of the rta-clade, however, diverge from this standard scheme and have a mixture of  <dig> × 6mer and  <dig> × 6mer hemocyanins  <cit> .
this hemocyanin type is built by six distinct g-type subunits, with subunit csahc <dig>  forming the inter-hexamer bridge within the  <dig> × 6mer molecules  <cit> .
this suggests a loss of the other six subunit types  during evolution.
thus the ancestor of rta-clade spiders most likely had a simple hexameric hemocyanin, exclusively built by g-type subunits.
some  <dig> ma, the reconstruction of a more complex hemocyanin type commenced.
this might be explained by physiological and behavioral changes that e.g.
required a higher oxygen capacity in the hemolymph and/or a hemocyanin with a higher cooperativity.
there is little information about the subunit agohcx, which is only known from the ests of a. gomesiana  <cit> .
phylogenetic analyses place the protein with a long branch at the base of the arachnid d/f-subunit clade.
the fact that the agohcx sequence was found in ests and does not display any nonsense mutation suggests that this unique hemocyanin-like protein is translated into a functional protein.
however, it is neither known whether hcx is restricted to certain taxa  nor whether it is component of the hemocyanin oligomer.
the latter seems to be unlikely because of its derived sequence.
in addition, no evidence for hcx was found in the hemocyanin of the closely related tarantula e. californicum, despite more than  <dig> years of research.
absence of hemocyanin in some chelicerate taxa
notably, some arachnids do not have hemocyanin or any other o2-transport protein in their hemolymph.
despite the large number of ests obtained, no hemocyanin sequences were detected in the harvestman p. opilio , the pseudoscorpion c. cancroides  and the sun spider g. dorsalis .
this observation essentially agrees with previous findings  <cit> , with the exception that kempter et al.
<cit>  suggested the presence of a dodecameric hemocyanin in the harvestman leiobunum limbatum.
however, re-evaluation of the original data and new experiments suggest that the protein in question may actually be a vitellogenin-like, di-tetrameric protein , similar to those found in other arachnids  <cit> .
in this context, we would also like to note that in contrast to previous suggestions  <cit> , the haplogyne spider dysdera does not possess a  <dig> × 6mer hemocyanin.
a recent reinvestigation of a number of individuals demonstrated that dysdera lacks any hemocyanin, but express a tetrameric non-respiratory protein .
we formally cannot exclude that hemocyanins are present in other species of these taxa or are expressed only under certain environmental conditions.
however, we consider such scenario unlikely because such specific expression is not observed in other chelicerates.
in addition, no evidence for hemocyanin was found in the  <dig>  ests or the available genomic sequences of acari.
thus mites and ticks most likely lack hemocyanin as well.
opiliones, pseudoscorpiones and solifugae are apulmonate arachnids.
the absence of hemocyanin may be a synapomorphic character and an indication for a close relationship of these taxa, which agrees with some character-based phylogenetic studies  <cit>  .
morphological and/or physiological characteristics may have rendered a respiratory protein unnecessary.
notably, apulmonate arachnids do not breathe through book lungs, but possess trachea, which may be sufficient to support the aerobic metabolism.
the hemocyanin-less acari  are usually small and also have trachea  <cit> .
by contrast, spiders, scorpions, whip spiders, whip scorpions have book lungs, which are filled with hemolymph and which may limit o <dig> consumption.
here, hemocyanin may be required for efficient o <dig> uptake and distribution.
implications for chelicerate phylogeny
hemocyanin sequences have been successfully used to infer arthropod phylogeny  <cit> .
traditionally, chelicerata were considered as the sister group of the mandibulata, a taxon that comprises myriapoda, crustacea and hexapoda  <cit> .
however, several molecular phylogenetic studies have provided evidence for a common clade of myriapoda and chelicerata .
morphological evidence is poor and restricted to similarities of neurogenesis  <cit> , which may, however, also represent a plesiomorphic state.
in our study, we received some support for the monophyly of mandibulata, which agrees with previous studies employing hemocyanin sequences  <cit> , as well as other molecular approaches  <cit> .
the relationship among the major chelicerate lineages is controversial.
phylogenetic trees derived from the hemocyanin subunit sequences  or the concatenated alignment  can also be used to deduce the relative position of some chelicerate taxa, while others cannot be considered due to the lack of hemocyanin .
notably, weygoldt and paulus  <cit>  suggested a taxon "apulmonata", which joins solifugae, opiliones, pseudoscorpiones, acari, ricinulei , and palpigradi .
palpigradi and ricinulei were not available for our studies.
however, the absence of hemocyanin in acari, opiliones, pseudoscorpiones and solifugae tentatively supports monophyletic "apulmonata".
morphological and molecular studies have placed the pycnogonids  either as sister group of the euchelicerata  <cit> , nested within the chelicerata  <cit> , or considered them as the sister group of all other euarthropoda .
our phylogenetic tree  strongly supports the inclusion of the pycnogonida in the chelicerata as sister group of the euchelicerata.
this position is also tentatively supported by the hemocyanin mono-hexamer and is in line with recent neuroanatomical studies, which demonstrated the homology of deuterocerebral appendages of pycnogonida and euchelicerata  <cit> .
an ingroup position of the pycnogonida with the arachnida, as deduced from complete mitochondrial dna sequences  <cit>  is not supported by our data and should be considered unlikely.
in agreement with morphological considerations and most previous molecular phylogenetic studies, the xiphosura form the sister group of the arachnida.
within the arachnids, the relative positions of scorpiones, araneae, uropygi and amblypygi are controversial  <cit> .
we found monophyletic tetrapulmonata , which is the sister group of the scorpions.
in previous phylogenetic analyses, the relative position of the taxa araneae, amblypygi and uropygi has been controversial.
while some morphological studies favor a sister group relationship between araneae and amblypygi, forming the taxon labellata  <cit> , others support a common taxon referred to as pedipalpi, which comprises the uropygi and amblypygi  <cit> .
the latter view is supported by our molecular phylogenetic trees.
this finding holds for the tree derived from the concatenated alignment  as well as for most analyses of single subunits , with the exception of an unusual position of e. bacillifer subunit c and unresolved relationships among the subunits e.
the fossil record of chelicerates is far from being complete, but still allows the estimation of the evolutionary history of this taxon  <cit> .
the true origin of the chelicerates is currently uncertain, but dates back at least to the early cambrian period  <cit> .
we calculated that the first split within the chelicerates occurred  <dig> ma .
this slightly predates the earliest stemline chelicerates, which derive from the lower cambrian maotianshan shale some  <dig> ma  <cit> .
the first putative pycnogonid derives from the orsten fauna ~ <dig> ma  <cit> , the oldest xiphosuran fossil was found in a late ordovician lagerstätte and dates ~ <dig> ma  <cit>  and the first unambiguous arachnid is a ~ <dig> ma old silurian scorpion  <cit> .
these dates are actually close to our molecular clock estimates .
the first fossils of true spiders, whip spiders and whip scorpions were found in carboniferous strata, which are 310- <dig> ma old, although the tetrapulmonata are probably of devonian origin  <cit> .
our calculations confirm this notion and date the origin of the clade leading to tetrapulmonata  <dig> ma.
the oldest opisthothele fossil  is a mygalomorph spider dating  <dig> ma, while the oldest representative of the sistergroup araneomorphae  is of early cretaceous origin  <cit> .
we calculated the earliest divergence within the araneae  <dig> ma, which is somewhat older.
the most successful subgroup within the araneomorphae are the entelegynae, which are subdivided into the orbicularidae and the spiders of the rta clade.
the lower bound of divergence of the orbicularidae  and the rta-clade  is a net from an orbicularian spider from the early cretaceous period, some  <dig> ma.
we calculated that the formation of the cupiennius-type hemocyanin commenced about  <dig> ma, which should be considered as the lower bound for the time of emergence of the rta-clade.
CONCLUSIONS
our results clearly demonstrate that chelicerate hemocyanin structure is conservative, but also allows innovations.
there is little doubt that hemocyanin evolution commenced as a hexamer with a single subunit type, as present today in the sea spider.
the first hemocyanin oligo-hexamer emerged early in euchelicerate evolution, probably associated with the demand for better oxygen supply.
gene duplications led to the formation of a  <dig> × 6mer hemocyanin in early euchelicerates, which was structurally retained in the arachnids.
in xiphosurans, however, an  <dig> × 6mer hemocyanin built from two identical  <dig> × 6mers emerged.
although in both arachnids and xiphosurans at least two additional but independent subunit duplications occurred, the architecture of the  <dig> × 6mer has remained conserved in most taxa for more than  <dig> ma.
only in the spiders of the rta-clade, gene losses and independent duplications gave rise to a novel hemocyanin version, as exemplified by the  <dig> × 6mer hemocyanin of c. salei.
again changing physiological demands may have been the cause for these events.
the conservative structure of hemocyanins makes them an excellent marker to trace chelicerate evolution, which is only limited by the absence of hemocyanin in some taxa.
abbreviations
cir: cox-ingersoll-ross process; ests: expressed sequence tags; log: lognormal autocorrelated clock model; ma: million years ago; rta: retrolateral tibial apophysis; ugm: uncorrelated gamma multipliers.
authors' contributions
tb conceived the study and carried out the phylogenetic analyses.
pr, cp and jb provided and analyzed sequence data.
pr performed the molecular clock analyses.
jm drafted the model of hemocyanin evolution.
pr, jm and tb drafted the manuscript.
all authors read and approved the final version of the manuscript.
supplementary material
additional file 1
list of sequences used in this study.
the accession numbers of the cdna sequences are given, except , which has been derived by conventional protein sequencing.
su = subunit.
click here for file
additional file 2
multiple sequence alignment of chelicerate, crustacean, myriapod and insect hemocyanins, and selected arthropod phenoloxidases.
click here for file
additional file 3
molecular properties of chelicerate hemocyanin cdna and the deduced amino acid.
the asterisks  denote incomplete n-terminal sequences of p. imperator hemocyanins subunits, with  <dig> and  <dig> amino acids missing.
click here for file
additional file 4
phylogenetic tree of the chelicerate hemocyanin subunits excluding agohcx.
the numbers at the nodes represent bayesian posterior probabilities estimated with the wag model of amino acid substitution.
the species abbreviations are: aau, androctonus australis; ago, acanthoscurria gomesiana; cro, carcinoscorpius rotundicauda; csa, cupiennius salei; eba, euphrynichus bacillifer; eca, eurypelma californicum; esp, endeis spinosa; lpo, limulus polyphemus; mgi, mastigoproctus giganteus; nin, nephila inaurata; pin, pandinus imperator; ttr, tachypleus tridentatus.
see additional file  <dig> for abbreviations of the proteins.
click here for file
additional file 5
divergence times of chelicerate hemocyanin subunit types .
rates across sites were modeled assuming a gamma distribution  or with a dirichlet process .
divergence time priors were either uniform or modeled with a birth death process.
hard or soft bounds were applied.
divergence times are given in ma.
click here for file
additional file 6
divergence times of chelicerate taxa, as estimated from the hemocyanin sequences .
rates across sites were modeled assuming a gamma distribution  or with a dirichlet process .
divergence time priors were either uniform or modeled with a birth death process.
hard or soft bounds were applied.
divergence times are given in ma.
click here for file
acknowledgements
we thank falko roeding, carina edel and natalija winokurowa for their help with cloning and sequencing, and even solem for his support with the experiments with dysdera.
we are grateful to jason dunlop for checking the chelicerate calibration points.
the work was supported by grants of the deutsche forschungsgemeinschaft .

 hemocyanins are large copper-proteins that transport o in the hemolymph of many arthropods and mollusks. however, the hemocyanins of these two phyla are structurally different and emerged independently. hemocyanins evolved early in the arthropod stem lineage from the phenoloxidases, which are o2-consuming enzymes involved in the melanin pathway. other members of the arthropod hemocyanin superfamily have lost the ability to bind copper and thus o and gave rise to the non-respiratory pseudo-hemocyanins in decapod crustaceans and the hexamerins in hexapods, which serve as storage proteins. arthropod hemocyanins form hexamers or oligo-hexamers of identical or related subunits with a molecular mass of about kda. in each subunit, o2-binding is mediated by two cu+ ions, which are coordinated by six histidine residues. based on biochemical, immunochemical and molecular phylogenetic analyses, distinct hemocyanin subunit types have been identified in chelicerata, myriapoda, crustacea and hexapoda. these subunits experienced an independent evolution within each of these taxa, with the exception of a more complex pattern within the pancrustacea. within the chelicerates, biochemical analyses have demonstrated the presence of hemocyanins in xiphosura, scorpiones, uropygi, amblypygi, and araneae, but failed to identify these respiratory proteins in pycnogonida, solifugae and acari. complete sets of hemocyanin subunit sequences are available from the tarantula eurypelma californicum, the hunting spider cupiennius salei, the golden orb web spider nephila inaurata and the horse shoe crab carcinoscorpius rotundicauda. chelicerate hemocyanins are composed of up to eight distinct subunit types and form either × × × or × 6mers. each subunit type occupies a distinct position within the native oligomer. the subunits have similar oxygen binding properties, but different physico-chemical characteristics and evolutionary origins. xiphosurans have the largest hemocyanin molecules known, consisting of subunits and up to eight distinct subunits types in limulus polyphemus. scorpiones, amblypygi, uropygi, and some araneae have × 6mer hemocyanins. the × 6mer hemocyanin of the tarantula e. californicum is the best studied example and comprises seven distinct subunit types. a similar subunit composition was found in many other araneae, the amblypygi and the uropygi. among the araneae, variations from this "standard" scheme have been found in the entelegyne spiders of the rta-clade. in this large taxon, × 6mer and × 6mer hemocyanins occur that consist only of g-type subunits, but have lost the other six subunit types present in other araneae. scorpion hemocyanins are composed of eight subunit types named 3a, 3b, 3c, 5a and 5b. immunological and structural studies have suggested the correspondence between the distinct araneaen, scorpion and limulus hemocyanin subunits. this indicates that the evolution of distinct subunit types preceded the separation of xiphosura and arachnida. to understand the evolution of the complex oligomeric structures of chelicerate hemocyanins, we have obtained novel cdna sequences of hemocyanin subunits from a sea spider, a horseshoe crab, a whip scorpion and a whip spider. together with the previously sequenced hemocyanins from horseshoe crabs, scorpions and spiders, and those assembled from expressed sequence tags, full length chelicerate subunit sequences are available. these data allow us i. to trace the hemocyanin subunit evolution, ii. to reconstruct the emergence of the hemocyanin oligomers, and iii. to infer chelicerate phylogeny and divergence times. methods sequencing of chelicerate hemocyanin cdna full length coding sequences were available for the hemocyanins of the tarantula e. californicum, the mangrove horseshoe crab c. rotundicauda, the golden orb web spider n. inaurata and the hunting spider c. salei. in addition, hemocyanin primary structures had been obtained by conventional amino acid sequencing from the horseshoe crab tachypleus tridentatus and the scorpion androctonus australis. ests were generated from total rna of the sea spider endeis spinosa as described before. gene ontology assessment and blast searches identified six ests with significant similarities to arthropod hemocyanins. assembly of the ests resulted in a single cdna sequence. two cdna clones from the original library were selected for sequencing by a commercial service, which both yielded identical sequences. a cdna library form the horseshoe crab l. polyphemus total rna was prepared and screened with specific anti-limulus-hemocyanin antibodies. four complete hemocyanin cdnas sequences were obtained by primer walking. the cdnas were assigned after translation to distinct subunits on the basis of known n-termini, identifying subunits ii, iiia, iv and vi. an additional hemocyanin cdna was identified from a set of ests and the complete coding sequence was obtained by primer walking. a cloneminer cdna library was constructed from total rna of a female emperor scorpion pandinus imperator and submitted to pyrosequencing, resulting in high-quality reads. hemocyanin subunit sequences were deduced from the assembled contigs. a single whip spider euphrynichus bacillifer was purchased from a commercial pet supplier. a cdna library was constructed employing the mint universal kit. reads were obtained from the cdna by pyrosequencing and the hemocyanin sequences were deduced from the assembled contigs. a single whip scorpion mastigoproctus giganteus was purchased from a commercial pet supplier. total rna was extracted and converted into a mint universal cdna library, from which reads were obtained by pyrosequencing. full length coding sequences of hemocyanin subunits a, d, e, f, and g, as well as partial sequences from subunits b and c were deduced from the assembled contigs. the missing fragments of subunits b and c were obtained by rt-rcr employing gene specific primers. the cdna fragments were cloned into pgem and sequenced. the final subunit sequences have been submitted to the databases under the accession numbers fr865920-fr sequence assembly and analyses the web-based tools provided by the expasy molecular biology server of the swiss institute of bioinformatics were used for cdna translation and the analyses of amino acid sequences. a multiple sequence alignment of the amino acid sequences of all available arthropod hemocyanins and selected arthropod phenoloxidases was constructed employing mafft with the g-ins-i routine and the blosum matrix at http://mafft.cbrc.jp/alignment/server/. a complete list of sequences used in this study is provided in additional file for the phylogenetic inferences, signal peptides as well as the n- and c-terminal extensions of some phenoloxidases and hemocyanins were excluded from the multiple sequence alignment. the final alignment comprised sequences and positions. phylogenetic analyses the most appropriate model of amino acid sequence evolution was selected with prottest using the akaike information criterion. bayesian phylogenetic analysis was performed using mrbayes.. we assumed the wag model with a gamma distribution of substitution rates. metropolis-coupled markov chain monte carlo sampling was performed with one cold and three heated chains. two independent runs were performed in parallel for million generations until the average standard deviation of split frequencies was <. starting trees were random and the trees were sampled every 100th generation. the program tracer was used to examine log-likelihood plots and markov chain monte carlo summaries for all parameters. posterior probabilities were estimated on the final trees. trees were displayed using the arthropod phenoloxidases as an outgroup. molecular clock calculations the program phylobayes was used for molecular clock estimates, employing the mrbayes consensus tree as input. first three relaxed clock models, the lognormal autocorrelated clock model, the cox-ingersoll-ross process and uncorrelated gamma multipliers were compared by tenfold cross-validation with eight replicates, as specified in phylobayes. rates across sites were modeled assuming a discrete gamma distribution with four categories. divergence time priors were either uniform or modeled with a birth death process. node ages were calculated using either hard constrains, which do not allow calibrated nodes to fall outside the calibration dates or soft bounds, which allows for divergence times outside the calibration interval. for each of these settings rates across sites were modeled assuming a discrete gamma distribution with four categories and a dirichlet process. all calculations were run for cycles. the tree was calibrated with fossil constraints. the maximum age for the separation of arthropod subphyla and thus the maximum age of the origin of the chelicerata was set to the base of the cambrian period ma ago. stratigraphic information was obtained from http://www.fossilrecord.net. numerical ages derive from the "international stratigraphic chart".  hemocyanin sequences a putative hemocyanin of the sea spider e. spinosa was identified in the ests of this species. the esphc cdna measures bp and translates into a protein of amino acids with a predicted molecular mass of kda. full length cdnas of five hemocyanin subunits of the atlantic horseshoe crab l. polyphemus were obtained. the predicted l. polyphemus hemocyanin subunits measure between and amino acids, with molecular masses of - kda. subunits i, iia and v could not identified in the cdna library or in the ests. seven hemocyanin subunit are available from the mangrove horseshoe crab c. rotundicauda, which corresponds to the subunits i, ii, iiia, iiib, iv, v, and vi. thus seven of the eight hemocyanin subunit types of xiphosura could be included in our analyses. the nature of an eighth type identified exclusively in l. polyphemus and its exact topological position within the × 6mer remains unclear. its n-terminal sequence resembles that of subunit iiia and it may occupy a homologous position in some hexamers. pyrosequencing was employed to obtain ests from the scorpion p. imperator, the whip spider e. bacillifer, the whip scorpion m. giganteus, the pseudoscorpion chelifer cancroides, the sun spider gluvia dorsalis and the harvestman phalangium opilio. hemocyanin sequences were identified in e. bacillifer and m. giganteus, but not in p. opilio, g. dorsalis or c. cancroides. in both, e. bacillifer and m. giganteus seven hemocyanin sequences were found, which are orthologous to e. californicum hemocyanin subunits a-g. the predicted proteins measure between and amino acids, with molecular masses of - kda. the publicly available chelicerate ests were searched for hemocyanin sequences using the tblastn algorithm. we identified hemocyanin sequences in the ests from the tarantula acanthoscurria gomesiana, the hunting spider c. salei, the orb web spider nephila antipodiana and the tarantula aphonopelma sp.. a total of ests from the house spider parasteatoda tepidariorum display significant similarities to hemocyanin and represent all seven subunits found in other araneae. however, none of the ests from c. salei, n. antipodiana, aphonopelma sp. or p. tepidariorum could be assembled to a complete coding sequence and therefore these sequences were not included in our analyses. lorenzini and colleagues obtained ests from a hemocyte library of a. gomesiana. a total of ests display significant similarities with hemocyanin, as identified by blast searches. these sequences were assembled into eight contigs. on the basis of similarity searches, seven subunits were assigned to arachnid hemocyanin subunits a-g. the eighth putative hemocyanin sequence has no obvious ortholog among the chelicerate hemocyanin subunits. phylogeny of chelicerate hemocyanin subunits bayesian phylogenetic analyses show that the hemocyanins tree basically follows the accepted arthropod relationships on the level of the subphyla. the euarthropods split in the two sister groups chelicerata and mandibulata; though, the clade comprising the hemocyanins of the mandibulata is poorly supported. myriapod and pancrustacean hemocyanins each are monophyletic, whereas the crustacean hemocyanins are not because the remipede hemocyanins were found more closely related to those of the insects than to the malacostracan hemocyanins. chelicerate hemocyanins are monophyletic. the lineage leading to the hemocyanin subunit from the sea spider e. spinosa diverges first. within the euchelicerate hemocyanins, four well-supported clades of distinct subunit types were identified. clade is the first branch within the euchelicerate hemocyanins and is formed by the arachnid b/c-type and xiphosuran v/vi-type subunits. these subunits facilitate the inter-hexamer contacts within the × 6mer, as well as between the × 6mer half-structures. while the xiphosuran subunits v and vi form a common monophyletic clade, the arachnid hemocyanin subunits c do not. here, the sequence of the putative c-subunit of the whip spider e. bacillifer was found more diverged, mirrored by its basal position within the b/c clade. clade comprises the arachnid a-type and the xiphosuran type ii subunits. in agreement with immunological studies with the scorpion a. australis, the scorpion p. imperator has two a-type subunits, of which pimhc3a groups with a basal position to the other arachnid a-type subunits, while the position of pimhc3b is not well resolved. note that subunit 3b represents a unique feature of scorpions that does not occur in the other arachnid hemocyanins. the common clade that includes the remaining euchelicerate subunits received bayesian support. this clade splits into two sub-clades, leading to arachnid subunits d and f, and xiphosuran subunits iiib and iv, on the one hand, and arachnid subunits e and g, and xiphosuran subunits i and iiia, on the other. within the arachnid d/f and e/g subtrees, respectively, the observed subunit relationships essentially mirror the expected phylogeny of the species, although the positions of the scorpion subunits are somewhat ambiguous. the hcx sequence, which had only been identified in a. gomesiana, is included in the d/f-clade, a position that is confirmed by pairwise comparisons, revealing that agohcx displays the highest sequence similarity to the arachnid f-subunits. exclusion of agohcx from the phylogenetic analyses results in the same topology as described here, but bayesian support slightly increased throughout the tree. entelegyne spiders of the rta-clade, such as c. salei, are known to diverge from the arachnid standard scheme of × 6mer hemocyanins and to possess a mixture of × 6mers and × 6mers. in agreement with previous studies, we found that the c. salei hemocyanin subunits all belong to the arachnid g-type, whereas the other six types appear to have been lost during the evolution of this taxon. the c. salei hemocyanin subunits consistently group with the subunit g of n. inaurata, an entelegyne spider with a × 6mer hemocyanin. molecular clock analyses of chelicerate hemocyanins a timescale of chelicerate hemocyanin evolution was inferred on the basis of the chelicerate hemocyanins described above. due to the divergent evolutionary rates, we excluded the hemocyanin sequences of onychophora, myriapoda, crustacea and hexapoda, as well as the phenoloxidases. cross-validation shows that that the uncorrelated gamma clock model fits better than the log-normal model or the cir process. the cross-validation score of ugm vs. log-normal was +/- and of ugm vs. cir +/-. further support comes from a comparison of the calculated divergence times of different orthologous subunit pairs, which are most similar under ugm. thus ugm was applied in our calculations. the time estimates with soft bounds and hard bounds were similar, with estimates being on average ~5% older when using hard bounds. no notable difference was observed between the results obtained with rates modeled with the dirichlet process or with a discrete gamma distribution. the divergence times resulting from the estimates obtained with soft bounds, a birth death process and rates modeled with the gamma distribution are displayed in figure. we applied five calibration points, which correspond to four upper and two lower bounds derived from the fossil record. subunits with uncertain orthology were ignored. we first calculated the divergence times of the distinct subunit types in arachnida and xiphosura. the earliest split of euchelicerate hemocyanins is formed by the clade of arachnid b/c and xiphosuran v/vi-type subunits, which occurred ~ ma. within this clade, the exact arrangement of the three clades consisting of i. arachnid b-type subunits, ii. arachnid c-type subunits and xiphosuran v/vi-type subunits is not well resolved, which is reflected by a rapid diversification 437- ma. xiphosuran subunits v and vi separated ~ ma. the a-type subunits, which include the xiphosuran subunit ii, split ~ ma, followed by the four clades defined above, consisting of arachnid subunits d/f, arachnid e/g, xiphosuran i/iiia, and xiphosuran iiib/iv subunits, which commenced to diversify ~ ma. arachnid d- and f-type subunits split ~ ma, arachnid e- and g-type subunits ~ ma. separation events within the xiphosuran subunits occurred and ma, respectively. in subsequent analyses, the hemocyanin sequences were employed to estimate chelicerate divergence times. hemocyanin-derived phylogeny of chelicerata we specifically studied the relationships among chelicerate taxa by concatenating the seven orthologous hemocyanin subunit sequences, representing spider subunits a-g, into a single sequence alignment. clear orthologs were available for e. californicum, a. gomesiana, n. inaurata, e. bacillifer, and m. giganteus, respectively. there is no c-type subunit in n. inaurata, which was coded as missing data. p. imperator has two a-type subunits, of which we selected pimhciiia on the basis of its closer relationship to araneaen a-type subunits. in case of c. rotundicauda, clear orthology assessment was only possible for subunit a. we assigned crohcv and vi to b and c, crohciiib and iv to d and f, and crohci and iiia to e and g-subunits. in a second approach, we exchanged these sequences. because the single e. spinosa hemocyanin subunit constitutes a conclusive outgroup for all seven subunit types, seven copies of this sequence were concatenated. bayesian phylogenetic analyses were applied and the sea spider e. spinosa was used as outgroup. all nodes displayed a bayesian support of and there was no effect of an exchange of the ambiguous subunits from c. rotundicauda. the xiphosura are the sister taxon of the arachnida. within the arachnida, the scorpion p. imperator diverged first. whip spiders and whip scorpions form a common clade. the pedipalpi form the sister taxon of the araneae, represented by the mygalomorph spiders e. californicum and a. gomesiana on the one hand, and the entelegyne spider n. inaurata on the other. according to the molecular clock calculations, the hemocyanins of pycnogonida and euchelicerata diverged ~ ma. considering the different orthologous subunits, xiphosuran and arachnid hemocyanins separated between to ma. we calculated that the orthologous hemocyanin subunits of scorpiones and tetrapulmonata split ~ ma. within the scorpions, p. imperator and a. australis separated ~ ma. the hemocyanins of pedipalpi and araneae diverged ~ ma, those of amblypygi and uropygi ma. within the araneae, hemocyanins of n. inaurata and the mygalomorphae diverged ~ ma. the hemocyanins of e. californicum and a. gomesiana split ma. n. inaurata subunit g and the c. salei hemocyanins separated ~ ma; the xiphosurans l. polyphemus and c. rotundicauda diverged ~ ma. discussion hemocyanin subunits assemble into hexamers, which may form quaternary structures comprising up to × 6mers. the evolutionary advantage of large oligomers presumably lies in a higher o2-carrying capacity per mol and higher cooperativity, which also enhances the o transport, combined with a low viscosity and a low colloid-osmotic pressure of the hemolymph. the phylogenetic tree permits inferring the origins and modifications of these complex protein structures in the chelicerates. the presence of only a single subunit in e. spinosa along with its basal position in the tree suggests that early chelicerate hemocyanins had a simple, homo-hexameric structure. this hypothesis is supported by the independent emergence of hemocyanin oligo-hexamers in the other arthropod subphyla, which hints to more simple hemocyanins in the last common arthropod ancestor. early emergence of euchelicerate hemocyanin oligomers the phylogenetic analyses demonstrate that the early euchelicerate hemocyanin, which was already used for o supply in the last common ancestor of the arachnids and the xiphosurans more than ma, was composed of at least four distinct subunit types. these subunits were the ancestors of the subunits represented in clades 1- respectively. according to our phylogenetic reconstruction, clade which includes arachnid b/c and xiphosuran v/vi subunits, diverged first. both b/c and v/vi-subunits form heterodimers, which are responsible for the contacts between the hexamers. thus the emergence of clade 1-subunits was most likely associated with the organization of the first oligo-hexameric hemocyanin. this event must have taken place very early in the evolution of euchelicerates, and may be associated with significant morphological and physiological changes. we calculated that the corresponding gene duplication occurred ~ ma. it may be speculated that first a × 6mer consisting of two distinct subunit types evolved. such dodecamers occur today in various crustaceans and in the spiders of the rta-clade. alternatively, this gene duplication has already resulted in a typical chelicerate × 6mer, which is stabilized exclusively by a central tetrameric ring of clade subunits. notably, stable × 6mers were obtained in hybrid reassembly experiments from mixtures of a clade heterodimer as "linker" and another subunit type as "hexamer former". these experiments convincingly worked with scorpion heterodimer 5b-3c plus l. polyphemus subunit ii, tarantula heterodimer b-c plus scorpion subunit and limulus heterodimer v-vi plus scorpion subunit the next step in evolution was the separation of clade from the remaining subunits ~ ma. clade includes arachnid subunit a and xiphosuran subunit ii. with the exception of pimhc3b, the a-subunits follow the expected phylogeny of the euchelicerates. in the hemocyanin quaternary structure, subunit a/ii is located at the inter-hexamer interface of the basic × 6mer. it may therefore be speculated that this step was required not only for stabilization of the × 6mer hemocyanin, but also for improving cooperativity. an additional gene duplication event, which may have taken place ~ ma, long before the arachnid-xiphosuran split, gave rise to both, clade consisting of arachnid subunits d and f, and xiphosuran subunits iiib and iv, and clade comprising arachnid subunits e and g, and xiphosuran subunits i and iiia. even though this largely is in accordance with previous comparative immunochemical studies, a common origin of subunits e and i, d and iv, f and iiib, and iiia and g, respectively, as proposed before, can now be excluded. the separation of the arachnid subunits d and f, and e and g, respectively, occurred before scorpions and spiders diverged. thus the last common ancestor of all arachnids had a × 6mer of seven subunits that were similar to subunit types a-g, demonstrating the evolutionary success of this conserved structure, which has remained essentially unchanged for more than ma. independent but parallel evolution of hemocyanin oligomers in xiphosura and arachnida in horseshoe crabs, a large × 6mer evolved. although cooperativity is not further enhanced by this step, it might have been required for reducing the osmotic pressure and the viscosity of the hemolymph. notably, the duplication of the subunits corresponding to arachnid b and c, d and f and e and g occurred independently and the separation of xiphosuran subunits occurred more than ma later. nevertheless, a comparable subunit diversity evolved: six distinct subunit types for each topological position within the hexamer, plus a heterogeneity of the central linker unit to form an asymmetric × 6mer. consequently, today both arachnid and xiphosuran hemocyanins consist of seven distinct subunits, but only subunits a and ii are one-to-one-orthologs. thus, there was an evolutionary pressure to maximize the distinctiveness of subunits within the hemocyanin hexamer, which may be explained by better regulatory properties. evolution of the arachnid hemocyanins as outlined above, the principal structure of the early arachnid hemocyanin was most likely a × 6mer. previous immunological and structural investigations identified orthologs between scorpions and spiders. the phylogenetic tree shows that these studies were essentially correct. we should also note that in our previous phylogenetic analyses, the subunit aauhc was assigned to the araneaen g-type subunits, while on the basis of the structural and immunological similarities aauhc was homologized with e-type subunits. the new tree, which includes more sequences, suggests that the protein-based studies were correct and aauhc is indeed an e-type subunit. the × 6mer structure is present in the mygalomorph spiders, and is found also in many entelegynae. the entelegyne spiders of the rta-clade, however, diverge from this standard scheme and have a mixture of × 6mer and × 6mer hemocyanins. this hemocyanin type is built by six distinct g-type subunits, with subunit csahc forming the inter-hexamer bridge within the × 6mer molecules. this suggests a loss of the other six subunit types during evolution. thus the ancestor of rta-clade spiders most likely had a simple hexameric hemocyanin, exclusively built by g-type subunits. some ma, the reconstruction of a more complex hemocyanin type commenced. this might be explained by physiological and behavioral changes that e.g. required a higher oxygen capacity in the hemolymph and/or a hemocyanin with a higher cooperativity. there is little information about the subunit agohcx, which is only known from the ests of a. gomesiana. phylogenetic analyses place the protein with a long branch at the base of the arachnid d/f-subunit clade. the fact that the agohcx sequence was found in ests and does not display any nonsense mutation suggests that this unique hemocyanin-like protein is translated into a functional protein. however, it is neither known whether hcx is restricted to certain taxa nor whether it is component of the hemocyanin oligomer. the latter seems to be unlikely because of its derived sequence. in addition, no evidence for hcx was found in the hemocyanin of the closely related tarantula e. californicum, despite more than years of research. absence of hemocyanin in some chelicerate taxa notably, some arachnids do not have hemocyanin or any other o2-transport protein in their hemolymph. despite the large number of ests obtained, no hemocyanin sequences were detected in the harvestman p. opilio, the pseudoscorpion c. cancroides and the sun spider g. dorsalis. this observation essentially agrees with previous findings, with the exception that kempter et al. suggested the presence of a dodecameric hemocyanin in the harvestman leiobunum limbatum. however, re-evaluation of the original data and new experiments suggest that the protein in question may actually be a vitellogenin-like, di-tetrameric protein, similar to those found in other arachnids. in this context, we would also like to note that in contrast to previous suggestions, the haplogyne spider dysdera does not possess a × 6mer hemocyanin. a recent reinvestigation of a number of individuals demonstrated that dysdera lacks any hemocyanin, but express a tetrameric non-respiratory protein. we formally cannot exclude that hemocyanins are present in other species of these taxa or are expressed only under certain environmental conditions. however, we consider such scenario unlikely because such specific expression is not observed in other chelicerates. in addition, no evidence for hemocyanin was found in the ests or the available genomic sequences of acari. thus mites and ticks most likely lack hemocyanin as well. opiliones, pseudoscorpiones and solifugae are apulmonate arachnids. the absence of hemocyanin may be a synapomorphic character and an indication for a close relationship of these taxa, which agrees with some character-based phylogenetic studies. morphological and/or physiological characteristics may have rendered a respiratory protein unnecessary. notably, apulmonate arachnids do not breathe through book lungs, but possess trachea, which may be sufficient to support the aerobic metabolism. the hemocyanin-less acari are usually small and also have trachea. by contrast, spiders, scorpions, whip spiders, whip scorpions have book lungs, which are filled with hemolymph and which may limit o consumption. here, hemocyanin may be required for efficient o uptake and distribution. implications for chelicerate phylogeny hemocyanin sequences have been successfully used to infer arthropod phylogeny. traditionally, chelicerata were considered as the sister group of the mandibulata, a taxon that comprises myriapoda, crustacea and hexapoda. however, several molecular phylogenetic studies have provided evidence for a common clade of myriapoda and chelicerata. morphological evidence is poor and restricted to similarities of neurogenesis, which may, however, also represent a plesiomorphic state. in our study, we received some support for the monophyly of mandibulata, which agrees with previous studies employing hemocyanin sequences, as well as other molecular approaches. the relationship among the major chelicerate lineages is controversial. phylogenetic trees derived from the hemocyanin subunit sequences or the concatenated alignment can also be used to deduce the relative position of some chelicerate taxa, while others cannot be considered due to the lack of hemocyanin. notably, weygoldt and paulus suggested a taxon "apulmonata", which joins solifugae, opiliones, pseudoscorpiones, acari, ricinulei, and palpigradi. palpigradi and ricinulei were not available for our studies. however, the absence of hemocyanin in acari, opiliones, pseudoscorpiones and solifugae tentatively supports monophyletic "apulmonata". morphological and molecular studies have placed the pycnogonids either as sister group of the euchelicerata, nested within the chelicerata, or considered them as the sister group of all other euarthropoda. our phylogenetic tree strongly supports the inclusion of the pycnogonida in the chelicerata as sister group of the euchelicerata. this position is also tentatively supported by the hemocyanin mono-hexamer and is in line with recent neuroanatomical studies, which demonstrated the homology of deuterocerebral appendages of pycnogonida and euchelicerata. an ingroup position of the pycnogonida with the arachnida, as deduced from complete mitochondrial dna sequences is not supported by our data and should be considered unlikely. in agreement with morphological considerations and most previous molecular phylogenetic studies, the xiphosura form the sister group of the arachnida. within the arachnids, the relative positions of scorpiones, araneae, uropygi and amblypygi are controversial. we found monophyletic tetrapulmonata, which is the sister group of the scorpions. in previous phylogenetic analyses, the relative position of the taxa araneae, amblypygi and uropygi has been controversial. while some morphological studies favor a sister group relationship between araneae and amblypygi, forming the taxon labellata, others support a common taxon referred to as pedipalpi, which comprises the uropygi and amblypygi. the latter view is supported by our molecular phylogenetic trees. this finding holds for the tree derived from the concatenated alignment as well as for most analyses of single subunits, with the exception of an unusual position of e. bacillifer subunit c and unresolved relationships among the subunits e. the fossil record of chelicerates is far from being complete, but still allows the estimation of the evolutionary history of this taxon. the true origin of the chelicerates is currently uncertain, but dates back at least to the early cambrian period. we calculated that the first split within the chelicerates occurred ma. this slightly predates the earliest stemline chelicerates, which derive from the lower cambrian maotianshan shale some ma. the first putative pycnogonid derives from the orsten fauna ~ ma, the oldest xiphosuran fossil was found in a late ordovician lagerstätte and dates ~ ma and the first unambiguous arachnid is a ~ ma old silurian scorpion. these dates are actually close to our molecular clock estimates. the first fossils of true spiders, whip spiders and whip scorpions were found in carboniferous strata, which are 310- ma old, although the tetrapulmonata are probably of devonian origin. our calculations confirm this notion and date the origin of the clade leading to tetrapulmonata ma. the oldest opisthothele fossil is a mygalomorph spider dating ma, while the oldest representative of the sistergroup araneomorphae is of early cretaceous origin. we calculated the earliest divergence within the araneae ma, which is somewhat older. the most successful subgroup within the araneomorphae are the entelegynae, which are subdivided into the orbicularidae and the spiders of the rta clade. the lower bound of divergence of the orbicularidae and the rta-clade is a net from an orbicularian spider from the early cretaceous period, some ma. we calculated that the formation of the cupiennius-type hemocyanin commenced about ma, which should be considered as the lower bound for the time of emergence of the rta-clade.  our results clearly demonstrate that chelicerate hemocyanin structure is conservative, but also allows innovations. there is little doubt that hemocyanin evolution commenced as a hexamer with a single subunit type, as present today in the sea spider. the first hemocyanin oligo-hexamer emerged early in euchelicerate evolution, probably associated with the demand for better oxygen supply. gene duplications led to the formation of a × 6mer hemocyanin in early euchelicerates, which was structurally retained in the arachnids. in xiphosurans, however, an × 6mer hemocyanin built from two identical × 6mers emerged. although in both arachnids and xiphosurans at least two additional but independent subunit duplications occurred, the architecture of the × 6mer has remained conserved in most taxa for more than ma. only in the spiders of the rta-clade, gene losses and independent duplications gave rise to a novel hemocyanin version, as exemplified by the × 6mer hemocyanin of c. salei. again changing physiological demands may have been the cause for these events. the conservative structure of hemocyanins makes them an excellent marker to trace chelicerate evolution, which is only limited by the absence of hemocyanin in some taxa. abbreviations cir: cox-ingersoll-ross process; ests: expressed sequence tags; log: lognormal autocorrelated clock model; ma: million years ago; rta: retrolateral tibial apophysis; ugm: uncorrelated gamma multipliers. authors' contributions tb conceived the study and carried out the phylogenetic analyses. pr, cp and jb provided and analyzed sequence data. pr performed the molecular clock analyses. jm drafted the model of hemocyanin evolution. pr, jm and tb drafted the manuscript. all authors read and approved the final version of the manuscript. supplementary material additional file 1 list of sequences used in this study. the accession numbers of the cdna sequences are given, except, which has been derived by conventional protein sequencing. su = subunit. click here for file additional file 2 multiple sequence alignment of chelicerate, crustacean, myriapod and insect hemocyanins, and selected arthropod phenoloxidases. click here for file additional file 3 molecular properties of chelicerate hemocyanin cdna and the deduced amino acid. the asterisks denote incomplete n-terminal sequences of p. imperator hemocyanins subunits, with and amino acids missing. click here for file additional file 4 phylogenetic tree of the chelicerate hemocyanin subunits excluding agohcx. the numbers at the nodes represent bayesian posterior probabilities estimated with the wag model of amino acid substitution. the species abbreviations are: aau, androctonus australis; ago, acanthoscurria gomesiana; cro, carcinoscorpius rotundicauda; csa, cupiennius salei; eba, euphrynichus bacillifer; eca, eurypelma californicum; esp, endeis spinosa; lpo, limulus polyphemus; mgi, mastigoproctus giganteus; nin, nephila inaurata; pin, pandinus imperator; ttr, tachypleus tridentatus. see additional file for abbreviations of the proteins. click here for file additional file 5 divergence times of chelicerate hemocyanin subunit types. rates across sites were modeled assuming a gamma distribution or with a dirichlet process. divergence time priors were either uniform or modeled with a birth death process. hard or soft bounds were applied. divergence times are given in ma. click here for file additional file 6 divergence times of chelicerate taxa, as estimated from the hemocyanin sequences. rates across sites were modeled assuming a gamma distribution or with a dirichlet process. divergence time priors were either uniform or modeled with a birth death process. hard or soft bounds were applied. divergence times are given in ma. click here for file acknowledgements we thank falko roeding, carina edel and natalija winokurowa for their help with cloning and sequencing, and even solem for his support with the experiments with dysdera. we are grateful to jason dunlop for checking the chelicerate calibration points. the work was supported by grants of the deutsche forschungsgemeinschaft.
hemocyanin subunit sequences had been available from horseshoe crabs and various spiders, but not from any other chelicerate taxon. subsequently, independent gene duplication events gave rise to the other distinct subunits in each of the × 6mer hemocyanin of xiphosura and the × 6mer of arachnida. the hemocyanin sequences were used to infer the evolutionary history of chelicerates. arthropod hemocyanins are hexamers or oligomers of hexamers, which are characterized by a high o transport capacity and a high cooperativity, thereby enhancing o supply. 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. 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. hemocyanin subunits from a sea spider, a scorpion, a whip scorpion and a whip spider were sequenced. while the sea spider has a simple hexameric hemocyanin, four distinct subunit types evolved before xiphosura and arachnida diverged around ma ago, suggesting the existence of a × 6mer at that time. the molecular clock calculations were in excellent agreement with the fossil record, also demonstrating the applicability of hemocyanins for such approach. 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. oxygen transport in the hemolymph of many arthropod species is facilitated by large copper-proteins referred to as hemocyanins. a × 6mer hemocyanin consisting of seven subunit types is conserved in most arachnids since more than ma, although some entelegyne spiders display selective subunit loss and independent oligomerization.
 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 o transport capacity and a high cooperativity, thereby enhancing o supply. hemocyanin subunit sequences had been available from horseshoe crabs and various spiders, 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.  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. 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 ma ago, suggesting the existence of a × 6mer at that time. subsequently, independent gene duplication events gave rise to the other distinct subunits in each of the × 6mer hemocyanin of xiphosura and the × 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.  formation of a complex hemocyanin oligomer commenced early in the evolution of euchelicerates. a × 6mer hemocyanin consisting of seven subunit types is conserved in most arachnids since more than 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. 