Lack of conventional oxygen-linked proton and anion binding sites does not impair allosteric regulation of oxygen binding in dwarf caiman hemoglobin

Evolution and molecular basis of a novel allosteric property of crocodilian hemoglobin

The extraordinary breath-hold diving capacity of crocodilians has been ascribed to a unique mode of allosterically regulating hemoglobin (Hb)-oxygenation in circulating red blood cells. We investigated the origin and mechanistic basis of this novel biochemical phenomenon by performing directed mutagenesis experiments on resurrected ancestral Hbs. Comparisons of Hb function between the common ancestor of archosaurs (the group that includes crocodilians and birds) and the last common ancestor of modern crocodilians revealed that regulation of Hb-O2 affinity via allosteric binding of bicarbonate ions represents a croc-specific innovation that evolved in combination with the loss of allosteric regulation by ATP binding. Mutagenesis experiments revealed that evolution of the novel allosteric function in crocodilians and the concomitant loss of ancestral function were not mechanistically coupled and were caused by different sets of substitutions. The gain of bicarbonate sensitivity in crocodilian Hb involved the direct effect of few amino acid substitutions at key sites in combination with indirect effects of numerous other substitutions at structurally disparate sites. Such indirect interaction effects suggest that evolution of the novel protein function was conditional on neutral mutations that produced no adaptive benefit when they first arose but that contributed to a permissive background for subsequent function-altering mutations at other sites. Due to the context dependence of causative substitutions, the unique allosteric properties of crocodilian Hb cannot be easily transplanted into divergent homologs of other species.

Changes in hemoglobin function and isoform expression during embryonic development in the American alligator, Alligator mississippiensis

In the developing embryos of egg-laying vertebrates, O₂ flux takes place across a fixed surface area of the eggshell and the chorioallantoic membrane. In the case of crocodilians, the developing embryo may experience a decrease in O₂ flux when the nest becomes hypoxic, which may cause compensatory adjustments in blood O₂ transport. However, whether the switch from embryonic to adult hemoglobin isoforms (isoHbs) plays some role in these adjustments is unknown. Here, we provide a detailed characterization of the developmental switch of isoHb synthesis in the American alligator, Alligator mississippiensis. We examined the in vitro functional properties and subunit composition of purified alligator isoHbs expressed during embryonic developmental stages in normoxia and hypoxia (10% O₂). We found distinct patterns of isoHb expression in alligator embryos at different stages of development, but these patterns were not affected by hypoxia. Specifically, alligator embryos expressed two main isoHbs: HbI, prevalent at early developmental stages, with a high O₂ affinity and high ATP sensitivity, and HbII, prevalent at later stages and identical to the adult protein, with a low O₂ affinity and high CO₂ sensitivity. These results indicate that whole blood O₂ affinity is mainly regulated by ATP in the early embryo and by CO₂ and bicarbonate from the late embryo until adult life, but the developmental regulation of isoHb expression is not affected by hypoxia exposure.

New insights into the allosteric effects of CO2 and bicarbonate on crocodilian hemoglobin

Crocodilians are unique among vertebrates in that their hemoglobin (Hb) O2 binding is allosterically regulated by bicarbonate, which forms in red blood cells upon hydration of CO2. Although known for decades, this remarkable mode of allosteric control has not yet been experimentally verified with direct evidence of bicarbonate binding to crocodilian Hb, probably because of confounding CO2-mediated effects. Here, we provide the first quantitative analysis of the separate allosteric effects of CO2 and bicarbonate on purified Hb of the spectacled caiman (Caiman crocodilus). Using thin-layer gas diffusion chamber and Tucker chamber techniques, we demonstrate that both CO2 and bicarbonate bind to Hb with high affinity and strongly decrease O2 saturation of Hb. We propose that both effectors bind to an unidentified positively charged site containing a reactive amino group in the low-O2 affinity T conformation of Hb. These results provide the first experimental evidence that bicarbonate binds directly to crocodilian Hb and promotes O2 delivery independently of CO2. Using the gas diffusion chamber, we observed similar effects in Hbs of a phylogenetically diverse set of other caiman, alligator and crocodile species, suggesting that the unique mode of allosteric regulation by CO2 and bicarbonate evolved >80-100 million years ago in the common ancestor of crocodilians. Our results show a tight and unusual linkage between O2 and CO2 transport in the blood of crocodilians, where the build-up of erytrocytic CO2 and bicarbonate ions during breath-hold diving or digestion facilitates O2 delivery, while Hb desaturation facilitates CO2 transport as protein-bound CO2 and bicarbonate.

Carbon dioxide and bicarbonate accumulation in caiman erythrocytes during diving

The ability of crocodilian haemoglobins to bind HCO3 - has been appreciated for more than half a century, but the functional implication of this is exceptional mechanism has not previously been assessed in vivo Therefore, the goal of the present study was to address the hypothesis that CO2 primarily binds to Hb, rather than being accumulated in plasma as in other vertebrates, during diving in caimans. Here, we demonstrate that CO2 primarily accumulates within the erythrocyte during diving and that most of the accumulated CO2 is bound to haemoglobin. Furthermore, we show that this HCO3 --binding is tightly associated with the progressive blood deoxygenation during diving, therefore, crocodilians differ from the classic vertebrate pattern, where HCO3 - accumulates in the plasma upon excretion from the erythrocytes by the Cl--HCO3 --exchanger.

Effect of NH2-terminal acetylation on the oxygenation properties of vertebrate haemoglobin

In vertebrate haemoglobin (Hb), the NH2-terminal residues of the α- and β-chain subunits are thought to play an important role in the allosteric binding of protons (Bohr effect), CO2 (as carbamino derivatives), chloride ions, and organic phosphates. Accordingly, acetylation of the α- and/or β-chain NH2-termini may have significant effects on the oxygenation properties of Hb. Here we investigate the effect of NH2-terminal acetylation by using a newly developed expression plasmid system that enables us to compare recombinantly expressed Hbs that are structurally identical except for the presence or absence of NH2-terminal acetyl groups. Experiments with native and recombinant Hbs of representative vertebrates reveal that NH2-terminal acetylation does not impair the Bohr effect, nor does it significantly diminish responsiveness to allosteric cofactors, such as chloride ions or organic phosphates. These results suggest that observed variation in the oxygenation properties of vertebrate Hbs is principally explained by amino acid divergence in the constituent globin chains rather than post-translational modifications of the globin chain NH2-termini.

Recent genome duplications facilitate the phenotypic diversity of Hb repertoire in the Cyprinidae

Whole-genome duplications (WGDs) are an important contributor to phenotypic innovations in evolutionary history. The diversity of blood oxygen transport traits is the perfect reflection of physiological versatility for evolutionary success among vertebrates. In this study, the evolutionary changes of hemoglobin (Hb) repertoire driven by the recent genome duplications were detected in representative Cyprinidae fish, including eight diploid and four tetraploid species. Comparative genomic analysis revealed a substantial variation in both membership composition and intragenomic organization of Hb genes in these species. Phylogenetic reconstruction analyses were conducted to characterize the evolutionary history of these genes. Data were integrated with the expression profiles of the genes during ontogeny. Our results indicated that genome duplications facilitated the phenotypic diversity of the Hb gene family; each was associated with species-specific changes in gene content via gene loss and fusion after genome duplications. This led to repeated evolutionary transitions in the ontogenic regulation of Hb gene expression. Our results revealed that genome duplications helped to generate phenotypic changes in Cyprinidae Hb systems.

Synthesis of Recombinant Human Hemoglobin With NH2 -Terminal Acetylation in Escherichia coli

The development of new technologies for the efficient expression of recombinant hemoglobin (rHb) is of interest for experimental studies of protein biochemistry and the development of cell-free blood substitutes in transfusion medicine. Expression of rHb in Escherichia coli host cells has numerous advantages, but one disadvantage of using prokaryotic systems to express eukaryotic proteins is that they are incapable of performing post-translational modifications such as NH₂ -terminal acetylation. One possible solution is to coexpress additional enzymes that can perform the necessary modifications in the host cells. Here, we report a new method for synthesizing human rHb with proper NH₂ -terminal acetylation. Mass spectrometry experiments involving native and recombinant human Hb confirmed the efficacy of the new technique in producing correctly acetylated globin chains. Finally, functional experiments provided insights into the effects of NH₂ -terminal acetylation on O₂ binding properties. © 2020 Wiley Periodicals LLC. Basic Protocol 1: Gene synthesis and cloning the cassette to the expression plasmid Basic Protocol 2: Selection of E. coli expression strains for coexpression Basic Protocol 3: Large-scale recombinant hemoglobin expression and purification Support Protocol 1: Measuring O₂ equilibration curves Support Protocol 2: Mass spectrometry to confirm NH₂ -terminal acetylation.

Structure and function of crocodilian hemoglobins and allosteric regulation by chloride, ATP, and CO2

Hemoglobins (Hbs) of crocodilians are reportedly characterized by unique mechanisms of allosteric regulatory control, but there are conflicting reports regarding the importance of different effectors, such as chloride ions, organic phosphates, and CO₂. Progress in understanding the unusual properties of crocodilian Hbs has also been hindered by a dearth of structural information. Here, we present the first comparative analysis of blood properties and Hb structure and function in a phylogenetically diverse set of crocodilian species. We examine mechanisms of allosteric regulation in the Hbs of 13 crocodilian species belonging to the families Crocodylidae and Alligatoridae. We also report new amino acid sequences for the α- and β-globins of these taxa, which, in combination with structural analyses, provide insights into molecular mechanisms of allosteric regulation. All crocodilian Hbs exhibited a remarkably strong sensitivity to CO₂, which would permit effective O₂ unloading to tissues in response to an increase in metabolism during intense activity and diving. Although the Hbs of all crocodilians exhibit similar intrinsic O₂-affinities, there is considerable variation in sensitivity to Cl^- ions and ATP, which appears to be at least partly attributable to variation in the extent of NH₂-terminal acetylation. Whereas chloride appears to be a potent allosteric effector of all crocodile Hbs, ATP has a strong, chloride-independent effect on Hb-O₂ affinity only in caimans. Modeling suggests that allosteric ATP binding has a somewhat different structural basis in crocodilian and mammalian Hbs.

Oxygenation properties of hemoglobin and the evolutionary origins of isoform multiplicity in an amphibious air-breathing fish, the blue-spotted mudskipper (Boleophthalmus pectinirostris)

Among the numerous lineages of teleost fish that have independently transitioned from obligate water breathing to facultative air breathing, evolved properties of hemoglobin (Hb)-O₂ transport may have been shaped by the prevalence and severity of aquatic hypoxia (which influences the extent to which fish are compelled to switch to aerial respiration) as well as the anatomical design of air-breathing structures and the cardiovascular system. Here, we examined the structure and function of Hbs in an amphibious, facultative air-breathing fish, the blue-spotted mudskipper (Boleophthalmus pectinirostris). We also characterized the genomic organization of the globin gene clusters of the species and we integrated phylogenetic and comparative genomic analyses to unravel the duplicative history of the genes that encode the subunits of structurally distinct mudskipper Hb isoforms (isoHbs). The B. pectinirostris isoHbs exhibit high intrinsic O₂ affinities, similar to those of hypoxia-tolerant, water-breathing teleosts, and remarkably large Bohr effects. Genomic analysis of conserved synteny revealed that the genes that encode the α-type subunits of the two main adult isoHbs are members of paralogous gene clusters that represent products of the teleost-specific whole-genome duplication. Experiments revealed no appreciable difference in the oxygenation properties of co-expressed isoHbs in spite of extensive amino acid divergence between the alternative α-chain subunit isoforms. It therefore appears that the ability to switch between aquatic and aerial respiration does not necessarily require a division of labor between functionally distinct isoHbs with specialized oxygenation properties.

Oxygenation properties and underlying molecular mechanisms of hemoglobins in plateau zokor (Eospalax baileyi)

The plateau zokor (Eospalax baileyi) is a species of subterranean rodent endemic to the Tibetan Plateau. It is well adapted to the cold and hypoxic and hypercapnic burrow. To study the oxygenation properties of plateau zokor hemoglobins (Hbs), we measured intrinsic Hb-O₂ affinities and their sensitivities to pH (Bohr effect); CO₂; Cl^-, 2,3-diphosphoglycerate (DPG); and temperature using purified Hbs from zokor and mouse. The optimal deoxyHb model of plateau zokor was constructed and used to study its structural characteristics by molecular dynamics simulations. O₂ binding results revealed that plateau zokor Hbs exhibit remarkably high intrinsic Hb-O₂ affinity, low CO₂ effects compared with human and the relatively low anion allosteric effector sensitivities (DPG and Cl^-) at normal temperature, which would safeguard the pulmonary Hb-O₂ loading under hypoxic and hypercapnic conditions. Furthermore, the high anion allosteric effector sensitivities at low temperature and low temperature sensitivities of plateau zokor Hbs would facilitate the releasing of O₂ in cold extremities and metabolic tissues. However, the high Hb-O₂ affinity of plateau zokor is not compensated by high pH sensitivity as the Bohr factors of plateau zokor Hbs were as low as those of mouse. The results of molecular dynamics simulations revealed the reduced hydrogen bonding between the α1β1- and α2β2-dimer interface of deoxyHb in zokor compared with mouse. It may be the primary mechanism of the high intrinsic Hb-O₂ affinities in zokor. Specifically, substitution of the 131Ser→Asn in the α2-chain weakened the connection between α1- and β2-subunit.

The role of mutation bias in adaptive molecular evolution: insights from convergent changes in protein function

An underexplored question in evolutionary genetics concerns the extent to which mutational bias in the production of genetic variation influences outcomes and pathways of adaptive molecular evolution. In the genomes of at least some vertebrate taxa, an important form of mutation bias involves changes at CpG dinucleotides: if the DNA nucleotide cytosine (C) is immediately 5' to guanine (G) on the same coding strand, then-depending on methylation status-point mutations at both sites occur at an elevated rate relative to mutations at non-CpG sites. Here, we examine experimental data from case studies in which it has been possible to identify the causative substitutions that are responsible for adaptive changes in the functional properties of vertebrate haemoglobin (Hb). Specifically, we examine the molecular basis of convergent increases in Hb-O₂ affinity in high-altitude birds. Using a dataset of experimentally verified, affinity-enhancing mutations in the Hbs of highland avian taxa, we tested whether causative changes are enriched for mutations at CpG dinucleotides relative to the frequency of CpG mutations among all possible missense mutations. The tests revealed that a disproportionate number of causative amino acid replacements were attributable to CpG mutations, suggesting that mutation bias can influence outcomes of molecular adaptation. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Gene Turnover and Diversification of the α- and β-Globin Gene Families in Sauropsid Vertebrates

The genes that encode the α- and β-chain subunits of vertebrate hemoglobin have served as a model system for elucidating general principles of gene family evolution, but little is known about patterns of evolution in amniotes other than mammals and birds. Here, we report a comparative genomic analysis of the α- and β-globin gene clusters in sauropsids (archosaurs and nonavian reptiles). The objectives were to characterize changes in the size and membership composition of the α- and β-globin gene families within and among the major sauropsid lineages, to reconstruct the evolutionary history of the sauropsid α- and β-globin genes, to resolve orthologous relationships, and to reconstruct evolutionary changes in the developmental regulation of gene expression. Our comparisons revealed contrasting patterns of evolution in the unlinked α- and β-globin gene clusters. In the α-globin gene cluster, which has remained in the ancestral chromosomal location, evolutionary changes in gene content are attributable to the differential retention of paralogous gene copies that were present in the common ancestor of tetrapods. In the β-globin gene cluster, which was translocated to a new chromosomal location, evolutionary changes in gene content are attributable to differential gene gains (via lineage-specific duplication events) and gene losses (via lineage-specific deletions and inactivations). Consequently, all major groups of amniotes possess unique repertoires of embryonic and postnatally expressed β-type globin genes that diversified independently in each lineage. These independently derived β-type globins descend from a pair of tandemly linked paralogs in the most recent common ancestor of sauropsids.

Allosteric mechanisms underlying the adaptive increase in hemoglobin-oxygen affinity of the bar-headed goose

The high blood-O2 affinity of the bar-headed goose (Anser indicus) is an integral component of the biochemical and physiological adaptations that allow this hypoxia-tolerant species to undertake migratory flights over the Himalayas. The high blood-O2 affinity of this species was originally attributed to a single amino acid substitution of the major hemoglobin (Hb) isoform, HbA, which was thought to destabilize the low-affinity T state, thereby shifting the T-R allosteric equilibrium towards the high-affinity R state. Surprisingly, this mechanistic hypothesis has never been addressed using native proteins purified from blood. Here, we report a detailed analysis of O2 equilibria and kinetics of native major HbA and minor HbD isoforms from bar-headed goose and greylag goose (Anser anser), a strictly lowland species, to identify and characterize the mechanistic basis for the adaptive change in Hb function. We find that HbA and HbD of bar-headed goose have consistently higher O2 affinities than those of the greylag goose. The corresponding Hb isoforms of the two species are equally responsive to physiological allosteric cofactors and have similar Bohr effects. Thermodynamic analyses of O2 equilibrium curves according to the two-state Monod-Wyman-Changeaux model revealed higher R-state O2 affinities in the bar-headed goose Hbs, associated with lower O2 dissociation rates, compared with the greylag goose. Conversely, the T state was not destabilized and the T-R allosteric equilibrium was unaltered in bar-headed goose Hbs. The physiological implication of these results is that increased R-state affinity allows for enhanced O2 saturation in the lungs during hypoxia, but without impairing O2 delivery to tissues.

Compensatory mutations and epistasis for protein function

Adaptive protein evolution may be facilitated by neutral amino acid mutations that confer no benefit when they first arise but which potentiate subsequent function-altering mutations via direct or indirect structural mechanisms. Theoretical and empirical results indicate that such compensatory interactions (intramolecular epistasis) can exert a strong influence on trajectories of protein evolution. For this reason, assessing the form and prevalence of intramolecular epistasis and characterizing biophysical mechanisms of compensatory interaction are important research goals at the nexus of structural biology and molecular evolution. Here I review recent insights derived from protein-engineering studies, and I describe an approach for identifying and characterizing mechanisms of epistasis that integrates experimental data on structure-function relationships with analyses of comparative sequence data.

Molecular basis of hemoglobin adaptation in the high-flying bar-headed goose

During the adaptive evolution of a particular trait, some selectively fixed mutations may be directly causative and others may be purely compensatory. The relative contribution of these two classes of mutation to adaptive phenotypic evolution depends on the form and prevalence of mutational pleiotropy. To investigate the nature of adaptive substitutions and their pleiotropic effects, we used a protein engineering approach to characterize the molecular basis of hemoglobin (Hb) adaptation in the high-flying bar-headed goose (Anser indicus), a hypoxia-tolerant species renowned for its trans-Himalayan migratory flights. To test the effects of observed substitutions on evolutionarily relevant genetic backgrounds, we synthesized all possible genotypic intermediates in the line of descent connecting the wildtype bar-headed goose genotype with the most recent common ancestor of bar-headed goose and its lowland relatives. Site-directed mutagenesis experiments revealed one major-effect mutation that significantly increased Hb-O2 affinity on all possible genetic backgrounds. Two other mutations exhibited smaller average effect sizes and less additivity across backgrounds. One of the latter mutations produced a concomitant increase in the autoxidation rate, a deleterious side-effect that was fully compensated by a second-site mutation at a spatially proximal residue. The experiments revealed three key insights: (i) subtle, localized structural changes can produce large functional effects; (ii) relative effect sizes of function-altering mutations may depend on the sequential order in which they occur; and (iii) compensation of deleterious pleiotropic effects may play an important role in the adaptive evolution of protein function.

Divergent and parallel routes of biochemical adaptation in high-altitude passerine birds from the Qinghai-Tibet Plateau

When different species experience similar selection pressures, the probability of evolving similar adaptive solutions may be influenced by legacies of evolutionary history, such as lineage-specific changes in genetic background. Here we test for adaptive convergence in hemoglobin (Hb) function among high-altitude passerine birds that are native to the Qinghai-Tibet Plateau, and we examine whether convergent increases in Hb-O₂ affinity have a similar molecular basis in different species. We documented that high-altitude parid and aegithalid species from the Qinghai-Tibet Plateau have evolved derived increases in Hb-O₂ affinity in comparison with their closest lowland relatives in East Asia. However, convergent increases in Hb-O₂ affinity and convergence in underlying functional mechanisms were seldom attributable to the same amino acid substitutions in different species. Using ancestral protein resurrection and site-directed mutagenesis, we experimentally confirmed two cases in which parallel substitutions contributed to convergent increases in Hb-O₂ affinity in codistributed high-altitude species. In one case involving the ground tit (Parus humilis) and gray-crested tit (Lophophanes dichrous), parallel amino acid replacements with affinity-enhancing effects were attributable to nonsynonymous substitutions at a CpG dinucleotide, suggesting a possible role for mutation bias in promoting recurrent changes at the same site. Overall, most altitude-related changes in Hb function were caused by divergent amino acid substitutions, and a select few were caused by parallel substitutions that produced similar phenotypic effects on the divergent genetic backgrounds of different species.

Hemoglobin polymerization via disulfide bond formation in the hypoxia-tolerant turtle Trachemys scripta: implications for antioxidant defense and O2 transport

The ability of many reptilian hemoglobins (Hbs) to form high-molecular weight polymers, albeit known for decades, has not been investigated in detail. Given that turtle Hbs often contain a high number of cysteine (Cys), potentially contributing to the red blood cell defense against reactive oxygen species, we have examined whether polymerization of Hb could occur via intermolecular disulfide bonds in red blood cells of freshwater turtle Trachemys scripta, a species that is highly tolerant of hypoxia and oxidative stress. We find that one of the two Hb isoforms of the hemolysate HbA is prone to polymerization in vitro into linear flexible chains of different size that are visible by electron microscopy but not the HbD isoform. Polymerization of purified HbA is favored by hydrogen peroxide, a main cellular reactive oxygen species and a thiol oxidant, and inhibited by thiol reduction and alkylation, indicating that HbA polymerization is due to disulfide bonds. By using mass spectrometry, we identify Cys5 of the α^A-subunit of HbA as specifically responsible for forming disulfide bonds between adjacent HbA tetramers. Polymerization of HbA does not affect oxygen affinity, cooperativity, and sensitivity to the allosteric cofactor ATP, indicating that HbA is still fully functional. Polymers also form in T. scripta blood after exposure to anoxia but not normoxia, indicating that they are of physiological relevance. Taken together, these results show that HbA polymers may form during oxidative stress and that Cys5α^A of HbA is a key element of the antioxidant capacity of turtle red blood cells.

O2 binding and CO2 sensitivity in haemoglobins of subterranean African mole rats

Inhabiting deep and sealed subterranean burrows, mole rats exhibit a remarkable suite of specializations, including eusociality (living in colonies with single breeding queens), extraordinary longevity, cancer immunity and poikilothermy, and extreme tolerance of hypoxia and hypercapnia. With little information available on adjustments in haemoglobin (Hb) function that may mitigate the impact of exogenous and endogenous constraints on the uptake and internal transport of O₂, we measured haematological characteristics, as well as Hb-O₂ binding affinity and sensitivity to pH (Bohr effect), CO₂, temperature and 2,3-diphosphoglycerate (DPG, the major allosteric modulator of Hb-O₂ affinity in red blood cells) in four social and two solitary species of African mole rats (family Bathyergidae) originating from different biomes and soil types across Central and Southern Africa. We found no consistent patterns in haematocrit (Hct) and blood and red cell DPG and Hb concentrations or in intrinsic Hb-O₂ affinity and its sensitivity to pH and DPG that correlate with burrowing, sociality and soil type. However, the results reveal low specific (pH independent) effects of CO₂ on Hb-O₂ affinity compared with humans that predictably safeguard pulmonary loading under hypoxic and hypercapnic burrow conditions. The O₂ binding characteristics are discussed in relation to available information on the primary structure of Hbs from adult and developmental stages of mammals subjected to hypoxia and hypercapnia and the molecular mechanisms underlying functional variation in rodent Hbs.

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.

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.

Interspecific variation and plasticity in hemoglobin nitrite reductase activity and its correlation with oxygen affinity in vertebrates

Deoxygenated hemoglobin (Hb) is a nitrite reductase that reduces naturally occurring nitrite to nitric oxide (NO), supplying physiological relevant NO under hypoxic conditions. The nitrite reductase activity is modulated by the allosteric equilibrium between the R and T structures of Hb that also determines oxygen affinity. In the present study we investigated nitrite reductase activity and O₂ affinity in Hbs from ten different vertebrate species under identical conditions to disclose interspecific variations and allow an extended test for a correlation between the rate constant for nitrite reduction and O₂ affinity. We also tested plastic changes in Hb properties via addition of T-structure-stabilizing organic phosphates (ATP and GTP). The decay in deoxyHb during its reaction with nitrite was exponential-like in ectotherms (Atlantic hagfish, carp, crucian carp, brown trout, rainbow trout, cane toad, Indian python and red-eared slider turtle), while it was sigmoid in mammals (harbor porpoise and rabbit). Typically, hypoxia-tolerant species showed a faster reaction than intolerant species. Addition of ATP and GTP decreased O₂ affinity and slowed the rate of nitrite reduction in a concentration-dependent manner. The initial second order rate constant of the deoxyHb-mediated nitrite reduction showed a strong curvilinear correlation with oxygen affinity among all ectothermic vertebrates, and the relationship also applied to plastic variations of Hb properties via organic phosphates. The relationship predicts high nitrite reductase activity in hypoxic tolerant species with high Hb-O₂ affinity and reveals that the decrease in erythrocyte ATP and/or GTP during acclimation to hypoxia in ectotherms increases the erythrocyte NO generating capacity.

Causes of molecular convergence and parallelism in protein evolution

To what extent is the convergent evolution of protein function attributable to convergent or parallel changes at the amino acid level? The mutations that contribute to adaptive protein evolution may represent a biased subset of all possible beneficial mutations owing to mutation bias and/or variation in the magnitude of deleterious pleiotropy. A key finding is that the fitness effects of amino acid mutations are often conditional on genetic background. This context dependence (epistasis) can reduce the probability of convergence and parallelism because it reduces the number of possible mutations that are unconditionally acceptable in divergent genetic backgrounds. Here, I review factors that influence the probability of replicated evolution at the molecular level.

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.

Contribution of a mutational hot spot to hemoglobin adaptation in high-altitude Andean house wrens

A key question in evolutionary genetics is why certain mutations or certain types of mutation make disproportionate contributions to adaptive phenotypic evolution. In principle, the preferential fixation of particular mutations could stem directly from variation in the underlying rate of mutation to function-altering alleles. However, the influence of mutation bias on the genetic architecture of phenotypic evolution is difficult to evaluate because data on rates of mutation to function-altering alleles are seldom available. Here, we report the discovery that a single point mutation at a highly mutable site in the β(A)-globin gene has contributed to an evolutionary change in hemoglobin (Hb) function in high-altitude Andean house wrens (Troglodytes aedon). Results of experiments on native Hb variants and engineered, recombinant Hb mutants demonstrate that a nonsynonymous mutation at a CpG dinucleotide in the β(A)-globin gene is responsible for an evolved difference in Hb-O2 affinity between high- and low-altitude house wren populations. Moreover, patterns of genomic differentiation between high- and low-altitude populations suggest that altitudinal differentiation in allele frequencies at the causal amino acid polymorphism reflects a history of spatially varying selection. The experimental results highlight the influence of mutation rate on the genetic basis of phenotypic evolution by demonstrating that a large-effect allele at a highly mutable CpG site has promoted physiological differentiation in blood O2 transport capacity between house wren populations that are native to different elevations.

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.

Differences in Hematological Traits between High- and Low-Altitude Lizards (Genus Phrynocephalus)

Phrynocephalus erythrurus (Lacertilia: Agamidae) is considered to be the highest living reptile in the world (about 4500-5000 m above sea level), whereas Phrynocephalus przewalskii inhabits low altitudes (about 1000-1500 m above sea level). Here, we report the differences in hematological traits between these two different Phrynocephalus species. Compared with P. przewalskii, the results indicated that P. erythrurus own higher oxygen carrying capacity by increasing red blood cell count (RBC), hemoglobin concentration ([Hb]) and hematocrit (Hct) and these elevations could promote oxygen carrying capacity without disadvantage of high viscosity. The lower partial pressure of oxygen in arterial blood (PaO2) of P. erythrurus did not cause the secondary alkalosis, which may be attributed to an efficient pulmonary system for oxygen (O2) loading. The elevated blood-O2 affinity in P. erythrurus may be achieved by increasing intrinsic O2 affinity of isoHbs and balancing the independent effects of potential heterotropic ligands. We detected one α-globin gene and three β-globin genes with 1 and 33 amino acid substitutions between these two species, respectively. Molecular dynamics simulation results showed that amino acids substitutions in β-globin chains could lead to the elimination of hydrogen bonds in T-state Hb models of P. erythrurus. Based on the present data, we suggest that P. erythrurus have evolved an efficient oxygen transport system under the unremitting hypobaric hypoxia.

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.