Figures
Abstract
Background
The evolution of oxygenic photosynthesis during Precambrian times entailed the diversification of strategies minimizing reactive oxygen species-associated damage. Four families of oxygen-carrier proteins (hemoglobin, hemerythrin and the two non-homologous families of arthropodan and molluscan hemocyanins) are known to have evolved independently the capacity to bind oxygen reversibly, providing cells with strategies to cope with the evolutionary pressure of oxygen accumulation. Oxygen-binding hemerythrin was first studied in marine invertebrates but further research has made it clear that it is present in the three domains of life, strongly suggesting that its origin predated the emergence of eukaryotes.
Results
Oxygen-binding hemerythrins are a monophyletic sub-group of the hemerythrin/HHE (histidine, histidine, glutamic acid) cation-binding domain. Oxygen-binding hemerythrin homologs were unambiguously identified in 367/2236 bacterial, 21/150 archaeal and 4/135 eukaryotic genomes. Overall, oxygen-binding hemerythrin homologues were found in the same proportion as single-domain and as long protein sequences. The associated functions of protein domains in long hemerythrin sequences can be classified in three major groups: signal transduction, phosphorelay response regulation, and protein binding. This suggests that in many organisms the reversible oxygen-binding capacity was incorporated in signaling pathways. A maximum-likelihood tree of oxygen-binding hemerythrin homologues revealed a complex evolutionary history in which lateral gene transfer, duplications and gene losses appear to have played an important role.
Conclusions
Hemerythrin is an ancient protein domain with a complex evolutionary history. The distinctive iron-binding coordination site of oxygen-binding hemerythrins evolved first in prokaryotes, very likely prior to the divergence of Firmicutes and Proteobacteria, and spread into many bacterial, archaeal and eukaryotic species. The later evolution of the oxygen-binding hemerythrin domain in both prokaryotes and eukaryotes led to a wide variety of functions, ranging from protection against oxidative damage in anaerobic and microaerophilic organisms, to oxygen supplying to particular enzymes and pathways in aerobic and facultative species.
Citation: Alvarez-Carreño C, Becerra A, Lazcano A (2016) Molecular Evolution of the Oxygen-Binding Hemerythrin Domain. PLoS ONE 11(6): e0157904. https://doi.org/10.1371/journal.pone.0157904
Editor: Nikolas Nikolaidis, California State University Fullerton, UNITED STATES
Received: April 5, 2016; Accepted: June 7, 2016; Published: June 23, 2016
Copyright: © 2016 Alvarez-Carreño et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: Claudia Alvarez is a doctoral student in the Programa de Doctorado en Ciencias Biomédicas, Universidad Nacional Autónoma de México (UNAM), and received fellowship 290175 from CONACYT.
Competing interests: The authors have declared that no competing interests exist.
Abbreviations: aLRT, approximate likelihood-ratio test; AMP, Adenosine monophosphate; FBXL5, F-box and leucine-rich repeat protein 5; GGDEF, diguanylate cyclase; HAMP, histidine kinases, adenyl cyclases, methyl-accepting proteins and phosphatases; MCP signal, methyl accepting chemotaxis protein; McHr, Methylococcus capsulatus hemerythrin; ML, maximum-likelihood; O2-binding hemerythrin, oxygen-binding hemerythrin; PDB, Protein Data Bank; pHMM, profile Hidden Markov Model
Introduction
Before the evolution of oxygenic photosynthesizers, sources of free oxygen were scarce [1]. Free molecular oxygen constitutes 21% of present-day terrestrial atmosphere and its main source is, essentially, oxygenic photosynthesis. Accumulation of free atmospheric oxygen during the Precambrian [1–3] is, undoubtedly, one of the major changes in the history of the planet and may be considered the most significant biogeochemical process after the origin of life itself. Oxygen-dependent metabolism evolved first in bacteria and is pervasive in contemporary eukaryotes.
The evolution of metazoans was constrained by the oxygen requirements of tissues [4–6]. Therefore, oxygen-carrier proteins that maintain a continuous delivery of oxygen while avoiding autoxidation as well as the formation and accumulation of reactive oxygen species [4] became essential for the development of animals. Four evolutionarily unrelated families of oxygen-carrier proteins are known: hemoglobin, hemerythrin and the two non-homologous families of molluscan and arthropodan hemocyanins. Such diverse assortment can only be understood in terms of the selective pressure imposed on the biosphere by free oxygen.
Hemerythrin is a small 118 amino acid protein classically studied in four metazoan phyla: Sipuncula, Brachiopoda, Priapulida and Annelida [7,8]; it has also been identified in other eukaryotes as well as in bacterial and archaeal genomes [9–12]. However, the molecular evolution of hemerythrin and its relationship to the geochemical history of the Earth has been largely ignored. In this work, we studied the phylogenetic distribution of hemerythrin-like sequences in 2521 completely sequenced bacterial, archaeal and eukaryotic genomes, and correlated the possible evolutionary scenarios with what is currently known about the structure and function of the hemerythrin domain in both prokaryotes and eukaryotes.
Results
Sequence similarity search
Oxygen-binding hemerythrin (O2-binding Hr) homologs were identified based on the statistical significance of aligning O2-binding Hr sequences from two annelid species, Phascolopsis gouldii [13] and Themiste hennahi [14], with non-mutated protein sequences annotated as hemerythrin or hemerythrin-like proteins in the PDB. The sequences of two O2-binding Hrs, from the proteobacteria Methylococcus capsulatus [15] and Desulfovibrio vulgaris [16], were statistically similar to the annelid sequences (Table 1), and were used as queries in subsequent sequence similarity searches. In contrast, the N-terminal domain of the human F-box and leucine-rich repeat protein 5 (FBXL5), which possess an iron-coordination site that does not bind oxygen [17–19], was not significantly similar to any of the annelid hemerythrin sequences used here as reference (Table 1). We constructed a profile Hidden Markov Model (S1 Fig) with 148 nodes, exclusively based on O2-binding Hr homologues, which constitute a divergent monophyletic sub-group of the hemerythrin/HHE (histidine, histidine, glutamic acid) cation-binding domain (Fig 1).
Midpoint-rooted maximum-likelihood tree of single domain hemerythrin/HHE cation-binding domain sequences. Internal nodes with approximate likelihood-ratio test lower than 0.6 were collapsed. Each one of the nodes represents a sequence identified by the Pfam-A hemerythrin/HHE cation-binding domain profile in a database of completely sequenced cellular genomes. Sequences names appear at the tips of the branches. Names in purple indicate sequences also identified by a hand-curated oxygen-binding hemerythrin profile.
Architecture of hemerythrin-containing sequences
O2-binding Hr sequences were identified in 367/2236 bacterial, 21/150 archaeal and 4/135 eukaryotic genomes (Fig 2). Sequences with subject coverage lower than 85% were considered long sequences and the presence of additional protein domains was investigated using the Pfam-A database (Figs 3 and 4). O2-binding Hr homologues were found in the same proportion as single-domain and as long protein sequences in archaeal, bacterial and eukaryotic genomes overall. A total of 56 different domain architectures were identified in long O2-binding Hr sequences (Fig 4 and S2 Fig). The largest group of long O2-binding Hr sequence architectures (56.9%) contained N-terminal or C-terminal regions of variable size with no clear homologs in Pfam-A. The most frequent location of the O2-binding Hr domain in long protein sequences is the protein termini (117 proteins at the C-terminus, 52 proteins at the N-terminus), suggesting a later incorporation by gene fusion events. The most frequent architectures were: 1) N-terminal methyl accepting chemotaxis protein domain (MCP signal) with a C-terminal O2-binding Hr domain (33 protein sequences); 2) N-terminal O2-binding Hr domain with a C-terminal diguanylate cyclase (GGDEF) domain (22 protein sequences); and 3) N-terminal histidine kinases, adenyl cyclases, methyl-accepting proteins and phosphatases (HAMP) domain with a MCP signal and a C-terminal O2-binding Hr domain (14 protein sequences). The Gene Ontology terms associated to the Pfam-A domains identified corresponded to: 1) signal transduction (43%); 2) phosphorelay response regulation (7%); and 3) protein binding (6%).
Phylogenetic tree based on a small subunit rRNA guide tree containing only completely sequenced species. Bacterial and archaeal species are collapsed on the phylum and group level. Eukaryotic species are collapsed together. n: number of species contained within each collapsed branch. The red bar is proportional to the number of species with at least one hemerythrin sequence in each collapsed branch. The total number of genomes with at least one O2-binding hemerytrhin is indicated by a red number next to the bar.
Phylogenetic tree based on a small subunit rRNA guide tree. Branch lengths are arbitrary. Each node corresponds to a completely sequenced species with at least one O2-binding Hr sequence homolog. Species names were replaced by their unique KEGG Organisms code [20]. The height of the bar charts is proportional to the absolute number of O2-binding Hr copies in each category: single domain O2-binding Hr sequences (red), long O2-binding Hr sequences (pink).
Protein domain architectures were obtained as specified in the Methods section. Defined Pfam-A domains are represented by boxes. Only architectures where three or more protein sequences were identified are depicted here. The size of the boxes representing each domian and spacing between contiguous boxes are arbitrary. Hemerythrin domain overlapping other protein domains is represented as half-filled blue boxes.
Phylogenetic distribution of oxygen-binding hemerythrin sequences
Of the 2521 cellular genomes that were analyzed in this work, a total of 392 (367 bacterial, 21 archaeal and 4 eukaryotic genomes) encoded for O2-binding Hr sequences (Fig 2). The number of O2-binding Hr copies in a genome varies widely, in particular, species of Spirochaetes, Chlorobi, Acidobacteria, Firmicutes-Clostridia, and α-, β- and δ- Proteobacteria present several copies of single-domain O2-binding Hr sequences (Fig 3). The greatest number of single-domain O2-binding Hr copies was found in Magnetospirillum magneticum AMB-1, a facultative α-Proteobacteria encoding 15 paralogous sequences (Fig 3).
The sample studied here included archaeal genomes only from the Crenarchaeota and Euryarchaeota phyla due to an under-representation of Archaea in genome databases and the low agreement between the databases consulted (See the Methods section). Within Crenarchaeota, O2-binding Hr was found only in five species of the order Desulfurococcales within long protein sequences. In Euryarchaeota, five species of the Methanomicrobia class contained single-domain O2-binding Hr, while six species of the Methanococci class had either single-domain O2-binding Hr or long sequences consisting of an O2-binding Hr domain with an AMP-binding domain or with a short orphan elongation. Eukaryotic O2-binding Hr sequences were found in Micromonas pusilla CCMP1545 and Acanthamoeba castellanii Neff only as single-domain sequences, and in Nematostella vectensis and Naegleria gruberi as both single-domain and long protein sequences. The profile Hidden Markov Model used in this work identified a negligible number of O2-binding Hr sequences in eight bacterial phyla (Thermotogae, Actinobacteria, Tenericutes, Chlamydiae, Planctomycetes, Bacteroidetes and Fusobacteria) and in the class Bacilli of Firmicutes (Fig 2). With the exception of Fusobacteria and Tenericutes, the homologous hemerythrin/HHE cation-binding domain was found in species from the phylogenetic groups where O2-binding Hr was not present (S3 Fig). In early-branching Bacteria (Aquificae, Chloroflexi, Deinococcus-Thermus and Cyanobacteria) and Firmicutes, single-domain O2-binding Hr accounted for 91.6% of the sequences identified. In contrast, in later bacterial groups (Acidobacteria, Chlorobi, Spirochaetes, and Proteobacteria) only 31.9% of the O2-binding Hr homologues were single-domain sequences.
As shown in Fig 2, the sequences of single-domain HHE cation-binding hemerythrin and single-domain O2-binding Hr, which we have analyzed in this work, are the outcome of an ancient gene duplication predating the explosive divergence of Bacteria during Precambrian times. The maximum-likelihood tree in Fig 2 also shows that the single-domain O2-binding Hr is a well-defined, monophyletic, highly divergent group distinct from the hemerytrin/HHE cation-binding domain.
The tree in Fig 5 shows two well-defined single-domain O2-binding Hr derived groups, a and b, formed by sequences encoded mostly by anaerobic species. The derived clade a includes sequences from Deinococcus-Thermus; Cyanobacteria; Firmicutes; Spirochaetes; and δ-Proteobacteria. The derived clade b is formed by sequences from Acidobacteria; Chlorobi; α-, β-, γ- and δ-Proteobacteria, as well as from the archaeal Methanospirillum hungatei JF-1, where it is present most likely because of a lateral gene event. Sequences in both derived clusters diverged from an early gene duplication of the single-domain O2-binding Hr gene, and also appear to have been subject of lateral gene transfer, gene loss, gene duplication and orthologous replacement events.
Internal nodes with approximate likelihood-ratio test lower than 0.6 were collapsed. Names of the sequences appear at the tips of the branches. Sub-clusters are named by Greek letters on the base of the tree and by letters and numbers in the direved clusters. Node names appear in color when there is at least another copy of single-domain O2-binding Hr in a separate cluster of the tree. Oxygen requirement of the species source, according to the Genomes OnLine Database [21], is indicated by color bars at the tips of the nodes.
There are several cases in which more than one O2-binding Hr sequence can be identified in a single cellular genome. This may be due to lateral gene transfer events or, most likely, to paralogous duplication events (Table 2).
The many cases in which highly similar single-domain O2-binding Hr copies are closely grouped in the same cluster in the phylogenetic tree suggest recent paralogous duplication events. For instance, Magnetospirillum magneticum AMB-1, a facultative proteobacteria, encodes for 15 single-domain O2-binding Hr gene copies, some of which are located in the same cluster, while others are found in distant clusters: amb4136 is represented by a single branch; amb4265, amb2654, amb4171, amb2239 and amb1415 are in cluster ε; amb1987, amb0226, amb1549 and amb2249, in cluster ι; amb3418, amb4296 and amb3966 in cluster η; amb0569 in cluster ζ; and amb1952 in cluster b1, indicating duplication events that occurred at different times during evolution of M. magneticum AMB-1 (Fig 5). Alternatively, the peculiarities of the distribution of one or more O2-binding Hr gene copies can be explained by lateral gene transfer events.
The overall topology of the single-domain O2-binding Hr phylogenetic tree is in clear disagreement with the 16/18S rRNA reference tree (Fig 5). This indicates that the significance of O2-binding Hr as a phylogenetic marker is hindered by a complex history of gene losses, gene duplications, paralogous replacements and lateral gene transfer events. The base of the tree is characterized by a highly branching pattern of single-domain O2-binding Hr sequences encoded by aerobic organisms including Bacteria, Archaea and Eukarya (Fig 5 and Table 3). The ample distribution of highly divergent O2-binding Hr sequences among aerobic organisms is probably best understood by the adaptive value of oxygen-binding proteins in a Precambrian environment that was becoming increasingly oxidizing as time went by.
Agreement with the topology of the species tree was observed mainly at the sub-group level of the O2-binding Hr sequence tree. For instance, as shown in Fig 5 and Table 3, sequences from Aquificae form a minor sub-group (sub-group α), suggesting that a single-domain O2-binding hemerythrin was present in their common ancestor. A comparable situation may be seen in the θ sub-group, which can be interpreted as the outcome of vertically inherited single-domain O2-binding Hr in species of Cyanobacteria subclass Oscillatoriophycideae. In the case of archaeal sequences, Methanococci and Methanomicrobia form two separate groups. Sequences from Methanococci, and bacterial sequences from Anaerolinea thermophila UNI-1, Persephonella marina EX-H1, Dechlorosoma suillum PS and Acidovorax ebreus TPSY form cluster μ. A distinct independent sub-group ν is formed by sequences from Methanomicrobia, and bacterial sequences from Melioribacter roseus P3M-2, Geobacter lovleyi SZ and Pelobacter carbinolicus DSM 2380. The topology of the single domain O2-binding Hr tree, together with the absence of single-domain O2-binding Hr in most archaeal groups, could indicate the acquisition by Methanomicrobia and Methanococci of the single-domain O2-binding Hr sequences in two independent events of lateral gene transfer. Eukaryotic sequences from Nematostella vectensis and Naegleria gruberi cluster together in sub-group λ.
Discussion
The HHE cation-binding domain was first predicted by bioinformatics methods as a domain composed of two helical regions and a conserved HHE cation-binding site [22]. The hemerythrin-like domain family is a repetition of the HHE cation-binding domain, which folds into an up-and-down bundle of four left-handed helices [14,23]. The molecular function of proteins containing the hemerythrin/HHE cation-binding domain, including metazoan oxygen-carrier hemerythrins, is often related to O2 or reactive oxygen species responses [10,16,24–29].
The search for O2-binding Hr domain homologs was performed using a manually curated profile Hidden Markov Model, and resulted in the identification of sequences of 86 to 2425 amino acids length, using a profile that contained 148 nodes (S1 Fig), highly weighting the position of the iron-coordinating amino acids in X-ray solved structures of O2-binding Hrs. Sequences with subject coverage lower than 85% were considered long sequences. To identify the presence of possible additional domains, we searched the protein profile database Pfam-A, which confirmed known additional domains in 37.2% of long sequences. In agreement with previous bioinformatics searches [9,11], the most frequent domain found in long O2-binding Hr sequences was the MCP signal domain (Fig 4). The characterization of the hemerythrin domain of the methyl-accepting chemotaxis protein dcrH from Desulfovibrio vulgaris has shown that the four-helix bundle fold and the amino acids of the active site are conserved [24]. It has been proposed that the hemerythrin domain in this structure could have a role in signal transduction [30].
O2-binding Hr sequences were identified in two archaeal groups (Euryarchaeota and Crenarchaeota), ten bacterial groups (Aquificae, Deinococcus-Thermus, Cyanobacteria, Chloroflexi, Firmicutes, Acidobacteria, Chlorobi, Spirochaetes and Proteobacteria), and four eukaryotic species (Naegleria gruberi, Micromonas pusilla CCMP1545, Nematostella vectensis and Acanthamoeba castellanii). It is also known to be present in marine invertebrates (Sipuncula, Brachiopoda, Priapulida and Annelida) [8]. By far, the highest number of O2-binding Hr sequences is found in Proteobacteria (Figs 2 and 3).
Single-domain and long O2-binding Hr-containing sequences are not randomly distributed in the phylogenetic tree, but exhibit a differential distribution among taxonomic groups, particularly in Bacteria, where two separate groups can be distinguished. In the first bacterial group, that includes deep-branching bacteria and the Firmicutes-Clostridia clade, O2-binding Hr homologues are predominantly single-domain sequences (91.6%). The expression of the single-domain HerA hemerythrin in the microaerophilic bacteria Campylobacter jejuni [26] reduces its susceptibility to oxygen and hydrogen peroxide-mediated damage of two iron-sulphur cluster enzymes (pyruvate:acceptor oxydoreductase and 2-oxoglutarate:acceptor oxidoreductase). This suggest that single-domain O2-binding Hr may also be involved in the prevention of oxygen-mediated damage in microaerophilic aquificales and in anaerobic clostridia, and probably reflect an evolutionary adaptation to the presence of free molecular oxygen. The absence of O2-binding Hr homologues in Mollicutes and the class Bacilli of Firmicutes, two bacterial groups closely related to clostridial species, is probably due to secondary loss.
The second bacterial group includes Acidobacteria, Chlorobi, Spirochaetes and Proteobacteria. Within this group, long O2-binding Hr sequences (68.1%) are more frequent than single-domain (31.9%) O2-binding Hr. In Proteobacteria, which is the most diverse bacterial phylum, the study of the O2-binding Hr domain led to the identification of three different functions: a) as a 135-amino acid domain in a chemotactic protein of Desulfovubrio vulgaris [24]; b) as a single-domain protein McHr in Methylococcus capsulatus, that supplies oxygen to the membrane-bound methane monooxygenase [10]; and c) as a single-domain protein HerA in Campylobacter jejuni, which protects iron-sulphur cluster enzymes from oxidative damage [26]. This exemplifies how the O2-binding Hr domain may have been coopted into specific physiological functions by different species, particularly in later-branching bacteria, where it was incorporated into a wide variety of sequences. In some cases, the additional sequences within long O2-binding Hr sequences are group-specific. For instance, long O2-binding Hr sequences from Acidobacteria and Chlorobi exhibit elongations that do not match domains in the Pfam-A database, and two architectures of the long O2-binding Hr sequences from Spirochaetes genomes contain a domain called Cache_1, which is an acronym for calcium channels and chemotaxis receptor.
Variations at the iron-coordinating amino acid residues in 106 single-domain O2-binding Hr sequences were observed (S4 Fig). This is quite evident in sequences from the Deinococcus-Thermus, Cyanobacteria, Firmicutes-Clostridia and from the α-, β- and δ- Proteobacteria clades. Amino acid substitutions may modify the reversible oxygen-binding ability of hemerythrins [31], or even the native metal-ion preference. Molecular characterization of the variants could clarify whether sequence variations at the iron-coordination site produce loss of function, neo-functionalization, or alternative iron-binding mechanisms.
Hemerythrin homologues have an ample biological distribution, and are present in the three domains of life. However, as shown here, their distribution is not universal, which is consistent with previous genomic searches, and may be explained by frequent gene losses [11,12]. Although the topology of the O2-binding Hr tree indicates intense horizontal gene transfers, gene duplications and differential gene loss (Fig 5), it is clear that O2-binding Hr is widely distributed in Clostridia and all five classes of Proteobacteria, suggesting that it is an ancient Precambrian trait that was already present before the divergence of Firmicutes and Proteobacteria. Alternatively, the evolution of O2-binding Hr may have occurred in a divergent phylogenetic group, and was horizontally transferred, through independent events, to the individual orders of Clostridia and/or Proteobacteria. This alternative scenario represents a less parsimonious explanation to the overall O2-binding Hr distribution. Instead, horizontal gene transfer could explain the cases where the presence of O2-binding Hr is circumscribed to a particular family or genus, as appears to be the case for archaean O2-binding Hr.
Conclusion
The four known evolutionary unrelated oxygen-carrier protein families (hemoglobin, hemerytrhin and the two non-homologous molluscan and arthropodan hemocyanins) represent a polyphyletic response to the selective pressure imposed by the accumulation of free oxygen during the Precambrian. Like other molecular mechanisms involved in the protection and repair of oxygen-induced damage that evolved during Precambrian times [32], the ample biological distribution of O2-binding Hr probably reflects its adaptive significance in an environment became increasingly oxidizing. The biological group where hemerythrin first evolved remains unknown, but the distribution of O2-binding Hr sequences suggest that its capacity to reversibly bind oxygen was exploited first by prokaryotic species to fulfill a range of different needs, ranging from protection against oxidative damage to oxygen supply to particular enzymes and pathways. More specifically, the available data suggests that it may have originated prior to the divergence of Firmicutes and Proteobacteria. Thereafter, the incorporation of the O2-binding Hr domain into preexisting proteins, combined with other mechanisms involved in the protection against oxidative damage, may have allowed a functional diversification of the protein repertoire particularly in the proteobacteria, one of the most diverse groups of bacteria. The evolution of O2-binding Hr sequences appears to parallel the evolution of strategies allowing the incorporation of oxygen to biological processes as a consequence of accumulation of free oxygen in the Precambrian oceans and atmosphere.
Methods
Sequence similarity searches
Eleven non-mutated hemerythrin and hemerythrin-like amino acid sequences retrieved from the Protein Data Bank (PDB) database [33] were compared to the hemerythrin sequence from Themiste hennahi (PDB codes 2MHR) and from Phascolopsis gouldii (PDB code 1I4Y) with the program PRSS version 36.3.8c (number of shuffles: 200; scoring matrix: Blosum50; open gap penalty: -10; extension gap penalty: -2) from UVa FASTA Server [34,35]. Expect value lower than 1e-3 of one-to-one alignments was used as cutoff for homologous sequences.
An initial BLAST search against the Kegg Genes database release 75.1[20,36] identified 365 protein sequences with expect value < 1e-5 and subject coverage > 95%. The sequences were aligned with MAFFT G-INS-I algorithm [37] with default parameters (gap opening penalty: 1.53; offset: 0.0). A profile Hidden Markov Model (pHMM) was constructed with HMMER 3.1b1 [38]. The Pfam hemerythrin profile (http://pfam.xfam.org/family/hemerythrin) [39] has a low restriction to include sequences with changes at the iron coordination amino acids observed in oxygen-binding hemerythrins (S1 Fig). The distribution of the hemerythrin domain homologues found by the Pfam profile is shown in S3 Fig. The hand curated hemerythrin pHMM was scanned against 2521 genomes from the Kegg Genes database release 75.1 that were unequivocally traced to an entry of the SSU RefNR99 guide tree from the SILVA SSU database release 119 [40] based on taxonomy identifiers from both databases.
Domain annotation and associated biological processes
Each sequence was scanned with hmmscan from the HMMER 3.1b1 software using the database of protein families Pfam-A version 28.0 [39]. The hemerythrin profile in Pfam-A was substituted with the pHMM that we generated. Domain overlapping was solved with a pearl script preferring domains with lower expect values. Gene Ontology information on molecular function of Pfam domains was obtained from the Pfam web page [41].
Phylogenetic analysis
The maximum-likelihood (ML) tree of single-domain oxygen-binding hemerythrin sequences and of hemerythrin HHE cation-binding domain sequences were constructed with PHYML [42]. The parameters for the construction of the trees (Table 4) were automatically selected with SMS: Smart Model Selection using the Akaike Information Criterion at the South of France bioinformatics platform (http://www.atgc-montpellier.fr). Branch support was given by the approximate likelihood-ratio test (aLRT) for branches [43]. Internal nodes with low support (< 0.6) were collapsed. Additional bootstrap support is provided by a 100-replicates as shown in S5 Fig. The multiple sequence alignments were constructed with MAFFT, L-INS-i algorithm using the default parameters (gap opening penalty: 1.53; offset: 0.0). The maximum-likelihood tree of hemerythrin HHE cation-binding domain sequences was midpoint-rooted; the maximum-likelihood of single-domain oxygen-binding hemerythrin sequences was rooted following the topology of the hemerythrin HHE cation-binding domain tree. Tree visualization was made with Interactive Tree Of Life V1.0 [44].
Supporting Information
S1 Fig. Conservation of the iron-coordinating amino acids is reflected on the calculated profile Hidden Markov Model.
Logo representation of the multiple sequence alignment used as seed to calculate the profile Hidden Markov models. The vertical axis indicates the information content of a sequence position. The one-letter notation for amino acid sequences was used. Glutamic and aspartic acid (purple), histidine (cyan). (A) O2-binding Hr model. Positions of the iron-coordinating amino acids are indicated by arrows. Position and length of helical structures were predicted by Ali2D. (B) Pfam-A hemerythrin model.
https://doi.org/10.1371/journal.pone.0157904.s001
(TIF)
S2 Fig. Schematic reperesentation of the domain architectures of long O2-binding hemerythrin sequences.
Protein domains identified in long O2-binding hemerythrin sequences are designated by their short Pfam-A family name. The number of sequences showing a particular architecture is indicated after a tabular space.
https://doi.org/10.1371/journal.pone.0157904.s002
(TIF)
S3 Fig. Proportion of genomes with hemerythrin HHE cation-binding domain.
Phylogenetic tree based on a small subunit rRNA guide tree containing only completely sequenced species. Bacterial and archaeal species are collapsed on the phylum level. Eukaryotic species are collapsed together. n: number of species within the collapsed branch. The red bar is proportional to the number of species with at least one hemerythrin HHE cation-binding domain sequence in each group. The total number of genomes with at least one HHE cation-binding domain sequence is indicated by a purple number next to the bar.
https://doi.org/10.1371/journal.pone.0157904.s003
(TIF)
S4 Fig. Number of single-domain hemerythrin copies with modifications of the iron-coordinating amino acids within each phylogenetic group.
https://doi.org/10.1371/journal.pone.0157904.s004
(TIF)
S5 Fig. Hemerythrin /HHE cation.binding phylogeny (100 bootstrap replicates).
Sequences names are the same as in Fig 1.
https://doi.org/10.1371/journal.pone.0157904.s005
(TIF)
Acknowledgments
Claudia Alvarez is a doctoral student in the Programa de Doctorado en Ciencias Biomédicas, Universidad Nacional Autónoma de México (UNAM), and received fellowship 290175 from CONACYT. We are indebted to J. Peter Gogarten for insightful discussions and to Luis Montaño, Sara Islas and Ricardo Hernández for their help with the manuscript. Support from DGAPA-PAPIIT(IN223916) is gratefully acknowledged.
Author Contributions
Conceived and designed the experiments: AL. Performed the experiments: CAC AB AL. Analyzed the data: CAC AB AL. Contributed reagents/materials/analysis tools: CAC AB AL. Wrote the paper: CAC AB AL.
References
- 1. Canfield DE. The early history of atmospheric oxygen: Homage to Robert M. Garrels. Annu Rev Earth Planet Sci. 2005;33:1–36.
- 2. Bekker A, Holland HD, Wang P-L, Rumble D, Stein HJ, Hannah JL, et al. Dating the rise of atmospheric oxygen. Nature. 2004;427(6970):117–20. pmid:14712267
- 3. Stüeken EE, Buick R, Anbar AD. Selenium isotopes support free O2 in the latest Archean. Geology [Internet]. 2015;43(3):259–62. Available: http://geology.gsapubs.org/content/43/3/259.abstract?rss=1\nhttp://geology.gsapubs.org/cgi/doi/10.1130/G36218.1
- 4. Kurtz DM. Oxygen-carrying proteins: three solutions to a common problem. Essays Biochem. 1999;34:85–100. pmid:10730190
- 5. Terwilliger NB. Functional adaptations of oxygen-transport proteins. J Exp Biol. 1998 Apr;201(Pt 8):1085–98. pmid:9510522
- 6. Planavsky NJ, Reinhard CT, Wang X, Thomson D, McGoldrick P, Rainbird RH, et al. Low Mid-Proterozoic atmospheric oxygen levels and the delayed rise of animals. Science (80-) [Internet]. 2014;346(6209):635–8. Available: http://www.sciencemag.org/content/346/6209/635.abstract\nhttp://www.sciencemag.org/content/346/6209/635.full.pdf pmid:25359975
- 7. Cook SF. The action of potassium cyanide and potassium ferricyanide on certain respiratory pigments. J Gen Physiol. 1928;11(4):339–48. pmid:19872401
- 8. Vanin S, Negrisolo E, Bailly X, Bubacco L, Beltramini M, Salvato B. Molecular evolution and phylogeny of sipunculan hemerythrins. J Mol Evol. 2006 Jan;62(1):32–41. pmid:16362484
- 9. Bailly X, Vanin S, Chabasse C, Mizuguchi K, Vinogradov SN. A phylogenomic profile of hemerythrins, the nonheme diiron binding respiratory proteins. BMC Evol Biol. 2008 Jan;8:244. pmid:18764950
- 10. Karlsen OA, Ramsevik L, Bruseth LJ, Larsen Ø, Brenner A, Berven FS, et al. Characterization of a prokaryotic haemerythrin from the methanotrophic bacterium Methylococcus capsulatus (Bath). FEBS J. 2005 May;272(10):2428–40. pmid:15885093
- 11. Martín-Durán JM, De Mendoza A, Sebé-Pedrós A, Ruiz-Trillo I, Hejnol A. A broad genomic survey reveals multiple origins and frequent losses in the evolution of respiratory hemerythrins and hemocyanins. Genome Biol Evol. 2013;5:1435–42. pmid:23843190
- 12. French CE, Bell JML, Ward FB. Diversity and distribution of hemerythrin-like proteins in prokaryotes. FEMS Microbiol Lett. 2008 Feb;279(2):131–45. pmid:18081840
- 13. Farmer CS, Kurtz DM Jr, Liu ZJ, Wang BC, Rose J, Ai J, et al. The crystal structures of Phascolopsis gouldii wild type and L98Y methemerythrins: structural and functional alterations of the O2 binding pocket. J Biol Inorg Chem [Internet]. 2001;6(4):418–29. Available: http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&dopt=Citation&list_uids=11372200 pmid:11372200
- 14. Sheriff S, Hendrickson WA, Smith JL. Structure of myohemerythrin in the azidomet state at 1.7/1.3 A resolution. J Mol Biol [Internet]. 1987;197(2):273–96. Available: http://www.ncbi.nlm.nih.gov/pubmed/3681996 pmid:3681996
- 15. Chen KH-C, Chuankhayan P, Wu H-H, Chen C-J, Fukuda M, Yu SS-F, et al. The bacteriohemerythrin from Methylococcus capsulatus (Bath): Crystal structures reveal that Leu114 regulates a water tunnel. J Inorg Biochem. Elsevier Inc.; 2015;15–7.
- 16. Isaza CE, Silaghi-Dumitrescu R, Iyer RB, Kurtz DM, Chan MK. Structural basis for O2 sensing by the hemerythrin-like domain of a bacterial chemotaxis protein: Substrate tunnel and fluxional N terminus. Biochemistry. 2006;45(30):9023–31. pmid:16866347
- 17. Rouault TA. An ancient gauge for iron leaps in translational elongation. Science (80-). 2009;326(October):676–7.
- 18. Shu C, Sung MW, Stewart MD, Igumenova TI, Tan X, Li P. The structural basis of iron sensing by the human F-box protein FBXL5. ChemBioChem. 2012;13(6):788–91. pmid:22492618
- 19. Thompson JW, Salahudeen A, Chollangi S, Ruiz JC, Brautigam C, Makris TM, et al. Structural and molecular characterization of iron-sensing hemerythrin-like domain within F-box and leucine-rich repeat protein 5 (FBXL5). J Biol Chem [Internet]. 2012;287(10):7357–65. Available: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3293540&tool=pmcentrez&rendertype=abstract pmid:22253436
- 20. Kanehisa M, Sato Y, Kawashima M, Furumichi M, Tanabe M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res [Internet]. 2016;44(D1):D457–62. Available: pmid:26476454
- 21. Reddy TBK, Thomas AD, Stamatis D, Bertsch J, Isbandi M, Jansson J, et al. The Genomes OnLine Database (GOLD) v.5: A metadata management system based on a four level (meta)genome project classification. Nucleic Acids Res. 2015;43(D1):D1099–106.
- 22. Yeats C, Bentley S, Bateman A. New knowledge from old: in silico discovery of novel protein domains in Streptomyces coelicolor. BMC Microbiol [Internet]. 2003;3:3. Available: http://www.ncbi.nlm.nih.gov/pubmed/12625841 pmid:12625841
- 23.
Hendrickson WA, Smit JL, Sheriff S. Structure and Function of Hemerythrins. Respiratory pigments in animals [Internet]. Springer Berlin Heidelberg; 1985. p. 1–8. Available: https://books.google.com.mx/books?id=AljmCAAAQBAJ
- 24. Xiong J, Kurtz DM, Ai J, Sanders-loehr J. A hemerythrin-like domain in a bacterial chemotaxis protein. Biochemistry. 2000 May 2;39(17):5117–25. pmid:10819979
- 25. Vashisht AA, Zumbrennen KB, Huang X, Powers DN, Durazo A, Sun D, et al. Control of iron homeostasis by an iron-regulated ubiquitin ligase. Science (80-). 2009;326(5953):718–21. pmid:19762596
- 26. Kendall JJ, Barrero-Tobon AM, Hendrixson DR, Kelly DJ. Hemerythrins in the microaerophilic bacterium Campylobacter jejuni help protect key iron-sulphur cluster enzymes from oxidative damage. Environ Microbiol. 2014 Apr;16(4):1105–21. pmid:24245612
- 27. Li X, Tao J, Hu X, Chan J, Xiao J, Mi K. A bacterial hemerythrin-like protein MsmHr inhibits the SigF-dependent hydrogen peroxide response in mycobacteria. Front Microbiol. 2015;5(January):1–11.
- 28. Chow ED, Liu OW, O’Brien S, Madhani HD. Exploration of whole-genome responses of the human AIDS-associated yeast pathogen Cryptococcus neoformans var grubii: Nitric oxide stress and body temperature. Curr Genet. 2007;52(3–4):137–48. pmid:17661046
- 29. Li X, Li J, Hu X, Huang L, Xiao J, Chan J, et al. Differential roles of the hemerythrin-like proteins of Mycobacterium smegmatis in hydrogen peroxide and erythromycin susceptibility. Sci Rep [Internet]. Nature Publishing Group; 2015;5(August):16130. Available: pmid:26607739
- 30. Onoda A, Okamoto Y, Sugimoto H, Shiro Y, Hayashi T. Crystal structure and spectroscopic studies of a stable mixed-valent state of the hemerythrin-like domain of a bacterial chemotaxis protein. Inorg Chem. 2011 Jun 6;50(11):4892–9. pmid:21528842
- 31. Okamoto Y, Onoda A, Sugimoto H, Takano Y, Hirota S, Kurtz DM, et al. Crystal structure, exogenous ligand binding, and redox properties of an engineered diiron active site in a bacterial hemerythrin. Inorg Chem. 2013;52(22):13014–20. pmid:24187962
- 32. Delaye L, Becerra A, Orgel L, Lazcano A. Molecular evolution of peptide methionine sulfoxide reductases (MsrA and MsrB): On the early development of a mechanism that protects against oxidative damage. J Mol Evol. 2007;64(1):15–32. pmid:17180746
- 33. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, et al. The Protein Data Bank. Nucleic Acids Res. 2000;28(1):235–42. pmid:10592235
- 34. Smith T, Waterman M. Identification of common molecular subsequences. J Mol Biol [Internet]. 1981;147(1):195–7. Available: http://dx.doi.org/10.1016/0022-2836(81)90087-5 pmid:7265238
- 35. Pearson WR. Searching protein sequence libraries: Comparison of the sensitivity and selectivity of the Smith-Waterman and FASTA algorithms. Genomics. 1991;11(3):635–50. pmid:1774068
- 36. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic Local Alignment Search Tool. J Mol Biol. 1990;403–10.
- 37. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–80. pmid:23329690
- 38. Eddy SR. Accelerated profile HMM searches. PLoS Comput Biol. 2011 Oct;7(10):e1002195. pmid:22039361
- 39. Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, et al. Pfam: The protein families database. Nucleic Acids Res. 2014;42(D1):222–30.
- 40. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucleic Acids Res. 2013;41(D1):590–6.
- 41. Mitchell A, Chang H-Y, Daugherty L, Fraser M atthe., Hunter S, Lopez R, et al. The InterPro protein families database: the classification resource after 15 years. Nucleic Acids Res. 2014;43(D1):D213–21.
- 42. Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52(5):696–704. pmid:14530136
- 43. Anisimova M, Gascuel O. Approximate likelihood-ratio test for branches: a fast, accurate, and powerful alternative. 2006;55(4):539–52. pmid:16785212
- 44. Letunic I, Bork P. Interactive Tree of Life v2: Online annotation and display of phylogenetic trees made easy. Nucleic Acids Res. 2011;39(SUPPL. 2):475–8.