Signatures of selection and gene conversion associated with human color vision variation

Trichromatic color vision in humans results from the combination of red, green, and blue photopigment opsins. Although color vision genes have been the targets of active molecular and psychophysical research on color vision abnormalities, little is known about patterns of normal genetic variation in these genes among global human populations. The current study presents nucleotide sequence analyses and tests of neutrality for a 5.5-kb region of the X-linked long-wave "red" opsin gene (OPN1LW) in 236 individuals from ethnically diverse human populations. Our analysis of the recombination landscape across OPN1LW reveals an unusual haplotype structure associated with amino acid replacement variation in exon 3 that is consistent with gene conversion. Compared with the absence of OPN1LW amino acid replacement fixation since divergence from chimpanzee, the human population exhibits a significant excess of high-frequency OPN1LW replacements. Our results suggest that subtle changes in L-cone opsin wavelength absorption may have been adaptive during human evolution.     

Diversifying selection of the tumor-growth promoter angiogenin in primate evolution

Diversifying selection drives the rapid differentiation of gene sequences and is one of the main forces behind adaptive evolution. Most genes known to be shaped by diversifying selection are those involved in host-pathogen or male-female interactions characterized as molecular "arms races." Here we report the unexpected detection of diversifying selection in the evolution of a tumor-growth promoter, angiogenin (ANG). A comparison among 11 primate species demonstrates that ANG has a significantly higher rate of nucleotide substitution at nonsynonymous sites than at synonymous sites, a hallmark of positive selection acting at the molecular level. Furthermore, we observed significant charge diversity at the molecular surface, suggesting the presence of selective pressures in the microenvironment of ANG, including its binding molecules. A population survey of ANG in chimpanzees, however, reveals no polymorphism, which may have resulted from a recent selective sweep of a charge-altering substitution in chimpanzee evolution. Functional assays of recombinant ANGs from the human and owl monkey indicate that primate ANGs retain angiogenic activity despite rapid evolution. Our study, together with findings of similar selection in the primate breast cancer suppressor gene, BRCA1, reveals an intriguing phenomenon of unusual selective pressures on, and adaptive evolution of, cancer-related genes in primate evolution.     

Deleterious mutations at the mitochondrial ND3 gene in South American marsh rats (Holochilus).

Statistical analyses of DNA sequences have revealed patterns of nonneutral evolution in mitochondrial DNA of mice, humans, and Drosophila. Here we report patterns of mitochondrial sequence evolution in South American marsh rats (genus Holochilus). We sequenced the complete mitochondrial ND3 gene in 82 Holochilus brasiliensis and 21 H. vulpinus to test the neutral prediction that the ratio of nonsynonymous to synonymous nucleotide changes is the same within and between species. Within H. brasiliensis we observed a greater number of amino acid polymorphisms than expected based on interspecific comparisons. This contingency table analysis suggests that many amino acid polymorphisms are mildly deleterious. Several tests of the frequency distribution also revealed departures from a neutral, equilibrium model, and these departures were observed for both nonsynonymous and synonymous sites. In general, an excess of rare sites was observed, consistent with either a recent selective sweep or with populations not at mutation-drift equilibrium.     

Evidence for balancing selection from nucleotide sequence analyses of human G6PD

Glucose-6-phosphate dehydrogenase (G6PD) mutations that result in reduced enzyme activity have been implicated in malarial resistance and constitute one of the best examples of selection in the human genome. In the present study, we characterize the nucleotide diversity across a 5.2-kb region of G6PD in a sample of 160 Africans and 56 non-Africans, to determine how selection has shaped patterns of DNA variation at this gene. Our global sample of enzymatically normal B alleles and A, A-, and Med alleles with reduced enzyme activities reveals many previously uncharacterized silent-site polymorphisms. In comparison with the absence of amino acid divergence between human and chimpanzee G6PD sequences, we find that the number of G6PD amino acid polymorphisms in human populations is significantly high. Unlike many other G6PD-activity alleles with reduced activity, we find that the age of the A variant, which is common in Africa, may not be consistent with the recent emergence of severe malaria and therefore may have originally had a historically different adaptive function. Overall, our observations strongly support previous genotype-phenotype association studies that proposed that balancing selection maintains G6PD deficiencies within human populations. The present study demonstrates that nucleotide sequence analyses can reveal signatures of both historical and recent selection in the genome and may elucidate the impact that infectious disease has had during human evolution.     

Structural diversity of the classical H-2 genes: K, D, and L.

Twenty-three class I DNA sequences, representing alleles of the H-2K, D, and L loci, were analyzed to assess patterns of nucleotide and amino acid diversity. Comparisons of the allelic and nonallelic sequences revealed locus specificity in regions encoding the leader peptides and the carboxyl-terminal segments of the Ag presenting molecules. Analyses focusing on the sequences that determine the Ag binding domains revealed weak or insignificant allelic associations, a finding that is in sharp contrast to previously observed relationships among the homologous human sequences. The amino acid positions exhibiting high diversity in the encoded glycoproteins in both mice and humans are localized primarily to the Ag binding site. In the mouse, diverse amino acids were positioned similarly in the K and D/L glycoproteins, although in humans, the A and B glycoproteins exhibit distinctive differences in their locations within the Ag binding site. The absence of locus specificity among the sequences that determine the Ag binding domains of the mouse is consistent with the hypothesis that ectopic gene conversion leads to interlocus exchange of class I sequences. Comparable interlocus exchanges among human class I genes have not played a similar role in shaping human A and B sequences. The basis of this difference between mice and humans is not clear. The nature of amino acid substitutions distinguishing class I loci in mice and humans are comparable, and the role of natural selection in determining diversity appears to be similar in the two species.     

Positive selection driving the evolution of a gene of male reproduction, Acp26Aa, of Drosophila: II. Divergence versus polymorphism

The evolution of the gene for a male ejaculatory protein, Acp26Aa, has been shown to be driven by positive selection when nonsibling species in the Drosophila melanogaster subgroup are compared. To know if selection has been operating in the recent past and to understand the details of its dynamics, we obtained DNA sequences of Acp26Aa and the nearby Acp26Ab gene from 39 D. melanogaster chromosomes. Together with the 10 published sequences, we analyzed 49 sequences from five populations in four continents. The southern African population is somewhat differentiated from all other populations, but its nucleotide diversity is lower at these two loci. We find the following results for Acp26Aa: (1) The R: S (replacement : silent changes) ratio is significantly higher in the between-species comparisons than in the within-species data by the McDonald and Kreitman test. Positive selection is probably responsible for the excess of amino acid replacements between species. (2) However, within-species nucleotide diversity is high. Neither the Tajima test nor the Fu and Li test indicates a reduction in nucleotide diversity due to positive selection in the recent past. (3) The newly derived nucleotides in D. melanogaster are at high frequency significantly more often than predicted by the neutral equilibrium. Since the nearby Acp26Ab gene does not show these patterns, these observations cannot be attributed to the characteristics of this chromosomal region. We suggest that positive selection is active, but may be weak, for each amino acid change in the Acp26Aa gene.      

Identification of Single Nucleotide Polymorphisms and Analysis of Linkage Disequilibrium in Different Bamboo Species Using the Candidate Gene Approach

Bamboos are one of the most beautiful and useful plants on Earth. The genetic background and population structure of bamboos are well known, which helps accelerate the process of artificial domestication of bamboo. Partial sequences of six genes involved in nitrogen use efficiency in 32 different bamboo species were analyzed for occurrence of single nucleotide polymorphisms (SNPs). The nucleotide diversity θw and total nucleotide polymorphisms π_T of the sequenced DNA regions was 0.05137 and 0.03332, respectively. Both πnonsyn /πsyn and Ka/Ks values were <1. The nucleotide sequences of these six genes were inferred to be relatively conserved, and the haplotype diversity was relatively high. The results of evolutionary neutrality tests showed that the six genes were in line with neutral evolution, and that the NRT2.1 and AMT2.1 gene sequences may have experienced negative selection. An inter-SNP recombination event at the NRT2.1 gene in the all pooled sample, of all 32 bamboo species was the lowest at 0.0645, whereas the AMT gene recombination events were all >0.1. Estimation and analysis of linkage disequilibrium of five genes revealed that with the increase in nucleotide sequence length, the degree of SNP linkage disequilibrium decreased rapidly. We inferred the population genetic structure of 32 bamboo species based on the SNP loci of six genes with frequencies >18%. 32 bamboo species were divided into five categories, which indicated that the combined population of all bamboo species had obvious multivariate characteristics and was heterogeneous; red (Group 1) and green (Group 2) were the main groups.     

Signatures of functional constraint at aye-aye opsin genes: the potential of adaptive color vision in a nocturnal primate

While color vision perception is thought to be adaptively correlated with foraging efficiency for diurnal mammals, those that forage exclusively at night may not need color vision nor have the capacity for it. Indeed, although the basic condition for mammals is dichromacy, diverse nocturnal mammals have only monochromatic vision, resulting from functional loss of the short-wavelength sensitive opsin gene. However, many nocturnal primates maintain intact two opsin genes and thus have dichromatic capacity. The evolutionary significance of this surprising observation has not yet been elucidated. We used a molecular population genetics approach to test evolutionary hypotheses for the two intact opsin genes of the fully nocturnal aye-aye (Daubentonia madagascariensis), a highly unusual and endangered Madagascar primate. No evidence of gene degradation in either opsin gene was observed for any of 8 aye-aye individuals examined. Furthermore, levels of nucleotide diversity for opsin gene functional sites were lower than those for 15 neutrally evolving intergenic regions (>25 kb in total), which is consistent with a history of purifying selection on aye-aye opsin genes. The most likely explanation for these findings is that dichromacy is advantageous for aye-ayes despite their nocturnal activity pattern. We speculate that dichromatic nocturnal primates may be able to perceive color while foraging under moonlight conditions, and suggest that behavioral and ecological comparisons among dichromatic and monochromatic nocturnal primates will help to elucidate the specific activities for which color vision perception is advantageous.     

Multilocus patterns of polymorphism and selection across the X chromosome of Caenorhabditis remanei

Natural selection and neutral processes such as demography, mutation, and gene conversion all contribute to patterns of polymorphism within genomes. Identifying the relative importance of these varied components in evolution provides the principal challenge for population genetics. To address this issue in the nematode Caenorhabditis remanei, I sampled nucleotide polymorphism at 40 loci across the X chromosome. The site-frequency spectrum for these loci provides no evidence for population size change, and one locus presents a candidate for linkage to a target of balancing selection. Selection for codon usage bias leads to the non-neutrality of synonymous sites, and despite its weak magnitude of effect (N(e)s approximately 0.1), is responsible for profound patterns of diversity and divergence in the C. remanei genome. Although gene conversion is evident for many loci, biased gene conversion is not identified as a significant evolutionary process in this sample. No consistent association is observed between synonymous-site diversity and linkage-disequilibrium-based estimators of the population recombination parameter, despite theoretical predictions about background selection or widespread genetic hitchhiking, but genetic map-based estimates of recombination are needed to rigorously test for a diversity-recombination relationship. Coalescent simulations also illustrate how a spurious correlation between diversity and linkage-disequilibrium-based estimators of recombination can occur, due in part to the presence of unbiased gene conversion. These results illustrate the influence that subtle natural selection can exert on polymorphism and divergence, in the form of codon usage bias, and demonstrate the potential of C. remanei for detecting natural selection from genomic scans of polymorphism.     

Recombination and selection at Brassica self-incompatibility loci

In Brassica species, self-incompatibility is controlled genetically by haplotypes involving two known genes, SLG and SRK, and possibly an as yet unknown gene controlling pollen incompatibility types. Alleles at the incompatibility loci are maintained by frequency-dependent selection, and diversity at SLG and SRK appears to be very ancient, with high diversity at silent and replacement sites, particularly in certain "hypervariable" portions of the genes. It is important to test whether recombination occurs in these genes before inferences about function of different parts of the genes can be made from patterns of diversity within their sequences. In addition, it has been suggested that, to maintain the relationship between alleles within a given S-haplotype, recombination is suppressed in the S-locus region. The high diversity makes many population genetic measures of recombination inapplicable. We have analyzed linkage disequilibrium within the SLG gene of two Brassica species, using published coding sequences. The results suggest that intragenic recombination has occurred in the evolutionary history of these alleles. This is supported by patterns of synonymous nucleotide diversity within both the SLG and SRK genes, and between domains of the SRK gene. Finally, clusters of linkage disequilibrium within the SLG gene suggest that hypervariable regions are under balancing selection, and are not merely regions of relaxed selective constraint.     

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.     

DNA variation at the rp49 gene region of Drosophila simulans: evolutionary inferences from an unusual haplotype structure.

An approximately 1.3-kb region including the rp49 gene plus its 5' and 3' flanking regions was sequenced in 24 lines of Drosophila simulans (10 from Spain and 14 from Mozambique). Fifty-four nucleotide and 8 length polymorphisms were detected. All nucleotide polymorphisms were silent: 52 in noncoding regions and 2 at synonymous sites in the coding region. Estimated silent nucleotide diversity was similar in both populations (pi = 0.016, for the total sample). Nucleotide variation revealed an unusual haplotype structure showing a subset of 11 sequences with a single polymorphism. This haplotype was present at intermediate frequencies in both the European and the African samples. The presence of such a major haplotype in a highly recombining region is incompatible with the neutral equilibrium model. This haplotype structure in both a derived and a putatively ancestral population can be most parsimoniously explained by positive selection. As the rate of recombination in the rp49 region is high, the target of selection should be close to or within the region studied.     

Mosaic nature of the wolbachia surface protein

Lateral gene transfer and recombination play important roles in the evolution of many parasitic bacteria. Here we investigate intragenic recombination in Wolbachia bacteria, considered among the most abundant intracellular bacteria on earth. We conduct a detailed analysis of the patterns of variation and recombination within the Wolbachia surface protein, utilizing an extensive set of published and new sequences from five main supergroups of Wolbachia. Analysis of nucleotide and amino acid sequence variations confirms four hypervariable regions (HVRs), separated by regions under strong conservation. Comparison of shared polymorphisms reveals a complex mosaic structure of the gene, characterized by a clear intragenic recombining of segments among several distinct strains, whose major recombination effect is shuffling of a relatively conserved set of amino acid motifs within each of the four HVRs. Exchanges occurred both within and between the arthropod supergroups. Analyses based on phylogenetic methods and a specific recombination detection program (MAXCHI) significantly support this complex partitioning of the gene, indicating a chimeric origin of wsp. Although wsp has been widely used to define macro- and microtaxonomy among Wolbachia strains, these results clearly show that it is not suitable for this purpose. The role of wsp in bacterium-host interactions is currently unknown, but results presented here indicate that exchanges of HVR motifs are favored by natural selection. Identifying host proteins that interact with wsp variants should help reveal how these widespread bacterial parasites affect and evolve in response to the cellular environments of their invertebrate hosts.     

Positive selection and interallelic recombination at the merozoite surface antigen-1 (MSA-1) locus of Plasmodium falciparum.

DNA sequences of alleles at the merozoite surface antigen-1 (MSA-1) gene locus of the malaria parasite Plasmodium falciparum show evidence of repeated past recombination events between alleles. These include both (1) nonreciprocal recombination events that have homogenized certain gene regions among alleles and (2) reciprocal recombination events that have combined allelic segments with divergent evolutionary histories, thereby enhancing allelic diversity. In three different gene regions, the rate of nonsynonymous nucleotide substitution significantly exceeds that of synonymous nucleotide substitution, implying that positive Darwinian selection has acted to diversify alleles at the amino acid level. The MSA-1 polymorphism seems to be quite ancient; the two major allelic types have been maintained for approximately 35 Myr.     

Population genetic analysis of the N-acylsphingosine amidohydrolase gene associated with mental activity in humans

To understand the evolution of human mental activity, we performed population genetic analyses of nucleotide sequences ( approximately 11 kb) from a worldwide sample of 60 chromosomes of the N-acylsphingosine amidohydrolase (ASAH1) gene. ASAH1 hydrolyzes ceramides and regulates neuronal development, and its deficiency often results in mental retardation. In the region ( approximately 4.4 kb) encompassing exons 3 and 4 of this gene, two distinct lineages (V and M) have been segregating in the human population for 2.4 +/- 0.4 million years (MY). The persistence of these two lineages is attributed to ancient population structure of humans in Africa. However, all haplotypes belonging to the V lineage exhibit strong linkage disequilibrium, a high frequency (62%), and small nucleotide diversity (pi = 0.05%). These features indicate a signature of positive Darwinian selection for the V lineage. Compared with the orthologs in mammals and birds, it is only Val at amino acid site 72 that is found exclusively in the V lineage in humans, suggesting that this Val is a likely target of positive selection. Computer simulation confirms that demographic models of modern humans except for the ancient population structure cannot explain the presence of two distinct lineages, and neutrality is incompatible with the observed small genetic variation of the V lineage at ASAH1. On the basis of the above observations, it is argued that positive selection is possibly operating on ASAH1 in the modern human population.     

Population genetics of Neisseria gonorrhoeae in a high-prevalence community using a hypervariable outer membrane porB and 13 slowly evolving housekeeping genes

Baltimore, Md., is an urban community with a high prevalence of Neisseria gonorrhoeae. Due to partially protective immune responses, introduction of new strains from other host populations, and exposure of N. gonorrhoeae to antibiotics, the phenotypic and genotypic characteristics of the circulating strains can fluctuate over time. Understanding the overall genetic diversity and population structure of N. gonorrhoeae is essential for informing public health interventions to eliminate this pathogen. We studied gonococci population genetics in Baltimore by analyzing a hypervariable and strongly selected outer membrane porB gene and 13 slowly evolving and presumably neutral housekeeping genes (abcZ, adk, aroE, fumC, gdh, glnA, gnd, pdhC, pgm, pilA, ppk, pyrD, and serC) in 204 isolates collected in 1991, 1996, and 2001 from male and female patients of two public sexually transmitted diseases clinics. Genetic diversity (), recombination (C), growth (g), population structure, and adaptive selection under codon-substitution and amino acid property models were estimated and compared between these two gene classes. Estimates of the F(ST) fixation index and the chi(2) test of sequence absolute frequencies revealed significant temporal substructuring for both gene types. Baltimore's N. gonorrhoeae populations have increased since 1991 as indicated by consistent positive values of g. Female patients showed similar or lower levels of and C than male patients. Within the MLST housekeeping genes, levels of and C ranged from 0.001-0.013 and 0.000-0.018, respectively. Overall recombination seems to be the dominant force driving evolution in these populations. All loci showed amino acid sites and physicochemical properties under adaptive (or positive-destabilizing) selection, rejecting the generally assumed hypothesis of stabilizing selection for these MLST genes. Within the porB gene, protein I B showed higher and C values than protein I A. Directional positive selection possibly mediated by the immune system operates to a significant extent in the protein I sequences, as indicated by the distribution of the positively selected sites in the surface-exposed loops. Thirteen amino acid physicochemical properties seem to drive protein evolution of the PI porins in N. gonorrhoeae.     

Adaptive and Slightly Deleterious Evolution in a Conifer, <Emphasis Type="Italic">Cryptomeria japonica</Emphasis>

In order to evaluate effects of the population structure and natural selection on organisms having long generation times, we surveyed DNA polymorphisms at five loci encoding 9-cis-epoxycarotenoid dioxygenase (NCED), ammonium transporter, calmodulin, aquaporin, and the second major allergen with polymethylgalacturonase enzyme activity in the pollen (Cryj2) in a conifer, Cryptomeria japonica. The average nucleotide diversity at silent sites across 12 loci including the previously analyzed seven loci was 0.0044. The population recombination rate (4Nr, where N and r are the effective population size and recombination rate per base per generation, respectively) was estimated as 0.00046 and a slow reduction in the population size was indicated, according to the maximum likelihood method implemented in LAMARC. At NCED, the McDonald-Kreitman (MK) test revealed an excess of replacement polymorphisms, suggesting contributions of slightly deleterious mutations. In contrast, the MK test revealed an excess of replacement divergence at Cryj2 and a maximum likelihood approach using the PAML package revealed that certain amino acid sites had a nonsynonymous/synonymous substitution rate ratio (omega) > 4.0, indicating adaptive evolution at this locus. The overall analysis of the 12 loci suggested that adaptive, neutral, and slightly deleterious evolution played important roles in the evolution of C. japonica.     

Increased constraints on MC4R during primate and human evolution

The melanocortin 4 receptor (MC4R) is routinely investigated for the role it plays in human obesity, as mutations in MC4R are the most common dominantly inherited form of the disease. As little is known about the evolutionary history of this locus, we investigated patterns of variation at MC4R in a worldwide sample of 1,015 humans from 51 populations, and in 8 central chimpanzees. There is a significant paucity of diversity at MC4R in humans, but not in chimpanzees. The spectrum of mutations in humans, combined with the overall low level of diversity, suggests that most (if not all) of the observed non-synonymous polymorphisms are likely to be transient deleterious mutations. The MC4R coding region was resequenced in 12 primate species and sequences from an additional 29 vertebrates were included in molecular evolutionary analyses. MC4R is highly conserved throughout vertebrate evolution, and has apparently been subject to high levels of continuous purifying selection that increased approximately threefold during primate evolution. Furthermore, the strong selection extends to codon usage bias, where most silent mutations are expected to be either quickly fixed or removed from the population, which may help explain the unusually low levels of silent polymorphisms in humans. Finally, there is a significant tendency for non-synonymous mutations that impact MC4R function to occur preferentially at sites that are identified by evolutionary analyses as being subject to very strong purifying selection. The information from this study should help inform future epidemiological investigations of MC4R. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Mitochondrial-nuclear interactions and accelerated compensatory evolution: evidence from the primate cytochrome C oxidase complex

Accelerated rates of mitochondrial protein evolution have been proposed to reflect Darwinian coadaptation for efficient energy production for mammalian flight and brain activity. However, several features of mammalian mtDNA (absence of recombination, small effective population size, and high mutation rate) promote genome degradation through the accumulation of weakly deleterious mutations. Here, we present evidence for "compensatory" adaptive substitutions in nuclear DNA- (nDNA) encoded mitochondrial proteins to prevent fitness decline in primate mitochondrial protein complexes. We show that high mutation rate and small effective population size, key features of primate mitochondrial genomes, can accelerate compensatory adaptive evolution in nDNA-encoded genes. We combine phylogenetic information and the 3D structure of the cytochrome c oxidase (COX) complex to test for accelerated compensatory changes among interacting sites. Physical interactions among mtDNA- and nDNA-encoded components are critical in COX evolution; amino acids in close physical proximity in the 3D structure show a strong tendency for correlated evolution among lineages. Only nuclear-encoded components of COX show evidence for positive selection and adaptive nDNA-encoded changes tend to follow mtDNA-encoded amino acid changes at nearby sites in the 3D structure. This bias in the temporal order of substitutions supports compensatory weak selection as a major factor in accelerated primate COX evolution.