Evolution of duplicated beta-globin genes and the structural basis of hemoglobin isoform differentiation in Mus

Single-cell RNA sequencing identifies Fgf23-expressing osteocytes in response to 1,25-dihydroxyvitamin D3 treatment

Fibroblast growth factor 23 (FGF23), a hormone, mainly produced by osteocytes, regulates phosphate and vitamin D metabolism. By contrast, 1,25-dihydroxyvitamin D₃, the active form of vitamin D, has been shown to enhance FGF23 production. While it is likely that osteocytes are heterogenous in terms of gene expression profiles, specific subpopulations of Fgf23-expressing osteocytes have not been identified. Single-cell RNA sequencing (scRNA-seq) technology can characterize the transcriptome of an individual cell. Recently, scRNA-seq has been used for bone tissue analysis. However, owing to technical difficulties associated with isolation of osteocytes, studies using scRNA-seq analysis to characterize FGF23-producing osteocytes are lacking. In this study, we characterized osteocytes secreting FGF23 from murine femurs in response to calcitriol (1,25-dihydroxyvitamin D₃) using scRNA-seq. We first detected Dmp1, Mepe, and Phex expression in murine osteocytes by in situ hybridization and used these as marker genes of osteocytes. After decalcification, enzyme digestion, and removal of CD45⁺ cells, femoral bone cells were subjected to scRNA-seq. We identified cell clusters containing osteocytes using marker gene expression. While Fgf23 expression was observed in some osteocytes isolated from femurs of calcitriol-injected mice, no Fgf23 expression was detected in untreated mice. In addition, the expression of several genes which are known to be changed after 1,25-dihydroxyvitamin D₃ treatment such as Ccnd2, Fn1, Igfbp7, Pdgfa, and Timp1 was also affected by calcitriol treatment in Fgf23-expressing osteocytes, but not in those lacking Fgf23 expression, even after calcitriol administration. Furthermore, box-and-whisker plots indicated that Fgf23 expression was observed in osteocytes with higher expression levels of the Fam20c, Dmp1, and Phex genes, whose inactivating mutations have been shown to cause FGF23-related hypophosphatemic diseases. These results indicate that osteocytes are heterogeneous with respect to their responsiveness to 1,25-dihydroxyvitamin D₃, and sensitivity to 1,25-dihydroxyvitamin D₃ is one of the characteristics of osteocytes with Fgf23 expression. It is likely that there is a subpopulation of osteocytes expressing several genes, including Fgf23, involved in phosphate metabolism.

The rise and fall of globins in the amphibia

The globin gene repertoire of gnathostome vertebrates is dictated by differential retention and loss of nine paralogous genes: androglobin, neuroglobin, globin X, cytoglobin, globin Y, myoglobin, globin E, and the α- and β-globins. We report the globin gene repertoire of three orders of modern amphibians: Anura, Caudata, and Gymnophiona. Combining phylogenetic and conserved synteny analysis, we show that myoglobin and globin E were lost only in the Batrachia clade, but retained in Gymnophiona. The major amphibian groups also retained different paralogous copies of globin X. None of the amphibian presented α^D-globin gene. Nevertheless, two clades of β-globins are present in all amphibians, indicating that the amphibian ancestor possessed two paralogous proto β-globins. We also show that orthologs of the gene coding for the monomeric hemoglobin found in the heart of Rana catesbeiana are present in Neobatrachia and Pelobatoidea species we analyzed. We suggest that these genes might perform myoglobin- and globin E-related functions. We conclude that the repertoire of globin genes in amphibians is dictated by both retention and loss of the paralogous genes cited above and the rise of a new globin gene through co-option of an α-globin, possibly facilitated by a prior event of transposition.

PML Nuclear Body Biogenesis, Carcinogenesis, and Targeted Therapy

Targeted therapy has become increasingly important in cancer therapy. For example, targeting the promyelocytic leukemia PML protein in leukemia has proved to be an effective treatment. PML is the core component of super-assembled structures called PML nuclear bodies (NBs). Although this nuclear megaDalton complex was first observed in the 1960s, the mechanism of its assembly remains poorly understood. We review recent breakthroughs in the PML field ranging from a revised assembly mechanism to PML-driven genome organization and carcinogenesis. In addition, we highlight that oncogenic oligomerization might also represent a promising target in the treatment of leukemias and solid tumors.

Oxygenation properties of hemoglobin and the evolutionary origins of isoform multiplicity in an amphibious air-breathing fish, the blue-spotted mudskipper (Boleophthalmus pectinirostris)

Among the numerous lineages of teleost fish that have independently transitioned from obligate water breathing to facultative air breathing, evolved properties of hemoglobin (Hb)-O₂ transport may have been shaped by the prevalence and severity of aquatic hypoxia (which influences the extent to which fish are compelled to switch to aerial respiration) as well as the anatomical design of air-breathing structures and the cardiovascular system. Here, we examined the structure and function of Hbs in an amphibious, facultative air-breathing fish, the blue-spotted mudskipper (Boleophthalmus pectinirostris). We also characterized the genomic organization of the globin gene clusters of the species and we integrated phylogenetic and comparative genomic analyses to unravel the duplicative history of the genes that encode the subunits of structurally distinct mudskipper Hb isoforms (isoHbs). The B. pectinirostris isoHbs exhibit high intrinsic O₂ affinities, similar to those of hypoxia-tolerant, water-breathing teleosts, and remarkably large Bohr effects. Genomic analysis of conserved synteny revealed that the genes that encode the α-type subunits of the two main adult isoHbs are members of paralogous gene clusters that represent products of the teleost-specific whole-genome duplication. Experiments revealed no appreciable difference in the oxygenation properties of co-expressed isoHbs in spite of extensive amino acid divergence between the alternative α-chain subunit isoforms. It therefore appears that the ability to switch between aquatic and aerial respiration does not necessarily require a division of labor between functionally distinct isoHbs with specialized oxygenation properties.

Gene Turnover and Diversification of the α- and β-Globin Gene Families in Sauropsid Vertebrates

The genes that encode the α- and β-chain subunits of vertebrate hemoglobin have served as a model system for elucidating general principles of gene family evolution, but little is known about patterns of evolution in amniotes other than mammals and birds. Here, we report a comparative genomic analysis of the α- and β-globin gene clusters in sauropsids (archosaurs and nonavian reptiles). The objectives were to characterize changes in the size and membership composition of the α- and β-globin gene families within and among the major sauropsid lineages, to reconstruct the evolutionary history of the sauropsid α- and β-globin genes, to resolve orthologous relationships, and to reconstruct evolutionary changes in the developmental regulation of gene expression. Our comparisons revealed contrasting patterns of evolution in the unlinked α- and β-globin gene clusters. In the α-globin gene cluster, which has remained in the ancestral chromosomal location, evolutionary changes in gene content are attributable to the differential retention of paralogous gene copies that were present in the common ancestor of tetrapods. In the β-globin gene cluster, which was translocated to a new chromosomal location, evolutionary changes in gene content are attributable to differential gene gains (via lineage-specific duplication events) and gene losses (via lineage-specific deletions and inactivations). Consequently, all major groups of amniotes possess unique repertoires of embryonic and postnatally expressed β-type globin genes that diversified independently in each lineage. These independently derived β-type globins descend from a pair of tandemly linked paralogs in the most recent common ancestor of sauropsids.

Rapid and Sensitive Assessment of Globin Chains for Gene and Cell Therapy of Hemoglobinopathies

The β-hemoglobinopathies sickle cell anemia and β-thalassemia are the focus of many gene-therapy studies. A key disease parameter is the abundance of globin chains because it indicates the level of anemia, likely toxicity of excess or aberrant globins, and therapeutic potential of induced or exogenous β-like globins. Reversed-phase high-performance liquid chromatography (HPLC) allows versatile and inexpensive globin quantification, but commonly applied protocols suffer from long run times, high sample requirements, or inability to separate murine from human β-globin chains. The latter point is problematic for in vivo studies with gene-addition vectors in murine disease models and mouse/human chimeras. This study demonstrates HPLC-based measurements of globin expression (1) after differentiation of the commonly applied human umbilical cord blood-derived erythroid progenitor-2 cell line, (2) in erythroid progeny of CD34+ cells for the analysis of clustered regularly interspaced short palindromic repeats/Cas9-mediated disruption of the globin regulator BCL11A, and (3) of transgenic mice holding the human β-globin locus. At run times of 8 min for separation of murine and human β-globin chains as well as of human γ-globin chains, and with routine measurement of globin-chain ratios for 12 nL of blood (tested for down to 0.75 nL) or of 300,000 in vitro differentiated cells, the methods presented here and any variant-specific adaptations thereof will greatly facilitate evaluation of novel therapy applications for β-hemoglobinopathies.

Gene Duplication and Evolutionary Innovations in Hemoglobin-Oxygen Transport

During vertebrate evolution, duplicated hemoglobin (Hb) genes diverged with respect to functional properties as well as the developmental timing of expression. For example, the subfamilies of genes that encode the different subunit chains of Hb are ontogenetically regulated such that functionally distinct Hb isoforms are expressed during different developmental stages. In some vertebrate taxa, functional differentiation between co-expressed Hb isoforms may also contribute to physiologically important divisions of labor.

Hemoglobin-oxygen affinity in high-altitude vertebrates: is there evidence for an adaptive trend?

In air-breathing vertebrates at high altitude, fine-tuned adjustments in hemoglobin (Hb)-O2 affinity provide an energetically efficient means of mitigating the effects of arterial hypoxemia. However, it is not always clear whether an increased or decreased Hb-O2 affinity should be expected to improve tissue O2 delivery under different degrees of hypoxia, due to the inherent trade-off between arterial O2 loading and peripheral O2 unloading. Theoretical results indicate that the optimal Hb-O2 affinity varies as a non-linear function of environmental O2 availability, and the threshold elevation at which an increased Hb-O2 affinity becomes advantageous depends on the magnitude of diffusion limitation (the extent to which O2 equilibration at the blood-gas interface is limited by the kinetics of O2 exchange). This body of theory provides a framework for interpreting the possible adaptive significance of evolved changes in Hb-O2 affinity in vertebrates that have colonized high-altitude environments. To evaluate the evidence for an empirical generalization and to test theoretical predictions, I synthesized comparative data in a phylogenetic framework to assess the strength of the relationship between Hb-O2 affinity and native elevation in mammals and birds. Evidence for a general trend in mammals is equivocal, but there is a remarkably strong positive relationship between Hb-O2 affinity and native elevation in birds. Evolved changes in Hb function in high-altitude birds provide one of the most compelling examples of convergent biochemical adaptation in vertebrates.

Tracing the eastward dispersal of the house mouse, Mus musculus

Here we describe recent advances in our understanding of the natural history of the house mouse, Mus musculus, with a focus on the genetic characteristics of the home territories and how this relates to prehistoric eastward movements from the predicted source areas. Recent studies of mitochondrial and nuclear gene sequences provide insight into the ancient divergence of the three subspecies groups, M. m. castaneus (CAS), M. m. domesticus (DOM), and M. m. musculus (MUS), with inferred natural habits (homelands) in central (Iran, Afghanistan, Pakistan, and India), western (western Iran), and northern (central Asia) areas, respectively. Our mitochondrial DNA and nuclear gene analyses indicate that only one local lineage of CAS extended its range via historical rapid expansion at two different times to Southeast Asia and East Asia, including Japan and southern Sakhalin. This is suggestive of a rapid range expansion of CAS out of its homeland, perhaps associated with the spread of agricultural practices in Asia. The subspecies group MUS now occurs in a large portion of northern Eurasia from eastern Europe in the West to the Japanese Islands in the East, including Uzbekistan, Kazakhstan, southern Siberia, northern China, and Korea, showing divergent patterns in terms of Mus musculus genetics, particularly in relation to nuclear gene sequences, allozymes (e.g., hemoglobin), morphological characteristics, and cytogenetic C-banding patterns. In this review article, we explain the complex spatial patterns of MUS. We postulate that two historical dispersal events took place, from two different source areas, and tentatively assign the taxon names "musculus" and "wagneri" to the two populations, which are associated with distinct genetic modules.

Gene turnover in the avian globin gene families and evolutionary changes in hemoglobin isoform expression

The apparent stasis in the evolution of avian chromosomes suggests that birds may have experienced relatively low rates of gene gain and loss in multigene families. To investigate this possibility and to explore the phenotypic consequences of variation in gene copy number, we examined evolutionary changes in the families of genes that encode the α- and β-type subunits of hemoglobin (Hb), the tetrameric α2β2 protein responsible for blood-O2 transport. A comparative genomic analysis of 52 bird species revealed that the size and membership composition of the α- and β-globin gene families have remained remarkably constant during approximately 100 My of avian evolution. Most interspecific variation in gene content is attributable to multiple independent inactivations of the α(D)-globin gene, which encodes the α-chain subunit of a functionally distinct Hb isoform (HbD) that is expressed in both embryonic and definitive erythrocytes. Due to consistent differences in O2-binding properties between HbD and the major adult-expressed Hb isoform, HbA (which incorporates products of the α(A)-globin gene), recurrent losses of α(D)-globin contribute to among-species variation in blood-O2 affinity. Analysis of HbA/HbD expression levels in the red blood cells of 122 bird species revealed high variability among lineages and strong phylogenetic signal. In comparison with the homologous gene clusters in mammals, the low retention rate for lineage-specific gene duplicates in the avian globin gene clusters suggests that the developmental regulation of Hb synthesis in birds may be more highly conserved, with orthologous genes having similar stage-specific expression profiles and similar functional properties in disparate taxa.

Adaptive phylogeography: functional divergence between haemoglobins derived from different glacial refugia in the bank vole

Over the years, researchers have used presumptively neutral molecular variation to infer the origins of current species' distributions in northern latitudes (especially Europe). However, several reported examples of genic and chromosomal replacements suggest that end-glacial colonizations of particular northern areas may have involved genetic input from different source populations at different times, coupled with competition and selection. We investigate the functional consequences of differences between two bank vole (Clethrionomys glareolus) haemoglobins deriving from different glacial refugia, one of which partially replaced the other in Britain during end-glacial climate warming. This allows us to examine their adaptive divergence and hence a possible role of selection in the replacement. We determine the amino acid substitution Ser52Cys in the major expressed β-globin gene as the allelic difference. We use structural modelling to reveal that the protein environment renders the 52Cys thiol a highly reactive functional group and we show its reactivity in vitro. We demonstrate that possessing the reactive thiol in haemoglobin increases the resistance of bank vole erythrocytes to oxidative stress. Our study thus provides striking evidence for physiological differences between products of genic variants that spread at the expense of one another during colonization of an area from different glacial refugia.

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.

Repeated evolution of chimeric fusion genes in the β-globin gene family of laurasiatherian mammals

The evolutionary fate of chimeric fusion genes may be strongly influenced by their recombinational mode of origin and the nature of functional divergence between the parental genes. In the β-globin gene family of placental mammals, the two postnatally expressed δ- and β-globin genes (HBD and HBB, respectively) have a propensity for recombinational exchange via gene conversion and unequal crossing-over. In the latter case, there are good reasons to expect differences in retention rates for the reciprocal HBB/HBD and HBD/HBB fusion genes due to thalassemia pathologies associated with the HBD/HBB "Lepore" deletion mutant in humans. Here, we report a comparative genomic analysis of the mammalian β-globin gene cluster, which revealed that chimeric HBB/HBD fusion genes originated independently in four separate lineages of laurasiatherian mammals: Eulipotyphlans (shrews, moles, and hedgehogs), carnivores, microchiropteran bats, and cetaceans. In cases where an independently derived "anti-Lepore" duplication mutant has become fixed, the parental HBD and/or HBB genes have typically been inactivated or deleted, so that the newly created HBB/HBD fusion gene is primarily responsible for synthesizing the β-type subunits of adult and fetal hemoglobin (Hb). Contrary to conventional wisdom that the HBD gene is a vestigial relict that is typically inactivated or expressed at negligible levels, we show that HBD-like genes often encode a substantial fraction (20-100%) of β-chain Hbs in laurasiatherian taxa. Our results indicate that the ascendancy or resuscitation of genes with HBD-like coding sequence requires the secondary acquisition of HBB-like promoter sequence via unequal crossing-over or interparalog gene conversion.

The molecular evolution of cytochrome P450 genes within and between drosophila species

We map 114 gene gains and 74 gene losses in the P450 gene family across the phylogeny of 12 Drosophila species by examining the congruence of gene trees and species trees. Although the number of P450 genes varies from 74 to 94 in the species examined, we infer that there were at least 77 P450 genes in the ancestral Drosophila genome. One of the most striking observations in the data set is the elevated loss of P450 genes in the Drosophila sechellia lineage. The gain and loss events are not evenly distributed among the P450 genes-with 30 genes showing no gene gains or losses whereas others show as many as 20 copy number changes among the species examined. The P450 gene clades showing the fewest number of gene gain and loss events tend to be those evolving with the most purifying selection acting on the protein sequences, although there are exceptions, such as the rapid rate of amino acid replacement observed in the single copy phantom (Cyp306a1) gene. Within D. melanogaster, we observe gene copy number polymorphism in ten P450 genes including multiple cases of interparalog chimeras. Nonallelic homologous recombination (NAHR) has been associated with deleterious mutations in humans, but here we provide a second possible example of an NAHR event in insect P450s being adaptive. Specifically, we find that a polymorphic Cyp12a4/Cyp12a5 chimera correlates with resistance to an insecticide. Although we observe such interparalog exchange in our within-species data sets, we have little evidence of it between species, raising the possibility that such events may occur more frequently than appreciated but are masked by subsequent sequence change.

Deer mouse hemoglobin exhibits a lowered oxygen affinity owing to mobility of the E helix

The deer mouse, Peromyscus maniculatus, exhibits altitude-associated variation in hemoglobin oxygen affinity. To examine the structural basis of this functional variation, the structure of the hemoglobin was solved. Recombinant hemoglobin was expressed in Escherichia coli and was purified by ion-exchange chromatography. Recombinant hemoglobin was crystallized by the hanging-drop vapor-diffusion method using polyethylene glycol as a precipitant. The obtained orthorhombic crystal contained two subunits in the asymmetric unit. The refined structure was interpreted as the aquo-met form. Structural comparisons were performed among hemoglobins from deer mouse, house mouse and human. In contrast to human hemoglobin, deer mouse hemoglobin lacks the hydrogen bond between α1Trp14 in the A helix and α1Thr67 in the E helix owing to the Thr67Ala substitution. In addition, deer mouse hemoglobin has a unique hydrogen bond at the α1β1 interface between residues α1Cys34 and β1Ser128.

Synchronous evolution of an odor biosynthesis pathway and behavioral response

BACKGROUND

Rodents use olfactory cues for species-specific behaviors. For example, mice emit odors to attract mates of the same species, but not competitors of closely related species. This implies rapid evolution of olfactory signaling, although odors and chemosensory receptors involved are unknown.

RESULTS

Here, we identify a mouse chemosignal, trimethylamine, and its olfactory receptor, trace amine-associated receptor 5 (TAAR5), to be involved in species-specific social communication. Abundant (>1,000-fold increased) and sex-dependent trimethylamine production arose de novo along the Mus lineage after divergence from Mus caroli. The two-step trimethylamine biosynthesis pathway involves synergy between commensal microflora and a sex-dependent liver enzyme, flavin-containing monooxygenase 3 (FMO3), which oxidizes trimethylamine. One key evolutionary alteration in this pathway is the recent acquisition in Mus of male-specific Fmo3 gene repression. Coincident with its evolving biosynthesis, trimethylamine evokes species-specific behaviors, attracting mice, but repelling rats. Attraction to trimethylamine is abolished in TAAR5 knockout mice, and furthermore, attraction to mouse scent is impaired by enzymatic depletion of trimethylamine or TAAR5 knockout.

CONCLUSIONS

TAAR5 is an evolutionarily conserved olfactory receptor required for a species-specific behavior. Synchronized changes in odor biosynthesis pathways and odor-evoked behaviors could ensure species-appropriate social interactions.

Whole-genome duplication and the functional diversification of teleost fish hemoglobins

Subsequent to the two rounds of whole-genome duplication that occurred in the common ancestor of vertebrates, a third genome duplication occurred in the stem lineage of teleost fishes. This teleost-specific genome duplication (TGD) is thought to have provided genetic raw materials for the physiological, morphological, and behavioral diversification of this highly speciose group. The extreme physiological versatility of teleost fish is manifest in their diversity of blood–gas transport traits, which reflects the myriad solutions that have evolved to maintain tissue O₂ delivery in the face of changing metabolic demands and environmental O₂ availability during different ontogenetic stages. During the course of development, regulatory changes in blood–O₂ transport are mediated by the expression of multiple, functionally distinct hemoglobin (Hb) isoforms that meet the particular O₂-transport challenges encountered by the developing embryo or fetus (in viviparous or oviparous species) and in free-swimming larvae and adults. The main objective of the present study was to assess the relative contributions of whole-genome duplication, large-scale segmental duplication, and small-scale gene duplication in producing the extraordinary functional diversity of teleost Hbs. To accomplish this, we integrated phylogenetic reconstructions with analyses of conserved synteny to characterize the genomic organization and evolutionary history of the globin gene clusters of teleosts. These results were then integrated with available experimental data on functional properties and developmental patterns of stage-specific gene expression. Our results indicate that multiple α- and β-globin genes were present in the common ancestor of gars (order Lepisoteiformes) and teleosts. The comparative genomic analysis revealed that teleosts possess a dual set of TGD-derived globin gene clusters, each of which has undergone lineage-specific changes in gene content via repeated duplication and deletion events. Phylogenetic reconstructions revealed that paralogous genes convergently evolved similar functional properties in different teleost lineages. Consistent with other recent studies of globin gene family evolution in vertebrates, our results revealed evidence for repeated evolutionary transitions in the developmental regulation of Hb synthesis.

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.

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.

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.

A tale of two genomes: contrasting patterns of phylogeographic structure in a widely distributed bat

One of the most widely distributed bats in the New World, the big brown bat (Eptesicus fuscus) exhibits well-documented geographic variation in morphology and life history traits, suggesting the potential for significant phylogeographic structure as well as adaptive differentiation among populations. In a pattern broadly consistent with morphologically defined subspecies, we found deeply divergent mitochondrial lineages restricted to different geographic regions. In contrast, sequence data from two nuclear loci suggest a general lack of regional genetic structure except for peripheral populations in the Caribbean and Mexico/South America. Coalescent analyses suggest that the striking difference in population structure between genomes cannot be attributed solely to different rates of lineage sorting, but is likely due to male-mediated gene flow homogenizing nuclear genetic diversity across most of the continental range. Despite this ongoing gene flow, selection has apparently been effective in producing and maintaining adaptive differentiation among populations, while strong female site fidelity, maintained over the course of millions of years, has produced remarkably deep divergence among geographically isolated matrilines. Our results highlight the importance of evaluating multiple genetic markers for a more complete understanding of population structure and history.

Evolutionary and functional properties of a two-locus beta-globin polymorphism in Indian house mice

Electrophoretic surveys of hemoglobin (Hb) polymorphism in house mice from South Asia and the Middle East have revealed that two alternative beta-globin haplotypes, Hbb(d) and Hbb(p), are often present at intermediate frequencies in geographically disparate populations. Both haplotypes harbor two functionally distinct beta-globin paralogs, HBB-T1 (which encodes the beta-chain subunits of the major Hb isoform) and HBB-T2 (which encodes the beta-chains of the minor Hb isoform). The Hbb(d) and Hbb(p) haplotypes share identical HBB-T1 alleles, but products of the alternative HBB-T2 alleles (d(minor) and p(minor)) are distinguished by two amino acid substitutions. To investigate the possible adaptive significance of the Hbb(d)/Hbb(p) polymorphism we conducted a population genetic analysis of the duplicated beta-globin genes of Indian house mice (Mus castaneus) in conjunction with experimental studies of Hb function in inbred strains of mice that carry the alternative Hbb(d) and Hbb(p) haplotypes. The main objectives of this study were (i) to characterize patterns of nucleotide polymorphism and linkage disequilibrium in the duplicated beta-globin genes of M. castaneus, (ii) to test the hypothesis that the Hbb(d) and Hbb(p) haplotypes are maintained as a balanced polymorphism, and (iii) to assess whether allelic differences in the alternative minor Hb isoforms (d(minor) and p(minor)) are associated with different O(2)-binding properties. A multilocus analysis of polymorphism and divergence revealed that levels of diversity at the HBB-T2 gene exceeded neutral expectations, and reconstructed haplotype networks for both beta-globin paralogs revealed extensive allele sharing with several other closely related species of Mus. However, despite this suggestive evidence for balancing selection, O(2)-equilibrium curves revealed no discernible functional differences between red cell lysates containing the d(minor) and p(minor) Hb isoforms. If the d(minor) and p(minor) alleles are maintained as a balanced polymorphism, our results indicate that the associated fitness variance is not directly related to respiratory functions of Hb.

Lineage-specific patterns of functional diversification in the alpha- and beta-globin gene families of tetrapod vertebrates

The alpha- and beta-globin gene families of jawed vertebrates have diversified with respect to both gene function and the developmental timing of gene expression. Phylogenetic reconstructions of globin gene family evolution have provided suggestive evidence that the developmental regulation of hemoglobin synthesis has evolved independently in multiple vertebrate lineages. For example, the embryonic beta-like globin genes of birds and placental mammals are not 1:1 orthologs. Despite the similarity in developmental expression profiles, the genes are independently derived from lineage-specific duplications of a beta-globin pro-ortholog. This suggests the possibility that other vertebrate taxa may also possess distinct repertoires of globin genes that were produced by repeated rounds of lineage-specific gene duplication and divergence. Until recently, investigations into this possibility have been hindered by the dearth of genomic sequence data from nonmammalian vertebrates. Here, we report new insights into globin gene family evolution that were provided by a phylogenetic analysis of vertebrate globins combined with a comparative genomic analysis of three key sauropsid taxa: a squamate reptile (anole lizard, Anolis carolinensis), a passeriform bird (zebra finch, Taeniopygia guttata), and a galliform bird (chicken, Gallus gallus). The main objectives of this study were 1) to characterize evolutionary changes in the size and membership composition of the alpha- and beta-globin gene families of tetrapod vertebrates and 2) to test whether functional diversification of the globin gene clusters occurred independently in different tetrapod lineages. Results of our comparative genomic analysis revealed several intriguing patterns of gene turnover in the globin gene clusters of different taxa. Lineage-specific differences in gene content were especially pronounced in the beta-globin gene family, as phylogenetic reconstructions revealed that amphibians, lepidosaurs (as represented by anole lizard), archosaurs (as represented by zebra finch and chicken), and mammals each possess a distinct independently derived repertoire of beta-like globin genes. In contrast to the ancient functional diversification of the alpha-globin gene cluster in the stem lineage of tetrapods, the physiological division of labor between early- and late-expressed genes in the beta-globin gene cluster appears to have evolved independently in several tetrapod lineages.