The alphaD-globin gene originated via duplication of an embryonic alpha-like globin gene in the ancestor of tetrapod vertebrates

Structural and functional characterization of haemoglobin genes in Labeo catla: Insights into hypoxic adaptation and survival

Haemoglobin (HB) protein comprises four subunits: two identical α-subunits (HBA) and two identical β-subunits (HBB), encoded by the HBA and HBB genes. In this investigation, 5'/3' RACE PCR (Rapid Amplification of cDNA Ends) was used to obtain complete coding sequences (CDSs) of both the genes from farmed Labeo catla. The resulting CDSs were 432 base pairs and 447 base pairs for HBA and HBB, respectively, corresponding to 143 and 148 amino acids. Phylogenetic analysis revealed close relationships with other cyprinids, with Labeo rohita being the closest relative. Functional analysis and protein structure prediction were conducted using bioinformatics tools. Expression profiling of both genes was checked in various tissues under control (C) and hypoxic (H) conditions. Notably, under hypoxia, HBA and HBB genes were significantly upregulated (P < 0.05) initially, followed by a return to normal expression levels. Similar trends were observed for Hif1α (Hypoxia-inducible factor one alpha) and EPO (Erythropoietin) genes. Additionally, haematological indices also significantly increased corresponding to the gene expressions. However, with the decrease in the expression of these genes an onset of mortality was observed in the hypoxia (H) treated groups. The results of the current study explored the role of haemoglobin genes in adaptation to the hypoxic condition.

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.

Whole genome sequencing of the fast-swimming Southern bluefin tuna (Thunnus maccoyii)

The economically important Southern bluefin tuna (Thunnus maccoyii) is a world-famous fast-swimming fish, but its genomic information is limited. Here, we performed whole genome sequencing and assembled a draft genome for Southern bluefin tuna, aiming to generate useful genetic data for comparative functional prediction. The final genome assembly is 806.54 Mb, with scaffold and contig N50 values of 3.31 Mb and 67.38 kb, respectively. Genome completeness was evaluated to be 95.8%. The assembled genome contained 23,403 protein-coding genes and 236.1 Mb of repeat sequences (accounting for 29.27% of the entire assembly). Comparative genomics analyses of this fast-swimming tuna revealed that it had more than twice as many hemoglobin genes (18) as other relatively slow-moving fishes (such as seahorse, sunfish, and tongue sole). These hemoglobin genes are mainly localized in two big clusters (termed as "MNˮ and "LAˮ respectively), which is consistent with other reported fishes. However, Thr39 of beta-hemoglobin in the MN cluster, conserved in other fishes, was mutated as cysteine in tunas including the Southern bluefin tuna. Since hemoglobins are reported to transport oxygen efficiently for aerobic respiration, our genomic data suggest that both high copy numbers of hemoglobin genes and an adjusted function of the beta-hemoglobin may support the fast-swimming activity of tunas. In summary, we produced a primary genome assembly and predicted hemoglobin-related roles for the fast-swimming Southern bluefin tuna.

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.

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.

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.

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.

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.

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.

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.

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.

Distinctive patterns of evolution of the δ-globin gene (HBD) in primates

In most vertebrates, hemoglobin (Hb) is a heterotetramer composed of two dissimilar globin chains, which change during development according to the patterns of expression of α- and β-globin family members. In placental mammals, the β-globin cluster includes three early-expressed genes, ε(HBE)-γ(HBG)-ψβ(HBBP1), and the late expressed genes, δ (HBD) and β (HBB). While HBB encodes the major adult β-globin chain, HBD is weakly expressed or totally silent. Paradoxically, in human populations HBD shows high levels of conservation typical of genes under strong evolutionary constraints, possibly due to a regulatory role in the fetal-to-adult switch unique of Anthropoid primates. In this study, we have performed a comprehensive phylogenetic and comparative analysis of the two adult β-like globin genes in a set of diverse mammalian taxa, focusing on the evolution and functional divergence of HBD in primates. Our analysis revealed that anthropoids are an exception to a general pattern of concerted evolution in placental mammals, showing a high level of sequence conservation at HBD, less frequent and shorter gene conversion events. Moreover, this lineage is unique in the retention of a functional GATA-1 motif, known to be involved in the control of the developmental expression of the β-like globin genes. We further show that not only the mode but also the rate of evolution of the δ-globin gene in higher primates are strictly associated with the fetal/adult β-cluster developmental switch. To gain further insight into the possible functional constraints that have been shaping the evolutionary history of HBD in primates, we calculated dN/dS (ω) ratios under alternative models of gene evolution. Although our results indicate that HBD might have experienced different selective pressures throughout primate evolution, as shown by different ω values between apes and Old World Monkeys + New World Monkeys (0.06 versus 0.43, respectively), these estimates corroborated a constrained evolution for HBD in Anthropoid lineages, which is unlikely to be related to protein function. Collectively, these findings suggest that sequence change at the δ-globin gene has been under strong selective constraints over 65 Myr of primate evolution, likely due to a regulatory role in ontogenic switches of gene expression.

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.

Evolution of α- and β-globin genes and their regulatory systems in light of the hypothesis of domain organization of the genome

The α- and β-globin gene domains are a traditional model for study of the domain organization of the eucaryotic genome because these genes encode hemoglobin, a physiologically important protein. The α-globin and β-globin gene domains are organized in completely different ways, while the expression of globin genes is tightly coordinated, which makes it extremely interesting to study the origin of these genes and the evolution of their regulatory systems. In this review, the organization of the α- and β-globin gene domains and their genomic environment in different taxonomic groups are comparatively analyzed. A new hypothesis of possible evolutionary pathways for segregated α- and β-globin gene domains of warm-blooded animals is proposed.

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.

Function and evolution of vertebrate globins

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

Repeated evolution of chimeric fusion genes in the β-globin gene family of laurasiatherian mammals

The evolutionary fate of chimeric fusion genes may be strongly influenced by their recombinational mode of origin and the nature of functional divergence between the parental genes. In the β-globin gene family of placental mammals, the two postnatally expressed δ- and β-globin genes (HBD and HBB, respectively) have a propensity for recombinational exchange via gene conversion and unequal crossing-over. In the latter case, there are good reasons to expect differences in retention rates for the reciprocal HBB/HBD and HBD/HBB fusion genes due to thalassemia pathologies associated with the HBD/HBB "Lepore" deletion mutant in humans. Here, we report a comparative genomic analysis of the mammalian β-globin gene cluster, which revealed that chimeric HBB/HBD fusion genes originated independently in four separate lineages of laurasiatherian mammals: Eulipotyphlans (shrews, moles, and hedgehogs), carnivores, microchiropteran bats, and cetaceans. In cases where an independently derived "anti-Lepore" duplication mutant has become fixed, the parental HBD and/or HBB genes have typically been inactivated or deleted, so that the newly created HBB/HBD fusion gene is primarily responsible for synthesizing the β-type subunits of adult and fetal hemoglobin (Hb). Contrary to conventional wisdom that the HBD gene is a vestigial relict that is typically inactivated or expressed at negligible levels, we show that HBD-like genes often encode a substantial fraction (20-100%) of β-chain Hbs in laurasiatherian taxa. Our results indicate that the ascendancy or resuscitation of genes with HBD-like coding sequence requires the secondary acquisition of HBB-like promoter sequence via unequal crossing-over or interparalog gene conversion.

Relaxed functional constraints on triplicate α-globin gene in the bank vole suggest a different evolutionary history from other rodents

Gene duplication plays an important role in the origin of evolutionary novelties, but the mechanisms responsible for the retention and functional divergence of the duplicated copy are not fully understood. The α-globin genes provide an example of a gene family with different numbers of gene duplicates among rodents. Whereas Rattus and Peromyscus each have three adult α-globin genes (HBA-T1, HBA-T2 and HBA-T3), Mus has only two copies. High rates of amino acid evolution in the independently derived HBA-T3 genes of Peromyscus and Rattus have been attributed to positive selection. Using RACE PCR, reverse transcription-PCR (RT-PCR) and RNA-seq, we show that another rodent, the bank vole Clethrionomys glareolus, possesses three transcriptionally active α-globin genes. The bank vole HBA-T3 gene is distinguished from each HBA-T1 and HBA-T2 by 20 amino acids and is transcribed 23- and 4-fold lower than HBA-T1 and HBA-T2, respectively. Polypeptides corresponding to all three genes are detected by electrophoresis, demonstrating that the translated products of HBA-T3 are present in adult erythrocytes. Patterns of codon substitution and the presence of low-frequency null alleles suggest a postduplication relaxation of purifying selection on bank vole HBA-T3.

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.

Evolution of hemoglobin and its genes

Insights into the evolution of hemoglobins and their genes are an abundant source of ideas regarding hemoglobin function and regulation of globin gene expression. This article presents the multiple genes and gene families encoding human globins, summarizes major events in the evolution of the hemoglobin gene clusters, and discusses how these studies provide insights into regulation of globin genes. Although the genes in and around the α-like globin gene complex are relatively stable, the β-like globin gene clusters are more dynamic, showing evidence of transposition to a new locus and frequent lineage-specific expansions and deletions. The cis-regulatory modules controlling levels and timing of gene expression are a mix of conserved and lineage-specific DNA, perhaps reflecting evolutionary constraint on core regulatory functions shared broadly in mammals and adaptive fine-tuning in different orders of mammals.

Gene duplication and the evolution of hemoglobin isoform differentiation in birds

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

Developmental regulation of hemoglobin synthesis in the green anole lizard Anolis carolinensis

Tetrapod vertebrates possess multiple α- and β-like globin genes that are ontogenetically regulated, such that functionally distinct hemoglobin (Hb) isoforms are synthesized during different stages of development. The α- and β-like globin genes of amphibians, birds and mammals are differentially expressed during embryonic development and postnatal life, but little is known about the developmental regulation of globin gene expression in non-avian reptiles. Here we report an investigation into the developmental regulation of Hb synthesis in the green anole lizard Anolis carolinensis. We tested two hypotheses derived from comparative genomic studies of the globin gene clusters in tetrapod vertebrates. First, we tested whether the product of the Anolis α(D)-globin gene is incorporated into embryonic Hb, thereby performing the role that would normally be performed by the embyronic α(E)-globin gene (which has been deleted from the green anole genome). Second, we tested whether two 'lizard-specific' β-globin paralogs have independently evolved a division of labor between an early-expressed embryonic gene and a later-expressed adult gene. Results of a proteomic analysis revealed that α- and β-like globin genes of the anole are differentially expressed during embryonic development. However, the same repertoire of α- and β-chain Hb isoforms was expressed during all stages of development and postnatal life, and the ontogenetic shifts in isoform composition were relatively subtle. In contrast to the pattern that has been documented in other tetrapod vertebrates, it appears that the developmental regulation of Hb synthesis in the green anole lizard does not involve discrete, stage-specific switches in gene activation and gene silencing.

Phylogenetic diversification of the globin gene superfamily in chordates

Phylogenetic reconstructions provide a means of inferring the branching relationships among members of multigene families that have diversified via successive rounds of gene duplication and divergence. Such reconstructions can illuminate the pathways by which particular expression patterns and protein functions evolved. For example, phylogenetic analyses can reveal cases in which similar expression patterns or functional properties evolved independently in different lineages, either through convergence, parallelism, or evolutionary reversals. The purpose of this article is to provide a robust phylogenetic framework for interpreting experimental data and for generating hypotheses about the functional evolution of globin proteins in chordate animals. To do this, we present a consensus phylogeny of the chordate globin gene superfamily. We document the relative roles of gene duplication and whole-genome duplication in fueling the functional diversification of vertebrate globins, and we unravel patterns of shared ancestry among globin genes from representatives of the three chordate subphyla (Craniata, Urochordata, and Cephalochordata). Our results demonstrate the value of integrating phylogenetic analyses with genomic analyses of conserved synteny to infer the duplicative origins and evolutionary histories of globin genes. We also discuss a number of case studies that illustrate the importance of phylogenetic information when making inferences about the evolution of globin gene expression and protein function. Finally, we discuss why the globin gene superfamily presents special challenges for phylogenetic analysis, and we describe methodological approaches that can be used to meet those challenges.

Indefinite noncooperative self-association of chicken deoxy hemoglobin D

The minor tetrameric hemoglobin (Hb), Hb D, of chicken red blood cells self-associates upon deoxygenation. This self-association enhances the cooperativity of oxygen binding. The maximal Hill coefficient is greater than 4 at high Hb concentrations. Previous measurements at low Hb concentrations were consistent with a monomer-to-dimer equilibrium and an association constant of ∼1.3-1.6 × 10(4) M(-1). Here, the Hb tetramer is considered as the monomer. However, new results indicate that the association extends beyond the dimer. We show by combination of Hb oligomer modeling and sedimentation velocity analyses that the data can be well described by an indefinite noncooperative or isodesmic association model. In this model, the deoxy Hb D associates noncooperatively to give a linear oligomeric chain with an equilibrium association constant of 1.42 × 10(4) M(-1) at 20°C for each step. The data are also well described by a monomer-dimer-tetramer equilibrium model with monomer-to-dimer and dimer-to-tetramer association constants of 1.87 and 1.03 × 10(4) M(-1) at 20°C, respectively. A hybrid recombinant Hb D was prepared with recombinant α(D)-globin and native β-globin to give a Hb D tetramer (α(2)(D)β(2)). This rHb D undergoes decreased deoxygenation-dependent self-association compared with the native Hb D. Residue glutamate 138 has previously been proposed to influence intertetramer interactions. Our results with recombinant Hb D show that Glu138 plays no role in deoxy Hb D intertetramer interactions.

Globin gene structure in a reptile supports the transpositional model for amniote α- and β-globin gene evolution

The haemoglobin protein, required for oxygen transportation in the body, is encoded by α- and β-globin genes that are arranged in clusters. The transpositional model for the evolution of distinct α-globin and β-globin clusters in amniotes is much simpler than the previously proposed whole genome duplication model. According to this model, all jawed vertebrates share one ancient region containing α- and β-globin genes and several flanking genes in the order MPG-C16orf35-(α-β)-GBY-LUC7L that has been conserved for more than 410 million years, whereas amniotes evolved a distinct β-globin cluster by insertion of a transposed β-globin gene from this ancient region into a cluster of olfactory receptors flanked by CCKBR and RRM1. It could not be determined whether this organisation is conserved in all amniotes because of the paucity of information from non-avian reptiles. To fill in this gap, we examined globin gene organisation in a squamate reptile, the Australian bearded dragon lizard, Pogona vitticeps (Agamidae). We report here that the α-globin cluster (HBK, HBA) is flanked by C16orf35 and GBY and is located on a pair of microchromosomes, whereas the β-globin cluster is flanked by RRM1 on the 3' end and is located on the long arm of chromosome 3. However, the CCKBR gene that flanks the β-globin cluster on the 5' end in other amniotes is located on the short arm of chromosome 5 in P. vitticeps, indicating that a chromosomal break between the β-globin cluster and CCKBR occurred at least in the agamid lineage. Our data from a reptile species provide further evidence to support the transpositional model for the evolution of β-globin gene cluster in amniotes.

Genomic organization and evolution of the Atlantic salmon hemoglobin repertoire

BACKGROUND

The genomes of salmonids are considered pseudo-tetraploid undergoing reversion to a stable diploid state. Given the genome duplication and extensive biological data available for salmonids, they are excellent model organisms for studying comparative genomics, evolutionary processes, fates of duplicated genes and the genetic and physiological processes associated with complex behavioral phenotypes. The evolution of the tetrapod hemoglobin genes is well studied; however, little is known about the genomic organization and evolution of teleost hemoglobin genes, particularly those of salmonids. The Atlantic salmon serves as a representative salmonid species for genomics studies. Given the well documented role of hemoglobin in adaptation to varied environmental conditions as well as its use as a model protein for evolutionary analyses, an understanding of the genomic structure and organization of the Atlantic salmon α and β hemoglobin genes is of great interest.

RESULTS

We identified four bacterial artificial chromosomes (BACs) comprising two hemoglobin gene clusters spanning the entire α and β hemoglobin gene repertoire of the Atlantic salmon genome. Their chromosomal locations were established using fluorescence in situ hybridization (FISH) analysis and linkage mapping, demonstrating that the two clusters are located on separate chromosomes. The BACs were sequenced and assembled into scaffolds, which were annotated for putatively functional and pseudogenized hemoglobin-like genes. This revealed that the tail-to-tail organization and alternating pattern of the α and β hemoglobin genes are well conserved in both clusters, as well as that the Atlantic salmon genome houses substantially more hemoglobin genes, including non-Bohr β globin genes, than the genomes of other teleosts that have been sequenced.

CONCLUSIONS

We suggest that the most parsimonious evolutionary path leading to the present organization of the Atlantic salmon hemoglobin genes involves the loss of a single hemoglobin gene cluster after the whole genome duplication (WGD) at the base of the teleost radiation but prior to the salmonid-specific WGD, which then produced the duplicated copies seen today. We also propose that the relatively high number of hemoglobin genes as well as the presence of non-Bohr β hemoglobin genes may be due to the dynamic life history of salmon and the diverse environmental conditions that the species encounters.Data deposition: BACs S0155C07 and S0079J05 (fps135): GenBank GQ898924; BACs S0055H05 and S0014B03 (fps1046): GenBank GQ898925.

Lineage-specific patterns of functional diversification in the alpha- and beta-globin gene families of tetrapod vertebrates

The alpha- and beta-globin gene families of jawed vertebrates have diversified with respect to both gene function and the developmental timing of gene expression. Phylogenetic reconstructions of globin gene family evolution have provided suggestive evidence that the developmental regulation of hemoglobin synthesis has evolved independently in multiple vertebrate lineages. For example, the embryonic beta-like globin genes of birds and placental mammals are not 1:1 orthologs. Despite the similarity in developmental expression profiles, the genes are independently derived from lineage-specific duplications of a beta-globin pro-ortholog. This suggests the possibility that other vertebrate taxa may also possess distinct repertoires of globin genes that were produced by repeated rounds of lineage-specific gene duplication and divergence. Until recently, investigations into this possibility have been hindered by the dearth of genomic sequence data from nonmammalian vertebrates. Here, we report new insights into globin gene family evolution that were provided by a phylogenetic analysis of vertebrate globins combined with a comparative genomic analysis of three key sauropsid taxa: a squamate reptile (anole lizard, Anolis carolinensis), a passeriform bird (zebra finch, Taeniopygia guttata), and a galliform bird (chicken, Gallus gallus). The main objectives of this study were 1) to characterize evolutionary changes in the size and membership composition of the alpha- and beta-globin gene families of tetrapod vertebrates and 2) to test whether functional diversification of the globin gene clusters occurred independently in different tetrapod lineages. Results of our comparative genomic analysis revealed several intriguing patterns of gene turnover in the globin gene clusters of different taxa. Lineage-specific differences in gene content were especially pronounced in the beta-globin gene family, as phylogenetic reconstructions revealed that amphibians, lepidosaurs (as represented by anole lizard), archosaurs (as represented by zebra finch and chicken), and mammals each possess a distinct independently derived repertoire of beta-like globin genes. In contrast to the ancient functional diversification of the alpha-globin gene cluster in the stem lineage of tetrapods, the physiological division of labor between early- and late-expressed genes in the beta-globin gene cluster appears to have evolved independently in several tetrapod lineages.

Differential loss of embryonic globin genes during the radiation of placental mammals

The differential gain and loss of genes from homologous gene families represents an important source of functional variation among the genomes of different species. Differences in gene content between species are primarily attributable to lineage-specific gene gains via duplication and lineage-specific losses via deletion or inactivation. Here, we use a comparative genomic approach to investigate this process of gene turnover in the beta-globin gene family of placental mammals. By analyzing genomic sequence data from representatives of each of the main superordinal clades of placental mammals, we were able to reconstruct pathways of gene family evolution during the basal radiation of this physiologically and morphologically diverse vertebrate group. Our analysis revealed that an initial expansion of the nonadult portion of the beta-globin gene cluster in the ancestor of placental mammals was followed by the differential loss and retention of ancestral gene lineages, thereby generating variation in the complement of embryonic globin genes among contemporary species. The sorting of epsilon-, gamma-, and eta-globin gene lineages among the basal clades of placental mammals has produced species differences in the functional types of hemoglobin isoforms that can be synthesized during the course of embryonic development.

Platypus globin genes and flanking loci suggest a new insertional model for beta-globin evolution in birds and mammals

BACKGROUND

Vertebrate alpha (alpha)- and beta (beta)-globin gene families exemplify the way in which genomes evolve to produce functional complexity. From tandem duplication of a single globin locus, the alpha- and beta-globin clusters expanded, and then were separated onto different chromosomes. The previous finding of a fossil beta-globin gene (omega) in the marsupial alpha-cluster, however, suggested that duplication of the alpha-beta cluster onto two chromosomes, followed by lineage-specific gene loss and duplication, produced paralogous alpha- and beta-globin clusters in birds and mammals. Here we analyse genomic data from an egg-laying monotreme mammal, the platypus (Ornithorhynchus anatinus), to explore haemoglobin evolution at the stem of the mammalian radiation.

RESULTS

The platypus alpha-globin cluster (chromosome 21) contains embryonic and adult alpha- globin genes, a beta-like omega-globin gene, and the GBY globin gene with homology to cytoglobin, arranged as 5'-zeta-zeta'-alphaD-alpha3-alpha2-alpha1-omega-GBY-3'. The platypus beta-globin cluster (chromosome 2) contains single embryonic and adult globin genes arranged as 5'-epsilon-beta-3'. Surprisingly, all of these globin genes were expressed in some adult tissues. Comparison of flanking sequences revealed that all jawed vertebrate alpha-globin clusters are flanked by MPG-C16orf35 and LUC7L, whereas all bird and mammal beta-globin clusters are embedded in olfactory genes. Thus, the mammalian alpha- and beta-globin clusters are orthologous to the bird alpha- and beta-globin clusters respectively.

CONCLUSION

We propose that alpha- and beta-globin clusters evolved from an ancient MPG-C16orf35-alpha-beta-GBY-LUC7L arrangement 410 million years ago. A copy of the original beta (represented by omega in marsupials and monotremes) was inserted into an array of olfactory genes before the amniote radiation (>315 million years ago), then duplicated and diverged to form orthologous clusters of beta-globin genes with different expression profiles in different lineages.

Adaptive functional divergence among triplicated alpha-globin genes in rodents

The functional divergence of duplicated genes is thought to play an important role in the evolution of new developmental and physiological pathways, but the role of positive selection in driving this process remains controversial. The objective of this study was to test whether amino acid differences among triplicated alpha-globin paralogs of the Norway rat (Rattus norvegicus) and the deer mouse (Peromyscus maniculatus) are attributable to a relaxation of purifying selection or to a history of positive selection that has adapted the gene products to new or modified physiological tasks. In each rodent species, the two paralogs at the 5'-end of the alpha-globin gene cluster (HBA-T1 and HBA-T2) are evolving in concert and are therefore identical or nearly identical in sequence. However, in each case, the HBA-T1 and HBA-T2 paralogs are distinguished from the third paralog at the 3'-end of the gene cluster (HBA-T3) by multiple amino acid substitutions. An analysis of genomic sequence data from several rodent species revealed that the HBA-T3 genes of Rattus and Peromyscus originated via independent, lineage-specific duplication events. In the independently derived HBA-T3 genes of both species, a likelihood analysis based on a codon-substitution model revealed that accelerated rates of amino acid substitution are attributable to positive directional selection, not to a relaxation of purifying selection. As a result of functional divergence among the triplicated alpha-globin genes in Rattus and Peromyscus, the red blood cells of both rodent species contain a mixture of functionally distinct alpha-chain hemoglobin isoforms that are predicted to have different oxygen-binding affinities. In P. maniculatus, a species that is able to sustain physiological function under conditions of chronic hypoxia at high altitude, the coexpression of distinct hemoglobin isoforms with graded oxygen affinities is expected to broaden the permissible range of arterial oxygen tensions for pulmonary/tissue oxygen transport.

Rapid rates of lineage-specific gene duplication and deletion in the alpha-globin gene family

Phylogeny reconstructions of the globin gene families have revealed that paralogous genes within species are often more similar to one another than they are to their orthologous counterparts in closely related species. This pattern has been previously attributed to mechanisms of concerted evolution such as interparalog gene conversion that homogenize sequence variation between tandemly duplicated genes and therefore create the appearance of recent common ancestry. Here we report a comparative genomic analysis of the alpha-globin gene family in mammals that reveal a surprisingly high rate of lineage-specific gene duplication and deletion via unequal crossing-over. Results of our analysis reveal that patterns of sequence similarity between paralogous alpha-like globin genes from the same species are only partly explained by concerted evolution between preexisting gene duplicates. In a number of cases, sequence similarity between paralogous sequences from the same species is attributable to recent ancestry between the products of de novo gene duplications. As a result of this surprisingly rapid rate of gene gain and loss, many mammals possess alpha-like globin genes that have no orthologous counterparts in closely related species. The resultant variation in gene copy number among species may represent an important source of regulatory variation that affects physiologically important aspects of blood oxygen transport and aerobic energy metabolism.