The Mini-hemoglobins in Neural and Body Wall Tissue of the Nemertean Worm, Cerebratulus lacteus*

Thomas L. VandergonDagger §, Claire K. RiggsDagger , Thomas A. GorrDagger , James M. Colacinoparallel , and Austen F. RiggsDagger **

From the Dagger  Department of Zoology, University of Texas, Austin, Texas 78712-1064 and the parallel  Department of Biological Sciences, Clemson University, Clemson, South Carolina 29631

    ABSTRACT
Top
Abstract
Introduction
Procedures
Results & Discussion
References

Hemoglobin (Hb) occurs in circulating red blood cells, neural tissue, and body wall muscle tissue of the nemertean worm, Cerebratulus lacteus. The neural and body wall tissue each express single major Hb components for which the amino acid sequences have been deduced from cDNA and genomic DNA. These 109-residue globins form the smallest stable Hbs known. The globin genes have three exons and two introns with splice sites in the highly conserved positions of most globin genes. Alignment of the sequences with those of other globins indicates that the A, B, and H helices are about one-half the typical length. Phylogenetic analysis indicates that shortening results in a small tendency of globins to group together regardless of their actual relationships. The neural and body wall Hbs in situ are half-saturated with O2 at 2.9 and 4.1 torr, respectively. The Hill coefficient for the neural Hb in situ, ~2.9, suggests that the neural Hb self-associates in the deoxy state at least to tetramers at the 2-3 mM (heme) concentration estimated in the cells. The Hb must dissociate upon oxygenation and dilution because the weight-average molecular mass of the HbO2 in vitro is only about 18 kDa at 2-3 µM heme concentration. Calculations suggest that the Hb can function as an O2 store capable of extending neuronal activity in an anoxic environment for 5-30 min.

    INTRODUCTION
Top
Abstract
Introduction
Procedures
Results & Discussion
References

Hemoglobin was first observed in neural tissue of invertebrates by Lankester (3) who noted the brilliant crimson color of the ganglia of the polychaetous annelid Aphrodite aculeata. Neural Hbs have since been recorded in or associated with nervous tissue of other annelids (4-7), molluscs (5, 8-12), arthropods (13, 14), and nemerteans (15, 16). The functional significance of Hb1 in neural tissue has been addressed in very few animals. Chalazonitis et al. (17) correlated the oxygenation state of neural Hb in Aplysia dipilans with the electrical activity of the neural ganglia and found that firing activity was proportional to the degree of oxygenation of the Hb. The neural tissue Hbs of the clams Tellina alternata and Spisula solidissima can extend the time of O2 delivery to nerves during anoxic periods by as much as 30 min by acting as an O2 store (12, 18).

Recently, DeWilde et al. (19) isolated and determined the amino acid sequence of the Hb from the neural tissue of A. aculeata. They found that the Hb was dimeric and had a relatively high O2 affinity and that the sequence and gene structure clearly showed it to be a member of the globin family.

Many nemertean worms express intracellular Hbs in red blood cells, body wall muscle tissue, and neural tissue. We report here the amino acid sequences, gene structure, and oxygen equilibria in situ for both the body wall and the neural tissue Hbs of the marine nemertean, Cerebratulus lacteus. We also propose a role of the neural hemoglobin as an oxygen store. We have used the amino acid sequences of the globins in maximum parsimony analyses to address possible phylogenetic relationships. These Hbs are very small, yet stable, unlike the artificially truncated mini-Mbs (20-22). This finding should make the C. lacteus Hbs particularly valuable for studies of folding and stability.

    EXPERIMENTAL PROCEDURES
Top
Abstract
Introduction
Procedures
Results & Discussion
References

Animals and Tissue Preparation

Eighty-four specimens of C. lacteus were purchased from the Department of Marine Resources, Marine Biological Laboratories, Woods Hole, MA, and maintained at 4 °C in the laboratory in 0.45 µm of filtered seawater. Animals used for oxygen-binding experiments were kept a maximum of 2 weeks. Tissue samples needed for protein, RNA, and DNA experiments were dissected from animals within 72 h of delivery, immediately frozen in a dry ice/ethanol bath, and stored at -80 °C. Tissues were prepared for microscopy by relaxing the animals in MgCl2 isotonic to seawater and fixing the head and trunk segments in 2.5% glutaraldehyde in 200 mM phosphate buffer, pH 7.4, made isotonic (0.58 M) by addition of NaCl. Segments were postfixed in 1.0% osmium tetroxide in the same buffer, dehydrated in an ethanol series, and embedded in Epon 812 resin. Ultrathin sections were stained with methylene blue and examined by light microscopy.

Isolation and Purification of Hemoglobins

Hb-containing tissues were ground in liquid N2, then transferred to 50 mM Tris acetate, 10 mM EDTA, pH 7.5, at 0 °C, followed by centrifugation at 12,000 × g for 10 min at 4 °C. Brain and lateral nerve tissue extracts are referred to as neural Hb, and body wall muscle extracts are referred to as body wall Hb. Tissue extracts were stored at 0 °C if used immediately or refrozen in liquid N2 and stored at -80 °C. Hbs were purified first by size-exclusion chromatography with G-75-SF Sephadex (Sigma) in a column (12 × 450 mm) equilibrated at 4 °C with 50 mM Tris acetate, 10 mM EDTA, pH 7.4 (measured at room temperature). Hb-containing fractions from this column were pooled, dialyzed against 10 mM Tris acetate, pH 8.0, and concentrated. The Hb was further purified by HPLC anion-exchange chromatography on a Synchropak AX300 column (4.6 × 250 mm), (SynChrom, Inc., Lafayette, IN), initially equilibrated with 10 mM Tris acetate, pH 8.0, at room temperature. Hbs were eluted from the column with a linear gradient from 0 to 50% of 10 mM Tris acetate, 500 mM sodium acetate, pH 8.0, at a flow rate of 0.5 ml/min. Hb fractions, identified by absorbance at 415 nm, were pooled, dialyzed against 10 mM Tris, pH 8.0, and lyophilized. These Hb samples were resuspended in 0.1% trifluoroacetic acid in water and further purified by RP-HPLC on a Synchropak RP300 C-18 column (4.6 × 250 mm) equilibrated with 0.1% trifluoroacetic acid in water at room temperature. Globins were eluted by establishing a series of linear gradients from 0 to 80% acetonitrile containing 0.1% trifluoroacetic acid over 70 min at a flow rate of 0.8 ml/min.

Neural tissue extracts from four animals were analyzed on a size-exclusion column, TSK-30 (4.6 × 300 mm, Bio-Rad), on an LDC Milton Roy CCM automated HPLC system. Extracts and molecular size standards (Sigma, catalog number MW-GF-70) were run at room temperature on the column equilibrated with 50 mM Tris acetate, pH 7.5, over 15 min at a flow rate of 1.0 ml/min. The standards were bovine serum albumin, ovalbumin, beta -lactoglobulin, horse heart myoglobin, and cytochrome c. Protein standards and Hb samples were monitored at 235 and 415 nm, respectively.

Amino Acid Composition and Sequence Analyses

Purified globins were hydrolyzed in 6 N HCl, sealed in vacuo in glass tubes, and hydrolyzed at 110 °C for 24 or 72 h. Amino acid compositions were determined with a Beckman model 121MB analyzer. Amino-terminal sequences of the globins were determined with a model 477A protein sequencer (Applied Biosystems, Inc.). Both instruments are at the University of Texas Microanalysis Facility.

Preparation and Amplification of Neural Globin mRNA

Neural tissue was ground in liquid N2 and dissolved in a solution containing 3 M LiCl, 6 M urea, and 0.2% sodium dodecyl sulfate at 0 °C. Brain homogenate was extracted with phenol/chloroform/isoamyl alcohol (25:24:1 v/v/v), and the RNA was precipitated with ethanol and resuspended in diethylpyrocarbonate-treated, sterile, deionized distilled water and stored at -80 °C. Poly(A)+ RNA was isolated from the total RNA by oligo(dT) affinity chromatography and reverse-transcribed using Moloney-murine leukemia virus reverse transcriptase (Promega). This was followed by extraction with phenol/chloroform/isoamyl alcohol, ethanol precipitation of the cDNA, resuspension in diethylpyrocarbonate-treated water, and storage at -80 °C. The cDNA was used as a template to amplify globin cDNA in the polymerase chain reaction (PCRTM, Hoffmann-La Roche). Primers were oligo(dT) and a redundant oligonucleotide based on the NH2-terminal residues of the globin (VNWAAV = 5'-GTXAA(T/C)TGGGCXGCXGT-3'). The reaction conditions were 100 mM Tris-HCl (pH 8.3 at 25 °C), 50 mM KCl, 2 mM MgCl2, 1 µg/µl bovine serum albumin, 0.25 µM oligo(dT), 1 µM redundant primer, 200 µM dNTPs, and 2.5 units of AmpliTaq DNA polymerase (Perkin-Elmer) in a volume of 100 µl. Reactants were heated to 98 °C for 2 min prior to the addition of enzyme, and amplification was performed for 30 cycles (94 °C, 1.5 min; 45 °C, 1 min; 72 °C, 1.5 min (cycles 1-10), and then 2.5 min (cycles 11-30)). PCR products were cloned into the pCR-1000 vector (TA Cloning System, Invitrogen, San Diego).

Genomic DNA Isolation and Sequencing

A genomic DNA library was constructed in the Lambda Dash II vector (Stratagene) from BamHI-restricted, size-selected genomic DNA that came from the sperm of a single animal. The library was screened with a 345-bp AccI fragment from the neural globin cDNA PCR product (see "Results"). A genomic clone containing a 25-kb insert that hybridized to the globin cDNA probe was isolated from the library, and several enzyme restriction fragments were subcloned into the SmaI site of pBluescript (Stratagene) for sequencing. Fig. 1 shows the positions of the various clones used to determine the organization and sequence of the genes for neural and body wall globins.


View larger version (13K):
[in this window]
[in a new window]
 
Fig. 1.   Comparison of the number of nucleotides (nt) in exons, introns, and adjacent 5' and 3' non-coding regions in C. lacteus globin genes. Boxes represent exons, and horizontal lines represent 5' and 3' non-coding regions and introns. The body wall globin sequence was determined from the p1-3 and TA9 PCR-generated constructs and from genomic clone fragments ClaI-3, ClaI-7, HpaI-6 and SphI-7. The neural globin sequence was determined from the pP8 and p4-7 PCR-generated constructs and from the p42 4.5-kb HpaI fragment of a genomic DNA clone. The diagram is not drawn to scale.

Body Wall Globin-- The 6-kb HpaI clone (HpaI-6) contained 5.1 kb of a body wall globin gene that was truncated in the first intron of the gene and lacked exon I. Other clones (SphI-7, 7 kb; ClaI-3, 3 kb; and ClaI-7, 7 kb) provided additional sequence information (see Fig. 1). The missing exon I region was isolated from genomic DNA by the use of two PCR amplification protocols. The first utilized a redundant primer within exon I and a non-redundant primer from within intron I of the body wall globin gene to amplify the intervening sequence from a template of XhoI-digested genomic DNA yielding fragment TA9 (see Fig. 1). The reaction mixture consisted of 20 mM Tris-HCl (pH 8.8 at 25 °C), 10 mM KCl, 10 mM (NH4)2SO4, 2 mM MgSO4, 0.1% Triton X-100, 2.4 µM redundant primer, 0.2 µM non-redundant primer, 400 µM dNTPs, and 2.5 units of AmpliTaq DNA polymerase. The region from intron I into the region 5' to exon I of the body wall globin gene was amplified in a second PCR reaction using a set of two nested primers from intron I and two primers from the PCR in vitro Single Site Amplification and Cloning System (catalog number TAK R015, Panvera Corp., Madison, WI) in a reaction mix given by the manufacturer. HindIII-digested genomic DNA was used as the template for the reaction which yielded fragment p1-3 (see Fig. 1). All PCR products were subcloned into pUC19 or pCRTMII vector (TA cloning system, Invitrogen, San Diego, CA) for sequencing performed by the dideoxy method using Sequenase 2.0 (U. S. Biochemical Corp.).

Neural Globin-- A second screening of the genomic DNA library, using the 581-bp cDNA for neural globin, yielded a clone containing an HpaI digestion fragment of 4.5 kb that contained the complete sequence of the neural globin gene except for the first two codons of exon 1 and the adjoining 5' upstream sequence (clone p42, see Fig. 1). The missing region was obtained with the PCR in vitro Single Site Amplification and Cloning System with a pair of non-redundant, nested primers located in the region of exon 2 and a sample of PstI-digested genomic DNA yielding fragment pP8 (see Fig. 1). A second pair of non-redundant, nested primers from the exon 3 sequence was used to amplify exon 3 and intron 2 of the neuronal globin gene from a sample of XbaI-digested genomic DNA yielding fragment p4-7 (see Fig. 1). The PCR products and fragment p42 were cloned and sequenced as described above.

Phylogenetic Analysis of Globin Sequences

Alignment-- Amino acid sequences of globins were obtained from the literature. Sequences of 22 globins from vertebrate and invertebrate animals (including Cerebratulus), bacteria, protists, algae, plants, and fungi were aligned by eye in a series of steps as follows. First, the 84 conserved residue positions identified by Kapp et al. (23) as shared by all members of the globin family were used for a core alignment that included the conserved CD1 Phe and F8 His. The template of Vinogradov et al. (24) was used to facilitate the identification of conserved hydrophilic or hydrophobic residues at external and internal locations, respectively. This template is a version of template I of Bashford et al. (25) extended to include invertebrates. This core comprises the shared parts of helices A, B, C, E, F, G, and H. The second step was to fill in the interhelical regions by using known x-ray structures of individual globins where possible (10 out of 22 globins). The helical regions, as identified in the x-ray analyses, were aligned with a minimum number of insertions or deletions in order to maintain the overall tertiary structure. Gaps were confined almost entirely to the loops between helices. Finally, the external/buried/internal classification for individual residues of Fermi and Perutz (26), originally developed for vertebrate myoglobins and hemoglobins, was taken as an additional template. This template includes both helical and interhelical regions, unlike those of Bashford and Vinogradov (24, 25), which carry only helical information. The Fermi template helps to identify and preserve the chemical patterns such as the arrangement of polar and non-polar amino acids in interhelical loops. All interhelical regions were expanded just enough to accommodate the longest sequence in these regions. The resulting alignment needed a total of 179 positions for the 22 globins with ~110 positions common to a majority of sequences.

Analysis-- The aligned sequences were analyzed for relationships by employing the maximum parsimony method. All programs involved in this analysis are part of the PHYLIP Program Package, version 3.5, of Felsenstein (27). The original sequences were bootstrapped either 100 or 1,000 times using the program SEQBOOT (28). The replicate data sets so obtained were analyzed by the program PROTPARS which calculated a corresponding number of unrooted parsimony trees with the use of both the global rearrangement option and the option for a randomized input order of sequences. Finally, the CONSENSE program was used to generate a single majority-rule consensus tree from that population of 100 or 1,000 parsimony trees, with calculated bootstrap support values at each node. The human alpha -globin sequence was taken as outgroup for rooting the tree because the monophyly of vertebrate myoglobins and hemoglobins is established.

Oxygen Equilibria

The oxygen binding properties of the neural and body wall Hbs were determined in situ by thin layer microspectrophotometry at 15 °C (29, 30). Hb solutions for in vitro experiments were prepared by homogenizing tissues on ice in 50 mM Tris buffers at various pH values followed by centrifugation at 12,000 × g at 4 °C for 10 min to remove cellular debris. Tissue for in situ experiments was prepared by cutting a piece of the brain or body wall muscle and rinsing it in a solution of 0.45 µm of filtered seawater at 0 °C buffered with 50 mM Tris at the appropriate pH and containing 20 mM KCN to inhibit cellular respiration. Thin sections were removed from the pieces with a razor blade and placed in a fresh buffer at 0 °C. All oxygen binding experiments were performed at 15 °C. All pH measurements were made at room temperature. Absorbance data from 12 wavelength (nm) pairs (410/420, 425, 430, 435, 440; 415/425, 430, 435, 440; 420/430, 435; 425/435) were used in a two-wavelength method for determining fractional saturation values. The tissue slice pathlength, measured optically, varied between 80 and 240 µm. Mean absorbance changes were 0.050 and 0.183 at 415 and 430 nm, respectively. Between 5 and 7 experiments were performed under each set of conditions. Typically, 40-60 data points were averaged for each of 6 oxygen pressures in each experiment. The combined data for the experiments under each set of conditions were fitted by nonlinear least squares (31) directly to the Hill equation, y = KPn/(1 + KPn) to estimate n and P50 (= nth root of K-1). Comparison of the absorption spectra before and after an experiment suggests that Met-Hb formation during an experiment is no more than about 5%.

    RESULTS AND DISCUSSION
Top
Abstract
Introduction
Procedures
Results & Discussion
References

Hemoglobin Expression-- Examination of C. lacteus neural ganglia and lateral nerve cords in fresh tissue shows that the cells surrounding axon bundles are uniformly bright red (HC in Fig. 2, A and B). Microspectrophotometry of this tissue gives the characteristic spectrum of HbO2, which is not detectable within the axonal bundles (AX), although some Hb-containing cells appear to have processes that enter the axonal bundles. Cells within the body wall also express Hb, but at differing concentrations, being most concentrated in the trunk region surrounding the gonadal sacs where it may aid in the delivery of oxygen to the developing gametes.


View larger version (108K):
[in this window]
[in a new window]
 
Fig. 2.   A, transverse section through the cerebral ganglia of C. lacteus; B, cross-section of the lateral nerve cord in the mid-trunk region. Both photomicrographs show central axonal material (AX) surrounded by nerve and glial cell bodies. The latter regions have a dark red color in fresh sections and are the hemoglobin-containing (HC) regions. PR, proboscis; CG, cerebral ganglia; RH, rhynchocoel; FG, foregut.

The concentration of Hb in the neural ganglia of the brain was estimated by dissecting out the brains from two animals and analyzing the entire extracts by size-exclusion HPLC. Integration of the absorption peaks at 415 nm yielded 1.1 and 1.7 nmol of Hb for the two extracts. We assumed a millimolar extinction coefficient of 131 at 415 nm, the value for human HbO2 (32). The volume of the brain shown in Fig. 2A would be completely contained within a sphere of radius 600 µm and volume 0.9 µl. This, then, provides an upper limit to the brain volume and provides a minimal concentration of 1.2-1.9 mM. The actual concentration must be substantially higher because the volume of 0.9 µl includes parts of the foregut, rhynchocoel, and the axonal tissue itself, none of which have Hb. We conclude that the concentration is at least 2-3 mM and could be significantly higher. This estimate is comparable to the 2.1-2.8 mM values for myoglobin found in the muscles of seals, porpoises, and tuna (33-35).

Structure of Globins and Globin Genes-- A single major globin protein fraction was isolated from neural tissue extracts in a three-step chromatographic procedure (Fig. 3, A-C). The body wall Hb also gave a single major peak in the first two chromatographic steps, but RP-HPLC of the Hb fraction showed several minor peaks (Fig. 3D), raising the possibility of additional Hbs in this tissue. NH2-terminal sequences of neural and body wall globins were determined (43 and 45 residues, respectively) and found to share 79% amino acid sequence identity (Fig. 4). The first eight residues of both globins are identical (VNWAAVVD), allowing us to construct a redundant oligonucleotide primer directed at both globin genes (see "Experimental Procedures"). Use of this primer and oligo(dT) permitted the amplification and cloning of a 581-bp cDNA fragment that encodes the complete sequence of a globin from the neural tissue (Fig. 5). Five clones were completely sequenced and found to be identical in the coding region, but they possessed some differences in the 3' non-coding region. The amino acid sequence determined by direct sequencing was DFY at residues 9-11 but was EL-STOP in the PCR-generated cDNA clones (Figs. 4 and 5). We conclude that the EL-STOP sequence is the result of errors in amplification and/or cloning because the neural globin genomic DNA sequence (Fig. 6) matches the direct determination of residues 1-42. An additional difference is that the cDNA sequence has Leu at position 56, whereas the genomic DNA gives Ile for this position.


View larger version (15K):
[in this window]
[in a new window]
 
Fig. 3.   A, elution profile of C. lacteus whole neural tissue lysate separated on a G-75-SF column equilibrated with 50 mM Tris acetate, 10 mM EDTA, pH 7.4, at 4 °C. Two equivalent runs are shown, monitored at 280 or 415 nm. B, elution profile of the neural Hb gel filtration fraction (A) separated by anion-exchange HPLC on a Synchropak AX300 column initially equilibrated with 10 mM Tris acetate, pH 8.0. Hemoglobins were eluted from the column by establishing a linear gradient from 0 to 50% 10 mM Tris acetate, 0.5 M sodium acetate, pH 8.0, at a flow rate of 0.5 ml/min. C, elution profile of the neural Hb anion-exchange fraction (B) separated by RP-HPLC on a Synchropak RP300 C-18 column equilibrated with 0.1% trifluoroacetic acid in water. Globins were eluted from the column by establishing a series of linear gradients from 0 to 80% acetonitrile containing 0.1% trifluoroacetic acid at a flow rate of 0.8 ml/min. D, elution profile of C. lacteus body wall Hb anion-exchange fraction separated by RP-HPLC. All Hb- and globin-containing chromatographic fractions are marked with a bar.


View larger version (14K):
[in this window]
[in a new window]
 
Fig. 4.   Comparison of C. lacteus globin sequences. The underlined sequences were determined both by direct sequence analysis and from cDNA (neural) or genomic DNA (neural and body wall). Boxed residues 9-11 (neural) were not correct in the cDNA (see text) but were obtained both by direct sequencing and from genomic DNA. Vertical bars mark the differences between the two globins. Asterisks mark the universally conserved CD1 (F), distal E7 (Q), and proximal F8 (H) residues.


View larger version (34K):
[in this window]
[in a new window]
 
Fig. 5.   Sequence of the 581-bp cDNA encoding the C. lacteus neural globin. The first 17 nucleotides (underlined) correspond to the redundant primer (see "Experimental Procedures"). Splice junction positions predicted from sequence alignment and subsequently verified from genomic clones are marked with arrows. The single-letter amino acid designations for residues of the globin chain are listed below their corresponding codons. A putative polyadenylation signal in the 3' non-coding region is boxed. Asterisks mark nucleotides that are C or T in genomic DNA. Codons 9-11 code for DFY in the genomic DNA and are in agreement with direct sequencing (see text).


View larger version (69K):
[in this window]
[in a new window]
 
Fig. 6.   Genomic DNA sequence encoding the three exons of C. lacteus neural globin. Splice junction sequences are underlined. Numbers to the right are the corresponding nucleotide numbers in the GenBankTM file. Coding regions are in capital letters, and single-letter amino acid designations are used for the globin sequence.

The amino acid sequence of the body wall globin was deduced from the genomic DNA (Fig. 7). This sequence matches the first 45 NH2-terminal residues determined directly (Fig. 4). The gene contains three exons separated by two introns. Intron 1 separates exons 1 and 2, splitting a Glu codon (residue 13) corresponding to position B12 in both neural and body wall globins. Intron 2 separates exons 2 and 3 between codons for His-80 (G6) and Asn-81 (G7) in neural globin and Thr-80 (G6) and His-81 (G7) in body wall globin. The positions of the 5' and 3' splice junctions are identical in both genes (Figs. 6 and 7) and match those found for typical globins (36). The two globin sequences are 78% identical and are clearly derived by duplication of an ancestral gene.


View larger version (70K):
[in this window]
[in a new window]
 
Fig. 7.   Genomic DNA sequence encoding the three exons of C. lacteus body wall globin. Splice junction sequences and a putative polyadenylation signal sequence are underlined. Numbers to the right are the corresponding nucleotide numbers in the GenBankTM file. Coding regions are in capital letters, and single-letter amino acid designations are used for the globin sequence.

The amino acid compositions of neural and body wall globins are compared in Table I with those deduced from DNA. The agreement is within 5-7% for the neural globin except for Cys and Trp, which are partially destroyed by acid hydrolysis, and Val. The low Val value presumably results from difficulty in cleaving the Val-Val bond between residues 6 and 7. The analysis of body wall globin shows greater deviations with unexplained high Ser and low Ala values.

                              
View this table:
[in this window]
[in a new window]
 
Table I
Amino acid composition of C. lacteus globins
Data are the averages of molar ratios of amino acids recovered after hydrolysis in 6 N HCl in vacuo for 24 and 72 h. Values in parentheses are derived from the nucleotide sequence.

Protein Structure-- The chains of neural and body wall Hbs have only 109 residues which makes them the smallest naturally occurring Hbs known. Unlike the proteolytically generated mini-Mb of 108 residues (20-22), the mini-Hbs of C. lacteus are stable. Paramecium Hb (37) is slightly larger (116 residues) than the nemertean Hbs and autoxidizes (38) 10-fold faster at pH 8 than sperm whale Mb (153 residues). The C. lacteus Hbs are sufficiently stable to permit the HbO2 to survive the lengthy purification procedures and to make possible the measurement of O2 equilibria over a period of several hours. In contrast, the mini-myoglobins (20-22) are rapidly oxidized in air, and O2 equilibrium measurements are impossible. We examine here how the C. lacteus chains differ from those of other Hbs and how these chains may achieve stability despite their short length.

The Alignment Reveals Amino- and Carboxyl-terminal Deletions-- The alignment of 22 globins (Fig. 8) shows an impressive conservation of hydrophobic residues at 37 internal positions throughout the globin family (i.e. the 34 columns of residues with blue backgrounds plus the CD1, E7, and F8 columns in red in Fig. 8) This pattern of 37 non-polar amino acids in globins at solvent-inaccessible locations, identified by Gerstein et al. (68), was applied to a sequence alignment by Kapp et al. (23). It allows a much improved refinement of the alignment of all globins. These residues at solvent-inaccessible positions are a major factor for maintaining the universal globin fold. The C. lacteus sequences have 34 of these 37 positions retained and lack only A8, B6, and H19 because of deletions (see Fig. 8). Kapp et al. (23) have also identified 84 positions conserved in all 700 globin sequences examined (boxed regions in Fig. 8). According to our alignment, the C. lacteus globins have all of these conserved positions except for 7 residues deleted in regions B5-B7 and H16-H19. Despite their shortness, the globins of this nemertean have clearly retained most of the 37 solvent-inaccessible and the 84 core residue positions. However, large deletions appear at both ends of C. lacteus globins. The alignment shows that greater parts of the A and B helix, as well as from the H helix, are missing. Still, all six A-helical core residues are maintained in the nemertine globins. A tryptophan occupying position 3 of this conserved core forms a hydrophobic anchor that makes contacts with the E, G, and H helices. This position, Trp-3 in C. lacteus globins, corresponds to human alpha chain residue Trp-14 (A12) which forms contacts with Ala-63 (E12), Leu-66 (E15), Thr-67 (E16), Val-70 (E19), Leu-105 (G12), and Leu-109 (G16) (see Ref. 69). The alignment suggests that very similar contacts are likely in C. lacteus neural Hb: Val-49 (E12), Ile-52 (E15), Asn-53 (E16), Leu-56 (E19), Leu-86 (G12), and Cys-90 (G16). The alignment preserves the propensity for a hydrogen bond-forming residue at position E16 in 15/22 of the chains shown in Fig. 8. The six-residue core of the A helix identified by Kapp et al. (23) is part of the 13-residue core found in six high resolution x-ray studies (70). The six-residue core is conserved in all globins including that of C. lacteus. We conclude that C. lacteus globins have retained the minimal core of residues in helix A required for the Mb fold.


View larger version (72K):
[in this window]
[in a new window]
 


View larger version (83K):
[in this window]
[in a new window]
 


View larger version (58K):
[in this window]
[in a new window]
 
Fig. 8.   Sequence alignment. Sequences: data on hemoglobins and genes from the following species were used. a, human: alpha -globin (26); b, Escherichia coli: hmp gene, globin domain (39); c, Vitreoscilla stercocaria: VtHb gene (40-42); d, Saccharomyces cerevisiae: flavo-Hb, globin domain (43); e, Candida norvegensis: flavo-Hb, globin domain (44); f, Glycine max (soybean): Lb a (45, 46); g, Parasponia andersonii: Lb I (47); h, Chlamydomonas eugametos: LI 637 gene (48); i, Nostoc commune: glbN gene (49); j, Tetrahymena pyriformis: Hb I (50, 51); k, Paramecium caudatum (37, 52); l, Cerebratulus lacteus: neural Hb (this study); m, Pseudoterranova decipiens: repeat 1 (53); n, Scapharca inaequivalvis: homodimeric Hb I (54, 55); o, Lumbricus terrestris: Hb chain b (splice junction positions determined by Donahue (85) (56); p, Glycera dibranchiata: monomeric (M) Hb II (all residues are as shown in Ref. 57 except residue 57 which is Asp in chemically determined sequence; no direct evidence exists for Ala in this position.) (57); q, Urechis caupo: Hb F1 (58, 59); r, Chironomus thummi thummi: monomeric Hb IV (60-62); s, Artemia spec.: T4 domain (63, 64); t, Caudina (Molpadia) arenicola: monomeric Hb C (structure of gene determined by Thomas (86) (65); u, Petromyzon marinus; Hb II (87) (major component; 66, 67); v, Physeter catodon (sperm whale): myoglobin (26). Templates: The templates of Fermi and Perutz (26) and Vinogradov et al. (24) are given above and below the sequences, respectively (see text). Positions are designated as follows: e, external; b, buried or surface crevice (Fermi only); i, internal; i, i, i, internal with volume variability low (>100 Å3), medium (~40-100 Å3), and severe (<40 Å3), respectively (Vinogradov only). Classification of size restrictions follow Table IV of Bashford et al. (25); a,= any residue (Vinogradov only). The 84 positions always conserved (Kapp et al. (23)) are boxed. Dots above each sequence mark every 10th residue. Color code: red bars, CD1, E7 (distal), and F8 (proximal); blue bars, the 34 conserved, solvent-inaccessible positions; brown background, exon/intron splice junctions of the genes that occur either within a single codon (phase 1 or 2) or between two codons (phase 0); white background, intronless genes as in bacteria, yeasts, or Chironomus Hb IV; yellow, splice junctions expected on the basis of homology but that have yet to be determined. Lines above sequences mark helical regions identified by x-ray crystallography in the following Hbs: human alpha , Vitreoscilla, Glycine, Scapharca, Glycera, Urechis, Chironomus, Caudina, Petromyzon, and Physeter (sperm whale). The x-ray structure of the monomeric Chironomus Hb III (62) was used for the Hb IV sequence shown. Ambiguous assignment to one or another of two adjacent helices or between a helix and an interhelical region occurs in the following positions: Ukechis, residues 37, 58, 78, and 81, and Chironomus residues 31, 72, and 76.

The truncated A helix (residues 1-7) may be particularly important for the stability of C. lacteus Hbs. The artificial mini-myoglobin (20-22) is unstable and rapidly autoxidizes. It comprises residues 32-139 of sperm whale Mb and so lacks the entire A helix and has a shortened COOH-terminal H helix. The C. lacteus Hbs have a similar deletion from the H helix. The COOH-terminal residue of C. lacteus Hbs corresponds to residue 138 of sperm whale Mb; the mini-Mb extends only one residue further. It appears, therefore, that the presence of an A helix, albeit shortened, in C. lacteus Hb and its absence in mini-myoglobin may be correlated with the difference in stability.

Most globins have a B helix that is almost twice as long as that in C. lacteus. How might the apparent gap between helices A and B in C. lacteus globins be accommodated? The gap between residues 8 and 9 of the C. lacteus globins corresponds to residues B1 to B7 in the human alpha  chain. Chironomus IV globin has a gap that corresponds to human alpha  chain helical residues B1 and B4. The high resolution x-ray structure of the very similar Chironomus III globin (61) shows one way in which such apparent gaps can be readily accommodated. Residues 12-19 in Chironomus III Hb form a 310 helix. Eight residues of a 310 helix would be about 4 Å longer than an eight-residue alpha  helix. Thus, helical segments could "stretch" to maintain essential contacts. There is no necessary correlation between globin chain length and stability, provided the essential features of the myoglobin fold are conserved.

The alignment of the C. lacteus sequence (Fig. 8) places some limitations on the minimal requirements for folding. Jennings and Wright (71) have identified 20 residues in helices A, B, G, and H in sperm whale apoMb for which the amide protons become protected within 5 ms during the initial folding and helix formation. C. lacteus globin lacks residues in five of these helical positions: A7-A9 and H19-H20. Therefore, residues at these positions do not appear to be required for folding.

Heme Pocket-- The structural stability of Hbs depends strongly on the affinity for heme. Ferrous Mbs and Hbs show little tendency to lose heme because of the covalent bond between the proximal His and the Fe2+ atom. This bond is much weaker in ferric Hb (Met-Hb). The lower retention of the hemin in Met-Mb and Met-Hb then depends more on the hydrophobicity of the heme pocket and, in aquo-Met-Hb, on the water molecule that is anchored to the iron and to distal residues (72). Additional stabilization is contributed by electrostatic contacts with the propionic acid groups on the heme. Table II shows the residues of the heme pocket in sperm whale Mb and in the beta  chain of human Hb. Without exception, all residues in the positions that are hydrophobic in both these chains are also hydrophobic in C. lacteus globins.

                              
View this table:
[in this window]
[in a new window]
 
Table II
Comparison of heme pocket residues of C. lacteus neural and body wall globins with those determined for sperm whale myoglobin and human beta -globin

The presence of Thr at E11 in C. lacteus Hbs is unusual and may contribute substantially to the stability of the Hbs. The mutation Val-68 right-arrow Thr (E11) in sperm whale Mb (72) results in a very large increase in stability by slowing the rate of hemin loss from the ferric Mb because the water molecule bound to the Fe3+ of the hemin is stabilized by hydrogen bonds to both the distal His-64 (E7) and Thr-68 (E11). The distal His at E7 of most vertebrate Hbs is replaced by Gln in the Hbs of C. lacteus. This substitution is not, however, associated with any significant increase in the rate of hemin loss in sperm whale Mb (72). So, the corresponding substitution in C. lacteus globins should not be expected to alter their rate of hemin loss significantly. Glutamine is found in the distal E7 position in virtually all of the prokaryotic and protist globins as well as many invertebrate globins. This suggests that Gln might be the ancestral residue at this position.

Additional stabilization of the hemin in C. lacteus globins is provided by the interaction of residue Lys-47 (E10) corresponding to human beta  Lys-66 (E10) which is known to form an electrostatic link with the propionic acid group of pyrrole IV of the heme (26). The Arg at F7 and Lys or Arg at C7 also should provide stabilization in the globins by propionic acid group interaction.

Hydrophilicity-- The hydrophilicity patterns of Fig. 9 show the close similarity of the body wall and neural Hbs of C. lacteus. The neural Hb has, however, substantially more hydrophilicity than the body wall Hb. This difference can be quantified by examining the 24 positions where the amino acid residues differ (Fig. 4). The sums of the hydrophilicity values at these positions for body wall and neural Hbs are 7.9 and -39.4, respectively.2 Why the Hbs should differ in this way is puzzling. Perhaps it reflects differences in the nature of the subunit interfaces involved in self-association.


View larger version (28K):
[in this window]
[in a new window]
 
Fig. 9.   Comparison of the hydrophilicity patterns of neural and body wall globins of C. lacteus. Values were determined using the Kyte-Doolittle procedure (73) calculated by MacVector 3.5 (International Biotechnologies, Inc., New Haven, CT) and a hydrophilicity window size of 5 residues. Hydrophilic residues have positive values. Helical designations are those shown in Fig. 8.

Fig. 10 shows a residue-by-residue comparison of the hydrophobicity of body wall Hb and sperm whale Mb. The alignment is identical to that in Fig. 8. The high degree of correspondence evident between the two Hbs provides substantial support for the overall alignment. The least satisfactory correspondence is in helix A where Trp of position 14 (sperm whale Mb) is used as an anchor (see arrow). Comparison of the two patterns in the EF region suggests an alternative to that in Fig. 8: the Ser-Ala-Asp- sequence (residues 59-61 of C. lacteus neural or body wall globin) corresponds better in hydrophilicity to sperm whale Mb if it is moved two residues to the right (see Fig. 8).


View larger version (40K):
[in this window]
[in a new window]
 
Fig. 10.   Comparison of the hydrophobicity patterns of C. lacteus body wall Hb (solid line) and sperm whale Mb (dashed line). The alignment is identical to that in Fig. 8. Arrows mark anchoring alignment residues at positions 14 (Trp), 43 (CD1 Phe), 64 (distal heme residue, His, or Gln), and 93 (proximal His). Bull's eye" circles mark positions with the same hydrophobicity in each chain. Gaps where the alignment provides no residue in C. lacteus globin that corresponds to a residue in sperm whale Mb occur at positions 19-26, 45, 50, 79-80, 84-85, 96, and 124-127. Hydrophilic residues are negative. The use of a moving average of a 5-residue window in Fig. 9 accounts for the apparent qualitative differences from this figure where hydrophobicities of individual residues are plotted.

Phylogenetic Analysis of Globin Sequences-- The 22 aligned globin sequences (Fig. 8) were analyzed as described under "Experimental Procedures." Fig. 11 summarizes the results. Consensus tree A is based on 1,000 bootstrap replicates, and the consensus trees B-D each are based on 100 replicates of the original data set of the 22 sequences.


View larger version (34K):
[in this window]
[in a new window]
 
Fig. 11.   The results of the phylogenetic analysis described in the text. A, all sequences are full-length; 1,000 bootstrap replicates. B, all sequences are shortened to include only the 93 positions E1-H19 (100 bootstrap replicates). C, only the 84 common-core positions identified by Kapp et al. (23) are included (100 replicates). D, all sequences are full-length except that of Glycera which has been shortened by deleting positions that correspond to those missing in C. lacteus and in the other four truncated sequences of the naturally short globins (Glycera*).10

Tree A, rooted with human alpha -globin, shows the consensus topology for all sequences (179 positions, as in Fig. 8). This tree has several important features as follows. 1) Three major clusters of globins are present (besides vertebrate globins), and their monophyletic origins are well supported (nodal bootstrap confidence is shown in bold, >= 95%, yeasts and bacteria (99% confidence), plants (99%), and the truncated globins of protozoans, the green alga Chlamydomonas, and the cyanobacterium Nostoc (100%). This finding of three clades is in excellent agreement with the results of a phylogenetic analysis of Moens et al. (74). The surprising similarities between yeast and bacterial globins and between all four truncated sequences were explained by these authors with two horizontal gene transfers as follows: one from a common ancestor of both yeasts to the prokaryotic ancestor of Escherichia, Vitreoscilla, and Alcaligenes; and the other from the common ancestor of ciliate protozoans and the green algae to the cyanobacterial lineage that evolved into Nostoc. The latter hypothesis is consistent with the extensive intracellular symbiosis between ciliates and algae that occur. 2) A gradient of decreasing vertebrate similarities from human alpha  chain to Petromyzon, the sperm whale Mb, and Caudina globins exists (see also trees B-D) that leads into the pool of all other globins consistent with the present understanding of the relationship of echinoderms to vertebrates. 3) A complete lack of resolution for the interrelationships of the major clusters and of invertebrate globins (extremely low bootstrap confidences) is evident. 4) Cerebratulus globin (boxed) is apparently grouped to the clade of the other truncated globins with a bootstrap value at 40%.

Recent data (75) on RNA sequences show that nemerteans should be classified with protostomous coelomates and are evolutionarily distant from protists, algae, and bacteria. How, then, can one explain the grouping of C. lacteus globins with those of more primitive organisms? Is it possible that the association is caused by the common shortness of their globins? Several modifiers of parsimony topologies are known. Factors such as the input order of sequences, the use of different outgroups or ingroup representatives (sampling bias), the selective attraction of long branches to each other, and the different evolutionary rates among the branches can all result in false groupings. The length of sequences is certainly another variable that can give biased trees. Gaps are treated as the 21st character state in the program PROTPARS that we used, with each deletion being three changes away from all other character states (i.e. codons for all 20 amino acids). Therefore, short sequences with more matches of gaps in homologous positions to each other than to longer, non-deleted proteins, have a surplus of unique 21st state identities that might contribute to group them together to the exclusion of longer sequences. The effect of sequence length on the apparent topology is examined in trees B-D.

Tree B shows the result of shortening all sequences to 93 positions (E1-H19). This truncation removes length-variable NH2- and COOH termini as well as the AB, CD, and D helix regions and minimizes the number of gaps involved in the analysis. The strength of the monophyletic support of the yeast/bacteria, the plant, and the "truncation" clades was not affected by this procedure, but the bootstrap value for the Cerebratulus truncation linkage increased from ~40 to ~59%. This finding shows that both the bootstrap confidence and the relative arrangement of lineages can be manipulated by artificial changes in sequence length.

Tree C takes this approach one step further. Here, the topology was determined by using only the 84 common core positions described by Kapp et al. (23) (see boxes in Fig. 8). This procedure confines the sequence data to the parts A10-A15, B5-C7, CD1, CD2, E1-E20, F1-F19, G4-G17, and H6-H19. This set of positions completely lacks any interhelical sequences except for CD1 and CD2. The result of this is that the bootstrap confidence for two major units (plants and bacteria/yeasts) drops below the 95% significance cutoff (bootstrap values in lightface), and the apparent relationship of Cerebratulus to other short sequences disappears. The sequences have become too short and too conservative to reveal any phylogenetic differentiation. This finding suggests that the most phylogenetically informative positions are apparently located within interhelical and not in helical regions, especially in the EF, FG, and GH loops (compare trees B and C). This conclusion makes clear that phylogenetic differentiation using globins relies on the interhelical segments rather than on the functionally important but conservative blocks of helical residues.

Tree D displays the result of artificially shortening Glycera globin from 147 to 109 residues with no change in the other sequences. The shortening of Glycera globin was accomplished by introducing deletions in positions homologous to those found in the five short globins. The result of this procedure is to move the Glycera globin even closer to the other short globins except for Cerebratulus which is now displaced to an insignificant level of correspondence. A similar result was obtained by artificial truncation of Urechis globin which also displaced Cerebratulus and gave the following grouping3: (Cerebratulus (Urechis*, short sequences)45)40. Shortening Chironomus yielded a ((Cerebratulus, Chironomus*)43 short sequences)49 topology (data not shown). These results demonstrate a difficulty in comparing homologous sequences of significantly different lengths. We conclude that the apparent relationship of the Cerebratulus sequence to other short sequences is a result of a common shortness, not common ancestry. Furthermore, the alignment in Fig. 8 makes clear that deletions in four miniglobins (Chlamydomonas, Nostoc, Tetrahymena, and Paramecium) occurred independently of those in Cerebratulus because the positions of the deletions differ.

Hemoglobin Function-- Data on oxygen equilibria measured in situ by microspectrophotometry (Fig. 12, A and B, and Table III) show that the neural tissue Hb has a higher oxygen binding affinity (P50 approx 2.9 torr) than that of the body wall Hb (P50 approx 4.1 torr) at 15 °C. Oxygen binding by the neural Hb in situ exhibits little or no dependence upon extracellular pH within the narrow range of pH 7.3-7.9 (Delta log P50/Delta pH < 0.1). Intracellular pH was not measured. Therefore, these results could mean that the intracellular pH was relatively constant and independent of extracellular pH.


View larger version (14K):
[in this window]
[in a new window]
 
Fig. 12.   Hill plots of oxygen equilibrium data of (A), neural Hb in situ, pH 7.3 and (B), body wall Hb, open circles, in vitro, pH 7.4 (B); filled circles, in situ, pH 7.3. The solid lines are those calculated as described in the text.

                              
View this table:
[in this window]
[in a new window]
 
Table III
Summary of O2 equilibrium data for hemoglobins of C. lacteus determined by microspectrophotometry

Our results show that the oxygen equilibria for body wall Hb in vivo and in vitro at pH 7.4 are identical within the errors of measurement. The spectrum of the body wall Hb in vitro indicates that Met-Hb formation is no more than about 5% during the long 2-5-h O2 binding experiments which indicates substantial stability. Under similar conditions, the truncated mini-Mbs would be 100% oxidized. Some Met-Hb was present initially, but we cannot accurately quantify this in the tissues because the extinction coefficients are unknown and because other absorbing substances are present. Although no in vitro O2 binding measurements were made with the neural Hb, it proved to be relatively stable throughout the lengthy chromatographic preparation.

Oxygen binding by both Hbs is clearly cooperative. The mean Hill coefficient, 2.9, for O2 binding by neural Hb at an estimated concentration of ~2-3 mM (heme) in situ is only slightly lower than that for human Hb. This means that the neural Hb must self-associate, probably at least to tetramers. The estimated molecular weight obtained by size-exclusion HPLC for HbO2 is only 17,700 ± 2,000 obtained at a mean concentration across the peak of 2-5 µM. This value is the weight-average molecular weight expected in a mixture of monomers (Mr approx  12,000) and dimers if the weight fraction of dimers were 0.48. These data can be interpreted in terms of an oxygenation-linked and concentration-dependent dissociation of presumed Hb tetramers: Hb4 right-left-arrows 2 Hb2 right-left-arrows 4 Hb. Oxygenation-linked dissociation processes are widespread in Hbs of both vertebrates and invertebrates (76-79) where the oligomeric Hbs are always more tightly associated in the deoxygenated state.

The cooperative behavior of the body wall Hb is similar to that of the neural Hb, but the data in Fig. 12B show a distinct difference. The slopes of the Hill plots for the body wall Hb in situ and in vitro show a pronounced decrease to 1.0-1.4 above 50% saturation (log (y/(1-y)) = 0). Such a drop may result from dissociation of the oligomer, Hb4 right-left-harpoons  2 Hb2 right-left-harpoons  4 Hb, at earlier stages of oxygenation than occurs with the neural Hb. The maximal Hill coefficient approaches 2.9 for both in situ and in vitro body wall Hbs between y = 0.3 and y = 0.5.

Cerebratulus lives in intertidal and subtidal sediment in shallow water and mud flats along the North Atlantic coast, quite often deep below the surface in the anoxic zone where it burrows in search of prey. The relatively high concentration of high affinity Hb within the neural tissue suggests that the Hb also may function as an O2 store when the animal encounters anoxic conditions during burrowing. Oxygen storage has been demonstrated for the neural Hb of the clams T. alternata and S. solidissima, which live in very similar habitats (12, 18), and has been suggested for the neural Hb of A. aculeata (19). An oxygen store could also be useful under normoxic conditions to augment the oxygen supply to nervous tissues (29, 80). Such tissues may exhibit bursts of activity during which oxygen demand could temporarily exceed oxygen transport from the environment. A nearby auxiliary oxygen source would be valuable under such circumstances. We used a simple model to calculate the capacity of the potential O2 store in the lateral nerves of Cerebratulus. Storage capacity was then divided by the potential O2 demand to give estimates of the maximum time the stored O2 would last. We assumed that no O2 was lost by diffusion (see last paragraph) and that the O2 consumption rate was constant and homogeneous throughout the nerves.

The lateral nerves of Cerebratulus are roughly cylindrical with nerve cell bodies and Hb restricted to an outer shell and axons to a central core (Fig. 2B). The volume of a cylindrical shell, 1.0 cm long, determined from measurements on calibrated photomicrographs, was approximately 6.3 × 10-4 cm3. Multiplying the shell volume by the heme concentration, 2 to 3 µmol/ml, gives the oxygen storage capacity 1.3-1.9 nmol or 29.1-42.6 nl of O2 STPD. Oxygen demand was estimated from literature values for O2 consumption rates of nemertean whole animals and resting or stimulated nerves (5, 7, 81, 82). These values range from 1,500 nl of O2 STPD/g-min (resting) to 6,000 nl of O2 STPD/g-min (stimulated). The volume of the entire nerve cord (shell plus core) was 9.5 × 10-4 cm3, giving O2 consumption rates between 1.4 and 5.7 nl of O2 STPD/min for the 1.0-cm model. The store could, therefore, theoretically supply O2 for periods ranging from 5 min (small O2 store and high O2 demand) to 31 min (large store and low demand). These durations are comparable to those found for T. alternata neural Hb which may extend normal neuronal activity for up to 30 min during exposure to anoxic conditions (12, 18).

The large Hill coefficient of both neural and body wall Hb will cause them to unload a majority of bound O2 within a narrow range of pO2 near their P50 values (12, 83, 84). This would tend to stabilize the tissue pO2 during unloading and could be useful in providing optimal conditions for nerve and muscle cell activity. These factors could also serve to prevent or reduce diffusive O2 loss from the neural tissue O2 store, because the surrounding body wall tissue, with its lower affinity Hb, would tend to act as an additional O2 source for the neural tissue. Thus, our simplifying assumption for the model calculations, that no O2 is lost by diffusion, may have some basis in fact. These considerations and the computations suggest that Cerebratulus neural Hb could function as a useful O2 store to extend or augment neuronal activity under both anoxic and normoxic conditions. The presence of high affinity, low molecular weight Hbs in the body wall and nervous tissues could also serve to facilitate O2 transport to and within these tissues.

    ACKNOWLEDGEMENTS

We thank Peter Smith, Elizabeth Balser, and Pat Behrens for assistance with the experiments; Alfred Wheeler and Edward Ruppert for the use of equipment; and Janet Young, Gwen Gage, Kristina Schlegel, and David Green for assistance with graphics. We thank Sandra Smith of the University of Texas Microsequencing Facility for automated sequencing of proteins.

    FOOTNOTES

* This research was supported in part by Welch Foundation Grant F-213 (to A. F. R.), National Science Foundation Grants DMB-8600614 (to J. M. C.), MCB 9205764, 9511759 and 972385 (to A. F. R.), and National Institutes of Health Grant GM35847 (to A. F. R.). Preliminary accounts of this work have been presented (1, 2).The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

The nucleotide sequence(s) reported in this paper has been submitted to the GenBankTM/EMBL Data Bank with accession number(s) AF033001 and AF033002.

§ Present address: Natural Sciences Division, Pepperdine University, 24255 Pacific Coast Hwy., Malibu, CA 90263. Tel.: 310-456-4326; Fax: 310-456-4758; E-mail: tvanderg{at}pepperdine.edu.

Present address: Hematology-Oncology Division, LMRC-2, Brigham and Women's Hospital, 221 Longwood Ave., Boston, MA 02115.

** To whom correspondence should be addressed.

1 The abbreviations used are: Hb, hemoglobin; Mb, myoglobin; HPLC, high performance liquid chromatography; RP-HPLC, reverse phase HPLC; STPD, standard temperature and pressure, dry; PCR, polymerase chain reaction; bp, base pair; kb, kilobase pair.

2 We use the Kyte-Doolittle "hydropathic index" for these values. Hydrophilic residues are negative in this scale (73). The MacVector Program reverses the sign to provide a hydrophilicity scale (Fig. 9). The original Kyte-Doolittle scale (here called hydrophobicity) is used in Fig. 10.

3 Only the starred sequences are shortened. Bootstrap confidences are given as superscripts.

    REFERENCES
Top
Abstract
Introduction
Procedures
Results & Discussion
References

  1. Vandergon, T. L., and Colacino, J. M. (1989) Am. Zool. 29, 59
  2. Vandergon, T. L., Riggs, C. K., Colacino, J. M., and Riggs, A. (1990) Am. Zool. 31, 72
  3. Lankester, E. R. (1872) Proc. R. Soc. Lond. B Biol. Sci. 21, 70-81
  4. Manwell, C. (1960) Arch. Biochem. Biophys. 89, 194-201[Medline] [Order article via Infotrieve]
  5. Wittenberg, J. B., Brown, P. K., and Wittenberg, B. A. (1965) Biochim. Biophys. Acta 109, 518-529[Medline] [Order article via Infotrieve]
  6. Wittenberg, B. A., Briehl, R. W., and Wittenberg, J. B. (1965) Biochem. J. 96, 363-371[Medline] [Order article via Infotrieve]
  7. Weber, R. E. (1978) in Physiology of Annelids (Mills, P. J., ed), pp. 393-446, Academic Press, London
  8. Kennedy, D. (1960) J. Gen. Physiol. 44, 277-299[Abstract/Free Full Text]
  9. Arvanitaki, A., and Chalazonitis, N. (1960) Bull. Inst. Oceanogr. 57, 1-83
  10. Strittmatter, P., and Burch, H. B. (1963) Biochim. Biophys. Acta 78, 562-563[CrossRef][Medline] [Order article via Infotrieve]
  11. Schindelmeiser, I., Kuhlmann, E., and Nolte, A. (1979) Comp. Biochem. Physiol. 64, 149-154
  12. Kraus, D. W., and Colacino, J. M. (1986) Science 232, 90-92
  13. Fox, H. M. (1955) Proc. R. Soc. Lond. B Biol. Sci. 143, 203-214[Medline] [Order article via Infotrieve]
  14. Fox, H. M. (1957) Nature 179, 148[Medline] [Order article via Infotrieve]
  15. Vernet, G. (1973) Comp. Biochem. Physiol. 44, 1053-1058
  16. Varndell, I. M. (1980) J. Exp. Mar. Biol. Ecol. 45, 157-172
  17. Chalazonitis, N., Gola, M., and Arvanitaki, A. (1965) C. R. Soc. Biol. 159, 240-244
  18. Doeller, J. E., and Kraus, D. W. (1988) Biol. Bull. 174, 67-76
  19. Dewilde, S., Blaxter, M., Van Hauwaert, M. L., Vanfleteren, J., Esmans, E. L., Marden, M., Griffon, N., and Moens, L. (1996) J. Biol. Chem. 271, 19865-19870[Abstract/Free Full Text]
  20. De Sanctis, G., Falcioni, G., Giardina, B., Ascoli, F., and Brunori, M. (1988) J. Mol. Biol. 200, 725-733[Medline] [Order article via Infotrieve]
  21. De Sanctis, G., Ascoli, F., and Brunori, M. (1994) Proc. Natl. Acad. Sci. U. S. A. 91, 11507-11511[Abstract/Free Full Text]
  22. DeSanctis, G., Falcioni, G., Grenolli, F., Desideri, A., Polizio, F., Giardina, B., Ascoli, F., and Brunori, M. (1992) J. Mol. Biol. 222, 637-643
  23. Kapp, O. H., Moens, L., Vanfleteren, J., Trotman, C. N. A., Suzuki, T., and Vinogradov, S. N. (1995) Protein Sci. 4, 2179-2190[Abstract/Free Full Text]
  24. Vinogradov, S. N., Walz, D. A., Pohajdak, B., Moens, L., Kapp, O. H., Suzuki, T., and Trotman, C. N. A. (1993) Comp. Biochem. Physiol. 106, 1-2[CrossRef]
  25. Bashford, D., Chothia, C., and Lesk, A. M. (1987) J. Mol. Biol. 196, 199-216[Medline] [Order article via Infotrieve]
  26. Fermi, G., and Perutz, M. F. (1981) in Atlas of Molecular Structures 2: Haemoglobin and Myoglobin (Phillips, D. C., and Richards, F. M., eds), p. 26, Clarendon Press, Oxford
  27. Felsenstein, J. (1993) PHYLIP, version 3.5, distributed by author, Dept. of Genetics, Univ. of Washington, Seattle
  28. Felsenstein, J. (1985) Evolution 39, 783-791
  29. Colacino, J. M., and Kraus, D. W. (1984) Comp. Biochem. Physiol. 79, 363-369
  30. Vandergon, T. L., and Colacino, J. M. (1989) Comp. Biochem. Physiol. 94, 31-39
  31. Bevington, P. R. (1969) Data Reduction and Error Analysis for the Physical Sciences, McGraw-Hill Inc., New York
  32. Van Assendelft, O. W. (1970) Spectrophotometry of Haemoglobin Derivatives, p. 57, Royal VanGorcum, Assen
  33. Andersen, H. T. (1996) Physiol. Rev. 46, 212-243
  34. Blessing, M. H. (1971) Comp. Biochem. Physiol. 41, 475-480
  35. Wittenberg, B. A., and Wittenberg, J. B. (1990) Annu. Rev. Biophys. Biophys. Chem. 19, 217-241[CrossRef][Medline] [Order article via Infotrieve]
  36. Hardison, R. C. (1991) in Evolution at the Molecular Level (Selander, R. K., Clark, A. G., and Whittam, T. S., eds), pp. 272-333, Sinauer Associates, Inc., Sunderland, MA
  37. Iwaasa, H., Takagi, T., and Shikama, K. (1989) J. Mol. Biol. 208, 355-358[Medline] [Order article via Infotrieve]
  38. Tsubamoto, Y., Matsuoka, A., Yusa, K., and Shikama, K. (1990) Eur. J. Biochem. 193, 55-59[Abstract]
  39. Vasudevan, S. G., Armarego, W. L. F., Shaw, D. C., Lilley, P. E., Dixon, N. E., and Poole, R. K. (1991) Mol. Gen. Genet. 226, 49-58[Medline] [Order article via Infotrieve]
  40. Wakabayashi, S., Matsubara, H., and Webster, D. A. (1986) Nature 322, 481-483[Medline] [Order article via Infotrieve]
  41. Khosla, C., and Bailey, J. E. (1988) Mol. Gen. Genet. 214, 158-161[Medline] [Order article via Infotrieve]
  42. Tarricone, C., Galizzi, A., Ascenzi, P., and Bolognesi, M. (1997) Structure 5, 497-507[Medline] [Order article via Infotrieve]
  43. Zhu, H., and Riggs, A. F. (1992) Proc. Natl Acad. Sci. U. S. A. 89, 5015-5019[Abstract]
  44. Iwaasa, H., Takagi, T., and Shikama, K. (1992) J. Mol. Biol. 227, 948-954[CrossRef][Medline] [Order article via Infotrieve]
  45. Hyldig-Nielsen, J. J., Jenssen, E. O., Paludan, K., Wiborg, O., Garret, R., Jørgensen, G. P., and Marcker, K. A. (1982) Nucleic Acids Res. 10, 689-701[Abstract]
  46. Ollis, D. L., Appleby, C. A., Colman, P. M., Cutten, A. E., Guss, J. M., Venkatappa, M. P., and Freeman, H. C. (1983) Aust. J. Chem. 36, 451-468
  47. Landsmann, J., Dennis, E. S., Higgins, T. J. V., Appleby, C. A., Kortt, A. A., and Peacock, W. J. (1986) Nature 324, 166-168
  48. Couture, M., Chamberland, H., St-Pierre, B., Lafontaine, J., and Guertin, M. (1994) Mol. Gen. Genet. 243, 184-197
  49. Potts, M., Angeloni, S. V., Ebel, R. E., and Bassam, D. (1992) Science 256, 1690-1692[Medline] [Order article via Infotrieve]
  50. Iwaasa, H., Takagi, T., and Shikama, K. (1990) J. Biol. Chem. 265, 8603-8609[Abstract/Free Full Text]
  51. Takagi, T., Iwaasa, H., Yuasa, H., Shikama, K., Takemasa, T., and Watanabe, Y. (1993) Biochim. Biophys. Acta 1173, 75-78[Medline] [Order article via Infotrieve]
  52. Yamauchi, K., Ochiai, T., and Usuki, I. (1992) Biochim. Biophys. Acta 1171, 81-87[Medline] [Order article via Infotrieve]
  53. Dixon, B., Walker, B., Kimmins, W., and Pohajdak, B. (1992) J. Mol. Evol. 35, 131-136[Medline] [Order article via Infotrieve]
  54. Petruzzelli, R., Goffredo, B. M., Barra, D., Bossa, F., Boffi, A., Verzili, D., Ascoli, F., and Chiancone, E. (1985) FEBS Lett. 184, 328-332[CrossRef][Medline] [Order article via Infotrieve]
  55. Royer, W. E., Jr., Hendrickson, W. A., and Chiancone, E. (1989) J. Biol. Chem. 264, 21052-21061[Abstract/Free Full Text]
  56. Fushitani, K., Matsuura, M. S. A., and Riggs, A. F. (1988) J. Biol. Chem. 263, 6502-6517[Abstract/Free Full Text]
  57. Arents, G., and Love, W. E. (1989) J. Mol. Biol. 210, 149-161[Medline] [Order article via Infotrieve]
  58. Garey, J. R., and Riggs, A. F. (1986) J. Biol. Chem. 261, 16446-16450[Abstract/Free Full Text]
  59. Kolatkar, P. R., Hackert, M. L., and Riggs, A. F. (1994) J. Mol. Biol. 237, 87-97[CrossRef][Medline] [Order article via Infotrieve]
  60. Kleinschmidt, T., and Sgouros, J. G. (1987) Biol. Chem. Hoppe-Seyler 368, 579-615[Medline] [Order article via Infotrieve]
  61. Steigemann, W., and Weber, E. (1979) J. Mol. Biol. 127, 309-338[Medline] [Order article via Infotrieve]
  62. Antoine, M., and Niessing, J. (1984) Nature 310, 795-798
  63. Trotman, C. N. A., Manning, A. M., Moens, L., and Tate, W. P. (1991) J. Biol. Chem. 266, 13789-13795[Abstract/Free Full Text]
  64. Moens, L., Vanfleteren, J., De Baere, I., Jellie, A. M., Tate, W., and Trotman, C. N. A. (1992) FEBS Lett. 312, 105-109[CrossRef][Medline] [Order article via Infotrieve]
  65. Mitchell, D. T., Kitto, G. B., and Hackert, M. L. (1995) J. Mol. Biol. 251, 421-431[CrossRef][Medline] [Order article via Infotrieve]
  66. Hombrados, I., Rodewald, K., Neuzil, E., and Braunitzer, G. (1983) Biochimie (Paris) 65, 247-257[Medline] [Order article via Infotrieve]
  67. Honzatko, R. B., Hendrickson, W. A., and Love, W. E. (1985) J. Mol. Biol. 184, 147-164[Medline] [Order article via Infotrieve]
  68. Gerstein, M., Sonnhammer, E. L. L., and Chothia, C. (1994) J. Mol. Biol. 236, 1067-1078[Medline] [Order article via Infotrieve]
  69. Sack, J. S., Andrews, L. C., Magnus, K. A., Hanson, J. C., Rubin, J., and Love, W. E. (1978) Hemoglobin 2, 153-169[Medline] [Order article via Infotrieve]
  70. Aronson, H.-E. G., Royer, W. E., and Hendrickson, W. A. (1994) Protein Sci. 3, 1706-1711[Abstract/Free Full Text]
  71. Jennings, P. A., and Wright, P. E. (1993) Science 262, 892-896[Medline] [Order article via Infotrieve]
  72. Hargrove, M. S., Wilkinson, A. J., and Olson, J. S. (1996) Biochemistry 35, 11300-11309[CrossRef][Medline] [Order article via Infotrieve]
  73. Kyte, J., and Doolittle, R. F. (1982) J. Mol. Biol. 157, 105-132[Medline] [Order article via Infotrieve]
  74. Moens, L., Vanfleteren, J., Van de Peer, Y., Peeters, K., Kapp, O., Czelusniak, J., Goodman, M., Blaxter, M., and Vinogradov, S. N. (1996) Mol. Biol. Evol. 13, 324-333[Abstract]
  75. Turbeville, J. M., Field, K. G., and Raff, R. A. (1992) Mol. Biol. Evol. 9, 235-249[Abstract]
  76. Cobb, J. A., Manning, D., Kolatkar, P. R., Cox, D. J., and Riggs, A. F. (1992) J. Biol. Chem. 267, 1183-1189[Abstract/Free Full Text]
  77. Chiancone, E., Vecchini, P., Verzili, D., Ascoli, F., and Antonini, E. (1981) J. Mol. Biol. 152, 577-592[Medline] [Order article via Infotrieve]
  78. Roberts, M. S., Terwilliger, R. C., and Terwilliger, N. B. (1984) Comp. Biochem. Physiol. 77, 237-243
  79. Bonaventura, J., and Kitto, G. B. (1973) in Comparative Physiology (Bolis, L., Schmidt-Nielsen, K., and Maddresll, S. H. P., eds), pp. 493-507, Elsevier/North-Holland Biomedical Press, Amsterdam
  80. Atkinson, H. J. (1975) J. Exp. Biol. 62, 1-9[Abstract]
  81. Poluhowich, J. J. (1970) Comp. Biochem. Physiol. 36, 817-821
  82. Ritchie, J. M., and Straub, R. W. (1979) J. Physiol. (Lond.) 287, 315-327[Medline] [Order article via Infotrieve]
  83. Vandergon, T. L., and Colacino, J. M. (1991) Physiol. Zool. 64, 1561-1577
  84. Colacino, J. M., Brannan, J. E., and Fields, M. A. (1987) J. Theor. Biol. 126, 437-447[Medline] [Order article via Infotrieve]
  85. Donahue, R. A. (1992) The Determination of the cDNA Sequences for chains b and d of the Hemoglobin of Lumbricus terrestris, Master's thesis, University of Texas, Austin
  86. Thomas, P.-W. (1994) Genomic Analysis and Evolutionary Studies on an Echinoderm Hemoglobin Family.Ph.D. dissertation, University of Texas, Austin
  87. Qiu, Y. (1997) Genes and Subunit Association of Hemoglobins in the Lamprey, Petromyzon marinus.Ph.D. dissertation, University of Texas, Austin


Copyright © 1998 by The American Society for Biochemistry and Molecular Biology, Inc.