New genes originated via multiple recombinational pathways in the beta-globin gene family of rodents

A New Homotetramer Hemoglobin in the Pulmonary Surfactant of Plateau Zokors (Myospalax Baileyi)

The plateau zokor (Myospalax baileyi) is a native species to the Qinghai-Tibetan Plateau, inhabiting hypoxia and hypercapnia sealed subterranean burrows that pose several unique physiological challenges. In this study, we observed a novel heme-containing protein in the pulmonary surfactant (PS) of plateau zokor, identified the encoding gene of the protein, predicted its origination and structure, verified its expression in alveolar epithelial cells, and determined the protein’s affinity to oxygen and its effect on the oxygen-dissolving capability in the PS of plateau zokors. The protein is an unusual homotetramer hemoglobin consisting of four γ-like subunits, and the subunit is encoded by a paralog gene of γ, that is γ-like. The divergence time of γ-like from γ is estimated by the molecular clock to be about 2.45 Mya. The generation of γ-like in plateau zokors might well relate to long-time stress of the high land hypoxia. Unlike γ, the γ-like has a hypoxia response element (HRE) and a lung tissue-specific enhancer in its upstream region, and it is expressed specifically in lung tissues and up-regulated by hypoxia. The protein is named as γ₄-like which is expressed specifically in Alveolar epithelial type II (ATII) cells and secreted into the alveolar cavities through the osmiophilic multilamellar body (LBs). The γ₄-like has a higher affinity to oxygen, and that increases significantly oxygen-dissolving capability in the PS of plateau zokors by its oxygenation function, which might be beneficial for the plateau zokors to obtain oxygen from the severe hypoxia environments by facilitating oxygen diffusion from alveoli to blood.

Genetic variation in haemoglobin is associated with evolved changes in breathing in high-altitude deer mice

Physiological systems often have emergent properties but the effects of genetic variation on physiology are often unknown, which presents a major challenge to understanding the mechanisms of phenotypic evolution. We investigated whether genetic variants in haemoglobin (Hb) that contribute to high-altitude adaptation in deer mice (Peromyscus maniculatus) are associated with evolved changes in the control of breathing. We created F2 inter-population hybrids of highland and lowland deer mice to test for phenotypic associations of α- and β-globin variants on a mixed genetic background. Hb genotype had expected effects on Hb-O2 affinity that were associated with differences in arterial O2 saturation in hypoxia. However, high-altitude genotypes were also associated with breathing phenotypes that should contribute to enhancing O2 uptake in hypoxia. Mice with highland α-globin exhibited a more effective breathing pattern, with highland homozygotes breathing deeper but less frequently across a range of inspired O2, and this difference was comparable to the evolved changes in breathing pattern in deer mouse populations native to high altitude. The ventilatory response to hypoxia was augmented in mice that were homozygous for highland β-globin. The association of globin variants with variation in breathing phenotypes could not be recapitulated by acute manipulation of Hb-O2 affinity, because treatment with efaproxiral (a synthetic drug that acutely reduces Hb-O2 affinity) had no effect on breathing in normoxia or hypoxia. Therefore, adaptive variation in Hb may have unexpected effects on physiology in addition to the canonical function of this protein in circulatory O2 transport.

The adaptive benefit of evolved increases in hemoglobin-O2 affinity is contingent on tissue O2 diffusing capacity in high-altitude deer mice

BACKGROUND

Complex organismal traits are often the result of multiple interacting genes and sub-organismal phenotypes, but how these interactions shape the evolutionary trajectories of adaptive traits is poorly understood. We examined how functional interactions between cardiorespiratory traits contribute to adaptive increases in the capacity for aerobic thermogenesis (maximal O₂ consumption, V̇O₂max, during acute cold exposure) in high-altitude deer mice (Peromyscus maniculatus). We crossed highland and lowland deer mice to produce F₂ inter-population hybrids, which expressed genetically based variation in hemoglobin (Hb) O₂ affinity on a mixed genetic background. We then combined physiological experiments and mathematical modeling of the O₂ transport pathway to examine the links between cardiorespiratory traits and V̇O₂max.

RESULTS

Physiological experiments revealed that increases in Hb-O₂ affinity of red blood cells improved blood oxygenation in hypoxia but were not associated with an enhancement in V̇O₂max. Sensitivity analyses performed using mathematical modeling showed that the influence of Hb-O₂ affinity on V̇O₂max in hypoxia was contingent on the capacity for O₂ diffusion in active tissues.

CONCLUSIONS

These results suggest that increases in Hb-O₂ affinity would only have adaptive value in hypoxic conditions if concurrent with or preceded by increases in tissue O₂ diffusing capacity. In high-altitude deer mice, the adaptive benefit of increasing Hb-O₂ affinity is contingent on the capacity to extract O₂ from the blood, which helps resolve controversies about the general role of hemoglobin function in hypoxia tolerance.

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.

Evolution of Hemoglobin Genes in a Subterranean Rodent Species (Lasiopodomys mandarinus)

The Mandarin vole (Lasiopodomys mandarinus), a typical subterranean rodent, has undergone hematological adaptations to tolerate the hypoxic/hypercapnic underground environment. Hemoglobin (Hb) genes encode respiratory proteins functioning principally in oxygen binding and transport to various tissues and organs. To investigate the evolution of α- and β-hemoglobin (Hb) in subterranean rodent species, we sequenced Hb genes of the Mandarin vole and the related aboveground Brandt's vole (L. brandtii). Sequencing showed that in both voles, α-globin was encoded by a cluster of five functional genes in the following linkage order: HBZ, HBA-T1, HBQ-T1, HBA-T2, and HBQ-T2; among these, HBQ-T2 is a pseudogene in both voles. The β-globin gene cluster in both voles also included five functional genes in the following linkage order: HBE, HBE/HBG, HBG, HBB-T1, and HBB-T2. Phylogenetic analysis revealed that the Mandarin vole underwent convergent evolution with its related aboveground species (Brandt's vole) but not with other subterranean rodent species. Selection pressure analyses revealed that α- and β-globin genes are under strong purifying selection (ω < 1), and branch-site analyses identified positive selection sites on HBAT-T1 and HBB-T1 in different subterranean rodent species. This suggests that the adaptive evolution of these genes enhanced the ability of Hb to store and transport oxygen in subterranean rodent species. Our findings highlight the critical roles of Hb genes in the evolution of hypoxia tolerance in subterranean rodent species.

Ontogenesis of evolved changes in respiratory physiology in deer mice native to high altitude

High-altitude environments are cold and hypoxic, and many high-altitude natives have evolved changes in respiratory physiology that improve O2 uptake in hypoxia as adults. Altricial mammals undergo a dramatic metabolic transition from ectothermy to endothermy in early post-natal life, which may influence the ontogenetic development of respiratory traits at high altitude. We examined the developmental changes in respiratory and haematological traits in deer mice (Peromyscus maniculatus) native to high altitude, comparing the respiratory responses to progressive hypoxia between highland and lowland deer mice. Among adults, highlanders exhibited higher total ventilation and a more effective breathing pattern (relatively deeper tidal volumes), for mice that were caught and tested at their native altitudes and those lab-raised in normoxia. Lab-raised progeny of each population were also tested at post-natal day (P)7, 14, 21 and 30. Highlanders developed an enhanced hypoxic ventilatory response by P21, concurrent with the full maturation of the carotid bodies, and their more effective breathing pattern arose by P14; these ages correspond to critical benchmarks in the full development of homeothermy in highlanders. However, highlanders exhibited developmental delays in ventilatory sensitivity to hypoxia, hyperplasia of type I cells in the carotid body and increases in blood haemoglobin content compared with lowland mice. Nevertheless, highlanders maintained consistently higher arterial O2 saturation in hypoxia across development, in association with increases in blood-O2 affinity that were apparent from birth. We conclude that evolved changes in respiratory physiology in high-altitude deer mice become expressed in association with the post-natal development of endothermy.

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.

Adaptive Changes in Hemoglobin Function in High-Altitude Tibetan Canids Were Derived via Gene Conversion and Introgression

A key question in evolutionary biology concerns the relative importance of different sources of adaptive genetic variation, such as de novo mutations, standing variation, and introgressive hybridization. A corollary question concerns how allelic variants derived from these different sources may influence the molecular basis of phenotypic adaptation. Here, we use a protein-engineering approach to examine the phenotypic effect of putatively adaptive hemoglobin (Hb) mutations in the high-altitude Tibetan wolf that were selectively introgressed into the Tibetan mastiff, a high-altitude dog breed that is renowned for its hypoxia tolerance. Experiments revealed that the introgressed coding variants confer an increased Hb–O₂ affinity in conjunction with an enhanced Bohr effect. We also document that affinity-enhancing mutations in the β-globin gene of Tibetan wolf were originally derived via interparalog gene conversion from a tandemly linked β-globin pseudogene. Thus, affinity-enhancing mutations were introduced into the β-globin gene of Tibetan wolf via one form of intragenomic lateral transfer (ectopic gene conversion) and were subsequently introduced into the Tibetan mastiff genome via a second form of lateral transfer (introgression). Site-directed mutagenesis experiments revealed that the increased Hb–O₂ affinity requires a specific two-site combination of amino acid replacements, suggesting that the molecular underpinnings of Hb adaptation in Tibetan mastiff (involving mutations that arose in a nonexpressed gene and which originally fixed in Tibetan wolf) may be qualitatively distinct from functionally similar changes in protein function that could have evolved via sequential fixation of de novo mutations during the breed’s relatively short duration of residency at high altitude.

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.

PKR activation and eIF2α phosphorylation mediate human globin mRNA splicing at spliceosome assembly

Short elements in mammalian mRNA can control gene expression by activating the RNA-dependent protein kinase PKR that attenuates translation by phosphorylating cytoplasmic eukaryotic initiation factor 2α (eIF2α). We demonstrate a novel, positive role for PKR activation and eIF2α phosphorylation in human globin mRNA splicing. PKR localizes in splicing complexes and associates with splicing factor SC35. Splicing and early-stage spliceosome assembly on β-globin pre-mRNA depend strictly on activation of PKR by a codon-containing RNA fragment within exon 1 and on phosphorylation of nuclear eIF2α on Serine 51. Nonphosphorylatable mutant eIF2αS51A blocked β-globin mRNA splicing in cells and nuclear extract. Mutations of the β-globin RNA activator abrogated PKR activation and profoundly affected mRNA splicing efficiency. PKR depletion abrogated splicing and spliceosome assembly; recombinant PKR effectively restored splicing. Excision of the first intron of β-globin induces strand displacement within the RNA activator of PKR by a sequence from exon 2, a structural rearrangement that silences the ability of spliced β-globin mRNA to activate PKR. Thus, the ability to activate PKR is transient, serving solely to enable splicing. α-Globin pre-mRNA splicing is controlled likewise but positions of PKR activator and silencer are reversed, demonstrating evolutionary flexibility in how PKR activation regulates globin mRNA splicing through eIF2α phosphorylation.

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.

Genetically based low oxygen affinities of felid hemoglobins: lack of biochemical adaptation to high-altitude hypoxia in the snow leopard

Genetically based modifications of hemoglobin (Hb) function that increase blood-O2 affinity are hallmarks of hypoxia adaptation in vertebrates. Among mammals, felid Hbs are unusual in that they have low intrinsic O2 affinities and reduced sensitivities to the allosteric cofactor 2,3-diphosphoglycerate (DPG). This combination of features compromises the acclimatization capacity of blood-O2 affinity and has led to the hypothesis that felids have a restricted physiological niche breadth relative to other mammals. In seeming defiance of this conjecture, the snow leopard (Panthera uncia) has an extraordinarily broad elevational distribution and occurs at elevations above 6000 m in the Himalayas. Here, we characterized structural and functional variation of big cat Hbs and investigated molecular mechanisms of Hb adaptation and allosteric regulation that may contribute to the extreme hypoxia tolerance of the snow leopard. Experiments revealed that purified Hbs from snow leopard and African lion exhibited equally low O2 affinities and DPG sensitivities. Both properties are primarily attributable to a single amino acid substitution, β2His→Phe, which occurred in the common ancestor of Felidae. Given the low O2 affinity and reduced regulatory capacity of feline Hbs, the extreme hypoxia tolerance of snow leopards must be attributable to compensatory modifications of other steps in the O2-transport pathway.

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.

Intraspecific polymorphism, interspecific divergence, and the origins of function-altering mutations in deer mouse hemoglobin

Major challenges for illuminating the genetic basis of phenotypic evolution are to identify causative mutations, to quantify their functional effects, to trace their origins as new or preexisting variants, and to assess the manner in which segregating variation is transduced into species differences. Here, we report an experimental analysis of genetic variation in hemoglobin (Hb) function within and among species of Peromyscus mice that are native to different elevations. A multilocus survey of sequence variation in the duplicated HBA and HBB genes in Peromyscus maniculatus revealed that function-altering amino acid variants are widely shared among geographically disparate populations from different elevations, and numerous amino acid polymorphisms are also shared with closely related species. Variation in Hb-O₂ affinity within and among populations of P. maniculatus is attributable to numerous amino acid mutations that have individually small effects. One especially surprising feature of the Hb polymorphism in P. maniculatus is that an appreciable fraction of functional standing variation in the two transcriptionally active HBA paralogs is attributable to recurrent gene conversion from a tandemly linked HBA pseudogene. Moreover, transpecific polymorphism in the duplicated HBA genes is not solely attributable to incomplete lineage sorting or introgressive hybridization; instead, it is mainly attributable to recurrent interparalog gene conversion that has occurred independently in different species. Partly as a result of concerted evolution between tandemly duplicated globin genes, the same amino acid changes that contribute to variation in Hb function within P. maniculatus also contribute to divergence in Hb function among different species of Peromyscus. In the case of function-altering Hb mutations in Peromyscus, there is no qualitative or quantitative distinction between segregating variants within species and fixed differences between species.

Hemoglobin function and allosteric regulation in semi-fossorial rodents (family Sciuridae) with different altitudinal ranges

Semi-fossorial ground squirrels face challenges to respiratory gas transport associated with the chronic hypoxia and hypercapnia of underground burrows, and such challenges are compounded in species that are native to high altitude. During hibernation, such species must also contend with vicissitudes of blood gas concentrations and plasma pH caused by episodic breathing. Here, we report an analysis of hemoglobin (Hb) function in six species of marmotine ground squirrels with different altitudinal distributions. Regardless of their native altitude, all species have high Hb-O2 affinities, mainly due to suppressed sensitivities to allosteric effectors [2,3-diphosphoglycerate (DPG) and chloride ions]. This suppressed anion sensitivity is surprising given that all canonical anion-binding sites are conserved. Two sciurid species, the golden-mantled and thirteen-lined ground squirrel, have Hb-O2 affinities that are characterized by high pH sensitivity and low thermal sensitivity relative to the Hbs of humans and other mammals. The pronounced Bohr effect is surprising in light of highly unusual amino acid substitutions at the C-termini that are known to abolish the Bohr effect in human HbA. Taken together, the high O2 affinity of sciurid Hbs suggests an enhanced capacity for pulmonary O2 loading under hypoxic and hypercapnic conditions, while the large Bohr effect should help to ensure efficient O2 unloading in tissue capillaries. In spite of the relatively low thermal sensitivities of the sciurid Hbs, our results indicate that the effect of hypothermia on Hb oxygenation is the main factor contributing to the increased blood-O2 affinity in hibernating ground squirrels.

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.

Whole-genome duplication and the functional diversification of teleost fish hemoglobins

Subsequent to the two rounds of whole-genome duplication that occurred in the common ancestor of vertebrates, a third genome duplication occurred in the stem lineage of teleost fishes. This teleost-specific genome duplication (TGD) is thought to have provided genetic raw materials for the physiological, morphological, and behavioral diversification of this highly speciose group. The extreme physiological versatility of teleost fish is manifest in their diversity of blood–gas transport traits, which reflects the myriad solutions that have evolved to maintain tissue O₂ delivery in the face of changing metabolic demands and environmental O₂ availability during different ontogenetic stages. During the course of development, regulatory changes in blood–O₂ transport are mediated by the expression of multiple, functionally distinct hemoglobin (Hb) isoforms that meet the particular O₂-transport challenges encountered by the developing embryo or fetus (in viviparous or oviparous species) and in free-swimming larvae and adults. The main objective of the present study was to assess the relative contributions of whole-genome duplication, large-scale segmental duplication, and small-scale gene duplication in producing the extraordinary functional diversity of teleost Hbs. To accomplish this, we integrated phylogenetic reconstructions with analyses of conserved synteny to characterize the genomic organization and evolutionary history of the globin gene clusters of teleosts. These results were then integrated with available experimental data on functional properties and developmental patterns of stage-specific gene expression. Our results indicate that multiple α- and β-globin genes were present in the common ancestor of gars (order Lepisoteiformes) and teleosts. The comparative genomic analysis revealed that teleosts possess a dual set of TGD-derived globin gene clusters, each of which has undergone lineage-specific changes in gene content via repeated duplication and deletion events. Phylogenetic reconstructions revealed that paralogous genes convergently evolved similar functional properties in different teleost lineages. Consistent with other recent studies of globin gene family evolution in vertebrates, our results revealed evidence for repeated evolutionary transitions in the developmental regulation of Hb synthesis.

Gene duplication, genome duplication, and the functional diversification of vertebrate globins

The functional diversification of the vertebrate globin gene superfamily provides an especially vivid illustration of the role of gene duplication and whole-genome duplication in promoting evolutionary innovation. For example, key globin proteins that evolved specialized functions in various aspects of oxidative metabolism and oxygen signaling pathways (hemoglobin [Hb], myoglobin [Mb], and cytoglobin [Cygb]) trace their origins to two whole-genome duplication events in the stem lineage of vertebrates. The retention of the proto-Hb and Mb genes in the ancestor of jawed vertebrates permitted a physiological division of labor between the oxygen-carrier function of Hb and the oxygen-storage function of Mb. In the Hb gene lineage, a subsequent tandem gene duplication gave rise to the proto α- and β-globin genes, which permitted the formation of multimeric Hbs composed of unlike subunits (α(2)β(2)). The evolution of this heteromeric quaternary structure was central to the emergence of Hb as a specialized oxygen-transport protein because it provided a mechanism for cooperative oxygen-binding and allosteric regulatory control. Subsequent rounds of duplication and divergence have produced diverse repertoires of α- and β-like globin genes that are ontogenetically regulated such that functionally distinct Hb isoforms are expressed during different stages of prenatal development and postnatal life. In the ancestor of jawless fishes, the proto Mb and Hb genes appear to have been secondarily lost, and the Cygb homolog evolved a specialized respiratory function in blood-oxygen transport. Phylogenetic and comparative genomic analyses of the vertebrate globin gene superfamily have revealed numerous instances in which paralogous globins have convergently evolved similar expression patterns and/or similar functional specializations in different organismal lineages.

Copy number polymorphism in the α-globin gene cluster of European rabbit (Oryctolagus cuniculus)

Comparative genomic studies have revealed that mammals typically possess two or more tandemly duplicated copies of the α-globin (HBA) gene. The domestic rabbit represents an exception to this general rule, as this species was found to possess a single HBA gene. Previous electrophoretic surveys of HBA polymorphism in natural populations of the European rabbit (Oryctolagus cuniculus) revealed extensive geographic variation in the frequencies of three main electromorphs. The variation in frequency of two electromorphs is mainly partitioned between two distinct subspecies of European rabbit, and a third is restricted to the hybrid zone between the two rabbit subspecies in Iberia. Here we report the results of a survey of nucleotide polymorphism, which revealed HBA copy number polymorphism in Iberian populations of the European rabbit. By characterizing patterns of HBA polymorphism in populations from the native range of the European rabbit, we were able to identify the specific amino-acid substitutions that distinguish the previously characterized electromorphs. Within the hybrid zone, we observed the existence of a second HBA gene duplicate, named HBA2, that mostly represents a novel sequence haplotype, which occurs in higher frequency within the hybrid zone, and thus appears to have arisen in hybrids of the two distinct subspecies. Although this novel gene is also present in other wild Iberian populations, it is almost absent from French populations, which suggest a recent ancestry, associated with the establishment of the post-Pleistocene contact zone between the two European rabbit subspecies.

Oxygenation properties and oxidation rates of mouse hemoglobins that differ in reactive cysteine content

House mice (genus Mus) harbor extensive allelic variation at two tandemly duplicated genes that encode the β-chain subunits of adult hemoglobin (Hb). Alternative haplotypes differ in the level of sequence divergence between the two β-globin gene duplicates: the Hbb(d) and Hbb(p) haplotypes harbor two structurally distinct β-globin genes, whereas the Hbb(s) haplotype harbors two β-globin duplicates that are identical in sequence. One especially salient difference between the s-type Hbs relative to the d- and p-type Hbs relates to the number of reactive β-chain cysteine residues. In addition to the highly conserved cysteine residue at β93, the d- and p-type Hbs contain an additional reactive cysteine residue at β13. To assess the functional consequences of allelic variation in β-globin cysteine content, we measured O(2)-binding properties and H(2)O(2)-induced oxidation rates of mono- and dicysteinyl β-Hbs from 4 different inbred strains of mice: C57BL/6J, BALB/cJ, MSM/Ms, and CAROLI/EiJ. The experiments revealed that purified Hbs from the various mouse strains did not exhibit substantial variation in O(2)-binding properties, but s-type Hb (which contains a single reactive β-chain cysteine residue) was far more readily oxidized to Fe(3+) metHb by H(2)O(2) than other mouse Hbs that contain two reactive β-chain cysteine residues. These results suggest that the possession of an additional reactive cysteine residue may protect against metHb formation under oxidizing conditions. The allelic differences in β-globin cysteine content could affect aspects of redox signaling and oxidative/nitrosative stress responses that are mediated by Hb-S-nitrosylation and Hb-S-glutathionylation pathways.

Altitudinal variation at duplicated β-globin genes in deer mice: effects of selection, recombination, and gene conversion

Spatially varying selection on a given polymorphism is expected to produce a localized peak in the between-population component of nucleotide diversity, and theory suggests that the chromosomal extent of elevated differentiation may be enhanced in cases where tandemly linked genes contribute to fitness variation. An intriguing example is provided by the tandemly duplicated β-globin genes of deer mice (Peromyscus maniculatus), which contribute to adaptive differentiation in blood-oxygen affinity between high- and low-altitude populations. Remarkably, the two β-globin genes segregate the same pair of functionally distinct alleles due to a history of interparalog gene conversion and alleles of the same functional type are in perfect coupling-phase linkage disequilibrium (LD). Here we report a multilocus analysis of nucleotide polymorphism and LD in highland and lowland mice with different genetic backgrounds at the β-globin genes. The analysis of haplotype structure revealed a paradoxical pattern whereby perfect LD between the two β-globin paralogs (which are separated by 16.2 kb) is maintained in spite of the fact that LD within both paralogs decays to background levels over physical distances of less than 1 kb. The survey of nucleotide polymorphism revealed that elevated levels of altitudinal differentiation at each of the β-globin genes drop away quite rapidly in the external flanking regions (upstream of the 5' paralog and downstream of the 3' paralog), but the level of differentiation remains unexpectedly high across the intergenic region. Observed patterns of diversity and haplotype structure are difficult to reconcile with expectations of a two-locus selection model with multiplicative fitness.

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

If rates of postduplication gene retention are positively correlated with levels of functional constraint, then gene duplicates that have been retained in a restricted number of taxonomic lineages would be expected to exhibit relatively low levels of sequence conservation. Paradoxical patterns are presented by gene duplicates that have been retained in a small number of taxa but which are nonetheless subject to strong purifying selection relative to paralogous members of the same multigene family. This pattern suggests that such genes may have been co-opted for novel, lineage-specific functions. One possible example involves the enigmatic globin-E gene (GbE), which appears to be exclusively restricted to birds. Available data indicate that this gene is expressed exclusively in the avian eye, but its physiological function remains a mystery. In contrast to the highly restricted phyletic distribution of GbE, the overwhelming majority of jawed vertebrates (gnathostomes) possess copies of the related cytoglobin (Cygb) and myoglobin (Mb) genes. The purpose of the present study was 1) to assess the phyletic distribution of the Cygb, Mb, and GbE genes among vertebrates, 2) to elucidate the duplicative origins and evolutionary histories of these three genes, and 3) to evaluate the relative levels of functional constraint of these genes based on comparative sequence analysis. To accomplish these objectives, we conducted a combined phylogenetic and comparative genomic analysis involving taxa that represent each of the major lineages of gnathostome vertebrates. Results of synteny comparisons and phylogenetic topology tests revealed that GbE is clearly not the product of a recent, bird-specific duplication event. Instead, GbE originated via duplication of a proto-Mb gene in the stem lineage of gnathostomes. Unlike the Mb gene, which has been retained in all major gnathostome lineages other than amphibians, the GbE gene has been retained only in the lineage leading to modern birds and has been independently lost in at least four major lineages: teleost fish, amphibians, mammals, and nonavian reptiles. Despite the restricted phyletic distribution of this gene, our results indicate that GbE is one of the most highly conserved globins in the avian genome.

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.

A tale of two genomes: contrasting patterns of phylogeographic structure in a widely distributed bat

One of the most widely distributed bats in the New World, the big brown bat (Eptesicus fuscus) exhibits well-documented geographic variation in morphology and life history traits, suggesting the potential for significant phylogeographic structure as well as adaptive differentiation among populations. In a pattern broadly consistent with morphologically defined subspecies, we found deeply divergent mitochondrial lineages restricted to different geographic regions. In contrast, sequence data from two nuclear loci suggest a general lack of regional genetic structure except for peripheral populations in the Caribbean and Mexico/South America. Coalescent analyses suggest that the striking difference in population structure between genomes cannot be attributed solely to different rates of lineage sorting, but is likely due to male-mediated gene flow homogenizing nuclear genetic diversity across most of the continental range. Despite this ongoing gene flow, selection has apparently been effective in producing and maintaining adaptive differentiation among populations, while strong female site fidelity, maintained over the course of millions of years, has produced remarkably deep divergence among geographically isolated matrilines. Our results highlight the importance of evaluating multiple genetic markers for a more complete understanding of population structure and history.

Genetic differences in hemoglobin function between highland and lowland deer mice

In high-altitude vertebrates, adaptive changes in blood-O(2) affinity may be mediated by modifications of hemoglobin (Hb) structure that affect intrinsic O(2) affinity and/or responsiveness to allosteric effectors that modulate Hb-O(2) affinity. This mode of genotypic specialization is considered typical of mammalian species that are high-altitude natives. Here we investigated genetically based differences in Hb-O(2) affinity between highland and lowland populations of the deer mouse (Peromyscus maniculatus), a generalist species that has the broadest altitudinal distribution of any North American mammal. The results of a combined genetic and proteomic analysis revealed that deer mice harbor a high level of Hb isoform diversity that is attributable to allelic polymorphism at two tandemly duplicated alpha-globin genes and two tandemly duplicated beta-globin genes. This high level of isoHb diversity translates into a correspondingly high level of interindividual variation in Hb functional properties. O(2) equilibrium experiments revealed that the Hbs of highland mice exhibit slightly higher intrinsic O(2) affinities and significantly lower Cl(-) sensitivities relative to the Hbs of lowland mice. The experiments also revealed distinct biochemical properties of deer mouse Hb related to the anion-dependent allosteric regulation of O(2) affinity. In conjunction with previous findings, our results demonstrate that modifications of Hb structure that alter allosteric anion sensitivity play an important role in the adaptive fine-tuning of blood-O(2) affinity.

Evolutionary and functional properties of a two-locus beta-globin polymorphism in Indian house mice

Electrophoretic surveys of hemoglobin (Hb) polymorphism in house mice from South Asia and the Middle East have revealed that two alternative beta-globin haplotypes, Hbb(d) and Hbb(p), are often present at intermediate frequencies in geographically disparate populations. Both haplotypes harbor two functionally distinct beta-globin paralogs, HBB-T1 (which encodes the beta-chain subunits of the major Hb isoform) and HBB-T2 (which encodes the beta-chains of the minor Hb isoform). The Hbb(d) and Hbb(p) haplotypes share identical HBB-T1 alleles, but products of the alternative HBB-T2 alleles (d(minor) and p(minor)) are distinguished by two amino acid substitutions. To investigate the possible adaptive significance of the Hbb(d)/Hbb(p) polymorphism we conducted a population genetic analysis of the duplicated beta-globin genes of Indian house mice (Mus castaneus) in conjunction with experimental studies of Hb function in inbred strains of mice that carry the alternative Hbb(d) and Hbb(p) haplotypes. The main objectives of this study were (i) to characterize patterns of nucleotide polymorphism and linkage disequilibrium in the duplicated beta-globin genes of M. castaneus, (ii) to test the hypothesis that the Hbb(d) and Hbb(p) haplotypes are maintained as a balanced polymorphism, and (iii) to assess whether allelic differences in the alternative minor Hb isoforms (d(minor) and p(minor)) are associated with different O(2)-binding properties. A multilocus analysis of polymorphism and divergence revealed that levels of diversity at the HBB-T2 gene exceeded neutral expectations, and reconstructed haplotype networks for both beta-globin paralogs revealed extensive allele sharing with several other closely related species of Mus. However, despite this suggestive evidence for balancing selection, O(2)-equilibrium curves revealed no discernible functional differences between red cell lysates containing the d(minor) and p(minor) Hb isoforms. If the d(minor) and p(minor) alleles are maintained as a balanced polymorphism, our results indicate that the associated fitness variance is not directly related to respiratory functions of Hb.

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.

Evolution of duplicated beta-globin genes and the structural basis of hemoglobin isoform differentiation in Mus

The functional diversification of multigene families may be strongly influenced by mechanisms of concerted evolution such as interparalog gene conversion. The beta-globin gene family of house mice (genus Mus) represents an especially promising system for evaluating the effects of gene conversion on the functional divergence of duplicated genes. Whereas the majority of mammalian species possess tandemly duplicated copies of the adult beta-globin gene that are identical in sequence, natural populations of house mice are often polymorphic for distinct two-locus haplotypes that differ in levels of functional divergence between duplicated beta-globin genes, HBB-T1 and HBB-T2. Here, we use a phylogenetic approach to unravel the complex evolutionary history of the HBB-T1 and HBB-T2 paralogs in a taxonomically diverse set of species in the genus Mus. The main objectives of this study were 1) to reconstruct the evolutionary history of the different HBB haplotypes of house mice, 2) to assess the role of recombinational exchange between HBB-T1 and HBB-T2 in promoting concerted evolution, 3) to assess the role of recombinational exchange between HBB-T1 and HBB-T2 in creating chimeric genes, and 4) to assess the structural basis of hemoglobin isoform differentiation in species that possess distinct HBB paralogs. Results of our phylogenetic survey revealed that the HBB-T1 and HBB-T2 genes in different species of Mus exhibit the full range of evolutionary outcomes with respect to levels of interparalog divergence. At one end of the spectrum, the two identical HBB paralogs on the Hbb(s) haplotype (shared by Mus domesticus, Mus musculus, and Mus spretus) represent a classic example of concerted evolution. At the other end of the spectrum, the two distinct HBB paralogs on the Hbb(d), Hbb(p), Hbb(w1), and Hbb(w2) haplotypes (shared by multiple species in the subgenus Mus) show no trace of gene conversion and are distinguished by a number of functionally important amino acid substitutions. Because the possession of distinct HBB paralogs expands the repertoire of functionally distinct hemoglobin isoforms that can be synthesized during fetal development and postnatal life, variation in the level of functional divergence between HBB-T1 and HBB-T2 may underlie important physiological variation within and among species.

Evolutionary and functional insights into the mechanism underlying high-altitude adaptation of deer mouse hemoglobin

Adaptive modifications of heteromeric proteins may involve genetically based changes in single subunit polypeptides or parallel changes in multiple genes that encode distinct, interacting subunits. Here we investigate these possibilities by conducting a combined evolutionary and functional analysis of duplicated globin genes in natural populations of deer mice (Peromyscus maniculatus) that are adapted to different elevational zones. A multilocus analysis of nucleotide polymorphism and linkage disequilibrium revealed that high-altitude adaptation of deer mouse hemoglobin involves parallel functional differentiation at multiple unlinked gene duplicates: two alpha-globin paralogs on chromosome 8 and two beta-globin paralogs on chromosome 1. Differences in O(2)-binding affinity of the alternative beta-chain hemoglobin isoforms were entirely attributable to allelic differences in sensitivity to 2,3-diphosphoglycerate (DPG), an allosteric cofactor that stabilizes the low-affinity, deoxygenated conformation of the hemoglobin tetramer. The two-locus beta-globin haplotype that predominates at high altitude is associated with suppressed DPG-sensitivity (and hence, increased hemoglobin-O(2) affinity), which enhances pulmonary O(2) loading under hypoxia. The discovery that allelic differences in DPG-sensitivity contribute to adaptive variation in hemoglobin-O(2) affinity illustrates the value of integrating evolutionary analyses of sequence variation with mechanistic appraisals of protein function. Investigation into the functional significance of the deer mouse beta-globin polymorphism was motivated by the results of population genetic analyses which revealed evidence for a history of divergent selection between elevational zones. The experimental measures of O(2)-binding properties corroborated the tests of selection by demonstrating a functional difference between the products of alternative alleles.

Origin and ascendancy of a chimeric fusion gene: the beta/delta-globin gene of paenungulate mammals

The delta-globin gene (HBD) of eutherian mammals exhibits a propensity for recombinational exchange with the closely linked beta-globin gene (HBB) and has been independently converted by the HBB gene in multiple lineages. Here we report the presence of a chimeric beta/delta fusion gene in the African elephant (Loxodonta africana) that was created by unequal crossing-over between misaligned HBD and HBB paralogs. The recombinant chromosome that harbors the beta/delta fusion gene in elephants is structurally similar to the "anti-Lepore" duplication mutant of humans (the reciprocal exchange product of the hemoglobin Lepore deletion mutant). However, the situation in the African elephant is unique in that the chimeric beta/delta fusion gene supplanted the parental HBB gene and is therefore solely responsible for synthesizing the beta-chain subunits of adult hemoglobin. A phylogenetic survey of beta-like globin genes in afrotherian and xenarthran mammals revealed that the origin of the chimeric beta/delta fusion gene and the concomitant inactivation of the HBB gene predated the radiation of "Paenungulata," a clade of afrotherian mammals that includes three orders: Proboscidea (elephants), Sirenia (dugongs and manatees), and Hyracoidea (hyraxes). The reduced fitness of the human Hb Lepore deletion mutant helps to explain why independently derived beta/delta fusion genes (which occur on an anti-Lepore chromosome) have been fixed in a number of mammalian lineages, whereas the reciprocal delta/beta fusion gene (which occurs on a Lepore chromosome) has yet to be documented in any nonhuman mammal. This illustrates how the evolutionary fates of chimeric fusion genes can be strongly influenced by their recombinational mode of origin.