Comparison of DNA and protein polymorphisms between humans and chimpanzees

The Impact of Linked Selection in Chimpanzees: A Comparative Study

Levels of nucleotide diversity vary greatly across the genomes of most species owing to multiple factors. These include variation in the underlying mutation rates, as well as the effects of both direct and linked selection. Fundamental to interpreting the relative importance of these forces is the common observation of a strong positive correlation between nucleotide diversity and recombination rate. While indeed observed in humans, the interpretation of this pattern has been difficult in the absence of high-quality polymorphism data and recombination maps in closely related species. Here, we characterize genetic features driving nucleotide diversity in Western chimpanzees using a recently generated whole genome polymorphism data set. Our results suggest that recombination rate is the primary predictor of nucleotide variation with a strongly positive correlation. In addition, telomeric distance, regional GC-content, and regional CpG-island content are strongly negatively correlated with variation. These results are compared with humans, with both similarities and differences interpreted in the light of the estimated effective population sizes of the two species as well as their strongly differing recent demographic histories.

A human-specific allelic group of the MHC DRB1 gene in primates

Background

Diversity among human leukocyte antigen (HLA) molecules has been maintained by host-pathogen coevolution over a long period of time. Reflecting this diversity, the HLA loci are the most polymorphic in the human genome. One characteristic of HLA diversity is long-term persistence of allelic lineages, which causes trans-species polymorphisms to be shared among closely related species. Modern humans have disseminated across the world after their exodus from Africa, while chimpanzees have remained in Africa since the speciation event between humans and chimpanzees. It is thought that modern humans have recently acquired resistance to novel pathogens outside Africa. In the present study, we investigated HLA alleles that could contribute to this local adaptation in humans and also studied the contribution of natural selection to human evolution by using molecular data.

Results

Phylogenetic analysis of HLA-DRB1 genes identified two major groups, HLA Groups A and B. Group A formed a monophyletic clade distinct from DRB1 alleles in other Catarrhini, suggesting that Group A is a human-specific allelic group. Our estimates of divergence time suggested that seven HLA-DRB1 Group A allelic lineages in humans have been maintained since before the speciation event between humans and chimpanzees, while chimpanzees possess only one DRB1 allelic lineage (Patr-DRB1*03), which is a sister group to Group A. Experimental data showed that some Group A alleles bound to peptides derived from human-specific pathogens. Of the Group A alleles, three exist at high frequencies in several local populations outside Africa.

Conclusions

HLA Group A alleles are likely to have been retained in human lineages for a long period of time and have not expanded since the divergence of humans and chimpanzees. On the other hand, most orthologs of HLA Group A alleles may have been lost in the chimpanzee due to differences in selective pressures. The presence of alleles with high frequency outside of Africa suggests these HLA molecules result from the local adaptations of humans. Our study helps elucidate the mechanism by which the human adaptive immune system has coevolved with pathogens over a long period of time.

Weak selection and protein evolution

The "nearly neutral" theory of molecular evolution proposes that many features of genomes arise from the interaction of three weak evolutionary forces: mutation, genetic drift, and natural selection acting at its limit of efficacy. Such forces generally have little impact on allele frequencies within populations from generation to generation but can have substantial effects on long-term evolution. The evolutionary dynamics of weakly selected mutations are highly sensitive to population size, and near neutrality was initially proposed as an adjustment to the neutral theory to account for general patterns in available protein and DNA variation data. Here, we review the motivation for the nearly neutral theory, discuss the structure of the model and its predictions, and evaluate current empirical support for interactions among weak evolutionary forces in protein evolution. Near neutrality may be a prevalent mode of evolution across a range of functional categories of mutations and taxa. However, multiple evolutionary mechanisms (including adaptive evolution, linked selection, changes in fitness-effect distributions, and weak selection) can often explain the same patterns of genome variation. Strong parameter sensitivity remains a limitation of the nearly neutral model, and we discuss concave fitness functions as a plausible underlying basis for weak selection.

Parasite-mediated evolution of the functional part of the MHC in primates

The major histocompatibility complex (MHC) is a key model of genetic polymorphism, but the mechanisms underlying its extreme variability are debated. Most hypotheses for MHC diversity focus on pathogen-driven selection and predict that MHC polymorphism evolves under the pressure of a diverse parasite fauna. Several studies reported that certain alleles offer protection against certain parasites, yet it remains unclear whether variation in parasite pressure more generally covaries with allelic diversity and rates of molecular evolution of MHC across species. We tested this prediction in a comparative study of 41 primate species. We characterized polymorphism of the exon 2 of DRB region of the MHC class II. Our phylogenetic analyses controlled for the potential effects of neutral mutation rate, population size, geographic origin and body mass and revealed that nematode species richness associates positively with nonsynonymous nucleotide substitution rate at the functional part of the molecule. We failed to find evidence for allelic diversity being strongly related to parasite species richness. Continental distribution was a strong predictor of both allelic diversity and substitution rate, with higher values in Malagasy and Neotropical primates. These results indicate that parasite pressure can influence the different estimates of MHC polymorphism, whereas geography plays an independent role in the natural history of MHC.

Gene flow between species of Lake Victoria haplochromine fishes

The haplochromine cichlid fishes of Lake Victoria (LV), East Africa, are a textbook example of adaptive radiation-a rapid divergence of multiple morphologically distinguishable forms from a few founding lineages. The forms are generally believed to constitute a "flock" of several hundred reproductively isolated species in a dozen or so genera. This belief has, until now, not been subjected to a test, however. Here, we compare genetic variation at 11 loci in 10 haplochromine populations of 6 different species. Although the genetic diversity in the populations is quite high, using a variety of statistical tests, we find no evidence of genetic differentiation among the populations of LV haplochromines. On genetic distance trees, populations of the same species intermingle with those of different species. At the molecular level, the species are indistinguishable from one another. Genetic comparisons with closely related species in 2 crater lakes indicate that the species within LV continue exchanging genes. These observations have important implications for phylogenetic reconstruction. The approach used in this study is applicable to other instances of adaptive radiation.

Ancestral population sizes and species divergence times in the primate lineage on the basis of intron and BAC end sequences

The effective sizes of ancestral populations and species divergence times of six primate species (humans, chimpanzees, gorillas, orangutans, and representatives of Old World monkeys and New World monkeys) are estimated by applying the two-species maximum likelihood (ML) method to intron sequences of 20 different loci. Examination of rate heterogeneity of nucleotide substitutions and intragenic recombination identifies five outrageous loci (ODC1, GHR, HBE, INS, and HBG). The estimated ancestral polymorphism ranges from 0.21 to 0.96% at major divergences in primate evolution. One exceptionally low polymorphism occurs when African and Asian apes diverged. However, taking into consideration the possible short generation times in primate ancestors, it is concluded that the ancestral population size in the primate lineage was no smaller than that of extant humans. Furthermore, under the assumption of 6 million years (myr) divergence between humans and chimpanzees, the divergence time of humans from gorillas, orangutans. Old World monkeys, and New World monkeys is estimated as 7.2, 18, 34, and 65 myr ago, respectively, which are generally older than traditional estimates. Beside the intron sequences, three other data sets of orthologous sequences are used between the human and the chimpanzee comparison. The ML application to these data sets including 58,156 random BAC end sequences (BES) shows that the nucleotide substitution rate is as low as 0.6-0.8 x 10(-9) per site per year and the extent of ancestral polymorphism is 0.33-0.51%. With such a low substitution rate and short generation time, the relatively high extent of polymorphism suggests a fairly large effective population size in the ancestral lineage common to humans and chimpanzees.

Differential alu mobilization and polymorphism among the human and chimpanzee lineages

Alu elements are primate-specific members of the SINE (short interspersed element) retroposon family, which comprise approximately 10% of the human genome. Here we report the first chromosomal-level comparison examining the Alu retroposition dynamics following the divergence of humans and chimpanzees. We find a twofold increase in Alu insertions in humans in comparison to the common chimpanzee (Pan troglodytes). The genomic diversity (polymorphism for presence or absence of the Alu insertion) associated with these inserts indicates that, analogous to recent nucleotide diversity studies, the level of chimpanzee Alu diversity is approximately 1.7 times higher than that of humans. Evolutionarily recent Alu subfamily structure differs markedly between the human and chimpanzee lineages, with the major human subfamilies remaining largely inactive in the chimpanzee lineage. We propose a population-based model to account for the observed fluctuation in Alu retroposition rates across primate taxa.

A prevalent POLG CAG microsatellite length allele in humans and African great apes

The human nuclear gene for the catalytic subunit of mitochondrial DNA polymerase gamma ( POLG) contains within its coding region a CAG microsatellite encoding a polyglutamine repeat. Previous studies demonstrated an association between length variation at this repeat and male infertility, suggesting a mechanism whereby the prevalent (CAG)(10) allele, which occurs at a frequency of >80% in different populations, could be maintained by selection. Sequence analysis of the POLG CAG microsatellite region of more than 1000 human chromosomes reveals that virtually all allelic variation at the locus is accounted for by length variation of the CAG repeat. Analysis of POLG from African great apes shows that a prevalent length allele is present in each species, although its exact length is species-specific. In common chimpanzee ( Pan troglodytes) a number of different sequence variants contribute to the prevalent length allele, strongly supporting the idea that the length of the POLG microsatellite region, rather than its exact nucleotide or amino acid sequence, is what is maintained. Analysis of POLG in other primates indicates that the repeat has expanded from a shorter, glutamine-rich sequence, present in the common ancestor of Old and New World monkeys.

Comparative linkage-disequilibrium analysis of the beta-globin hotspot in primates

Recombination rates vary both across the genome and between different species, but little information is available about the temporal and physical scales over which such rates change. To shed light on these questions, we performed a high-resolution analysis of a genomic region within the beta-globin gene cluster that is known to experience elevated recombination rates in humans. For this purpose, we developed new linkage disequilibrium-based methods that thoroughly search for subsets of the data with unusually high or unusually low estimated values of the population-recombination parameter (4Nr, where N is the effective population size and r is the crossover rate between adjacent base pairs). By resequencing a 15-kb segment in a human population sample, we were able to narrow the recombinational hotspot to a segment <2 kb in length that coincides with the beta-globin replication origin. In addition, we analyzed the orthologous region in samples of rhesus macaques and common chimpanzees. Whereas the analysis of the chimpanzee data is complicated by the sample structure, the macaque data imply that this region may not be a hotspot in that species. These results suggest a time scale for the evolution of hotspots in primates. Furthermore, they allow us to propose diverged sequence elements that may contribute to the differences in the recombinational landscape in the two species.

Low nucleotide diversity in chimpanzees and bonobos

Comparison of the levels of nucleotide diversity in humans and apes may provide much insight into the mechanisms of maintenance of DNA polymorphism and the demographic history of these organisms. In the past, abundant mitochondrial DNA (mtDNA) polymorphism data indicated that nucleotide diversity (pi) is more than threefold higher in chimpanzees than in humans. Furthermore, it has recently been claimed, on the basis of limited data, that this is also true for nuclear DNA. In this study we sequenced 50 noncoding, nonrepetitive DNA segments randomly chosen from the nuclear genome in 9 bonobos and 17 chimpanzees. Surprisingly, the pi value for bonobos is only 0.078%, even somewhat lower than that (0.088%) for humans for the same 50 segments. The pi values are 0.092, 0.130, and 0.082% for East, Central, and West African chimpanzees, respectively, and 0.132% for all chimpanzees. These values are similar to or at most only 1.5 times higher than that for humans. The much larger difference in mtDNA diversity than in nuclear DNA diversity between humans and chimpanzees is puzzling. We speculate that it is due mainly to a reduction in effective population size (N(e)) in the human lineage after the human-chimpanzee divergence, because a reduction in N(e) has a stronger effect on mtDNA diversity than on nuclear DNA diversity. Sequence data from this article have been deposited with the GenBank Data libraries under accession nos. AY 275957-AY 277244.