The transcriptome analysis of the whole-body of the gastropod mollusk Limax flavus and screening of putative antimicrobial peptide and protein genes.

The gastropod mollusk Limax flavus, one of the most widespread pests in China, is used to treat infectious diseases in traditional Chinese medicine. However, little genomic information is available for this non-model species. In this study, the whole-body transcriptome of L. flavus was sequenced using next generation sequencing technology. A total of 6.81 Gb clean reads were obtained, which were assembled into 150,766 transcripts with 132,206 annotated unigenes. Functionally classification assigned 30,542 unigenes to 56 Gene Ontology terms, 16,745 unigenes were divided into 26 euKaryotic Ortholog Groups of proteins categories, and 13,854 unigenes were assigned to 230 Kyoto Encyclopedia of Genes and Genomes pathways. Furthermore, we identified 17,251 simple sequence repeats and several kinds of antimicrobial peptide and protein (AMPs) genes. The transcriptome data of L. flavus will provide a valuable genomic resource for further studies on this species, and the AMPs identified in L. flavus will support its medical potential.

Processive Recoding and Metazoan Evolution of Selenoprotein P: Up to 132 UGAs in Molluscs.

Selenoproteins typically contain a single selenocysteine, the 21st amino acid, encoded by a context-redefined UGA. However, human selenoprotein P (SelenoP) has a redox-functioning selenocysteine in its N-terminal domain and nine selenium transporter-functioning selenocysteines in its C-terminal domain. Here we show that diverse SelenoP genes are present across metazoa with highly variable numbers of Sec-UGAs, ranging from a single UGA in certain insects, to 9 in common spider, and up to 132 in bivalve molluscs. SelenoP genes were shaped by a dynamic evolutionary process linked to selenium usage. Gene evolution featured modular expansions of an ancestral multi-Sec domain, which led to particularly Sec-rich SelenoP proteins in many aquatic organisms. We focused on molluscs, and chose Pacific oyster Magallana gigas as experimental model. We show that oyster SelenoP mRNA with 46 UGAs is translated full-length in vivo. Ribosome profiling indicates that selenocysteine specification occurs with ∼5% efficiency at UGA1 and approaches 100% efficiency at distal 3' UGAs. We report genetic elements relevant to its expression, including a leader open reading frame and an RNA structure overlapping the initiation codon that modulates ribosome progression in a selenium-dependent manner. Unlike their mammalian counterparts, the two SECIS elements in oyster SelenoP (3'UTR recoding elements) do not show functional differentiation in vitro. Oysters can increase their tissue selenium level up to 50-fold upon supplementation, which also results in extensive changes in selenoprotein expression.

Contrasting patterns of population connectivity between regions in a commercially important mollusc Haliotis rubra: integrating population genetics, genomics and marine LiDAR data.

Estimating contemporary genetic structure and population connectivity in marine species is challenging, often compromised by genetic markers that lack adequate sensitivity, and unstructured sampling regimes. We show how these limitations can be overcome via the integration of modern genotyping methods and sampling designs guided by LiDAR and SONAR data sets. Here we explore patterns of gene flow and local genetic structure in a commercially harvested abalone species (Haliotis rubra) from southeastern Australia, where the viability of fishing stocks is believed to be dictated by recruitment from local sources. Using a panel of microsatellite and genomewide SNP markers, we compare allele frequencies across a replicated hierarchical sampling area guided by bathymetric LiDAR imagery. Results indicate high levels of gene flow and no significant genetic structure within or between benthic reef habitats across 1400 km of coastline. These findings differ to those reported for other regions of the fishery indicating that larval supply is likely to be spatially variable, with implications for management and long-term recovery from stock depletion. The study highlights the utility of suitably designed genetic markers and spatially informed sampling strategies for gaining insights into recruitment patterns in benthic marine species, assisting in conservation planning and sustainable management of fisheries.

The molecular evolution of the pif family proteins in various species of mollusks.

Various novel proteins have been identified from many kinds of mollusk shells. Although such matrix proteins are believed to play important roles in the calcium carbonate crystal formation of shells, no common proteins that interact with calcium carbonate or that are involved in the molecular mechanisms behind shell formation have been identified. Pif consists of two proteins, Pif 80 and Pif 97, which are encoded by a single mRNA. Pif 80 was identified as a key acidic protein that regulates the formation of the nacreous layer in Pinctada fucata, while Pif 97 has von Willebrand factor type A (VWA) and chitin-binding domains. In this study, we identified Pif homologues from Pinctada margaritifera, Pinctada maxima, Pteria penguin, Mytilus galloprovincialis, and in the genome database of Lottia gigantea in order to compare their primary protein sequences. The VWA and chitin-binding domains are conserved in all Pif 97 homologues, whereas the amino acid sequences of the Pif 80 regions differ markedly among the species. Sequence alignment revealed the presence of a novel significantly conserved sequence between the chitin-binding domain and the C-terminus of Pif 97. Further examination of the Pif 80 regions suggested that they share a sequence that is similar to the laminin G domain. These results indicate that all Pif molecules in bivalves and gastropods may be derived from a common ancestral gene. These comparisons may shed light on the correlation between molecular evolution and morphology in mollusk shell microstructure.

First molluscan theta-class Glutathione S-Transferase: identification, cloning, characterization and transcriptional analysis post immune challenges.

Glutathione S-Transferases (GSTs) are multifunctional cytosolic isoenzymes, distinctly known as phase II detoxification enzymes. GSTs play a significant role in cellular defense against toxicity and have been identified in nearly all organisms studied to date, from bacteria to mammals. In this study, we have identified a full-length cDNA of the theta class GST from Ruditapes philippinarum (RpGSTθ), an important commercial edible molluscan species. RpGSTθ was cloned and the recombinant protein expressed, in order to study its biochemical characteristics and determine its physiological activities. The cDNA comprised an ORF of 693 bp, encoding 231 amino acids with a predicted molecular mass of 27 kDa and an isoelectric point of 8.2. Sequence analysis revealed that RpGSTθ possessed characteristic conserved domains of the GST_N family, Class Theta subfamily (PSSM: cd03050) and GST_C_family Super family (PSSM: cl02776). Phylogenetic analysis showed that RpGSTθ evolutionarily linked with other theta class homologues. The recombinant protein was expressed in Escherichia coli BL21(DE3) cells and the purified enzyme showed high activity with GST substrates like CDNB and 4-NBC. Glutathione dependent peroxidase activity of GST, investigated with cumene hydroperoxide as substrate affirmed the antioxidant property of rRpGSTθ. By quantitative PCR, RpGSTθ was found to be ubiquitously expressed in all tissues examined, with the highest levels occurring in gills, mantle, and hemocytes. Since GSTs may act as detoxification enzymes to mediate immune defense, the effects of pathogen associated molecular pattern, lipopolysaccharide and intact Vibrio tapetis bacteria challenge on RpGSTθ gene transcription were studied. Furthermore, the RpGSTθ expression changes induced by immune challenges were similar to those of the antioxidant defense enzyme manganese superoxide dismutase (RpMnSOD). To our knowledge, RpGSTθ is the first molluscan theta class GST reported, and its immune-related role in Manila clam may provide insights into potential therapeutic targets for protecting this important aquaculture species.

Effects of dietary fish oil replacement with flaxseed oil on tissue fatty acid composition and expression of desaturase and elongase genes.

BACKGROUND: This study was conducted to evaluate the effects of total or partial substitution of dietary fish oil (FO) by flaxseed oil (FlaxO) in Jade Tiger hybrid abalone on fatty acid composition of muscle, gonad and digestive gland, and the expression of desaturase and elongase genes. Abalone were fed five different experimental diets in which FO (control diet) was serially replaced by 25%, 50%, 75% and 100% FlaxO respectively. RESULTS: Muscle, gonad and digestive gland of abalone fed the control diet and the diets containing 25%, 50% and 75% FlaxO showed significantly higher (P < 0.05) levels of eicosapentaenoic acid (EPA), docosapentaenoic acid (DPA) and docosahexaenoic acid (DHA) compared to those fed the 100% FlaxO. The results also showed that Δ-6 desaturase and elongase gene expression in muscle was increased in a graded manner by increasing dietary FlaxO. The expression of both genes was higher in abalone fed the FlaxO-substituted diets compared to the abalone fed FO. CONCLUSION: The replacement of FO with FlaxO in commercial abalone diets at levels of 25-75% can improve the composition of health-benefiting n-3 polyunsturated fatty acids in tissues of cultured hybrid abalone, and achieve similar outcomes to FO supplementation.

Gene cloning and biochemical characterization of the BMP-2 of Pinctada fucata.

The bone morphogenetic proteins (BMPs) constitute a subfamily of the transforming growth factor type beta (TGF-beta) supergene family. BMP-2 plays an important role not only in osteoblast differentiation but also in pattern formation during development. To determine the function of BMP-2 in Pinctada fucata development and hard tissue formation, we isolated a BMP-2 genomic DNA clone and the BMP-2 cDNA. The deduced BMP-2 sequence consisted of 447 amino acids. The BMP-2 gene was composed of three exons. The C-terminal portion (149 amino acids) had 86% and 66% identity to the Crassostrea gigas and the human BMP-2 sequence respectively. The 5' flanking promoter region contained putative glioma transcription factor (Gli) and retinoic acid receptor (RAR) responding elements. The BMP-2 gene was expressed strongly in the inner part of the mantle tissue, corresponding to the nacreous aragonite shell layer. This finding suggests that BMP-2 has a key role in nacreous layer formation.

A novel oligoalginate lyase from abalone, Haliotis discus hannai, that releases disaccharide from alginate polymer in an exolytic manner.

We previously reported the isolation and cDNA cloning of an endolytic alginate lyase, HdAly, from abalone Haliotis discus hannai [Carbohydr. Res.2003, 338, 2841-2852]. Although HdAly preferentially degraded mannuronate-rich substrates, it was incapable of degrading unsaturated oligomannuronates smaller than tetrasaccharide. In the present study, we used conventional chromatographic techniques to isolate a novel unsaturated-trisaccharide-degrading enzyme, named HdAlex, from the digestive fluid of the abalone. The HdAlex showed a molecular weight of 32,000 on SDS-PAGE and could degrade not only unsaturated trisaccharide but also alginate and mannuronate-rich polymers at an optimal pH and temperature of 7.1 and 42 degrees C, respectively. Upon digestion of alginate polymer, HdAlex decreased the viscosity of the alginate at a slower rate than did HdAly, producing only unsaturated disaccharide without any intermediate oligosaccharides. These results indicate that HdAlex degrades the alginate polymer in an exolytic manner. Because HdAlex split saturated trisaccharide producing unsaturated disaccharide, we considered that this enzyme cleaved the alginate at the second glycoside linkage from the reducing terminus. The primary structure of HdAlex was deduced with cDNAs amplified from an abalone hepatopancreas cDNA library by the polymerase chain reaction. The translational region of 822 bp in the total 887-bp sequence of HdAlex cDNA encoded an amino-acid sequence of 273 residues. The N-terminal sequence of 16 residues, excluding the initiation methionine, was regarded as the signal peptide of this enzyme. The amino-acid sequence of the remaining 256 residues shared 62-67% identities with those of the polysaccharide lyase family-14 (PL14) enzymes such as HdAly and turban-shell alginate lyase SP2. To our knowledge, HdAlex is the first exolytic oligoalginate lyase belonging to PL14.

Genes expressed in a turrid venom duct: divergence and similarity to conotoxins.

The toxoglossate mollusks are a large group of venomous animals (>10,000 species) conventionally divided into three groups, the cone snails, the auger snails, and the turrid snails; turrids account for >90% of the biodiversity of toxoglossans. Only the venoms of cone snails have been intensively investigated, with little work focused on turrids. We report the first broad characterization of genes expressed in venom ducts of any turrid species. Twenty-three different cDNA clones encoding putative toxins were characterized from the venom duct of the turrine species Lophiotoma olangoensis Olivera 2002 and belong to 16 different gene families. Of the 16 different Lophiotoma olangoensis gene families that encode putative toxins, for only 1 was there clear evidence of sequence similarity with any conotoxin gene family. The I-like gene family of Lophiotoma olangoensis was found to be related to the K channel-targeted I(2) conotoxin superfamily. Most putative Lophiotoma toxins are cysteine-rich polypeptides, with a significant fraction much larger (>80 amino acids) than the toxins from cone snails. A small number were not cysteine-rich but had hydrophobic amino acid clusters interspersed with arginine residues. This is only 1 of >10,000 different turrid venoms that needs to be characterized. From this study, a common origin with Conus for one family of putative turrid toxins is indicated.

Developmental expression of two Haliotis asinina hemocyanin isoforms.

Hemocyanins are large copper-containing respiratory proteins that play a role in oxygen transport in many molluscs. In some species only one hemocyanin isoform is present while in others two are expressed. The physiological relevance of these isoforms is unclear and the developmental and tissue-specific expression of hemocyanin genes is largely unknown. Here we show that two hemocyanin genes in the gastropod Haliotis asinina, which encode H. asinina hemocyanin (HaH1) and HaH2 isoforms, are developmentally expressed. These genes initially are expressed in a small number of mesenchyme cells at trochophore and pre-torsional veliger stages, with HaH1 expression slightly preceding HaH2. These cells largely are localized to the visceral mass, although a small number of cells are present in head and foot regions. Following metamorphosis the isoforms show overlapping as well as isoform-specific expression profiles, suggesting some degree of isoform-specific function.

Precursors of novel Gla-containing conotoxins contain a carboxy-terminal recognition site that directs gamma-carboxylation.

Vitamin K-dependent gamma-glutamyl carboxylase catalyzes the conversion of glutamyl residues to gamma-carboxyglutamate. Its substrates include vertebrate proteins involved in blood coagulation, bone mineralization, and signal transduction and invertebrate ion channel blockers known as conotoxins. Substrate recognition involves a recognition element, the gamma-carboxylation recognition site, typically located within a cleavable propeptide preceding the targeted glutamyl residues. We have purified two novel gamma-carboxyglutamate-containing conotoxins, Gla-TxX and Gla-TxXI, from the venom of Conus textile. Their cDNA-deduced precursors have a signal peptide but no apparent propeptide. Instead, they contain a C-terminal extension that directs gamma-carboxylation but is not found on the mature conotoxin. A synthetic 13-residue "postpeptide" from the Gla-TxXI precursor reduced the K(m) for the reaction of the Conus gamma-carboxylase with peptide substrates, including FLEEL and conantokin-G, by up to 440-fold, regardless of whether it was positioned at the N- or C-terminal end of the mature toxin. Comparison of the postpeptides to propeptides from other conotoxins suggested some common elements, and amino acid substitutions of these residues perturbed gamma-carboxylation of the Gla-TxXI peptide. The demonstration of a functional and transferable C-terminal postpeptide in these conotoxins indicates the presence of the gamma-carboxylation recognition site within the postpeptide and defines a novel precursor structure for vitamin K-dependent polypeptides. It also provides the first formal evidence to prove that gamma-carboxylation occurs as a post-translational rather than a cotranslational process.

Scallop DMT functions as a Ca2+ transporter.

We identified a DMT (divalent metal transporter) homologous protein that functions as a Ca(2+) transporter. Scallop DMT cDNA encodes a 539-amino-acid protein with 12 putative membrane-spanning domains and has a consensus transport motif in the fourth extracellular loop. Since its mRNA is significantly expressed in the gill and intestine, it is assumed that scallop DMT transports Ca(2+) from seawater by the gill and from food by the intestine. Scallop DMT lacks the iron-responsive element commonly found in iron-regulatory proteins, suggesting that it is free of the post-transcriptional regulation from intracellular Fe(2+) concentration. Scallop DMT distinctly functions as a Ca(2+) transporter unlike other DMTs, however, it also transports Fe(2+) and Cd(2+) similar to them.

cDNA cloning of an alginate lyase from abalone, Haliotis discus hannai.

An alginate lyase, termed HdAly in the present paper, was isolated from the hepatopancreas of abalone, Haliotis discus hannai, by ammonium sulfate fractionation, followed by TOYOPEARL CM-650M column chromatography. Enzymatic properties of HdAly were similar to those of previously reported Haliotis and Turbo poly(M) lyases, e.g., it preferentially degraded a poly(beta-D-mannuronate)-rich substrate with an optimal pH and temperature at pH 8.0 and 45 degrees C, respectively. In order to determine the primary structure of abalone lyase that is still poorly understood, cDNAs for HdAly were cloned by PCR from the abalone hepatopancreas cDNA library and sequenced. From the nucleotide sequences of the cDNAs, the sequence of 909 bp in total was determined, and the amino acid sequence of 273 residues was deduced from the translational region of 822 bp locating at nucleotide positions 27-848. The N-terminal region of 16 residues, except for the initiation Met in the deduced sequence, was regarded as the signal peptide since it was absent in the HdAly protein and showed high similarity to the consensus sequence for signal peptides of eukaryote secretary proteins. This suggests that HdAly is initially produced as a precursor possessing the signal peptide in hepatopancreatic cells and then secreted into digestive tract as the mature form. Thus, the mature HdAly was regarded to consist of 256 residues with the calculated molecular mass of 28895.5 Da. The amino acid sequence of HdAly showed 85 and 28% identity to those of Turbo cornutus alginate lyase SP2 and the C-terminal region of Chlorella virus lyase-like protein CL2, respectively, while it showed no significant identity to those of any bacterial alginate lyases. In order to provide the basis for the structure-function studies and various applications of the abalone lyase, a bacterial expression system was constructed by means of the HdAly-cDNA and pET-3a expression plasmid. Although the active recombinant HdAly was hardly produced at a cultivation temperature 37 degrees C in Escherichia coli BL21 (DE3), a small amount of soluble and active enzyme could be produced when the temperature was lowered to 19 degrees C.

cDNA cloning of two A-superfamily conotoxins from Conus striatus.

The full-length cDNAs of two A-superfamily conotoxins, kappaA-SIVA and alpha-SII, were respectively cloned and sequenced from Conus striatus using 3' RACE and 5' RACE. The cDNA of kappaA-SIVA encodes a precursor of 68 residues, including a signal peptide of 21 residues, a pro-peptide of 17 residues, and a mature peptide of 30 residues with an additional residue Gly which is prerequisite for the amidation of the preceding C-terminal Cys. The cDNA-deduced sequence of alpha-SII is composed of a signal peptide of 21 residues, a pro-peptide of 29 residues, a mature peptide of 19 residues and three additional residues Arg-Thr-Ile at the C-terminus. This tripeptide might be cleaved off by proteolytic processing. Although these two conotoxins belong to different families and target voltage-gated potassium channel and nicotinic acetylcholine receptor, respectively, they share the same signal sequence, and both are processed at the common signal site -X-Arg- immediately before the mature peptide sequences. The length of 3' untranslational region of alpha-conotoxin SII was extraordinarily large about 10 times longer than that of kappaA-SIVA with 770 and 75 bp, respectively. The elucidated cDNAs of these two toxins will facilitate a better understanding of the process of their post-translational modifications.

A new metallothionein gene from the giant keyhole limpet Megathura crenulata.

Metallothioneins (MTs) are small soluble proteins ubiquitously expressed in animals and plants. Different isoforms are present in deuterostomes and protostomes. They do not differ greatly in primary structure, but are clearly distinguishable. Here, I present the gene and the complete cDNA of a novel MT from the mollusk Megathura crenulata. This protein is closely related to the Cu-inducible MTs of the vineyard snail Helix pomatia, but has also some minor sequence features typical of Cd-inducible isoforms of H. pomatia and other molluscs. Overall, the deduced primary structure is similar to the known molluscan MTs, but in addition possesses an insertion of 5 amino acids not found in any other molluscan MTs, protostomic or deuterostomic MTs. In addition, a pentapeptide insertion, characteristic of mammalian MT-3 is present but it lacks the functional tetrapeptide CPCP within the beta-region of those MT-3 proteins that are known to suppress neuronal growth processes. The M. crenulata MT is a novel form of MT in comparison to all other known MTs. Possible functional aspects for this new MT are discussed.

Traceability of aquatic animals.

Effective methods of traceability are urgently required for use in research as well as in different types of aquaculture operations and to control trade in aquatic animals and products. In regard to the marking of fish, many different tagging methods have been described and the method to be used depends on the purpose and need for tagging. In contrast, for molluscs and crustaceans, only a few methods of marking such animals have been described, due to the practical difficulties. The authors first describe the different methods for tracing fish and fishery products, by means of external tags, such as Floy tags, Carlin tags and passive integrated transponder tags; chemical marking using inorganic substances such as silver nitrate or potassium nitrate, pigments, oxytetracycline, etc.; and several different types of electronic devices in which basic information such as the strain of fish, farm of origin or weight can be stored. Genetic traceability using deoxyribonucleic acid profiling is developing quite rapidly for cultured brood stocks and wild populations. This technique may be used with very high degrees of confidence to assign to or exclude animals or products from their claimed origin, paternity or strain, and may be used as evidence in court proceedings. The second section of this paper describes the traceability of live molluscs for restocking and for human consumption. In these applications, genetic markers have been demonstrated to be suitable. Mechanical tagging on a small scale for research purposes has also been used. Otherwise, the only means of tracing live molluscs are the movement documents and the labelling on boxes that certifies the origin of the commodity. The third section describes the methods available for tracing live and dead crustaceans. A large variety of physical tagging methods for decapod crustaceans is described, such as the injection of biological stains (fast green, Niagara sky blue, trypan red and blue) and external tags such as coloured streamer tags, wire tags and a variety of anchor tags. Furthermore, a number of different internal coding methods, such as the coded micro-wire tags and injected elastomer tags are discussed in detail. As is the case for fish, genetic molecular techniques are also applied in population studies of crustaceans; some of the molecular genetic methods are described. Prawns for human consumption are most frequently packed whole or as tails after the necessary sorting, washing and freezing and the only way of performing a traceback is through documents relating to movement, invoices, health certificates and labelling of the boxes. The minimum requirements for labelling would be the content of the packages, i.e. species, quantity, identification of the manufacturer (name and address), packing place, importer/exporter or vendor of the product, in addition to the loading bill number.

A comparative study of asymmetric migration events across a marine biogeographic boundary.

In many nonclonal, benthic marine species, geographic distribution is mediated by the dispersal of their larvae. The dispersal and recruitment of marine larvae may be limited by temperature gradients that can affect mortality or by ocean currents that can directly affect the movements of pelagic larvae. We focus on Point Conception, a well-known biogeographic boundary between the Californian and Oregonian biogeographic provinces, to investigate whether ocean currents affect patterns of gene flow in intertidal marine invertebrates. The predominance of pelagically dispersing species with northern range limits at Point Conception suggests that ocean currents can affect species distributions by erecting barriers to the dispersal of planktonic larvae. In this paper, we investigate whether the predominantly southward currents have left a recognizable genetic signature in species with pelagically dispersing larvae whose ranges span Point Conception. We use patterns of genetic diversity and a new method for inferring cladistic migration events to test the hypothesis that southward currents increase southward gene flow for species with pelagically dispersing larvae. We collected mitochondrial DNA (mtDNA) sequence data for the barnacles Balanus glandula and Chthamalus fissus and also reanalyzed a previously published mtDNA dataset (Strongylocentrotus purpuratus, Edmands et al. 1996). For all three species, our cladistic approach identified an excess of southward migration events across Point Conception. In data from a fourth species with nondispersing larvae (Nucella emarginata, Marko 1998), our method suggests that ocean currents have not played a role in generating genetic structure.

Amino acid sequence of a Ca(2+)-transporting ATPase from the sarcoplasmic reticulum of the cross-striated part of the adductor muscle of the deep sea scallop: comparison to serca enzymes of other animals.

The RT PCR approach was used to obtain the nucleotide sequence of the mRNA of a sarco/endoplasmic reticulum calcium transporting ATPase (SERCA) from the cross-striated (phasic) part of the adductor muscle of the deep sea scallop. Initially, degenerate primers based on consensus sequences among SERCAs and tryptic fragments of the scallop Ca-ATPase were used. The sequence was then extended using homologous primers and the 5' and 3' ends of the transcript determined by 5' and 3' RACE. The mRNA codes for a polypeptide chain 994 amino acid residues long (coded for by 2982 nucleotides) and has a 195 bp 5' untranslated region, with a 697 bp 3' untranslated region. The scallop enzyme shows strongest amino acid similarity to the SERCA enzyme of Loligo, followed by those of Drosophila and Artemia. It resembles the vertebrate SERCA3 in that it does not possess the phospholamban binding motif and so is unlikely to be regulated by protein kinase A mediated signals.

Sequence of the Octopus dofleini hemocyanin subunit: structural and evolutionary implications.

Sequencing of the subunit of the hemocyanin of Octopus dofleini has been completed from a cDNA library. This represents the first molluscan hemocyanin to be completely sequenced. The sequence determined is for one of the two distinguishable cDNAs which have been recognized for this protein. The protein subunit has 2896 amino acids and contains seven functional units, each carrying two sets of three invariant histidine residues constituting the binding sites (A and B) for two copper atoms. The accompanying paper identifies this site in the C-terminal functional unit (Odg). Differences in sequence for the two cDNAs, for the region in which both are available, are concentrated in the "linker regions" between functional units. The sequences of the seven units exhibit high similarity, averaging about 40% identity, with a concentration of conserved sequences in the region surrounding the copper binding sites. The sequences around the B-site show significant homology to the sequences of arthropod hemocyanins. Comparison of the functional unit sequences in terms of hydrophobicity and surface exposure profiles, as well as regions of probable secondary structure, indicate that all functional units probably have a common tertiary folding; the protein subunit is a string of similarly folded beads. A number of putative N-linked carbohydrate binding sites can be recognized in the sequence; one of these corresponds to the carbohydrate observed in the X-ray diffraction study of functional unit Odg as disclosed in the accompying paper. Phylogenetic analysis of the sequences of the O. dofleini functional units, and comparison with other available molluscan sequences indicates that the multi-domain subunit structure must have arisen over a relatively brief period, preceeding the differentiation of major molluscan types.