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-O2 affinity in multiple high-altitude taxa. In six of the eight taxon pairs, high-altitude taxa evolved derived increases in Hb-O2 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-O2 affinity, indicating that the repeatability of phenotypic evolution does not require parallelism at the molecular level.

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

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

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.

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.

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.

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.

Decoding the architecture and origins of mechanisms for developmental polyphenism.

Developmental polyphenism affords a single genotype multiple solutions to match an organism to its environment. Because polyphenism is the extreme example of how development deviates from a linear genetic blueprint, it demands a genetic explanation for how environmental cues shunt development to hypothetically alternative modules. We highlight several recent advances that have begun to illuminate genetic mechanisms for polyphenism and how this recurring developmental novelty may arise. An emerging genetic knowledge of polyphenism is providing precise targets for testing hypotheses of how switch mechanisms are built-out of olfactory, nutrient-sensing, hormone-reception, and developmental and genetic buffering systems-to accommodate plasticity. Moreover, classic and new model systems are testing the genetic basis of polyphenism's proposed causal roles in evolutionary change.

Selective forces acting during multi-domain protein evolution: the case of multi-domain globins.

Multi-domain proteins form the majority of proteins in eukaryotes. During their formation by tandem duplication or gene fusion, new interactions between domains may arise as a result of the structurally-forced proximity of domains. The proper function of the formed proteins likely required the molecular adjustment of these stress zones by specific amino acid replacements, which should be detectable by the molecular signature of selection that governed their changes. We used multi-domain globins from three different invertebrate lineages to investigate the selective forces that acted throughout the evolution of these molecules. In the youngest of these molecules [Branchipolynoe scaleworm; original duplication ca. 60 million years (Ma)], we were able to detect some amino acids under positive selection corresponding to the initial duplication event. In older lineages (didomain globin from bivalve mollusks and nematodes), there was no evidence of amino acid positions under positive selection, possibly the result of accumulated non-adaptative mutations since the original duplication event (165 and 245 Ma, respectively). Some amino acids under positive selection were sometimes detected in later branches, either after speciation events, or after the initial duplication event. In Branchipolynoe, the position of the amino acids under positive selection on a 3D model suggests some of them are located at the interface between two domains; while others are locate in the heme pocket.