The amino acid sequence of hemoglobin III from the symbiont-harboring clamLucina pectinata

The cytoplasmic hemoglobin III from the gill of the symbiont-harboring clamLucina pectinata consists of 152 amino acid residues, has a calculated M_m of 18,068, including heme, and has N-acetyl-serine as the N-terminal residue. Based on the alignment of its sequence with other vertebrate and nonvertebrate globins, it retains the invariant residues Phe45 at position CD1 and His98 at the proximal position F8, as well as the highly conserved Trp16 and Pro39 at positions A12 and C2, respectively. The most likely candidate for the distal residue at position E7 is Gln66.Lucina hemoglobin III shares 95 identical residues with hemoglobin II (J. D. Hockenhull-Johnsonet al., J. Prot. Chem. 10, 609–622, 1991), including Tyr at position B10, which has been shown to be capable of entering the distal heme cavity and placing its hydroxyl group within a 2.8 Å of the water molecule occupying the distal ligand position, by modeling the hemoglobin II sequence using the crystal structure of sperm whale metmyoglobin. The amino acid sequences of the twoLucina globins are compared in detail with the known sequences of mollusc globins, including seven cytoplasmic and 11 intracellular globins. Relative to 75% homology between the twoLucina globins (counting identical and conserved residues), both sequences have percent homology scores ranging from 36–49% when compared to the two groups of mollusc globins. The highest homology appears to exist between theLucina globins and the cytoplasmic hemoglobin ofBusycon canaliculatum.     

The amino acid sequence of hemoglobin II from the symbiont-harboring clamLucina pectinata

The cytoplasmic hemoglobin II from the gill of the clamLucina pectinata consists of 150 amino acid residues, has a calculatedM _m of 17,476, including heme and an acetylated N-terminal residue. It retains the invariant residues Phe 44 at position CD1 and His 65 at the proximal position F8, as well as the highly conserved Trp 15 at position A12 and Pro 38 at position C2. The most likely candidate for the distal residue at position E7, based on the alignment with other globins, is Gln 65. However, optical and EPR spectroscopic studies of the ferri Hb II (Kraus, D. W., Wittenberg, J. B., Lu, J. F., and Peisach, J.,J. Biol. Chem. 265, 16054–16059, 1990) have implicated a tyrosinate oxygen as the distal ligand. Modeling of theLucina Hb II sequence, using the crystal structure of sperm whale aquometmyoglobin, showed that Tyr 30 substituting for the Leu located at position B10 can place its oxygen within 2.8 Å of the water molecule occupying the distal ligand position. This structural alteration is facilitated by the coordinate mutation of the residue at position CD4, from Phe 46 in the sperm whale myoglobin sequence to Leu 47 inLucina Hb II.     

The amino acid sequence of hemoglobin II from the symbiont-harboring clamLucina pectinata

The cytoplasmic hemoglobin II from the gill of the clamLucina pectinata consists of 150 amino acid residues, has a calculatedM _m of 17,476, including heme and an acetylated N-terminal residue. It retains the invariant residues Phe 44 at position CD1 and His 65 at the proximal position F8, as well as the highly conserved Trp 15 at position A12 and Pro 38 at position C2. The most likely candidate for the distal residue at position E7, based on the alignment with other globins, is Gln 65. However, optical and EPR spectroscopic studies of the ferri Hb II (Kraus, D. W., Wittenberg, J. B., Lu, J. F., and Peisach, J.,J. Biol. Chem. 265, 16054–16059, 1990) have implicated a tyrosinate oxygen as the distal ligand. Modeling of theLucina Hb II sequence, using the crystal structure of sperm whale aquometmyoglobin, showed that Tyr 30 substituting for the Leu located at position B10 can place its oxygen within 2.8 Å of the water molecule occupying the distal ligand position. This structural alteration is facilitated by the coordinate mutation of the residue at position CD4, from Phe 46 in the sperm whale myoglobin sequence to Leu 47 inLucina Hb II.     

Amino acid sequence of myoglobin from the chitonLiolophura japonica and a phylogenetic tree for molluscan globins

Myoglobin was isolated from the radular muscle of the chitonLiolophura japonica, a primitive archigastropodic mollusc.Liolophura contains three monomeric myoglobins (I, II, and III), and the complete amino acid sequence of myoglobin I has been determined. It is composed of 145 amino acid residues, and the molecular mass was calculated to be 16,070 D. The E7 distal histidine, which is replaced by valine or glutamine in several molluscan globins, is conserved inLiolophura myoglobin. The autoxidation rate at physiological conditions indicated thatLiolophura oxymyoglobin is fairly stable when compared with other molluscan myoglobins. The amino acid sequence ofLiolophura myoglobin shows low homology (11–21%) with molluscan dimeric myoglobins and hemoglobins, but shows higher homology (26–29%) with monomeric myoglobins from the gastropodic molluscsAplysia, Dolabella, andBursatella. A phylogenetic tree was constructed from 19 molluscan globin sequences. The tree separated them into two distinct clusters, a cluster for muscle myoglobins and a cluster for erythrocyte or gill hemoglobins. The myoglobin cluster is divided further into two subclusters, corresponding to monomeric and dimeric myoglobins, respectively.Liolophura myoglobin was placed on the branch of monomeric myoglobin lineage, showing that it diverged earlier from other monomeric myoglobins. The hemoglobin cluster is also divided into two subclusters. One cluster contains homodimeric, heterodimeric, tetrameric, and didomain chains of erythrocyte hemoglobins of the blood clamsAnadara, Scapharca, andBarbatia. Of special interest is the other subcluster. It consists of three hemoglobin chains derived from the bacterial symbiont-harboring clamsCalyptogena andLucina, in which hemoglobins are supposed to play an important role in maintaining the symbiosis with sulfide bacteria.     

Amino acid sequence of myoglobin from the chitonLiolophura japonica and a phylogenetic tree for molluscan globins

Myoglobin was isolated from the radular muscle of the chitonLiolophura japonica, a primitive archigastropodic mollusc.Liolophura contains three monomeric myoglobins (I, II, and III), and the complete amino acid sequence of myoglobin I has been determined. It is composed of 145 amino acid residues, and the molecular mass was calculated to be 16,070 D. The E7 distal histidine, which is replaced by valine or glutamine in several molluscan globins, is conserved inLiolophura myoglobin. The autoxidation rate at physiological conditions indicated thatLiolophura oxymyoglobin is fairly stable when compared with other molluscan myoglobins. The amino acid sequence ofLiolophura myoglobin shows low homology (11–21%) with molluscan dimeric myoglobins and hemoglobins, but shows higher homology (26–29%) with monomeric myoglobins from the gastropodic molluscsAplysia, Dolabella, andBursatella. A phylogenetic tree was constructed from 19 molluscan globin sequences. The tree separated them into two distinct clusters, a cluster for muscle myoglobins and a cluster for erythrocyte or gill hemoglobins. The myoglobin cluster is divided further into two subclusters, corresponding to monomeric and dimeric myoglobins, respectively.Liolophura myoglobin was placed on the branch of monomeric myoglobin lineage, showing that it diverged earlier from other monomeric myoglobins. The hemoglobin cluster is also divided into two subclusters. One cluster contains homodimeric, heterodimeric, tetrameric, and didomain chains of erythrocyte hemoglobins of the blood clamsAnadara, Scapharca, andBarbatia. Of special interest is the other subcluster. It consists of three hemoglobin chains derived from the bacterial symbiont-harboring clamsCalyptogena andLucina, in which hemoglobins are supposed to play an important role in maintaining the symbiosis with sulfide bacteria.     

Amino Acid Sequence, Spectral, Oxygen-Binding, and Autoxidation Properties of Indoleamine Dioxygenase-Like Myoglobin from the Gastropod Mollusc Turbo cornutus

Myoglobin was isolated from the radular muscle of the archaeogastropod mollusc Turbo cornutus (Turbinidae). This myoglobin is a monomer carrying one protoheme group; the molecular mass was estimated by SDS–PAGE to be about 40 kDa, 2.5 times larger than that of usual myoglobin. The cDNA-derived amino acid sequence of 375 residues was determined, of which 327 residues were identified directly by chemical sequencing of internal peptides. The amino acid sequence of Turbo myoglobin showed no significant homology with any other usual 16-kDa globins, but showed 36% identity with the myoglobin from Sulculus diversicolor (Haliotiidae) and 27% identity with human indoleamine 2,3-dioxygenase, a tryptophan-degrading enzyme containing heme. Thus, the Turbo myoglobin can be counted among the myoglobins which evolved from the same ancestor as that of indoleamine 2,3-dioxygenase. The absorbance ratio of to CT maximum (/CT) of Turbo metmyoglobin was 17.8, indicating that this myoglobin probably possesses a histidine residue near the sixth coordination position of heme iron. The Turbo myoglobin binds oxygen reversibly. Its oxygen equilibrium properties are similar to those of Sulculus myoglobin, giving P ₅₀ = 3.5 mm Hg at pH 7.4 and 20°C. The pH dependence of autoxidation of Turbo oxymyoglobin was quite different from that of mammalian myoglobin, suggesting a unique protein folding around the heme cavity of Turbo myoglobin. A kinetic analysis of autoxidation indicates that the amino acid residue with pK _a = 5.4 is involved in the reaction. The autoxidation reaction was enhanced markedly at pH 7.6, but not at pH 5.5 and 6.3 in the presence of tryptophan. We suggest that a noncatalytic binding site for tryptophan, in which several dissociation groups with pK _a 7.6 are involved, remains in Turbo myoglobin as a relic of molecular evolution.     

Crystal structure of the carbon monoxide complex of human cytoglobin

Cytoglobin (Cgb) is a vertebrate heme‐containing globin‐protein expressed in a broad range of mammalian tissues. Unlike myoglobin, Cgb displays a hexa‐coordinated (bis‐hystidyl) heme iron atom, having the heme distal His81(E7) residue as the endogenous sixth ligand. In the present study, we crystallized human Cgb in the presence of a reductant Na₂S₂O₄ under a carbon monoxide (CO) atmosphere, and determined the crystal structure at 2.6 Å resolution. The CO ligand occupies the sixth axial position of the heme ferrous iron. Eventually, the imidazole group of His81(E7) is expelled from the sixth position and swings out of the distal heme pocket. The flipping motion of the His81 imidazole group accompanies structural readjustments of some residues (Gln62, Phe63, Gln72, and Ser75) in both the CD‐corner and D‐helix regions of Cgb. On the other hand, no significant structural changes were observed in other Cgb regions, for example, on the proximal side. These structural alterations that occurred as a result of exogenous ligand (CO) binding are clearly different from those observed in other vertebrate hexa‐coordinated globins (mouse neuroglobin, Drosophila melanogaster hemoglobin) and penta‐coordinated sperm whale myoglobin. The present study provides the structural basis for further discussion of the unique ligand‐binding properties of Cgb. Proteins 2011. © 2011 Wiley‐Liss, Inc.     

Amino acid sequence of myoglobin from the molluscBursatella leachii

The complete amino acid sequence of myoglobin from the triturative stomach of gastropodic molluscBursatella leachii has been determined. It is composed of 146 amino acid residues, is acetylated at the N-terminus, and contains a single histidine residue at position 95 which corresponds to the heme-binding proximal histidine. The E7 distal histidine, which is conserved widely in myoglobins and hemoglobins, is replaced by valine inBursatella myoglobin. The amino acid sequence ofBursatella myoglobin shows strong homology (73–84%) with those ofAplysia andDolabella myoglobins.     

X-Ray crystallographic structural characteristics of Arabidopsis hemoglobin I and their functional implications

Genome of the model dicot flowering plant, Arabidopsis thaliana, a popular tool for understanding molecular biology of plant physiology, encodes all three classes of plant hemoglobins that differ in their sequence, ligand binding and spectral properties. As such these globins are of considerable attention. Crystal structures of few members of plant class I nonsymbiotic hemoglobin have been described earlier. Here we report the crystal structure of Arabidopsis class I hemoglobin (AHb1) to 2.2 ? and compare its key features with the structures of similar nonsymbiotic hemoglobin from other species. Crystal structure of AHb1 is homologous to the related members with similar globin fold and heme pocket architecture. The structure is homodimeric in the asymmetric unit with both distal and proximal histidines coordinating to the heme iron atom. Residues lining the dimeric interface are also conserved in AHb1 with the exception of additional electrostatic interaction between H112 and E113 of each subunit and that involving Y119 through two water molecules. In addition, differences in heme pocket non-covalent interactions, a novel Ser residue at F7 position, Xe binding site variability, internal cavity topology differences, CD loop conformation and stability and other such properties might explain kinetic variability in AHb1. Detailed cavity analysis of AHb1 showed the presence of a novel long tunnel connecting the distal pockets of both the monomers. Presence of such tunnel, along with conformational heterogeneity observed in the two chains, might suggest cooperative ligand binding and support its role in NO scavenging. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.     

Derivation of the globins from type b cytochromes

Summary Similarities in the amino acid sequences of vertebrate and invertebrate globins, b₅ and b₂ cytochromes and chicken sulfite oxidase point to a common ancestry for all of these proteins. The distal heme ligand (histidine or its equivalent) is common to both sets of proteins, but the proximal histidine ligand of the cytochromes is replaced by another histidine residue in the globins. This explains why the heme is reversed between globins and b₅ cytochromes. It seems likely that the genes for primitive globins contained three exons, the first two of which were derived from a cytochromelike DNA sequence. A model is presented to show how globins may have evolved from a pre-existing type bcytochrome; the complexity of the required changes is an indication that all globins are monophyletic.     

Spectroscopic characterization of a truncated hemoglobin from the nitrogen-fixing bacterium Herbaspirillum seropedicae

The Herbaspirillum seropedicae genome sequence encodes a truncated hemoglobin typical of group II (Hs-trHb1) members of this family. We show that His-tagged recombinant Hs-trHb1 is monomeric in solution, and its optical spectrum resembles those of previously reported globins. NMR analysis allowed us to assign heme substituents. All data suggest that Hs-trHb1 undergoes a transition from an aquomet form in the ferric state to a hexacoordinate low-spin form in the ferrous state. The close positions of Ser-E7, Lys-E10, Tyr-B10, and His-CD1 in the distal pocket place them as candidates for heme coordination and ligand regulation. Peroxide degradation kinetics suggests an easy access to the heme pocket, as the protein offered no protection against peroxide degradation when compared with free heme. The high solvent exposure of the heme may be due to the presence of a flexible loop in the access pocket, as suggested by a structural model obtained by using homologous globins as templates. The truncated hemoglobin described here has unique features among truncated hemoglobins and may function in the facilitation of O₂ transfer and scavenging, playing an important role in the nitrogen-fixation mechanism. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Hemoglobins,XLVIIII

Summary Hagfish hemoglobin has three main components, one of which is Hb III. It is monomeric and consists of 148 amino acid residues (M = 17 350). Its complete primary structure, previously published, is discussed here. The proximal amino acid (F8) of the heme linkage is histidine as always in the hemoglobins, but the regularly expected distal histidine E7 is substituted by glutamine. This substitution, leading to a new kind of heme linkage, has hitherto only been demonstrated in opossum hemoglobin. It is suggested that E7, Gln, is directed out of the heme pocket, and that the adjacent Ell, Ile, is directed toward the inside of the pocket, giving the distal heme contact instead of histidine.Myxine Hb III has an additional tail of 9 amino acid residues at its N-terminal end, as has the hemoglobin ofLampetra fluviatilis. The genetic codes ofMyxine andLampetra hemoglobins show 117 differences, in spite of many morphological resemblances between hagfish and lamprey. Their primary hemoglobin structures show differences substantial enough to bo compatible with the divergence of the two families some 400–500 million years ago.     

The polymeric hemoglobin molecule of Artemia. Interpretation of translated cDNA sequence of nine domains.

Translated cDNA for Artemia hemoglobin provided sequence data for almost nine domains, from the fourth residue of the A helix of one domain through 1405 residues to a stop codon after the ninth domain. The domain sequences were all different (homology between pairs 17-38%) but aligned well with each other and with conventional globins, satisfying the requirements for Phe at CD1, His at F8 and most other highly conserved features of globins including His at E7. Features found to be characteristic of Artemia globin and present in all nine domains were Phe at B10, Tyr at C4, Gly at F5, Phe at G5 and Gly at H22. Approximately 14 residues including a consensus -Val-Asp-Pro-Val-Thr-Gly-Leu- were available to form the linker between each pair of domains. The Artemia sequence data were compared with the crystal structures of Chironomus thummi thummi erythrocruorin III and sperm whale myoglobin in order to identify features of structural similarity and to examine the consequences of the differences. The Artemia sequences were compatible with the main helices and critical features of the globin fold. Possible modifications to the C helix, FG turn, and GH turn were studied in terms of molecular coordinates.     

Protein fold and structure in the truncated (2/2) globin family

Analysis of amino acids sequences and protein folds has recently unraveled the structural bases and details of several proteins from the recently discovered "truncated hemoglobin" family. The analysis here presented, in agreement with previous surveys, shows that truncated hemoglobins can be classified in three main groups, based on their structural properties. Crystallographic analyses have shown that all three groups adopt a 2-on-2 alpha-helical sandwich fold, resulting from apparent editing of the classical 3-on-3 alpha-helical sandwich of vertebrate and invertebrate conventional globins. Specific structural features distinguish each of the three groups. Among these, a protein matrix tunnel system is typical of group I, a Trp residue at the G8 topological site is conserved in groups II and III, and TyrB10 is almost invariant through the three groups. A strongly intertwined network of hydrogen bonds stabilizes the heme bound ligand, despite variability of the heme distal residues observed in the different proteins considered. Details of ligand recognition in the three groups are discussed at the light of residue conservation and of differing ligand diffusion pathways to the heme. Based on structural analyses of the family-specific fold, we endorse a recent proposal of leaving the "truncated hemoglobins" term, that does not represent properly the observed 2-on-2 alpha-helical sandwich fold, and adopting the simple "2/2Hb" term to concisely address this protein family.     

Self-incompatibility in Petunia inflata: isolation and characterization of cDNAs encoding three S-allele-associated proteins

Summary We have identified three alleles of the S-locus controlling self-incompatibility and their associated pistil proteins in Petunia inflata, a species that displays monofactorial gametophytic self-incompatibility. These S-allele-associated proteins (S-proteins) are pistil specific, and their levels are developmentally regulated. The amino-terminal sequences determined for the three S-proteins are highly conserved and show considerable homology to those of S-proteins from Petunia hybrida, Nicotiana alata and Lycopersicon peruvianum, three other species of the Solanaceae that also exhibit gametophytic self-incompatibility. cDNA clones encoding the three S-proteins were isolated and sequenced. Comparison of their deduced amino acid sequences reveals an average homology of 75.6%, with conserved and variable residue interspersed throughout the protein. Of the 137 conserved residues, 53 are also conserved in the N. alata S-proteins studies so far; of the 64 variable residues, 29 were identified as hypervariable based on calculation of the Similarity Index. There is only one hypervariable region of significant length, and it consists of eight consecutive hypervariable residues. This region correspond approximately to the hypervariable region HV2 identified in N. alata S-proteins. Of the two classes of N. alata S-proteins previously identified, one class exhibits greater homology to the three P. inflata S-proteins reported here than to the other class of N. alata S-proteins.     

Amino acid sequence of the dimeric hemoglobin (Hb I) from the deep-sea cold-seep clam Calyptogena soyoae and the phylogenetic relationship with other molluscan globins

The deep-sea cold-seep clam Calyptogena soyoae has two homodimeric hemoglobins (Hbs I and II) in erythrocytes. The complete amino acid sequence of Hb I has been determined. It is composed of 144 amino acid residues, has a high content of hydrophobic residues, and a calculated molecular weight of 16,350 including a heme group. The sequence of Calyptogena Hb I showed high homology (42% identity) with that of Calyptogena Hb II (Suzuki, T., Takagi T. and Ohta, S. (1989) Biochem. J. 260, 177-182), although it has a long insertion of seven residues in the C-terminal region compared with Hb II. On the other hand, it showed low homology (12-20% identity) with other molluscan globins. As well as Hb II, Calyptogena Hb I lacked the N-terminal extension of 7-9 residues characteristic of molluscan intracellular hemoglobins, and the distal (E7) histidine was replaced by glutamine. A phylogenetic tree was constructed from 13 molluscan globins belonging to the five families Aplysiidae, Galeodidae, Potamididae, Arcidae and Vesicomyidae. The globin sequences of Calyptogena (Vesicomyidae) were found to be rather distant from other globin sequences, suggesting that they might conserve a primitive form of molluscan globins.     

Analysis of conserved glutamate residues in Porphyromonas gingivalis outer membrane receptor HmuR: toward a further understanding of heme uptake

The aim of this study was to broaden the current knowledge about the Porphyromonas gingivalis heme receptor HmuR. Site-directed mutagenesis was employed to replace Glu427, Glu448, Glu458 and Glu503 by alanines and to construct a triple Glu427Ala/Glu448Ala/Glu 458Ala mutant. All iron/heme-starved P. gingivalis mutants showed decreased growth recovery when human serum as the iron/heme source was used, hmuR::ermF, hmuR ^E503A and hmuR ^{E427A,E448A,E458A} mutant strains being the most affected. E. coli cells expressing HmuR with mutated glutamate residues bound hemin, hemoglobin and hemin–serum albumin complex with the same efficiency as did the wild-type recombinant protein, suggesting that the residues were not directly involved in heme binding. These data indicate that in addition to two conserved histidine residues (His95 and His434), NPDL and YRAP motifs, conserved glutamate residues are important for HmuR to utilize heme present in serum hemoproteins.     

The cDNA-derived amino acid sequence of indoleamine dioxygenase like-myoglobin from the gastropod molluscOmphalius pfeifferi

Myoglobin was isolated from the radular muscle of the archaegastropod molluscOmphalius pfeifferi (Trochidae). The molecular mass was estimated by SDS-PAGE to be about 40 kDa, 2.5 times larger than that of usual myoglobin. The cDNA forOmphalius myoglobin was amplified by polymerase chain reaction, and the cDNA-derived amino acid sequence of 375 residues was determined, of which 73 residues were identified directly by the chemical sequencing of internal peptides. The amino acid sequence ofOmphalius myoglobin showed no significant homology with any other usual 16-kDa globins, but showed 84% and 36% identities with indoleamine dioxygenase-like myoglobins fromBattilus (Turbinidae) andSulculus (Haliotiidae), respectively. It also shows significant homology (26% identity) with human indoleamine 2,3-dioxygenase, a tryptophan-degrading enzyme containing heme. The distribution of indoleamine dioxygenase-like myoglobins suggests that they must have arisen exclusively along the specified lineage including the three families Haliotiidae, Turbinidae, and Trochidae of Archaegastropoda in molluscan evolution.     

The cDNA-derived amino acid sequence of indoleamine dioxygenase like-myoglobin from the gastropod molluscOmphalius pfeifferi

Myoglobin was isolated from the radular muscle of the archaegastropod molluscOmphalius pfeifferi (Trochidae). The molecular mass was estimated by SDS-PAGE to be about 40 kDa, 2.5 times larger than that of usual myoglobin. The cDNA forOmphalius myoglobin was amplified by polymerase chain reaction, and the cDNA-derived amino acid sequence of 375 residues was determined, of which 73 residues were identified directly by the chemical sequencing of internal peptides. The amino acid sequence ofOmphalius myoglobin showed no significant homology with any other usual 16-kDa globins, but showed 84% and 36% identities with indoleamine dioxygenase-like myoglobins fromBattilus (Turbinidae) andSulculus (Haliotiidae), respectively. It also shows significant homology (26% identity) with human indoleamine 2,3-dioxygenase, a tryptophan-degrading enzyme containing heme. The distribution of indoleamine dioxygenase-like myoglobins suggests that they must have arisen exclusively along the specified lineage including the three families Haliotiidae, Turbinidae, and Trochidae of Archaegastropoda in molluscan evolution.     

The complete amino acid sequences of four globins from the land leech Haemadipsa zeylanica var. japonica

The amino acid sequences of four globins from the land leech, Haemadipsa zeylanica var. japonica, were determined using nucleotide sequencing and protein sequencing. The mature globin-molecules were composed of 146 amino acid residues for M-1 globin, 156 for M-2 globin, 143 for D-1 globin, and 149 for D-2 globin. Alignment of the four kinds of globins by Clustal X revealed 22 invariant amino acids. The four globins were 26–33% identical. A striking feature of amino acid alteration was: the replacement of the E7 distal-His of D-1 globin by phenylalanine because histidine is conserved among the rest of the globins of H. zeylanica, those of other representative species (Lumbricus and Tylorrhynchus) of Annelida and most other hemoglobins. A phylogenetic tree constructed of 18 globin structures including two species of leeches, H. zeylanica (a land leech) and Macrobdella decora (a freshwater leech), T. heterochaetus (a representative species of polychaetes), L. terrestris (a representative species of oligochaetes), and human α and β globins strongly indicated that the leech globins first separated from globin lineage of annelids.