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.

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-O 2 affinity for O 2 loading in the pulmonary circulation and O 2 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-O 2 affinity and native elevation. Protein engineering experiments and ancestral-state reconstructions revealed that this correlation is attributable to derived increases in Hb-O 2 affinity in highland lineages, as well as derived reductions in Hb-O 2 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-O 2 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.high-altitude adaptation | hypoxia | parallel evolution | protein evolution | epistasis

Background:The functional significance of hemoglobin heterogeneity remains a mystery. Results: In adult birds, the HbD isoform (related to embryonic hemoglobin) exhibits distinct oxygenation properties relative to the major HbA isoform. Conclusion: Substitutions that distinguish HbD from HbA are not shared with embryonic hemoglobin. Significance: Differences between isoforms stem from derived (nonancestral) changes in duplicated genes, not from the retention of an ancestral condition.

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.

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.