Low nucleotide diversity in chimpanzees and bonobos

Morphological variation of the Pan talus relative to that of Gorilla

OBJECTIVES

Differences in talar articular morphology relative to locomotion have recently been found within Pan and Gorilla. Whole-bone talar morphology within, and shared variation among, Pan and Gorilla (sub)species, however, has yet to be investigated. Here we separately analyze talar external shape within Pan (P. t. troglodytes, P. t. schweinfurthii, P. t. verus, P. paniscus) and Gorilla (G. g. gorilla, G. b. beringei, G. b. graueri) relative to degree of arboreality and body size. Pan and Gorilla are additionally analyzed together to determine if consistent shape differences exist within the genera.

MATERIALS AND METHODS

Talar external shape was quantified using a weighted spherical harmonic analysis. Shape variation both within and among Pan and Gorilla was described using principal component analyses. Root mean square distances were calculated between taxon averages, and resampling statistics conducted to test for pairwise differences.

RESULTS

P. t. verus (most arboreal Pan) talar shape significantly differs from other Pan taxa (p < 0.05 for pairwise comparisons) driven by more asymmetrical trochlear rims and a medially-set talar head. P. t. troglodytes, P. t. schweinfurthii, and P. paniscus do not significantly differ (p > 0.05 for pairwise comparisons). All gorilla taxa exhibit significantly different talar morphologies (p < 0.007 for pairwise comparisons). The more terrestrial subspecies of G. beringei and P. troglodytes exhibit a superoinferiorly taller talar head/neck complex.

DISCUSSION

P. t. verus exhibits talar morphologies that have been previously related to more frequent arboreality. The adaptations in the more terrestrial G. beringei and P. troglodytes subspecies may serve to facilitate load transmission.

Carpal allometry of African apes among mammals

OBJECTIVES

Morphological variation in African ape carpals has been used to support the idea that Pan and Gorilla evolved knuckle-walking independently. Little work, however, has focused on the effect of body mass on carpal morphology. Here, we compare carpal allometry in Pan and Gorilla to that of other quadrupedal mammals with similar body mass differences. If allometric trends in Pan and Gorilla carpals mirror those of other mammals with similar body mass variation, then body mass differences may provide a more parsimonious explanation for African ape carpal variation than the independent evolution of knuckle-walking.

MATERIALS AND METHODS

Three linear measurements were collected on the capitate, hamate, lunate, and scaphoid (or scapholunate) of 39 quadrupedal species from six mammalian families/subfamilies. Relationships between linear measurements and estimated body mass were analyzed using reduced major axis regression. Slopes were compared to 0.33 for isometry.

RESULTS

Within Hominidae, higher body mass taxa (Gorilla) have relatively anteroposteriorly wider, mediolaterally wider, and/or proximodistally shorter capitates, hamates, and scaphoids than low body mass taxa (Pan). These allometric relationships are mirrored in most, but not all, mammalian families/subfamilies included in the analysis.

CONCLUSIONS

Within most mammalian families/subfamilies, carpals of high body mass taxa are proximodistally shorter, anteroposteriorly wider, and mediolaterally wider than those of low body mass taxa. These distinctions may be caused by the need to accommodate relatively higher forelimb loading associated with greater body mass. Because these trends occur within multiple mammalian families/subfamilies, some carpal variation in Pan and Gorilla is consistent with body mass differences.

Ampliconic Genes on the Great Ape Y Chromosomes: Rapid Evolution of Copy Number but Conservation of Expression Levels

Multicopy ampliconic gene families on the Y chromosome play an important role in spermatogenesis. Thus, studying their genetic variation in endangered great ape species is critical. We estimated the sizes (copy number) of nine Y ampliconic gene families in population samples of chimpanzee, bonobo, and orangutan with droplet digital polymerase chain reaction, combined these estimates with published data for human and gorilla, and produced genome-wide testis gene expression data for great apes. Analyzing this comprehensive data set within an evolutionary framework, we, first, found high inter- and intraspecific variation in gene family size, with larger families exhibiting higher variation as compared with smaller families, a pattern consistent with random genetic drift. Second, for four gene families, we observed significant interspecific size differences, sometimes even between sister species—chimpanzee and bonobo. Third, despite substantial variation in copy number, Y ampliconic gene families’ expression levels did not differ significantly among species, suggesting dosage regulation. Fourth, for three gene families, size was positively correlated with gene expression levels across species, suggesting that, given sufficient evolutionary time, copy number influences gene expression. Our results indicate high variability in size but conservation in gene expression levels in Y ampliconic gene families, significantly advancing our understanding of Y-chromosome evolution in great apes.

The population genetics of nonmigratory Allen's Hummingbird (Selasphorus sasin sedentarius) following a recent mainland colonization

Allen's Hummingbird comprises two subspecies, one migratory (Selasphorus sasin sasin) and one nonmigratory (S. s. sedentarius). The nonmigratory subspecies, previously endemic to the California Channel Islands, apparently colonized the California mainland on the Palos Verdes Peninsula some time before 1970 and now breeds throughout coastal southern California. We sequenced and compared populations of mainland nonmigratory Allen's Hummingbird to Channel Island populations from Santa Catalina, San Clemente, and Santa Cruz Island. We found no evidence of founder effects on the mainland population. Values of nucleotide diversity on the mainland were higher than on the Channel Islands. There were low levels of divergence between the Channel Islands and the mainland, and Santa Cruz Island was the most genetically distinct. Ecological niche models showed that rainfall and temperature variables on the Channel Islands are similar in the Los Angeles basin and predicted continued expansion of nonmigratory Allen's Hummingbird north along the coast and inland. We also reviewed previous genetic studies of vertebrate species found on the Channel Islands and mainland and showed that broad conclusions regarding island-mainland patterns remain elusive. Challenges include the idiosyncratic nature of colonization itself as well as the lack of a comprehensive approach that incorporates similar markers and sampling strategies across taxa, which, within the context of a comparative study of island-mainland relationships, may lead to inconsistent results.

Similar patterns of genetic diversity and linkage disequilibrium in Western chimpanzees (Pan troglodytes verus) and humans indicate highly conserved mechanisms of MHC molecular evolution

BACKGROUND

Many species are threatened with extinction as their population sizes decrease with changing environments or face novel pathogenic threats. A reduction of genetic diversity at major histocompatibility complex (MHC) genes may have dramatic effects on populations' survival, as these genes play a key role in adaptive immunity. This might be the case for chimpanzees, the MHC genes of which reveal signatures of an ancient selective sweep likely due to a viral epidemic that reduced their population size a few million years ago. To better assess how this past event affected MHC variation in chimpanzees compared to humans, we analysed several indexes of genetic diversity and linkage disequilibrium across seven MHC genes on four cohorts of chimpanzees and we compared them to those estimated at orthologous HLA genes in a large set of human populations.

RESULTS

Interestingly, the analyses uncovered similar patterns of both molecular diversity and linkage disequilibrium across the seven MHC genes in chimpanzees and humans. Indeed, in both species the greatest allelic richness and heterozygosity were found at loci A, B, C and DRB1, the greatest nucleotide diversity at loci DRB1, DQA1 and DQB1, and both significant global linkage disequilibrium and the greatest proportions of haplotypes in linkage disequilibrium were observed at pairs DQA1 ~ DQB1, DQA1 ~ DRB1, DQB1 ~ DRB1 and B ~ C. Our results also showed that, despite some differences among loci, the levels of genetic diversity and linkage disequilibrium observed in contemporary chimpanzees were globally similar to those estimated in small isolated human populations, in contrast to significant differences compared to large populations.

CONCLUSIONS

We conclude, first, that highly conserved mechanisms shaped the diversity of orthologous MHC genes in chimpanzees and humans. Furthermore, our findings support the hypothesis that an ancient demographic decline affecting the chimpanzee populations - like that ascribed to a viral epidemic - exerted a substantial effect on the molecular diversity of their MHC genes, albeit not more pronounced than that experienced by HLA genes in human populations that underwent rapid genetic drift during humans' peopling history. We thus propose a model where chimpanzees' MHC genes regenerated molecular variation through recombination/gene conversion and/or balancing selection after the selective sweep.

Adapting Biased Gene Conversion theory to account for intensive GC-content deterioration in the human genome by novel mutations

We examined seventy million well-characterized human mutations, and their impact on G+C-compositional dynamics, in order to understand the formation and maintenance of major genomic nucleotide sequence patterns. Among novel mutations, those that change a strong (S) base pair G:C/C:G to a weak (W) pair A:T/T:A occur at nearly twice the frequency of the opposite mutations. Such imbalance puts strong downward pressure on overall GC-content. However, along protracted paths to fixation, S→W mutations are much less likely to propagate than W→S mutations. The magnitude of relative propagation disadvantages for S→W mutations is inexplicable by any currently-accepted model. This fact forced us to re-examine the quantitative features of Biased Gene Conversion (BGC) theory. Revised parameters of BGC that, per average individual, convert 7–14 W base pairs into S pairs, would account for the S-content turnover differences between new and old mutations, and make BGC an instrumental force for nucleotide dynamics and evolution. BGC should thus be considered seriously in both theories and biomedical practice. In particular, BGC should be taken into account during allele imputations, where missing SNP alleles are computationally predicted based on the information about several neighboring alleles. Finally, we analyzed the effect of neighboring nucleotide context on the mutation frequencies, dynamics, and GC-composition turnover. For this purpose, we examined genomic regions having extremely biased nucleotide compositions (enriched for S-, W-, purine/pyrimidine strand asymmetry, or AC/GT-strand asymmetry). It was found that point mutations in these regions preferentially degrade the nucleotide inhomogeneities, decreasing the sequence biases. Degradation of sequence bias is highest for novel mutations, and considerably lower for older mutations (those widespread across populations). Besides BGC, there may be additional, still uncharacterized molecular mechanisms that either preserve genomic regions with biased nucleotide compositions from mutational degradation or fail to degrade such inhomogeneities in specific chromosomal regions.

Evolution of the Human Nervous System Function, Structure, and Development

The nervous system-in particular, the brain and its cognitive abilities-is among humans' most distinctive and impressive attributes. How the nervous system has changed in the human lineage and how it differs from that of closely related primates is not well understood. Here, we consider recent comparative analyses of extant species that are uncovering new evidence for evolutionary changes in the size and the number of neurons in the human nervous system, as well as the cellular and molecular reorganization of its neural circuits. We also discuss the developmental mechanisms and underlying genetic and molecular changes that generate these structural and functional differences. As relevant new information and tools materialize at an unprecedented pace, the field is now ripe for systematic and functionally relevant studies of the development and evolution of human nervous system specializations.

Genetic diversity in humans and non-human primates and its evolutionary consequences

Genetic diversity is a key parameter in population genetics and is important for understanding the process of evolution and for the development of appropriate conservation strategies. Recent advances in sequencing technology have enabled the measurement of genetic diversity of various organisms at the nucleotide level and on a genome-wide scale, yielding more precise estimates than were previously achievable. In this review, I have compiled and summarized the estimates of genetic diversity in humans and non-human primates based on recent genome-wide studies. Although studies on population genetics demonstrated fluctuations in population sizes over time, general patterns have emerged. As shown previously, genetic diversity in humans is one of the lowest among primates; however, certain other primate species exhibit genetic diversity that is comparable to or even lower than that in humans. There exists greater than 10-fold variation in genetic diversity among primate species, and I found weak correlation with species fecundity but not with body or propagule size. I further discuss the potential evolutionary consequences of population size decline on the evolution of primate species. The level of genetic diversity negatively correlates with the ratio of non-synonymous to synonymous polymorphisms in a population, suggesting that proportionally greater numbers of slightly deleterious mutations segregate in small rather than large populations. Although population size decline is likely to promote the fixation of slightly deleterious mutations, there are molecular mechanisms, such as compensatory mutations at various molecular levels, which may prevent fitness decline at the population level. The effects of slightly deleterious mutations from theoretical and empirical studies and their relevance to conservation biology are also discussed in this review.

How did bonobos come to range south of the congo river? Reconsideration of the divergence of Pan paniscus from other Pan populations

While investigating the genetic structure in wild bonobos,(1) we realized that the widely accepted scenario positing that the Pleistocene appearance of the Congo River separated the common ancestor of chimpanzees (Pan troglodytes) and bonobos (P. paniscus) into two species is not supported by recent geographical knowledge about the formation of the Congo River. We explored the origin of bonobos using a broader biogeographical perspective by examining local faunas in the central African region. The submarine Congo River sediments and paleotopography of central Africa show that the Congo River has functioned as a geographical barrier for the last 34 million years. This evidence allows us to hypothesize that when the river was first formed, the ancestor of bonobos did not inhabit the current range of the species on the left bank of the Congo River but that, during rare times when the Congo River discharge decreased during the Pleistocene, one or more founder populations of ancestral Pan paniscus crossed the river to its left bank. The proposed scenario for formation of the Congo River and the corridor hypothesis for an ancestral bonobo population is key to understanding the distribution of great apes and their evolution.

IMa2p--parallel MCMC and inference of ancient demography under the Isolation with migration (IM) model

IMa2 and related programs are used to study the divergence of closely related species and of populations within species. These methods are based on the sampling of genealogies using MCMC, and they can proceed quite slowly for larger data sets. We describe a parallel implementation, called IMa2p, that provides a nearly linear increase in genealogy sampling rate with the number of processors in use. IMa2p is written in OpenMPI and C++, and scales well for demographic analyses of a large number of loci and populations, which are difficult to study using the serial version of the program.

Neutral nuclear variation in Baboons (genus Papio) provides insights into their evolutionary and demographic histories

Baboons (genus Papio) are distributed over most of sub-Saharan Africa and in the southern portion of the Arabian Peninsula. Six distinct morphotypes, with clearly defined geographic distributions, are recognized (the olive, chacma, yellow, Guinea, Kinda, and hamadryas baboons). The evolutionary relationships among baboon forms have long been a controversial issue. Phylogenetic analyses based on mitochondrial DNA sequences revealed that the modern baboon morphotypes are mitochondrially paraphyletic or polyphyletic. The discordance between mitochondrial lineages and morphology is indicative of extensive introgressive hybridization between ancestral baboon populations. To gain insights into the evolutionary relationships among morphotypes and their demographic history, we performed an analysis of nuclear variation in baboons. We sequenced 13 noncoding, putatively neutral, nuclear regions, and scored the presence/absence of 18 polymorphic transposable elements in a sample of 45 baboons belonging to five of the six recognized baboon forms. We found that the chacma baboon is the sister-taxon to all other baboons and the yellow baboon is the sister-taxon to an unresolved northern clade containing the olive, Guinea, and hamadryas baboons. We estimated that the diversification of baboons occurred entirely in the Pleistocene, the earliest split dating ∼1.5 million years ago, and that baboons have experienced relatively large and constant effective population sizes for most of their evolutionary history (∼30,000 to 95,000 individuals).

Speciation, population structure, and demographic history of the Mojave Fringe-toed Lizard (Uma scoparia), a species of conservation concern

The North American deserts were impacted by both Neogene plate tectonics and Quaternary climatic fluctuations, yet it remains unclear how these events influenced speciation in this region. We tested published hypotheses regarding the timing and mode of speciation, population structure, and demographic history of the Mojave Fringe-toed Lizard (Uma scoparia), a sand dune specialist endemic to the Mojave Desert of California and Arizona. We sampled 109 individual lizards representing 22 insular dune localities, obtained DNA sequences for 14 nuclear loci, and found that U. scoparia has low genetic diversity relative to the U. notata species complex, comparable to that of chimpanzees and southern elephant seals. Analyses of genotypes using Bayesian clustering algorithms did not identify discrete populations within U. scoparia. Using isolation-with-migration (IM) models and a novel coalescent-based hypothesis testing approach, we estimated that U. scoparia diverged from U. notata in the Pleistocene epoch. The likelihood ratio test and the Akaike Information Criterion consistently rejected nested speciation models that included parameters for migration and population growth of U. scoparia. We reject the Neogene vicariance hypothesis for the speciation of U. scoparia and define this species as a single evolutionarily significant unit for conservation purposes.

The apolipoprotein E (APOE) gene appears functionally monomorphic in chimpanzees (Pan troglodytes)

Background

The human apolipoprotein E (APOE) gene is polymorphic, with three primary alleles (E2, E3, E4) that differ at two key non-synonymous sites. These alleles are functionally different in how they bind to lipoproteins, and this genetic variation is associated with phenotypic variation for several medical traits, including cholesterol levels, cardiovascular health, Alzheimer’s disease risk, and longevity. The relative frequencies of these alleles vary across human populations, and the evolution and maintenance of this diversity is much debated. Previous studies comparing human and chimpanzee APOE sequences found that the chimpanzee sequence is most similar to the human E4 allele, although the resulting chimpanzee protein might function like the protein coded for by the human E3 allele. However, these studies have used sequence data from a single chimpanzee and do not consider whether chimpanzees, like humans, show intra-specific and subspecific variation at this locus.

Methodology and Principal Findings

To examine potential intraspecific variation, we sequenced the APOE gene of 32 chimpanzees. This sample included 20 captive individuals representing the western subspecies (P. troglodytes verus) and 12 wild individuals representing the eastern subspecies (P. t. schweinfurthii). Variation in our resulting sequences was limited to one non-coding, intronic SNP, which showed fixed differences between the two subspecies. We also compared APOE sequences for all available ape genera and fossil hominins. The bonobo APOE protein is identical to that of the chimpanzee, and the Denisovan APOE exhibits all four human-specific, non-synonymous changes and appears functionally similar to the human E4 allele.

Conclusions

We found no coding variation within and between chimpanzee populations, suggesting that the maintenance of functionally diverse APOE polymorphisms is a unique feature of human evolution.

Contrasting patterns of X/Y polymorphism distinguish Carica papaya from other sex chromosome systems

The sex chromosomes of the tropical crop papaya (Carica papaya) are evolutionarily young and consequently allow for the examination of evolutionary mechanisms that drive early sex chromosome divergence. We conducted a molecular population genetic analysis of four X/Y gene pairs from a collection of 45 wild papaya accessions. These population genetic analyses reveal striking differences in the patterns of polymorphism between the X and Y chromosomes that distinguish them from other sex chromosome systems. In most sex chromosome systems, the Y chromosome displays significantly reduced polymorphism levels, whereas the X chromosome maintains a level of polymorphism that is comparable to autosomal loci. However, the four papaya sex-linked loci that we examined display diversity patterns that are opposite this trend: the papaya X alleles exhibit significantly reduced polymorphism levels, whereas the papaya Y alleles maintain greater than expected levels of diversity. Our analyses suggest that selective sweeps in the regions of the X have contributed to this pattern while also revealing geographically restricted haplogroups on the Y. We discuss the possible role sexual selection and/or genomic conflict have played in shaping the contrasting patterns of polymorphism found for the papaya X and Y chromosomes.

Genetic analysis of cytomegalovirus in malignant gliomas

Human cytomegalovirus (HCMV) has been found in malignant gliomas at variable frequencies with efforts to date focused on characterizing the role(s) of single gene products in disease. Here, we reexamined the HCMV prevalence in malignant gliomas using different methods and began to dissect the genetics of HCMV in tumors. HCMV DNA was found in 16/17 (94%) tumor specimens. Viral DNA copy numbers were found to be low and variable, ranging from 10(2) to 10(6) copies/500 ng of total DNA. The tumor tissues had incongruences between viral DNA copy numbers and protein levels. However, nonlatent protein expression was detected in many tumors. The viral UL83 gene, encoding pp65, was found to segregate into five cancer-associated genotypes with a bias for amino acid changes in glioblastoma multiforme (GBM) in comparison to the low-grade tumors. Deep sequencing of a GBM-associated viral population resulted in 81,224 bp of genome coverage. Sequence analysis revealed the presence of intact open reading frames and higher numbers of high-frequency variations within the repeat long region compared to the unique long region, which harbors many core genes, and the unique short region (P = 0.001). This observation was in congruence with phylogenetic analyses across replication-competent viral strains in databases. The tumor-associated viral population was less variable (π = 0.1% and π(AA) = 0.08%) than that observed in other clinical infections. Moreover, 42/46 (91.3%) viral genes analyzed had dN/dS scores of <1, which is indicative of high amino acid sequence conservation. Taken together, these findings raise the possibility that replication-competent HCMV may exist in malignant gliomas.

A selection index for gene expression evolution and its application to the divergence between humans and chimpanzees

The importance of gene regulation in animal evolution is a matter of long-standing interest, but measuring the impact of selection on gene expression has proven a challenge. Here, we propose a selection index of gene expression as a straightforward method for assessing the mode and strength of selection operating on gene expression levels. The index is based on the widely used McDonald-Kreitman test and requires the estimation of four quantities: the within-species and between-species expression variances as well as the sequence heterozygosity and divergence of neutrally evolving sequences. We apply the method to data from human and chimpanzee lymphoblastoid cell lines and show that gene expression is in general under strong stabilizing selection. We also demonstrate how the same framework can be used to estimate the proportion of adaptive gene expression evolution.

A fine-scale chimpanzee genetic map from population sequencing

To study the evolution of recombination rates in apes, we developed methodology to construct a fine-scale genetic map from high-throughput sequence data from 10 Western chimpanzees, Pan troglodytes verus. Compared to the human genetic map, broad-scale recombination rates tend to be conserved, but with exceptions, particularly in regions of chromosomal rearrangements and around the site of ancestral fusion in human chromosome 2. At fine scales, chimpanzee recombination is dominated by hotspots, which show no overlap with those of humans even though rates are similarly elevated around CpG islands and decreased within genes. The hotspot-specifying protein PRDM9 shows extensive variation among Western chimpanzees, and there is little evidence that any sequence motifs are enriched in hotspots. The contrasting locations of hotspots provide a natural experiment, which demonstrates the impact of recombination on base composition.

A trans-specific polymorphism in ZC3HAV1 is maintained by long-standing balancing selection and may confer susceptibility to multiple sclerosis

The human ZC3HAV1 gene encodes an antiviral protein. The longest splicing isoform of ZC3HAV1 contains a C-terminal PARP-like domain, which has evolved under positive selection in primates. We analyzed the evolutionary history of this same domain in humans and in Pan troglodytes. We identified two variants that segregate in both humans and chimpanzees; one of them (rs3735007) does not occur at a hypermutable site and accounts for a nonsynonymous substitution (Thr851Ile). The probability that the two trans-specific polymorphisms have occurred independently in the two lineages was estimated to be low (P = 0.0054), suggesting that at least one of them has arisen before speciation and has been maintained by selection. Population genetic analyses in humans indicated that the region surrounding the shared variants displays strong evidences of long-standing balancing selection. Selection signatures were also observed in a chimpanzee population sample. Inspection of 1000 Genomes data confirmed these findings but indicated that search for selection signatures using low-coverage whole-genome data may need masking of repetitive sequences. A case–control study of more than 1,000 individuals from mainland Italy indicated that the Thr851Ile SNP is significantly associated with susceptibility to multiple sclerosis (MS) (odds ratio [OR] = 1.47, 95% confidence intervals [CI]: 1.08–1.99, P = 0.011). This finding was confirmed in a larger sample of 4,416 Sardinians cases/controls (OR = 1.18, 95% CI: 1.037–1.344, P = 0.011), but not in a population from Belgium. We provide one of the first instances of human/chimpanzee trans-specific coding variant located outside the major histocompatibility complex region. The selective pressure is likely to be virus driven; in modern populations, this variant associates with susceptibility to MS, possibly via the interaction with environmental factors.

Comparative RNA sequencing reveals substantial genetic variation in endangered primates

Comparative genomic studies in primates have yielded important insights into the evolutionary forces that shape genetic diversity and revealed the likely genetic basis for certain species-specific adaptations. To date, however, these studies have focused on only a small number of species. For the majority of nonhuman primates, including some of the most critically endangered, genome-level data are not yet available. In this study, we have taken the first steps toward addressing this gap by sequencing RNA from the livers of multiple individuals from each of 16 mammalian species, including humans and 11 nonhuman primates. Of the nonhuman primate species, five are lemurs and two are lorisoids, for which little or no genomic data were previously available. To analyze these data, we developed a method for de novo assembly and alignment of orthologous gene sequences across species. We assembled an average of 5721 gene sequences per species and characterized diversity and divergence of both gene sequences and gene expression levels. We identified patterns of variation that are consistent with the action of positive or directional selection, including an 18-fold enrichment of peroxisomal genes among genes whose regulation likely evolved under directional selection in the ancestral primate lineage. Importantly, we found no relationship between genetic diversity and endangered status, with the two most endangered species in our study, the black and white ruffed lemur and the Coquerel's sifaka, having the highest genetic diversity among all primates. Our observations imply that many endangered lemur populations still harbor considerable genetic variation. Timely efforts to conserve these species alongside their habitats have, therefore, strong potential to achieve long-term success.

A genome sequence resource for the aye-aye (Daubentonia madagascariensis), a nocturnal lemur from Madagascar

We present a high-coverage draft genome assembly of the aye-aye (Daubentonia madagascariensis), a highly unusual nocturnal primate from Madagascar. Our assembly totals ~3.0 billion bp (3.0 Gb), roughly the size of the human genome, comprised of ~2.6 million scaffolds (N50 scaffold size = 13,597 bp) based on short paired-end sequencing reads. We compared the aye-aye genome sequence data with four other published primate genomes (human, chimpanzee, orangutan, and rhesus macaque) as well as with the mouse and dog genomes as nonprimate outgroups. Unexpectedly, we observed strong evidence for a relatively slow substitution rate in the aye-aye lineage compared with these and other primates. In fact, the aye-aye branch length is estimated to be ~10% shorter than that of the human lineage, which is known for its low substitution rate. This finding may be explained, in part, by the protracted aye-aye life-history pattern, including late weaning and age of first reproduction relative to other lemurs. Additionally, the availability of this draft lemur genome sequence allowed us to polarize nucleotide and protein sequence changes to the ancestral primate lineage-a critical period in primate evolution, for which the relevant fossil record is sparse. Finally, we identified 293,800 high-confidence single nucleotide polymorphisms in the donor individual for our aye-aye genome sequence, a captive-born individual from two wild-born parents. The resulting heterozygosity estimate of 0.051% is the lowest of any primate studied to date, which is understandable considering the aye-aye's extensive home-range size and relatively low population densities. Yet this level of genetic diversity also suggests that conservation efforts benefiting this unusual species should be prioritized, especially in the face of the accelerating degradation and fragmentation of Madagascar's forests.

Long-term balancing selection at the antiviral gene OAS1 in Central African chimpanzees

Oligoadenylate synthetases (OAS) are interferon-induced enzymes that participate in the first line of defense against a wide range of viral infection in animals. Upon activation by viral double-stranded RNA, OAS synthesizes (2-5) oligoadenylates, which activate RNase L, leading to the nonspecific degradation of cellular and viral RNA. Some association studies in humans suggest that variation at one of the OAS genes, OAS1, could be influencing host susceptibility to viral infection. We assessed the diversity of OAS1 in hominoid primates with a focus on chimpanzees. We found that the OAS1 gene is extremely polymorphic in Central African chimpanzee and exhibits levels of silent and replacement diversity much higher than neutral regions of the chimpanzee genome. This level of variation strongly suggests that balancing selection is acting on OAS1, and indeed, this conclusion was validated by several tests of neutrality. We further demonstrated that balancing selection has been acting at this locus since the split between chimpanzees, humans, and gorillas (~8.6 Ma) and caused the persistence of two deeply divergent allelic lineages in Central African chimpanzees. These two groups of OAS1 alleles differ by a large number of amino acids (a.a.), including several a.a. putatively involved in RNA binding. It is therefore very likely that variation at the OAS1 locus affects the innate immune response of individuals to specific viral infection. Our data strongly suggest that interactions between viral RNA and OAS1 are responsible for the maintenance of ancestral polymorphisms at this locus for at least 13.2 My.

Copy number variation analysis in the great apes reveals species-specific patterns of structural variation

Copy number variants (CNVs) are increasingly acknowledged as an important source of evolutionary novelties in the human lineage. However, our understanding of their significance is still hindered by the lack of primate CNV data. We performed intraspecific comparative genomic hybridizations to identify loci harboring copy number variants in each of the four great apes: bonobos, chimpanzees, gorillas, and orangutans. For the first time, we could analyze differences in CNV location and frequency in these four species, and compare them with human CNVs and primate segmental duplication (SD) maps. In addition, for bonobo and gorilla, patterns of CNV and nucleotide diversity were studied in the same individuals. We show that CNVs have been subject to different selective pressures in different lineages. Evidence for purifying selection is stronger in gorilla CNVs overlapping genes, while positive selection appears to have driven the fixation of structural variants in the orangutan lineage. In contrast, chimpanzees and bonobos present high levels of common structural polymorphism, which is indicative of relaxed purifying selection together with the higher mutation rates induced by the known burst of segmental duplication in the ancestor of the African apes. Indeed, the impact of the duplication burst is noticeable by the fact that bonobo and chimpanzee share more CNVs with gorilla than expected. Finally, we identified a number of interesting genomic regions that present high-frequency CNVs in all great apes, while containing only very rare or even pathogenic structural variants in humans.

Bonobos fall within the genomic variation of chimpanzees

To gain insight into the patterns of genetic variation and evolutionary relationships within and between bonobos and chimpanzees, we sequenced 150,000 base pairs of nuclear DNA divided among 15 autosomal regions as well as the complete mitochondrial genomes from 20 bonobos and 58 chimpanzees. Except for western chimpanzees, we found poor genetic separation of chimpanzees based on sample locality. In contrast, bonobos consistently cluster together but fall as a group within the variation of chimpanzees for many of the regions. Thus, while chimpanzees retain genomic variation that predates bonobo-chimpanzee speciation, extensive lineage sorting has occurred within bonobos such that much of their genome traces its ancestry back to a single common ancestor that postdates their origin as a group separate from chimpanzees.

Comparative and demographic analysis of orang-utan genomes

“Orangutan” is derived from the Malay term “man of the forest” and aptly describes the Southeast Asian great apes native to Sumatra and Borneo. The orangutan species, Pongo abelii (Sumatran) and Pongo pygmaeus (Bornean), are the most phylogenetically distant great apes from humans, thereby providing an informative perspective on hominid evolution. Here we present a Sumatran orangutan draft genome assembly and short read sequence data from five Sumatran and five Bornean orangutan genomes. Our analyses reveal that, compared to other primates, the orangutan genome has many unique features. Structural evolution of the orangutan genome has proceeded much more slowly than other great apes, evidenced by fewer rearrangements, less segmental duplication, a lower rate of gene family turnover and surprisingly quiescent Alu repeats, which have played a major role in restructuring other primate genomes. We also describe the first primate polymorphic neocentromere, found in both Pongo species, emphasizing the gradual evolution of orangutan genome structure. Orangutans have extremely low energy usage for a eutherian mammal1, far lower than their hominid relatives. Adding their genome to the repertoire of sequenced primates illuminates new signals of positive selection in several pathways including glycolipid metabolism. From the population perspective, both Pongo species are deeply diverse; however, Sumatran individuals possess greater diversity than their Bornean counterparts, and more species-specific variation. Our estimate of Bornean/Sumatran speciation time, 400k years ago (ya), is more recent than most previous studies and underscores the complexity of the orangutan speciation process. Despite a smaller modern census population size, the Sumatran effective population size (N_e) expanded exponentially relative to the ancestral N_e after the split, while Bornean N_e declined over the same period. Overall, the resources and analyses presented here offer new opportunities in evolutionary genomics, insights into hominid biology, and an extensive database of variation for conservation efforts.

More reliable estimates of divergence times in Pan using complete mtDNA sequences and accounting for population structure

Here, we report the sequencing and analysis of eight complete mitochondrial genomes of chimpanzees (Pan troglodytes) from each of the three established subspecies (P. t. troglodytes, P. t. schweinfurthii and P. t. verus) and the proposed fourth subspecies (P. t. ellioti). Our population genetic analyses are consistent with neutral patterns of evolution that have been shaped by demography. The high levels of mtDNA diversity in western chimpanzees are unlike those seen at nuclear loci, which may reflect a demographic history of greater female to male effective population sizes possibly owing to the characteristics of the founding population. By using relaxed-clock methods, we have inferred a timetree of chimpanzee species and subspecies. The absolute divergence times vary based on the methods and calibration used, but relative divergence times show extensive uniformity. Overall, mtDNA produces consistently older times than those known from nuclear markers, a discrepancy that is reduced significantly by explicitly accounting for chimpanzee population structures in time estimation. Assuming the human-chimpanzee split to be between 7 and 5 Ma, chimpanzee time estimates are 2.1-1.5, 1.1-0.76 and 0.25-0.18 Ma for the chimpanzee/bonobo, western/(eastern + central) and eastern/central chimpanzee divergences, respectively.

Genomic-scale capture and sequencing of endogenous DNA from feces

Genomic-level analyses of DNA from non-invasive sources would facilitate powerful conservation and evolutionary studies in natural populations of endangered and otherwise elusive species. However, the typical low quantity and poor quality of DNA that is extracted from non-invasive samples have generally precluded such work. Here we apply a modified DNA capture protocol that, when used in combination with massively-parallel sequencing technology, facilitates efficient and highly-accurate resequencing of megabases of specified nuclear genomic regions from fecal DNA samples. We validated our approach by comparing genetic variants identified from corresponding fecal and blood DNA samples of six western chimpanzees (Pan troglodytes verus) across more than 1.5 megabases of chromosome 21, chromosome X, and the complete mitochondrial genome. Our results suggest that it is now feasible to conduct genomic studies in natural populations for which constraints on invasive sampling have otherwise long been a barrier. The data we collected also provided an opportunity to examine western chimpanzee genetic diversity at unprecedented scale. Despite high mitochondrial genome diversity (π = 0.585%), western chimpanzees have a low ratio (0.42) of X chromosomal (π = 0.034%) to autosomal (chromosome 21 π = 0.081%) sequence diversity, a pattern that may reflect an unusual demographic history of this subspecies.

Diversification of bitter taste receptor gene family in western chimpanzees

In mammals, bitter taste is mediated by T2R genes, which belong to the large family of seven transmembrane G protein-coupled receptors. Because T2Rs are directly involved in the interaction between mammals and their dietary sources, it is likely that these genes evolved to reflect species' specific diets during mammalian evolution. Here, we investigated the sequences of all 28 putative functional chimpanzee T2R genes (cT2Rs) in 46 western chimpanzees to compare the intraspecies variations in chimpanzees to those already known for all 25 human functional T2R genes (hT2Rs). The numbers of functional genes varied among individuals in western chimpanzees, and most chimpanzees had two or three more functional genes than humans. Similarly to hT2Rs, cT2Rs showed high nucleotide diversity along with a large number of amino acid substitutions. Comparison of the nucleotide substitution patterns in cT2Rs with those in five cT2R pseudogenes and 14 autosomal intergenic noncoding regions among the same individuals revealed that the evolution of cT2R genes was almost identical to that of putative neutral regions with slight but significantly positive Tajima's D values, suggesting that selective constraint on these genes was relaxed with weak balancing selection. These trends have resulted in the occurrence of various divergent alleles of T2Rs within the western chimpanzee populations and in heterozygous individuals who might have the ability to taste a broader range of substances.

Y chromosomal variation tracks the evolution of mating systems in chimpanzee and bonobo

The male-specific regions of the Y chromosome (MSY) of the human and the chimpanzee (Pan troglodytes) are fully sequenced. The most striking difference is the dramatic rearrangement of large parts of their respective MSYs. These non-recombining regions include ampliconic gene families that are known to be important for male reproduction,and are consequently under significant selective pressure. However, whether the published Y-chromosomal pattern of ampliconic fertility genes is invariable within P. troglodytes is an open but fundamental question pertinent to discussions of the evolutionary fate of the Y chromosome in different primate mating systems. To solve this question we applied fluorescence in situ hybridisation (FISH) of testis-specific expressed ampliconic fertility genes to metaphase Y chromosomes of 17 chimpanzees derived from 11 wild-born males and 16 bonobos representing seven wild-born males. We show that of eleven P. troglodytes Y-chromosomal lines, ten Y-chromosomal variants were detected based on the number and arrangement of the ampliconic fertility genes DAZ (deleted in azoospermia) and CDY (chromodomain protein Y)-a so-far never-described variation of a species' Y chromosome. In marked contrast, no variation was evident among seven Y-chromosomal lines of the bonobo, P. paniscus, the chimpanzee's closest living relative. Although, loss of variation of the Y chromosome in the bonobo by a founder effect or genetic drift cannot be excluded, these contrasting patterns might be explained in the context of the species' markedly different social and mating behaviour. In chimpanzees, multiple males copulate with a receptive female during a short period of visible anogenital swelling, and this may place significant selection on fertility genes. In bonobos, however, female mate choice may make sperm competition redundant (leading to monomorphism of fertility genes), since ovulation in this species is concealed by the prolonged anogenital swelling, and because female bonobos can occupy high-ranking positions in the group and are thus able to determine mate choice more freely.

Evolutionary history of chimpanzees inferred from complete mitochondrial genomes

Investigations into the evolutionary history of the common chimpanzee, Pan troglodytes, have produced inconsistent results due to differences in the types of molecular data considered, the model assumptions employed, and the quantity and geographical range of samples used. We amplified and sequenced 24 complete P. troglodytes mitochondrial genomes from fecal samples collected at multiple study sites throughout sub-Saharan Africa. Using a "relaxed molecular clock," fossil calibrations, and 12 additional complete primate mitochondrial genomes, we analyzed the pattern and timing of primate diversification in a Bayesian framework. Our results support the recognition of four chimpanzee subspecies. Within P. troglodytes, we report a mean (95% highest posterior density [HPD]) time since most recent common ancestor (tMRCA) of 1.026 (0.811-1.263) Ma for the four proposed subspecies, with two major lineages. One of these lineages (tMRCA = 0.510 [0.387-0.650] Ma) contains P. t. verus (tMRCA = 0.155 [0.101-0.213] Ma) and P. t. ellioti (formerly P. t. vellerosus; tMRCA = 0.157 [0.102-0.215] Ma), both of which are monophyletic. The other major lineage contains P. t. schweinfurthii (tMRCA = 0.111 [0.077-0.146] Ma), a monophyletic clade nested within the P. t. troglodytes lineage (tMRCA = 0.380 [0.296-0.476] Ma). We utilized two analysis techniques that may be of widespread interest. First, we implemented a Yule speciation prior across the entire primate tree with separate coalescent priors on each of the chimpanzee subspecies. The validity of this approach was confirmed by estimates based on more traditional techniques. We also suggest that accurate tMRCA estimates from large computationally difficult sequence alignments may be obtained by implementing our novel method of bootstrapping smaller randomly subsampled alignments.

Evolutionary genetics in wild primates: combining genetic approaches with field studies of natural populations

Ecological and evolutionary studies of wild primates hold important keys to understanding both the shared characteristics of primate biology and the genetic and phenotypic differences that make specific lineages, including our own, unique. Although complementary genetic research on nonhuman primates has long been of interest, recent technological and methodological advances now enable functional and population genetic studies in an unprecedented manner. In the past several years, novel genetic data sets have revealed new information about the demographic history of primate populations and the genetics of adaptively important traits. In combination with the rich history of behavioral, ecological, and physiological work on natural primate populations, genetic approaches promise to provide a compelling picture of primate evolution in the past and in the present day.

Can intra-Y gene conversion oppose the degeneration of the human Y chromosome? A simulation study

The human Y is a genetically degenerate chromosome, which has lost about 97% of the genes originally present. Most of the remaining human Y genes are in large duplicated segments (ampliconic regions) undergoing intense Y–Y gene conversion. It has been suggested that Y–Y gene conversion may help these genes getting rid of deleterious mutations that would inactivate them otherwise. Here, we tested this idea by simulating the evolution of degenerating Y chromosomes with or without gene conversion using the most up-to-date population genetics parameters for humans. We followed the fate of a variant with Y–Y gene conversion in a population of Y chromosomes where Y–Y gene conversion is originally absent. We found that this variant gets fixed more frequently than the neutral expectation, which supports the idea that gene conversion is beneficial for a degenerating Y chromosome. Interestingly, a very high rate of gene conversion is needed for an effect of gene conversion to be observed. This suggests that high levels of Y-Y gene conversion observed in humans may have been selected to oppose the Y degeneration. We also studied with a similar approach the evolution of ampliconic regions on the Y chromosomes and found that the fixation of many copies at once is unlikely, which suggest these regions probably evolved gradually unless selection for increased dosage favored large-scale duplication events. Exploring the parameter space showed that Y–Y gene conversion may be beneficial in most mammalian species, which is consistent with recent data in chimpanzees and mice.

Dental and phylogeographic patterns of variation in gorillas

Gorilla patterns of variation have great relevance for studies of human evolution. In this study, molar morphometrics were used to evaluate patterns of geographic variation in gorillas. Dental specimens of 323 adult individuals, drawn from the current distribution of gorillas in equatorial Africa were divided into 14 populations. Discriminant analyses and Mahalanobis distances were used to study population structure. Results reveal that: 1) the West and East African gorillas form distinct clusters, 2) the Cross River gorillas are well separated from the rest of the western populations, 3) gorillas from the Virunga mountains and the Bwindi Forest can be differentiated from the lowland gorillas of Utu and Mwenga-Fizi, 4) the Tshiaberimu gorillas are distinct from other eastern gorillas, and the Kahuzi-Biega gorillas are affiliated with them. These findings provide support for a species distinction between Gorilla gorilla and Gorilla beringei, with subspecies G. g. diehli, G. g. gorilla, G. b. graueri, G. b. beringei, and possibly, G. b. rex-pygmaeorum. Clear correspondence between dental and other patterns of taxonomic diversity demonstrates that dental data reveal underlying genetic patterns of differentiation. Dental distances increased predictably with altitude but not with geographic distances, indicating that altitudinal segregation explains gorilla patterns of population divergence better than isolation-by-distance. The phylogeographic pattern of gorilla dental metric variation supports the idea that Plio-Pleistocene climatic fluctuations and local mountain building activity in Africa affected gorilla phylogeography. I propose that West Africa comprised the historic center of gorilla distribution and experienced drift-gene flow equilibrium, whereas Nigeria and East Africa were at the periphery, where climatic instability and altitudinal variation promoted drift and genetic differentiation. This understanding of gorilla population structure has implications for gorilla conservation, and for understanding the distribution of sympatric chimpanzees and Plio-Pleistocene hominins.

Adaptation and constraint at Toll-like receptors in primates

Frequent positive selection is a hallmark of genes involved in the adaptive immune system of vertebrates, but the incidence of positive selection for genes underlying innate immunity in vertebrates has not been well studied. The toll-like receptors (TLRs) of the innate immune system represent the first line of defense against pathogens. TLRs lie directly at the host-environment interface, and they target microbial molecules. Because of this, they might be subject to frequent positive selection due to coevolutionary dynamics with their microbial counterparts. However, they also recognize conserved molecular motifs, and this might constrain their evolution. Here, we investigate the evolution of the ten human TLRs in the framework of these competing ideas. We studied rates of protein evolution among primate species and we analyzed patterns of polymorphism in humans and chimpanzees. This provides a window into TLR evolution at both long and short timescales. We found a clear signature of positive selection in the rates of substitution across primates in most TLRs. Some of the implicated sites fall in structurally important protein domains, involve radical amino acid changes, or overlap with polymorphisms with known clinical associations in humans. However, within species, patterns of nucleotide variation were generally compatible with purifying selection, and these patterns differed between humans and chimpanzees and between viral and nonviral TLRs. Thus, adaptive evolution at TLRs does not appear to reflect a constant turnover of alleles and instead might be more episodic in nature. This pattern is consistent with more ephemeral pathogen-host associations rather than with long-term coevolution.

On the diversity of malaria parasites in African apes and the origin of Plasmodium falciparum from Bonobos

The origin of Plasmodium falciparum, the etiological agent of the most dangerous forms of human malaria, remains controversial. Although investigations of homologous parasites in African Apes are crucial to resolve this issue, studies have been restricted to a chimpanzee parasite related to P. falciparum, P. reichenowi, for which a single isolate was available until very recently. Using PCR amplification, we detected Plasmodium parasites in blood samples from 18 of 91 individuals of the genus Pan, including six chimpanzees (three Pan troglodytes troglodytes, three Pan t. schweinfurthii) and twelve bonobos (Pan paniscus). We obtained sequences of the parasites' mitochondrial genomes and/or from two nuclear genes from 14 samples. In addition to P. reichenowi, three other hitherto unknown lineages were found in the chimpanzees. One is related to P. vivax and two to P. falciparum that are likely to belong to distinct species. In the bonobos we found P. falciparum parasites whose mitochondrial genomes indicated that they were distinct from those present in humans, and another parasite lineage related to P. malariae. Phylogenetic analyses based on this diverse set of Plasmodium parasites in African Apes shed new light on the evolutionary history of P. falciparum. The data suggested that P. falciparum did not originate from P. reichenowi of chimpanzees (Pan troglodytes), but rather evolved in bonobos (Pan paniscus), from which it subsequently colonized humans by a host-switch. Finally, our data and that of others indicated that chimpanzees and bonobos maintain malaria parasites, to which humans are susceptible, a factor of some relevance to the renewed efforts to eradicate malaria.

Elucidating population histories using genomic DNA sequences

In 1993, Cliff Jolly suggested that rather than debating species definitions and classifications, energy would be better spent investigating multidimensional patterns of variation and gene flow among populations. Until now, however, genetic studies of wild primate populations have been limited to very small portions of the genome. Access to complete genome sequences of humans, chimpanzees, macaques, and other primates makes it possible to design studies surveying substantial amounts of DNA sequence variation at multiple genetic loci in representatives of closely related but distinct wild primate populations. Such data can be analyzed with new approaches that estimate not only when populations diverged but also the relative amounts and directions of subsequent gene flow. These analyses will reemphasize the difficulty of achieving consistent species and subspecies definitions by revealing the extent of variation in the amount and duration of gene flow accompanying population divergences.

The divergence of chimpanzee species and subspecies as revealed in multipopulation isolation-with-migration analyses

The divergence of bonobos and three subspecies of the common chimpanzee was examined under a multipopulation isolation-with-migration (IM) model with data from 73 loci drawn from the literature. A benefit of having a full multipopulation model, relative to conducting multiple pairwise analyses between sampled populations, is that a full model can reveal historical gene flow involving ancestral populations. An example of this was found in which gene flow is indicated between the western common chimpanzee subspecies and the ancestor of the central and the eastern common chimpanzee subspecies. The results of a full analysis on all four populations are strongly consistent with analyses on pairs of populations and generally similar to results from previous studies. The basal split between bonobos and common chimpanzees was estimated at 0.93 Ma (0.68-1.54 Ma, 95% highest posterior density interval), with the split among the ancestor of three common chimpanzee populations at 0.46 Ma (0.35-0.65), and the most recent split between central and eastern common chimpanzee populations at 0.093 Ma (0.041-0.157). Population size estimates mostly fell in the range from 5,000 to 10,000 individuals. The exceptions are the size of the ancestor of the common chimpanzee and the bonobo, at 17,000 (8,000-28,000) individuals, and the central common chimpanzee and its immediate ancestor with the eastern common chimpanzee, which have effective size estimates at 27,000 (16,000-44,000) and 32,000 (19,000-54,000) individuals, respectively.

Sequencing primate genomes: what have we learned?

We summarize the progress in whole-genome sequencing and analyses of primate genomes. These emerging genome datasets have broadened our understanding of primate genome evolution revealing unexpected and complex patterns of evolutionary change. This includes the characterization of genome structural variation, episodic changes in the repeat landscape, differences in gene expression, new models regarding speciation, and the ephemeral nature of the recombination landscape. The functional characterization of genomic differences important in primate speciation and adaptation remains a significant challenge. Limited access to biological materials, the lack of detailed phenotypic data and the endangered status of many critical primate species have significantly attenuated research into the genetic basis of primate evolution. Next-generation sequencing technologies promise to greatly expand the number of available primate genome sequences; however, such draft genome sequences will likely miss critical genetic differences within complex genomic regions unless dedicated efforts are put forward to understand the full spectrum of genetic variation.

The action of falciparum malaria on the human and chimpanzee genomes compared: absence of evidence for a genomic signature of malaria at HBB and G6PD in three subspecies of chimpanzee

The historical association between Plasmodium and primates has meant that many Plasmodium species have coevolved with specific primate hosts. However, unlike humans that are infected by species such as P. falciparum that cause severe malaria, many non-human primates are infected by Plasmodium species that only cause mild disease. Here we investigate whether the genomic signatures of plasmodial infection found in humans are also present in chimpanzees. We find no evidence of the major deleterious mutations at HBB (beta-globin) and G6PD in chimpanzees that confer resistance to malaria caused by P. falciparum nor evidence of long-term balancing selection at these loci. Our knowledge of malaria prevalence and pathogenesis in wild chimpanzees is severely limited, but it may be the case that beta-globin and G6PD variation are not adaptive in chimpanzees because malaria is rare and/or less detrimental in this species. Alternatively, chimpanzees may utilise mechanisms that are different from those of humans to protect against malaria.

A comparative approach shows differences in patterns of numt insertion during hominoid evolution

Nuclear integrations of mitochondrial DNA (numts) are widespread among eukaryotes, although their prevalence differs greatly among taxa. Most knowledge of numt evolution comes from analyses of whole-genome sequences of single species or, more recently, from genomic comparisons across vast phylogenetic distances. Here we employ a comparative approach using human and chimpanzee genome sequence data to infer differences in the patterns and processes underlying numt integrations. We identified 66 numts that have integrated into the chimpanzee nuclear genome since the human-chimp divergence, which is significantly greater than the 37 numts observed in humans. By comparing these closely related species, we accurately reconstructed the preintegration target site sequence and deduced nucleotide changes associated with numt integration. From >100 species-specific numts, we quantified the frequency of small insertions, deletions, duplications, and instances of microhomology. Most human and chimpanzee numt integrations were accompanied by microhomology and short indels of the kind typically observed in the nonhomologous end-joining pathway of DNA double-strand break repair. Human-specific numts have integrated into regions with a significant deficit of transposable elements; however, the same was not seen in chimpanzees. From a separate data set, we also found evidence for an apparent increase in the rate of numt insertions in the last common ancestor of humans and the great apes using a polymerase chain reaction-based screen. Last, phylogenetic analyses indicate that mitochondrial-numt alignments must be at least 500 bp, and preferably >1 kb in length, to accurately reconstruct hominoid phylogeny and recover the correct point of numt insertion.

Evolutionary toggling of the MAPT 17q21.31 inversion region

Using comparative sequencing approaches, we investigated the evolutionary history of the European-enriched 17q21.31 MAPT inversion polymorphism. We present a detailed, BAC-based sequence assembly of the inverted human H2 haplotype and compare it to the sequence structure and genetic variation of the corresponding 1.5-Mb region for the noninverted H1 human haplotype and that of chimpanzee and orangutan. We found that inversion of the MAPT region is similarly polymorphic in other great ape species, and we present evidence that the inversions occurred independently in chimpanzees and humans. In humans, the inversion breakpoints correspond to core duplications with the LRRC37 gene family. Our analysis favors the H2 configuration and sequence haplotype as the likely great ape and human ancestral state, with inversion recurrences during primate evolution. We show that the H2 architecture has evolved more extensive sequence homology, perhaps explaining its tendency to undergo microdeletion associated with mental retardation in European populations.

A history of recurrent positive selection at the toll-like receptor 5 in primates

Many genes involved in immunity evolve rapidly. It remains unclear, however, to what extent pattern-recognition receptors (PRRs) of the innate immune system in vertebrates are subject to recurrent positive selection imposed by pathogens, as suggested by studies in Drosophila, or whether they are evolutionarily constrained. Here, we show that Toll-like receptor 5 (TLR5), a member of the Toll-like receptor family of innate immunity genes that responds to bacterial flagellin, has undergone a history of adaptive evolution in primates. We have identified specific residues that have changed multiple times, sometimes in parallel in primates, and are thus likely candidates for selection. Most of these changes map to the extracellular leucine-rich repeats involved in pathogen recognition, and some are likely to have an effect on protein function due to the radical nature of the amino acid substitutions that are involved. These findings suggest that vertebrate PRRs might show similar patterns of evolution to Drosophila PRRs, in spite of the acquisition of the more complex and specific vertebrate adaptive immune system. At shorter timescales, however, we found no evidence of adaptive evolution in either humans or chimpanzees. In fact, we found that one mutation that abolishes TLR5 function is present at high frequencies in many human populations. Patterns of variation indicate that this mutation is not young, and its high frequency suggests some functional redundancy for this PRR in humans.

Hominoid seminal protein evolution and ancestral mating behavior

Hominoid mating systems show extensive variation among species. The degree of sexual dimorphism in body size and canine size varies among primates in accordance with their mating system, as does the testes size and the consistency of ejaculated semen, in response to differing levels of sperm competition. To investigate patterns of evolution at hominoid seminal proteins and to make inferences regarding the mating systems of extinct taxa, we sequenced the entire coding region of the prostate-specific transglutaminase (TGM4) gene in human, chimpanzee, bonobo, western lowland gorilla, eastern lowland gorilla, orangutan, and siamang, including multiple humans, chimps, and gorillas. Partial DNA sequence of the coding regions was also obtained for one eastern lowland gorilla at the semenogelin genes (SEMG1 and SEMG2), which code for the predominant proteins in semen. Patterns of nucleotide variation and inferred protein sequence change were evaluated within and between species. Combining the present data with previous studies demonstrates a high rate of amino acid substitutions, and low intraspecific variation, at seminal proteins in Pan, presumably driven by strong sperm competition. Both gorilla species apparently possess nonfunctional TGM4, SEMG1, and SEMG2 genes, suggesting that gorillas have had low sperm competition, and therefore their current polygynous mating system, for a long time before their divergence. Similarly, orangutans show longstanding stasis at TGM4, which may be interpreted as evidence for an unchanging mating system for most of their evolution after their divergence from African apes. In contrast to the great apes, the data from humans could be interpreted as evidence of fluctuations between different mating systems or alternatively as a relaxed functional constraint in these proteins. It is our hope that this study is a first step toward developing a model to predict ancestral mating systems from extant molecular data to complement interpretations from the fossil record.

Strong signature of natural selection within an FHIT intron implicated in prostate cancer risk

Previously, a candidate gene linkage approach on brother pairs affected with prostate cancer identified a locus of prostate cancer susceptibility at D3S1234 within the fragile histidine triad gene (FHIT), a tumor suppressor that induces apoptosis. Subsequent association tests on 16 SNPs spanning approximately 381 kb surrounding D3S1234 in Americans of European descent revealed significant evidence of association for a single SNP within intron 5 of FHIT. In the current study, re-sequencing and genotyping within a 28.5 kb region surrounding this SNP further delineated the association with prostate cancer risk to a 15 kb region. Multiple SNPs in sequences under evolutionary constraint within intron 5 of FHIT defined several related haplotypes with an increased risk of prostate cancer in European-Americans. Strong associations were detected for a risk haplotype defined by SNPs 138543, 142413, and 152494 in all cases (Pearson's chi(2) = 12.34, df 1, P = 0.00045) and for the homozygous risk haplotype defined by SNPs 144716, 142413, and 148444 in cases that shared 2 alleles identical by descent with their affected brothers (Pearson's chi(2) = 11.50, df 1, P = 0.00070). In addition to highly conserved sequences encompassing SNPs 148444 and 152413, population studies revealed strong signatures of natural selection for a 1 kb window covering the SNP 144716 in two human populations, the European American (pi = 0.0072, Tajima's D = 3.31, 14 SNPs) and the Japanese (pi = 0.0049, Fay & Wu's H = 8.05, 14 SNPs), as well as in chimpanzees (Fay & Wu's H = 8.62, 12 SNPs). These results strongly support the involvement of the FHIT intronic region in an increased risk of prostate cancer.