Genetic differences in hemoglobin function between highland and lowland deer mice

Experimental approaches to evaluate the contributions of candidate protein-coding mutations to phenotypic evolution

Identifying mechanisms of molecular adaptation can provide important insights into the process of phenotypic evolution, but it can be exceedingly difficult to quantify the phenotypic effects of specific mutational changes. To verify the adaptive significance of genetically based changes in protein function, it is necessary to document functional differences between the products of derived and wild-type alleles and to demonstrate that such differences impinge on higher-level physiological processes (and ultimately, fitness). In the case of metabolic enzymes, this requires documenting in vivo differences in reaction rate that give rise to differences in flux through the pathway in which the enzymes function. These measured differences in pathway flux should then give rise to differences in cellular or systemic physiology that affect fitness-related variation in whole-organism performance. Efforts to establish these causal connections between genotype, phenotype, and fitness require experiments that carefully control for environmental variation and background genetic variation. Here, we discuss experimental approaches to evaluate the contributions of amino-acid mutations to adaptive phenotypic change. We discuss conceptual and methodological issues associated with in vitro and in vivo studies of protein function, and the evolutionary insights that can be gleaned from such studies. We also discuss the importance of isolating the effects of individual mutations to distinguish between positively selected substitutions that directly contribute to improvements in protein function versus positively selected, compensatory substitutions that mitigate negative pleiotropic effects of antecedent changes.

Copy number polymorphism in the α-globin gene cluster of European rabbit (Oryctolagus cuniculus)

Comparative genomic studies have revealed that mammals typically possess two or more tandemly duplicated copies of the α-globin (HBA) gene. The domestic rabbit represents an exception to this general rule, as this species was found to possess a single HBA gene. Previous electrophoretic surveys of HBA polymorphism in natural populations of the European rabbit (Oryctolagus cuniculus) revealed extensive geographic variation in the frequencies of three main electromorphs. The variation in frequency of two electromorphs is mainly partitioned between two distinct subspecies of European rabbit, and a third is restricted to the hybrid zone between the two rabbit subspecies in Iberia. Here we report the results of a survey of nucleotide polymorphism, which revealed HBA copy number polymorphism in Iberian populations of the European rabbit. By characterizing patterns of HBA polymorphism in populations from the native range of the European rabbit, we were able to identify the specific amino-acid substitutions that distinguish the previously characterized electromorphs. Within the hybrid zone, we observed the existence of a second HBA gene duplicate, named HBA2, that mostly represents a novel sequence haplotype, which occurs in higher frequency within the hybrid zone, and thus appears to have arisen in hybrids of the two distinct subspecies. Although this novel gene is also present in other wild Iberian populations, it is almost absent from French populations, which suggest a recent ancestry, associated with the establishment of the post-Pleistocene contact zone between the two European rabbit subspecies.

Oxygenation properties and oxidation rates of mouse hemoglobins that differ in reactive cysteine content

House mice (genus Mus) harbor extensive allelic variation at two tandemly duplicated genes that encode the β-chain subunits of adult hemoglobin (Hb). Alternative haplotypes differ in the level of sequence divergence between the two β-globin gene duplicates: the Hbb(d) and Hbb(p) haplotypes harbor two structurally distinct β-globin genes, whereas the Hbb(s) haplotype harbors two β-globin duplicates that are identical in sequence. One especially salient difference between the s-type Hbs relative to the d- and p-type Hbs relates to the number of reactive β-chain cysteine residues. In addition to the highly conserved cysteine residue at β93, the d- and p-type Hbs contain an additional reactive cysteine residue at β13. To assess the functional consequences of allelic variation in β-globin cysteine content, we measured O(2)-binding properties and H(2)O(2)-induced oxidation rates of mono- and dicysteinyl β-Hbs from 4 different inbred strains of mice: C57BL/6J, BALB/cJ, MSM/Ms, and CAROLI/EiJ. The experiments revealed that purified Hbs from the various mouse strains did not exhibit substantial variation in O(2)-binding properties, but s-type Hb (which contains a single reactive β-chain cysteine residue) was far more readily oxidized to Fe(3+) metHb by H(2)O(2) than other mouse Hbs that contain two reactive β-chain cysteine residues. These results suggest that the possession of an additional reactive cysteine residue may protect against metHb formation under oxidizing conditions. The allelic differences in β-globin cysteine content could affect aspects of redox signaling and oxidative/nitrosative stress responses that are mediated by Hb-S-nitrosylation and Hb-S-glutathionylation pathways.

Altitudinal variation at duplicated β-globin genes in deer mice: effects of selection, recombination, and gene conversion

Spatially varying selection on a given polymorphism is expected to produce a localized peak in the between-population component of nucleotide diversity, and theory suggests that the chromosomal extent of elevated differentiation may be enhanced in cases where tandemly linked genes contribute to fitness variation. An intriguing example is provided by the tandemly duplicated β-globin genes of deer mice (Peromyscus maniculatus), which contribute to adaptive differentiation in blood-oxygen affinity between high- and low-altitude populations. Remarkably, the two β-globin genes segregate the same pair of functionally distinct alleles due to a history of interparalog gene conversion and alleles of the same functional type are in perfect coupling-phase linkage disequilibrium (LD). Here we report a multilocus analysis of nucleotide polymorphism and LD in highland and lowland mice with different genetic backgrounds at the β-globin genes. The analysis of haplotype structure revealed a paradoxical pattern whereby perfect LD between the two β-globin paralogs (which are separated by 16.2 kb) is maintained in spite of the fact that LD within both paralogs decays to background levels over physical distances of less than 1 kb. The survey of nucleotide polymorphism revealed that elevated levels of altitudinal differentiation at each of the β-globin genes drop away quite rapidly in the external flanking regions (upstream of the 5' paralog and downstream of the 3' paralog), but the level of differentiation remains unexpectedly high across the intergenic region. Observed patterns of diversity and haplotype structure are difficult to reconcile with expectations of a two-locus selection model with multiplicative fitness.

Molecular spandrels: tests of adaptation at the genetic level

Although much progress has been made in identifying the genes (and, in rare cases, mutations) that contribute to phenotypic variation, less is known about the effects that these genes have on fitness. Nonetheless, genes are commonly labelled as 'adaptive' if an allele has been shown to affect a phenotype with known or suspected functional importance or if patterns of nucleotide variation at the locus are consistent with positive selection. In these cases, the 'adaptive' designation may be premature and may lead to incorrect conclusions about the relationships between gene function and fitness. Experiments to test targets and agents of natural selection within a genomic context are necessary for identifying the adaptive consequences of individual alleles.

High altitude adaptation in Daghestani populations from the Caucasus

We have surveyed 15 high-altitude adaptation candidate genes for signals of positive selection in North Caucasian highlanders using targeted re-sequencing. A total of 49 unrelated Daghestani from three ethnic groups (Avars, Kubachians, and Laks) living in ancient villages located at around 2,000 m above sea level were chosen as the study population. Caucasian (Adygei living at sea level, N = 20) and CEU (CEPH Utah residents with ancestry from northern and western Europe; N = 20) were used as controls. Candidate genes were compared with 20 putatively neutral control regions resequenced in the same individuals. The regions of interest were amplified by long-PCR, pooled according to individual, indexed by adding an eight-nucleotide tag, and sequenced using the Illumina GAII platform. 1,066 SNPs were called using false discovery and false negative thresholds of ~6%. The neutral regions provided an empirical null distribution to compare with the candidate genes for signals of selection. Two genes stood out. In Laks, a non-synonymous variant within HIF1A already known to be associated with improvement in oxygen metabolism was rediscovered, and in Kubachians a cluster of 13 SNPs located in a conserved intronic region within EGLN1 showing high population differentiation was found. These variants illustrate both the common pathways of adaptation to high altitude in different populations and features specific to the Daghestani populations, showing how even a mildly hypoxic environment can lead to genetic adaptation.

Developmental regulation of hemoglobin synthesis in the green anole lizard Anolis carolinensis

Tetrapod vertebrates possess multiple α- and β-like globin genes that are ontogenetically regulated, such that functionally distinct hemoglobin (Hb) isoforms are synthesized during different stages of development. The α- and β-like globin genes of amphibians, birds and mammals are differentially expressed during embryonic development and postnatal life, but little is known about the developmental regulation of globin gene expression in non-avian reptiles. Here we report an investigation into the developmental regulation of Hb synthesis in the green anole lizard Anolis carolinensis. We tested two hypotheses derived from comparative genomic studies of the globin gene clusters in tetrapod vertebrates. First, we tested whether the product of the Anolis α(D)-globin gene is incorporated into embryonic Hb, thereby performing the role that would normally be performed by the embyronic α(E)-globin gene (which has been deleted from the green anole genome). Second, we tested whether two 'lizard-specific' β-globin paralogs have independently evolved a division of labor between an early-expressed embryonic gene and a later-expressed adult gene. Results of a proteomic analysis revealed that α- and β-like globin genes of the anole are differentially expressed during embryonic development. However, the same repertoire of α- and β-chain Hb isoforms was expressed during all stages of development and postnatal life, and the ontogenetic shifts in isoform composition were relatively subtle. In contrast to the pattern that has been documented in other tetrapod vertebrates, it appears that the developmental regulation of Hb synthesis in the green anole lizard does not involve discrete, stage-specific switches in gene activation and gene silencing.

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.

Phylogenetic diversification of the globin gene superfamily in chordates

Phylogenetic reconstructions provide a means of inferring the branching relationships among members of multigene families that have diversified via successive rounds of gene duplication and divergence. Such reconstructions can illuminate the pathways by which particular expression patterns and protein functions evolved. For example, phylogenetic analyses can reveal cases in which similar expression patterns or functional properties evolved independently in different lineages, either through convergence, parallelism, or evolutionary reversals. The purpose of this article is to provide a robust phylogenetic framework for interpreting experimental data and for generating hypotheses about the functional evolution of globin proteins in chordate animals. To do this, we present a consensus phylogeny of the chordate globin gene superfamily. We document the relative roles of gene duplication and whole-genome duplication in fueling the functional diversification of vertebrate globins, and we unravel patterns of shared ancestry among globin genes from representatives of the three chordate subphyla (Craniata, Urochordata, and Cephalochordata). Our results demonstrate the value of integrating phylogenetic analyses with genomic analyses of conserved synteny to infer the duplicative origins and evolutionary histories of globin genes. We also discuss a number of case studies that illustrate the importance of phylogenetic information when making inferences about the evolution of globin gene expression and protein function. Finally, we discuss why the globin gene superfamily presents special challenges for phylogenetic analysis, and we describe methodological approaches that can be used to meet those challenges.

Integrating evolutionary and functional approaches to infer adaptation at specific loci

Inferences about adaptation at specific loci are often exclusively based on the static analysis of DNA sequence variation. Ideally,population-genetic evidence for positive selection serves as a stepping-off point for experimental studies to elucidate the functional significance of the putatively adaptive variation. We argue that inferences about adaptation at specific loci are best achieved by integrating the indirect, retrospective insights provided by population-genetic analyses with the more direct, mechanistic insights provided by functional experiments. Integrative studies of adaptive genetic variation may sometimes be motivated by experimental insights into molecular function, which then provide the impetus to perform population genetic tests to evaluate whether the functional variation is of adaptive significance. In other cases, studies may be initiated by genome scans of DNA variation to identify candidate loci for recent adaptation. Results of such analyses can then motivate experimental efforts to test whether the identified candidate loci do in fact contribute to functional variation in some fitness-related phenotype. Functional studies can provide corroborative evidence for positive selection at particular loci, and can potentially reveal specific molecular mechanisms of adaptation.

Phenotypic plasticity and genetic adaptation to high-altitude hypoxia in vertebrates

High-altitude environments provide ideal testing grounds for investigations of mechanism and process in physiological adaptation. In vertebrates, much of our understanding of the acclimatization response to high-altitude hypoxia derives from studies of animal species that are native to lowland environments. Such studies can indicate whether phenotypic plasticity will generally facilitate or impede adaptation to high altitude. Here, we review general mechanisms of physiological acclimatization and genetic adaptation to high-altitude hypoxia in birds and mammals. We evaluate whether the acclimatization response to environmental hypoxia can be regarded generally as a mechanism of adaptive phenotypic plasticity, or whether it might sometimes represent a misdirected response that acts as a hindrance to genetic adaptation. In cases in which the acclimatization response to hypoxia is maladaptive, selection will favor an attenuation of the induced phenotypic change. This can result in a form of cryptic adaptive evolution in which phenotypic similarity between high- and low-altitude populations is attributable to directional selection on genetically based trait variation that offsets environmentally induced changes. The blunted erythropoietic and pulmonary vasoconstriction responses to hypoxia in Tibetan humans and numerous high-altitude birds and mammals provide possible examples of this phenomenon. When lowland animals colonize high-altitude environments, adaptive phenotypic plasticity can mitigate the costs of selection, thereby enhancing prospects for population establishment and persistence. By contrast, maladaptive plasticity has the opposite effect. Thus, insights into the acclimatization response of lowland animals to high-altitude hypoxia can provide a basis for predicting how altitudinal range limits might shift in response to climate change.