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 α-globin paralogs on chromosome 8 and two β-globin paralogs on chromosome 1. Differences in O 2 -binding affinity of the alternative β-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 β-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 β-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.