Bacterial and archaeal globins - a revised perspective.

A bioinformatics survey of putative globins in over 2200 bacterial and some 140 archaeal genomes revealed that over half the bacterial and approximately one fifth of archaeal genomes contain genes encoding globins that were classified into three families: the M (myoglobin-like), and S (sensor) families all exhibiting the canonical 3/3 myoglobin fold, and the T family (truncated myoglobin fold). Although the M family comprises 2 subfamilies, flavohemoglobins (FHbs) and single domain globins (SDgbs), the S family encompasses chimeric globin-coupled sensors (GCSs), single domain Pgbs (protoglobins) and SSDgbs (sensor single domain globins). The T family comprises three classes TrHb1s, TrHb2s and TrHb3s, characterized by the abbreviated 2/2 myoglobin fold. The Archaea contain only Pgbs, GCSs and TrHb1s. The smallest globin-bearing genomes are the streamlined genomes (~1.3Mbp) of the SAR11 clade of alphaproteobacteria and the slightly larger (ca. 1.7Mbp) genomes of Aquificae. The smallest genome with members of all three families is the 2.3Mbp genome of the extremophile Methylacidiphilum infernorum (Verrumicrobia). Of the 147 possible combinations of the eight globin subfamilies, only 83 are observed. Although binary combinations are infrequent and ternary combinations are rare, the FHb+TrHb2 combination is the most commonly observed. Of the possible functions of bacterial globins we discuss the two principal ones - nitric oxide detoxification via the NO dioxygenase or denitrosylase activities and the sensing of oxygen concentration in the environmental niche. In only few cases has a physiological role been demonstrated in vivo. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Evolution of the globin gene family in deuterostomes: lineage-specific patterns of diversification and attrition.

In the Metazoa, globin proteins display an underlying unity in tertiary structure that belies an extraordinary diversity in primary structures, biochemical properties, and physiological functions. Phylogenetic reconstructions can reveal which of these functions represent novel, lineage-specific innovations, and which represent ancestral functions that are shared with homologous globin proteins in other eukaryotes and even prokaryotes. To date, our understanding of globin diversity in deuterostomes has been hindered by a dearth of genomic sequence data from the Ambulacraria (echinoderms + hemichordates), the sister group of chordates, and the phylum Xenacoelomorpha, which includes xenoturbellids, acoelomorphs, and nemertodermatids. Here, we report the results of a phylogenetic and comparative genomic analysis of the globin gene repertoire of deuterostomes. We first characterized the globin genes of the acorn worm, Saccoglossus kowalevskii, a representative of the phylum Hemichordata. We then integrated genomic sequence data from the acorn worm into a comprehensive analysis of conserved synteny and phylogenetic relationships among globin genes from representatives of the eight lineages that comprise the superphylum Deuterostomia. The primary aims were 1) to unravel the evolutionary history of the globin gene superfamily in deuterostomes and 2) to use the estimated phylogeny to gain insights into the functional evolution of deuterostome globins. Results of our analyses indicate that the deuterostome common ancestor possessed a repertoire of at least four distinct globin paralogs and that different subsets of these ancestral genes have been retained in each of the descendant organismal lineages. In each major deuterostome group, a different subset of ancestral precursor genes underwent lineage-specific expansions of functional diversity through repeated rounds of gene duplication and divergence. By integrating results of the phylogenetic analysis with available functional data, we discovered that circulating oxygen-transport hemoglobins evolved independently in several deuterostome lineages and that intracellular nerve globins evolved independently in chordates and acoelomorph worms.

Androglobin: a chimeric globin in metazoans that is preferentially expressed in Mammalian testes.

Comparative genomic studies have led to the recent identification of several novel globin types in the Metazoa. They have revealed a surprising evolutionary diversity of functions beyond the familiar O(2) supply roles of hemoglobin and myoglobin. Here we report the discovery of a hitherto unrecognized family of proteins with a unique modular architecture, possessing an N-terminal calpain-like domain, an internal, circular permuted globin domain, and an IQ calmodulin-binding motif. Putative orthologs are present in the genomes of many metazoan taxa, including vertebrates. The calpain-like region is homologous to the catalytic domain II of the large subunit of human calpain-7. The globin domain satisfies the criteria of a myoglobin-like fold but is rearranged and split into two parts. The recombinantly expressed human globin domain exhibits an absorption spectrum characteristic of hexacoordination of the heme iron atom. Molecular evolutionary analyses indicate that this chimeric globin family is phylogenetically ancient and originated in the common ancestor to animals and choanoflagellates. In humans and mice, the gene is predominantly expressed in testis tissue, and we propose the name "androglobin" (Adgb). Expression is associated with postmeiotic stages of spermatogenesis and is insensitive to experimental hypoxia. Evidence exists for increased gene expression in fertile compared with infertile males.

The globin gene family of the cephalochordate amphioxus: implications for chordate globin evolution.

BACKGROUND: The lancelet amphioxus (Cephalochordata) is a close relative of vertebrates and thus may enhance our understanding of vertebrate gene and genome evolution. In this context, the globins are one of the best studied models for gene family evolution. Previous biochemical studies have demonstrated the presence of an intracellular globin in notochord tissue and myotome of amphioxus, but the corresponding gene has not yet been identified. Genomic resources of Branchiostoma floridae now facilitate the identification, experimental confirmation and molecular evolutionary analysis of its globin gene repertoire. RESULTS: We show that B. floridae harbors at least fifteen paralogous globin genes, all of which reveal evidence of gene expression. The protein sequences of twelve globins display the conserved characteristics of a functional globin fold. In phylogenetic analyses, the amphioxus globin BflGb4 forms a common clade with vertebrate neuroglobins, indicating the presence of this nerve globin in cephalochordates. Orthology is corroborated by conserved syntenic linkage of BflGb4 and flanking genes. The kinetics of ligand binding of recombinantly expressed BflGb4 reveals that this globin is hexacoordinated with a high oxygen association rate, thus strongly resembling vertebrate neuroglobin. In addition, possible amphioxus orthologs of the vertebrate globin X lineage and of the myoglobin/cytoglobin/hemoglobin lineage can be identified, including one gene as a candidate for being expressed in notochord tissue. Genomic analyses identify conserved synteny between amphioxus globin-containing regions and the vertebrate ß-globin locus, possibly arguing against a late transpositional origin of the ß-globin cluster in vertebrates. Some amphioxus globin gene structures exhibit minisatellite-like tandem duplications of intron-exon boundaries ("mirages"), which may serve to explain the creation of novel intron positions within the globin genes. CONCLUSIONS: The identification of putative orthologs of vertebrate globin variants in the B. floridae genome underlines the importance of cephalochordates for elucidating vertebrate genome evolution. The present study facilitates detailed functional studies of the amphioxus globins in order to trace conserved properties and specific adaptations of respiratory proteins at the base of chordate evolution.

Structural determinants for the formation of sulfhemeprotein complexes.

Several hemoglobins were explored by UV-Vis and resonance Raman spectroscopy to define sulfheme complex formation. Evaluation of these proteins upon the reaction with H(2)O(2) or O(2) in the presence of H(2)S suggest: (a) the formation of the sulfheme derivate requires a HisE7 residue in the heme distal site with an adequate orientation to form an active ternary complex; (b) that the ternary complex intermediate involves the HisE7, the peroxo or ferryl species, and the H(2)S molecule. This moiety precedes and triggers the sulfheme formation.

Lessons from the post-genomic era: Globin diversity beyond oxygen binding and transport.

Vertebrate hemoglobin (Hb) and myoglobin (Mb) were among the first proteins whose structures and sequences were determined over 50 years ago. In the subsequent pregenomic period, numerous related proteins came to light in plants, invertebrates and bacteria, that shared the myoglobin fold, a signature sequence motif characteristic of a 3-on-3 α-helical sandwich. Concomitantly, eukaryote and bacterial globins with a truncated 2-on-2 α-helical fold were discovered. Genomic information over the last 20 years has dramatically expanded the list of known globins, demonstrating their existence in a limited number of archaeal genomes, a majority of bacterial genomes and an overwhelming majority of eukaryote genomes. In vertebrates, 6 additional globin types were identified, namely neuroglobin (Ngb), cytoglobin (Cygb), globin E (GbE), globin X (GbX), globin Y (GbY) and androglobin (Adgb). Furthermore, functions beyond the familiar oxygen transport and storage have been discovered within the vertebrate globin family, including NO metabolism, peroxidase activity, scavenging of free radicals, and signaling functions. The extension of the knowledge on globin functions suggests that the original roles of bacterial globins must have been enzymatic, involved in defense against NO toxicity, and perhaps also as sensors of O2, regulating taxis away or towards high O2 concentrations. In this review, we aimed to discuss the evolution and remarkable functional diversity of vertebrate globins with particular focus on the variety of non-canonical expression sites of mammalian globins and their according impressive variability of atypical functions.

A phylogenomic profile of hemerythrins, the nonheme diiron binding respiratory proteins.

BACKGROUND: Hemerythrins, are the non-heme, diiron binding respiratory proteins of brachiopods, priapulids and sipunculans; they are also found in annelids and bacteria, where their functions have not been fully elucidated. RESULTS: A search for putative Hrs in the genomes of 43 archaea, 444 bacteria and 135 eukaryotes, revealed their presence in 3 archaea, 118 bacteria, several fungi, one apicomplexan, a heterolobosan, a cnidarian and several annelids. About a fourth of the Hr sequences were identified as N- or C-terminal domains of chimeric, chemotactic gene regulators. The function of the remaining single domain bacterial Hrs remains to be determined. In addition to oxygen transport, the possible functions in annelids have been proposed to include cadmium-binding, antibacterial action and immunoprotection. A Bayesian phylogenetic tree revealed a split into two clades, one encompassing archaea, bacteria and fungi, and the other comprising the remaining eukaryotes. The annelid and sipunculan Hrs share the same intron-exon structure, different from that of the cnidarian Hr. CONCLUSION: The phylogenomic profile of Hrs demonstrated a limited occurrence in bacteria and archaea and a marked absence in the vast majority of multicellular organisms. Among the metazoa, Hrs have survived in a cnidarian and in a few protostome groups; hence, it appears that in metazoans the Hr gene was lost in deuterostome ancestor(s) after the radiata/bilateria split. Signal peptide sequences in several Hirudinea Hrs suggest for the first time, the possibility of extracellular localization. Since the alpha-helical bundle is likely to have been among the earliest protein folds, Hrs represent an ancient family of iron-binding proteins, whose primary function in bacteria may have been that of an oxygen sensor, enabling aerophilic or aerophobic responses. Although Hrs evolved to function as O2 transporters in brachiopods, priapulids and sipunculans, their function in annelids remains to be elucidated. Overall Hrs exhibit a considerable lack of evolutionary success in metazoans.

The Caenorhabditis globin gene family reveals extensive nematode-specific radiation and diversification.

BACKGROUND: Globin isoforms with variant properties and functions have been found in the pseudocoel, body wall and cuticle of various nematode species and even in the eyespots of the insect-parasite Mermis nigrescens. In fact, much higher levels of complexity exist, as shown by recent whole genome analysis studies. In silico analysis of the genome of Caenorhabditis elegans revealed an unexpectedly high number of globin genes featuring a remarkable diversity in gene structure, amino acid sequence and expression profiles. RESULTS: In the present study we have analyzed whole genomic data from C. briggsae, C. remanei, Pristionchus pacificus and Brugia malayi and EST data from several other nematode species to study the evolutionary history of the nematode globin gene family. We find a high level of conservation of the C. elegans globin complement, with even distantly related nematodes harboring orthologs to many Caenorhabditis globins. Bayesian phylogenetic analysis resolves all nematode globins into two distinct globin classes. Analysis of the globin intron-exon structures suggests extensive loss of ancestral introns and gain of new positions in deep nematode ancestors, and mainly loss in the Caenorhabditis lineage. We also show that the Caenorhabditis globin genes are expressed in distinct, mostly non-overlapping, sets of cells and that they are all under strong purifying selection. CONCLUSION: Our results enable reconstruction of the evolutionary history of the globin gene family in the nematode phylum. A duplication of an ancestral globin gene occurred before the divergence of the Platyhelminthes and the Nematoda and one of the duplicated genes radiated further in the nematode phylum before the split of the Spirurina and Rhabditina and was followed by further radiation in the lineage leading to Caenorhabditis. The resulting globin genes were subject to processes of subfunctionalization and diversification leading to cell-specific expression patterns. Strong purifying selection subsequently dampened further evolution and facilitated fixation of the duplicated genes in the genome.

Tracing globin phylogeny using PSIBLAST searches based on groups of sequences.

PSIBLAST search of a protein database with a query sequence is a widely used tool for the detection of related but evolutionarily distant sequences. Iterations carried out until convergence (i.e., until the majority of the sequences most similar to the query sequence are extracted from the database) also produces an ordered list of more distantly related (i.e., less similar sequences) and false positives belonging to other protein families. Thus, a PSIBLAST search based on one group of globin sequences will provide sequences left in the database that are more distantly related (i.e., belong to other groups of globins) ordered according to the E value or bit score. The relative order of the scores should yield a clue to the relative similarity of the query group to the other groups of globins. Histograms of E values or bit scores from PSIBLAST searches using vertebrate myoglobins, cytoglobins, alpha- and beta-globins, and neuroglobins as query groups show distributions that are congruent with the accepted phylogenetic tree. An illustration of more distant relationships is demonstrated by the results using neuroglobins as the query group, which show a striking similarity to bacterial single-domain globins and flavohemoglobins from bacteria and eukaryotes. Furthermore, it is observed that sequences belonging to the three undoubtedly ancient globin lineages form very broad distributions, while recently evolved groups such as cytoglobins have narrow distributions. Thus, the breadth of a distribution of E values or bit scores for a query group may be related to its evolutionary age.

Mass mapping of large globin complexes by scanning transmission electron microscopy.

Scanning transmission electron microscopy (STEM) of unstained, freeze-dried biological macromolecules in the dark-field mode provides an image based on the number of electrons elastically scattered by the constituent atoms of the macromolecule. The image of each isolated particle provides information about the projected structure of the latter, and its integrated intensity is directly related to the mass of the selected particle. Particle images can be sorted by shape, providing independent histograms of mass to study assembly/disassembly intermediates. STEM is optimized for low-dose imaging and is suitable for accurate measurement of particle masses over the range from about 30 kDa to 1,000 MDa. This article describes the details of the method developed at the Brookhaven National Laboratory STEM facility and illustrates its application to the mass mapping of large globin complexes.

Expression and in silico structural analysis of a rice (Oryza sativa) hemoglobin 5.

This work reports the analysis of an additional hemoglobin (hb) gene copy, hb5, in the genome of rice. The amino acid sequence of Hb5 differs from the previously determined rice Hbs 1-4 in missing 11 residues in helix E. Transcripts of hb5 were found to be ubiquitous in rice organs, and hormone- and stress-response promoters exist upstream of the rice hb5 gene. Furthermore, the modeled structure of Hb5 based on the known crystal structure of rice Hb1 is unusual in that the putative distal His is distant from the heme Fe. This observation suggests that Hb5 binds and releases O(2) easily and thus that it functions as an O(2)-carrier or in some aspects of the O(2) metabolism.

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.

A phylogenomic profile of globins.

BACKGROUND: Globins occur in all three kingdoms of life: they can be classified into single-domain globins and chimeric globins. The latter comprise the flavohemoglobins with a C-terminal FAD-binding domain and the gene-regulating globin coupled sensors, with variable C-terminal domains. The single-domain globins encompass sequences related to chimeric globins and "truncated" hemoglobins with a 2-over-2 instead of the canonical 3-over-3 alpha-helical fold. RESULTS: A census of globins in 26 archaeal, 245 bacterial and 49 eukaryote genomes was carried out. Only approximately 25% of archaea have globins, including globin coupled sensors, related single domain globins and 2-over-2 globins. From one to seven globins per genome were found in approximately 65% of the bacterial genomes: the presence and number of globins are positively correlated with genome size. Globins appear to be mostly absent in Bacteroidetes/Chlorobi, Chlamydia, Lactobacillales, Mollicutes, Rickettsiales, Pastorellales and Spirochaetes. Single domain globins occur in metazoans and flavohemoglobins are found in fungi, diplomonads and mycetozoans. Although red algae have single domain globins, including 2-over-2 globins, the green algae and ciliates have only 2-over-2 globins. Plants have symbiotic and nonsymbiotic single domain hemoglobins and 2-over-2 hemoglobins. Over 90% of eukaryotes have globins: the nematode Caenorhabditis has the most putative globins, approximately 33. No globins occur in the parasitic, unicellular eukaryotes such as Encephalitozoon, Entamoeba, Plasmodium and Trypanosoma. CONCLUSION: Although Bacteria have all three types of globins, Archaeado not have flavohemoglobins and Eukaryotes lack globin coupled sensors. Since the hemoglobins in organisms other than animals are enzymes or sensors, it is likely that the evolution of an oxygen transport function accompanied the emergence of multicellular animals.

Phylogeny of Echinoderm Hemoglobins.

BACKGROUND: Recent genomic information has revealed that neuroglobin and cytoglobin are the two principal lineages of vertebrate hemoglobins, with the latter encompassing the familiar myoglobin and α-globin/ß-globin tetramer hemoglobin, and several minor groups. In contrast, very little is known about hemoglobins in echinoderms, a phylum of exclusively marine organisms closely related to vertebrates, beyond the presence of coelomic hemoglobins in sea cucumbers and brittle stars. We identified about 50 hemoglobins in sea urchin, starfish and sea cucumber genomes and transcriptomes, and used Bayesian inference to carry out a molecular phylogenetic analysis of their relationship to vertebrate sequences, specifically, to assess the hypothesis that the neuroglobin and cytoglobin lineages are also present in echinoderms. RESULTS: The genome of the sea urchin Strongylocentrotus purpuratus encodes several hemoglobins, including a unique chimeric 14-domain globin, 2 androglobin isoforms and a unique single androglobin domain protein. Other strongylocentrotid genomes appear to have similar repertoires of globin genes. We carried out molecular phylogenetic analyses of 52 hemoglobins identified in sea urchin, brittle star and sea cucumber genomes and transcriptomes, using different multiple sequence alignment methods coupled with Bayesian and maximum likelihood approaches. The results demonstrate that there are two major globin lineages in echinoderms, which are related to the vertebrate neuroglobin and cytoglobin lineages. Furthermore, the brittle star and sea cucumber coelomic hemoglobins appear to have evolved independently from the cytoglobin lineage, similar to the evolution of erythroid oxygen binding globins in cyclostomes and vertebrates. CONCLUSION: The presence of echinoderm globins related to the vertebrate neuroglobin and cytoglobin lineages suggests that the split between neuroglobins and cytoglobins occurred in the deuterostome ancestor shared by echinoderms and vertebrates.

An electrospray ionization mass spectrometric study of the subunit structure of the giant hemoglobin from the leech Nephelopsis oscura.

The subunit structure of the giant, extracellular hexagonal bilayer (HBL) hemoglobin (Hb) from the leech Nephelopsis oscura was investigated by electrospray ionization mass spectrometry (ESI-MS) employing a maximum entropy deconvolution of its complex, multiply charged ESI spectra. The denatured unreduced Hb consisted of three monomer globin chains (M), a1 = 16535 Da, a2 = 17171 Da and a3 = 17315 Da, five nonglobin linker chains, L1 = 24512 Da, L2 = 24586 Da, L3 = 24979 Da, L4 = 25006 Da, and L5 = 25566 Da and two subunits of 32950 Da and 33125 Da. ESI-MS of the denatured, reduced Hb showed that the latter were disulfide-bonded heterodimers (D) of globin chains b1 = 16322 Da and b2 = 16499 Da with chain c = 16632 Da. Time-of-flight ESI-MS of the Hb at pH 3.8, 4.5, 5.0, 5.8 and 7.0 revealed a distribution of charge states from 32(+) to 37(+) with masses decreasing from 211 to 208.5 kDa with increase in cone voltage from 60 to 160 V, indicating the presence of a subassembly comprising 12 globin chains. The subunit composition 6M + 3D + 12h, where M = 16993 Da and D = 33004 Da are the weighted masses and h = 616.5 Da, provides a calculated mass, 208.37 kDa that is closest to 208.5 kDa. Our experimental findings are consistent with the bracelet model of HBL Hbs, verified by the recent low-resolution crystal structure of Lumbricus Hb, wherein an HBL arrangement of 12 globin dodecamer subassemblies is tethered to a central complex of 36 linker chains for a total mass of 208.37 x 12 + 24.94 x 36 = 3398 kDa.

The stoichiometry of the four linker subunits of Lumbricus terrestris hemoglobin suggests an asymmetric distribution.

The extracellular, giant ( approximately 3.6 MDa) hexagonal bilayer hemoglobin of the earthworm Lumbricus terrestris consists of 12 dodecamers of globin chains tethered to a central complex of 36 non-globin, linker chains (24-32 kDa). Four types of linker chains L1-L4 have been detected by electrospray ionization (ESI) and by matrix-assisted laser desorption ionization (MALDI) mass spectrometry (MS) and isolated by reversed phase high pressure liquid chromatography (HPLC). Deconvolution of the HPLC elution profile and of the MS spectra provided the following individual linker contents, expressed as percent of the sum of the four linker peak areas: HPLC-21% L1, 37% L2, 23% L3 and 19% L4, MALDI-47% L1, 29% L2, 16% L3 and 8% L4; ESI-24% L1, 16% L2, 40% L3 and 20% L4; respectively. Comparison with electrophoretic results revealed a surprising lack of overall agreement between all the methods. The calculated mean values of the available linker contents were found to be 32+/-12% L1, 28+/-9% L2, 27+/-10% L3 and 13+/-7% L4, suggesting the following relative stoichiometry: L1: L2: L3: L4 approximately 1: 1: 1: 0.5. With a total of 36 linkers, a hexagonally symmetric distribution of each of the four linker chains is impossible. Thus, the asymmetric linker distribution provides an explanation for the existence of a large dipole moment of Lumbricus terrestris hemoglobin, 17,300+/-2300 Da (Takashima et al., 1999).

Hemoglobins in the genome of the cryptomonad Guillardia theta.

Cryptomonads, are a lineage of unicellular and mostly photosynthetic algae, that acquired their plastids through the "secondary" endosymbiosis of a red alga - and still retain the nuclear genome (nucleomorph) of the latter. We find that the genome of the cryptomonad Guillardia theta comprises genes coding for 13 globin domains, of which 6 occur within two large chimeric proteins. All the sequences adhere to the vertebrate 3/3 myoglobin fold. Although several globins have no introns, the remainder have atypical intron locations. Bayesian phylogenetic analyses suggest that the G. theta Hbs are related to the stramenopile and chlorophyte single domain globins.

Microbial eukaryote globins.

A bioinformatics survey of about 120 protist and 240 fungal genomes and transcriptomes revealed a broad array of globins, representing five of the eight subfamilies identified in bacteria. Most conspicuous is the absence of protoglobins and globin-coupled sensors, except for a two-domain globin in Leishmanias, that comprises a nucleotidyl cyclase domain, and the virtual absence of truncated group 3 globins. In contrast to bacteria, co-occurrence of more than two globin subfamilies appears to be rare in protists. Although globins were lacking in the Apicomplexa and the Microsporidia intracellular pathogens, they occurred in the pathogenic Trypanosomatidae, Stramenopiles and certain fungi. Flavohaemoglobins (FHbs) and related single-domain globins occur across the protist groups. Fungi are unique in having FHbs co-occurring with sensor single-domain globins (SSDgbs). Obligately biotrophic fungi covered in our analysis lack globins. Furthermore, SSDgbs occur only in a heterolobosean amoeba, Naegleria and the stramenopile Hyphochytrium. Of the three subfamilies of truncated Mb-fold globins, TrHb1s appear to be the most widespread, occurring as multiple copies in chlorophyte and ciliophora genomes, many as multidomain proteins. Although the ciliates appear to have only TrHb1s, the chlorophytes have Mb-like globins and TrHb2s, both closely related to the corresponding plant globins. The presently available number of protist genomes is inadequate to provide a definitive census of their globins. Bayesian molecular analyses of single-domain 3/3 Mb-fold globins suggest a close relationship of chlorophyte and haptophyte globins, including choanoflagellate and Capsaspora globins to land plant symbiotic and non-symbiotic haemoglobins and to vertebrate neuroglobins.

A phylogenetic analysis of the globins in fungi.

BACKGROUND: All globins belong to one of three families: the F (flavohemoglobin) and S (sensor) families that exhibit the canonical 3/3 α-helical fold, and the T (truncated 3/3 fold) globins characterized by a shortened 2/2 α-helical fold. All eukaryote 3/3 hemoglobins are related to the bacterial single domain F globins. It is known that Fungi contain flavohemoglobins and single domain S globins. Our aims are to provide a census of fungal globins and to examine their relationships to bacterial globins. RESULTS: Examination of 165 genomes revealed that globins are present in >90% of Ascomycota and ~60% of Basidiomycota genomes. The S globins occur in Blastocladiomycota and Chytridiomycota in addition to the phyla that have FHbs. Unexpectedly, group 1 T globins were found in one Blastocladiomycota and one Chytridiomycota genome. Phylogenetic analyses were carried out on the fungal globins, alone and aligned with representative bacterial globins. The Saccharomycetes and Sordariomycetes with two FHbs form two widely divergent clusters separated by the remaining fungal sequences. One of the Saccharomycete groups represents a new subfamily of FHbs, comprising a previously unknown N-terminal and a FHb missing the C-terminal moiety of its reductase domain. The two Saccharomycete groups also form two clusters in the presence of bacterial FHbs; the surrounding bacterial sequences are dominated by Proteobacteria and Bacilli (Firmicutes). The remaining fungal FHbs cluster with Proteobacteria and Actinobacteria. The Sgbs cluster separately from their bacterial counterparts, except for the intercalation of two Planctomycetes and a Proteobacterium between the Fungi incertae sedis and the Blastocladiomycota and Chytridiomycota. CONCLUSION: Our results are compatible with a model of globin evolution put forward earlier, which proposed that eukaryote F, S and T globins originated via horizontal gene transfer of their bacterial counterparts to the eukaryote ancestor, resulting from the endosymbiotic events responsible for the origin of mitochondria and chloroplasts.

The evolution of land plant hemoglobins.

This review discusses the evolution of land plant hemoglobins within the broader context of eukaryote hemoglobins and the three families of bacterial globins. Most eukaryote hemoglobins, including metazoan globins and the symbiotic and non-symbiotic plant hemoglobins, are homologous to the bacterial 3/3-fold flavohemoglobins. The remaining plant hemoglobins are homologous to the bacterial 2/2-fold group 2 hemoglobins. We have proposed that all eukaryote globins were acquired via horizontal gene transfer concomitant with the endosymbiotic events responsible for the origin of mitochondria and chloroplasts. Although the 3/3 hemoglobins originated in the ancestor of green algae and plants prior to the emergence of embryophytes at about 450 mya, the 2/2 hemoglobins appear to have originated via horizontal gene transfer from a bacterium ancestral to present day Chloroflexi. Unlike the 2/2 hemoglobins, the evolution of the 3/3 hemoglobins was accompanied by duplication, diversification, and functional adaptations. Duplication of the ancestral plant nshb gene into the nshb-1 and nshb-2 lineages occurred prior to the monocot-dicot divergence at ca. 140 mya. It was followed by the emergence of symbiotic hemoglobins from a non-symbiotic hemoglobin precursor and further specialization, leading to leghemoglobins in N2-fixing legume nodules concomitant with the origin of nodulation at ca. 60 mya. The transition of non-symbiotic to symbiotic hemoglobins (including to leghemoglobins) was accompanied by the alteration of heme-Fe coordination from hexa- to penta-coordination. Additional genomic information about Charophyte algae, the sister group to land plants, is required for the further clarification of plant globin phylogeny.