Low nucleotide diversity in chimpanzees and bonobos

Different selective pressures shape the molecular evolution of color vision in chimpanzee and human populations

A population genetic analysis of the long-wavelength opsin (OPN1LW, "red") color vision gene in a global sample of 236 human nucleotide sequences had previously discovered nine amino acid replacement single nucleotide polymorphisms, which were found at high frequencies in both African and non-African populations and associated with an unusual haplotype diversity. Although this pattern of nucleotide diversity is consistent with balancing selection, it has been argued that a recombination "hot spot" or gene conversion within and between X-linked color vision genes alone may explain these patterns. The current analysis investigates a closely related primate with trichromatism to determine whether color vision gene amino acid polymorphism and signatures of adaptive evolution are characteristic of humans alone. Our population sample of 56 chimpanzee (Pan troglodytes) OPN1LW sequences shows three singleton amino acid polymorphisms and no unusual recombination or linkage disequilibrium patterns across the approximately 5.5-kb region analyzed. Our comparative population genetic approach shows that the patterns of OPN1LW variation in humans and chimpanzees are consistent with positive and purifying selection within the two lineages, respectively. Although the complex role of color vision has been greatly documented in primate evolution in general, it is surprising that trichromatism has followed very different selective trajectories even between humans and our closest relatives.

Genetic subspecies diversity of the chimpanzee CD4 virus-receptor gene

Chimpanzees are naturally and asymptomatically infected by simian immunodeficiency virus (SIV). Pathogenic properties of SIV/HIV vary and differences in susceptibility and pathogenicity of SIV/HIV depend in part on host-specific factors such as virus-receptor/co-receptor interactions. Since CD4 plays a primary role in virus binding and since SIVcpz have been found only in two African chimpanzee subspecies, we characterized the genetic diversity of CD4 receptors in all four recognized subspecies of chimpanzees. We found noticeable variation in the first variable region V1 of CD4 and in intron six among the subspecies of chimpanzees. We found the CD4 receptor to be conserved in individuals belonging to the P. t. verus subspecies and divergent from the other three subspecies, which harbored highly variable CD4 receptors. The CD4 receptor of chimpanzees differed from that of humans. We question whether the observed diversity can explain the species-specific differences in susceptibility to and pathogenicity of SIV/HIV.

Near neutrality: leading edge of the neutral theory of molecular evolution

The nearly neutral theory represents a development of Kimura's neutral theory of molecular evolution that makes testable predictions that go beyond a mere null model. Recent evidence has strongly supported several of these predictions, including the prediction that slightly deleterious variants will accumulate in a species that has undergone a severe bottleneck or in cases where recombination is reduced or absent. Because bottlenecks often occur in speciation and slightly deleterious mutations in coding regions will usually be nonsynonymous, we should expect that the ratio of nonsynonymous to synonymous fixed differences between species should often exceed the ratio of nonsynonymous to synonymous polymorphisms within species. Many data support this prediction, although they have often been wrongly interpreted as evidence for positive Darwinian selection. The use of conceptually flawed tests for positive selection has become widespread in recent years, seriously harming the quest for an understanding of genome evolution. When properly analyzed, many (probably most) claimed cases of positive selection will turn out to involve the fixation of slightly deleterious mutations by genetic drift in bottlenecked populations. Slightly deleterious variants are a transient feature of evolution in the long term, but they have substantially affected contemporary species, including our own.

Adaptive and slightly deleterious evolution in a conifer, Cryptomeria japonica

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.

Analysis of chimpanzee history based on genome sequence alignments

Population geneticists often study small numbers of carefully chosen loci, but it has become possible to obtain orders of magnitude for more data from overlaps of genome sequences. Here, we generate tens of millions of base pairs of multiple sequence alignments from combinations of three western chimpanzees, three central chimpanzees, an eastern chimpanzee, a bonobo, a human, an orangutan, and a macaque. Analysis provides a more precise understanding of demographic history than was previously available. We show that bonobos and common chimpanzees were separated ∼1,290,000 years ago, western and other common chimpanzees ∼510,000 years ago, and eastern and central chimpanzees at least 50,000 years ago. We infer that the central chimpanzee population size increased by at least a factor of 4 since its separation from western chimpanzees, while the western chimpanzee effective population size decreased. Surprisingly, in about one percent of the genome, the genetic relationships between humans, chimpanzees, and bonobos appear to be different from the species relationships. We used PCR-based resequencing to confirm 11 regions where chimpanzees and bonobos are not most closely related. Study of such loci should provide information about the period of time 5–7 million years ago when the ancestors of humans separated from those of the chimpanzees.

Studies of population history traditionally examine a small number of genetic regions in many individuals; however, with genome sequencing technologies it is possible to assemble data sets with thousands more aligned sequences albeit in fewer individuals. To explore whether such data can provide useful insights about population history, we assembled large-scale data sets consisting of overlaps of random genome sequencing reads from chimpanzees and bonobos. Analysis of these data finds that bonobos and chimpanzees split from each other about 1.29 million years ago, western and central chimpanzees about 0.51 million years ago, and eastern and central chimpanzees at least 50,000 years ago. We find that the chimpanzee population has fluctuated significantly in size over the past half million years, with the central chimpanzee population size expanding dramatically, and the western chimpanzee population size contracting. Surprisingly, we also find that there are widespread regions of the genome where chimpanzees and bonobos are less closely related to each other than any of them are to humans. In these regions, chimpanzees and bonobos share a common genetic ancestor dating back to speciation from humans, providing a new source of information about that evolutionary event.

Life-history traits drive the evolutionary rates of mammalian coding and noncoding genomic elements

A comprehensive phylogenetic framework is indispensable for investigating the evolution of genomic features in mammals as a whole, and particularly in humans. Using the ENCODE sequence data, we estimated mammalian neutral evolutionary rates and selective pressures acting on conserved coding and noncoding elements. We show that neutral evolutionary rates can be explained by the generation time (GT) hypothesis. Accordingly, primates (especially humans), having longer GTs than other mammals, display slower rates of neutral evolution. The evolution of constrained elements, particularly of nonsynonymous sites, is in agreement with the expectations of the nearly neutral theory of molecular evolution. We show that rates of nonsynonymous substitutions (dN) depend on the population size of a species. The results are robust to the exclusion of hypermutable CpG prone sites. The average rate of evolution in conserved noncoding sequences (CNCs) is 1.7 times higher than in nonsynonymous sites. Despite this, CNCs evolve at similar or even lower rates than nonsynonymous sites in the majority of basal branches of the eutherian tree. This observation could be the result of an overall gradual or, alternatively, lineage-specific relaxation of CNCs. The latter hypothesis was supported by the finding that 3 of the 20 longest CNCs displayed significant relaxation of individual branches. This observation may explain why the evolution of CNCs fits the expectations of the nearly neutral theory less well than the evolution of nonsynonymous sites.

Statistical evaluation of alternative models of human evolution

An appropriate model of recent human evolution is not only important to understand our own history, but it is necessary to disentangle the effects of demography and selection on genome diversity. Although most genetic data support the view that our species originated recently in Africa, it is still unclear if it completely replaced former members of the Homo genus, or if some interbreeding occurred during its range expansion. Several scenarios of modern human evolution have been proposed on the basis of molecular and paleontological data, but their likelihood has never been statistically assessed. Using DNA data from 50 nuclear loci sequenced in African, Asian and Native American samples, we show here by extensive simulations that a simple African replacement model with exponential growth has a higher probability (78%) as compared with alternative multiregional evolution or assimilation scenarios. A Bayesian analysis of the data under this best supported model points to an origin of our species approximately 141 thousand years ago (Kya), an exit out-of-Africa approximately 51 Kya, and a recent colonization of the Americas approximately 10.5 Kya. We also find that the African replacement model explains not only the shallow ancestry of mtDNA or Y-chromosomes but also the occurrence of deep lineages at some autosomal loci, which has been formerly interpreted as a sign of interbreeding with Homo erectus.

Genomic Comparisons of Humans and Chimpanzees

The genome consists of the entire DNA present in the nucleus of the fertilized embryo, which is then duplicated in every cell in the body. A draft sequence of the chimpanzee genome is now available, providing opportunities to better understand genetic contributions to human evolution, development, and disease. Sequence differences from the human genome were confirmed to be ∼1% in areas that can be precisely aligned, representing ∼35 million single base-pair differences. Some 45 million nucleotides of insertions and deletions unique to each lineage were also discovered, making the actual difference between the two genomes ∼4%. We discuss the opportunities and challenges that arise from this information and the need for comparison with additional species, as well as population genetic studies. Finally, we present a few examples of interesting findings resulting from genome-wide analyses, candidate gene studies, and combined approaches, emphasizing the pros and cons of each approach.

A new approach to estimate parameters of speciation models with application to apes

How populations diverge and give rise to distinct species remains a fundamental question in evolutionary biology, with important implications for a wide range of fields, from conservation genetics to human evolution. A promising approach is to estimate parameters of simple speciation models using polymorphism data from multiple loci. Existing methods, however, make a number of assumptions that severely limit their applicability, notably, no gene flow after the populations split and no intralocus recombination. To overcome these limitations, we developed a new Markov chain Monte Carlo method to estimate parameters of an isolation-migration model. The approach uses summaries of polymorphism data at multiple loci surveyed in a pair of diverging populations or closely related species and, importantly, allows for intralocus recombination. To illustrate its potential, we applied it to extensive polymorphism data from populations and species of apes, whose demographic histories are largely unknown. The isolation-migration model appears to provide a reasonable fit to the data. It suggests that the two chimpanzee species became reproductively isolated in allopatry approximately 850 Kya, while Western and Central chimpanzee populations split approximately 440 Kya but continued to exchange migrants. Similarly, Eastern and Western gorillas and Sumatran and Bornean orangutans appear to have experienced gene flow since their splits approximately 90 and over 250 Kya, respectively.

Variation and diversity in Homo erectus: a 3D geometric morphometric analysis of the temporal bone

Although the level of taxonomic diversity within the fossil hominin species Homo erectus (sensu lato) is continually debated, there have been relatively few studies aiming to quantify the morphology of this species. Instead, most researchers have relied on qualitative descriptions or the evaluation of nonmetric characters, which in many cases display continuous variation. Also, only a few studies have used quantitative data to formally test hypotheses regarding the taxonomic composition of the "erectus" hypodigm. Despite these previous analyses, however, and perhaps in part due to these varied approaches for assessing variation within specimens typically referred to H. erectus (sensu lato) and the general lack of rigorous statistical testing of how variation within this taxon is partitioned, there is currently little consensus regarding whether this group is a single species, or whether it should instead be split into separate temporal or geographically delimited taxa. In order to evaluate possible explanations for variation within H. erectus, we tested the general hypothesis that variation within the temporal bone morphology of H. erectus is consistent with that of a single species, using great apes and humans as comparative taxa. Eighteen three-dimensional (3D) landmarks of the temporal bone were digitized on a total of 520 extant and fossil hominid crania. Landmarks were registered by Generalized Procrustes Analysis, and Procrustes distances were calculated for comparisons of individuals within and between the extant taxa. Distances between fossil specimens and between a priori groupings of fossils were then compared to the distances calculated within the extant taxa to assess the variation within the H. erectus sample relative to that of known species, subspecies, and populations. Results of these analyses indicate that shape variation within the entire H. erectus sample is generally higher than extant hominid intraspecific variation, and putative H. ergaster specimens are significantly different from other specimens in H. erectus (sensu lato). However, shape distances within geographical groups of H. erectus are also high, and OH 9 and Dmanisi 2280 are morphologically distinct from the Koobi Fora specimens that are sometimes classified as H. ergaster. These findings suggest that, although H. erectus may be composed of multiple species, the differentiation is complex, and specimens cannot easily be grouped geographically or chronologically. Consequently, more complicated scenarios seeking to explain the observed variation within H. erectus must be considered.

Episodes of Natural Selection Shaped the Interactions of IgA-Fc with FcαRI and Bacterial Decoy Proteins

FcalphaRI, a receptor for IgA-Fc, recruits myeloid cells to attack IgA-coated pathogens. By competing with FcalphaRI for IgA, bacterial decoys, like SSL7 of Staphylococcus aureus, subvert this defense. We examined how pathogen selection has driven the diversification and coevolution of IgA and FcalphaRI. In higher primates, the IgA binding site of FcalphaRI diversified under positive selection, a strong episode occurring in hominoid ancestors about the time of the IgA gene duplication. The differential binding of SSL7 to IgA-Fc of different species correlates with substitution at seven positions in IgA-Fc, two of which were positively selected in higher primates. Two others, which reduce SSL7 binding, emerged during episodes of positive selection in the rabbit and rodent lineages. The FcalphaRI-IgA interaction evolves episodically under two types of positive selection: pressure from pathogen decoys selects for IgA escape variants which, in turn, selects for FcalphaRI variants to keep up with the novel IgA. When FcalphaRI cannot keep up, its function is lost and the gene becomes susceptible to elimination, as occurred in the mouse genome, either by chance or selection on one of the many linked, variable immune system genes. A cluster of positively selected residues presents a putative binding site for unknown IgA-binding factors.

Genetic structure of chimpanzee populations

Little is known about the history and population structure of our closest living relatives, the chimpanzees, in part because of an extremely poor fossil record. To address this, we report the largest genetic study of the chimpanzees to date, examining 310 microsatellites in 84 common chimpanzees and bonobos. We infer three common chimpanzee populations, which correspond to the previously defined labels of "western," "central," and "eastern," and find little evidence of gene flow between them. There is tentative evidence for structure within western chimpanzees, but we do not detect distinct additional populations. The data also provide historical insights, demonstrating that the western chimpanzee population diverged first, and that the eastern and central populations are more closely related in time.

Understanding the recent evolution of the human genome: insights from human-chimpanzee genome comparisons

The sequencing of the chimpanzee genome and the comparison with its human counterpart have begun to reveal the spectrum of genetic changes that has accompanied human evolution. In addition to gross karyotypic rearrangements such as the fusion that formed human chromosome 2 and the human-specific pericentric inversions of chromosomes 1 and 18, there is considerable submicroscopic structural variation involving deletions, duplications, and inversions. Lineage-specific segmental duplications, detected by array comparative genomic hybridization and direct sequence comparison, have made a very significant contribution to this structural divergence, which is at least three-fold greater than that due to nucleotide substitutions. Since structural genomic changes may have given rise to irreversible functional differences between the diverging species, their detailed analysis could help to identify the biological processes that have accompanied speciation. To this end, interspecies comparisons have revealed numerous human-specific gains and losses of genes as well as changes in gene expression. The very considerable structural diversity (polymorphism) evident within both lineages has, however, hampered the analysis of the structural divergence between the human and chimpanzee genomes. The concomitant evaluation of genetic divergence and diversity at the nucleotide level has nevertheless served to identify many genes that have evolved under positive selection and may thus have been involved in the development of human lineage-specific traits. Genes that display signs of weak negative selection have also been identified and could represent candidate loci for complex genomic disorders. Here, we review recent progress in comparing the human and chimpanzee genomes and discuss how the differences detected have improved our understanding of the evolution of the human genome.

Accelerated molecular evolution in Microtus (Rodentia) as assessed via complete mitochondrial genome sequences

Microtus is one of the most taxonomically diverse mammalian genera, including over 60 extant species. These rodents have evolved rapidly, as the genus originated less than 2 million years ago. If these numbers are taken at face value, then an average of 30 microtine speciation events have occurred every million years. One explanation for the rapid rate of cladogenesis in Microtus could be the karyotypic differentiation exhibited across the genus: diploid numbers range from 17 to 64. Despite the striking chromosomal variability within Microtus, phenotypic variation is unremarkable. To determine whether nucleotide substitution rates are also elevated in voles, we sequenced the entire mitochondrial DNA (mtDNA) genome of the Eurasian sibling vole (Microtus rossiaemeridionalis). We compared this genome to another previously sequenced vole mtDNA genome (Microtus kikuchii) and performed pairwise sequence comparisons with the mtDNA genomes of ten additional mammalian genera. We found that microtine mtDNA genomes are evolving more rapidly than any other mammalian lineage we sampled, as gauged by the rate of nucleotide substitution across the entire mtDNA genome as well as at each individual protein-coding gene. Additionally, we compared substitution rates within the cytochrome b gene to seven other rodent genera and found that Microtus mtDNA is evolving fastest. The root cause of accelerated evolution in Microtus remains uncertain, but merits further investigation.

Comment on "Population size does not influence mitochondrial genetic diversity in animals"

Bazin et al. (Reports, 28 April, 2006, p. 570) found no relationship between mitochondrial DNA (mtDNA) diversity and population size when comparing across large groups of animals. We show empirically that species with smaller populations, as represented by eutherian mammals, exhibit a positive correlation between mtDNA and allozyme variation, suggesting that mtDNA diversity may correlate with population size in these animals.

Ancient and recent positive selection transformed opioid cis-regulation in humans

Changes in the cis-regulation of neural genes likely contributed to the evolution of our species' unique attributes, but evidence of a role for natural selection has been lacking. We found that positive natural selection altered the cis-regulation of human prodynorphin, the precursor molecule for a suite of endogenous opioids and neuropeptides with critical roles in regulating perception, behavior, and memory. Independent lines of phylogenetic and population genetic evidence support a history of selective sweeps driving the evolution of the human prodynorphin promoter. In experimental assays of chimpanzee-human hybrid promoters, the selected sequence increases transcriptional inducibility. The evidence for a change in the response of the brain's natural opioids to inductive stimuli points to potential human-specific characteristics favored during evolution. In addition, the pattern of linked nucleotide and microsatellite variation among and within modern human populations suggests that recent selection, subsequent to the fixation of the human-specific mutations and the peopling of the globe, has favored different prodynorphin cis-regulatory alleles in different parts of the world.

Combining sperm typing and linkage disequilibrium analyses reveals differences in selective pressures or recombination rates across human populations

A previous polymorphism survey of the type 2 diabetes gene CAPN10 identified a segment showing an excess of polymorphism levels in all population samples, coinciding with localized breakdown of linkage disequilibrium (LD) in a sample of Hausa from Cameroon, but not in non-African samples. This raised the possibility that a recombination hotspot is present in all populations and we had insufficient power to detect it in the non-African data. To test this possibility, we estimated the crossover rate by sperm typing in five non-African men; these estimates were consistent with the LD decay in the non-African, but not in the Hausa data. Moreover, resequencing the orthologous region in a sample of Western chimpanzees did not show either an excess of polymorphism level or rapid LD decay, suggesting that the processes underlying the patterns observed in humans operated only on the human lineage. These results suggest that a hotspot of recombination has recently arisen in humans and has reached higher frequency in the Hausa than in non-Africans, or that there is no elevation in crossover rate in any human population, and the observed variation results from long-standing balancing selection.

Reconstructing human origins in the genomic era

Analyses of recently acquired genomic sequence data are leading to important insights into the early evolution of anatomically modern humans, as well as into the more recent demographic processes that accompanied the global radiation of Homo sapiens. Some of the new results contradict early, but still influential, conclusions that were based on analyses of gene trees from mitochondrial DNA and Y-chromosome sequences. In this review, we discuss the different genetic and statistical methods that are available for studying human population history, and identify the most plausible models of human evolution that can accommodate the contrasting patterns observed at different loci throughout the genome.

Subspecies composition and founder contribution of the captive U.S. chimpanzee (Pan troglodytes) population

Chimpanzees are presently classified into three subspecies: Pan troglodytes verus from west Africa, P.t. troglodytes from central Africa, and P.t. schweinfurthii from east Africa. A fourth subspecies (P.t. vellerosus), from Cameroon and northern Nigeria, has been proposed. These taxonomic designations are based on geographical origins and are reflected in sequence variation in the first hypervariable region (HVR-I) of the mtDNA D-loop. Although advances have been made in our understanding of chimpanzee phylogenetics, little has been known regarding the subspecies composition of captive chimpanzees. We sequenced part of the mtDNA HVR-I region in 218 African-born population founders and performed a phylogenetic analysis with previously characterized African sequences of known provenance to infer subspecies affiliations. Most founders were P.t. verus (95.0%), distantly followed by the troglodytes schweinfurthii clade (4.6%), and a single P.t. vellerosus (0.4%). Pedigree-based estimates of genomic representation in the descendant population revealed that troglodytes schweinfurthii founder representation was reduced in captivity, vellerosus representation increased due to prolific breeding by a single male, and reproductive variance resulted in uneven representation among male P.t.verus founders. No increase in mortality was evident from between-subspecies interbreeding, indicating a lack of outbreeding depression. Knowledge of subspecies and their genomic representation can form the basis for phylogenetically informed genetic management of extant chimpanzees to preserve rare genetic variation for research, conservation, or possible future breeding.

Demographic history and genetic differentiation in apes

Comparisons of genetic variation between humans and great apes are hampered by the fact that we still know little about the demographics and evolutionary history of the latter species. In addition, characterizing ape genetic variation is important because they are threatened with extinction, and knowledge about genetic differentiation among groups may guide conservation efforts. We sequenced multiple intergenic autosomal regions totaling 22,400 base pairs (bp) in ten individuals each from western, central, and eastern chimpanzee groups and in nine bonobos, and 16,000 bp in ten Bornean and six Sumatran orangutans. These regions are analyzed together with homologous information from three human populations and gorillas. We find that whereas orangutans have the highest diversity, western chimpanzees have the lowest, and that the demographic histories of most groups differ drastically. Special attention should therefore be paid to sampling strategies and the statistics chosen when comparing levels of variation within and among groups. Finally, we find that the extent of genetic differentiation among "subspecies" of chimpanzees and orangutans is comparable to that seen among human populations, calling the validity of the "subspecies" concept in apes into question.

Molecular taxonomy in the dark: evolutionary history, phylogeography, and diversity of cave crayfish in the subgenus Aviticambarus, genus Cambarus

Freshwater crayfish species in the subgenus Aviticambarus (Cambaridae: Cambarus) are restricted to caves along the Cumberland Plateau, the Sequatchie Valley, and the Highland Rim which extend along the Tennessee River in southcentral Tennessee and northern Alabama. Historically, three stygobitic species, Cambarus jonesi, Cambarus hamulatus, and Cambarus veitchorum, comprise this subgenus. We examine species' boundaries and phylogeographic structure of this imperiled Southern Appalachian assemblage to shed light on patterns of cave colonization. We also provide estimates of genetic diversity for conservation status assessment. Using geologic evidence, phylogeographic analyses, and sequence data from five gene regions (two nuclear: Histone H3 and GAPDH and three mitochondrial: 12S, 16S, and CO1 totaling almost 2700 base pairs), we identify two well-supported cryptic species in addition to the three currently recognized taxa. Four of these taxa exhibit low levels of genetic variation both currently and historically, which may indicate local extirpation events associated with geological and river basin changes. Our results also support other recent findings that pre-Pleistocene paleodrainages may best explain the current patterns of aquatic faunal biodiversity in the Southern Appalachians.

Subterranean phylogeography of freshwater crayfishes shows extensive gene flow and surprisingly large population sizes

Subterranean animals are currently viewed as highly imperiled, precariously avoiding extinction in an extreme environment of darkness. This assumption is based on a hypothesis that the reduction in visual systems and morphology common in cave faunas reflects a genetic inability to adapt and persist coupled with the perception of a habitat that is limited, disconnected, and fragile. Accordingly, 95% of cave fauna in the United States are presumed endangered due to surface environmental degradation and limited geographic distributions. Our study explores the subterranean phylogeography of stygobitic crayfishes in the southeastern United States, a global hotspot of groundwater biodiversity, using extensive geographic sampling and molecular data. Despite their endangered status, our results show that subterranean crayfish species have attained moderate to high levels of genetic diversity over their evolutionary histories with large population sizes and extensive gene flow among karst systems. We then compare the subterranean population histories to those of common surface stream-dwelling crayfishes. Our results show recent drastic declines in genetic variability in the surface crayfish and suggest that these species also warrant conservation attention.

Contrasting histories of G6PD molecular evolution and malarial resistance in humans and chimpanzees

Although mutations in the glucose-6-phosphate dehydrogenase (G6PD) gene result in several blood-related diseases in humans, they also confer resistance to malarial infection. This association between G6PD and malaria was supported by population genetic analyses of the G6PD locus, which indicated that these mutations may have recently risen in frequency in certain geographic regions as a result of positive selection. Here we characterize nucleotide sequence variation in a 5.2-kb region of the G6PD locus in a population sample of 56 chimpanzees, as well as among 7 other nonhuman primates, to compare with that in humans in determining whether other primates that are impacted by malaria also exhibit patterns of G6PD polymorphism or divergence consistent with positive selection. We find that chimpanzees have several amino acid variants but that the overall pattern at G6PD in chimpanzees, as well as in Old and New World primates in general, can be explained by recent purifying selection as well as strong functional constraint dating back to at least 30-40 MYA. These comparative analyses suggest that the recent signature of positive selection at G6PD in humans is unique.

Hotspots for copy number variation in chimpanzees and humans

Copy number variation is surprisingly common among humans and can be involved in phenotypic diversity and variable susceptibility to complex diseases, but little is known of the extent of copy number variation in nonhuman primates. We have used two array-based comparative genomic hybridization platforms to identify a total of 355 copy number variants (CNVs) in the genomes of 20 wild-born chimpanzees (Pan troglodytes) and have compared the identified chimpanzee CNVs to known human CNVs from previous studies. Many CNVs were observed in the corresponding regions in both chimpanzees and humans; especially those CNVs of higher frequency. Strikingly, these loci are enriched 20-fold for ancestral segmental duplications, which may facilitate CNV formation through nonallelic homologous recombination mechanisms. Therefore, some of these regions may be unstable "hotspots" for the genesis of copy number variation, with recurrent duplications and deletions occurring across and within species.

Independent evolution of bitter-taste sensitivity in humans and chimpanzees

It was reported over 65 years ago that chimpanzees, like humans, vary in taste sensitivity to the bitter compound phenylthiocarbamide (PTC). This was suggested to be the result of a shared balanced polymorphism, defining the first, and now classic, example of the effects of balancing selection in great apes. In humans, variable PTC sensitivity is largely controlled by the segregation of two common alleles at the TAS2R38 locus, which encode receptor variants with different ligand affinities. Here we show that PTC taste sensitivity in chimpanzees is also controlled by two common alleles of TAS2R38; however, neither of these alleles is shared with humans. Instead, a mutation of the initiation codon results in the use of an alternative downstream start codon and production of a truncated receptor variant that fails to respond to PTC in vitro. Association testing of PTC sensitivity in a cohort of captive chimpanzees confirmed that chimpanzee TAS2R38 genotype accurately predicts taster status in vivo. Therefore, although Fisher et al.'s observations were accurate, their explanation was wrong. Humans and chimpanzees share variable taste sensitivity to bitter compounds mediated by PTC receptor variants, but the molecular basis of this variation has arisen twice, independently, in the two species.

GC content evolution of the human and mouse genomes: insights from the study of processed pseudogenes in regions of different recombination rates

Processed pseudogenes are generated by reverse transcription of a functional gene. They are generally nonfunctional after their insertion and, as a consequence, are no longer subjected to the selective constraints associated with functional genes. Because of this property they can be used as neutral markers in molecular evolution. In this work, we investigated the relationship between the evolution of GC content in recently inserted processed pseudogenes and the local recombination pattern in two mammalian genomes (human and mouse). We confirmed, using original markers, that recombination drives GC content in the human genome and we demonstrated that this is also true for the mouse genome despite lower recombination rates. Finally, we discussed the consequences on isochores evolution and the contrast between the human and the mouse pattern.

The chimpanzee-specific pericentric inversions that distinguish humans and chimpanzees have identical breakpoints in Pan troglodytes and Pan paniscus

Seven of nine pericentric inversions that distinguish human (HSA) and chimpanzee karyotypes are chimpanzee-specific. In this study we investigated whether the two extant chimpanzee species, Pan troglodytes (common chimpanzee) and Pan paniscus (bonobo), share exactly the same pericentric inversions. The methods applied were FISH with breakpoint-spanning BAC/PAC clones and PCR analyses of the breakpoint junction sequences. Our findings for the homologues to HSA 4, 5, 9, 12, 16, and 17 confirm for the first time at the sequence level that these pericentric inversions have identical breakpoints in the common chimpanzee and the bonobo. Therefore, these inversions predate the separation of the two chimpanzee species 0.86-2 Mya. Further, the inversions distinguishing human and chimpanzee karyotypes may be regarded as early acquisitions, such that they are likely to have been present at the time of human/chimpanzee divergence. According to the chromosomal speciation theory the inversions themselves could have promoted human speciation.

Comparing the human and chimpanzee genomes: searching for needles in a haystack

The chimpanzee genome sequence is a long-awaited milestone, providing opportunities to explore primate evolution and genetic contributions to human physiology and disease. Humans and chimpanzees shared a common ancestor approximately 5-7 million years ago (Mya). The difference between the two genomes is actually not approximately 1%, but approximately 4%--comprising approximately 35 million single nucleotide differences and approximately 90 Mb of insertions and deletions. The challenge is to identify the many evolutionarily, physiologically, and biomedically important differences scattered throughout these genomes while integrating these data with emerging knowledge about the corresponding "phenomes" and the relevant environmental influences. It is logical to tackle the genetic aspects via both genome-wide analyses and candidate gene studies. Genome-wide surveys could eliminate the majority of genomic sequence differences from consideration, while simultaneously identifying potential targets of opportunity. Meanwhile, candidate gene approaches can be based on such genomic surveys, on genes that may contribute to known differences in phenotypes or disease incidence/severity, or on mutations in the human population that impact unique aspects of the human condition. These two approaches will intersect at many levels and should be considered complementary. We also cite some known genetic differences between humans and great apes, realizing that these likely represent only the tip of the iceberg.

Population history, biogeography, and taxonomy of orangutans (Genus: Pongo) based on a population genetic meta-analysis of multiple loci

This paper examines orangutan population history and evolution through a meta-analysis of seven loci collected from both Sumatran and Bornean orangutans. Within orangutans, most loci show that the Sumatran population is about twice as diverse as the Bornean population. Orangutans are more diverse than African apes and humans. Sumatran and Bornean populations show significant genetic differentiation from one another and their history does not differ significantly from an 'island model' (population splitting without gene flow). Two different methods support a divergence of Bornean and Sumatran orangutans at 2.7-5 million years ago. This suggests that Pleistocene events, such as the cyclical exposure of the Sunda shelf and the Toba volcanic eruption, did not have a major impact on the divergence of Bornean and Sumatran orangutans. Pairwise mismatch analyses, however, suggest that Bornean orangutans have undergone a recent population expansion (beginning 39,000-64,000 years ago), while Sumatran orangutan populations were stable. Pleistocene events may have contributed to these aspects of orangutan population history. These conclusions are applied to the debate on orangutan taxonomy.

Variable molecular clocks in hominoids

Generation time is an important determinant of a neutral molecular clock. There are several human-specific life history traits that led to a substantially longer generation time in humans than in other hominoids. Indeed, a long generation time is considered an important trait that distinguishes humans from their closest relatives. Therefore, humans may exhibit a significantly slower molecular clock as compared to other hominoids. To investigate this hypothesis, we performed a large-scale analysis of lineage-specific rates of single-nucleotide substitutions among hominoids. We found that humans indeed exhibit a significant slowdown of molecular evolution compared to chimpanzees and other hominoids. However, the amount of fixed differences between humans and chimpanzees appears extremely small, suggesting a very recent evolution of human-specific life history traits. Notably, chimpanzees also exhibit a slower rate of molecular evolution compared to gorillas and orangutans in the regions analyzed.

Should the draft chimpanzee sequence be finished?

Owing to the availability of genome working drafts (WDs), current comparative-sequence studies are frequently performed on a genome-wide scale. In this article, we appraise the utility of WD sequences in the detection of genomic differences in closely related species. We compared human DNA sequences with draft and high-quality versions of the corresponding chimpanzee loci to reveal the overall high quality of the chimp WD sequence. Nevertheless, a significant proportion of the differences between WD and high-quality sequences we observed can be attributed to sequencing errors in the draft. Although we suggest methods to reduce the number of such false positives efficiently, our study emphasizes the benefit expected from finishing the chimpanzee genome sequence.

Polymorphic micro-inversions contribute to the genomic variability of humans and chimpanzees

A combination of inter- and intra-species genome comparisons is required to identify and classify the full spectrum of genetic changes, both subtle and gross, that have accompanied the evolutionary divergence of humans and other primates. In this study, gene order comparisons of 11,518 human and chimpanzee orthologous gene pairs were performed to detect regions of inverted gene order that are potentially indicative of small-scale rearrangements such as inversions. By these means, a total of 71 potential micro-rearrangements were detected, nine of which were considered to represent micro-inversions encompassing more than three genes. These putative inversions were then investigated by FISH and/or PCR analyses and the authenticity of five of the nine inversions, ranging in size from approximately 800 kb to approximately 4.4 Mb, was confirmed. These inversions mapped to 1p13.2-13.3, 7p22.1, 7p13-14.1, 18p11.21-11.22 and 19q13.12 and encompass 50, 14, 16, 7 and 16 known genes, respectively. Intriguingly, four of the confirmed inversions turned out to be polymorphic: three were polymorphic in the chimpanzee and one in humans. It is concluded that micro-inversions make a significant contribution to genomic variability in both humans and chimpanzees and inversion polymorphisms may be more frequent than previously realized.

Selection in favor of nucleotides G and C diversifies evolution rates and levels of polymorphism at mammalian synonymous sites

The impact of synonymous nucleotide substitutions on fitness in mammals remains controversial. Despite some indications of selective constraint, synonymous sites are often assumed to be neutral, and the rate of their evolution is used as a proxy for mutation rate. We subdivide all sites into four classes in terms of the mutable CpG context, nonCpG, postC, preG, and postCpreG, and compare four-fold synonymous sites and intron sites residing outside transposable elements. The distribution of the rate of evolution across all synonymous sites is trimodal. Rate of evolution at nonCpG synonymous sites, not preceded by C and not followed by G, is approximately 10% below that at such intron sites. In contrast, rate of evolution at postCpreG synonymous sites is approximately 30% above that at such intron sites. Finally, synonymous and intron postC and preG sites evolve at similar rates. The relationship between the levels of polymorphism at the corresponding synonymous and intron sites is very similar to that between their rates of evolution. Within every class, synonymous sites are occupied by G or C much more often than intron sites, whose nucleotide composition is consistent with neutral mutation-drift equilibrium. These patterns suggest that synonymous sites are under weak selection in favor of G and C, with the average coefficient s approximately 0.25/Ne approximately 10(-5), where Ne is the effective population size. Such selection decelerates evolution and reduces variability at sites with symmetric mutation, but has the opposite effects at sites where the favored nucleotides are more mutable. The amino-acid composition of proteins dictates that many synonymous sites are CpGprone, which causes them, on average, to evolve faster and to be more polymorphic than intron sites. An average genotype carries approximately 10(7) suboptimal nucleotides at synonymous sites, implying synergistic epistasis in selection against them.

Strong and weak male mutation bias at different sites in the primate genomes: insights from the human-chimpanzee comparison

Male mutation bias is a higher mutation rate in males than in females thought to result from the greater number of germ line cell divisions in males. If errors in DNA replication cause most mutations, then the magnitude of male mutation bias, measured as the male-to-female mutation rate ratio (alpha), should reflect the relative excess of male versus female germ line cell divisions. Evolutionary rates averaged among all sites in a sequence and compared between mammalian sex chromosomes were shown to be indeed higher in males than in females. However, it is presently unknown whether individual classes of substitutions exhibit such bias. To address this issue, we investigated male mutation bias separately at non-CpG and CpG sites using human-chimpanzee whole-genome alignments. We observed strong male mutation bias at non-CpG sites: alpha in the X-autosome comparison was approximately 6-7, which was similar to the male-to-female ratio in the number of germ line cell divisions. In contrast, mutations at CpG sites exhibited weak male mutation bias: alpha in the X-autosome comparison was only approximately 2-3. This is consistent with the methylation-induced and replication-independent mechanism of CpG transitions, which constitute the majority of mutations at CpG sites. Interestingly, our study also indicated weak male mutation bias for transversions at CpG sites, implying a spontaneous mechanism largely not associated with replication. Male mutation bias was equally strong at CpG and non-CpG sites located within unmethylated "CpG islands," suggesting the replication-dependent origin of these mutations. Thus, we found that the strength of male mutation bias is nonuniform in the primate genomes. Importantly, we discovered that male mutation bias depends on the proportion of CpG sites in the loci compared. This might explain the differences in the magnitude of primate male mutation bias observed among studies.

Role of selection in fixation of gene duplications

New genes commonly appear through complete or partial duplications of pre-existing genes. Duplications of long DNA segments are constantly produced by rare mutations, may become fixed in a population by selection or random drift, and are subject to divergent evolution of the paralogous sequences after fixation, although gene conversion can impede this process. New data shed some light on each of these processes. Mutations which involve duplications can occur through at least two different mechanisms, backward strand slippage during DNA replication and unequal crossing-over. The background rate of duplication of a complete gene in humans is 10(-9)-10(-10) per generation, although many genes located within hot-spots of large-scale mutation are duplicated much more often. Many gene duplications affect fitness strongly, and are responsible, through gene dosage effects, for a number of genetic diseases. However, high levels of intrapopulation polymorphism caused by presence or absence of long, gene-containing DNA segments imply that some duplications are not under strong selection. The polymorphism to fixation ratios appear to be approximately the same for gene duplications and for presumably selectively neutral nucleotide substitutions, which, according to the McDonald-Kreitman test, is consistent with selective neutrality of duplications. However, this pattern can also be due to negative selection against most of segregating duplications and positive selection for at least some duplications which become fixed. Patterns in post-fixation evolution of duplicated genes do not easily reveal the causes of fixations. Many gene duplications which became fixed recently in a variety of organisms were positively selected because the increased expression of the corresponding genes was beneficial. The effects of gene dosage provide a unified framework for studying all phases of the life history of a gene duplication. Application of well-known methods of evolutionary genetics to accumulating data on new, polymorphic, and fixed duplication will enhance our understanding of the role of natural selection in the evolution by gene duplication.

The use of racial, ethnic, and ancestral categories in human genetics research

The global dispersal of anatomically modern humans over the past 100,000 years has produced patterns of phenotypic variation that have exerted--and continue to exert--powerful influences on the lives of individuals and the experiences of groups. The recency of our common ancestry and continued gene flow among populations have resulted in less genetic differentiation among geographically distributed human populations than is observed in many other mammalian species. Nevertheless, differences in appearance have contributed to the development of ideas about "race" and "ethnicity" that often include the belief that significant inherited differences distinguish humans. The use of racial, ethnic, and ancestral categories in genetics research can imply that group differences arise directly through differing allele frequencies, with little influence from socially mediated mechanisms. At the same time, careful investigations of the biological, environmental, social, and psychological attributes associated with these categories will be an essential component of cross-disciplinary research into the origins, prevention, and treatment of common diseases, including those diseases that differ in prevalence among groups.

Distribution of the strength of selection against amino acid replacements in human proteins

The impact of an amino acid replacement on the organism's fitness can vary from lethal to selectively neutral and even, in rare cases, beneficial. Substantial data are available on either pathogenic or acceptable replacements. However, the whole distribution of coefficients of selection against individual replacements is not known for any organism. To ascertain this distribution for human proteins, we combined data on pathogenic missense mutations, on human non-synonymous SNPs and on human-chimpanzee divergence of orthologous proteins. Fractions of amino acid replacements which reduce fitness by >10(-2), 10(-2)-10(-4), 10(-4)-10(-5) and <10(-5) are 25, 49, 14 and 12%, respectively. On average, the strength of selection against a replacement is substantially higher when chemically dissimilar amino acids are involved, and the Grantham's index of a replacement explains 35% of variance in the average logarithm of selection coefficients associated with different replacements. Still, the impact of a replacement depends on its context within the protein more than on its own nature. Reciprocal replacements are often associated with rather different selection coefficients, in particular, replacements of non-polar amino acids with polar ones are typically much more deleterious than replacements in the opposite direction. However, differences between evolutionary fluxes of reciprocal replacements are only weakly correlated with the differences between the corresponding selection coefficients.

Initial sequence of the chimpanzee genome and comparison with the human genome

Here we present a draft genome sequence of the common chimpanzee (Pan troglodytes). Through comparison with the human genome, we have generated a largely complete catalogue of the genetic differences that have accumulated since the human and chimpanzee species diverged from our common ancestor, constituting approximately thirty-five million single-nucleotide changes, five million insertion/deletion events, and various chromosomal rearrangements. We use this catalogue to explore the magnitude and regional variation of mutational forces shaping these two genomes, and the strength of positive and negative selection acting on their genes. In particular, we find that the patterns of evolution in human and chimpanzee protein-coding genes are highly correlated and dominated by the fixation of neutral and slightly deleterious alleles. We also use the chimpanzee genome as an outgroup to investigate human population genetics and identify signatures of selective sweeps in recent human evolution.

Alu retrotransposition-mediated deletion

Alu repeats contribute to genomic instability in primates via insertional and recombinational mutagenesis. Here, we report an analysis of Alu element-induced genomic instability through a novel mechanism termed retrotransposition-mediated deletion, and assess its impact on the integrity of primate genomes. For human and chimpanzee genomes, we find evidence of 33 retrotransposition-mediated deletion events that have eliminated approximately 9000 nucleotides of genomic DNA. Our data suggest that, during the course of primate evolution, Alu retrotransposition may have contributed to over 3000 deletion events, eliminating approximately 900 kb of DNA in the process. Potential mechanisms for the creation of Alu retrotransposition-mediated deletions include L1 endonuclease-dependent retrotransposition, L1 endonuclease-independent retrotransposition, internal priming on DNA breaks, and promiscuous target primed reverse transcription. A comprehensive analysis of the collateral effects by Alu mobilization on all primate genomes will require sequenced genomes from representatives of the entire order.

A structured ancestral population for the evolution of modern humans

The view that modern humans evolved through a bottleneck from a single founding group of archaic Homo is being challenged by new analyses of contemporary genetic variation. A wide range of middle to late Pleistocene ages for gene genealogies and evidence for early population structures point to a diverse and scattered ancestry associated with a metapopulation history of local extinctions, re-colonization and admixture. A different balance of the same processes has shaped chimpanzee diversity.

Nuclear insertions help and hinder inference of the evolutionary history of gorilla mtDNA

Numts are fragments of mitochondrial DNA (mtDNA) that have been translocated to the nucleus, where they can persist while their mitochondrial counterparts continue to rapidly evolve. Thus, numts represent 'molecular fossils' useful for comparison with mitochondrial variation, and are particularly suited for studies of the fast-evolving hypervariable segment of the mitochondrial control region (HV1). Here we used information from numts found in western gorillas (Gorilla gorilla) and eastern gorillas (Gorilla beringei) to estimate that these two species diverged about 1.3 million years ago (Ma), an estimate similar to recent calculations for the divergence of chimpanzee and bonobo. We also describe the sequence of a gorilla numt still possessing a segment lost from all contemporary gorilla mtDNAs. In contrast to that sequence, many numts of the HV1 are highly similar to authentic mitochondrial organellar sequences, making it difficult to determine whether purported mitochondrial sequences truly derive from that genome. We used all available organellar HV1 and corresponding numt sequences from gorillas in a phylogenetic analysis aimed at distinguishing these two types of sequences. Numts were found in several clades in the tree. This, in combination with the fact that only a limited amount of the extant variation in gorillas has been sampled, suggests that categorization of new sequences by the indirect means of phylogenetic comparison would be prone to uncertainty. We conclude that for taxa such as gorillas that contain numerous numts, direct approaches to the authentication of HV1 sequences, such as amplification strategies relying upon the circularity of the mtDNA molecule, remain necessary.

Reduced MIC gene repertoire variation in West African chimpanzees as compared to humans

The human major histocompatibility complex class I chain-related (MIC) genes are members of a multicopy family showing similarity to the classical HLA-A, HLA-B, and HLA-C genes. Only the MICA and MICB genes produce functional transcripts. In chimpanzees, however, only one MIC gene is expressed, showing an intermediate character, resulting from a deletion fusing the MICA and MICB gene segments together. The present population study illustrates that all chimpanzee haplotypes sampled possess the hybrid MICA/B gene. In contrast to the human situation this gene displays reduced allelic variation. The observed repertoire reduction of the chimpanzee MICA/B gene is in conformity with the severe repertoire condensation documented for Patr-B locus lineages, probably due to the close proximity of both genes.

Rivers influence the population genetic structure of bonobos (Pan paniscus)

Bonobos are large, highly mobile primates living in the relatively undisturbed, contiguous forest south of the Congo River. Accordingly, gene flow among populations is assumed to be extensive, but may be impeded by large, impassable rivers. We examined mitochondrial DNA control region sequence variation in individuals from five distinct localities separated by rivers in order to estimate relative levels of genetic diversity and assess the extent and pattern of population genetic structure in the bonobo. Diversity estimates for the bonobo exceed those for humans, but are less than those found for the chimpanzee. All regions sampled are significantly differentiated from one another, according to genetic distances estimated as pairwise FSTs, with the greatest differentiation existing between region East and each of the two Northern populations (N and NE) and the least differentiation between regions Central and South. The distribution of nucleotide diversity shows a clear signal of population structure, with some 30% of the variance occurring among geographical regions. However, a geographical patterning of the population structure is not obvious. Namely, mitochondrial haplotypes were shared among all regions excepting the most eastern locality and the phylogenetic analysis revealed a tree in which haplotypes were intermixed with little regard to geographical origin, with the notable exception of the close relationships among the haplotypes found in the east. Nonetheless, genetic distances correlated with geographical distances when the intervening distances were measured around rivers presenting effective current-day barriers, but not when straight-line distances were used, suggesting that rivers are indeed a hindrance to gene flow in this species.