Evolution of physiological performance capacities and environmental adaptation: insights from high-elevation deer mice (Peromyscus maniculatus)

Analysis of variation in whole-animal performance can shed light on causal connections between specific traits, integrated physiological capacities, and Darwinian fitness. Here, we review and synthesize information on naturally occurring variation in physiological performance capacities and how it relates to environmental adaptation in deer mice (Peromyscus maniculatus). We discuss how evolved changes in aerobic exercise capacity and thermogenic capacity have contributed to adaptation to high elevations. Comparative work on deer mice at high and low elevations has revealed evolved differences in aerobic performance capacities in hypoxia. Highland deer mice have consistently higher aerobic performance capacities under hypoxia relative to lowland natives, consistent with the idea that it is beneficial to have a higher maximal metabolic rate (as measured by the maximal rate of O₂ consumption, VO_2max) in an environment characterized by lower air temperatures and lower O₂ availability. Observed differences in aerobic performance capacities between highland and lowland deer mice stem from changes in numerous subordinate traits that alter the flux capacity of the O₂-transport system, the oxidative capacity of tissue mitochondria, and the relationship between O₂ consumption and ATP synthesis. Many such changes in physiological phenotype are associated with hypoxia-induced changes in gene expression. Research on natural variation in whole-animal performance forms a nexus between physiological ecology and evolutionary biology that requires insight into the natural history of the study species.

Remodeling Ancestral Phenotypic Plasticity in Local Adaptation: A New Framework to Explore the Role of Genetic Compensation in the Evolution of Homeostasis

Phenotypic plasticity is not universally adaptive. In certain cases, plasticity can result in phenotypic shifts that reduce fitness relative to the un-induced state. A common cause of such maladaptive plasticity is the co-option of ancestral developmental and physiological response systems to meet novel challenges. Because these systems evolved to meet specific challenges in an ancestral environment (e.g., localized and transient hypoxia), their co-option to meet a similar, but novel, stressor (e.g., reductions in ambient pO2 at high elevation) can lead to misdirected responses that reduce fitness. In such cases, natural selection should act to remodel phenotypic plasticity to suppress the expression of these maladaptive responses. Because these maladaptive responses reduce the fitness of colonizers in new environments, this remodeling of ancestral plasticity may be among the earliest steps in adaptive walks toward new local optima. Genetic compensation has been proposed as a general form of adaptive evolution that leads to the suppression of maladaptive plasticity to restore the ancestral trait value in the face of novel stimuli. Given their central role in the regulation of basic physiological functions, we argue that genetic compensation may often be achieved by modifications of homeostatic regulatory systems. We further suggest that genetic compensation to modify homeostatic systems can be achieved by two alternative strategies that differ in their mechanistic underpinnings; to our knowledge, these strategies have not been formally recognized by previous workers. We then consider how the mechanistic details of these alternative strategies may constrain their evolution. These considerations lead us to argue that genetic compensation is most likely to evolve by compensatory physiological changes that safeguard internal homeostatic conditions to prevent the expression of maladaptive portions of conserved reaction norms, rather than direct evolution of plasticity itself. Finally, we outline a simple experimental framework to test this hypothesis. Our goal is to stimulate research aimed at providing a deeper mechanistic understanding of whether and how phenotypic plasticity can be remodeled following environmental shifts that render ancestral responses maladaptive, an issue with increasing importance in our current era of rapid environmental change.

Divergent mitochondrial lineages arose within a large, panmictic population of the Savannah sparrow (Passerculus sandwichensis)

Unusual patterns of mtDNA diversity can reveal interesting aspects of a species' biology. However, making such inferences requires discerning among the many alternative scenarios that could underlie any given mtDNA pattern. Next-generation sequencing methods provide large, multilocus data sets with increased power to resolve unusual mtDNA patterns. A mtDNA-based phylogeography of the Savannah sparrow (Passerculus sandwichensis) previously identified two sympatric, but divergent (~2%) clades within the nominate subspecies group and a third clade that consisted of birds sampled from northwest Mexico. We revisited the phylogeography of this species using a population genomic data set to resolve the processes leading to the evolution of sympatric and divergent mtDNA lineages. We identified two genetic clusters in the genomic data set corresponding to (a) the nominate subspecies group and (b) northwestern Mexico birds. Following divergence, the nominate clade maintained a large, stable population, indicating that divergent mitochondrial lineages arose within a panmictic population. Simulations based on parameter estimates from this model further confirmed that this demographic history could produce observed levels of mtDNA diversity. Patterns of divergent, sympatric mtDNA lineages are frequently interpreted as admixture of historically isolated lineages. Our analyses reject this interpretation for Savannah sparrows and underscore the need for genomic data sets to resolve the evolutionary mechanisms behind anomalous, locus-specific patterns.

Maladaptive phenotypic plasticity in cardiac muscle growth is suppressed in high-altitude deer mice

How often phenotypic plasticity acts to promote or inhibit adaptive evolution is an ongoing debate among biologists. Recent work suggests that adaptive phenotypic plasticity promotes evolutionary divergence, though several studies have also suggested that maladaptive plasticity can potentiate adaptation. The role of phenotypic plasticity, adaptive, or maladaptive, in evolutionary divergence remains controversial. We examined the role of plasticity in evolutionary divergence between two species of Peromyscus mice that differ in native elevations. We used cardiac mass as a model phenotype, since ancestral hypoxia-induced responses of the heart may be both adaptive and maladaptive at high-altitude. While left ventricle growth should enhance oxygen delivery to tissues, hypertrophy of the right ventricle can lead to heart failure and death. We compared left- and right-ventricle plasticity in response to hypoxia between captive-bred P. leucopus (representing the ancestral lowland condition) and P. maniculatus from high-altitude. We found that maladaptive ancestral plasticity in right ventricle hypertrophy is reduced in high-altitude deer mice. Analysis of the heart transcriptome suggests that changes in expression of inflammatory signaling genes, particularly interferon regulatory factors, contribute to the suppression of right ventricle hypertrophy. We found weak evidence that adaptive plasticity of left ventricle mass contributes to evolution. Our results suggest that selection to suppress ancestral maladaptive plasticity plays a role in adaptation.

High-altitude shorebird migration in the absence of topographical barriers: avoiding high air temperatures and searching for profitable winds

Nearly 20% of all bird species migrate between breeding and nonbreeding sites annually. Their migrations include storied feats of endurance and physiology, from non-stop trans-Pacific crossings to flights at the cruising altitudes of jetliners. Despite intense interest in these performances, there remains great uncertainty about which factors most directly influence bird behaviour during migratory flights. We used GPS trackers that measure an individual's altitude and wingbeat frequency to track the migration of black-tailed godwits (Limosa limosa) and identify the abiotic factors influencing their in-flight migratory behaviour. We found that godwits flew at altitudes above 5000 m during 21% of all migratory flights, and reached maximum flight altitudes of nearly 6000 m. The partial pressure of oxygen at these altitudes is less than 50% of that at sea level, yet these extremely high flights occurred in the absence of topographical barriers. Instead, they were associated with high air temperatures at lower altitudes and increasing wind support at higher altitudes. Our results therefore suggest that wind, temperature and topography all play a role in determining migratory behaviour, but that their relative importance is context dependent. Extremely high-altitude flights may thus not be especially rare, but they may only occur in very specific environmental contexts.

Understanding Evolutionary Impacts of Seasonality: An Introduction to the Symposium

Seasonality is a critically important aspect of environmental variability, and strongly shapes all aspects of life for organisms living in highly seasonal environments. Seasonality has played a key role in generating biodiversity, and has driven the evolution of extreme physiological adaptations and behaviors such as migration and hibernation. Fluctuating selection pressures on survival and fecundity between summer and winter provide a complex selective landscape, which can be met by a combination of three outcomes of adaptive evolution: genetic polymorphism, phenotypic plasticity, and bet-hedging. Here, we have identified four important research questions with the goal of advancing our understanding of evolutionary impacts of seasonality. First, we ask how characteristics of environments and species will determine which adaptive response occurs. Relevant characteristics include costs and limits of plasticity, predictability, and reliability of cues, and grain of environmental variation relative to generation time. A second important question is how phenological shifts will amplify or ameliorate selection on physiological hardiness. Shifts in phenology can preserve the thermal niche despite shifts in climate, but may fail to completely conserve the niche or may even expose life stages to conditions that cause mortality. Considering distinct environmental sensitivities of life history stages will be key to refining models that forecast susceptibility to climate change. Third, we must identify critical physiological phenotypes that underlie seasonal adaptation and work toward understanding the genetic architectures of these responses. These architectures are key for predicting evolutionary responses. Pleiotropic genes that regulate multiple responses to changing seasons may facilitate coordination among functionally related traits, or conversely may constrain the expression of optimal phenotypes. Finally, we must advance our understanding of how changes in seasonal fluctuations are impacting ecological interaction networks. We should move beyond simple dyadic interactions, such as predator prey dynamics, and understand how these interactions scale up to affect ecological interaction networks. As global climate change alters many aspects of seasonal variability, including extreme events and changes in mean conditions, organisms must respond appropriately or go extinct. The outcome of adaptation to seasonality will determine responses to climate change.

Comparative Transcriptomics of Seasonal Phenotypic Flexibility in Two North American Songbirds

Phenotypic flexibility allows organisms to reversibly alter their phenotypes to match the changing demands of seasonal environments. Because phenotypic flexibility is mediated, at least in part, by changes in gene regulation, comparative transcriptomic studies can provide insights into the mechanistic underpinnings of seasonal phenotypic flexibility, and the extent to which regulatory responses to changing seasons are conserved across species. To begin to address these questions, we sampled individuals of two resident North American songbird species, American goldfinch (Spinus tristis) and black-capped chickadee (Poecile atricapillus) in summer and winter to measure seasonal variation in pectoralis transcriptomic profiles and to identify conserved and species-specific elements of these seasonal profiles. We found that very few genes exhibited divergent responses to changes in season between species, and instead, a core set of over 1200 genes responded to season concordantly in both species. Moreover, several key metabolic pathways, regulatory networks, and gene functional classes were commonly recruited to induce seasonal phenotypic shifts in these species. The seasonal transcriptomic responses mirror winter increases in pectoralis mass and cellular metabolic intensity documented in previous studies of both species, suggesting that these seasonal phenotypic responses are due in part to changes in gene expression. Despite growing evidence of muscle nonshivering thermogenesis (NST) in young precocial birds, we did not find strong evidence of upregulation of genes putatively involved in NST during winter in either species, suggesting that seasonal modification of muscular NST is not a prominent contributor to winter increases in thermogenic capacity for adult passerine birds. Together, these results provide the first comprehensive overview of potential common regulatory mechanisms underlying seasonally flexible phenotypes in wild, free-ranging birds.

Life-history characteristics influence physiological strategies to cope with hypoxia in Himalayan birds

Hypobaric hypoxia at high elevation represents an important physiological stressor for montane organisms, but optimal physiological strategies to cope with hypoxia may vary among species with different life histories. Montane birds exhibit a range of migration patterns; elevational migrants breed at high elevations but winter at low elevations or migrate further south, while high-elevation residents inhabit the same elevation throughout the year. Optimal physiological strategies to cope with hypoxia might therefore differ between species that exhibit these two migratory patterns, because they differ in the amount time spent at high elevation. We examined physiological parameters associated with blood-oxygen transport (haemoglobin concentration and haematocrit, i.e. the proportion of red blood cells in blood) in nine species of elevational migrants and six species of high-elevation residents that were sampled along a 2200 m (1000-3200 m) elevational gradient. Haemoglobin concentration increased with elevation within species regardless of migratory strategy, but it was only significantly correlated with haematocrit in elevational migrants. Surprisingly, haemoglobin concentration was not correlated with haematocrit in high-elevation residents, and these species exhibited higher mean cellular haemoglobin concentration than elevational migrants. Thus, alternative physiological strategies to regulate haemoglobin concentration and blood O₂ carrying capacity appear to differ among birds with different annual elevational movement patterns.

Circulatory mechanisms underlying adaptive increases in thermogenic capacity in high-altitude deer mice

We examined the circulatory mechanisms underlying adaptive increases in thermogenic capacity in deer mice (Peromyscus maniculatus) native to the cold hypoxic environment at high altitudes. Deer mice from high- and low-altitude populations were born and raised in captivity to adulthood, and then acclimated to normoxia or hypobaric hypoxia (simulating hypoxia at ∼4300 m). Thermogenic capacity [maximal O₂ consumption (V̇_O2,max), during cold exposure] was measured in hypoxia, along with arterial O₂ saturation (Sa_O2 ) and heart rate (f_H). Hypoxia acclimation increased V̇_O2,max by a greater magnitude in highlanders than in lowlanders. Highlanders also had higher Sa_O2  and extracted more O₂ from the blood per heartbeat (O₂ pulse=V̇_O2,max/f_H). Hypoxia acclimation increased f_H, O₂ pulse and capillary density in the left ventricle of the heart. Our results suggest that adaptive increases in thermogenic capacity involve integrated functional changes across the O₂ cascade that augment O₂ circulation and extraction from the blood.

Acclimation to hypoxia increases carbohydrate use during exercise in high-altitude deer mice

The low O₂ experienced at high altitude is a significant challenge to effective aerobic locomotion, as it requires sustained tissue O₂ delivery in addition to the appropriate allocation of metabolic substrates. Here, we tested whether high- and low-altitude deer mice (Peromyscus maniculatus) have evolved different acclimation responses to hypoxia with respect to muscle metabolism and fuel use during submaximal exercise. Using F₁ generation high- and low-altitude deer mice that were born and raised in common conditions, we assessed 1) fuel use during exercise, 2) metabolic enzyme activities, and 3) gene expression for key transporters and enzymes in the gastrocnemius. After hypoxia acclimation, highland mice showed a significant increase in carbohydrate oxidation and higher relative reliance on this fuel during exercise at 75% maximal O₂ consumption. Compared with lowland mice, highland mice had consistently higher activities of oxidative and fatty acid oxidation enzymes in the gastrocnemius. In contrast, only after hypoxia acclimation did activities of hexokinase increase significantly in the muscle of highland mice to levels greater than lowland mice. Highland mice also responded to acclimation with increases in muscle gene expression for hexokinase 1 and 2 genes, whereas both populations increased mRNA expression for glucose transporters. Changes in skeletal muscle with acclimation suggest that highland mice had an increased capacity for the uptake and oxidation of circulatory glucose. Our results demonstrate that highland mice have evolved a distinct mode of hypoxia acclimation that involves an increase in carbohydrate use during exercise.

Transcriptomic plasticity in brown adipose tissue contributes to an enhanced capacity for nonshivering thermogenesis in deer mice

For small mammals living at high altitude, aerobic heat generation (thermogenesis) is essential for survival during prolonged periods of cold, but is severely impaired under conditions of hypobaric hypoxia. Recent studies in deer mice (Peromyscus maniculatus) reveal adaptive enhancement of thermogenesis in high- compared to low-altitude populations under hypoxic cold stress, an enhancement that is attributable to modifications in the aerobic metabolism of muscles used in shivering. However, because small mammals rely heavily on nonshivering mechanisms for cold acclimatization, we tested for evidence of adaptive divergence in nonshivering thermogenesis (NST) under hypoxia. To do so, we measured NST and characterized transcriptional profiles of brown adipose tissue (BAT) in high- and low-altitude deer mice that were (i) wild-caught and acclimatized to their native altitude, and (ii) born and reared under common garden conditions at low elevation. We found that NST performance under hypoxia is enhanced in wild-caught, high-altitude deer mice, a difference that is associated with increased expression of coregulated genes that influence several physiological traits. These traits include vascularization and O2 supply to BAT, brown adipocyte proliferation and the uncoupling of oxidative phosphorylation from ATP synthesis in the generation of heat. Our results suggest that acclimatization to hypoxic cold stress is facilitated by enhancement of nonshivering heat production, which is driven by regulatory plasticity in a suite of genes that influence intersecting physiological pathways.

Divergence and Functional Degradation of a Sex Chromosome-like Supergene

A major challenge in biology is to understand the genetic basis of adaptation. One compelling idea is that groups of tightly linked genes (i.e., "supergenes" [1, 2]) facilitate adaptation in suites of traits that determine fitness. Despite their likely importance, little is known about how alternate supergene alleles arise and become differentiated, nor their ultimate fate within species. Herein we address these questions by investigating the evolutionary history of a supergene in white-throated sparrows, Zonotrichia albicollis. This species comprises two morphs, tan and white, that differ in pigmentation and components of social behavior [3-5]. Morph is determined by alternative alleles at a balanced >100-Mb inversion-based supergene, providing a unique system for studying gene-behavior relationships. Using over two decades of field data, we document near-perfect disassortative mating among morphs, as well as the fitness consequences of rare assortative mating. We use de novo whole-genome sequencing coupled with population- and phylogenomic data to show that alternate supergene alleles are highly divergent at over 1,000 genes and that these alleles originated prior to the split of Z. albicollis from its sister species and may be polymorphic in Z. albicollis due to a past hybridization event. We provide evidence that the "white" allele may be degrading, similar to neo-Y/W sex chromosomes. We further show that the "tan" allele has surprisingly low levels of genetic diversity yet does not show several canonical signatures of recurrent positive selection. We discuss these results in the context of the origin, molecular evolution, and possible fate of this remarkable polymorphism.

Functional Genomic Insights into Regulatory Mechanisms of High-Altitude Adaptation

Recent studies of indigenous human populations at high altitude have provided proof-of-principle that genome scans of DNA polymorphism can be used to identify candidate loci for hypoxia adaptation. When integrated with experimental analyses of physiological phenotypes, genome-wide surveys of DNA polymorphism and tissue-specific transcriptional profiles can provide insights into actual mechanisms of adaptation. It has been suggested that adaptive phenotypic evolution is largely mediated by cis-regulatory changes in genes that are located at integrative control points in regulatory networks. This hypothesis can be tested by conducting transcriptomic analyses of hypoxic signaling pathways in conjunction with experimental measures of vascular oxygen supply and metabolic pathway flux. Such studies may reveal whether the architecture of gene regulatory networks can be used to predict which loci (and which types of loci) are likely to be "hot spots" for adaptive physiological evolution. Functional genomic studies of deer mice (Peromyscus maniculatus) demonstrate how the integrated analysis of variation in tissue-specific transcriptomes, whole-animal physiological performance, and various subordinate traits can yield insights into the mechanistic underpinnings of high-altitude adaptation.

Adaptive Modifications of Muscle Phenotype in High-Altitude Deer Mice Are Associated with Evolved Changes in Gene Regulation

At high-altitude, small mammals are faced with the energetic challenge of sustaining thermogenesis and aerobic exercise in spite of the reduced O2 availability. Under conditions of hypoxic cold stress, metabolic demands of shivering thermogenesis and locomotion may require enhancements in the oxidative capacity and O2 diffusion capacity of skeletal muscle to compensate for the diminished tissue O2 supply. We used common-garden experiments involving highland and lowland deer mice (Peromyscus maniculatus) to investigate the transcriptional underpinnings of genetically based population differences and plasticity in muscle phenotype. We tested highland and lowland mice that were sampled in their native environments as well as lab-raised F1 progeny of wild-caught mice. Experiments revealed that highland natives had consistently greater oxidative fiber density and capillarity in the gastrocnemius muscle. RNA sequencing analyses revealed population differences in transcript abundance for 68 genes that clustered into two discrete transcriptional modules, and a large suite of transcripts (589 genes) with plastic expression patterns that clustered into five modules. The expression of two transcriptional modules was correlated with the oxidative phenotype and capillarity of the muscle, and these phenotype-associated modules were enriched for genes involved in energy metabolism, muscle plasticity, vascular development, and cell stress response. Although most of the individual transcripts that were differentially expressed between populations were negatively correlated with muscle phenotype, several genes involved in energy metabolism (e.g., Ckmt1, Ehhadh, Acaa1a) and angiogenesis (Notch4) were more highly expressed in highlanders, and the regulators of mitochondrial biogenesis, PGC-1α (Ppargc1a) and mitochondrial transcription factor A (Tfam), were positively correlated with muscle oxidative phenotype. These results suggest that evolved population differences in the oxidative capacity and capillarity of skeletal muscle involved expression changes in a small suite of coregulated genes.

Regulatory mechanisms of metabolic flexibility in the dark-eyed junco (Junco hyemalis)

Small temperate birds reversibly modify their aerobic performance to maintain thermoregulatory homeostasis under seasonally changing environmental conditions and these physiological adjustments may be attributable to changes in the expression of genes in the underlying regulatory networks. Here, we report the results of an experimental procedure designed to gain insight into the fundamental mechanisms of metabolic flexibility in the dark-eyed junco (Junco hyemalis). We combined genomic transcriptional profiles with measures of metabolic enzyme activities and whole-animal thermogenic performance from juncos exposed to four 6-week acclimation treatments that varied in temperature (cold, 3°C; warm, 24°C) and photoperiod (short day, 8 h light:16 h dark; long day, 16 h light:8 h dark). Cold-acclimated birds increased thermogenic capacity compared with warm-acclimated birds, and this enhanced performance was associated with upregulation of genes involved in muscle hypertrophy, angiogenesis, and lipid transport and oxidation, as well as with catabolic enzyme activities. These physiological changes occurred over ecologically relevant timescales, suggesting that birds make regulatory adjustments to interacting, hierarchical pathways in order to seasonally enhance thermogenic capacity.

High-altitude ancestry and hypoxia acclimation have distinct effects on exercise capacity and muscle phenotype in deer mice

The hypoxic and cold environment at high altitudes requires that small mammals sustain high rates of O2 transport for exercise and thermogenesis while facing a diminished O2 availability. We used laboratory-born and -raised deer mice (Peromyscus maniculatus) from highland and lowland populations to determine the interactive effects of ancestry and hypoxia acclimation on exercise performance. Maximal O₂consumption (V̇o(2max)) during exercise in hypoxia increased after hypoxia acclimation (equivalent to the hypoxia at ∼4,300 m elevation for 6-8 wk) and was consistently greater in highlanders than in lowlanders. V̇o(2max) during exercise in normoxia was not affected by ancestry or acclimation. Highlanders also had consistently greater capillarity, oxidative fiber density, and maximal activities of oxidative enzymes (cytochrome c oxidase and citrate synthase) in the gastrocnemius muscle, lower lactate dehydrogenase activity in the gastrocnemius, and greater cytochrome c oxidase activity in the diaphragm. Hypoxia acclimation did not affect any of these muscle traits. The unique gastrocnemius phenotype of highlanders was associated with higher mRNA and protein abundances of peroxisome proliferator-activated receptor γ (PPARγ). Vascular endothelial growth factor (VEGFA) transcript abundance was lower in highlanders, and hypoxia acclimation reduced the expression of numerous genes that regulate angiogenesis and energy metabolism, in contrast to the observed population differences in muscle phenotype. Lowlanders exhibited greater increases in blood hemoglobin content, hematocrit, and wet lung mass (but not dry lung mass) than highlanders after hypoxia acclimation. Genotypic adaptation to high altitude, therefore, improves exercise performance in hypoxia by mechanisms that are at least partially distinct from those underlying hypoxia acclimation.

Spatially variable coevolution between a haemosporidian parasite and the MHC of a widely distributed passerine

The environment shapes host-parasite interactions, but how environmental variation affects the diversity and composition of parasite-defense genes of hosts is unresolved. In vertebrates, the highly variable major histocompatibility complex (MHC) gene family plays an essential role in the adaptive immune system by recognizing pathogen infection and initiating the cellular immune response. Investigating MHC-parasite associations across heterogeneous landscapes may elucidate the role of spatially fluctuating selection in the maintenance of high levels of genetic variation at the MHC. We studied patterns of association between an avian haemosporidian blood parasite and the MHC of rufous-collared sparrows (Zonotrichia capensis) that inhabit environments with widely varying haemosporidian infection prevalence in the Peruvian Andes. MHC diversity peaked in populations with high infection prevalence, although intra-individual MHC diversity was not associated with infection status. MHC nucleotide and protein sequences associated with infection absence tended to be rare, consistent with negative frequency-dependent selection. We found an MHC variant associated with a ∽26% decrease in infection probability at middle elevations (1501-3100 m) where prevalence was highest. Several other variants were associated with a significant increase in infection probability in low haemosporidian prevalence environments, which can be interpreted as susceptibility or quantitative resistance. Our study highlights important challenges in understanding MHC evolution in natural systems, but may point to a role of negative frequency-dependent selection and fluctuating spatial selection in the evolution of Z. capensisMHC.

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.

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.