Molecular Evolution of the Oxygen-Binding Hemerythrin Domain

Conserved C-Terminal Tail Is Responsible for Membrane Localization and Function of Pseudomonas aeruginosa Hemerythrin

Many bacteria have hemerythrin (Hr) proteins that bind O2, including Pseudomonas aeruginosa, in which microoxia-induced Hr (Mhr) provide fitness advantages under microoxic conditions. Mhr has a 23 amino-acid extension at its C-terminus relative to a well-characterized Hr from Methylococcus capsulatus, and similar extensions are also found in Hrs from other bacteria. The last 11 amino acids of this extended, C-terminal tail are highly conserved in gammaproteobacteria and predicted to form a helix with positively charged and hydrophobic faces. In cellular fractionation assays, wild-type (WT) Mhr was found in both membrane and cytosolic fractions, while a MhrW143* variant lacking the last 11 residues was largely in the cytosol and did not complement Mhr function in competition assays. MhrL112Y, a variant that has a much longer-lived O2-bound form, was fully functional and had a similar localization pattern to that of WT Mhr. Both MhrW143* and MhrL112Y had secondary structures, stabilities, and O2-binding kinetics similar to those of WT Mhr. Fluorescence studies revealed that the C-terminal tail, and particularly the fragment corresponding to its last 11 residues, was sufficient and necessary for association with lipid vesicles. Molecular dynamics simulations and subsequent cellular analysis of Mhr variants have demonstrated that conserved, positively charged residues in the tail are important for Mhr interactions with negatively charged membranes and the contribution of this protein to competitive fitness. Together, these data suggest that peripheral interactions of Mhr with membranes are guided by the C-terminal tail and are independent of O2-binding.

Identification of Type 4B Secretion System Substrates That Are Conserved among Coxiella burnetii Genomes and Promote Intracellular Growth

Coxiella burnetii is a Gram-negative pathogen that infects a variety of mammalian hosts. Infection of domesticated ewes can cause fetal abortion, whereas acute human infection normally manifests as the flu-like illness Q fever. Successful host infection requires replication of the pathogen within the lysosomal Coxiella-containing vacuole (CCV). The bacterium encodes a type 4B secretion system (T4BSS) that delivers effector proteins into the host cell. Disruption of C. burnetii T4BSS effector export abrogates CCV biogenesis and bacterial replication. Over 150 C. burnetii T4BSS substrates have been designated often based on heterologous protein translocation by the Legionella pneumophila T4BSS. Cross-genome comparisons predict that many of these T4BSS substrates are truncated or absent in the acute-disease reference strain C. burnetii Nine Mile. This study investigated the function of 32 proteins conserved among diverse C. burnetii genomes that are reported to be T4BSS substrates. Despite being previously designated T4BSS substrates, many of the proteins were not translocated by C. burnetii when expressed fused to the CyaA or BlaM reporter tags. CRISPR interference (CRISPRi) indicated that of the validated C. burnetii T4BSS substrates, CBU0122, CBU1752, CBU1825, and CBU2007 promote C. burnetii replication in THP-1 cells and CCV biogenesis in Vero cells. When expressed in HeLa cells tagged at its C or N terminus with mCherry, CBU0122 localized to the CCV membrane and the mitochondria, respectively. Collectively, these data further define the repertoire of bona fide C. burnetii T4BSS substrates. IMPORTANCE Coxiella burnetii secretes effector proteins via a T4BSS that are required for successful infection. Over 150 C. burnetii proteins are reported to be T4BSS substrates and often by default considered putative effectors, but few have assigned functions. Many C. burnetii proteins were designated T4BSS substrates using heterologous secretion assays in L. pneumophila and/or have coding sequences that are absent or pseudogenized in clinically relevant C. burnetii strains. This study examined 32 previously reported T4BSS substrates that are conserved among C. burnetii genomes. Of the proteins tested that were previously designated T4BSS substrates using L. pneumophila, most were not exported by C. burnetii. Several T4BSS substrates that were validated in C. burnetii also promoted pathogen intracellular replication and one trafficked to late endosomes and the mitochondria in a manner suggestive of effector activity. This study identified several bona fide C. burnetii T4BSS substrates and further refined the methodological criteria for their designation.

A Short Tale of the Origin of Proteins and Ribosome Evolution

Proteins are the workhorses of the cell and have been key players throughout the evolution of all organisms, from the origin of life to the present era. How might life have originated from the prebiotic chemistry of early Earth? This is one of the most intriguing unsolved questions in biology. Currently, however, it is generally accepted that amino acids, the building blocks of proteins, were abiotically available on primitive Earth, which would have made the formation of early peptides in a similar fashion possible. Peptides are likely to have coevolved with ancestral forms of RNA. The ribosome is the most evident product of this coevolution process, a sophisticated nanomachine that performs the synthesis of proteins codified in genomes. In this general review, we explore the evolution of proteins from their peptide origins to their folding and regulation based on the example of superoxide dismutase (SOD1), a key enzyme in oxygen metabolism on modern Earth.

Structure-Function Relationships of Oxygen Transport Proteins in Marine Invertebrates Enduring Higher Temperatures and Deoxygenation

AbstractPredictions for climate change-to lesser and greater extents-reveal a common scenario in which marine waters are characterized by a deadly trio of stressors: higher temperatures, lower oxygen levels, and acidification. Ectothermic taxa that inhabit coastal waters, such as shellfish, are vulnerable to rapid and prolonged environmental disturbances, such as heatwaves, pollution-induced eutrophication, and dysoxia. Oxygen transport capacity of the hemolymph (blood equivalent) is considered the proximal driver of thermotolerance and respiration in many invertebrates. Moreover, maintaining homeostasis under environmental duress is inextricably linked to the activities of the hemolymph-based oxygen transport or binding proteins. Several protein groups fulfill this role in marine invertebrates: copper-based extracellular hemocyanins, iron-based intracellular hemoglobins and hemerythrins, and giant extracellular hemoglobins. In this brief text, we revisit the distribution and multifunctional properties of oxygen transport proteins, notably hemocyanins, in the context of climate change, and the consequent physiological reprogramming of marine invertebrates.

Biochar and Manure Applications Differentially Altered the Class 1 Integrons, Antimicrobial Resistance, and Gene Cassettes Diversity in Paddy Soils

Integrons are genetic components that are critically involved in bacterial evolution and antimicrobial resistance by assisting in the propagation and expression of gene cassettes. In recent decades, biochar has been introduced as a fertilizer to enhance physiochemical properties and crop yield of soil, while manure has been used as a fertilizer for centuries. The current study aimed to investigate the impact of biochar, manure, and a combination of biochar and manure on integrons, their gene cassettes, and relative antimicrobial resistance in paddy soil. Field experiments revealed class 1 (CL1) integrons were prevalent in all samples, with higher concentration and abundance in manure-treated plots than in biochar-treated ones. The gene cassette arrays in the paddy featured a broad pool of cassettes with a total of 35% novel gene cassettes. A majority of gene cassettes encoded resistance to aminoglycosides, heat shock protein, heavy metals, pilus secretory proteins, and twin-arginine translocases (Tat), TatA, TatB, and TatC. Both in combination and solo treatments, the diversity of gene cassettes was increased in the manure-enriched soil, however, biochar reduced the gene cassettes’ diversity and their cassettes array. Manure considerably enhanced CL1 integrons abundance and antimicrobial resistance, whereas biochar amendments significantly reduced integrons and antimicrobial resistance. The results highlighted the differential effects of biochar and manure on integrons and its gene cassette arrays, showing increased abundance of integrons and antibiotic resistance upon manure application and decrease of the same with biochar. The use of biochar alone or in combination with manure could be a beneficial alternative to mitigate the spread of antimicrobial resistance and bacterial evolution in the environment, specifically in paddy soils.

The hemerythrin-like diiron protein from Mycobacterium kansasii is a nitric oxide peroxidase

The hemerythrin-like protein from Mycobacterium kansasii (Mka HLP) is a member of a distinct class of oxo-bridged diiron proteins that are found only in mycobacterial species that cause respiratory disorders in humans. Because it had been shown to exhibit weak catalase activity and a change in absorbance on exposure to nitric oxide (NO), the reactivity of Mka HLP toward NO was examined under a variety of conditions. Under anaerobic conditions, we found that NO was converted to nitrite (NO₂⁻) via an intermediate, which absorbed light at 520 nm. Under aerobic conditions NO was converted to nitrate (NO₃⁻). In each of these two cases, the maximum amount of nitrite or nitrate formed was at best stoichiometric with the concentration of Mka HLP. When incubated with NO and H₂O₂, we observed NO peroxidase activity yielding nitrite and water as reaction products. Steady-state kinetic analysis of NO consumption during this reaction yielded a K_m for NO of 0.44 μM and a k_cat/K_m of 2.3 × 10⁵ M⁻¹s⁻¹. This high affinity for NO is consistent with a physiological role for Mka HLP in deterring nitrosative stress. This is the first example of a peroxidase that uses an oxo-bridged diiron center and a rare example of a peroxidase utilizing NO as an electron donor and cosubstrate. This activity provides a mechanism by which the infectious Mycobacterium may combat against the cocktail of NO and superoxide (O₂^•−) generated by macrophages to defend against bacteria, as well as to produce NO₂⁻ to adapt to hypoxic conditions.

Diversity of structures and functions of oxo-bridged non-heme diiron proteins

Oxo-bridged diiron proteins are a distinct class of non-heme iron proteins. Their active sites are composed of two irons that are coordinated by amino acid side chains, and a bridging oxygen that interacts with each iron. These proteins are members of the ferritin superfamily and share the structural feature of a four α-helix bundle that provides the residues that coordinate the irons. The different proteins also display a wide range of structures and functions. A prototype of this family is hemerythrin, which functions as an oxygen transporter. Several other hemerythrin-like proteins have been described with a diversity of functions including oxygen and iron sensing, and catalytic activities. Rubrerythrins react with hydrogen peroxide and rubrerythrin-like proteins possess a rubredoxin domain, in addition to the oxo-bridged diiron center. Other redox enzymes with oxo-bridged irons include flavodiiron proteins that act as O₂ or NO reductases, ribonucleotide reductase and methane monooxygenase. Ferritins have an oxo-bridged diiron in the ferroxidase center of the protein, which plays a role in the iron storage function of these proteins. There are also bacterial ferritins that exhibit catalytic activities. The structures and functions of this broad class of oxo-bridged diiron proteins are described and compared in this review.

Ferric iron reductases and their contribution to unicellular ferrous iron uptake

Iron is a necessary element for nearly all forms of life, and the ability to acquire this trace nutrient has been identified as a key virulence factor for the establishment of infection by unicellular pathogens. In the presence of O2, iron typically exists in the ferric (Fe3+) oxidation state, which is highly unstable in aqueous conditions, necessitating its sequestration into cofactors and/or host proteins to remain soluble. To counter this insolubility, and to compete with host sequestration mechanisms, many unicellular pathogens will secrete low molecular weight, high-affinity Fe3+ chelators known as siderophores. Once acquired, unicellular pathogens must liberate the siderophore-bound Fe3+ in order to assimilate this nutrient into metabolic pathways. While these organisms may hydrolyze the siderophore backbone to release the chelated Fe3+, this approach is energetically costly. Instead, iron may be liberated from the Fe3+-siderophore complex through reduction to Fe2+, which produces a lower-affinity form of iron that is highly soluble. This reduction is performed by a class of enzymes known as ferric reductases. Ferric reductases are broadly-distributed electron-transport proteins that are expressed by numerous infectious organisms and are connected to the virulence of unicellular pathogens. Despite this importance, ferric reductases remain poorly understood. This review provides an overview of our current understanding of unicellular ferric reductases (both soluble and membrane-bound), with an emphasis on the important but underappreciated connection between ferric-reductase mediated Fe3+ reduction and the transport of Fe2+ via ferrous iron transporters.

A Novel Marine Pathogen Isolated from Wild Cunners (Tautogolabrus adspersus): Comparative Genomics and Transcriptome Profiling of Pseudomonas sp. Strain J380

Cunner (Tautogolabrus adspersus) is a cleaner fish being considered for utilized in the North Atlantic salmon (Salmo salar) aquaculture industry to biocontrol sea lice infestations. However, bacterial diseases due to natural infections in wild cunners have yet to be described. This study reports the isolation of Pseudomonas sp. J380 from infected wild cunners and its phenotypic, genomic, and transcriptomic characterization. This Gram-negative motile rod-shaped bacterium showed a mesophilic (4-28 °C) and halotolerant growth. Under iron-limited conditions, Pseudomonas sp. J380 produced pyoverdine-type fluorescent siderophore. Koch's postulates were verified in wild cunners by intraperitoneally (i.p.) injecting Pseudomonas sp. J380 at 4 × 103, 4 × 105, and 4 × 107 colony forming units (CFU)/dose. Host-range and comparative virulence were also investigated in lumpfish and Atlantic salmon i.p. injected with ~106 CFU/dose. Lumpfish were more susceptible compared to cunners, and Atlantic salmon was resistant to Pseudomonas sp. J380 infection. Cunner tissues were heavily colonized by Pseudomonas sp. J380 compared to lumpfish and Atlantic salmon suggesting that it might be an opportunistic pathogen in cunners. The genome of Pseudomonas sp. J380 was 6.26 megabases (Mb) with a guanine-cytosine (GC) content of 59.7%. Biochemical profiles, as well as comparative and phylogenomic analyses, suggested that Pseudomonas sp. J380 belongs to the P. fluorescens species complex. Transcriptome profiling under iron-limited vs. iron-enriched conditions identified 1159 differentially expressed genes (DEGs). Cellular metabolic processes, such as ribosomal and energy production, and protein synthesis, were impeded by iron limitation. In contrast, genes involved in environmental adaptation mechanisms including two-component systems, histidine catabolism, and redox balance were transcriptionally up-regulated. Furthermore, iron limitation triggered the differential expression of genes encoding proteins associated with iron homeostasis. As the first report on a bacterial infection in cunners, the current study provides an overview of a new marine pathogen, Pseudomonas sp. J380.

Screening a Spliced Leader-Based Symbiodinium microadriaticum cDNA Library Using the Yeast-Two Hybrid System Reveals a Hemerythrin-Like Protein as a Putative SmicRACK1 Ligand

The dinoflagellate Symbiodiniaceae family plays a central role in the health of the coral reef ecosystem via the symbiosis that establishes with its inhabiting cnidarians and supports the host metabolism. In the last few decades, coral reefs have been threatened by pollution and rising temperatures which have led to coral loss. These events have raised interest in studying Symbiodiniaceae and their hosts; however, progress in understanding their metabolism, signal transduction pathways, and physiology in general, has been slow because dinoflagellates present peculiar characteristics. We took advantage of one of these peculiarities; namely, the post-transcriptional addition of a Dino Spliced Leader (Dino-SL) to the 5' end of the nuclear mRNAs, and used it to generate cDNA libraries from Symbiodinium microadriaticum. We compared sequences from two Yeast-Two Hybrid System cDNA Libraries, one based on the Dino-SL sequence, and the other based on the SMART technology (Switching Mechanism at 5' end of RNA Transcript) which exploits the template switching function of the reverse transcriptase. Upon comparison of the performance of both libraries, we obtained a significantly higher yield, number and length of sequences, number of transcripts, and better 5' representation from the Dino-SL based library than from the SMART library. In addition, we confirmed that the cDNAs from the Dino-SL library were adequately expressed in the yeast cells used for the Yeast-Two Hybrid System which resulted in successful screening for putative SmicRACK1 ligands, which yielded a putative hemerythrin-like protein.

Correlation of Conservation of Sequence and Structures of Mycobacterial Hemerythrin-like Proteins with Evolutionary Relationship and Host Pathogenicity

The Rv2633c gene of Mycobacterium tuberculosis, which plays a role in infection, encodes a hemerythrin-like protein (HLP). The crystal structure of an orthologue of Rv2633c, the HLP from Mycobacterium kansasii, revealed that it possessed structural features that were distinct from other hemerythrins and HLPs. These and other orthologous proteins comprise a distinct class of non-heme di-iron HLPs that are only found in mycobacteria. This study presents an analysis and comparison of protein sequences, putative structures, and evolutionary relationship of HLPs from 20 mycobacterial species that are known to cause tuberculosis or pulmonary disorders in humans. The results of this analysis allowed correlation of the physicochemical characteristics of amino acid residues that are substituted in these highly conserved sequences with their position in structures, possible effects on function, and evolutionary relationships. The sequences of the proteins from M. tuberculosis, Mycobacterium bovis, and other members of the M. tuberculosis complex, which cause tuberculosis, have substitutions not seen in the other non-tuberculous mycobacteria. Furthermore, groups of species that are closely related, based on phylogenetic analysis, possess substitutions of otherwise conserved residues not seen in other species that are less related. This information is correlated with the occurrence and clinical presentations of these groups of mycobacterial species. The results of this study provide a framework for structure-function studies to determine how subtle differences in the primary sequences of members of this family of proteins correlate with their structures and activities and how this may influence the infectious properties of the host species.

Crystal structure of a hemerythrin-like protein from Mycobacterium kansasii and homology model of the orthologous Rv2633c protein of M. tuberculosis

Pathogenic and opportunistic mycobacteria have a distinct class of non-heme di-iron hemerythrin-like proteins (HLPs). The first to be isolated was the Rv2633c protein, which plays a role in infection by Mycobacterium tuberculosis (Mtb), but could not be crystallized. This work presents the first crystal structure of an ortholog of Rv2633c, the mycobacterial HLP from Mycobacterium kansasii (Mka). This structure differs from those of hemerythrins and other known HLPs. It consists of five α-helices, whereas all other HLP domains have four. In contrast with other HLPs, the HLP domain is not fused to an additional protein domain. The residues ligating and surrounding the di-iron site are also unique among HLPs. Notably, a tyrosine occupies the position normally held by one of the histidine ligands in hemerythrin. This structure was used to construct a homology model of Rv2633c. The structure of five α-helices is conserved and the di-iron site ligands are identical in Rv2633c. Two residues near the ends of helices in the Mka HLP structure are replaced with prolines in the Rv2633c model. This may account for structural perturbations that decrease the solubility of Rv2633c relative to Mka HLP. Clusters of residues that differ in charge or polarity between Rv2633c and Mka HLP that point outward from the helical core could reflect a specificity for potential differential interactions with other protein partners in vivo, which are related to function. The Mka HLP exhibited weaker catalase activity than Rv2633c. Evidence was obtained for the interaction of Mka HLP irons with nitric oxide.

Pseudomonas aeruginosa lasR mutant fitness in microoxia is supported by an Anr-regulated oxygen-binding hemerythrin

Pseudomonas aeruginosa strains with loss-of-function mutations in the transcription factor LasR are frequently encountered in the clinic and the environment. Among the characteristics common to LasR-defective (LasR-) strains is increased activity of the transcription factor Anr, relative to their LasR+ counterparts, in low-oxygen conditions. One of the Anr-regulated genes found to be highly induced in LasR- strains was PA14_42860 (PA1673), which we named mhr for microoxic hemerythrin. Purified P. aeruginosa Mhr protein contained the predicted di-iron center and bound molecular oxygen with an apparent K _d of ∼1 µM. Both Anr and Mhr were necessary for fitness in lasR+ and lasR mutant strains in colony biofilms grown in microoxic conditions, and the effects were more striking in the lasR mutant. Among genes in the Anr regulon, mhr was most closely coregulated with the Anr-controlled high-affinity cytochrome c oxidase genes. In the absence of high-affinity cytochrome c oxidases, deletion of mhr no longer caused a fitness disadvantage, suggesting that Mhr works in concert with microoxic respiration. We demonstrate that Anr and Mhr contribute to LasR- strain fitness even in biofilms grown in normoxic conditions. Furthermore, metabolomics data indicate that, in a lasR mutant, expression of Anr-regulated mhr leads to differences in metabolism in cells grown on lysogeny broth or artificial sputum medium. We propose that increased Anr activity leads to higher levels of the oxygen-binding protein Mhr, which confers an advantage to lasR mutants in microoxic conditions.

Recent Insights into the Diversity and Evolution of Invertebrate Hemerythrins and Extracellular Globins

There are three broad groups of oxygen-transport proteins found in the haemolymph (blood) of invertebrates, namely the hemocyanins, the hemerythrins and the globins. Both hemerythrins and extracellular globins are iron-based proteins that are understudied when compared to the copper-containing hemocyanins. Recent evidence suggests that hemerythrins and (giant) extracellular globins (and their linker chains) are more widely distributed than previously thought and may have biological functions beyond oxygen transport and storage. Herein, we review contemporary literature of these often-neglected proteins with respect to their structural configurations on formation and ancestral states.

Gene refashioning through innovative shifting of reading frames in mosses

Early-diverging land plants such as mosses are known for their outstanding abilities to grow in various terrestrial habitats, incorporating tremendous structural and physiological innovations, as well as many lineage-specific genes. How these genes and functional innovations evolved remains unclear. In this study, we show that a dual-coding gene YAN/AltYAN in the moss Physcomitrella patens evolved from a pre-existing hemerythrin gene. Experimental evidence indicates that YAN/AltYAN is involved in fatty acid and lipid metabolism, as well as oil body and wax formation. Strikingly, both the recently evolved dual-coding YAN/AltYAN and the pre-existing hemerythrin gene might have similar physiological effects on oil body biogenesis and dehydration resistance. These findings bear important implications in understanding the mechanisms of gene origination and the strategies of plants to fine-tune their adaptation to various habitats.

Extant representatives of the earliest land plant lineages adapt to various terrestrial habitats with structural and physiological innovations. Here the authors show a dual-coding gene in the moss Physcomitrella patens evolved from a hemerythrin gene, with effects on oil body biogenesis and dehydration resistance.

Structure, function and evolution of the hemerythrin-like domain superfamily

Hemerythrin-like proteins have generally been studied for their ability to reversibly bind oxygen through their binuclear nonheme iron centers. However, in recent years, it has become increasingly evident that some members of the hemerythrin-like superfamily also participate in many other biological processes. For instance, the binuclear nonheme iron site of YtfE, a hemerythrin-like protein involved in the repair of iron centers in Escherichia coli, catalyzes the reduction of nitric oxide to nitrous oxide, and the human F-box/LRR-repeat protein 5, which contains a hemerythrin-like domain, is involved in intracellular iron homeostasis. Furthermore, structural data on hemerythrin-like domains from two proteins of unknown function, PF0695 from Pyrococcus furiosus and NMB1532 from Neisseria meningitidis, show that the cation-binding sites, typical of hemerythrin, can be absent or be occupied by metal ions other than iron. To systematically investigate this functional and structural diversity of the hemerythrin-like superfamily, we have collected hemerythrin-like sequences from a database comprising fully sequenced proteomes and generated a cluster map based on their all-against-all pairwise sequence similarity. Our results show that the hemerythrin-like superfamily comprises a large number of protein families which can be classified into three broad groups on the basis of their cation-coordinating residues: (a) signal-transduction and oxygen-carrier hemerythrins (H-HxxxE-HxxxH-HxxxxD); (b) hemerythrin-like (H-HxxxE-H-HxxxE); and, (c) metazoan F-box proteins (H-HExxE-H-HxxxE). Interestingly, all but two hemerythrin-like families exhibit internal sequence and structural symmetry, suggesting that a duplication event may have led to the origin of the hemerythrin domain.

Broad Phylogenetic Occurrence of the Oxygen-Binding Hemerythrins in Bilaterians

Animal tissues need to be properly oxygenated for carrying out catabolic respiration and, as such, natural selection has presumably favored special molecules that can reversibly bind and transport oxygen. Hemoglobins, hemocyanins, and hemerythrins (Hrs) fulfill this role, with Hrs being the least studied. Knowledge of oxygen-binding proteins is crucial for understanding animal physiology. Hr genes are present in the three domains of life, Archaea, Bacteria, and Eukaryota; however, within Animalia, Hrs has been reported only in marine species in six phyla (Annelida, Brachiopoda, Priapulida, Bryozoa, Cnidaria, and Arthropoda). Given this observed Hr distribution, whether all metazoan Hrs share a common origin is circumspect. We investigated Hr diversity and evolution in metazoans, by employing in silico approaches to survey for Hrs from of 120 metazoan transcriptomes and genomes. We found 58 candidate Hr genes actively transcribed in 36 species distributed in 11 animal phyla, with new records in Echinodermata, Hemichordata, Mollusca, Nemertea, Phoronida, and Platyhelminthes. Moreover, we found that “Hrs” reported from Cnidaria and Arthropoda were not consistent with that of other metazoan Hrs. Contrary to previous suggestions that Hr genes were absent in deuterostomes, we find Hr genes present in deuterostomes and were likely present in early bilaterians, but not in nonbilaterian animal lineages. As expected, the Hr gene tree did not mirror metazoan phylogeny, suggesting that Hrs evolutionary history was complex and besides the oxygen carrying capacity, the drivers of Hr evolution may also consist of secondary functional specializations of the proteins, like immunological functions.

Discovery and evolution of novel hemerythrin genes in annelid worms

BACKGROUND

Despite extensive study on hemoglobins and hemocyanins, little is known about hemerythrin (Hr) evolutionary history. Four subgroups of Hrs have been documented, including: circulating Hr (cHr), myohemerythrin (myoHr), ovohemerythrin (ovoHr), and neurohemerythrin (nHr). Annelids have the greatest diversity of oxygen carrying proteins among animals and are the only phylum in which all Hr subgroups have been documented. To examine Hr diversity in annelids and to further understand evolution of Hrs, we employed approaches to survey annelid transcriptomes in silico.

RESULTS

Sequences of 214 putative Hr genes were identified from 44 annelid species in 40 different families and Bayesian inference revealed two major clades with strong statistical support. Notably, the topology of the Hr gene tree did not mirror the phylogeny of Annelida as presently understood, and we found evidence of extensive Hr gene duplication and loss in annelids. Gene tree topology supported monophyly of cHrs and a myoHr clade that included nHrs sequences, indicating these designations are functional rather than evolutionary.

CONCLUSIONS

The presence of several cHrs in early branching taxa suggests that a variety of Hrs were present in the common ancestor of extant annelids. Although our analysis was limited to expressed-coding regions, our findings demonstrate a greater diversity of Hrs among annelids than previously reported.