Gene duplication and the evolution of hemoglobin isoform differentiation in birds

Convergent High O2 Affinity but Distinct ATP-Mediated Allosteric Regulation of Hemoglobins in Oviparous and Viviparous Eremias Lizards from the Qinghai-Tibet Plateau

The functional adaptation and underlying molecular mechanisms of hemoglobins (Hbs) have primarily concentrated on mammals and birds, with few reports on reptiles. This study aimed to investigate the convergent and species-specific high-altitude adaptation mechanisms of Hbs in two Eremias lizards from the Qinghai-Tibet Plateau. The Hbs of high-altitude E. argus and E. multiocellata were characterized by significantly high overall and intrinsic Hb-O2 affinity compared to their low-altitude populations. Despite the similarly low Cl- sensitivities, the Hbs of high-altitude E. argus exhibited higher ATP sensitivity and ATP-dependent Bohr effects than that of E. multiocellata, which could facilitate O2 unloading in respiring tissues. Eremias lizards Hbs exhibited similarly low temperature sensitivities and relatively high Bohr effects at lower temperatures, which could help to stably deliver and release O2 to cold extremities at low temperatures. The oxygenation properties of Hbs in high-altitude populations might be attributed to varying ratios of β2/β1 globin and substitutions on the β2-type globin. Notably, the Asn12Ala in lowland E. argus could cause localized destabilization of the E-helix in the tetrameric Hb by elimination of hydrogen bonds, thereby resulting in its lowest O2 affinity. This study provides a valuable reference for the high-altitude adaptation mechanisms of hemoglobins in reptiles.

Divergent Reactivity of Hemoglobin Isoforms α2Aβ2 and α2Dβ2 from Meleagris gallopavo in Native and Fatty Acid-Modified States

Meleagris gallopavo (turkey) coexpresses distinct hemoglobin (Hb) isoforms, Hb α2Aβ2 (HbA) and α2Dβ2 (HbD), at a ratio of ∼3:1 (HbA:HbD). Herein, the reactivities of HbA and HbD were investigated in their native and free fatty acid (FFA)-modified states. Results indicated that HbD displays elevated autoxidation (kox) and an increased propensity to oxidize lipids in its reduced (oxy) and oxidized (met) forms. Interestingly, metHbD displayed less heme-globin cross-linking compared to HbA. Regarding FFA-modified Hb, we found that an FFA mixture and linoleic acid (LA) produced a bis-histidyl ferric (Bis-His) Hb species, decreasing the ability of Hb to oxidize lipids. Using molecular docking, we found LA to hydrogen bond with β Arg C6, found at the α1β2 interface, but the extent of Bis-His formation differs between HbA and HbD. Our findings suggest HbA displays elevated oxidative stability compared to HbD and that FFA may act as allosteric effectors of metHb.

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.

Resistance to glycation in the zebra finch: Mass spectrometry-based analysis and its perspectives for evolutionary studies of aging

In humans, hyperglycemia is associated with protein glycation, which may contribute to aging. Strikingly, birds usually outlive mammals of the same body mass, while exhibiting high plasma glucose levels. However, how birds succeed in escaping pro-aging effects of glycation remains unknown. Using a specific mass spectrometry-based approach in captive zebra finches of known age, we recorded high glycaemia values but no glycated hemoglobin form was found. Still, we showed that zebra finch hemoglobin can be glycated in vitro, albeit only to a limited extent compared to its human homologue. This may be due to peculiar structural features, as supported by the unusual presence of three different tetramer populations with balanced proportions and a still bound cofactor that could be inositol pentaphosphate. High levels of the glycated forms of zebra finch plasma serotransferrin, carbonic anhydrase 2, and albumin were measured. Glucose, age or body mass correlations with either plasma glycated proteins or hemoglobin isoforms suggest that those variables may be future molecular tools of choice to monitor glycation and its link with individual fitness. Our molecular advance may help determine how evolution succeeded in associating flying ability, high blood glucose and long lifespan in birds.

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.

Changes in Hemoglobin Isoforms in the Peripheral Blood of Rats with Experimental Posthemorrhagic Anemia

Six hemoglobin isoforms were identified in the peripheral blood of male Wistar rats by PAAG electrophoresis. Erythrocytes can be separated into 6 fractions by fractional centrifugation; the fractions differ by the content of reticulocytes and cells with fetal hemoglobin. Cells in each fraction contain one light and one heavy hemoglobin isoforms, and their combination is specific to each fraction. Changes in the hemoglobin isoforms in the whole blood in case of blood loss can be associated with changes in physical and chemical properties of circulating erythrocytes, as well as with the formation of various cell populations resulting from activated erythropoiesis.

Genomic Insights into the Origin and Evolution of Molluscan Red-Bloodedness in the Blood Clam Tegillarca granosa

Blood clams differ from their molluscan kins by exhibiting a unique red-blood (RB) phenotype; however, the genetic basis and biochemical machinery subserving this evolutionary innovation remain unclear. As a fundamental step toward resolving this mystery, we presented the first chromosome-level genome and comprehensive transcriptomes of the blood clam Tegillarca granosa for an integrated genomic, evolutionary, and functional analyses of clam RB phenotype. We identified blood clam-specific and expanded gene families, as well as gene pathways that are of RB relevant. Clam-specific RB-related hemoglobins (Hbs) showed close phylogenetic relationships with myoglobins (Mbs) of blood clam and other molluscs without the RB phenotype, indicating that clam-specific Hbs were likely evolutionarily derived from the Mb lineage. Strikingly, similar to vertebrate Hbs, blood clam Hbs were present in a form of gene cluster. Despite the convergent evolution of Hb clusters in blood clam and vertebrates, their Hb clusters may have originated from a single ancestral Mb-like gene as evidenced by gene phylogeny and synteny analysis. A full suite of enzyme-encoding genes for heme synthesis was identified in blood clam, with prominent expression in hemolymph and resembling those in vertebrates, suggesting a convergence of both RB-related Hb and heme functions in vertebrates and blood clam. RNA interference experiments confirmed the functional roles of Hbs and key enzyme of heme synthesis in the maintenance of clam RB phenotype. The high-quality genome assembly and comprehensive transcriptomes presented herein serve new genomic resources for the super-diverse phylum Mollusca, and provide deep insights into the origin and evolution of invertebrate RB.

Evolved increases in hemoglobin-oxygen affinity and the Bohr effect coincided with the aquatic specialization of penguins

In diving birds like penguins, physiologic considerations suggest that increased hemoglobin (Hb)-O₂ affinity may improve pulmonary O₂ extraction and enhance dive capacity. We integrated experimental tests on whole-blood and native Hbs of penguins with protein engineering experiments on reconstructed ancestral Hbs. The experiments involving ancestral protein resurrection enabled us to test for evolved changes in Hb function in the stem lineage of penguins after divergence from their closest nondiving relatives. We demonstrate that penguins evolved an increased Hb-O₂ affinity in conjunction with a greatly augmented Bohr effect (i.e., reduction in Hb-O₂ affinity at low pH) that should maximize pulmonary O₂ extraction without compromising O₂ delivery at systemic capillaries.

Dive capacities of air-breathing vertebrates are dictated by onboard O₂ stores, suggesting that physiologic specialization of diving birds such as penguins may have involved adaptive changes in convective O₂ transport. It has been hypothesized that increased hemoglobin (Hb)-O₂ affinity improves pulmonary O₂ extraction and enhances the capacity for breath-hold diving. To investigate evolved changes in Hb function associated with the aquatic specialization of penguins, we integrated comparative measurements of whole-blood and purified native Hb with protein engineering experiments based on site-directed mutagenesis. We reconstructed and resurrected ancestral Hb representing the common ancestor of penguins and the more ancient ancestor shared by penguins and their closest nondiving relatives (order Procellariiformes, which includes albatrosses, shearwaters, petrels, and storm petrels). These two ancestors bracket the phylogenetic interval in which penguin-specific changes in Hb function would have evolved. The experiments revealed that penguins evolved a derived increase in Hb-O₂ affinity and a greatly augmented Bohr effect (i.e., reduced Hb-O₂ affinity at low pH). Although an increased Hb-O₂ affinity reduces the gradient for O₂ diffusion from systemic capillaries to metabolizing cells, this can be compensated by a concomitant enhancement of the Bohr effect, thereby promoting O₂ unloading in acidified tissues. We suggest that the evolved increase in Hb-O₂ affinity in combination with the augmented Bohr effect maximizes both O₂ extraction from the lungs and O₂ unloading from the blood, allowing penguins to fully utilize their onboard O₂ stores and maximize underwater foraging time.

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.

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.

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.

Biochemical pedomorphosis and genetic assimilation in the hypoxia adaptation of Tibetan antelope

Developmental shifts in stage-specific gene expression can provide a ready mechanism of phenotypic change by altering the rate or timing of ontogenetic events. We found that the high-altitude Tibetan antelope (Panthelops hodgsonii) has evolved an adaptive increase in blood-O2 affinity by truncating the ancestral ontogeny of globin gene expression such that a high-affinity juvenile hemoglobin isoform (isoHb) completely supplants the lower-affinity isoHb that is expressed in the adult red blood cells of other bovids. This juvenilization of blood properties represents a canalization of an acclimatization response to hypoxia that has been well documented in adult goats and sheep. We also found the genomic mechanism underlying this regulatory isoHb switch, revealing how a reversible acclimatization response became genetically assimilated as an irreversible adaptation to chronic hypoxia.

Origin of complexity in haemoglobin evolution

Most proteins associate into multimeric complexes with specific architectures1,2, which often have functional properties such as cooperative ligand binding or allosteric regulation3. No detailed knowledge is available about how any multimer and its functions arose during evolution. Here we use ancestral protein reconstruction and biophysical assays to elucidate the origins of vertebrate haemoglobin, a heterotetramer of paralogous α- and β-subunits that mediates respiratory oxygen transport and exchange by cooperatively binding oxygen with moderate affinity. We show that modern haemoglobin evolved from an ancient monomer and characterize the historical 'missing link' through which the modern tetramer evolved-a noncooperative homodimer with high oxygen affinity that existed before the gene duplication that generated distinct α- and β-subunits. Reintroducing just two post-duplication historical substitutions into the ancestral protein is sufficient to cause strong tetramerization by creating favourable contacts with more ancient residues on the opposing subunit. These surface substitutions markedly reduce oxygen affinity and even confer cooperativity, because an ancient linkage between the oxygen binding site and the multimerization interface was already an intrinsic feature of the protein's structure. Our findings establish that evolution can produce new complex molecular structures and functions via simple genetic mechanisms that recruit existing biophysical features into higher-level architectures.

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.

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'.

Very high blood oxygen affinity and large Bohr shift differentiates the air-breathing siamese fighting fish (Betta splendens) from the closely related anabantoid the blue gourami (Trichopodus trichopterus)

The Siamese fighting fish, Betta splendens, and the blue gourami, Trichopodus trichopterus, are two closely related air-breathing anabantoid fishes. B. splendens is a sedentary facultative air breather frequenting often hypoxic waters, while T. trichopterus is a more active obligatory air-breather inhabiting better oxygenated waters. Despite their close taxonomic relationship, previous studies have shown inter-specific differences in both physiological and morphological plasticity. Consequently, we hypothesized that B. splendens would have the higher blood oxygen affinity characteristics typical of more hypoxia-tolerant fishes. Whole blood oxygen equilibrium curves were determined at 27 °C and pHs of 7.62, 7.44 and 7.25. At a pH of 7.62, the blood O₂ affinity (P₅₀) of B. splendens was just 2.9 mmHg, while that of T. trichopterus was ~5 times higher at 14.7 mmHg. There were no significant differences in P₅₀ between males and females in either species. The Bohr coefficient in B. splendens and T. trichopterus was -1.79 and - 0.83, respectively. B. splendens, unlike T. trichopterus, showed a large Root effect. Hills cooperatively coefficient, n, was ~2 in both species, indicating a significant binding cooperative between oxygen and hemoglobin. Collectively, these differences in blood O₂ transport characteristics in these two closely related species are likely correlated with the differing habitats in which they breed and inhabit as adults, as well as different activity levels. Finally, the very high blood O₂ affinity of B. splendens is not extraordinary among air-breathing fish, as revealed by a review of the literature of blood oxygen affinity in air-breathing fishes.

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.

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.

Molecular evolution and expression of oxygen transport genes in livebearing fishes (Poeciliidae) from hydrogen sulfide rich springs

Hydrogen sulfide (H2S) is a natural toxicant in some aquatic environments that has diverse molecular targets. It binds to oxygen transport proteins, rendering them non-functional by reducing oxygen-binding affinity. Hence, organisms permanently inhabiting H2S-rich environments are predicted to exhibit adaptive modifications to compensate for the reduced capacity to transport oxygen. We investigated 10 lineages of fish of the family Poeciliidae that have colonized freshwater springs rich in H2S—along with related lineages from non-sulfidic environments—to test hypotheses about the expression and evolution of oxygen transport genes in a phylogenetic context. We predicted shifts in the expression of and signatures of positive selection on oxygen transport genes upon colonization of H2S-rich habitats. Our analyses indicated significant shifts in gene expression for multiple hemoglobin genes in lineages that have colonized H2S-rich environments, and three hemoglobin genes exhibited relaxed selection in sulfidic compared to non-sulfidic lineages. However, neither changes in gene expression nor signatures of selection were consistent among all lineages in H2S-rich environments. Oxygen transport genes may consequently be predictable targets of selection during adaptation to sulfidic environments, but changes in gene expression and molecular evolution of oxygen transport genes in H2S-rich environments are not necessarily repeatable across replicated lineages.

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.

Migration-Selection Balance Drives Genetic Differentiation in Genes Associated with High-Altitude Function in the Speckled Teal (Anas flavirostris) in the Andes

Local adaptation frequently occurs across populations as a result of migration-selection balance between divergent selective pressures and gene flow associated with life in heterogeneous landscapes. Studying the effects of selection and gene flow on the adaptation process can be achieved in systems that have recently colonized extreme environments. This study utilizes an endemic South American duck species, the speckled teal (Anas flavirostris), which has both high- and low-altitude populations. High-altitude speckled teal (A. f. oxyptera) are locally adapted to the Andean environment and mostly allopatric from low-altitude birds (A. f. flavirostris); however, there is occasional gene flow across altitudinal gradients. In this study, we used next-generation sequencing to explore genetic patterns associated with high-altitude adaptation in speckled teal populations, as well as the extent to which the balance between selection and migration have affected genetic architecture. We identified a set of loci with allele frequencies strongly correlated with altitude, including those involved in the insulin-like signaling pathway, bone morphogenesis, oxidative phosphorylation, responders to hypoxia-induced DNA damage, and feedback loops to the hypoxia-inducible factor pathway. These same outlier loci were found to have depressed gene flow estimates, as well as being highly concentrated on the Z-chromosome. Our results suggest a multifactorial response to life at high altitudes through an array of interconnected pathways that are likely under positive selection and whose genetic components seem to be providing an effective genomic barrier to interbreeding, potentially functioning as an avenue for population divergence and speciation.

Molecular basis of α-thalassemia

α-Thalassemia is an inherited, autosomal recessive, disorder characterized by a microcytic hypochromic anemia. It is one of the most common monogenic gene disorders in the world population. The clinical severity varies from almost asymptomatic, to mild microcytic hypochromic, and to a lethal hemolytic condition, called Hb Bart's Hydrops Foetalis Syndrome. The molecular basis are usually deletions and less frequently, point mutations affecting the expression of one or more of the duplicated α-genes. The clinical variation and increase in disease severity is directly related to the decreased expression of one, two, three or four copies of the α-globin genes. Deletions and point mutations in the α-globin genes and their regulatory elements have been studied extensively in carriers and patients and these studies have given insight into the α-globin genes are regulated. By looking at naturally occurring deletions and point mutations, our knowledge of globin-gene regulation and expression will continue to increase and will lead to new targets of therapy.

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.

Hemoglobin-oxygen affinity in high-altitude vertebrates: is there evidence for an adaptive trend?

In air-breathing vertebrates at high altitude, fine-tuned adjustments in hemoglobin (Hb)-O2 affinity provide an energetically efficient means of mitigating the effects of arterial hypoxemia. However, it is not always clear whether an increased or decreased Hb-O2 affinity should be expected to improve tissue O2 delivery under different degrees of hypoxia, due to the inherent trade-off between arterial O2 loading and peripheral O2 unloading. Theoretical results indicate that the optimal Hb-O2 affinity varies as a non-linear function of environmental O2 availability, and the threshold elevation at which an increased Hb-O2 affinity becomes advantageous depends on the magnitude of diffusion limitation (the extent to which O2 equilibration at the blood-gas interface is limited by the kinetics of O2 exchange). This body of theory provides a framework for interpreting the possible adaptive significance of evolved changes in Hb-O2 affinity in vertebrates that have colonized high-altitude environments. To evaluate the evidence for an empirical generalization and to test theoretical predictions, I synthesized comparative data in a phylogenetic framework to assess the strength of the relationship between Hb-O2 affinity and native elevation in mammals and birds. Evidence for a general trend in mammals is equivocal, but there is a remarkably strong positive relationship between Hb-O2 affinity and native elevation in birds. Evolved changes in Hb function in high-altitude birds provide one of the most compelling examples of convergent biochemical adaptation in vertebrates.

A preliminary study on the antibacterial mechanism of Tegillarca granosa hemoglobin by derived peptides and peroxidase activity

The blood clam, Tegillarca granosa, is one of the few bivalve molluscs containing hemoglobin (Hb). In the present study, we purified two types of T. granosa hemoglobin, Tg-HbI and Tg-HbII, using size exclusion chromatography and measured their antibacterial and peroxidase activities. We also tested antibacterial activities of peptides prepared by trypsin digestion of purified Tg-Hb and reversed-phase high-performance liquid chromatography purification. Purified Tg-HbI and Tg-HbII showed antibacterial activity against Escherichia coli, Pseudomonas putida, Bacillus subtilis, and Bacillus firmus, with differences in minimal inhibitory concentrations (MICs), but lacked antibacterial activity against Vibrio alginolyticus, Vibrio parahaemolyticus, Vibrio harveyi and Staphylococcus aureus. In contrast, 7 Tg-Hb derived peptides exhibited varying degrees of antibacterial activity against V. alginolyticus (MICs: 12-200 μg/ml), V. parahaemolyticus (11-100 μg/ml) and V. harveyi (1-200 μg/ml). The antibacterial activity of Hb derived peptides was confirmed by fluorescence microscopy. In addition, peroxidase activity was detected in Tg-HbI and Tg-HbII. The results indicated that in addition to functioning as a respiratory protein T. granosa hemoglobins likely play a role in host antibacterial defense probably via a peroxidase activity of native molecules and some internal peptides released from the proteins.

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.

Top-down-assisted bottom-up method for homologous protein sequencing: hemoglobin from 33 bird species

Ticks are vectors for disease transmission because they are indiscriminant in their feeding on multiple vertebrate hosts, transmitting pathogens between their hosts. Identifying the hosts on which ticks have fed is important for disease prevention and intervention. We have previously shown that hemoglobin (Hb) remnants from a host on which a tick fed can be used to reveal the host's identity. For the present research, blood was collected from 33 bird species that are common in the U.S. as hosts for ticks but that have unknown Hb sequences. A top-down-assisted bottom-up mass spectrometry approach with a customized searching database, based on variability in known bird hemoglobin sequences, has been devised to facilitate fast and complete sequencing of hemoglobin from birds with unknown sequences. These hemoglobin sequences will be added to a hemoglobin database and used for tick host identification. The general approach has the potential to sequence any set of homologous proteins completely in a rapid manner. Graphical Abstract ᅟ.

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.

Quantitative proteomics provides new insights into chicken eggshell matrix protein functions during the primary events of mineralisation and the active calcification phase

Eggshell is a bioceramic composed of 95% calcium carbonate mineral and 3.5% organic matrix. Its structural organisation is controlled by its organic matrix. We have used quantitative proteomics to study four key stages of shell mineralisation: 1) widespread deposition of amorphous calcium carbonate (ACC), 2) ACC transformation into crystalline calcite aggregates, 3) formation of larger calcite crystal units and 4) development of a columnar structure with preferential calcite crystal orientation. This approach explored the distribution of 216 shell matrix proteins found at the four stages. Variations in abundance according to these calcification events were observed for 175 proteins. A putative function related to the mineralisation process was predicted by bioinformatics for 77 of them and was further characterised. We confirmed the important role of lysozyme, ovotransferrin, ovocleidin-17 and ovocleidin-116 for shell calcification process, characterised major calcium binding proteins (EDIL3, ALB, MFGE8, NUCB2), and described novel proteoglycans core proteins (GPC4, HAPLN3). We suggest that OVAL and OC-17 play a role in the stabilisation of ACC. Finally, we report proteins involved in the regulation of proteins driving the mineralisation. They correspond to numerous molecular chaperones including CLU, PPIB and OCX21, protease and protease inhibitors including OVM and CST3, and regulators of phosphorylation.

Genetic approaches in comparative and evolutionary physiology

Whole animal physiological performance is highly polygenic and highly plastic, and the same is generally true for the many subordinate traits that underlie performance capacities. Quantitative genetics, therefore, provides an appropriate framework for the analysis of physiological phenotypes and can be used to infer the microevolutionary processes that have shaped patterns of trait variation within and among species. In cases where specific genes are known to contribute to variation in physiological traits, analyses of intraspecific polymorphism and interspecific divergence can reveal molecular mechanisms of functional evolution and can provide insights into the possible adaptive significance of observed sequence changes. In this review, we explain how the tools and theory of quantitative genetics, population genetics, and molecular evolution can inform our understanding of mechanism and process in physiological evolution. For example, lab-based studies of polygenic inheritance can be integrated with field-based studies of trait variation and survivorship to measure selection in the wild, thereby providing direct insights into the adaptive significance of physiological variation. Analyses of quantitative genetic variation in selection experiments can be used to probe interrelationships among traits and the genetic basis of physiological trade-offs and constraints. We review approaches for characterizing the genetic architecture of physiological traits, including linkage mapping and association mapping, and systems approaches for dissecting intermediary steps in the chain of causation between genotype and phenotype. We also discuss the promise and limitations of population genomic approaches for inferring adaptation at specific loci. We end by highlighting the role of organismal physiology in the functional synthesis of evolutionary biology.

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.

High affinity and temperature sensitivity of blood oxygen binding in Pangasianodon hypophthalmus due to lack of chloride-hemoglobin allosteric interaction

Air-breathing fishes represent interesting organisms in terms of understanding the physiological changes associated with the terrestrialization of vertebrates, and, further, are of great socio-economic importance for aquaculture in Southeast Asia. To understand how environmental factors, such as high temperature, affect O2 transport in air-breathing fishes, this study assessed the effects of temperature on O2 binding of blood and Hb in the economically important air-breathing fish Pangasianodon hypophthalmus. To determine blood O2 binding properties, blood was drawn from resting cannulated fishes and O2 binding curves made at 25°C and 35°C. To determine the allosteric regulation and thermodynamics of Hb O2 binding, Hb was purified, and O2 equilibria were recorded at five temperatures in the absence and presence of ATP and Cl(-). Whole blood had a high O2 affinity (O2 tension at half saturation P50 = 4.6 mmHg at extracellular pH 7.6 and 25°C), a high temperature sensitivity of O2 binding (apparent heat of oxygenation ΔH(app) = -28.3 kcal/mol), and lacked a Root effect. Further, the data on Hb revealed weak ATP binding and a complete lack of Cl(-) binding to Hb, which, in part, explains the high O2 affinity and high temperature sensitivity of blood O2 binding. This study demonstrates how a potent mechanism for increasing O2 affinity is linked to increased temperature sensitivity of O2 transport and provides a basic framework for a better understanding of how hypoxia-adapted species will react to increasing temperatures.

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.

Ancestral duplications and highly dynamic opsin gene evolution in percomorph fishes

Single-gene and whole-genome duplications are important evolutionary mechanisms that contribute to biological diversification by launching new genetic raw material. For example, the evolution of animal vision is tightly linked to the expansion of the opsin gene family encoding light-absorbing visual pigments. In teleost fishes, the most species-rich vertebrate group, opsins are particularly diverse and key to the successful colonization of habitats ranging from the bioluminescence-biased but basically dark deep sea to clear mountain streams. In this study, we report a previously unnoticed duplication of the violet-blue short wavelength-sensitive 2 (SWS2) opsin, which coincides with the radiation of highly diverse percomorph fishes, permitting us to reinterpret the evolution of this gene family. The inspection of close to 100 fish genomes revealed that, triggered by frequent gene conversion between duplicates, the evolutionary history of SWS2 is rather complex and difficult to predict. Coincidentally, we also report potential cases of gene resurrection in vertebrate opsins, whereby pseudogenized genes were found to convert with their functional paralogs. We then identify multiple novel amino acid substitutions that are likely to have contributed to the adaptive differentiation between SWS2 copies. Finally, using the dusky dottyback Pseudochromis fuscus, we show that the newly discovered SWS2A duplicates can contribute to visual adaptation in two ways: by gaining sensitivities to different wavelengths of light and by being differentially expressed between ontogenetic stages. Thus, our study highlights the importance of comparative approaches in gaining a comprehensive view of the dynamics underlying gene family evolution and ultimately, animal diversification.

Comparative genomics reveals insights into avian genome evolution and adaptation

Birds are the most species-rich class of tetrapod vertebrates and have wide relevance across many research fields. We explored bird macroevolution using full genomes from 48 avian species representing all major extant clades. The avian genome is principally characterized by its constrained size, which predominantly arose because of lineage-specific erosion of repetitive elements, large segmental deletions, and gene loss. Avian genomes furthermore show a remarkably high degree of evolutionary stasis at the levels of nucleotide sequence, gene synteny, and chromosomal structure. Despite this pattern of conservation, we detected many non-neutral evolutionary changes in protein-coding genes and noncoding regions. These analyses reveal that pan-avian genomic diversity covaries with adaptations to different lifestyles and convergent evolution of traits.

Epistasis constrains mutational pathways of hemoglobin adaptation in high-altitude pikas

A fundamental question in evolutionary genetics concerns the roles of mutational pleiotropy and epistasis in shaping trajectories of protein evolution. This question can be addressed most directly by using site-directed mutagenesis to explore the mutational landscape of protein function in experimentally defined regions of sequence space. Here, we evaluate how pleiotropic trade-offs and epistatic interactions influence the accessibility of alternative mutational pathways during the adaptive evolution of hemoglobin (Hb) function in high-altitude pikas (Mammalia: Lagomorpha). By combining ancestral protein resurrection with a combinatorial protein-engineering approach, we examined the functional effects of sequential mutational steps in all possible pathways that produced an increased Hb–O₂ affinity. These experiments revealed that the effects of mutations on Hb–O₂ affinity are highly dependent on the temporal order in which they occur: Each of three β-chain substitutions produced a significant increase in Hb–O₂ affinity on the ancestral genetic background, but two of these substitutions produced opposite effects when they occurred as later steps in the pathway. The experiments revealed pervasive epistasis for Hb–O₂ affinity, but affinity-altering mutations produced no significant pleiotropic trade-offs. These results provide insights into the properties of adaptive substitutions in naturally evolved proteins and suggest that the accessibility of alternative mutational pathways may be more strongly constrained by sign epistasis for positively selected biochemical phenotypes than by antagonistic pleiotropy.

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 conserved Phe GH5 of importance for hemoglobin intersubunit contact is mutated in gadoid fish

Background

Functionality of the tetrameric hemoglobin molecule seems to be determined by a few amino acids located in key positions. Oxygen binding encompasses structural changes at the interfaces between the α1β2 and α2β1 dimers, but also subunit interactions are important for the oxygen binding affinity and stability. The latter packing contacts include the conserved Arg B12 interacting with Phe GH5, which is replaced by Leu and Tyr in the α^A and α^D chains, respectively, of birds and reptiles.

Results

Searching all known hemoglobins from a variety of gnathostome species (jawed vertebrates) revealed the almost invariant Arg B12 coded by the AGG triplet positioned at an exon-intron boundary. Rare substitutions of Arg B12 in the gnathostome β globins were found in pig, tree shrew and scaled reptiles. Phe GH5 is also highly conserved in the β globins, except for the Leu replacement in the β1 globin of five marine gadoid species, gilthead seabream and the Comoran coelacanth, while Cys and Ile were found in burbot and yellow croaker, respectively. Atlantic cod β1 globin showed a Leu/Met polymorphism at position GH5 dominated by the Met variant in northwest-Atlantic populations that was rarely found in northeast-Atlantic cod. Site-specific analyses identified six consensus codons under positive selection, including 122β(GH5), indicating that the amino acid changes identified at this position may offer an adaptive advantage. In fact, computational mutation analysis showed that the replacement of Phe GH5 with Leu or Cys decreased the number of van der Waals contacts essentially in the deoxy form that probably causes a slight increase in the oxygen binding affinity.

Conclusions

The almost invariant Arg B12 and the AGG codon seem to be important for the packing contacts and pre-mRNA processing, respectively, but the rare mutations identified might be beneficial. The Leu122β1(GH5)Met and Met55β1(D6)Val polymorphisms in Atlantic cod hemoglobin modify the intradimer contacts B12-GH5 and H2-D6, while amino acid replacements at these positions in avian hemoglobin seem to be evolutionary adaptive in air-breathing vertebrates. The results support the theory that adaptive changes in hemoglobin functions are caused by a few substitutions at key positions.

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.

Hemoglobin isoform differentiation and allosteric regulation of oxygen binding in the turtle, Trachemys scripta

When freshwater turtles acclimatize to winter hibernation, there is a gradual transition from aerobic to anaerobic metabolism, which may require adjustments of blood O2 transport before turtles become anoxic. Here, we report the effects of protons, anionic cofactors, and temperature on the O2-binding properties of isolated hemoglobin (Hb) isoforms, HbA and HbD, in the turtle Trachemys scripta. We determined the primary structures of the constituent subunits of the two Hb isoforms, and we related the measured functional properties to differences in O2 affinity between untreated hemolysates from turtles that were acclimated to normoxia and anoxia. Our data show that HbD has a consistently higher O2 affinity compared with HbA, whereas Bohr and temperature effects, as well as thiol reactivity, are similar. Although sequence data show amino acid substitutions at two known β-chain ATP-binding site positions, we find high ATP affinities for both Hb isoforms, suggesting an alternative and stronger binding site for ATP. The high ATP affinities indicate that, although ATP levels decrease in red blood cells of turtles acclimating to anoxia, the O2 affinity would remain largely unchanged, as confirmed by O2-binding measurements of untreated hemolysates from normoxic and anoxic turtles. Thus, the increase in blood-O2 affinity that accompanies winter acclimation is mainly attributable to a decrease in temperature rather than in concentrations of organic phosphates. This is the first extensive study on freshwater turtle Hb isoforms, providing molecular evidence for adaptive changes in O2 transport associated with acclimation to severe hypoxia.

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

In contrast to other vertebrate hemoglobins (Hbs) whose high intrinsic O2 affinities are reduced by red cell allosteric effectors (mainly protons, CO2, organic phosphates, and chloride ions), crocodilian Hbs exhibit low sensitivity to organic phosphates and high sensitivity to bicarbonate (HCO3(-)), which is believed to augment Hb-O2 unloading during diving and postprandial alkaline tides when blood HCO3(-) levels and metabolic rates increase. Examination of α- and β-globin amino acid sequences of dwarf caiman (Paleosuchus palpebrosus) revealed a unique combination of substitutions at key effector binding sites compared with other vertebrate and crocodilian Hbs: β82Lys→Gln, β143His→Val, and β146His→Tyr. These substitutions delete positive charges and, along with other distinctive changes in residue charge and polarity, may be expected to disrupt allosteric regulation of Hb-O2 affinity. Strikingly, however, P. palpebrosus Hb shows a strong Bohr effect, and marked deoxygenation-linked binding of organic phosphates (ATP and DPG) and CO2 as carbamate (contrasting with HCO3(-) binding in other crocodilians). Unlike other Hbs, it polymerizes to large complexes in the oxygenated state. The highly unusual properties of P. palpebrosus Hb align with a high content of His residues (potential sites for oxygenation-linked proton binding) and distinctive surface Cys residues that may form intermolecular disulfide bridges upon polymerization. On the basis of its singular properties, P. palpebrosus Hb provides a unique opportunity for studies on structure-function coupling and the evolution of compensatory mechanisms for maintaining tissue O2 delivery in Hbs that lack conventional effector-binding residues.

High thermal sensitivity of blood enhances oxygen delivery in the high-flying bar-headed goose

The bar-headed goose (Anser indicus) crosses the Himalaya twice a year at altitudes where oxygen (O2) levels are less than half those at sea level and temperatures are below -20°C. Although it has been known for over three decades that the major hemoglobin (Hb) component of bar-headed geese has an increased affinity for O2, enhancing O2 uptake, the effects of temperature and interactions between temperature and pH on bar-headed goose Hb-O2 affinity have not previously been determined. An increase in breathing of the hypoxic and extremely cold air experienced by a bar-headed goose at altitude (due to the enhanced hypoxic ventilatory response in this species) could result in both reduced temperature and reduced levels of CO2 at the blood-gas interface in the lungs, enhancing O2 loading. In addition, given the strenuous nature of flapping flight, particularly in thin air, blood leaving the exercising muscle should be warm and acidotic, facilitating O2 unloading. To explore the possibility that features of blood biochemistry in this species could further enhance O2 delivery, we determined the P50 (the partial pressure of O2 at which Hb is 50% saturated) of whole blood from bar-headed geese under conditions of varying temperature and [CO2]. We found that blood-O2 affinity was highly temperature sensitive in bar-headed geese compared with other birds and mammals. Based on our analysis, temperature and pH effects acting on blood-O2 affinity (cold alkalotic lungs and warm acidotic muscle) could increase O2 delivery by twofold during sustained flapping flight at high altitudes compared with what would be delivered by blood at constant temperature and pH.