Differential loss and retention of cytoglobin, myoglobin, and globin-E during the radiation of vertebrates

Molecular evolution of the hemoglobin gene family across vertebrates

Adaptation to various altitudes and oxygen levels is a major aspect of vertebrate evolution. Hemoglobin is an erythrocyte protein belonging to the globin superfamily, and the α-, β-globin genes of jawed vertebrates encode tetrameric ((α₂β₂) hemoglobin, which contributes to aerobic metabolism by delivering oxygen from the respiratory exchange surfaces into cells. However, there are various gaps in knowledge regarding hemoglobin gene evolution, including patterns in cartilaginous fish and the roles of gene conversion in various taxa. Hence, we evaluated the evolutionary history of the vertebrate hemoglobin gene family by analyses of 97 species representing all classes of vertebrates. By genome-wide analyses, we extracted 879 hemoglobin sequences. Members of the hemoglobin gene family were conserved in birds and reptiles but variable in mammals, amphibians, and teleosts. Gene motifs, structures, and synteny were relatively well-conserved among vertebrates. Our results revealed that purifying selection contributed substantially to the evolution of all vertebrate hemoglobin genes, with mean d_N/d_S (ω) values ranging from 0.057 in teleosts to 0.359 in reptiles. In general, after the fish-specific genome duplication, the teleost hemoglobin genes showed variation in rates of evolution, and the β-globin genes showed relatively high ω values after a gene transposition event in amniotes. We also observed that the frequency of gene conversion was high in amniotes, with fewer hemoglobin genes and higher rates of evolution. Collectively, our findings provide detail insight into complex evolutionary processes shaping the vertebrate hemoglobin gene family, involving gene duplication, gene loss, purifying selection, and gene conversion.

Structural evolution of an amphibian-specific globin: A computational evolutionary biochemistry approach

Studies on the globin family are continuously revealing insights into the mechanisms of gene and protein evolution. The rise of a new globin gene type in Pelobatoidea and Neobatrachia (Amphibia:Anura) from an α-globin precursor provides the opportunity to investigate the genetic and physical mechanisms underlying the origin of new protein structural and functional properties. This amphibian-specific globin (globin A/GbA) discovered in the heart of Rana catesbeiana is a monomer. As the ancestral oligomeric state of α-globins is a homodimer, we inferred that the ancestral state was lost somewhere in the GbA lineage. Here, we combined computational molecular evolution with structural bioinformatics to determine the extent to which the loss of the homodimeric state is pervasive in the GbA clade. We also characterized the loci of GbA genes in Bufo bufo. We found two GbA clades in Neobatrachia. One was deleted in Ranidae, but retained and expanded to yield a new globin cluster in Bufonidae species. Loss of the ancestral oligomeric state seems to be pervasive in the GbA clade. However, a taxonomic sampling that includes more Pelobatoidea, as well as early Neobatrachia, lineages would be necessary to determine the oligomeric state of the last common ancestor of all GbA. The evidence presented here points out a possible loss of oligomerization in Pelobatoidea GbA as a result of amino acid substitutions that weaken the homodimeric state. In contrast, the loss of oligomerization in both Neobatrachia GbA clades was linked to independent deletions that disrupted many packing contacts at the homodimer interface.

The role of globins in cardiovascular physiology

Globin proteins exist in every cell type of the vasculature, from erythrocytes to endothelial cells, vascular smooth muscle cells, and peripheral nerve cells. Many globin subtypes are also expressed in muscle tissues (including cardiac and skeletal muscle), in other organ-specific cell types, and in cells of the central nervous system (CNS). The ability of each of these globins to interact with molecular oxygen (O2) and nitric oxide (NO) is preserved across these contexts. Endothelial α-globin is an example of extraerythrocytic globin expression. Other globins, including myoglobin, cytoglobin, and neuroglobin, are observed in other vascular tissues. Myoglobin is observed primarily in skeletal muscle and smooth muscle cells surrounding the aorta or other large arteries. Cytoglobin is found in vascular smooth muscle but can also be expressed in nonvascular cell types, especially in oxidative stress conditions after ischemic insult. Neuroglobin was first observed in neuronal cells, and its expression appears to be restricted mainly to the CNS and the peripheral nervous system. Brain and CNS neurons expressing neuroglobin are positioned close to many arteries within the brain parenchyma and can control smooth muscle contraction and thus tissue perfusion and vascular reactivity. Overall, reactions between NO and globin heme iron contribute to vascular homeostasis by regulating vasodilatory NO signals and scavenging reactive species in cells of the mammalian vascular system. Here, we discuss how globin proteins affect vascular physiology, with a focus on NO biology, and offer perspectives for future study of these functions.

Whole-Genome Duplications and the Diversification of the Globin-X Genes of Vertebrates

Globin-X (GbX) is an enigmatic member of the vertebrate globin gene family with a wide phyletic distribution that spans protostomes and deuterostomes. Unlike canonical globins such as hemoglobins and myoglobins, functional data suggest that GbX does not have a primary respiratory function. Instead, evidence suggests that the monomeric, membrane-bound GbX may play a role in cellular signaling or protection against the oxidation of membrane lipids. Recently released genomes from key vertebrates provide an excellent opportunity to address questions about the early stages of the evolution of GbX in vertebrates. We integrate bioinformatics, synteny, and phylogenetic analyses to characterize the diversity of GbX genes in nonteleost ray-finned fishes, resolve relationships between the GbX genes of cartilaginous fish and bony vertebrates, and demonstrate that the GbX genes of cyclostomes and gnathostomes derive from independent duplications. Our study highlights the role that whole-genome duplications (WGDs) have played in expanding the repertoire of genes in vertebrate genomes. Our results indicate that GbX paralogs have a remarkably high rate of retention following WGDs relative to other globin genes and provide an evolutionary framework for interpreting results of experiments that examine functional properties of GbX and patterns of tissue-specific expression. By identifying GbX paralogs that are products of different WGDs, our results can guide the design of experimental work to explore whether gene duplicates that originate via WGDs have evolved novel functional properties or expression profiles relative to singleton or tandemly duplicated copies of GbX.

The More, the Merrier? Multiple Myoglobin Genes in Fish Species, Especially in Gray Bichir (Polypterus senegalus) and Reedfish (Erpetoichthys calabaricus)

The members of the globin superfamily are a classical model system to investigate gene evolution and their fates as well as the diversity of protein function. One of the best-known globins is myoglobin (Mb), which is mainly expressed in heart muscle and transports oxygen from the sarcolemma to the mitochondria. Most vertebrates harbor a single copy of the myoglobin gene, but some fish species have multiple myoglobin genes. Phylogenetic analyses indicate an independent emergence of multiple myoglobin genes, whereby the origin is mostly the last common ancestor of each order. By analyzing different transcriptome data sets, we found at least 15 multiple myoglobin genes in the polypterid gray bichir (Polypterus senegalus) and reedfish (Erpetoichthys calabaricus). In reedfish, the myoglobin genes are expressed in a broad range of tissues but show very different expression values. In contrast, the Mb genes of the gray bichir show a rather scattered expression pattern; only a few Mb genes were found expressed in the analyzed tissues. Both, gray bichir and reedfish possess lungs which enable them to inhabit shallow and swampy waters throughout tropical Africa with frequently fluctuating and low oxygen concentrations. The myoglobin repertoire probably reflects the molecular adaptation to these conditions. The sequence divergence, the substitution rate, and the different expression pattern of multiple myoglobin genes in gray bichir and reedfish imply different functions, probably through sub- and neofunctionalization during evolution.

The rise and fall of globins in the amphibia

The globin gene repertoire of gnathostome vertebrates is dictated by differential retention and loss of nine paralogous genes: androglobin, neuroglobin, globin X, cytoglobin, globin Y, myoglobin, globin E, and the α- and β-globins. We report the globin gene repertoire of three orders of modern amphibians: Anura, Caudata, and Gymnophiona. Combining phylogenetic and conserved synteny analysis, we show that myoglobin and globin E were lost only in the Batrachia clade, but retained in Gymnophiona. The major amphibian groups also retained different paralogous copies of globin X. None of the amphibian presented α^D-globin gene. Nevertheless, two clades of β-globins are present in all amphibians, indicating that the amphibian ancestor possessed two paralogous proto β-globins. We also show that orthologs of the gene coding for the monomeric hemoglobin found in the heart of Rana catesbeiana are present in Neobatrachia and Pelobatoidea species we analyzed. We suggest that these genes might perform myoglobin- and globin E-related functions. We conclude that the repertoire of globin genes in amphibians is dictated by both retention and loss of the paralogous genes cited above and the rise of a new globin gene through co-option of an α-globin, possibly facilitated by a prior event of transposition.

Evolutionary History of the Globin Gene Family in Annelids

Animals depend on the sequential oxidation of organic molecules to survive; thus, oxygen-carrying/transporting proteins play a fundamental role in aerobic metabolism. Globins are the most common and widespread group of respiratory proteins. They can be divided into three types: circulating intracellular, noncirculating intracellular, and extracellular, all of which have been reported in annelids. The diversity of oxygen transport proteins has been underestimated across metazoans. We probed 250 annelid transcriptomes in search of globin diversity in order to elucidate the evolutionary history of this gene family within this phylum. We report two new globin types in annelids, namely androglobins and cytoglobins. Although cytoglobins and myoglobins from vertebrates and from invertebrates are referred to by the same name, our data show they are not genuine orthologs. Our phylogenetic analyses show that extracellular globins from annelids are more closely related to extracellular globins from other metazoans than to the intracellular globins of annelids. Broadly, our findings indicate that multiple gene duplication and neo-functionalization events shaped the evolutionary history of the globin family.

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 O₂, regulating taxis away or towards high O₂ 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.

Globins and nitric oxide homeostasis in fish embryonic development

Since the discovery of new members of the globin superfamily such as Cytoglobin, Neuroglobin and Globin X, in addition to the most well-known members, Hemoglobin and Myoglobin, different hypotheses have been suggested about their function in vertebrates. Globins are ubiquitously found in living organisms and can carry out different functions based on their ability to bind ligands such as O2, and nitric oxide (NO) and to catalyze reactions scavenging NO or generating NO by reducing nitrite. NO is a highly diffusible molecule with a central role in signaling important for egg maturation, fertilization and early embryonic development. The globins ability to scavenge or generate NO makes these proteins ideal candidates in regulating NO homeostasis depending on the micro environment and tissue NO demands. Different amounts of various globins have been found in zebrafish eggs and developing embryos where it's unlikely that they function as respiratory proteins and instead could play a role in maintaining embryonic NO homeostasis. Here we summarize the current knowledge concerning the role of NO in adult fish in comparison to mammals and we discuss NO function during embryonic development with possible implications for globins in maintaining embryonic NO homeostasis.

Adaptations to environmental change: Globin superfamily evolution in Antarctic fishes

The ancient origins and functional versatility of globins make them ideal subjects for studying physiological adaptation to environmental change. Our goals in this review are to describe the evolution of the vertebrate globin gene superfamily and to explore the structure/function relationships of hemoglobin, myoglobin, neuroglobin and cytoglobin in teleost fishes. We focus on the globins of Antarctic notothenioids, emphasizing their adaptive features as inferred from comparisons with human proteins. We dedicate this review to Guido di Prisco, our co-author, colleague, friend, and husband of C.V. Ever thoughtful, creative, and enthusiastic, Guido spearheaded study of the structure, function, and evolution of the hemoglobins of polar fishes - this review is testimony to his wide-ranging contributions. Throughout his career, Guido inspired younger scientists to embrace polar biological research, and he challenged researchers of all ages to explore evolutionary adaptation in the context of global climate change. Beyond his scientific contributions, we will miss his warmth, his culture, and his great intellect. Guido has left an outstanding legacy, one that will continue to inspire us and our research.

Genetic and functional diversity of the multiple lungfish myoglobins

It is known that the West African lungfish (Protopterus annectens) harbours multiple myoglobin (Mb) genes that are differentially expressed in various tissues and that the Mbs differ in their abilities to confer tolerance towards hypoxia. Here, we show that other lungfish species (Protopterus dolloi, Protopterus aethiopicus and Lepidosiren paradoxa) display a similar diversity of Mb genes and have orthologous Mb genes. To investigate the functional diversification of these genes, we studied the structures, O₂ binding properties and nitrite reductase enzymatic activities of recombinantly expressed P. annectens Mbs (PanMbs). CD spectroscopy and small-angle X-ray scattering revealed the typical globin-fold in all investigated recombinant Mbs, indicating a conserved structure. The highest O₂ affinity was measured for PanMb2 (P₅₀  = 0.88 Torr at 20 °C), which is mainly expressed in the brain, whereas the muscle-specific PanMb1 has the lowest O₂ affinity (P₅₀  = 3.78 Torr at 20 °C), suggesting that tissue-specific O₂ requirements have resulted in the emergence of distinct Mb types. Two of the mainly neuronally expressed Mbs (PanMb3 and PanMb4b) have the highest nitrite reductase rates. These data show different O₂ binding and enzymatic properties of lungfish Mbs, reflecting multiple subfunctionalisation and neofunctionalisation events that occurred early in the evolution of lungfish. Some Mbs may have also taken over the functions of neuroglobin and cytoglobin, which are widely expressed in vertebrates but appear to be missing in lungfish.

The Zebrafish Cytochrome b5/Cytochrome b5 Reductase/NADH System Efficiently Reduces Cytoglobins 1 and 2: Conserved Activity of Cytochrome b5/Cytochrome b5 Reductases during Vertebrate Evolution

Cytoglobin is a heme protein evolutionarily related to hemoglobin and myoglobin. Cytoglobin is expressed ubiquitously in mammalian tissues; however, its physiological functions are yet unclear. Phylogenetic analyses indicate that the cytoglobin gene is highly conserved in vertebrate clades, from fish to reptiles, amphibians, birds, and mammals. Most proposed roles for cytoglobin require the maintenance of a pool of reduced cytoglobin (Fe^II). We have shown previously that the human cytochrome b₅/cytochrome b₅ reductase system, considered a quintessential hemoglobin/myoglobin reductant, can reduce human and zebrafish cytoglobins ≤250-fold faster than human hemoglobin or myoglobin. It was unclear whether this reduction of zebrafish cytoglobins by mammalian proteins indicates a conserved pathway through vertebrate evolution. Here, we report the reduction of zebrafish cytoglobins 1 and 2 by the zebrafish cytochrome b₅ reductase and the two zebrafish cytochrome b₅ isoforms. In addition, the reducing system also supports reduction of Globin X, a conserved globin in fish and amphibians. Indeed, the zebrafish reducing system can maintain a fully reduced pool for both cytoglobins, and both cytochrome b₅ isoforms can support this process. We determined the P₅₀ for oxygen to be 0.5 Torr for cytoglobin 1 and 4.4 Torr for cytoglobin 2 at 25 °C. Thus, even at low oxygen tensions, the reduced cytoglobins may exist in a predominant oxygen-bound form. Under these conditions, the cytochrome b₅/cytochrome b₅ reductase system can support a conserved role for cytoglobins through evolution, providing electrons for redox signaling reactions such as nitric oxide dioxygenation, nitrite reduction, and phospholipid oxidation.

Myoglobin as a versatile peroxidase: Implications for a more important role for vertebrate striated muscle in antioxidant defense

Myoglobins (Mb) are ubiquitous proteins found in striated muscle of nearly all vertebrate taxa. Although their function is most commonly associated with facilitating oxygen storage and diffusion, Mb has also been implicated in cellular antioxidant defense. The oxidized (Fe³⁺) form of Mb (metMB) can react with hydrogen peroxide (H₂O₂) to produce ferrylMb. FerrylMb can be reduced back to metMb for another round of reaction with H₂O₂. In the present study, we have shown that horse skeletal muscle Mb displays peroxidase activity using 2,2'-azino-di-(3-ethylbenzothiazoline)-6-sulfonic acid (ABTS) and 3,3',5,5'-tetramethylbenzidine (TMB) as reducing substrates, as well as the biologically-relevant substrates NADH/NADPH, ascorbate, caffeic acid, and resveratrol. We have also shown that ferrylMb can be reduced by both ethanol and acetaldehyde, which are known to accumulate in some vertebrate tissues under anaerobic conditions, such as anoxic goldfish and crucian carp, implying a potential mechanism for ethanol detoxification in striated muscle. We found that metMb peroxidase activity is pH-dependent, increasing as pH decreases from 7.4 to 6.1, which is biologically relevant to anaerobic vertebrate muscle when incurring intracellular lactic acidosis. Finally, we found that metMb reacts with hypochlorite in a heme-dependent fashion, indicating that Mb could play a role in hypochlorite detoxification. Taken together, these data suggest that Mb peroxidase activity might be an important antioxidant mechanism in vertebrate cardiac and skeletal muscle under a variety of physiological conditions, such as those that might occur in contracting skeletal muscle or during hypoxia.

Globin E is a myoglobin-related, respiratory protein highly expressed in lungfish oocytes

Globins are a classical model system for the studies of protein evolution and function. Recent studies have shown that – besides the well-known haemoglobin and myoglobin – additional globin-types occur in vertebrates that serve different functions. Globin E (GbE) was originally identified as an eye-specific protein of birds that is distantly related to myoglobin. GbE is also present in turtles and the coelacanth but appeared to have been lost in other vertebrates. Here, we show that GbE additionally occurs in lungfish, the closest living relatives of the tetrapods. Each lungfish species harbours multiple (≥5) GbE gene copies. Surprisingly, GbE is exclusively and highly expressed in oocytes, with mRNA levels that exceed that of myoglobin in the heart. Thus, GbE is the first known oocyte-specific globin in vertebrates. No GbE transcripts were found in the ovary or egg transcriptomes of other vertebrates, suggesting a lungfish-specific function. Spectroscopic analysis and kinetic studies of recombinant GbE1 of the South American lungfish Lepidosiren paradoxa revealed a typical pentacoordinate globin with myoglobin-like O₂-binding kinetics, indicating similar functions. Our findings suggest that the multiple copies of GbE evolved to enhance O₂-supply in the developing embryo of lungfish, analogous to the embryonic and fetal haemoglobins of other vertebrates. In evolution, GbE must have changed its expression site from oocytes to eyes, or vice versa.

Enhancing Statistical Multiple Sequence Alignment and Tree Inference Using Structural Information

For highly divergent sequences, there is often insufficient information to reliably construct alignments and phylogenetic trees. Since protein structure may be strongly conserved despite large divergences in sequence, structural information can be used to help identify homology in such cases.While there exist well-studied models of sequence evolution, structurally informed alignment methods have typically made use of geometric measures of deviation that do not take into account the underlying mutational processes. In order to integrate structural information into sequence-based evolutionary models, we recently developed a stochastic model of structural evolution on a phylogenetic tree and implemented this as the StructAlign plugin for the StatAlign statistical alignment package.In this chapter, we will outline the types of analyses that can be carried out using StructAlign, illustrating how the inclusion of structural information can be used to inform joint estimation of alignments and trees. StructAlign can also be used to infer branch-specific rates of structural evolution, and analysis of an example globin dataset highlights strong variation in the inferred rate across the tree. While structure is more highly conserved within clades, the rate of structural divergence as a function of sequence variation is larger between functionally divergent proteins. Allowing for the rate of structural divergence to vary over the tree results in an improved fit to the empirically observed pairwise RMSD values.

Genome sequence of walking catfish (Clarias batrachus) provides insights into terrestrial adaptation

Background

Walking catfish (Clarias batrachus) is a freshwater fish capable of air-breathing and locomotion on land. It usually inhabits various low-oxygen habitats, burrows inside the mudflat, and sometimes “walks” to search for suitable environments during summer. It has evolved accessory air-breathing organs for respiring air and corresponding mechanisms to survive in such challenging environments. Thereby, it serves as a great model for understanding adaptations to terrestrial life.

Results

Comparative genomics with channel catfish (Ictalurus punctatus) revealed specific adaptations of C. batrachus in DNA repair, enzyme activator activity, and small GTPase regulator activity. Comparative analysis with 11 non-air-breathing fish species suggested adaptive evolution in gene expression and nitrogenous waste metabolic processes. Further, myoglobin, olfactory receptor related to class A G protein-coupled receptor 1, and sulfotransferase 6b1 genes were found to be expanded in the air-breathing walking catfish genome, with 15, 15, and 12 copies, respectively, compared to non-air-breathing fishes that possess only 1–2 copies of these genes. Additionally, we sequenced and compared the transcriptomes of the gill and the air-breathing organ to characterize the mechanism of aerial respiration involved in elastic fiber formation, oxygen binding and transport, angiogenesis, ion homeostasis and acid-base balance. The hemoglobin genes were expressed dramatically higher in the air-breathing organ than in the gill of walking catfish.

Conclusions

This study provides an important genomic resource for understanding the adaptive mechanisms of walking catfish to terrestrial environments. It is possible that the coupling of enhanced abilities for oxygen storage and oxygen transport through genomic expansion of myoglobin genes and transcriptomic up-regulation of hemoglobin and angiogenesis-related genes are important components of the molecular basis for adaptation of this aquatic species to terrestrial life.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5355-9) contains supplementary material, which is available to authorized users.

The Reprimo Gene Family: A Novel Gene Lineage in Gastric Cancer with Tumor Suppressive Properties

The reprimo (RPRM) gene family is a group of single exon genes present exclusively within the vertebrate lineage. Two out of three members of this family are present in humans: RPRM and RPRM-Like (RPRML). RPRM induces cell cycle arrest at G2/M in response to p53 expression. Loss-of-expression of RPRM is related to increased cell proliferation and growth in gastric cancer. This evidence suggests that RPRM has tumor suppressive properties. However, the molecular mechanisms and signaling partners by which RPRM exerts its functions remain unknown. Moreover, scarce studies have attempted to characterize RPRML, and its functionality is unclear. Herein, we highlight the role of the RPRM gene family in gastric carcinogenesis, as well as its potential applications in clinical settings. In addition, we summarize the current knowledge on the phylogeny and expression patterns of this family of genes in embryonic zebrafish and adult humans. Strikingly, in both species, RPRM is expressed primarily in the digestive tract, blood vessels and central nervous system, supporting the use of zebrafish for further functional characterization of RPRM. Finally, drawing on embryonic and adult expression patterns, we address the potential relevance of RPRM and RPRML in cancer. Active investigation or analytical research in the coming years should contribute to novel translational applications of this poorly understood gene family as potential biomarkers and development of novel cancer therapies.

Evolution and Expression of Tissue Globins in Ray-Finned Fishes

The globin gene family encodes oxygen-binding hemeproteins conserved across the major branches of multicellular life. The origins and evolutionary histories of complete globin repertoires have been established for many vertebrates, but there remain major knowledge gaps for ray-finned fish. Therefore, we used phylogenetic, comparative genomic and gene expression analyses to discover and characterize canonical “non-blood” globin family members (i.e., myoglobin, cytoglobin, neuroglobin, globin-X, and globin-Y) across multiple ray-finned fish lineages, revealing novel gene duplicates (paralogs) conserved from whole genome duplication (WGD) and small-scale duplication events. Our key findings were that: (1) globin-X paralogs in teleosts have been retained from the teleost-specific WGD, (2) functional paralogs of cytoglobin, neuroglobin, and globin-X, but not myoglobin, have been conserved from the salmonid-specific WGD, (3) triplicate lineage-specific myoglobin paralogs are conserved in arowanas (Osteoglossiformes), which arose by tandem duplication and diverged under positive selection, (4) globin-Y is retained in multiple early branching fish lineages that diverged before teleosts, and (5) marked variation in tissue-specific expression of globin gene repertoires exists across ray-finned fish evolution, including several previously uncharacterized sites of expression. In this respect, our data provide an interesting link between myoglobin expression and the evolution of air breathing in teleosts. Together, our findings demonstrate great-unrecognized diversity in the repertoire and expression of nonblood globins that has arisen during ray-finned fish evolution.

Efficient Reduction of Vertebrate Cytoglobins by the Cytochrome b5/Cytochrome b5 Reductase/NADH System

Cytoglobin is a heme-containing protein ubiquitous in mammalian tissues. Unlike the evolutionarily related proteins hemoglobin and myoglobin, cytoglobin shows a six-coordinated heme binding, with the heme iron coordinated by two histidine side chains. Cytoglobin is involved in cytoprotection pathways through yet undefined mechanisms, and it has recently been demonstrated that cytoglobin has redox signaling properties via nitric oxide (NO) and nitrite metabolism. The reduced, ferrous cytoglobin can bind oxygen and will react with NO in a dioxygenation reaction to form nitrate, which dampens NO signaling. When deoxygenated, cytoglobin can bind nitrite and reduce it to NO. This oxidoreductase activity could be catalytic if an effective reduction system exists to regenerate the reduced heme species. The nature of the physiological cytoglobin reducing system is unknown, although it has been proposed that ascorbate and cytochrome b₅ could fulfill this role. Here we describe that physiological concentrations of cytochrome b₅ and cytochrome b₅ reductase can reduce human and fish cytoglobins at rates up to 250-fold higher than those reported for their known physiological substrates, hemoglobin and myoglobin, and up to 100-fold faster than 5 mM ascorbate. These data suggest that the cytochrome b₅/cytochrome b₅ reductase system is a viable reductant for cytoglobin in vivo, allowing for catalytic oxidoreductase activity.

Stability-Mediated Epistasis Restricts Accessible Mutational Pathways in the Functional Evolution of Avian Hemoglobin

If the fitness effects of amino acid mutations are conditional on genetic background, then mutations can have different effects depending on the sequential order in which they occur during evolutionary transitions in protein function. A key question concerns the fraction of possible mutational pathways connecting alternative functional states that involve transient reductions in fitness. Here we examine the functional effects of multiple amino acid substitutions that contributed to an evolutionary transition in the oxygenation properties of avian hemoglobin (Hb). The set of causative changes included mutations at intradimer interfaces of the Hb tetramer. Replacements at such sites may be especially likely to have epistatic effects on Hb function since residues at intersubunit interfaces are enmeshed in networks of salt bridges and hydrogen bonds between like and unlike subunits; mutational reconfigurations of these atomic contacts can affect allosteric transitions in quaternary structure and the propensity for tetramer-dimer dissociation. We used ancestral protein resurrection in conjunction with a combinatorial protein engineering approach to synthesize genotypes representing the complete set of mutational intermediates in all possible forward pathways that connect functionally distinct ancestral and descendent genotypes. The experiments revealed that 1/2 of all possible forward pathways included mutational intermediates with aberrant functional properties because particular combinations of mutations promoted tetramer-dimer dissociation. The subset of mutational pathways with unstable intermediates may be selectively inaccessible, representing evolutionary roads not taken. The experimental results also demonstrate how epistasis for particular functional properties of proteins may be mediated indirectly by mutational effects on quaternary structural stability.

Predictable convergence in hemoglobin function has unpredictable molecular underpinnings

To investigate the predictability of genetic adaptation, we examined the molecular basis of convergence in hemoglobin function in comparisons involving 56 avian taxa that have contrasting altitudinal range limits. Convergent increases in hemoglobin-oxygen affinity were pervasive among high-altitude taxa, but few such changes were attributable to parallel amino acid substitutions at key residues. Thus, predictable changes in biochemical phenotype do not have a predictable molecular basis. Experiments involving resurrected ancestral proteins revealed that historical substitutions have context-dependent effects, indicating that possible adaptive solutions are contingent on prior history. Mutations that produce an adaptive change in one species may represent precluded possibilities in other species because of differences in genetic background.

Gene Duplication and Evolutionary Innovations in Hemoglobin-Oxygen Transport

During vertebrate evolution, duplicated hemoglobin (Hb) genes diverged with respect to functional properties as well as the developmental timing of expression. For example, the subfamilies of genes that encode the different subunit chains of Hb are ontogenetically regulated such that functionally distinct Hb isoforms are expressed during different developmental stages. In some vertebrate taxa, functional differentiation between co-expressed Hb isoforms may also contribute to physiologically important divisions of labor.

The Phenotypic Plasticity of Duplicated Genes in Saccharomyces cerevisiae and the Origin of Adaptations

Gene and genome duplication are the major sources of biological innovations in plants and animals. Functional and transcriptional divergence between the copies after gene duplication has been considered the main driver of innovations . However, here we show that increased phenotypic plasticity after duplication plays a more major role than thought before in the origin of adaptations. We perform an exhaustive analysis of the transcriptional alterations of duplicated genes in the unicellular eukaryote Saccharomyces cerevisiae when challenged with five different environmental stresses. Analysis of the transcriptomes of yeast shows that gene duplication increases the transcriptional response to environmental changes, with duplicated genes exhibiting signatures of adaptive transcriptional patterns in response to stress. The mechanism of duplication matters, with whole-genome duplicates being more transcriptionally altered than small-scale duplicates. The predominant transcriptional pattern follows the classic theory of evolution by gene duplication; with one gene copy remaining unaltered under stress, while its sister copy presents large transcriptional plasticity and a prominent role in adaptation. Moreover, we find additional transcriptional profiles that are suggestive of neo- and subfunctionalization of duplicate gene copies. These patterns are strongly correlated with the functional dependencies and sequence divergence profiles of gene copies. We show that, unlike singletons, duplicates respond more specifically to stress, supporting the role of natural selection in the transcriptional plasticity of duplicates. Our results reveal the underlying transcriptional complexity of duplicated genes and its role in the origin of adaptations.

Genomic conservation of erythropoietic microRNAs (erythromiRs) in white-blooded Antarctic icefish

White-blooded Antarctic crocodile icefish are the only vertebrates known to lack functional hemoglobin genes and red blood cells throughout their lives. We do not yet know, however, whether extinction of hemoglobin genes preceded loss of red blood cells or vice versa, nor whether erythropoiesis regulators disappeared along with hemoglobin genes in this erythrocyte-null clade. Several microRNAs, which we here call erythromiRs, are expressed primarily in developing red blood cells in zebrafish, mouse, and humans. Abrogating some erythromiRs, like mir144 and mir451a, leads to profound anemia, demonstrating a functional role in erythropoiesis. Here, we tested two not mutually exclusive hypotheses: 1) that the loss of one or more erythromiR genes extinguished the erythropoietic program of icefish and/or led to the loss of globin gene expression through pseudogenization; and 2) that some erythromiR genes were secondarily lost after the loss of functional hemoglobin and red blood cells in icefish. We explored small RNA transcriptomes generated from the hematopoietic kidney marrow of four Antarctic notothenioids: two red-blooded species (bullhead notothen Notothenia coriiceps and emerald notothen Trematomus bernacchii) and two white-blooded icefish (blackfin icefish Chaenocephalus aceratus and hooknose icefish Chionodraco hamatus). The N. coriiceps genome assembly anchored analyses. Results showed that, like the two red-blooded species, the blackfin icefish genome possessed and the marrow expressed all known erythromiRs. This result indicates that loss of hemoglobin and red blood cells in icefish was not caused by loss of known erythromiR genes. Furthermore, expression of only one erythromiR, mir96, appears to have been lost after the loss of red blood cells and hemoglobin-expression was not detected in the erythropoietic organ of hooknose icefish but was present in blackfin icefish. All other erythromiRs investigated, including mir144 and mir451a, were expressed by all four species and thus are present in the genomes of at least the two white-blooded icefish. Our results rule out the hypothesis that genomic loss of any known erythromiRs extinguished erythropoiesis in icefish, and suggest that after the loss of red blood cells, few erythromiRs experienced secondary loss. Results suggest that functions independent of erythropoiesis maintained erythromiRs, thereby highlighting the evolutionary resilience of miRNA genes in vertebrate genomes.

Insight into the phylogenetic relationship and structural features of vertebrate myoglobin family

Myoglobin (Mb) is studied to clarify the structure-function relationships in protein science. In this work, we report the results of a comparative analysis of amino acid sequences from 298 vertebrate Mbs. Forty-one high conserved residues were identified and seven of them were invariants [E18, G25, F43, V68, L72, H93 (proximal histidine) and H97]. E18 is the only invariant amino acid residue located out of the heme-pocket and Xe-cavities playing a role in interaction between the A and E-helices. A comparative analysis of several parameters related to amino acid composition shows an increase of average mass, accessible surface area and volume per residue from Actinopterygii to Mammalia and Aves. This may be due to an increased number of bulky residues reducing the non-specific cavities volume and thus improving the oxygen flow between the heme site and the outside of the protein. Finally, the phylogenetic analyses of Mb in vertebrates are consistent with an evolution that runs with the diversification of the species, but in which several episodes of gene duplication and lost have occurred, less frequently in the ancestors of great taxons, cartilaginous fishes and non-avian reptiles, most frequently in ray-finned fishes and mammals, and very frequently in birds.

Characterization of zebrafish neuroglobin and cytoglobins 1 and 2: Zebrafish cytoglobins provide insights into the transition from six-coordinate to five-coordinate globins

Neuroglobin (Ngb) and cytoglobin (Cygb) are two six-coordinate heme proteins of unknown physiological function. Although studies on the mammalian proteins have elucidated aspects of Ngb and Cygb biophysics and indicated potential functions, the properties of non-mammalian Ngbs and Cygbs are largely uncharacterized. We have expressed the recombinant zebrafish proteins Ngb, Cygb1, and Cygb2 in Escherichia coli and characterized their nitrite reduction rates, spectral properties, autoxidation rate constants, redox potentials and lipid binding properties. The three zebrafish proteins can catalyze the reduction of nitrite to nitric oxide with a broad range of reaction rate constants. (Ngb, 0.68 ± 0.04 M(-1) s(-1); Cygb1, 28.6 ± 3.1 M(-1) s(-1); Cygb2, 0.94 ± 0.18 M(-1) s(-1)). We observe that zebrafish Ngb and Cygb2 have comparable spectral features to those of human Ngb and Cygb, consistent with a six-coordinate heme, whereas unexpectedly Cygb1 has a five-coordinate heme, a slower autoxidation and in general has properties more akin to oxygen transport proteins. In agreement with a possible oxygen carrier and nitrite reductase role, we detect mRNA transcript for Cygb1 but not Cygb2 or Ngb in zebrafish blood. Unlike human Cygb, neither of the zebrafish globins binds oleic acid with high affinity. This finding suggests that lipid binding may be a trait acquired later during evolution and not an ancestral property of cytoglobins. Altogether, our results uncover unexpected properties of zebrafish globins and reveal the pivotal role of cytoglobins in the transition of heme globins from six-coordinate to five-coordinate oxygen carriers and nitrite reductases.

Convergent Evolution of Hemoglobin Function in High-Altitude Andean Waterfowl Involves Limited Parallelism at the Molecular Sequence Level

A fundamental question in evolutionary genetics concerns the extent to which adaptive phenotypic convergence is attributable to convergent or parallel changes at the molecular sequence level. Here we report a comparative analysis of hemoglobin (Hb) function in eight phylogenetically replicated pairs of high- and low-altitude waterfowl taxa to test for convergence in the oxygenation properties of Hb, and to assess the extent to which convergence in biochemical phenotype is attributable to repeated amino acid replacements. Functional experiments on native Hb variants and protein engineering experiments based on site-directed mutagenesis revealed the phenotypic effects of specific amino acid replacements that were responsible for convergent increases in Hb-O₂ affinity in multiple high-altitude taxa. In six of the eight taxon pairs, high-altitude taxa evolved derived increases in Hb-O₂ affinity that were caused by a combination of unique replacements, parallel replacements (involving identical-by-state variants with independent mutational origins in different lineages), and collateral replacements (involving shared, identical-by-descent variants derived via introgressive hybridization). In genome scans of nucleotide differentiation involving high- and low-altitude populations of three separate species, function-altering amino acid polymorphisms in the globin genes emerged as highly significant outliers, providing independent evidence for adaptive divergence in Hb function. The experimental results demonstrate that convergent changes in protein function can occur through multiple historical paths, and can involve multiple possible mutations. Most cases of convergence in Hb function did not involve parallel substitutions and most parallel substitutions did not affect Hb-O₂ affinity, indicating that the repeatability of phenotypic evolution does not require parallelism at the molecular level.

The convergent evolution of similar traits in different species could be due to repeated changes at the genetic level or different changes that produce the same phenotypic effect. To investigate the extent to which convergence in phenotype is caused by repeated mutations, we investigated the molecular basis of convergent changes in the oxygenation properties of hemoglobin (Hb) in eight pairs of high- and low-altitude waterfowl taxa from the Andes. The results revealed that convergent increases in Hb-O₂ affinity in highland taxa involved a combination of unique and repeated amino acid replacements. However, convergent changes in Hb function generally did not involve parallel substitutions, indicating that repeatability in the evolution of protein function does not require repeatability at the sequence level.

Oxygenation properties and isoform diversity of snake hemoglobins

Available data suggest that snake hemoglobins (Hbs) are characterized by a combination of unusual structural and functional properties relative to the Hbs of other amniote vertebrates, including oxygenation-linked tetramer-dimer dissociation. However, standardized comparative data are lacking for snake Hbs, and the Hb isoform composition of snake red blood cells has not been systematically characterized. Here we present the results of an integrated analysis of snake Hbs and the underlying α- and β-type globin genes to characterize 1) Hb isoform composition of definitive erythrocytes, and 2) the oxygenation properties of isolated isoforms as well as composite hemolysates. We used species from three families as subjects for experimental studies of Hb function: South American rattlesnake, Crotalus durissus (Viperidae); Indian python, Python molurus (Pythonidae); and yellow-bellied sea snake, Pelamis platura (Elapidae). We analyzed allosteric properties of snake Hbs in terms of the Monod-Wyman-Changeux model and Adair four-step thermodynamic model. Hbs from each of the three species exhibited high intrinsic O2 affinities, low cooperativities, small Bohr factors in the absence of phosphates, and high sensitivities to ATP. Oxygenation properties of the snake Hbs could be explained entirely by allosteric transitions in the quaternary structure of intact tetramers, suggesting that ligation-dependent dissociation of Hb tetramers into αβ-dimers is not a universal feature of snake Hbs. Surprisingly, the major Hb isoform of the South American rattlesnake is homologous to the minor HbD of other amniotes and, contrary to the pattern of Hb isoform differentiation in birds and turtles, exhibits a lower O2 affinity than the HbA isoform.

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.

The Full Globin Repertoire of Turtles Provides Insights into Vertebrate Globin Evolution and Functions

Globins are small heme proteins that play an important role in oxygen supply, but may also have other functions. Globins offer a unique opportunity to study the functional evolution of genes and proteins. We have characterized the globin repertoire of two different turtle species: the Chinese softshell turtle (Pelodiscus sinensis) and the western painted turtle (Chrysemys picta bellii). In the genomes of both species, we have identified eight distinct globin types: hemoglobin (Hb), myoglobin, neuroglobin, cytoglobin, globin E, globin X, globin Y, and androglobin. Therefore, along with the coelacanth, turtles are so far the only known vertebrates with a full globin repertoire. This fact allows for the first time a comparative analysis of the expression of all eight globins in a single species. Phylogenetic analysis showed an early divergence of neuroglobin and globin X before the radiation of vertebrates. Among the other globins, cytoglobin diverged first, and there is a close relationship between myoglobin and globin E; the position of globin Y is not resolved. The globin E gene was selectively lost in the green anole, and the genes coding for globin X and globin Y were deleted in chicken. Quantitative real-time reverse transcription polymerase chain reaction experiments revealed that myoglobin, neuroglobin, and globin E are highly expressed with tissue-specific patterns, which are in line with their roles in the oxidative metabolism of the striated muscles, the brain, and the retina, respectively. Histochemical analyses showed high levels of globin E in the pigment epithelium of the eye. Globin E probably has a myoglobin-like role in transporting O₂ across the pigment epithelium to supply in the metabolically highly active retina.

Gene turnover in the avian globin gene families and evolutionary changes in hemoglobin isoform expression

The apparent stasis in the evolution of avian chromosomes suggests that birds may have experienced relatively low rates of gene gain and loss in multigene families. To investigate this possibility and to explore the phenotypic consequences of variation in gene copy number, we examined evolutionary changes in the families of genes that encode the α- and β-type subunits of hemoglobin (Hb), the tetrameric α2β2 protein responsible for blood-O2 transport. A comparative genomic analysis of 52 bird species revealed that the size and membership composition of the α- and β-globin gene families have remained remarkably constant during approximately 100 My of avian evolution. Most interspecific variation in gene content is attributable to multiple independent inactivations of the α(D)-globin gene, which encodes the α-chain subunit of a functionally distinct Hb isoform (HbD) that is expressed in both embryonic and definitive erythrocytes. Due to consistent differences in O2-binding properties between HbD and the major adult-expressed Hb isoform, HbA (which incorporates products of the α(A)-globin gene), recurrent losses of α(D)-globin contribute to among-species variation in blood-O2 affinity. Analysis of HbA/HbD expression levels in the red blood cells of 122 bird species revealed high variability among lineages and strong phylogenetic signal. In comparison with the homologous gene clusters in mammals, the low retention rate for lineage-specific gene duplicates in the avian globin gene clusters suggests that the developmental regulation of Hb synthesis in birds may be more highly conserved, with orthologous genes having similar stage-specific expression profiles and similar functional properties in disparate taxa.

Lineage-specific loss of FGF17 within the avian orders Galliformes and Passeriformes

The genomic and developmental complexity of vertebrates is commonly attributed to two rounds of whole genome duplications which occurred at the base of the vertebrate radiation. These duplications led to the rise of several, multi-gene families of developmental proteins like the fibroblast growth factors (FGFs); a signaling protein family which functions at various stages of embryonic development. One of the major FGF assemblages arising from these duplications is the FGF8 subfamily, which includes FGF8, FGF17, and FGF18 in tetrapods. While FGF8 and FGF18 are found in all tetrapods and are critical for embryonic survival, genomic analyses suggest putative loss of FGF17 in various lineages ranging from frogs and fish, to the chicken. This study utilizes 27 avian genomes in conjunction with molecular analyses of chicken embryos to confirm the loss of FGF17 in chicken as a true, biological occurrence. FGF17 is also missing in the turkey, black grouse, Japanese quail and northern bobwhite genomes. These species, along with chicken, form a monophyletic clade in the order Galliformes. Four additional species, members of the clade Passeroidea, within the order Passeriformes, are also missing FGF17. Additionally, analysis of intact FGF17 in other avian lineages reveals that it is still under strong purifying selection, despite being seemingly dispensable. Thus, FGF17 likely represents a molecular spandrel arising from a genome duplication event and due to its high connectivity with FGF8/FGF18, and potential for interference with their function, is retained under strong purifying selection, despite itself not having a strong selective advantage.

Ancient Duplications and Expression Divergence in the Globin Gene Superfamily of Vertebrates: Insights from the Elephant Shark Genome and Transcriptome

Comparative analyses of vertebrate genomes continue to uncover a surprising diversity of genes in the globin gene superfamily, some of which have very restricted phyletic distributions despite their antiquity. Genomic analysis of the globin gene repertoire of cartilaginous fish (Chondrichthyes) should be especially informative about the duplicative origins and ancestral functions of vertebrate globins, as divergence between Chondrichthyes and bony vertebrates represents the most basal split within the jawed vertebrates. Here, we report a comparative genomic analysis of the vertebrate globin gene family that includes the complete globin gene repertoire of the elephant shark (Callorhinchus milii). Using genomic sequence data from representatives of all major vertebrate classes, integrated analyses of conserved synteny and phylogenetic relationships revealed that the last common ancestor of vertebrates possessed a repertoire of at least seven globin genes: single copies of androglobin and neuroglobin, four paralogous copies of globin X, and the single-copy progenitor of the entire set of vertebrate-specific globins. Combined with expression data, the genomic inventory of elephant shark globins yielded four especially surprising findings: 1) there is no trace of the neuroglobin gene (a highly conserved gene that is present in all other jawed vertebrates that have been examined to date), 2) myoglobin is highly expressed in heart, but not in skeletal muscle (reflecting a possible ancestral condition in vertebrates with single-circuit circulatory systems), 3) elephant shark possesses two highly divergent globin X paralogs, one of which is preferentially expressed in gonads, and 4) elephant shark possesses two structurally distinct α-globin paralogs, one of which is preferentially expressed in the brain. Expression profiles of elephant shark globin genes reveal distinct specializations of function relative to orthologs in bony vertebrates and suggest hypotheses about ancestral functions of vertebrate globins.

Molecular characterization, phylogenetic analysis and expression profiling of myoglobin and cytoglobin genes in response to heat stress in channel catfish Ictalurus punctatus

To understand the function of myoglobin (Mb) and cytoglobin (Cygb) in channel catfish Ictalurus punctatus in response to heat stress, mb and cygb genes were identified and characterized in this study. These genes were widely expressed in all the tested tissues, but strong tissue preferences were observed, with the mb gene being expressed most highly in the heart, cygb1 most highly expressed in the intestine and cygb2 most highly expressed in the brain. After heat-stress challenge, mb and cygb genes were up-regulated in almost all tested tissues. In general, such up-regulation was more dramatic in the tolerant group than in the intolerant group, suggesting that higher expression of mb and cygb genes contributed to greater tolerance of I. punctatus to heat stress.

Basal Gnathostomes provide unique insights into the evolution of vitamin B12 binders

The uptake and transport of vitamin B12 (cobalamin; Cbl) in mammals involves a refined system with three evolutionarily related transporters: transcobalamin 1 (Tcn1), transcobalamin 2 (Tcn2), and the gastric intrinsic factor (Gif). Teleosts have a single documented binder with intermediate features to the human counterparts. Consequently, it has been proposed that the expansion of Cbl binders occurred after the separation of Actinopterygians. Here, we demonstrate that the diversification of this gene family took place earlier in gnathostome ancestry. Our data indicates the presence of single copy orthologs of the Sarcopterygii/Tetrapoda duplicates Tcn1 and Gif, and Tcn2, in Chondrichthyes. In addition, a highly divergent Cbl binder was found in the Elasmobranchii. We unveil a complex scenario forged by genome, tandem duplications and lineage-specific gene loss. Our findings suggest that from an ancestral transporter, exhibiting large spectrum and high affinity binding, highly specific Cbl transporters emerged through gene duplication and mutations at the binding pocket.

Integrating evolutionary and functional tests of adaptive hypotheses: a case study of altitudinal differentiation in hemoglobin function in an Andean Sparrow, Zonotrichia capensis

In air-breathing vertebrates, the physiologically optimal blood-O2 affinity is jointly determined by the prevailing partial pressure of atmospheric O2, the efficacy of pulmonary O2 transfer, and internal metabolic demands. Consequently, genetic variation in the oxygenation properties of hemoglobin (Hb) may be subject to spatially varying selection in species with broad elevational distributions. Here we report the results of a combined functional and evolutionary analysis of Hb polymorphism in the rufous-collared sparrow (Zonotrichia capensis), a species that is continuously distributed across a steep elevational gradient on the Pacific slope of the Peruvian Andes. We integrated a population genomic analysis that included all postnatally expressed Hb genes with functional studies of naturally occurring Hb variants, as well as recombinant Hb (rHb) mutants that were engineered through site-directed mutagenesis. We identified three clinally varying amino acid polymorphisms: Two in the α(A)-globin gene, which encodes the α-chain subunits of the major HbA isoform, and one in the α(D)-globin gene, which encodes the α-chain subunits of the minor HbD isoform. We then constructed and experimentally tested single- and double-mutant rHbs representing each of the alternative α(A)-globin genotypes that predominate at different elevations. Although the locus-specific patterns of altitudinal differentiation suggested a history of spatially varying selection acting on Hb polymorphism, the experimental tests demonstrated that the observed amino acid mutations have no discernible effect on respiratory properties of the HbA or HbD isoforms. These results highlight the importance of experimentally validating the hypothesized effects of genetic changes in protein function to avoid the pitfalls of adaptive storytelling.

The globin gene repertoire of lampreys: convergent evolution of hemoglobin and myoglobin in jawed and jawless vertebrates

Agnathans (jawless vertebrates) occupy a key phylogenetic position for illuminating the evolution of vertebrate anatomy and physiology. Evaluation of the agnathan globin gene repertoire can thus aid efforts to reconstruct the origin and evolution of the globin genes of vertebrates, a superfamily that includes the well-known model proteins hemoglobin and myoglobin. Here, we report a comprehensive analysis of the genome of the sea lamprey (Petromyzon marinus) which revealed 23 intact globin genes and two hemoglobin pseudogenes. Analyses of the genome of the Arctic lamprey (Lethenteron camtschaticum) identified 18 full length and five partial globin gene sequences. The majority of the globin genes in both lamprey species correspond to the known agnathan hemoglobins. Both genomes harbor two copies of globin X, an ancient globin gene that has a broad phylogenetic distribution in the animal kingdom. Surprisingly, we found no evidence for an ortholog of neuroglobin in the lamprey genomes. Expression and phylogenetic analyses identified an ortholog of cytoglobin in the lampreys; in fact, our results indicate that cytoglobin is the only orthologous vertebrate-specific globin that has been retained in both gnathostomes and agnathans. Notably, we also found two globins that are highly expressed in the heart of P. marinus, thus representing functional myoglobins. Both genes have orthologs in L. camtschaticum. Phylogenetic analyses indicate that these heart-expressed globins are not orthologous to the myoglobins of jawed vertebrates (Gnathostomata), but originated independently within the agnathans. The agnathan myoglobin and hemoglobin proteins form a monophyletic group to the exclusion of functionally analogous myoglobins and hemoglobins of gnathostomes, indicating that specialized respiratory proteins for O2 transport in the blood and O2 storage in the striated muscles evolved independently in both lineages. This dual convergence of O2-transport and O2-storage proteins in agnathans and gnathostomes involved the convergent co-option of different precursor proteins in the ancestral globin repertoire of vertebrates.

Function and evolution of vertebrate globins

Globins are haem-proteins that bind O2 and thus play an important role in the animal's respiration and oxidative energy production. However, globins may also have other functions such as the decomposition or production of NO, the detoxification of reactive oxygen species or intracellular signalling. In addition to the well-investigated haemoglobins and myoglobins, genome sequence analyses have led to the identification of six further globin types in vertebrates: androglobin, cytoglobin, globin E, globin X, globin Y and neuroglobin. Here, we review the present state of knowledge on the functions, the taxonomic distribution and evolution of vertebrate globins, drawing conclusions about the functional changes underlying present-day globin diversity.

Cardiac myoglobin deficit has evolved repeatedly in teleost fishes

Myoglobin (Mb) is the classic vertebrate oxygen-binding protein present in aerobic striated muscles. It functions principally in oxygen delivery and provides muscle with its characteristic red colour. Members of the Antarctic icefish family (Channichthyidae) are widely thought to be extraordinary for lacking cardiac Mb expression, a fact that has been attributed to their low metabolic rate and unusual evolutionary history. Here, we report that cardiac Mb deficit, associated with pale heart colour, has evolved repeatedly during teleost evolution. This trait affects both gill- and air-breathing species from temperate to tropical habitats across a full range of salinities. Cardiac Mb deficit results from total pseudogenization in three-spined stickleback and is associated with a massive reduction in mRNA level in two species that evidently retain functional Mb. The results suggest that near or complete absence of Mb-assisted oxygen delivery to heart muscle is a common facet of teleost biodiversity, even affecting lineages with notable oxygen demands. We suggest that Mb deficit may affect how different teleost species deal with increased tissue oxygen demands arising under climate change.

Repeated elevational transitions in hemoglobin function during the evolution of Andean hummingbirds

Animals that sustain high levels of aerobic activity under hypoxic conditions (e.g., birds that fly at high altitude) face the physiological challenge of jointly optimizing blood-O2 affinity for O2 loading in the pulmonary circulation and O2 unloading in the systemic circulation. At high altitude, this challenge is especially acute for small endotherms like hummingbirds that have exceedingly high mass-specific metabolic rates. Here we report an experimental analysis of hemoglobin (Hb) function in South American hummingbirds that revealed a positive correlation between Hb-O2 affinity and native elevation. Protein engineering experiments and ancestral-state reconstructions revealed that this correlation is attributable to derived increases in Hb-O2 affinity in highland lineages, as well as derived reductions in Hb-O2 affinity in lowland lineages. Site-directed mutagenesis experiments demonstrated that repeated evolutionary transitions in biochemical phenotype are mainly attributable to repeated amino acid replacements at two epistatically interacting sites that alter the allosteric regulation of Hb-O2 affinity. These results demonstrate that repeated changes in biochemical phenotype involve parallelism at the molecular level, and that mutations with indirect, second-order effects on Hb allostery play key roles in biochemical adaptation.

Expression patterns and adaptive functional diversity of vertebrate myoglobins

Recent years have witnessed a new round of research on one of the most studied proteins - myoglobin (Mb), the oxygen (O2) carrier of skeletal and heart muscle. Two major discoveries have stimulated research in this field: 1) that Mb has additional protecting functions, such as the regulation of in vivo levels of the signaling molecule nitric oxide (NO) by scavenging and generating NO during normoxia and hypoxia, respectively; and 2) that Mb in vertebrates (particularly fish) is expressed as tissue-specific isoforms in other tissues than heart and skeletal muscle, such as vessel endothelium, liver and brain, as found in cyprinid fish. Furthermore, Mb has also been found to protect against oxidative stress after hypoxia and reoxygenation and to undergo allosteric, O2-linked S-nitrosation, as in rainbow trout. Overall, the emerging evidence, particularly from fish species, indicates that Mb fulfills a broader array of physiological functions in a wider range of different tissues than hitherto appreciated. This new knowledge helps to better understand how variations in Mb structure and function may correlate with differences in animals' lifestyles and hypoxia-tolerance. This review integrates old and new results on Mb expression patterns and functional properties amongst vertebrates and discusses how these may relate to adaptive variations in different species. This article is part of a special issue entitled: Oxygen Binding and Sensing Proteins.

Conservation of globin genes in the "living fossil" Latimeria chalumnae and reconstruction of the evolution of the vertebrate globin family

The (hemo-)globins are among the best-investigated proteins in biomedical sciences. These small heme-proteins play an important role in oxygen supply, but may also have other functions. In addition to well known hemoglobin and myoglobin, six other vertebrate globin types have been identified in recent years: neuroglobin, cytoglobin, globin E, globin X, globin Y, and androglobin. Analyses of the genome of the "living fossil" Latimeria chalumnae show that the coelacanth is the only known vertebrate that includes all eight globin types. Thus, Latimeria can also be considered as a "globin fossil". Analyses of gene synteny and phylogenetic reconstructions allow us to trace the evolution and the functional changes of the vertebrate globin family. Neuroglobin and globin X diverged from the other globin types before the separation of Protostomia and Deuterostomia. The cytoglobins, which are unlikely to be involved in O2 supply, form the earliest globin branch within the jawed vertebrates (Gnathostomata), but do not group with the agnathan hemoglobins, as it has been proposed before. There is strong evidence from phylogenetic reconstructions and gene synteny that the eye-specific globin E and muscle-specific myoglobin constitute a common clade, suggesting a similar role in intracellular O2 supply. Latimeria possesses two α- and two β-hemoglobin chains, of which one α-chain emerged prior to the divergence of Actinopterygii and Sarcopterygii, but has been retained only in the coelacanth. Notably, the embryonic hemoglobin α-chains of Gnathostomata derive from a common ancestor, while the embryonic β-chains - with the exception of a more complex pattern in the coelacanth and amphibians - display a clade-specific evolution. Globin Y is associated with the hemoglobin gene cluster, but its phylogenetic position is not resolved. Our data show an early divergence of distinct globin types in the vertebrate evolution before the emergence of tetrapods. The subsequent loss of globins in certain taxa may be associated with changes in the oxygen-dependent metabolism. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Gene duplication and the evolution of hemoglobin isoform differentiation in birds

The majority of bird species co-express two functionally distinct hemoglobin (Hb) isoforms in definitive erythrocytes as follows: HbA (the major adult Hb isoform, with α-chain subunits encoded by the α(A)-globin gene) and HbD (the minor adult Hb isoform, with α-chain subunits encoded by the α(D)-globin gene). The α(D)-globin gene originated via tandem duplication of an embryonic α-like globin gene in the stem lineage of tetrapod vertebrates, which suggests the possibility that functional differentiation between the HbA and HbD isoforms may be attributable to a retained ancestral character state in HbD that harkens back to a primordial, embryonic function. To investigate this possibility, we conducted a combined analysis of protein biochemistry and sequence evolution to characterize the structural and functional basis of Hb isoform differentiation in birds. Functional experiments involving purified HbA and HbD isoforms from 11 different bird species revealed that HbD is characterized by a consistently higher O(2) affinity in the presence of allosteric effectors such as organic phosphates and Cl(-) ions. In the case of both HbA and HbD, analyses of oxygenation properties under the two-state Monod-Wyman-Changeux allosteric model revealed that the pH dependence of Hb-O(2) affinity stems primarily from changes in the O(2) association constant of deoxy (T-state)-Hb. Ancestral sequence reconstructions revealed that the amino acid substitutions that distinguish the adult-expressed Hb isoforms are not attributable to the retention of an ancestral (pre-duplication) character state in the α(D)-globin gene that is shared with the embryonic α-like globin gene.

Comparative genomics of neuroglobin reveals its early origins

BACKGROUND

Neuroglobin (Ngb) is a hexacoordinated globin expressed mainly in the central and peripheral nervous system of vertebrates. Although several hypotheses have been put forward regarding the role of neuroglobin, its definite function remains uncertain. Ngb appears to have a neuro-protective role enhancing cell viability under hypoxia and other types of oxidative stress. Ngb is phylogenetically ancient and has a substitution rate nearly four times lower than that of other vertebrate globins, e.g. hemoglobin. Despite its high sequence conservation among vertebrates Ngb seems to be elusive in invertebrates.

PRINCIPAL FINDINGS

We determined candidate orthologs in invertebrates and identified a globin of the placozoan Trichoplax adhaerens that is most likely orthologous to vertebrate Ngb and confirmed the orthologous relationship of the polymeric globin of the sea urchin Strongylocentrotus purpuratus to Ngb. The putative orthologous globin genes are located next to genes orthologous to vertebrate POMT2 similarly to localization of vertebrate Ngb. The shared syntenic position of the globins from Trichoplax, the sea urchin and of vertebrate Ngb strongly suggests that they are orthologous. A search for conserved transcription factor binding sites (TFBSs) in the promoter regions of the Ngb genes of different vertebrates via phylogenetic footprinting revealed several TFBSs, which may contribute to the specific expression of Ngb, whereas a comparative analysis with myoglobin revealed several common TFBSs, suggestive of regulatory mechanisms common to globin genes.

SIGNIFICANCE

Identification of the placozoan and echinoderm genes orthologous to vertebrate neuroglobin strongly supports the hypothesis of the early evolutionary origin of this globin, as it shows that neuroglobin was already present in the placozoan-bilaterian last common ancestor. Computational determination of the transcription factor binding sites repertoire provides on the one hand a set of transcriptional factors that are responsible for the specific expression of the Ngb genes and on the other hand a set of factors potentially controlling expression of a couple of different globin genes.

Oxygen levels and the regulation of cell adhesion in the nervous system: a control point for morphogenesis in development, disease and evolution?

In this article, I discuss the hallmarks of hypoxia in vitro and in vivo and review work showing that many types of stem cell proliferate more robustly in lowered oxygen. I then discuss recent studies showing that alterations in the levels and the types of cell and substrate adhesion molecules are a notable response to reduced O(2) levels in both cultured primary neural stem cells and brain tissues in response to hypoxia in vivo. The ability of O(2) levels to regulate adhesion molecule expression is linked to the Wnt signaling pathway, which can control and be controlled by adhesion events. The ability of O(2) levels to influence cell adhesion also has far-reaching implications for development, ischemic trauma and neural regeneration, as well as for cancer and other diseases. Finally I discuss the possibility that the fluctuations in O(2) levels known to have occurred over evolutionary time could, by influencing adhesion systems, have contributed to early symbiotic events in unicellular organisms and to the emergence of multicellularity. It is not my intention to be exhaustive in these domains, which are far from my own field of study. Rather this article is meant to provoke and stimulate thinking about molecular evolution involving O(2) sensing and signaling during eras of geologic and atmospheric change that might inform modern studies on development and disease.