Complex signatures of selection and gene conversion in the duplicated globin genes of house mice

House Mice in the Atlantic Region: Genetic Signals of Their Human Transport

BACKGROUND/OBJECTIVES

The colonization history of house mice reflects the maritime history of humans that passively transported them worldwide. We investigated western house mouse colonization in the Atlantic region through studies of mitochondrial D-loop DNA sequences from modern specimens.

METHODS

We assembled a dataset of 758 haplotypes derived from 2765 mice from 47 countries/oceanic archipelagos (a combination of new and published data). Our maximum likelihood phylogeny recovered five previously identified clades, and we used the haplotype affinities within the phylogeny to infer house mouse colonization history, employing statistical tests and indices. From human history, we predefined four European source areas for mice in the Atlantic region (Northern Europe excluding Scandinavia, Southern Europe, Scandinavia, and Macaronesia) and we investigated the colonization from these source areas to different geographic areas in the Atlantic region.

RESULTS

Our inferences suggest mouse colonization of Scandinavia itself from Northern Europe, and Macaronesia from both Southern Europe and Scandinavia/Germany (the latter likely representing the transport of mice by Vikings). Mice on North Atlantic islands apparently derive primarily from Scandinavia, while for South Atlantic islands, North America, and Sub-Saharan Africa, the clearest source is Northern Europe, although mice on South Atlantic islands also had genetic inputs from Macaronesia and Southern Europe (for Tristan da Cunha). Macaronesia was a stopover for Atlantic voyages, creating an opportunity for mouse infestation. Mice in Latin America also apparently had multiple colonization sources, with a strong Southern European signal but also input from Northern Europe and/or Macaronesia.

CONCLUSIONS

D-loop sequences help discern the broad-scale colonization history of house mice and new perspectives on human history.

Varying strength of selection contributes to the intragenomic diversity of rRNA genes

Ribosome biogenesis in eukaryotes is supported by hundreds of ribosomal RNA (rRNA) gene copies that are encoded in the ribosomal DNA (rDNA). The multiple copies of rRNA genes are thought to have low sequence diversity within one species. Here, we present species-wide rDNA sequence analysis in Saccharomyces cerevisiae that challenges this view. We show that rDNA copies in this yeast are heterogeneous, both among and within isolates, and that many variants avoided fixation or elimination over evolutionary time. The sequence diversity landscape across the rDNA shows clear functional stratification, suggesting different copy-number thresholds for selection that contribute to rDNA diversity. Notably, nucleotide variants in the most conserved rDNA regions are sufficiently deleterious to exhibit signatures of purifying selection even when present in only a small fraction of rRNA gene copies. Our results portray a complex evolutionary landscape that shapes rDNA sequence diversity within a single species and reveal unexpectedly strong purifying selection of multi-copy genes.

The Bank Vole (Clethrionomys glareolus) as a Model System for Adaptive Phylogeography in the European Theater

The legacy of climatic changes during the Pleistocene glaciations allows inferences to be made about the patterns and processes associated with range expansion/colonization, including evolutionary adaptation. With the increasing availability of population genomic data, we have the opportunity to examine these questions in detail and in a variety of non-traditional model species. As an exemplar, here we review more than two decades of work by our group and others that illustrate the potential of a single “non-model model” mammal species - the bank vole (Clethrionomys glareolus), which is particularly well suited to illustrate the complexities that may be associated with range expansion and the power of genomics (and other datasets) to uncover them. We first summarize early phylogeographic work using mitochondrial DNA and then describe new phylogeographic insights gained from population genomic analysis of genome-wide SNP data to highlight the bank vole as one of the most compelling examples of a forest mammal, that survived in cryptic extra-Mediterranean (“northern”) glacial refugia in Europe, and as one of the species in which substantial replacement and mixing of lineages originating from different refugia occurred during end-glacial colonization. Our studies of bank vole hemoglobin structure and function, as well as our recent ecological niche modeling study examining differences among bank vole lineages, led us to develop the idea of “adaptive phylogeography.” This is what we call the study of the role of adaptive differences among populations in shaping phylogeographic patterns. Adaptive phylogeography provides a link between past population history and adaptation that can ultimately help predict the potential of future species responses to climate change. Because the bank vole is part of a community of organisms whose range has repeatedly contracted and then expanded in the past, what we learn from the bank vole will be useful for our understanding of a broad range of species.

The genomic basis of high-elevation adaptation in wild house mice (Mus musculus domesticus) from South America

Understanding the genetic basis of environmental adaptation in natural populations is a central goal in evolutionary biology. The conditions at high elevation, particularly the low oxygen available in the ambient air, impose a significant and chronic environmental challenge to metabolically active animals with lowland ancestry. To understand the process of adaptation to these novel conditions and to assess the repeatability of evolution over short timescales, we examined the signature of selection from complete exome sequences of house mice (Mus musculus domesticus) sampled across two elevational transects in the Andes of South America. Using phylogenetic analysis, we show that house mice colonized high elevations independently in Ecuador and Bolivia. Overall, we found distinct responses to selection in each transect and largely nonoverlapping sets of candidate genes, consistent with the complex nature of traits that underlie adaptation to low oxygen availability (hypoxia) in other species. Nonetheless, we also identified a small subset of the genome that appears to be under parallel selection at the gene and SNP levels. In particular, three genes (Col22a1, Fgf14, and srGAP1) bore strong signatures of selection in both transects. Finally, we observed several patterns that were common to both transects, including an excess of derived alleles at high elevation, and a number of hypoxia-associated genes exhibiting a threshold effect, with a large allele frequency change only at the highest elevations. This threshold effect suggests that selection pressures may increase disproportionately at high elevations in mammals, consistent with observations of some high-elevation diseases in humans.

High-altitude rodents have abundant collaterals that protect against tissue injury after cerebral, coronary and peripheral artery occlusion

Collateral number/density varies widely in brain and other tissues among strains of Mus musculus mice due to differences in genetic background. Recent studies have shown that prolonged exposure to reduced atmospheric oxygen induces additional collaterals to form, suggesting that natural selection may favor increased collaterals in populations native to high-altitude. High-altitude guinea pigs (Cavia) and deer mice (Peromyscus) were compared with lowland species of Peromyscus, Mus and Rattus (9 species/strains examined). Collateral density, diameter and other morphometrics were measured in brain where, importantly, collateral abundance reflects that in other tissues of the same individual. Guinea pigs and high-altitude deer mice had a greater density of pial collaterals than lowlanders. Consistent with this, guinea pigs and highlander mice evidenced complete and 80% protection against stroke, respectively. They also sustained significantly less ischemia in heart and lower extremities after arterial occlusion. Vessels of the circle of Willis, including the communicating collateral arteries, also exhibited unique features in the highland species. Our findings support the hypothesis that species native to high-altitude have undergone genetic selection for abundant collaterals, suggesting that besides providing protection in obstructive disease, collaterals serve a physiological function to optimize oxygen delivery to meet oxygen demand when oxygen is limiting.

Mice and Rats Display Different Ventilatory, Hematological, and Metabolic Features of Acclimatization to Hypoxia

Phylogeographic studies showed that house mice (Mus musculus) originated in the Himalayan region, while common rats (Rattus rattus and Rattus norvegicus) come from the lowlands of China and India. Accordingly, it has been proposed that its origins gave mice, but not rats, the ability to invade ecological niches at high altitudes (pre-adaptation). This proposal is strongly supported by the fact that house mice are distributed throughout the world, while common rats are practically absent above 2,500 m. Considering that the ability of mammals to colonize high-altitude environments (>2,500 m) is limited by their capability to tolerate reduced oxygen availability, in this work, we hypothesize that divergences in the ventilatory, hematological, and metabolic phenotypes of mice and rats establish during the process of acclimatization to hypoxia (Hx). To test this hypothesis male FVB mice and Sprague-Dawley (SD) rats were exposed to Hx (12% O₂) for 0 h (normoxic controls), 6 h, 1, 7, and 21 days. We assessed changes in ventilatory [minute ventilation (V_E), respiratory frequency (f_R), and tidal volume (V_T)], hematological (hematocrit and hemoglobin concentration), and metabolic [whole-body O₂ consumption (VO₂) and CO₂ production (VCO₂), and liver mitochondrial oxygen consumption rate (OCR) parameters]. Compared to rats, results in mice show increased ventilatory, metabolic, and mitochondrial response. In contrast, rats showed quicker and higher hematological response than mice and only minor ventilatory and metabolic adjustments. Our findings may explain, at least in part, why mice, but not rats, were able to colonize high-altitude habitats.

Balancing selection in Pattern Recognition Receptor signalling pathways is associated with gene function and pleiotropy in a wild rodent

Pathogen-mediated balancing selection is commonly considered to play an important role in the maintenance of genetic diversity, in particular in immune genes. However, the factors that may influence which immune genes are the targets of such selection are largely unknown. To address this, here we focus on Pattern Recognition Receptor (PRR) signalling pathways, which play a key role in innate immunity. We used whole-genome resequencing data from a population of bank voles (Myodes glareolus) to test for associations between balancing selection, pleiotropy and gene function in a set of 123 PRR signalling pathway genes. To investigate the effect of gene function, we compared genes encoding (a) receptors for microbial ligands versus downstream signalling proteins, and (b) receptors recognizing components of microbial cell walls, flagella and capsids versus receptors recognizing features of microbial nucleic acids. Analyses based on the nucleotide diversity of full coding sequences showed that balancing selection primarily targeted receptor genes with a low degree of pleiotropy. Moreover, genes encoding receptors recognizing components of microbial cell walls etc. were more important targets of balancing selection than receptors recognizing nucleic acids. Tests for localized signatures of balancing selection in coding and noncoding sequences showed that such signatures were mostly located in introns, and more evenly distributed among different functional categories of PRR pathway genes. The finding that signatures of balancing selection in full coding sequences primarily occur in receptor genes, in particular those encoding receptors for components of microbial cell walls etc., is consistent with the idea that coevolution between hosts and pathogens is an important cause of balancing selection on immune genes.

Emergence of a Chimeric Globin Pseudogene and Increased Hemoglobin Oxygen Affinity Underlie the Evolution of Aquatic Specializations in Sirenia

As limits on O₂ availability during submergence impose severe constraints on aerobic respiration, the oxygen binding globin proteins of marine mammals are expected to have evolved under strong evolutionary pressures during their land-to-sea transition. Here, we address this question for the order Sirenia by retrieving, annotating, and performing detailed selection analyses on the globin repertoire of the extinct Steller’s sea cow (Hydrodamalis gigas), dugong (Dugong dugon), and Florida manatee (Trichechus manatus latirostris) in relation to their closest living terrestrial relatives (elephants and hyraxes). These analyses indicate most loci experienced elevated nucleotide substitution rates during their transition to a fully aquatic lifestyle. While most of these genes evolved under neutrality or strong purifying selection, the rate of nonsynonymous/synonymous replacements increased in two genes (Hbz-T1 and Hba-T1) that encode the α-type chains of hemoglobin (Hb) during each stage of life. Notably, the relaxed evolution of Hba-T1 is temporally coupled with the emergence of a chimeric pseudogene (Hba-T2/Hbq-ps) that contributed to the tandemly linked Hba-T1 of stem sirenians via interparalog gene conversion. Functional tests on recombinant Hb proteins from extant and ancestral sirenians further revealed that the molecular remodeling of Hba-T1 coincided with increased Hb–O₂ affinity in early sirenians. Available evidence suggests that this trait evolved to maximize O₂ extraction from finite lung stores and suppress tissue O₂ offloading, thereby facilitating the low metabolic intensities of extant sirenians. In contrast, the derived reduction in Hb–O₂ affinity in (sub)Arctic Steller’s sea cows is consistent with fueling increased thermogenesis by these once colossal marine herbivores.

Contrasting patterns of coding and flanking region evolution in mammalian keratin associated protein-1 genes

Mammalian genomes contain a number of duplicated genes, and sequence identity between these duplicates can be maintained by purifying selection. However, between-duplicate recombination can also maintain sequence identity between copies, resulting in a pattern known as concerted evolution where within-genome repeats are more similar to each other than to orthologous repeats in related species. Here we investigated the tandemly-repeated keratin-associated protein 1 (KAP1) gene family, KRTAP1, which encodes proteins that are important components of hair and wool in mammals. Comparison of eutherian mammal KRTAP1 gene repeats within and between species shows a strong pattern of concerted evolution. However, in striking contrast to the coding regions of these genes, we find that the flanking regions have a divergent pattern of evolution. This contrast in evolutionary pattern transitions abruptly near the start and stop codons of the KRTAP1 genes. We reveal that this difference in evolutionary patterns is not explained by conventional purifying selection, nor is it likely a consequence of codon adaptation or reverse transcription of KRTAP1-n mRNA. Instead, the evidence suggests that these contrasting patterns result from short-tract gene conversion events that are biased to the KRTAP1 coding region by selection and/or differential sequence divergence. This work demonstrates the power that gene conversion has to finely shape the evolution of repetitive genes, and provides another distinctive pattern of contrasting evolutionary outcomes that results from gene conversion. A greater emphasis on exploring the evolution of multi-gene eukaryotic families will reveal how common different contrasting evolutionary patterns are in gene duplicates.

Altitudinal variation of the gut microbiota in wild house mice

The maintenance of oxygen homeostasis in the gut is critical for the maintenance of a healthy gut microbiota. However, few studies have explored how the concentration of atmospheric oxygen affects the gut microbiota in natural populations. High-altitude environments provide an opportunity to study the potential effects of atmospheric oxygen on the composition and function of the gut microbiota. Here, we characterized the caecal microbial communities of wild house mice (Mus musculus domesticus) in two independent altitudinal transects, one in Ecuador and one in Bolivia, from sea level to nearly 4,000 m. First, we found that differences in altitude were associated with differences in the gut microbial community after controlling for the effects of body mass, diet, reproductive status and population of origin. Second, obligate anaerobes tended to show a positive correlation with altitude, while all other microbes tended to show a negative correlation with altitude. These patterns were seen independently in both transects, consistent with the expected effects of atmospheric oxygen on gut microbes. Prevotella was the most-enriched genus at high elevations in both transects, consistent with observations in high-altitude populations of pikas, ruminants and humans, and also consistent with observations of laboratory mice exposed to hypoxic conditions. Lastly, the renin-angiotensin system, a recently proposed microbiota-mediated pathway of blood pressure regulation, was the top predicted metagenomic pathway enriched in high altitudes in both transects. These results suggest that high-altitude environments affect the composition and function of the gut microbiota in wild mammals.

Playing Hide-and-Seek in Beta-Globin Genes: Gene Conversion Transferring a Beneficial Mutation between Differentially Expressed Gene Duplicates

Increasing evidence suggests that adaptation to diverse environments often involves selection on existing variation rather than new mutations. A previous study identified a nonsynonymous single nucleotide polymorphism (SNP) in exon 2 of two paralogous β-globin genes of the bank vole (Clethrionomys glareolus) in Britain in which the ancestral serine (Ser) and the derived cysteine (Cys) allele represent geographically partitioned functional variation affecting the erythrocyte antioxidative capacity. Here we studied the geographical pattern of the two-locus Ser/Cys polymorphism throughout Europe and tested for the geographic correlation between environmental variables and allele frequency, expected if the polymorphism was under spatially heterogeneous environment-related selection. Although bank vole population history clearly is important in shaping the dispersal of the oxidative stress protective Cys allele, analyses correcting for population structure suggest the Europe-wide pattern is affected by geographical variation in environmental conditions. The β-globin phenotype is encoded by the major paralog HBB-T1 but we found evidence of bidirectional gene conversion of exon 2 with the low-expression paralog HBB-T2. Our data support the model where gene conversion reshuffling genotypes between high- and low- expressed paralogs enables tuning of erythrocyte thiol levels, which may help maintain intracellular redox balance under fluctuating environmental conditions. Therefore, our study suggests a possible role for gene conversion between differentially expressed gene duplicates as a mechanism of physiological adaptation of populations to new or changing environments.

Divergent Mitochondrial Antioxidant Activities and Lung Alveolar Architecture in the Lungs of Rats and Mice at High Altitude

Compared with mice, adult rats living at 3,600 m above sea level (SL—La Paz, Bolivia) have high hematocrit, signs of pulmonary hypertension, and low lung volume with reduced alveolar surface area. This phenotype is associated with chronic mountain sickness in humans living at high altitude (HA). We tested the hypothesis that this phenotype is associated with impaired gas exchange and oxidative stress in the lungs. We used rats and mice (3 months old) living at HA (La Paz) and SL (Quebec City, Canada) to measure arterial oxygen saturation under graded levels of hypoxia (by pulse oximetry), the alveolar surface area in lung slices and the activity of pro- (NADPH and xanthine oxidases—NOX and XO) and anti- (superoxide dismutase, and glutathione peroxidase—SOD and GPx) oxidant enzymes in cytosolic and mitochondrial lung protein extracts. HA rats have a lower arterial oxygen saturation and reduced alveolar surface area compared to HA mice and SL rats. Enzymatic activities (NOX, XO, SOD, and GPx) in the cytosol were similar between HA and SL animals, but SOD and GPx activities in the mitochondria were 2–3 times higher in HA vs. SL rats, and only marginally higher in HA mice vs. SL mice. Furthermore, the maximum activity of cytochrome oxidase-c (COX) measured in mitochondrial lung extracts was also 2 times higher in HA rats compared with SL rats, while there was only a small increase in HA mice vs. SL mice. Interestingly, compared with SL controls, alterations in lung morphology are not observed for young rats at HA (15 days after birth), and enzymatic activities are only slightly altered. These results suggest that rats living at HA have a gradual reduction of their alveolar surface area beyond the postnatal period. We can speculate that the elevation of SOD, GPx, and COX activities in the lung mitochondria are not sufficient to compensate for oxidative stress, leading to damage of the lung tissue in rats.

Panic-like escape response elicited in mice by exposure to CO2, but not hypoxia

Exposure to elevated concentrations of CO₂ or hypoxia has been widely used in psychiatric research as a panic provoking stimulus. However, the use of these respiratory challenges to model panic-like responses in experimental animals has been less straightforward. Little data is available, from behavioral and endocrine perspectives, to support the conclusion that a marked aversive situation, such as that experienced during panic attacks, was evoked in these animals. We here compared the behavioral responses of male CB57BL/6 mice during exposure to 20% CO₂ or 7% O₂ and its consequence on plasma levels of corticosterone. We also evaluated whether clinically-effective panicolytic drugs affect the behavioral responses expressed during CO₂ exposure. The results showed that whereas hypoxia caused a marked reduction in locomotion, inhalation of CO₂-enriched air evoked an active escape response, characterized by bouts of upward leaps directed to the border of the experimental cage, interpreted as escape attempts. Corticosterone levels were increased 30min after either of the respiratory challenges used, but it was higher in the hypoxia group. Chronic (21days), but not acute, treatment with fluoxetine or imipramine (5, 10 or 15mg/kg) or a single injection of alprazolam (0.025, 0.05 or 0.1mg/kg), but not of the anxiolytic diazepam (0.025, 0.05 or 0.1 and 1mg/kg) reduced the number of escape attempts, indicating a panicolytic-like effect. Altogether, the results suggest that whereas hypoxia increased anxiety, exposure to 20% CO₂ evoked a panic-like state. The latter condition/test protocol seems to be a simple and validated model for studying in mice pathophysiological mechanisms and the screening of novel drugs for panic disorder.

Frequent nonallelic gene conversion on the human lineage and its effect on the divergence of gene duplicates

Gene conversion is the copying of a genetic sequence from a "donor" region to an "acceptor." In nonallelic gene conversion (NAGC), the donor and the acceptor are at distinct genetic loci. Despite the role NAGC plays in various genetic diseases and the concerted evolution of gene families, the parameters that govern NAGC are not well characterized. Here, we survey duplicate gene families and identify converted tracts in 46% of them. These conversions reflect a large GC bias of NAGC. We develop a sequence evolution model that leverages substantially more information in duplicate sequences than used by previous methods and use it to estimate the parameters that govern NAGC in humans: a mean converted tract length of 250 bp and a probability of [Formula: see text] per generation for a nucleotide to be converted (an order of magnitude higher than the point mutation rate). Despite this high baseline rate, we show that NAGC slows down as duplicate sequences diverge-until an eventual "escape" of the sequences from its influence. As a result, NAGC has a small average effect on the sequence divergence of duplicates. This work improves our understanding of the NAGC mechanism and the role that it plays in the evolution of gene duplicates.

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.

HIF1α and physiological responses to hypoxia are correlated in mice but not in rats

We previously reported that rats and mice that have been raised for more than 30 generations in La Paz, Bolivia (3600m), display divergent physiological responses to high altitude (HA), including improved respiratory and metabolic control in mice. In the present study we asked whether these traits would also be present in response to hypoxia at sea level (SL). To answer this question, we exposed rats (SD) and mice (FVB) to normoxia (21% O2) or hypoxia (15 and 12% O2) for 6 hours and measured ventilation and metabolic rate (whole body plethysmography), and expression of the transcription factor HIF-1α (ELISA and Mass Spectrometry) and other proteins whose expression are regulated by hypoxia (Glucose Transporter 1, Pyruvate Dehydrogenase Kinase 1, and Angiopoietin 2 - Mass Spectrometry) in the brainstem. In response to hypoxia, compared with rats, mice had higher minute ventilation, lower metabolic rate, and higher expression of HIF-1α in the brainstem. In mice the expression level of HIF-1α was positively correlated with ventilation and negatively correlated with metabolic rate. In rats, the concentration of brainstem cytosolic protein decreased by 38% at 12% O2, while expression of the glucose transporter 1 increased. We conclude that mice and rats raised at sea level have divergent physiological and molecular responses to hypoxia, supporting the hypothesis that mice have innate traits that favor adaptation to altitude.

Genetically based low oxygen affinities of felid hemoglobins: lack of biochemical adaptation to high-altitude hypoxia in the snow leopard

Genetically based modifications of hemoglobin (Hb) function that increase blood-O2 affinity are hallmarks of hypoxia adaptation in vertebrates. Among mammals, felid Hbs are unusual in that they have low intrinsic O2 affinities and reduced sensitivities to the allosteric cofactor 2,3-diphosphoglycerate (DPG). This combination of features compromises the acclimatization capacity of blood-O2 affinity and has led to the hypothesis that felids have a restricted physiological niche breadth relative to other mammals. In seeming defiance of this conjecture, the snow leopard (Panthera uncia) has an extraordinarily broad elevational distribution and occurs at elevations above 6000 m in the Himalayas. Here, we characterized structural and functional variation of big cat Hbs and investigated molecular mechanisms of Hb adaptation and allosteric regulation that may contribute to the extreme hypoxia tolerance of the snow leopard. Experiments revealed that purified Hbs from snow leopard and African lion exhibited equally low O2 affinities and DPG sensitivities. Both properties are primarily attributable to a single amino acid substitution, β2His→Phe, which occurred in the common ancestor of Felidae. Given the low O2 affinity and reduced regulatory capacity of feline Hbs, the extreme hypoxia tolerance of snow leopards must be attributable to compensatory modifications of other steps in the O2-transport pathway.

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.

Divergent physiological responses in laboratory rats and mice raised at high altitude

Ecological studies show that mice can be found at high altitude (HA – up to 4000 m) while rats are absent at these altitudes, and there are no data to explain this discrepancy. We used adult laboratory rats and mice that have been raised for more than 30 generations in La Paz, Bolivia (3600 m), and compared their hematocrit levels, right ventricular hypertrophy (index of pulmonary hypertension) and alveolar surface area in the lungs. We also used whole-body plethysmography, indirect calorimetry and pulse oxymetry to measure ventilation, metabolic rate (O2 consumption and CO2 production), heart rate and pulse oxymetry oxygen saturation (pO2,sat) under ambient conditions, and in response to exposure to sea level PO2 (32% O2=160 mmHg for 10 min) and hypoxia (18% and 15% O2=90 and 75 mmHg for 10 min each). The variables used for comparisons between species were corrected for body mass using standard allometric equations, and are termed mass-corrected variables. Under baseline, compared with rats, adult mice had similar levels of pO2,sat, but lower hematocrit and hemoglobin levels, reduced right ventricular hypertrophy and higher mass-corrected alveolar surface area, tidal volume and metabolic rate. In response to sea level PO2 and hypoxia, mice and rats had similar changes of ventilation, but metabolic rate decreased much more in hypoxia in mice, while pO2,sat remained higher in mice. We conclude that laboratory mice and rats that have been raised at HA for >30 generations have different physiological responses to altitude. These differences might explain the different altitude distribution observed in wild rats and mice.

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.

Genomic organization and differential signature of positive selection in the alpha and beta globin gene clusters in two cetacean species

The hemoglobin of jawed vertebrates is a heterotetramer protein that contains two α- and two β-chains, which are encoded by members of α- and β-globin gene families. Given the hemoglobin role in mediating an adaptive response to chronic hypoxia, it is likely that this molecule may have experienced a selective pressure during the evolution of cetaceans, which have to deal with hypoxia tolerance during prolonged diving. This selective pressure could have generated a complex history of gene turnover in these clusters and/or changes in protein structure themselves. Accordingly, we aimed to characterize the genomic organization of α- and β-globin gene clusters in two cetacean species and to detect a possible role of positive selection on them using a phylogenetic framework. Maximum likelihood and Bayesian phylogeny reconstructions revealed that both cetacean species had retained a similar complement of putatively functional genes. For the α-globin gene cluster, the killer whale presents a complement of genes composed of HBZ, HBK, and two functional copies of HBA and HBQ genes, whereas the dolphin possesses HBZ, HBK, HBA and HBQ genes, and one HBA pseudogene. For the β-globin gene cluster, both species retained a complement of four genes, two early expressed genes—HBE and HBH—and two adult expressed genes—HBD and HBB. Our natural selection analysis detected two positively selected sites in the HBB gene (56 and 62) and four in HBA (15, 21, 49, 120). Interestingly, only the genes that are expressed during the adulthood showed the signature of positive selection.

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.

Relaxed functional constraints on triplicate α-globin gene in the bank vole suggest a different evolutionary history from other rodents

Gene duplication plays an important role in the origin of evolutionary novelties, but the mechanisms responsible for the retention and functional divergence of the duplicated copy are not fully understood. The α-globin genes provide an example of a gene family with different numbers of gene duplicates among rodents. Whereas Rattus and Peromyscus each have three adult α-globin genes (HBA-T1, HBA-T2 and HBA-T3), Mus has only two copies. High rates of amino acid evolution in the independently derived HBA-T3 genes of Peromyscus and Rattus have been attributed to positive selection. Using RACE PCR, reverse transcription-PCR (RT-PCR) and RNA-seq, we show that another rodent, the bank vole Clethrionomys glareolus, possesses three transcriptionally active α-globin genes. The bank vole HBA-T3 gene is distinguished from each HBA-T1 and HBA-T2 by 20 amino acids and is transcribed 23- and 4-fold lower than HBA-T1 and HBA-T2, respectively. Polypeptides corresponding to all three genes are detected by electrophoresis, demonstrating that the translated products of HBA-T3 are present in adult erythrocytes. Patterns of codon substitution and the presence of low-frequency null alleles suggest a postduplication relaxation of purifying selection on bank vole HBA-T3.

Both positive and negative selection pressures contribute to the polymorphism pattern of the duplicated human CYP21A2 gene

The human steroid 21-hydroxylase gene (CYP21A2) participates in cortisol and aldosterone biosynthesis, and resides together with its paralogous (duplicated) pseudogene in a multiallelic copy number variation (CNV), called RCCX CNV. Concerted evolution caused by non-allelic gene conversion has been described in great ape CYP21 genes, and the same conversion activity is responsible for a serious genetic disorder of CYP21A2, congenital adrenal hyperplasia (CAH). In the current study, 33 CYP21A2 haplotype variants encoding 6 protein variants were determined from a European population. CYP21A2 was shown to be one of the most diverse human genes (HHe=0.949), but the diversity of intron 2 was greater still. Contrary to previous findings, the evolution of intron 2 did not follow concerted evolution, although the remaining part of the gene did. Fixed sites (different fixed alleles of sites in human CYP21 paralogues) significantly accumulated in intron 2, indicating that the excess of fixed sites was connected to the lack of effective non-allelic conversion and concerted evolution. Furthermore, positive selection was presumably focused on intron 2, and possibly associated with the previous genetic features. However, the positive selection detected by several neutrality tests was discerned along the whole gene. In addition, the clear signature of negative selection was observed in the coding sequence. The maintenance of the CYP21 enzyme function is critical, and could lead to negative selection, whereas the presumed gene regulation altering steroid hormone levels via intron 2 might help fast adaptation, which broadly characterizes the genes of human CNVs responding to the environment.

Gene conversion yields novel gene combinations in paralogs of GOT1 in the copepod Tigriopus californicus

Background

Gene conversion of duplicated genes can slow the divergence of paralogous copies over time but can also result in other interesting evolutionary patterns. Islands of genetic divergence that persist in the face of gene conversion can point to gene regions undergoing selection for new functions. Novel combinations of genetic variation that differ greatly from the original sequence can result from the transfer of genetic variation between paralogous genes by rare gene conversion events. Genetically divergent populations of the copepod Tigriopus californicus provide an excellent model to look at the patterns of divergence among paralogs across multiple independent evolutionary lineages.

Results

In this study the evolution of a set of paralogous genes encoding putative aspartate transaminase proteins (called GOT1 here) are examined in populations of the copepod T. californicus. One pair of duplicated genes, GOT1p1 and GOT1p2, has regions of high divergence between the copies in the face of apparent on-going gene conversion. The GOT1p2 gene also has unique haplotypes in two populations that appear to have resulted from a transfer of genetic variation via inter-paralog gene conversion. A second pair of duplicated genes GOT1Sr and GOT1Sd also shows evidence of gene conversion, but this gene conversion does not appear to have maintained each as a functional copy in all populations.

Conclusions

The patterns of conservation and sequence divergence across this set of paralogous genes among populations of T. californicus suggest that some interesting evolutionary patterns are occurring at these loci. The results for the GOT1p1/GOT1p2 paralogs illustrate how gene conversion can factor in the creation of a mosaic pattern of regions of high divergence and low divergence. When coupled with rare gene conversion events of divergent regions, this pattern can result in the formation of novel proteins differing substantially from either original protein. The evolutionary patterns across these paralogs show how gene conversion can both constrain and facilitate diversification of genetic sequences.

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.

Genetic analysis of hematological parameters in incipient lines of the collaborative cross

Hematological parameters, including red and white blood cell counts and hemoglobin concentration, are widely used clinical indicators of health and disease. These traits are tightly regulated in healthy individuals and are under genetic control. Mutations in key genes that affect hematological parameters have important phenotypic consequences, including multiple variants that affect susceptibility to malarial disease. However, most variation in hematological traits is continuous and is presumably influenced by multiple loci and variants with small phenotypic effects. We used a newly developed mouse resource population, the Collaborative Cross (CC), to identify genetic determinants of hematological parameters. We surveyed the eight founder strains of the CC and performed a mapping study using 131 incipient lines of the CC. Genome scans identified quantitative trait loci for several hematological parameters, including mean red cell volume (Chr 7 and Chr 14), white blood cell count (Chr 18), percent neutrophils/lymphocytes (Chr 11), and monocyte number (Chr 1). We used evolutionary principles and unique bioinformatics resources to reduce the size of candidate intervals and to view functional variation in the context of phylogeny. Many quantitative trait loci regions could be narrowed sufficiently to identify a small number of promising candidate genes. This approach not only expands our knowledge about hematological traits but also demonstrates the unique ability of the CC to elucidate the genetic architecture of complex traits.

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.

Long-term balancing selection at the blood group-related gene B4galnt2 in the genus Mus (Rodentia; Muridae)

Recent surveys of the human genome have highlighted the significance of balancing selection in relation to understanding the evolutionary origins of disease-associated variation. Cis-regulatory variation at the blood group-related glycosyltransferase B4galnt2 is associated with a phenotype in mice that closely resembles a common human bleeding disorder, von Willebrand disease. In this study, we have performed a survey of the 5' flanking region of the B4galnt2 gene in several Mus musculus subspecies and Mus spretus. Our results reveal a clear pattern of trans-species polymorphism and indicate that allele classes conferring alternative tissue-specific expression patterns have been maintained for >2.8 My in the genus Mus. Furthermore, analysis of B4galnt2 expression patterns revealed the presence of an additional functional class of alleles, supporting a role for gastrointestinal phenotypes in the long-term maintenance of expression variation at this gene.

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.

Neutral and non-neutral evolution of duplicated genes with gene conversion

Gene conversion is one of the major mutational mechanisms involved in the DNA sequence evolution of duplicated genes. It contributes to create unique patters of DNA polymorphism within species and divergence between species. A typical pattern is so-called concerted evolution, in which the divergence between duplicates is maintained low for a long time because of frequent exchanges of DNA fragments. In addition, gene conversion affects the DNA evolution of duplicates in various ways especially when selection operates. Here, we review theoretical models to understand the evolution of duplicates in both neutral and non-neutral cases. We also explain how these theories contribute to interpreting real polymorphism and divergence data by using some intriguing examples.

Sequence variation at multiple loci influences red cell hemoglobin concentration

A substantial genetic contribution underlies variation in baseline peripheral blood counts. We performed quantitative trait locus/loci analyses to identify chromosome regions harboring genes influencing red cell hemoglobin concentration using the cell hemoglobin concentration mean (CHCM), a directly measured parameter analogous to the mean cell hemoglobin concentration. Fourteen significant loci (gene symbols Chcmq1-Chcmq14) were detected. Seven of these influenced CHCM in a sex-specific fashion, and 2 showed significant interactive effects (epistasis). For quantitative trait locus/loci detected in multiple crosses, confidence intervals were narrowed using statistical and bioinformatic approaches. Two strong candidate genes emerged and were further analyzed: adult β-globin (Hbb) for Chcmq3 on Chr 7, and transferrin (Trf) for Chcmq2 on Chr 9. High and low allele parental strains in crosses detecting Chcmq3 segregate 100% with the known ancestral haplotype blocks, hemoglobin (Hb) diffuse (Hbb(d)) and Hb single (Hbb(s)), respectively. Hbb(d) consists of nonidentical major and minor polypeptides and exhibits an increased positive charge relative to Hbb(s) due to the net loss of 2 negative residues in the Hbb(dminor) polypeptide, resulting in a pI of 7.85 versus 7.13. Thus, as shown in human erythrocytes, positively charged Hbs are associated with cell dehydration and increased CHCM in mouse erythrocytes.

Nonallelic gene conversion in the genus Drosophila

Nonallelic gene conversion has been proposed as a major force in homogenizing the sequences of paralogous genes. In this work, we investigate the extent and characteristics of gene conversion among gene families in nine species of the genus Drosophila. We carried out a genome-wide study of 2855 gene families (including 17,742 genes) and determined that conversion events involved 2628 genes. The proportion of converted genes ranged across species from 1 to 9% when paralogs of all ages were included. Although higher levels of gene conversion were found among young gene duplicates, at most 1-2% of the coding sequences of these duplicates were affected by conversion. Using a second approach relying on gene family size changes and gene-tree/species-tree reconciliation methods, we estimate that only 1-15% of gene trees are misled by gene conversion, depending on the lineage considered. Several features of paralogous genes correlate with gene conversion, such as intra-/interchromosomal location, level of nucleotide divergence, and GC content, although we found no definitive evidence for biased substitution patterns. After considering species-specific differences in the age and distance between paralogs, we found a highly significant difference in the amount of gene conversion among species. In particular, members of the melanogaster group showed the lowest proportion of converted genes. Our data therefore suggest underlying differences in the mechanistic basis of gene conversion among species.

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.

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

The functional diversification of multigene families may be strongly influenced by mechanisms of concerted evolution such as interparalog gene conversion. The beta-globin gene family of house mice (genus Mus) represents an especially promising system for evaluating the effects of gene conversion on the functional divergence of duplicated genes. Whereas the majority of mammalian species possess tandemly duplicated copies of the adult beta-globin gene that are identical in sequence, natural populations of house mice are often polymorphic for distinct two-locus haplotypes that differ in levels of functional divergence between duplicated beta-globin genes, HBB-T1 and HBB-T2. Here, we use a phylogenetic approach to unravel the complex evolutionary history of the HBB-T1 and HBB-T2 paralogs in a taxonomically diverse set of species in the genus Mus. The main objectives of this study were 1) to reconstruct the evolutionary history of the different HBB haplotypes of house mice, 2) to assess the role of recombinational exchange between HBB-T1 and HBB-T2 in promoting concerted evolution, 3) to assess the role of recombinational exchange between HBB-T1 and HBB-T2 in creating chimeric genes, and 4) to assess the structural basis of hemoglobin isoform differentiation in species that possess distinct HBB paralogs. Results of our phylogenetic survey revealed that the HBB-T1 and HBB-T2 genes in different species of Mus exhibit the full range of evolutionary outcomes with respect to levels of interparalog divergence. At one end of the spectrum, the two identical HBB paralogs on the Hbb(s) haplotype (shared by Mus domesticus, Mus musculus, and Mus spretus) represent a classic example of concerted evolution. At the other end of the spectrum, the two distinct HBB paralogs on the Hbb(d), Hbb(p), Hbb(w1), and Hbb(w2) haplotypes (shared by multiple species in the subgenus Mus) show no trace of gene conversion and are distinguished by a number of functionally important amino acid substitutions. Because the possession of distinct HBB paralogs expands the repertoire of functionally distinct hemoglobin isoforms that can be synthesized during fetal development and postnatal life, variation in the level of functional divergence between HBB-T1 and HBB-T2 may underlie important physiological variation within and among species.

Origin and ascendancy of a chimeric fusion gene: the beta/delta-globin gene of paenungulate mammals

The delta-globin gene (HBD) of eutherian mammals exhibits a propensity for recombinational exchange with the closely linked beta-globin gene (HBB) and has been independently converted by the HBB gene in multiple lineages. Here we report the presence of a chimeric beta/delta fusion gene in the African elephant (Loxodonta africana) that was created by unequal crossing-over between misaligned HBD and HBB paralogs. The recombinant chromosome that harbors the beta/delta fusion gene in elephants is structurally similar to the "anti-Lepore" duplication mutant of humans (the reciprocal exchange product of the hemoglobin Lepore deletion mutant). However, the situation in the African elephant is unique in that the chimeric beta/delta fusion gene supplanted the parental HBB gene and is therefore solely responsible for synthesizing the beta-chain subunits of adult hemoglobin. A phylogenetic survey of beta-like globin genes in afrotherian and xenarthran mammals revealed that the origin of the chimeric beta/delta fusion gene and the concomitant inactivation of the HBB gene predated the radiation of "Paenungulata," a clade of afrotherian mammals that includes three orders: Proboscidea (elephants), Sirenia (dugongs and manatees), and Hyracoidea (hyraxes). The reduced fitness of the human Hb Lepore deletion mutant helps to explain why independently derived beta/delta fusion genes (which occur on an anti-Lepore chromosome) have been fixed in a number of mammalian lineages, whereas the reciprocal delta/beta fusion gene (which occurs on a Lepore chromosome) has yet to be documented in any nonhuman mammal. This illustrates how the evolutionary fates of chimeric fusion genes can be strongly influenced by their recombinational mode of origin.

Species specificity in major urinary proteins by parallel evolution

Species-specific chemosignals, pheromones, regulate social behaviors such as aggression, mating, pup-suckling, territory establishment, and dominance. The identity of these cues remains mostly undetermined and few mammalian pheromones have been identified. Genetically-encoded pheromones are expected to exhibit several different mechanisms for coding 1) diversity, to enable the signaling of multiple behaviors, 2) dynamic regulation, to indicate age and dominance, and 3) species-specificity. Recently, the major urinary proteins (Mups) have been shown to function themselves as genetically-encoded pheromones to regulate species-specific behavior. Mups are multiple highly related proteins expressed in combinatorial patterns that differ between individuals, gender, and age; which are sufficient to fulfill the first two criteria. We have now characterized and fully annotated the mouse Mup gene content in detail. This has enabled us to further analyze the extent of Mup coding diversity and determine their potential to encode species-specific cues.

Our results show that the mouse Mup gene cluster is composed of two subgroups: an older, more divergent class of genes and pseudogenes, and a second class with high sequence identity formed by recent sequential duplications of a single gene/pseudogene pair. Previous work suggests that truncated Mup pseudogenes may encode a family of functional hexapeptides with the potential for pheromone activity. Sequence comparison, however, reveals that they have limited coding potential. Similar analyses of nine other completed genomes find Mup gene expansions in divergent lineages, including those of rat, horse and grey mouse lemur, occurring independently from a single ancestral Mup present in other placental mammals. Our findings illustrate that increasing genomic complexity of the Mup gene family is not evolutionarily isolated, but is instead a recurring mechanism of generating coding diversity consistent with a species-specific function in mammals.

Balancing selection at a frog antimicrobial peptide locus: fluctuating immune effector alleles?

Balancing selection is common on many defense genes, but it has rarely been reported for immune effector proteins such as antimicrobial peptides (AMPs). We describe genetic diversity at a brevinin-1 AMP locus in three species of leopard frogs (Rana pipiens, Rana blairi, and Rana palustris). Several highly divergent allelic lineages are segregating at this locus. That this unusual pattern results from balancing selection is demonstrated by multiple lines of evidence, including a ratio of nonsynonymous/synonymous polymorphism significantly higher than 1, the ZnS test, incongruence between the number of segregating sites and haplotype diversity, and significant Tajima's D values. Our data are more consistent with a model of fluctuating selection in which alleles change frequencies over time than with a model of stable balancing selection such as overdominance. Evidence for fluctuating selection includes skewed allele frequencies, low levels of synonymous variation, nonneutral values of Tajima's D within allelic lineages, an inverse relationship between the frequency of an allelic lineage and its degree of polymorphism, and divergent allele frequencies among populations. AMP loci could be important sites of adaptive genetic diversity, with consequences for host–pathogen coevolution and the ability of species to resist disease epidemics.

New genes originated via multiple recombinational pathways in the beta-globin gene family of rodents

Species differences in the size or membership composition of multigene families can be attributed to lineage-specific additions of new genes via duplication, losses of genes via deletion or inactivation, and the creation of chimeric genes via domain shuffling or gene fusion. In principle, it should be possible to infer the recombinational pathways responsible for each of these different types of genomic change by conducting detailed comparative analyses of genomic sequence data. Here, we report an attempt to unravel the complex evolutionary history of the beta-globin gene family in a taxonomically diverse set of rodent species. The main objectives were: 1) to characterize the genomic structure of the beta-globin gene cluster of rodents; 2) to assign orthologous and paralogous relationships among duplicate copies of beta-like globin genes; and 3) to infer the specific recombinational pathways responsible for gene duplications, gene deletions, and the creation of chimeric fusion genes. Results of our comparative genomic analyses revealed that variation in gene family size among rodent species is mainly attributable to the differential gain and loss of later expressed beta-globin genes via unequal crossing-over. However, two distinct recombinational mechanisms were implicated in the creation of chimeric fusion genes. In muroid rodents, a chimeric gamma/epsilon fusion gene was created by unequal crossing-over between the embryonic epsilon- and gamma-globin genes. Interestingly, this gamma/epsilon fusion gene was generated in the same fashion as the "anti-Lepore" 5'-delta-(beta/delta)-beta-3' duplication mutant in humans (the reciprocal exchange product of the pathological hemoglobin Lepore deletion mutant). By contrast, in the house mouse, Mus musculus, a chimeric beta/delta fusion pseudogene was created by a beta-globin --> delta-globin gene conversion event. Although the gamma/epsilon and beta/delta fusion genes share a similar chimeric gene structure, they originated via completely different recombinational pathways.

Differential loss of embryonic globin genes during the radiation of placental mammals

The differential gain and loss of genes from homologous gene families represents an important source of functional variation among the genomes of different species. Differences in gene content between species are primarily attributable to lineage-specific gene gains via duplication and lineage-specific losses via deletion or inactivation. Here, we use a comparative genomic approach to investigate this process of gene turnover in the beta-globin gene family of placental mammals. By analyzing genomic sequence data from representatives of each of the main superordinal clades of placental mammals, we were able to reconstruct pathways of gene family evolution during the basal radiation of this physiologically and morphologically diverse vertebrate group. Our analysis revealed that an initial expansion of the nonadult portion of the beta-globin gene cluster in the ancestor of placental mammals was followed by the differential loss and retention of ancestral gene lineages, thereby generating variation in the complement of embryonic globin genes among contemporary species. The sorting of epsilon-, gamma-, and eta-globin gene lineages among the basal clades of placental mammals has produced species differences in the functional types of hemoglobin isoforms that can be synthesized during the course of embryonic development.

Adaptive functional divergence among triplicated alpha-globin genes in rodents

The functional divergence of duplicated genes is thought to play an important role in the evolution of new developmental and physiological pathways, but the role of positive selection in driving this process remains controversial. The objective of this study was to test whether amino acid differences among triplicated alpha-globin paralogs of the Norway rat (Rattus norvegicus) and the deer mouse (Peromyscus maniculatus) are attributable to a relaxation of purifying selection or to a history of positive selection that has adapted the gene products to new or modified physiological tasks. In each rodent species, the two paralogs at the 5'-end of the alpha-globin gene cluster (HBA-T1 and HBA-T2) are evolving in concert and are therefore identical or nearly identical in sequence. However, in each case, the HBA-T1 and HBA-T2 paralogs are distinguished from the third paralog at the 3'-end of the gene cluster (HBA-T3) by multiple amino acid substitutions. An analysis of genomic sequence data from several rodent species revealed that the HBA-T3 genes of Rattus and Peromyscus originated via independent, lineage-specific duplication events. In the independently derived HBA-T3 genes of both species, a likelihood analysis based on a codon-substitution model revealed that accelerated rates of amino acid substitution are attributable to positive directional selection, not to a relaxation of purifying selection. As a result of functional divergence among the triplicated alpha-globin genes in Rattus and Peromyscus, the red blood cells of both rodent species contain a mixture of functionally distinct alpha-chain hemoglobin isoforms that are predicted to have different oxygen-binding affinities. In P. maniculatus, a species that is able to sustain physiological function under conditions of chronic hypoxia at high altitude, the coexpression of distinct hemoglobin isoforms with graded oxygen affinities is expected to broaden the permissible range of arterial oxygen tensions for pulmonary/tissue oxygen transport.

Genomic evidence for independent origins of beta-like globin genes in monotremes and therian mammals

Phylogenetic reconstructions of the beta-globin gene family in vertebrates have revealed that developmentally regulated systems of hemoglobin synthesis have been reinvented multiple times in independent lineages. For example, the functional differentiation of embryonic and adult beta-like globin genes occurred independently in birds and mammals. In both taxa, the embryonic beta-globin gene is exclusively expressed in primitive erythroid cells derived from the yolk sac. However, the "epsilon-globin" gene in birds is not orthologous to the epsilon-globin gene in mammals, because they are independently derived from lineage-specific duplications of a proto beta-globin gene. Here, we report evidence that the early and late expressed beta-like globin genes in monotremes and therian mammals (marsupials and placental mammals) are the products of independent duplications of a proto beta-globin gene in each of these two lineages. Results of our analysis of genomic sequence data from a large number of vertebrate taxa, including sequence from the recently completed platypus genome, reveal that the epsilon- and beta-globin genes of therian mammals arose via duplication of a proto beta-globin gene after the therian/monotreme split. Our analysis of genomic sequence from the platypus also revealed the presence of a duplicate pair of beta-like globin genes that originated via duplication of a proto beta-globin gene in the monotreme lineage. This discovery provides evidence that, in different lineages of mammals, descendent copies of the same proto beta-globin gene may have been independently neofunctionalized to perform physiological tasks associated with oxygen uptake and storage during embryonic development.

Rapid rates of lineage-specific gene duplication and deletion in the alpha-globin gene family

Phylogeny reconstructions of the globin gene families have revealed that paralogous genes within species are often more similar to one another than they are to their orthologous counterparts in closely related species. This pattern has been previously attributed to mechanisms of concerted evolution such as interparalog gene conversion that homogenize sequence variation between tandemly duplicated genes and therefore create the appearance of recent common ancestry. Here we report a comparative genomic analysis of the alpha-globin gene family in mammals that reveal a surprisingly high rate of lineage-specific gene duplication and deletion via unequal crossing-over. Results of our analysis reveal that patterns of sequence similarity between paralogous alpha-like globin genes from the same species are only partly explained by concerted evolution between preexisting gene duplicates. In a number of cases, sequence similarity between paralogous sequences from the same species is attributable to recent ancestry between the products of de novo gene duplications. As a result of this surprisingly rapid rate of gene gain and loss, many mammals possess alpha-like globin genes that have no orthologous counterparts in closely related species. The resultant variation in gene copy number among species may represent an important source of regulatory variation that affects physiologically important aspects of blood oxygen transport and aerobic energy metabolism.