Chandrasekhar Natarajan

Epistasis among adaptive mutations in deer mouse hemoglobin.

Holding Your Breath Hemoglobin and myoglobin are widely responsible for oxygen transport and storage (see the Perspective by Rezende). The ability of diving mammals to obtain enough oxygen to support extended dives and foraging is largely dependent on muscle myoglobin (Mb) content. Mirceta et al. (p. 1303) found that in mammalian lineages with an aquatic or semiaquatic lifestyle, Mb net charge increases, which may represent an adaptation to inhibit self-association of Mb at high intracellular concentrations. Epistasis results from nonadditive genetic interactions and can affect phenotypic evolution. Natarajan et al. (p. 1324) found that epistatic interactions were able to explain the increased hemoglobin oxygen-binding affinity observed in deer mice populations at high altitude. In mammals, the offloading of oxygen from hemoglobin is facilitated by a reduction in the bloods pH, driven by metabolically produced CO2. However, in fish, a reduction in blood pH reduces oxygen carrying capacity of hemoglobin. Rummer et al. (p. 1327) implanted fiber optic oxygen sensors within the muscles of rainbow trout and found that elevated CO2 levels in the water led to acidosis and elevated oxygen tensions. Deer mice have discovered that mutations distant from the oxygen-binding site help them live at high altitude. [Also see Perspective by Rezende] Epistatic interactions between mutant sites in the same protein can exert a strong influence on pathways of molecular evolution. We performed protein engineering experiments that revealed pervasive epistasis among segregating amino acid variants that contribute to adaptive functional variation in deer mouse hemoglobin (Hb). Amino acid mutations increased or decreased Hb-O2 affinity depending on the allelic state of other sites. Structural analysis revealed that epistasis for Hb-O2 affinity and allosteric regulatory control is attributable to indirect interactions between structurally remote sites. The prevalence of sign epistasis for fitness-related biochemical phenotypes has important implications for the evolutionary dynamics of protein polymorphism in natural populations.

Repeated elevational transitions in hemoglobin function during the evolution of Andean hummingbirds

Significance Hummingbirds have exceedingly high oxygen demands because of their elevated rates of aerobic metabolism, and yet they thrive in high-altitude environments in the Andes where oxygen is scarce. Here we report the finding that when hummingbird species colonized new elevational zones, evolutionary changes in the respiratory properties of hemoglobin were repeatedly mediated by the same amino acid replacements. Specifically, ancestral sequence reconstruction and protein engineering experiments revealed that parallel adaptation of hemoglobin function in multiple species is attributable to repeated amino acid replacements at a single pair of interacting sites. This striking parallelism at the molecular level suggests a surprising degree of reproducibility and predictability in adaptive protein evolution. 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.

Predictable convergence in hemoglobin function has unpredictable molecular underpinnings

Expect the unexpected In convergent evolution, similar environmental conditions produce similar sets of adaptations. Does similar convergence exist in the molecular underpinnings of such morphological changes? Natarajan et al. looked across more than 50 species of birds that have adapted to different elevations to identify patterns of similarity in hemoglobin-oxygen binding affinity (see the Perspective by Bridgham). Increases in hemoglobin-oxygen binding affinity occurred in alpine species, but the molecular changes underlying the hemoglobin changes were variable. Thus, even in cases where adaptive phenotypic change is predictable, the molecular pathways to these changes may not be. Science, this issue p. 336; see also p. 289 Improvements in hemoglobin-oxygen affinity are convergent across taxa, but the genetic pathways are different. To investigate the predictability of genetic adaptation, we examined the molecular basis of convergence in hemoglobin function in comparisons involving 56 avian taxa that have contrasting altitudinal range limits. Convergent increases in hemoglobin-oxygen affinity were pervasive among high-altitude taxa, but few such changes were attributable to parallel amino acid substitutions at key residues. Thus, predictable changes in biochemical phenotype do not have a predictable molecular basis. Experiments involving resurrected ancestral proteins revealed that historical substitutions have context-dependent effects, indicating that possible adaptive solutions are contingent on prior history. Mutations that produce an adaptive change in one species may represent precluded possibilities in other species because of differences in genetic background.

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.

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.

Gene Duplication and the Evolution of Hemoglobin Isoform Differentiation in Birds

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. 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 O2 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-O2 affinity stems primarily from changes in the O2 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.

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

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.

Contribution of a mutational hot spot to hemoglobin adaptation in high-altitude Andean house wrens

Significance Within a given gene, there may be many possible mutations that are capable of producing a particular change in phenotype. However, if some sites have especially high rates of mutation to function-altering alleles, then such mutations may make disproportionate contributions to phenotypic evolution. We report the discovery that a point mutation at a highly mutable site in the β-globin gene of Andean house wrens has produced a physiologically important change in the oxygenation properties of hemoglobin (Hb). The mutant allele that confers an increased Hb–O2 affinity is present at an unusually high frequency at high altitude. These findings suggest that site-specific variation in mutation rate may exert a strong influence on the genetic basis of phenotypic evolution. 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.

Expression and Purification of Recombinant Hemoglobin in Escherichia coli

Background Recombinant DNA technologies have played a pivotal role in the elucidation of structure-function relationships in hemoglobin (Hb) and other globin proteins. Here we describe the development of a plasmid expression system to synthesize recombinant Hbs in Escherichia coli, and we describe a protocol for expressing Hbs with low intrinsic solubilities. Since the α- and β-chain Hbs of different species span a broad range of solubilities, experimental protocols that have been optimized for expressing recombinant human HbA may often prove unsuitable for the recombinant expression of wildtype and mutant Hbs of other species. Methodology/Principal Findings As a test case for our expression system, we produced recombinant Hbs of the deer mouse (Peromyscus maniculatus), a species that has been the subject of research on mechanisms of Hb adaptation to hypoxia. By experimentally assessing the combined effects of induction temperature, induction time and E. coli expression strain on the solubility of recombinant deer mouse Hbs, we identified combinations of expression conditions that greatly enhanced the yield of recombinant protein and which also increased the efficiency of post-translational modifications. Conclusion/Significance Our protocol should prove useful for the experimental study of recombinant Hbs in many non-human animals. One of the chief advantages of our protocol is that we can express soluble recombinant Hb without co-expressing molecular chaperones, and without the need for additional reconstitution or heme-incorporation steps. Moreover, our plasmid construct contains a combination of unique restriction sites that allows us to produce recombinant Hbs with different α- and β-chain subunit combinations by means of cassette mutagenesis.