Globin gene family evolution and functional diversification in annelids

Globins are the most common type of oxygen-binding protein in annelids. In this paper, we show that circulating intracellular globin (Alvinella pompejana and Glycera dibranchiata), noncirculating intracellular globin (Arenicola marina myoglobin) and extracellular globin from various annelids share a similar gene structure, with two conserved introns at canonical positions B12.2 and G7.0. Despite sequence divergence between intracellular and extracellular globins, these data strongly suggest that these three globin types are derived from a common ancestral globin-like gene and evolved by duplication events leading to diversification of globin types and derived functions. A phylogenetic analysis shows a distinct evolutionary history of annelid extracellular hemoglobins with respect to intracellular annelid hemoglobins and mollusc and arthropod extracellular hemoglobins. In addition, dehaloperoxidase (DHP) from the annelid, Amphitrite ornata, surprisingly exhibits close phylogenetic relationships to some annelid intracellular globins. We have characterized the gene structure of A. ornata DHP to confirm assumptions about its homology with globins. It appears that it has the same intron position as in globin genes, suggesting a common ancestry with globins. In A. ornata, DHP may be a derived globin with an unusual enzymatic function.     

The multigenic family of the extracellular hemoglobin from the annelid polychaete Arenicola marina

The extracellular hemoglobin of the lugworm Arenicola marina which inhabits on the intertidal area, a sulfide-rich environment, comprises eight globin chains previously determined by mass spectrometry. We have cloned and sequenced five of the globin components. The deduced amino-acid sequences exhibit an extracellular signal peptide and two cysteine residues involved in an internal disulfide bond. The molecular weights calculated from the globin primary structures obtained from complete cDNA sequences are in good agreement with the mass spectrometry values obtained with the native hemoglobin. Phylogenetic analysis has allowed assigning the five A. marina sequences to the different globin sub-families. Two of the globins were found to be A2 globin chains lacking the cysteine residues proposed to be involved in the binding of hydrogen sulfide by such hemoglobin. We discuss the unusual absence of these cysteines in the light of their invariant occurrence in the A2 subfamily of hemoglobins from annelids inhabiting sulfide-rich environments.     

The sulfide binding function of annelid hemoglobins: relic of an old biosystem?

Invertebrates living in sulfide-rich environments have developed different strategies of coping with sulfide toxicity. Some bivalves and annelids have hemoglobins that are capable of binding sulfide for detoxification and/or transporting it to internal bacterial symbionts. Annelids living in the sulfide-rich environments have giant (approximately 3.6 MDa) hemoglobin, consisting of 144 globin chains arranged in a hexagonal bilayer structure held together by 36 nonglobin linker chains. Some globin chains contain either a free cysteine residue at positions Cys+1 or at position Cys+11 relative to the E7 distal residue in the E helix and EF interhelical region, respectively, which bind sulfide. The hexagonal bilayer hemoglobins of annelids living in environments lacking sulfide, do not have the corresponding free cysteine residues and cannot bind sulphide. Given that the early stages of life occurred under anoxic conditions in the presence of sulfide, it is possible that the sulfide binding function from modern annelid globins inhabiting sulphide rich habitats is an evolutionary relic. This proposal seems supported by the recent finding of "protoglobins" which also have a corresponding cysteine residue in Archea known to exist in hyperthermophilic and sulfide-rich environments.     

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.     

Amino acid sequence of the monomer subunit of the giant extracellular hemoglobin of the aquatic oligochaete, Tubifex tubifex

The extracellular hemoglobin of the aquatic oligochaete Tubifex tubifex consists of four subunits: a monomer of 16.5 kDa, a disulfide-bonded trimer of about 50 kDa and at least two subunits of about 30 kDa. The complete amino acid sequence of the monomeric subunit was determined: it consists of 141 amino acid residues and has a molecular mass of 16,286 Da including a heme group. 39 residues (28%) were found to be identical with those in the corresponding positions in the monomeric globin chains from Lumbricus terrestris, Pheretima sieboldi, and Tylorrhynchus heterochaetus. Tubifex and Lumbricus are most similar, with 75 amino acid identities (53%). There are eight invariant residues amongst these monomeric globins and the intracellular monomeric globin of Glycera and the human beta-globin. The monomeric globin from Tubifex aligns best with those of group A, globins which have a Cys in their second position and an invariant Lys-Val-Lys at positions 9-11 [Gotoh et al. (1987) Biochem. J. 241, 441-445]. The two cysteine residues, at positions 2 and 131, appear to be disulfide-bonded.     

Gene Structure and Molecular Phylogeny of the Linker Chains from the Giant Annelid Hexagonal Bilayer Hemoglobins

Giant extracellular hexagonal bilayer hemoglobin (HBL-Hb), found only in annelids, is an ∼3500-kDa heteropolymeric structure involved in oxygen transport. The HBL-Hbs are comprised of globin and linker chains, the latter being required for the assembly of the quaternary structure. The linker chains, varying in size from 225 to 283 amino acids, have a conserved cysteine-rich domain within their N-terminal moiety that is homologous to the cysteine-rich modules constituting the ligand binding domain of the low-density lipoprotein receptor (LDLR) protein family found in many metazoans. We have investigated the gene structure of linkers from Arenicola marina, Alvinella pompejana, Nereis diversicolor, Lumbricus terrestris, and Riftia pachyptila. We found, contrary to the results obtained earlier with linker genes from N. diversicolor and L. terrestris, that in all of the foregoing cases, the linker LDL-A module is flanked by two phase 1 introns, as in the human LDLR gene, with two more introns in the 3′ side whose positions varied with the species. In addition, we obtained 13 linker cDNAs that have been determined experimentally or found in the EST database LumbriBASE. A molecular phylogenetic analysis of the linker primary sequences demonstrated that they cluster into two distinct families of linker proteins. We propose that the common gene ancestor to annelid linker genes exhibited a four-intron and five-exon structure and gave rise to the two families subsequent to a duplication event.     

The cDNA sequences encoding two components of the polymeric fraction of the intracellular hemoglobin of Glycera dibranchiata

The intracellular hemoglobin of the polychaete Glycera dibranchiata consists of several components, some of which self-associate into a "polymeric" fraction. The cDNA library constructed from the poly(A+) mRNA of Glycera erythrocytes (Simons, P. C., and Satterlee, J. D. (1989) Biochemistry 28, 8525-8530) was screened with two oligodeoxynucleotide probes corresponding to the amino acid sequences MEEKVP and AMNSKV. Each of the two probes identified a full-length positive insert; these were sequenced using the dideoxynucleotide chain termination method. One clone was 630 bases long and contained 36 bases of 5'-untranslated RNA, a reading frame of 441 bases coding for the 147 amino acids of globin P2 including the residues MEEKVP, and a 3'-untranslated region of 153 bases. The other clone was 540 bases long and contained 24 bases of 5'-untranslated RNA, an open reading frame of 441 bases coding for globin P3 including the residues AMNSKV, and a 3'-untranslated region of 75 bases. The inferred amino acid sequences of the two globins were in agreement with the partial amino acid sequences obtained by chemical methods. The P2 and P3 globin sequences, together with the previously determined P1 sequence of a complete insert and partial sequences P4, P5, and P6 obtained from partial inserts (Zafar, R. S., Chow, L. H., Stern, M. S., Vinogradov, S. N., and Walz, D. A. (1990) Biochim. Biophys. Acta, in press) suggest that there are at least six components in the polymeric fraction of Glycera hemoglobin, which is in agreement with the results of polyacrylamide gel electrophoresis in Tris/glycine buffer, pH 8.3, 6 M urea. Nothern and dot blot analyses of Glycera erythrocyte poly(A+) mRNA using the foregoing two cDNA probes clearly demonstrated the presence of mature messages encoding both types of globins. Comparison of the polymeric sequences P1, P2, and P3 with the "monomeric" globins M-II and M-IV using the alignment and templates of Bashford et al. (Bashford, D., Chothia, C., and Lesk, A. M. (1987) J. Mol. Biol. 196, 199-216) showed that all five globins have identical residues at 39 positions. At 44 positions, the three polymeric globins share identical residues that differ from the identical residues at the corresponding locations in the monomeric sequences M-II and M-IV including position E7, where the latter have leucine instead of the distal histidine. At 15 positions, there occurs an alteration from polar to nonpolar or from a small nonpolar to a larger nonpolar residue in going from the monomeric to the polymeric globins.(ABSTRACT TRUNCATED AT 400 WORDS)     

The evolution of extracellular hemoglobins of annelids, vestimentiferans, and pogonophorans

The evolution of extracellular hemoglobins of annelids, vestimentiferans, and pogonophorans was investigated by applying cladistic and distance-based approaches to reconstruct the phylogenetic relationships of this group of respiratory pigments. We performed this study using the aligned sequences of globin and linker chains that are the constituents of these complex molecules. Three novel globin and two novel linker chains of Sabella spallanzanii described in an accompanying paper (Pallavicini, A., Negrisolo, E., Barbato, R., Dewilde, S., Ghiretti-Magaldi, A., Moens, L., and Lanfranchi, G. (2001) J. Biol. Chem. 276, 26384--26390) were also included. Our results allowed us to test previous hypotheses on the evolutionary pathways of these proteins and to formulate a new most parsimonious model of molecular evolution. According to this novel model, the genes coding for the polypeptides forming these composite molecules were already present in the common ancestor of annelids, vestimentiferans, and pogonophorans.     

Genetic variation in mouse beta globin cysteine content modifies glutathione metabolism: implications for the use of mouse models

Allelic variation in the mouse beta globin gene complex (Hbb) produces structurally different beta globins in different mouse strains. Like humans, mice with HbbS alleles produce a single beta globin with one reactive cysteine (beta Cys93). In contrast, mice with HbbD alleles produce two structurally different beta globins, each containing an additional cysteine (beta Cys13). beta Cys93 forms mixed disulfides with glutathione and plays a pivotal role in the activities of hemoglobin, glutathione, and nitric oxide. Similar roles for mouse beta Cys13 have not been described. We used capillary electrophoresis to compare reduced glutathione (GSH), glutathione disulfide (GSSG), and S-glutathionyl hemoglobin levels in erythrocytes from inbred C57BL/6J (homozygous HbbS/S) and 129S1/SvImJ (homozygous HbbD/D) mice and their homozygous and heterozygous B6129S/F2J hybrid offspring. S-glutathionyl hemoglobin was nearly undetectable in inbred or hybrid mice with only monocysteinyl beta globins (HbbS/S) but represented up to 10% of total hemoglobin in mice with polycysteinyl beta globins (HbbS/D or HbbD/D). The stepwise increase in beta globin sulfhydryl group concentration in HbbS/S, HbbS/D, and HbbD/D F2 mice was associated with increasing hemoglobin-bound glutathione and decreasing free glutathione (GSH + GSSG) concentrations. Total erythrocyte glutathione (GSH + GSSG + hemoglobin-bound) was not significantly different between groups. In vitro studies showed that beta Cys13 in mouse HbbD beta globins was more susceptible to disulfide exchange with GSSG than beta Cys93. We conclude that reactive beta globin sulfhydryl group concentration is genetically determined in mice, and that polycysteinyl beta globins markedly influence intraerythrocyte glutathione distribution between free and hemoglobin-bound compartments. Although Hbb heterozygosity and polycysteinyl beta globins are common in wild mouse populations, all common human beta globins contain only one reactive cysteine, and homozygosity is the norm. These fundamental differences in mouse and human beta globin genetics have important implications for the study of mouse biology and for the use of some mouse strains as models for humans.     

Structural characterization of hemoglobins from Monilifera and Frenulata tubeworms (Siboglinids): First discovery of giant hexagonal-bilayer hemoglobin in the former “Pogonophora” group

Siboglinids are symbiotic polychete annelids having hemoglobins as essential oxygen- and sulfide-carriers for their endosymbiotic bacteria. We analyzed the structure of the hemoglobins from two species of siboglinids: the monilifera Sclerolinum contortum and the frenulata Oligobrachia webbi (i.e. haakonmosbiensis) from Norwegian cold seeps. Measured by Multi-Angle Laser Light Scattering (MALLS), Sclerolinum shows a 3190 ± 50 kDa hexagonal bilayer hemoglobin (HBL-Hb) and a 461 ± 46 kDa ring-Hb, just as vestimentifera, whereas Oligobrachia has a 409 ± 3.7 kDa ring-Hb only. Electrospray Ionization-Mass Spectrometry (ESI-MS) showed Sclerolinum HBL-Hb composed of seven monomeric globins (15–16 kDa), three disulfide-bonded globin heterodimers and three linkers. The heterodimers always contain globin-b (15814.4 ± 1.5 Da). Sclerolinum ring-Hb is composed of globins and dimers with identical masses as its HBL-Hb, but lacks linkers. Oligobrachia ring-Hb has three globin monomers (14–15 kDa) only, with no disulfide-bonded dimers. Comparison of Sclerolinum hemoglobins between Storegga and Haakon Mosby Mud Volcano, using the normalized height of deconvoluted ESI-MS peaks, shows differences in globin monomers abundances that could reflect genetic differences or differential gene expression between distinct seep populations. The discovery of HBL-Hb in Sclerolinum is a new element supporting the hypothesis of monilifera being phylogenetically more closely related to vestimentifera, than to frenulata.     

Variation of folded polypeptide surface area with probe size

Three types of polypeptide surface area (contact, accessible, and molecular) have been studied as a function of the radius of a probe sphere used to map the surface. The surfaces are: (1) three alpha-helices, the H-helix of myoglobin, the E-helix of leghemoglobin, and an artificial polyalanine helix, each with 26 residues; (2) two globins, myoglobin and leghemoglobin, each with 153 residues; and (3) a two-center model system for which the three types of surface area have been calculated analytically. The two globin helices have almost identical surface areas as a function of probe size as do the two globins. The polyalanine helix surface area is smaller but similar in shape to the globin helix areas. All three helix contact areas tend to the same limit as the probe size increases, and the globin contact areas behave similarly. Fractal dimensions were calculated for the helix and globin contact and molecular surfaces. All fractal dimensions showed strong dependence on probe size. The contact fractal dimension peaks at larger values for both the helices and globins. Most residues do not make contact with large probes (15 A).     

Amino acid sequence of a major globin from the sea cucumber Paracaudina chilensis

The sea cucumber Paracaudina chilensis (Echinodermata) contains three major globins I, II and III in coelomic cells. The complete amino acid sequence of globin I has been determined. It is composed of 157 amino acid residues, is acetylated at the N-terminus, and has a characteristic N-terminal extension of 9-10 residues when compared with vertebrate globins. The sequence of Paracaudina globin I showed slightly higher homology with human alpha globin (25%) rather than with the invertebrate Anadara alpha globin (22%). Paracaudina globin I also showed strong homology (59%) with globin D from another sea cucumber, Molpadia arenicola (Mauri, F.C. (1985) Ph.D. dissertation, University of Texas). The globin sequences from the phylum Echinodermata have an important position in the molecular evolution of the globins, because they are the invertebrate group most closely related to the vertebrates.     

Non-functional conserved residues in globins and their possible role as a folding nucleus.

Structure-based sequence alignment of 728 sequences of different globin subfamilies shows that in each subfamily there are two clusters of consensually conserved residues. The first is the well-known "functional" cluster which includes six heme-binding conserved residues (Phe CD1, His F8; aliphatic E11, FG5; hydrophobic F4, G5) and seven other conserved residues (Pro C2; aliphatic H19; hydrophobic B10, B13, B14, CD4, E4) that do not bind the heme but belong to its immediate neighborhood. The second cluster revealed here (aliphatic A8, G16, G12; aromatic A12; hydrophobic H8 and possibly H12) is distant from the heme. It is entirely non-polar and includes one turn (i, i+4 positions) from each of helices A, G, and H. It is known that A, G, and H helices formed at the earliest stage of apomyoglobin folding remain relatively stable in the equilibrium molten globule state, and are likely to be tightly packed with each other in this state. We have shown the existence of two similar conserved clusters in c -type cytochromes, heme-binding and distal from the heme. The second cluster in c -cytochromes includes one turn from each of the N and C-terminal alpha-helices. These N and C-terminal helices in cytochrome c are formed at the earliest stage of protein folding, remain relatively stable in the molten globule state, and are tightly packed with each other in this state, similar to the observed behavior of the globins. At least these two large protein families (c -type cytochromes and globins) have a close similarity in the existence and mutual positions of non-functional conserved residues. We assume that non-functional conserved residues are requisite for the fast and correct folding of both of these protein families into their stable 3D structures.     

A model of globin evolution

Putative globins have been identified in 426 bacterial, 32 Archaeal and 67 eukaryote genomes. Among these sequences are the hitherto unsuspected presence of single domain sensor globins within Bacteria, Fungi, and a Euryarchaeote. Bayesian phylogenetic trees suggest that their occurrence in the latter two groups could be the result of lateral gene transfer from Bacteria. Iterated psiblast searches based on groups of globin sequences indicate that bacterial flavohemoglobins are closer to metazoan globins than to the other two lineages, the 2-over-2 globins and the globin-coupled sensors. Since Bacteria is the only kingdom to have all the subgroups of the three globin lineages, we propose a working model of globin evolution based on the assumption that all three lineages originated and evolved only in Bacteria. Although the 2-over-2 globins and the globin-coupled sensors recognize flavohemoglobins, there is little recognition between them. Thus, in the first stage of globin evolution, we favor a flavohemoglobin-like single domain protein as the ancestral globin. The next stage comprised the splitting off to single domain 2-over-2 and sensor-like globins, followed by the covalent addition of C-terminal domains resulting in the chimeric flavohemoglobins and globin-coupled sensors. The last stage encompassed the lateral gene transfers of some members of the three globin lineages to specific groups of Archaea and Eukaryotes.     

Protein and gene structure of a chlorocruorin chain of Eudistylia vancouverii

The polychaete annelid, Eudistylia vancouverii, contains as oxygen carrier a hexagonal bilayer (HBL) chlorocruorin. One of the globin chains, chain a1, has 142 amino acids (Mr 16,054.99) and its sequence deviates strongly from other nonvertebrate globin sequences. Unprecedented, it displays a Phe at the distal position E7 as well as at position B10, creating a very hydrophobic heme pocket probably responsible for the low oxygen affinity of the native molecule. Phylogenetic analysis of annelid globin chains clearly proves that globin chain a1 belongs to type I of globin chains having a pattern of 3 cysteine residues essential for the aggregation into a HBL structure. The gene coding for globin chain a1 is interrupted by 2 introns at the conserved positions B12.2 and G7.0. Based on protein and gene structure it can therefore be concluded that the globin chains of chlorocruorins are not fundamentally different from other annelid globin chains.     

Gene duplication,genome duplication,and the functional diversification of vertebrate globins

The functional diversification of the vertebrate globin gene superfamily provides an especially vivid illustration of the role of gene duplication and whole-genome duplication in promoting evolutionary innovation. For example, key globin proteins that evolved specialized functions in various aspects of oxidative metabolism and oxygen signaling pathways (hemoglobin [Hb], myoglobin [Mb], and cytoglobin [Cygb]) trace their origins to two whole-genome duplication events in the stem lineage of vertebrates. The retention of the proto-Hb and Mb genes in the ancestor of jawed vertebrates permitted a physiological division of labor between the oxygen-carrier function of Hb and the oxygen-storage function of Mb. In the Hb gene lineage, a subsequent tandem gene duplication gave rise to the proto α- and β-globin genes, which permitted the formation of multimeric Hbs composed of unlike subunits (α₂β₂). The evolution of this heteromeric quaternary structure was central to the emergence of Hb as a specialized oxygen-transport protein because it provided a mechanism for cooperative oxygen-binding and allosteric regulatory control. Subsequent rounds of duplication and divergence have produced diverse repertoires of α- and β-like globin genes that are ontogenetically regulated such that functionally distinct Hb isoforms are expressed during different stages of prenatal development and postnatal life. In the ancestor of jawless fishes, the proto Mb and Hb genes appear to have been secondarily lost, and the Cygb homolog evolved a specialized respiratory function in blood-oxygen transport. Phylogenetic and comparative genomic analyses of the vertebrate globin gene superfamily have revealed numerous instances in which paralogous globins have convergently evolved similar expression patterns and/or similar functional specializations in different organismal lineages.     

Globins in nonvertebrate species: dispersal by horizontal gene transfer and evolution of the structure-function relationships

Using a new template based on an alignment of 145 nonvertebrate globins we examined several recently determined sequences of putative globins and globin-like hemeproteins. We propose that all globins have evolved from a family of ancestral, approx. 17-kDa hemeproteins, which displayed the globin fold and functioned as redox proteins. Once atmospheric O2 became available the acquisition of oxygen-binding properties was initiated, culminating in the various highly specialized functions known as present. During this evolutionary process, we suggest that (1) high oxygen affinity may have been acquired repeatedly and (2) the formation of chimeric proteins containing both a globin and a flavin binding domain was an additional and distinct evolutionary trend. Furthermore, globin-like hemeproteins encompass hemeproteins produced through convergent evolution from nonglobin ancestral proteins to carry out O2-binding functions as well as hemeproteins whose sequences exhibit the loss of some or all of the structural determinants of the globin fold. We also propose that there occurred two cases of horizontal globin gene transfer, one from an ancestor common to the ciliates Paramecium and Tetrahymena and the green alga Chlamydomonas to a cyanobacterium ancestor and the other, from a eukaryote ancestor of the yeasts Saccharomyces and Candida to a bacterial ancestor of the proteobacterial genera Escherichia, Alcaligenes, and Vitreoscilla.      

Primary structure of the common polypeptide chain b from the multi-hemoglobin system of the hydrothermal vent tube worm Riftia pachyptila: An insight on the sulfide binding-site

The deep-sea tube worm Riftia pachyptila Jones possesses a multi-hemoglobin system with three different extracellular Hbs: two dissolved in the vascular blood, V1 (ca. 3,500 kDa) and V2 (ca. 400 kDa), and one in the coelomic fluid, C1 (ca. 400 kDa). V1 Hb consists of four heme-containing, globin chains (b–e) and four linker chains (L1–L4). V2 and C1 Hbs are exclusively built from globin chains, six for V2 (a–f) and five for C1 (a–e). The complete amino acid sequence of the isolated monomeric globin chain b, common to all Riftia Hbs, has been determined by automated Edman degradation sequencing of the peptides derived by digestion with trypsin, chymotrypsin, thermolysin, and CNBr. This polypeptide chain is composed of 144 amino acid residues, providing a Mr of 16, 135.0 Da. Moreover, the primary sequence of chain b revealed 3 Cys residues at position 4, 75, and 134. Cys-4 and Cys-134 are located at positions where an intra-chain disulfide bridge is formed in all annelid, vestimentiferan, or pogonophoran chains, but Cys-75 is located at a unique position only found in three globin chains belonging to Lamellibrachia and Oligobrachia, a vestimentiferan and a pogonophoran. In both groups, Hbs can bind sulfide reversibly to fuel the chemosynthetic process of the symbiotic bacteria they harbor. Sulfide-binding experiments performed on purified Hb fractions (i.e., V1, V2, and C1 Hbs) suggest that free Cys residues on globin chains, and the numerous Cys found in linker chains, as determined previously by ESI-MS, may be the sulfide binding-sites. Blocking the free Cys by N-ethylmaleimide, we confirmed that free cysteines were involved in sulfide-binding but did not account for the whole sulfide-binding capacity of V1 Hb. Furthermore, a phylogenetic tree was constructed from 18 globin-like chains of annelid, vetimentiferan, and pogonophoran extracellular Hbs to clarify the systematic position of tubeworms. Riftia chain b clearly belongs to the “strain A” family with 30 to 80% identity with the other sequences analyzed. Its position in the tree confirmed a close relationship between vestimentiferan, pogonophoran, and annelid Hbs. Proteins 29:562–574, 1997. © 1997 Wiley-Liss, Inc.     

Axolotl hemoglobin: cDNA-derived amino acid sequences of two alpha globins and a beta globin from an adult Ambystoma mexicanum

Erythrocytes of the adult axolotl, Ambystoma mexicanum, have multiple hemoglobins. We separated and purified two kinds of hemoglobin, termed major hemoglobin (Hb M) and minor hemoglobin (Hb m), from a five-year-old male by hydrophobic interaction column chromatography on Alkyl Superose. The hemoglobins have two distinct alpha type globin polypeptides (alphaM and alpham) and a common beta globin polypeptide, all of which were purified in FPLC on a reversed-phase column after S-pyridylethylation. The complete amino acid sequences of the three globin chains were determined separately using nucleotide sequencing with the assistance of protein sequencing. The mature globin molecules were composed of 141 amino acid residues for alphaM globin, 143 for alpham globin and 146 for beta globin. Comparing primary structures of the five kinds of axolotl globins, including two previously established alpha type globins from the same species, with other known globins of amphibians and representatives of other vertebrates, we constructed phylogenetic trees for amphibian hemoglobins and tetrapod hemoglobins. The molecular trees indicated that alphaM, alpham, beta and the previously known alpha major globin were adult types of globins and the other known alpha globin was a larval type. The existence of two to four more globins in the axolotl erythrocyte is predicted.     

The phylogeny of human globin genes investigated by the maximum parsimony method

Gene phylogenetic trees were constructed by the maximum parsimony method for various sets of ninety six globin chain amino acid sequences spanning plant and animal kingdoms. The method, executed by several computer programs, constructed ancestor and descendant globin messengers on tree topologies which required the least number of nucleotide replacements to account for the evolution of the globins. The human myoglobin-hemoglobin divergence was traced to a gene duplication which occurred either in the first vertebrates or earlier yet in the common ancestor of chordates and annelids, the alpha-beta divergence to a gene duplication in the common ancestor of teleosts and tetrapods, the gamma divergence from typical beta chains to a gene duplication in basal therian mammals, and the delta separation from beta to a duplication in the basal catarrhine primates. Evidence was provided by the globin phylogenies for the hominoid affinities of the gibbon and the close phyletic relationship of the African apes to man. Over the period of teleos-tetrapod divergence the globin messengers evolved at an average rate of 18.5 nucleotide replacements per 100 codons per 10⁸ years, a faster rate than most previous estimates. Very fast and very slow rates were encountered in different globin lineages and at different stages of descent, reducing the effectiveness of globins as molecular clocks. Rates increased with gene duplication and decreased after selection discovered useful specializations in the products of genes which had previously been freer to accept mutations. The early eutherian radiation was characterized by rapid rates of globin evolution, but the later hominoid radiation by extremely slow rates. This pattern was related to more complicated grades of internal organization evolving in human ancestors. The types of nucleotide replacements in the globin messengers over the long course of globin evolution did not seem indicative of any special mutational mechanisms.