Effects of spatially varying selection on nucleotide diversity and linkage disequilibrium: insights from deer mouse globin genes

An important goal of population genetics is to elucidate the effects of natural selection on patterns of DNA sequence variation. Here we report results of a study to assess the joint effects of selection, recombination, and gene flow in shaping patterns of nucleotide variation at genes involved in local adaptation. We first describe a new summary statistic, Z(g), that measures the between-sample component of linkage disequilibrium (LD). We then report results of a multilocus survey of nucleotide diversity and LD between high- and low-altitude populations of deer mice, Peromyscus maniculatus. The multilocus survey included two closely linked alpha-globin genes, HBA-T1 and HBA-T2, that underlie adaptation to different elevational zones. The primary goals were to assess whether the alpha-globin genes exhibit the hallmarks of spatially varying selection that are predicted by theory (i.e., sharply defined peaks in the between-population components of nucleotide diversity and LD) and to assess whether peaks in diversity and LD may be useful for identifying specific sites that distinguish selectively maintained alleles. Consistent with theoretical expectations, HBA-T1 and HBA-T2 were characterized by highly elevated levels of diversity between populations and between allele classes. Simulation and empirical results indicate that sliding-window analyses of Z(g) between allele classes may provide an effective means of pinpointing causal substitutions.     

Comparative analysis of testis protein evolution in rodents

Genes expressed in testes are critical to male reproductive success, affecting spermatogenesis, sperm competition, and sperm-egg interaction. Comparing the evolution of testis proteins at different taxonomic levels can reveal which genes and functional classes are targets of natural and sexual selection and whether the same genes are targets among taxa. Here we examine the evolution of testis-expressed proteins at different levels of divergence among three rodents, mouse (Mus musculus), rat (Rattus norvegicus), and deer mouse (Peromyscus maniculatus), to identify rapidly evolving genes. Comparison of expressed sequence tags (ESTs) from testes suggests that proteins with testis-specific expression evolve more rapidly on average than proteins with maximal expression in other tissues. Genes with the highest rates of evolution have a variety of functional roles including signal transduction, DNA binding, and egg-sperm interaction. Most of these rapidly evolving genes have not been identified previously as targets of selection in comparisons among more divergent mammals. To determine if these genes are evolving rapidly among closely related species, we sequenced 11 of these genes in six Peromyscus species and found evidence for positive selection in five of them. Together, these results demonstrate rapid evolution of functionally diverse testis-expressed proteins in rodents, including the identification of amino acids under lineage-specific selection in Peromyscus. Evidence for positive selection among closely related species suggests that changes in these proteins may have consequences for reproductive isolation.     

Peromyscus maniculatus--Mus musculus chromosome homology map derived from reciprocal cross species chromosome painting

The Mus musculus and Rattus norvegicus genomes have been extensively studied, yet despite the emergence of Peromyscus maniculatus as an NIH model for genome sequencing and biomedical research much remains unknown about the genome organization of Peromyscines. Contrary to their phylogenetic relationship, the genomes of Rattus and Peromyscus appear more similar at the gross karyotypic level than either does to Mus. We set out to define the chromosome homologies between Peromyscus, Mus and Rattus. Reciprocal cross-species chromosome painting and G-band homology assignments were used to delineate the conserved chromosome homology map between P. maniculatus and M. musculus. These data show that each species has undergone extensive chromosome rearrangements since they last shared a common ancestor 25 million years ago (mya). This analysis coupled with an inferred homology map with Rattus revealed a high level of chromosome conservation between Rattus and Peromyscus and indicated that the chromosomes of Mus are highly derived.     

Cross-transmission studies with Eimeria arizonensis, E. arizonensis-like oocysts and Eimeria langebarteli: host specificity at the genus and species level within the Muridae

Cross-transmission experiments were done using sporulated oocysts of Eimeria arizonensis from Peromyscus truei and Peromyscus maniculatus, and oocysts of 2 putative species that resemble E. arizonensis, i.e., Eimeria albigulae from Neotoma albigula, and Eimeria onychomysis from Onychomys leucogaster. Oocysts of each species were inoculated into representatives of P. maniculatus and the latter 2 rodent species. Other experiments were conducted wherein oocysts of Eimeria langebarteli from Peromyscus leucopus were given to P. truei and P. maniculatus. Oocysts of E. arizonensis from P. truei and P. maniculatus could be transmitted only to P. maniculatus; likewise, oocysts of E. albigulae and E. onychomysis produced patent infections only in N. albigula and O. leucogaster, respectively. Oocysts of E. langebarteli from P. leucopus could be transmitted to P. truei, but not P. maniculatus. These results indicate that E. arizonensis, and the morphologically similar E. albigulae and E. onychomysis, are distinct species that are not transmissible between the genera of their respective hosts (Peromyscus, Neotoma, Onychomys), and that some isolates of E. langebarteli, reported from 6 species of Peromyscus and Reithrodontomys megalotis, may not always be infective to P. maniculatus.     

Systematic implications of chromosomal data from two insular species of Peromyscus from the Gulf of California

G- and C-banded karyotypes for two insular species of deer mice, Peromyscus slevini and P. sejugis, are described and analyzed relative to the evolutionary relationship of these species to and their inclusion within the P. maniculatus species group. The chromosomal phenotype of P. slevini is unique among all banded karyotypes reported for Peromyscus, and comparison with published karyotypes suggests that P. slevini has systematic affinities with either the P. boylii or P. mexicanus species groups. The karyotypic data for P. sejugis clearly align these mice with P. maniculatus and provide a diagnostic character that supports the specific distinction between these taxa.     

MHC class II beta sequence diversity in the deer mouse (Peromyscus maniculatus): implications for models of balancing selection

We studied population polymorphism at a major histocompatibility complex (MHC) class II beta gene in the deer mouse (Peromyscus maniculatus). We found that: (i) a single population of P. maniculatus has significantly higher levels of DNA and protein sequence diversity than worldwide samples from homologous genes in other taxa, including humans and mice; and (ii) the genealogy of allelic sequences in P. maniculatus deviates significantly from theoretical expectation under a model of symmetric balancing selection, in that alleles are relatively more divergent than expected. We suggest that the observation of high levels of pairwise allelic sequence divergence and deviation of the genealogy from theoretical expectation in P. maniculatus together provide support for a divergent allele advantage model for the maintenance of MHC polymorphism.     

Selection on amino acid substitutions in Arabidopsis

Studies of nucleotide diversity have found an excess of low-frequency amino acid polymorphisms segregating in Arabidopsis thaliana, suggesting a predominance of weak purifying selection acting on amino acid polymorphism in this inbreeding species. Here, we investigate levels of diversity and divergence at synonymous and nonsynonymous sites in 6 circumpolar populations of the outbreeding Arabidopsis lyrata and compare these results with A. thaliana, to test for differences in mutation and selection parameters across genes, populations, and species. We find that A. lyrata shows an excess of low-frequency nonsynonymous polymorphisms both within populations and species wide, consistent with weak purifying selection similar to the patterns observed in A. thaliana. Furthermore, nonsynonymous polymorphisms tend to be more restricted in their population distribution in A. lyrata, consistent with purifying selection preventing their geographic spread. Highly expressed genes show a reduced ratio of amino acid to synonymous change for both polymorphism and fixed differences, suggesting a general pattern of stronger purifying selection on high-expression proteins.     

Positive selection in phytotoxic protein-encoding genes of Botrytis species

Evolutionary patterns of sequence divergence were analyzed in genes from the fungal genus Botrytis (Ascomycota), encoding phytotoxic proteins homologous to a necrosis and ethylene-inducing protein from Fusarium oxysporum. Fragments of two paralogous genes (designated NEP1 and NEP2) were amplified from all known Botrytis species and sequenced. NEP1 sequences of two Botrytis species contain premature stop codons, indicating that they may be non-functional. Both paralogs of all species encode proteins with a remarkably similar predicted secondary structure, however, they contain different types of post-translational modification motifs, which are conserved across the genus. While both NEP genes are, overall, under purifying selection, we identified a number of amino acids under positive selection based on inference using maximum likelihood models. Positively selected amino acids in NEP1 were not under selection in corresponding positions in NEP2. The biological significance of positively selected residues and the role of NEP proteins in pathogenesis remain to be resolved.     

Is promiscuity associated with enhanced selection on MHC-DQα in mice (genus Peromyscus)?

Reproductive behavior may play an important role in shaping selection on Major Histocompatibility Complex (MHC) genes. For example, the number of sexual partners that an individual has may affect exposure to sexually transmitted pathogens, with more partners leading to greater exposure and, hence, potentially greater selection for variation at MHC loci. To explore this hypothesis, we examined the strength of selection on exon 2 of the MHC-DQα locus in two species of Peromyscus. While the California mouse (P. californicus) is characterized by lifetime social and genetic monogamy, the deer mouse (P. maniculatus) is socially and genetically promiscuous; consistent with these differences in mating behavior, the diversity of bacteria present within the reproductive tracts of females is significantly greater for P. maniculatus. To test the prediction that more reproductive partners and exposure to a greater range of sexually transmitted pathogens are associated with enhanced diversifying selection on genes responsible for immune function, we compared patterns and levels of diversity at the Class II MHC-DQα locus in sympatric populations of P. maniculatus and P. californicus. Using likelihood based analyses, we show that selection is enhanced in the promiscuous P. maniculatus. This study is the first to compare the strength of selection in wild sympatric rodents with known differences in pathogen milieu.     

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.     

Anthropogenetical analysis of abnormal human alpha-globin gene cluster arrangement on chromosome 16

An earlier study of human globin gene polymorphism in two Adriatic islands of Olib and Silba showed an abnormal arrangement of alpha-globin genes in two different individuals. The next step was to determine the degree of the kinship relationship between the two probands, one with a deleted and another with triplicated alpha-globin gene on the island Silba, and to determine the stability of this disorder through generations. We reviewed the parish registers (Status Animarum) of the island of Silba, dating from the year 1527, and constructed family trees for the two probands. Restriction endonuclease mapping was performed to study the arrangement of the alpha-globin genes in the offspring of our probands. A total of 183 ancestors completed the two family trees. The kinship relationship between them was established in the 5th, 6th, and 7th generation. The analysis of alpha-globin genes in the offspring of our probands showed the triplicated alpha-globin genes in two persons. We also found alpha-globin gene triplication in other three relatives. We did not find any deleted alpha-globin genes. We determined the kinship relationship between the two probands, one with deleted and the other with triplicated alpha-globin genes. This finding enabled us to determine the stability of this gene disarrangement through generations. It also showed new possibilities in anthropogenetic research, by combining the analyses of parish registers with those of modern genetic methods, such as restriction endonuclease mapping.     

Evolution and diversity of axon guidance Robo receptor family genes

The diversity of axon guidance (AG) receptors reflects gains in complexity of the animal nervous system during evolution. Members of the Roundabout (Robo) family of receptors interact with Slit proteins and play important roles in many developmental processes, including AG and neural crest cell migration. There are four members of the Robo gene family. However, the evolutionary history of Robo family genes remain obscure. We analyzed the distribution of Robo family members in metazoan species ranging in complexity from hydras to humans. We undertook a phylogenetic analysis in metazoans, synteny analysis, and ancestral chromosome mapping in vertebrates, and detected selection pressure and functional divergence among four mammalian Robo paralogs. Based on our analysis, we proposed that the ancestral Robo gene could have undergone a tandem duplication in the vertebrate ancestor; then one round of whole genome duplication events occurred before the divergence of ancestral lamprey and gnathostome, generating four paralogs in early vertebrates. Robo4 paralog underwent segmental loss in the following evolutionary process. Our results showed that Robo3 paralog is under more powerful purifying selection pressure compared with other three paralogs, which could correlate with its unique expression pattern and function. Furthermore, we found four sites under positive selection pressure on the Ig1‐2 domains of Robo4 that might interfere with its binding to Slits ligand. Diverge analysis at the amino acid level showed that Robo4 paralog have relatively greater functional diversifications than other Robo paralogs. This coincides with the fact that Robo4 predominantly functions in vascular endothelial cells but not the nervous system.     

Kappa-chain constant-region gene sequences in genus Rattus: coding regions are diverging more rapidly than noncoding regions

We have determined the nucleotide sequence of a 1,200-base pair (bp) genomic fragment that includes the kappa-chain constant-region gene (C kappa) from two species of native Australian rodents, Rattus leucopus cooktownensis and Rattus colletti. Comparison of these sequences with each other and with other rodent C kappa genes shows three surprising features. First, the coding regions are diverging at a rate severalfold higher than that of the nearby noncoding regions. Second, replacement changes within the coding region are accumulating at a rate at least as great as that of silent changes. Third, most of the amino acid replacements are localized in one region of the C kappa domain--namely, the carboxy-terminal "bends" in the alpha-carbon backbone. These three features have previously been described from comparisons of the two allelic forms of C kappa genes in R. norvegicus. These data imply the existence of considerable evolutionary constraints on the noncoding regions (based on as yet undetermined functions) or powerful positive selection to diversify a portion of the constant-region domain (whose physiological significance is not known). These surprising features of C kappa evolution appear to be characteristic only of closely related C kappa genes, since comparison of rodent with human sequences shows the expected greater conservation of coding regions, as well as a predominance of silent nucleotide substitutions within the coding regions.      

ALPHA-CHAIN HEMOGLOBIN POLYMORPHISMS ARE CORRELATED WITH ALTITUDE IN THE DEER MOUSE,PEROMYSCUS MANICULATUS

In deer mouse (Peromyscus maniculatus) populations in the western United States, alpha-globin haplotype frequency, beta-globin haplotype frequency, and base-line blood oxygen affinity (measured after acclimation to low altitude) show strong correlations with native altitude. The correlations improve when an average regional altitude is substituted for the local altitude at collection sites. This substitution roughly compensates for the effects of gene exchange between populations in areas of highly variable topography. When subspecific effects are removed with covariate analyses a significant (P < 0.05) relationship remains only for alpha-globin haplotype frequency and altitude. Thus, alpha-globin haplotype frequency, beta-globin haplotype frequency, and base-line blood oxygen affinity may be explained by either subspecific or altitudinal effects, but subspecific effects explain a larger proportion of the variance. Part of the subspecific effect may be attributable to an underlying relationship of subspecies with altitude. The analyses for the alpha-globins in conjunction with other data on the effects of alpha-globins on blood oxygen affinity and whole-animal physiological performance are consistent with the hypothesis that the frequency of the alpha-globins evolved in response to selection resulting from the stress of high-altitude hypoxia.     

Serologic evidence of hantavirus infection in sigmodontine rodents in Mexico

Antibodies to hantaviruses in two species of sigmodontine rodents (Peromyscus maniculatus and Reithrodontomys sumichrasti) collected in central Mexico are reported. Peromyscus maniculatus, a common species throughout much of Mexico, is the reservoir of Sin Nombre virus (SNV), the etiologic agent of the great majority of cases of hantavirus pulmonary syndrome (HPS) in North America. Although the identity of the virus detected in P. maniculatus in Mexico could not be determined by these serologic results, our findings suggest that SNV may occur throughout the range of P. maniculatus in North America. If true, the failure to identify HPS in Mexico is not due to the absence of pathogenic hantaviruses in Mexico.     

Diverse selective modes among orthologs/paralogs of the chalcone synthase (Chs) gene family of Arabidopsis thaliana and its relative A. halleri ssp. gemmifera

As a model system, Arabidopsis thaliana and its wild relatives have played an important role in the study of genomics and evolution in plants. In this study, we examined the genetic diversity of the chalcone synthase (Chs) gene, which encodes a key enzyme of the flavonoid pathway and is located on chromosome five, as well as two Chs-like genes on the first and fourth chromosomes of Arabidopsis. The objectives of the study are to determine if natural selection operates differentially on the paralogs of the Chs gene family in A. thaliana and Arabidopsis halleri ssp. gemmifera. The mode of selection was inferred from Tajima's D values from noncoding and coding regions, as well as from the ratio of nonsynonymous to synonymous substitutions. Both McDonald-Kreitman and HKA tests revealed the effects of selection on the allelic distribution, except for the chromosome 1 paralog in ssp. gemmifera. The Chs gene on chromosome 5 was under purifying selection in both species. Significant, negative Tajima's D values at synonymous sites and positive Fay and Wu's H values within coding region, plus reduced genetic variability in introns, indicated effects of background selection in shaping the evolution of this gene region in A. thaliana. The Chs paralog on chromosome 1 was under positive selection in A. thaliana, while interspecific introgression and balancing selection determined the fates of the paralog and resulted in high heterogeneity in ssp. gemmifera. Local adaptation differentiated populations of Japan and China at the locus. In contrast, the other Chs-paralog of chromosome 4 was shaped by purifying selection in A. thaliana, while under positive selection in ssp. gemmifera, as indicated by dn/ds>1. Moreover, these contrasting patterns of selection have likely resulted in functional divergence in Arabidopsis, as indicated by radical amino acid substitutions at the chalcone synthase/stilbene synthase motif of the Chs genes. Unlike previous studies of the evolutionary history of A. thaliana, the high levels of genetic diversity in most gene regions of Chs paralogs and nonsignificant Tajima's D in the intron sequences of the Chs gene family in A. thaliana did not reflect the effects of a recent demographic expansion.     

Two types of triplicated alpha-globin loci in humans.

DNA from healthy Malaysian newborns was studied on gene maps after digestion with different restriction endonucleases. Of 65 newborns, two were found to be carriers of two different variants of triplicated alpha-globin loci. In variant no. 1, found in an Malay, the three alpha-globin genes are in an elongated DNA fragment on digestion with Eco RI and Bam HI. The third alpha-globin gene was found in a additional 3.7-kb fragment on digestion with Hpa I, Bgl II and Hind III. In variant no. 2, a new type of triplicated alpha-globin loci, found in a Chinese, the three alpha-globin genes reside in an elongated DNA fragment longer than that of variant no. 1 on digestion with Eco RI and Bam HI. The third alpha-globin gene was found in an additional 4.2-kb fragment on digestion with Hpa I and Hind III. Digestion of this variant DNA with Bg1 II produced an abnormal 16.7-kb fragment in addition to the normal 7.0-kb Bgl-II fragment. The locations of the restriction sites in the two types of triplicated alpha-globin loci are compatible with a mechanism of unequal crossing over following two different modes of misalignment.     

A Nonsynonymous/Synonymous Substitution Analysis of the B56 Gene Family Aids in Understanding B56 Isoform Diversity

Gene duplication leads to the formation of gene families, wherein purifying or neutral selection maintains the original gene function, while diversifying selection confers new functions onto duplicated genes. The B56 gene family is highly conserved; it is encoded by one gene in protists and fungi, and five genes in vertebrates. B56 regulates protein phosphatase 2A (PP2A), an abundant heterotrimeric serine/threonine phosphatase that functions as a tumor suppressor and consists of a scaffolding “A” and catalytic “C” subunit heterodimer bound to a regulatory “B” subunit. Individual regulatory B56 subunits confer disparate functions onto PP2A in various cell-cell signaling pathways. B56 proteins share a conserved central core domain, but have divergent N- and C-termini which play a role in isoform specificity. We carried out a nonsynonymous/synonymous substitution analysis to better understand the divergence of vertebrate B56 genes. When five B56 paralogs from ten vertebrate species were analyzed, the gene family displayed purifying selection; stronger purifying selection was revealed when individual B56 isoforms were analyzed separately. The B56 core experienced stronger purifying selection than the N- and C-termini, which correlates with the presence of several contacts between the core and the AC heterodimer. Indeed, the majority of the contact points that we analyzed between B56 and the AC heterodimer experienced strong purifying selection. B56 subfamilies showed distinct patterns of selection in their N- and C-termini. The C-terminus of the B56-1 subfamily and the N-terminus of the B56-2 subfamily exhibited strong purifying selection, suggesting that these termini carry out subfamily-specific functions, while the opposite termini exhibited diversifying selection and likely carry out isoform-specific functions. We also found reduced synonymous substitutions at the N- and C-termini when grouping B56 genes by species but not by isoform, suggesting species-specific codon bias may have a role in regulating B56 gene expression.     

Not born equal: increased rate asymmetry in relocated and retrotransposed rodent gene duplicates

Duplicated genes frequently evolve at different rates. This asymmetry is evidence of natural selection's ability to discriminate between the 2 copies, subjecting them to different levels of purifying selection or even permitting adaptive evolution of one or both copies. However, if gene duplication creates pairs of protein-coding sequences that are initially identical, this raises the question of how selection tells the 2 copies apart. Here, we investigated asymmetric sequence divergence of recently duplicated genes in rodents and related this to 2 possible sources of such asymmetry: gene relocation as a consequence of duplication and retrotransposition as a mechanism of gene duplication. We found that most young rodent duplicates that have been relocated were created by retrotransposition. The degree of rate asymmetry in gene pairs where one copy has been relocated (either by retrotransposition or DNA-based duplication) is greater than in pairs formed by local DNA-based duplication events. Furthermore, by considering the direction of transposition for distant duplicates, we found a consistent tendency for retrogenes to undergo accelerated protein evolution relative to their static paralogs, whereas DNA-based transpositions showed no such tendency. Finally, we demonstrate that the faster sequence evolution of retrogenes correlates with the profound alteration of their expression pattern that is precipitated by retrotransposition.