Globin and linker sequences of the giant extracellular hemoglobin from the leech Macrobdella decora

A detailed electrospray ionization mass spectrometric study of the 3.5-MDa hexagonal bilayer hemoglobin (HBL Hb) from the pond leech Macrobdella decora has shown it to consist of at least six 17-kDa globin chains, of which two are monomeric and the remaining four occur as disulfide-bonded heterodimers, and three 24-kDa nonglobin linker chains (Weber et al., J. Mol. Biol. 251: 703–720, 1995). The cDNA sequences of the five major constituent chains, globin chains IIA, IIB, B, and C and linker chain L1, are reported here. The globins and linkers share 30%–50% and 20%–30% identity, respectively, with other annelid sequences. Furthermore, IIB and C align with strain A of annelid sequences, whereas IIA and B align with the strain B sequences. Although chains B and C are monomeric, chains IIA and IIB form the main disulfide-bonded dimer. They also have some unusual features: the distal His (E7) is replaced by Phe in IIA, and the highly conserved CD1Phe is replaced by Leu in IIB. In spite of these unusual features, the functional properties of Macrobdella Hb are comparable to those of other HBL Hbs. A phylogenetic analysis of the globin sequences from Macrobdella, the polychaete Tylorrhynchus, the oligochaete Lumbricus, and the vestimentiferan Lamellibrachia, indicates that the two strains originated by gene duplication followed by additional duplication of each of the two strains. The mutation rate of the linkers appeared to be faster than that of the globin chains. The phylogenetic trees constructed using the Maximum Likelihood, Neighbor-Joining and Fitch methods showed the Macrobdella globin sequences to be closest to Lumbricus, in agreement with a view of annelid evolution in which the divergence of the polychaetes occurred before the divergence of the leeches from oligochaetes.