Demographic history and genetic differentiation in apes

Modeling of African population history using f-statistics is biased when applying all previously proposed SNP ascertainment schemes

f-statistics have emerged as a first line of analysis for making inferences about demographic history from genome-wide data. Not only are they guaranteed to allow robust tests of the fits of proposed models of population history to data when analyzing full genome sequencing data-that is, all single nucleotide polymorphisms (SNPs) in the individuals being analyzed-but they are also guaranteed to allow robust tests of models for SNPs ascertained as polymorphic in a population that is an outgroup in a phylogenetic sense to all groups being analyzed. True "outgroup ascertainment" is in practice impossible in humans because our species has arisen from a substructured ancestral population that does not descend from a homogeneous ancestral population going back many hundreds of thousands of years into the past. However, initial studies suggested that non-outgroup-ascertainment schemes might produce robust enough results using f-statistics, and that motivated widespread fitting of models to data using non-outgroup-ascertained SNP panels such as the "Affymetrix Human Origins array" which has been genotyped on thousands of modern individuals from hundreds of populations, or the "1240k" in-solution enrichment reagent which has been the source of about 70% of published genome-wide data for ancient humans. In this study, we show that while analyses of population history using such panels work well for studies of relationships among non-African populations and one African outgroup, when co-modeling more than one sub-Saharan African and/or archaic human groups (Neanderthals and Denisovans), fitting of f-statistics to such SNP sets is expected to frequently lead to false rejection of true demographic histories, and failure to reject incorrect models. Analyzing panels of SNPs polymorphic in archaic humans, which has been suggested as a solution for the ascertainment problem, has limited statistical power and retains important biases. However, by carrying out simulations of diverse demographic histories, we show that bias in inferences based on f-statistics can be minimized by ascertaining on variants common in a union of diverse African groups; such ascertainment retains high statistical power while allowing co-analysis of archaic and modern groups.

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.

Modeling of African population history using f -statistics can be highly biased and is not addressed by previously suggested SNP ascertainment schemes

f -statistics have emerged as a first line of analysis for making inferences about demographic history from genome-wide data. These statistics can provide strong evidence for either admixture or cladality, which can be robust to substantial rates of errors or missing data. f -statistics are guaranteed to be unbiased under "SNP ascertainment" (analyzing non-randomly chosen subsets of single nucleotide polymorphisms) only if it relies on a population that is an outgroup for all groups analyzed. However, ascertainment on a true outgroup that is not co-analyzed with other populations is often impractical and uncommon in the literature. In this study focused on practical rather than theoretical aspects of SNP ascertainment, we show that many non-outgroup ascertainment schemes lead to false rejection of true demographic histories, as well as to failure to reject incorrect models. But the bias introduced by common ascertainments such as the 1240K panel is mostly limited to situations when more than one sub-Saharan African and/or archaic human groups (Neanderthals and Denisovans) or non-human outgroups are co-modelled, for example, f 4 -statistics involving one non-African group, two African groups, and one archaic group. Analyzing panels of SNPs polymorphic in archaic humans, which has been suggested as a solution for the ascertainment problem, cannot fix all these problems since for some classes of f -statistics it is not a clean outgroup ascertainment, and in other cases it demonstrates relatively low power to reject incorrect demographic models since it provides a relatively small number of variants common in anatomically modern humans. And due to the paucity of high-coverage archaic genomes, archaic individuals used for ascertainment often act as sole representatives of the respective groups in an analysis, and we show that this approach is highly problematic. By carrying out large numbers of simulations of diverse demographic histories, we find that bias in inferences based on f -statistics introduced by non-outgroup ascertainment can be minimized if the derived allele frequency spectrum in the population used for ascertainment approaches the spectrum that existed at the root of all groups being co-analyzed. Ascertaining on sites with variants common in a diverse group of African individuals provides a good approximation to such a set of SNPs, addressing the great majority of biases and also retaining high statistical power for studying population history. Such a "pan-African" ascertainment, although not completely problem-free, allows unbiased exploration of demographic models for the widest set of archaic and modern human populations, as compared to the other ascertainment schemes we explored.

Population genomic monitoring provides insight into conservation status but no correlation with demographic estimates of extinction risk in a threatened trout

The current extinction crisis requires effective assessment and monitoring tools. Genetic approaches are appealing given the relative ease of field sampling required to estimate genetic diversity characteristics assumed related to population size, evolutionary potential, and extinction risk, and to evaluate hybridization with non‐native species simultaneously. However, linkages between population genetic metrics of diversity from survey‐style field collections and demographic estimates of population size and extinction risk are still in need of empirical examples, especially for remotely distributed species of conservation concern where the approach might be most beneficial. We capitalized on an exceptional opportunity to evaluate congruence between genetic diversity metrics and demographic‐based estimates of abundance and extinction risk from a comprehensive Multiple Population Viability Analysis (MPVA) in a threatened fish, the Lahontan cutthroat trout (LCT). We sequenced non‐native trout reference samples and recently collected and archived tissue samples of most remaining populations of LCT (N = 60) and estimated common genetic assessment metrics, predicting minimal hybridization with non‐native trout, low diversity, and declining diversity over time. We further hypothesized genetic metrics would correlate positively with MPVA‐estimated abundance and negatively with extinction probability. We uncovered several instances of hybridization that pointed to immediate management needs. After removing hybridized individuals, cautious interpretation of low effective population sizes (2–63) suggested reduced evolutionary potential for many LCT populations. Other genetic metrics did not decline over time nor correlate with MPVA‐based estimates of harmonic mean abundance or 30‐year extinction probability. Our results demonstrate benefits of genetic monitoring for efficiently detecting hybridization and, though genetic results were disconnected from demographic assessment of conservation status, they suggest reduced evolutionary potential and likely a higher conservation risk than currently recognized for this threatened fish. We emphasize that genetic information provides essential complementary insight, in addition to demographic information, for evaluating species status.

Population dynamics and genetic connectivity in recent chimpanzee history

Knowledge on the population history of endangered species is critical for conservation, but whole-genome data on chimpanzees (Pan troglodytes) is geographically sparse. Here, we produced the first non-invasive geolocalized catalog of genomic diversity by capturing chromosome 21 from 828 non-invasive samples collected at 48 sampling sites across Africa. The four recognized subspecies show clear genetic differentiation correlating with known barriers, while previously undescribed genetic exchange suggests that these have been permeable on a local scale. We obtained a detailed reconstruction of population stratification and fine-scale patterns of isolation, migration, and connectivity, including a comprehensive picture of admixture with bonobos (Pan paniscus). Unlike humans, chimpanzees did not experience extended episodes of long-distance migrations, which might have limited cultural transmission. Finally, based on local rare variation, we implement a fine-grained geolocalization approach demonstrating improved precision in determining the origin of confiscated chimpanzees.

Morphological variation of the maxilla in modern humans and African apes

Differences in morphology among modern humans and African apes are frequently used when assessing whether hominin fossils should be attributed to a single species or represent evidence for taxic diversity. A good understanding of the degree and structure of the intergeneric, interspecific, and intraspecific variation, including aspects such as sexual dimorphism and age, are key in this context. Here we explore the variation and differences shown by the maxilla of extant hominines, as maxillary morphology is central in the diagnosis of several hominin taxa. Our sample includes adults of all currently recognized hominine species and subspecies, with a balanced species sex ratio. In addition, we compared the adults with a small sample of late juveniles. The morphology of the maxillae was captured using three-dimensional landmarks, and the size and shape were analyzed using geometric morphometric methods. Key observations are that 1) the maxillae of all extant hominine species and subspecies show statistically significant differences, but complete separation in shape is only seen at the genus level; 2) the degree of variation is not consistent between genera, with subspecies of Gorilla being more different from each other than are species of Pan; 3) the pattern of sexual shape dimorphism is different in Pan, Gorilla, and Homo, often showing opposite trends; and 4) differentiation between maxillary shapes is increased after adjustment for static intraspecific allometry. These results provide a taxonomically up-to-date comparative morphological framework to help interpret the hominin fossil record, and we discuss the practical implications in that context.

How methodological changes have influenced our understanding of population structure in threatened species: insights from tiger populations across India

Unprecedented advances in sequencing technology in the past decade allow a better understanding of genetic variation and its partitioning in natural populations. Such inference is critical to conservation: to understand species biology and identify isolated populations. We review empirical population genetics studies of Endangered Bengal tigers within India, where 60–70% of wild tigers live. We assess how changes in marker type and sampling strategy have impacted inferences by reviewing past studies, and presenting three novel analyses including a single-nucleotide polymorphism (SNP) panel, genome-wide SNP markers, and a whole-mitochondrial genome network. At a broad spatial scale, less than 100 SNPs revealed the same patterns of population clustering as whole genomes (with the exception of one additional population sampled only in the SNP panel). Mitochondrial DNA indicates a strong structure between the northeast and other regions. Two studies with more populations sampled revealed further substructure within Central India. Overall, the comparison of studies with varied marker types and sample sets allows more rigorous inference of population structure. Yet sampling of some populations is limited across all studies, and these should be the focus of future sampling efforts. We discuss challenges in our understanding of population structure, and how to further address relevant questions in conservation genetics.

This article is part of the theme issue ‘Celebrating 50 years since Lewontin's apportionment of human diversity’.

Genetic Diversity in Chimpanzee Transcriptomics Does Not Represent Wild Populations

Chimpanzees (Pan troglodytes) are a genetically diverse species, consisting of four highly distinct subspecies. As humans' closest living relative, they have been a key model organism in the study of human evolution, and comparisons of human and chimpanzee transcriptomes have been widely used to characterize differences in gene expression levels that could underlie the phenotypic differences between the two species. However, the subspecies from which these transcriptomic data sets have been derived is not recorded in metadata available in the public NCBI Sequence Read Archive (SRA). Furthermore, labeling of RNA sequencing (RNA-seq) samples is for the most part inconsistent across studies, and the true number of individuals from whom transcriptomic data are available is difficult to ascertain. Thus, we have evaluated genetic diversity at the subspecies and individual level in 486 public RNA-seq samples available in the SRA, spanning the vast majority of public chimpanzee transcriptomic data. Using multiple population genetics approaches, we find that nearly all samples (96.6%) have some degree of Western chimpanzee ancestry. At the individual donor level, we identify multiple samples that have been repeatedly analyzed across different studies and identify a total of 135 genetically distinct individuals within our data, a number that falls to 89 when we exclude likely first- and second-degree relatives. Altogether, our results show that current transcriptomic data from chimpanzees are capturing low levels of genetic diversity relative to what exists in wild chimpanzee populations. These findings provide important context to current comparative transcriptomics research involving chimpanzees.

Cochlear shape distinguishes southern African early hominin taxa with unique auditory ecologies

Insights into potential differences among the bony labyrinths of Plio-Pleistocene hominins may inform their evolutionary histories and sensory ecologies. We use four recently-discovered bony labyrinths from the site of Kromdraai to significantly expand the sample for Paranthropus robustus. Diffeomorphometry, which provides detailed information about cochlear shape, reveals size-independent differences in cochlear shape between P. robustus and Australopithecus africanus that exceed those among modern humans and the African apes. The cochlea of P. robustus is distinctive and relatively invariant, whereas cochlear shape in A. africanus is more variable, resembles that of early Homo, and shows a degree of morphological polymorphism comparable to that evinced by modern species. The curvature of the P. robustus cochlea is uniquely derived and is consistent with enhanced sensitivity to low-frequency sounds. Combined with evidence for selection, our findings suggest that sound perception shaped distinct ecological adaptations among southern African early hominins.

The evolution of morality and the role of commitment

A considerable share of the literature on the evolution of human cooperation considers the question why we have not evolved to play the Nash equilibrium in prisoners' dilemmas or public goods games. In order to understand human morality and pro-social behaviour, we suggest that it would actually be more informative to investigate why we have not evolved to play the subgame perfect Nash equilibrium in sequential games, such as the ultimatum game and the trust game. The 'rationally irrational' behaviour that can evolve in such games gives a much better match with actual human behaviour, including elements of morality such as honesty, responsibility and sincerity, as well as the more hostile aspects of human nature, such as anger and vengefulness. The mechanism at work here is commitment, which does not need population structure, nor does it need interactions to be repeated. We argue that this shift in focus can not only help explain why humans have evolved to know wrong from right, but also why other animals, with similar population structures and similar rates of repetition, have not evolved similar moral sentiments. The suggestion that the evolutionary function of morality is to help us commit to otherwise irrational behaviour stems from the work of Robert Frank (American Economic Review, 77(4), 593-604, 1987; Passions within reason: The strategic role of the emotions, WW Norton, 1988), which has played a surprisingly modest role in the scientific debate to date.

Recent genetic connectivity and clinal variation in chimpanzees

Much like humans, chimpanzees occupy diverse habitats and exhibit extensive behavioural variability. However, chimpanzees are recognized as a discontinuous species, with four subspecies separated by historical geographic barriers. Nevertheless, their range-wide degree of genetic connectivity remains poorly resolved, mainly due to sampling limitations. By analyzing a geographically comprehensive sample set amplified at microsatellite markers that inform recent population history, we found that isolation by distance explains most of the range-wide genetic structure of chimpanzees. Furthermore, we did not identify spatial discontinuities corresponding with the recognized subspecies, suggesting that some of the subspecies-delineating geographic barriers were recently permeable to gene flow. Substantial range-wide genetic connectivity is consistent with the hypothesis that behavioural flexibility is a salient driver of chimpanzee responses to changing environmental conditions. Finally, our observation of strong local differentiation associated with recent anthropogenic pressures portends future loss of critical genetic diversity if habitat fragmentation and population isolation continue unabated.

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.

Population history of the golden eagle inferred from whole-genome sequencing of three of its subspecies

The application of evolutionary genetic research to investigate the potential for endangered species to adapt to changing environments is important for conservation biology. Effective population size (Ne) is informative for understanding adaptive potential as it refers to the genetic variation in breeding individuals who have contributed to contemporary and historic population diversity. We reconstruct fluctuations in Ne in three golden eagle subspecies (Japanese, Scottish, North American) using the pairwise sequential Markovian coalescent (PSMC) model based on whole-genome sequence data. Our results indicate the timing of subspeciation events and suggest significant ongoing demographic reductions since the start of the Last Glacial Period. Importantly, we find evidence for gene flow from continental populations into the ancestral Japanese population resulting in a short, sharp recovery in genetic diversity. Timing agrees with the palaeogeographic estimates of land bridge connections between the Japanese archipelago and Asian continent and matches a similar Ne spike in the Scottish population, but not in the North American population. Given contemporary declines in isolated Japanese and UK island populations, our study highlights a concerning loss of local genetic diversity, but also indicates the likely response of populations to genetic reinforcement from neighbouring subspecies, increasing management options and encouraging a range-wide species conservation approach.

Differences in MHC-B diversity and KIR epitopes in two populations of wild chimpanzees

The major histocompatibility complex (MHC) class I genes play a critical role within the immune system, both by the presentation of antigens from intracellular pathogens to immunocompetent cells and by the interaction with killer cell immunoglobulin-like receptors (KIR) on natural killer cells (NK cells). Genes of the MHC are highly diverse, and MHC variation can have effects on the immune functionality of individuals; hence, comparisons of MHC diversity among closely related phylogenetic taxa may give insight into the factors responsible for the shaping of its diversity. The four geographically separated chimpanzee subspecies differ in their overall genetic diversity, have different population histories, and are confronted with different pathogens in their natural habitat, all of which may affect MHC class I DNA sequence diversity. Here, we compare the MHC-B exon two DNA sequence diversity from 24 wild western and 46 wild eastern chimpanzees using necropsy and noninvasively collected fecal samples, respectively. We found a higher MHC-B exon two nucleotide diversity, in our western than eastern chimpanzees. The inclusion of previously published MHC-B exon two data from other western and eastern chimpanzees supported this finding. In addition, our results confirm and extend the finding of a very low C1 epitope frequency at eastern chimpanzee MHC-B molecules, which likely affects the ability of these molecules to interact with NK cells. While the understanding of the differing pathogen environments encountered by disparate populations of a species is a challenging endeavor, these findings highlight the potential for these pathogens to selectively shape immune system variation.

Humans and Chimpanzees Display Opposite Patterns of Diversity in Arylamine N-Acetyltransferase Genes

Among the many genes involved in the metabolism of therapeutic drugs, human arylamine N-acetyltransferases (NATs) genes have been extensively studied, due to their medical importance both in pharmacogenetics and disease epidemiology. One member of this small gene family, NAT2, is established as the locus of the classic human acetylation polymorphism in drug metabolism. Current hypotheses hold that selective processes favoring haplotypes conferring lower NAT2 activity have been operating in modern humans’ recent history as an adaptation to local chemical and dietary environments. To shed new light on such hypotheses, we investigated the genetic diversity of the three members of the NAT gene family in seven hominid species, including modern humans, Neanderthals and Denisovans. Little polymorphism sharing was found among hominids, yet all species displayed high NAT diversity, but distributed in an opposite fashion in chimpanzees and bonobos (Pan genus) compared to modern humans, with higher diversity in Pan species at NAT1 and lower at NAT2, while the reverse is observed in humans. This pattern was also reflected in the results returned by selective neutrality tests, which suggest, in agreement with the predicted functional impact of mutations detected in non-human primates, stronger directional selection, presumably purifying selection, at NAT1 in modern humans, and at NAT2 in chimpanzees. Overall, the results point to the evolution of divergent functions of these highly homologous genes in the different primate species, possibly related to their specific chemical/dietary environment (exposome) and we hypothesize that this is likely linked to the emergence of controlled fire use in the human lineage.

The effect of sexual selection on phenotypic diversification among human populations: A simulation study

Despite the generally low level of inter-population genetic differentiation in humans as compared with great apes, it has long been acknowledged that there is a considerable amount of geographic variations in human phenotypes, for example, skin pigmentation, cranial morphology, and soft-tissue facial morphology, to name but a few. Indeed, recent studies have suggested that the extent of inter-population diversity in some human phenotypes is greater than expected from random drift alone. Such an excess of phenotypic diversity is often attributed to adaptation to local environment. However, this account is valid only if populations differ in some ecological aspects that elicit differential selection acting on a given phenotypic feature. Another long-standing hypothesis is the sexual selection hypothesis, which claims that phenotypic diversity arises and/or is maintained owing to variations in preference for mating partners. In this paper, we explore the plausibility of the sexual selection hypothesis by means of computer simulations, in which the inter-population diversity of a quantitative trait is evaluated against the expectation from random drift, using the Q_ST-F_ST comparison. As possible driving factors of sexual selection, we consider two types of mate-choice preference: preference for the population average and preference for a culturally-transmitted arbitrary trend. Our simulations suggest that sexual selection can, under certain circumstances, maintain and/or generate a detectable amount of inter-population phenotypic diversity, even when populations are ecologically identical and loosely connected to each other by mutual migration. Since mating decisions in humans are considerably affected by social learning, human mate-choice preference may be more readily diversified between populations than in other animals. We suggest, therefore, that some of the observed human phenotypic variations may be better understood as a product of cultural, rather than ecological, diversification.

Reduced bonobo MHC class I diversity predicts a reduced viral peptide binding ability compared to chimpanzees

BACKGROUND

The highly polymorphic genes of the major histocompatibility complex (MHC) class I are involved in defense against viruses and other intracellular pathogens. Although several studies found reduced MHC class I diversity in bonobos in comparison to the closely related chimpanzee, it is unclear if this lower diversity also influences the functional ability of MHC class I molecules in bonobos. Here, we use a bioinformatic approach to analyze the viral peptide binding ability of all published bonobo MHC class I molecules (n = 58) in comparison to all published chimpanzee MHC class I molecules (n = 161) for the class I loci A, B, C and A-like.

RESULTS

We examined the peptide binding ability of all 219 different MHC class I molecules to 5,788,712 peptides derived from 1432 different primate viruses and analyzed the percentage of bound peptides and the overlap of the peptide binding repertoires of the two species. We conducted multiple levels of analysis on the "species"-, "population"- and "individual"-level to account for the characterization of MHC variation in a larger number of chimpanzees and their broader geographic distribution. We found a lower percentage of bound peptides in bonobos at the B locus in the "population"-level comparison and at the B and C loci in the "individual"-level comparison. Furthermore, we found evidence of a limited peptide binding repertoire in bonobos by tree-based visualization of functional clustering of MHC molecules, as well as an analysis of peptides bound by both species.

CONCLUSION

Our results suggest a reduced MHC class I viral peptide binding ability at the B and C loci in bonobos compared to chimpanzees. The effects of this finding on the immune defense against viruses in wild living bonobos are unclear. However, special caution is needed to prevent introduction and spread of new viruses to bonobos, as their defensive ability to cope with new viruses could be limited compared to chimpanzees.

Hominoid intraspecific cranial variation mirrors neutral genetic diversity

In humans, patterns of cranial variation mirror genetic diversity globally, implicating population history as a key driver of cranial disparity. Here, we demonstrate that the magnitude of genetic diversity within 12 extant ape taxa explains a large proportion of cranial shape variation. Taxa that are more genetically diverse tend to be more cranially diverse also. Our results suggest that neutral evolutionary processes such as mutation, genetic drift, and gene flow are reflected in both genetic and cranial diversity in apes. This work provides a perspective on intraspecific cranial variation in apes which has important implications for interpreting selective and developmental pressures on the cranium and for understanding shape variation in fossil hominin crania.

Natural selection, developmental constraint, and plasticity have all been invoked as explanations for intraspecific cranial variation in humans and apes. However, global patterns of human cranial variation are congruent with patterns of genetic variation, demonstrating that population history has influenced cranial variation in humans. Here we show that this finding is not unique to Homo sapiens but is also broadly evident across extant ape species. Specifically, taxa that exhibit greater intraspecific cranial shape variation also exhibit greater genetic diversity at neutral autosomal loci. Thus, cranial shape variation within hominoid taxa reflects the population history of each species. Our results suggest that neutral evolutionary processes such as mutation, gene flow, and genetic drift have played an important role in generating cranial variation within species. These findings are consistent with previous work on human cranial morphology and improve our understanding of the evolutionary processes that generate intraspecific cranial shape diversity within hominoids. This work has implications for the analysis of selective and developmental pressures on the cranium and for interpreting shape variation in fossil hominin crania.

Direct estimate of the spontaneous germ line mutation rate in African green monkeys

Here, I provide the first direct estimate of the spontaneous mutation rate in an Old World monkey, using a seven individual, three-generation pedigree of African green monkeys. Eight de novo mutations were identified within ∼1.5 Gbp of accessible genome, corresponding to an estimated point mutation rate of 0.94 × 10^-8 per site per generation, suggesting an effective population size of ∼12000 for the species. This estimation represents a significant improvement in our knowledge of the population genetics of the African green monkey, one of the most important nonhuman primate models in biomedical research. Furthermore, by comparing mutation rates in Old World monkeys with the only other direct estimates in primates to date-humans and chimpanzees-it is possible to uniquely address how mutation rates have evolved over longer time scales. While the estimated spontaneous mutation rate for African green monkeys is slightly lower than the rate of 1.2 × 10^-8 per base pair per generation reported in chimpanzees, it is similar to the lower range of rates of 0.96 × 10^-8 -1.28 × 10^-8 per base pair per generation recently estimated from whole genome pedigrees in humans. This result suggests a long-term constraint on mutation rate that is quite different from similar evidence pertaining to recombination rate evolution in primates.

The Demographic and Adaptive History of the African Green Monkey

Relatively little is known about the evolutionary history of the African green monkey (genus Chlorocebus) due to the lack of sampled polymorphism data from wild populations. Yet, this characterization of genetic diversity is not only critical for a better understanding of their own history, but also for human biomedical research given that they are one of the most widely used primate models. Here, I analyze the demographic and selective history of the African green monkey, utilizing one of the most comprehensive catalogs of wild genetic diversity to date, consisting of 1,795,643 autosomal single nucleotide polymorphisms in 25 individuals, representing all five major populations: C. a. aethiops, C. a. cynosurus, C. a. pygerythrus, C. a. sabaeus, and C. a tantalus. Assuming a mutation rate of 5.9 × 10-9 per base pair per generation and a generation time of 8.5 years, divergence time estimates range from 523 to 621 kya for the basal split of C. a. aethiops from the other four populations. Importantly, the resulting tree characterizing the relationship and split-times between these populations differs significantly from that presented in the original genome paper, owing to their neglect of within-population variation when calculating between population-divergence. In addition, I find that the demographic history of all five populations is well explained by a model of population fragmentation and isolation, rather than novel colonization events. Finally, utilizing these demographic models as a null, I investigate the selective history of the populations, identifying candidate regions potentially related to adaptation in response to pathogen exposure.

Teleosts as Model Organisms To Understand Host-Microbe Interactions

Host-microbe interactions are influenced by complex host genetics and environment. Studies across animal taxa have aided our understanding of how intestinal microbiota influence vertebrate development, disease, and physiology. However, traditional mammalian studies can be limited by the use of isogenic strains, husbandry constraints that result in small sample sizes and limited statistical power, reliance on indirect characterization of gut microbial communities from fecal samples, and concerns of whether observations in artificial conditions are actually reflective of what occurs in the wild. Fish models are able to overcome many of these limitations. The extensive variation in the physiology, ecology, and natural history of fish enriches studies of the evolution and ecology of host-microbe interactions. They share physiological and immunological features common among vertebrates, including humans, and harbor complex gut microbiota, which allows identification of the mechanisms driving microbial community assembly. Their accelerated life cycles and large clutch sizes and the ease of sampling both internal and external microbial communities make them particularly well suited for robust statistical studies of microbial diversity. Gnotobiotic techniques, genetic manipulation of the microbiota and host, and transparent juveniles enable novel insights into mechanisms underlying development of the digestive tract and disease states. Many diseases involve a complex combination of genes which are difficult to manipulate in homogeneous model organisms. By taking advantage of the natural genetic variation found in wild fish populations, as well as of the availability of powerful genetic tools, future studies should be able to identify conserved genes and pathways that contribute to human genetic diseases characterized by dysbiosis.

Significant differentiation in the apolipoprotein(a)/lipoprotein(a) trait between chimpanzees from Western and Central Africa

Elevated Lipoprotein(a) (Lp(a)) plasma concentrations are a risk factor for cardiovascular disease in humans, largely controlled by the LPA gene encoding apolipoprotein(a) (apo(a)). Lp(a) is composed of low-density lipoprotein (LDL) and apo(a) and restricted to Catarrhini. A variable number of kringle IV (KIV) domains in LPA lead to a size polymorphism of apo(a) that is inversely correlated with Lp(a) concentrations. Smaller apo(a) isoforms and higher Lp(a) levels in central chimpanzees (Pan troglodytes troglodytes [PTT]) compared to humans from Europe had been reported. We studied apo(a) isoforms and Lp(a) concentrations in 75 western (Pan troglodytes verus [PTV]) and 112 central chimpanzees, and 12 bonobos (Pan paniscus [PPA]), all wild born and living in sanctuaries in Sierra Leone, Republic of the Congo, and DR Congo, respectively, and 116 humans from Gabon. Lp(a) levels were severalfold higher in western than in central chimpanzees (181.0 ± 6.7 mg/dl vs. 56.5 ± 4.3 mg/dl), whereas bonobos showed intermediate levels (134.8 ± 33.4 mg/dl). Apo(a) isoform sizes differed significantly between subspecies (means 20.9 ± 2.2, 22.9 ± 4.4, and 23.8 ± 3.8 KIV repeats in PTV, PTT, and PPA, respectively). However, far higher isoform-associated Lp(a) concentrations for all isoform sizes in western chimpanzees offered the main explanation for the higher overall Lp(a) levels in this subspecies. Human Lp(a) concentrations (mean 47.9 ± 2.8 mg/dl) were similar to those in central chimpanzees despite larger isoforms (mean 27.1 ± 4.9 KIV). Lp(a) and LDL, apoB-100, and total cholesterol levels only correlated in PTV. This remarkable differentiation between chimpanzees from different African habitats and the trait's similarity in humans and chimpanzees from Central Africa poses the question of a possible impact of an environmental factor that has shaped the genetic architecture of LPA. Overall, studies on the cholesterol-containing particles of Lp(a) and LDL in chimpanzees should consider differentiation between subspecies.

MHC class I diversity in chimpanzees and bonobos

Major histocompatibility complex (MHC) class I genes are critically involved in the defense against intracellular pathogens. MHC diversity comparisons among samples of closely related taxa may reveal traces of past or ongoing selective processes. The bonobo and chimpanzee are the closest living evolutionary relatives of humans and last shared a common ancestor some 1 mya. However, little is known concerning MHC class I diversity in bonobos or in central chimpanzees, the most numerous and genetically diverse chimpanzee subspecies. Here, we used a long-read sequencing technology (PacBio) to sequence the classical MHC class I genes A, B, C, and A-like in 20 and 30 wild-born bonobos and chimpanzees, respectively, with a main focus on central chimpanzees to assess and compare diversity in those two species. We describe in total 21 and 42 novel coding region sequences for the two species, respectively. In addition, we found evidence for a reduced MHC class I diversity in bonobos as compared to central chimpanzees as well as to western chimpanzees and humans. The reduced bonobo MHC class I diversity may be the result of a selective process in their evolutionary past since their split from chimpanzees.

Unconstrained cranial evolution in Neandertals and modern humans compared to common chimpanzees

A variety of lines of evidence support the idea that neutral evolutionary processes (genetic drift, mutation) have been important in generating cranial differences between Neandertals and modern humans. But how do Neandertals and modern humans compare with other species? And how do these comparisons illuminate the evolutionary processes underlying cranial diversification? To address these questions, we used 27 standard cranial measurements collected on 2524 recent modern humans, 20 Neandertals and 237 common chimpanzees to estimate split times between Neandertals and modern humans, and between Pan troglodytes verus and two other subspecies of common chimpanzee. Consistent with a neutral divergence, the Neandertal versus modern human split-time estimates based on cranial measurements are similar to those based on DNA sequences. By contrast, the common chimpanzee cranial estimates are much lower than DNA-sequence estimates. Apparently, cranial evolution has been unconstrained in Neandertals and modern humans compared with common chimpanzees. Based on these and additional analyses, it appears that cranial differentiation in common chimpanzees has been restricted by stabilizing natural selection. Alternatively, this restriction could be due to genetic and/or developmental constraints on the amount of within-group variance (relative to effective population size) available for genetic drift to act on.

Great ape Y Chromosome and mitochondrial DNA phylogenies reflect subspecies structure and patterns of mating and dispersal

The distribution of genetic diversity in great ape species is likely to have been affected by patterns of dispersal and mating. This has previously been investigated by sequencing autosomal and mitochondrial DNA (mtDNA), but large-scale sequence analysis of the male-specific region of the Y Chromosome (MSY) has not yet been undertaken. Here, we use the human MSY reference sequence as a basis for sequence capture and read mapping in 19 great ape males, combining the data with sequences extracted from the published whole genomes of 24 additional males to yield a total sample of 19 chimpanzees, four bonobos, 14 gorillas, and six orangutans, in which interpretable MSY sequence ranges from 2.61 to 3.80 Mb. This analysis reveals thousands of novel MSY variants and defines unbiased phylogenies. We compare these with mtDNA-based trees in the same individuals, estimating time-to-most-recent common ancestor (TMRCA) for key nodes in both cases. The two loci show high topological concordance and are consistent with accepted (sub)species definitions, but time depths differ enormously between loci and (sub)species, likely reflecting different dispersal and mating patterns. Gorillas and chimpanzees/bonobos present generally low and high MSY diversity, respectively, reflecting polygyny versus multimale–multifemale mating. However, particularly marked differences exist among chimpanzee subspecies: The western chimpanzee MSY phylogeny has a TMRCA of only 13.2 (10.8–15.8) thousand years, but that for central chimpanzees exceeds 1 million years. Cross-species comparison within a single MSY phylogeny emphasizes the low human diversity, and reveals species-specific branch length variation that may reflect differences in long-term generation times.

Delimiting Species-Poor Data Sets using Single Molecular Markers: A Study of Barcode Gaps, Haplowebs and GMYC

Most single-locus molecular approaches to species delimitation available to date have been designed and tested on data sets comprising at least tens of species, whereas the opposite case (species-poor data sets for which the hypothesis that all individuals are conspecific cannot by rejected beforehand) has rarely been the focus of such attempts. Here we compare the performance of barcode gap detection, haplowebs and generalized mixed Yule-coalescent (GMYC) models to delineate chimpanzees and bonobos using nuclear sequence markers, then apply these single-locus species delimitation methods to data sets of one, three, or six species simulated under a wide range of population sizes, speciation rates, mutation rates and sampling efforts. Our results show that barcode gap detection and GMYC models are unable to delineate species properly in data sets composed of one or two species, two situations in which haplowebs outperform them. For data sets composed of three or six species, bGMYC and haplowebs outperform the single-threshold and multiple-threshold versions of GMYC, whereas a clear barcode gap is only observed when population sizes and speciation rates are both small. The latter conditions represent a "sweet spot" for molecular taxonomy where all the single-locus approaches tested work well; however, the performance of these methods decreases strongly when population sizes and speciation rates are high, suggesting that multilocus approaches may be necessary to tackle such cases.

Reconstructing the demographic history of orang-utans using Approximate Bayesian Computation

Investigating how different evolutionary forces have shaped patterns of DNA variation within and among species requires detailed knowledge of their demographic history. Orang-utans, whose distribution is currently restricted to the South-East Asian islands of Borneo (Pongo pygmaeus) and Sumatra (Pongo abelii), have likely experienced a complex demographic history, influenced by recurrent changes in climate and sea levels, volcanic activities and anthropogenic pressures. Using the most extensive sample set of wild orang-utans to date, we employed an Approximate Bayesian Computation (ABC) approach to test the fit of 12 different demographic scenarios to the observed patterns of variation in autosomal, X-chromosomal, mitochondrial and Y-chromosomal markers. In the best-fitting model, Sumatran orang-utans exhibit a deep split of populations north and south of Lake Toba, probably caused by multiple eruptions of the Toba volcano. In addition, we found signals for a strong decline in all Sumatran populations ~24 ka, probably associated with hunting by human colonizers. In contrast, Bornean orang-utans experienced a severe bottleneck ~135 ka, followed by a population expansion and substructuring starting ~82 ka, which we link to an expansion from a glacial refugium. We showed that orang-utans went through drastic changes in population size and connectedness, caused by recurrent contraction and expansion of rainforest habitat during Pleistocene glaciations and probably hunting by early humans. Our findings emphasize the fact that important aspects of the evolutionary past of species with complex demographic histories might remain obscured when applying overly simplified models.

Androgen receptor and monoamine oxidase polymorphism in wild bonobos

Androgen receptor gene (AR), monoamine oxidase A gene (MAOA) and monoamine oxidase B gene (MAOB) have been found to have associations with behavioral traits, such as aggressiveness, and disorders in humans. However, the extent to which similar genetic effects might influence the behavior of wild apes is unclear. We examined the loci AR glutamine repeat (ARQ), AR glycine repeat (ARG), MAOA intron 2 dinucleotide repeat (MAin2) and MAOB intron 2 dinucleotide repeat (MBin2) in 32 wild bonobos, Pan paniscus, and compared them with those of chimpanzees, Pan troglodytes, and humans. We found that bonobos were polymorphic on the four loci examined. Both loci MAin2 and MBin2 in bonobos showed a higher diversity than in chimpanzees. Because monoamine oxidase influences aggressiveness, the differences between the polymorphisms of MAin2 and MBin2 in bonobos and chimpanzees may be associated with the differences in aggression between the two species. In order to understand the evolution of these loci and AR, MAOA and MAOB in humans and non-human primates, it would be useful to conduct future studies focusing on the potential association between aggressiveness, and other personality traits, and polymorphisms documented in bonobos.

The sampling scheme matters: Pan troglodytes troglodytes and P. t. schweinfurthii are characterized by clinal genetic variation rather than a strong subspecies break

Populations of an organism living in marked geographical or evolutionary isolation from other populations of the same species are often termed subspecies and expected to show some degree of genetic distinctiveness. The common chimpanzee (Pan troglodytes) is currently described as four geographically delimited subspecies: the western (P. t. verus), the nigerian-cameroonian (P. t. ellioti), the central (P. t. troglodytes) and the eastern (P. t. schweinfurthii) chimpanzees. Although these taxa would be expected to be reciprocally monophyletic, studies have not always consistently resolved the central and eastern chimpanzee taxa. Most studies, however, used data from individuals of unknown or approximate geographic provenance. Thus, genetic data from samples of known origin may shed light on the evolutionary relationship of these subspecies. We generated microsatellite genotypes from noninvasively collected fecal samples of 185 central chimpanzees that were sampled across large parts of their range and analyzed them together with 283 published eastern chimpanzee genotypes from known localities. We observed a clear signal of isolation by distance across both subspecies. Further, we found that a large proportion of comparisons between groups taken from the same subspecies showed higher genetic differentiation than the least differentiated between-subspecies comparison. This proportion decreased substantially when we simulated a more clumped sampling scheme by including fewer groups. Our results support the general concept that the distribution of the sampled individuals can dramatically affect the inference of genetic population structure. With regard to chimpanzees, our results emphasize the close relationship of equatorial chimpanzees from central and eastern equatorial Africa and the difficult nature of subspecies definitions.

Great ape genomics

The great ape families are the species most closely related to our own, comprising chimpanzees, bonobos, gorillas, and orangutans. They live exclusively in tropical rainforests in Central Africa and the islands of Southeast Asia. Due to their close evolutionary relationship with humans, great apes share many cognitive, physiological, and morphological similarities with humans. The members of the great ape family make obvious models to facilitate the further understanding about humans' biology and history. This review will discuss how the recent addition of genome-wide data from great apes has furthered humans' understanding of these species and humanity, especially in the realm of evolutionary genetics.

Phenotypic variation in infants, not adults, reflects genotypic variation among chimpanzees and bonobos

Studies comparing phenotypic variation with neutral genetic variation in modern humans have shown that genetic drift is a main factor of evolutionary diversification among populations. The genetic population history of our closest living relatives, the chimpanzees and bonobos, is now equally well documented, but phenotypic variation among these taxa remains relatively unexplored, and phenotype-genotype correlations are not yet documented. Also, while the adult phenotype is typically used as a reference, it remains to be investigated how phenotype-genotye correlations change during development. Here we address these questions by analyzing phenotypic evolutionary and developmental diversification in the species and subspecies of the genus Pan. Our analyses focus on the morphology of the femoral diaphysis, which represents a functionally constrained element of the locomotor system. Results show that during infancy phenotypic distances between taxa are largely congruent with non-coding (neutral) genotypic distances. Later during ontogeny, however, phenotypic distances deviate from genotypic distances, mainly as an effect of heterochronic shifts between taxon-specific developmental programs. Early phenotypic differences between Pan taxa are thus likely brought about by genetic drift while late differences reflect taxon-specific adaptations.

Sequence diversity of Pan troglodytes subspecies and the impact of WFDC6 selective constraints in reproductive immunity

Recent efforts have attempted to describe the population structure of common chimpanzee, focusing on four subspecies: Pan troglodytes verus, P. t. ellioti, P. t. troglodytes, and P. t. schweinfurthii. However, few studies have pursued the effects of natural selection in shaping their response to pathogens and reproduction. Whey acidic protein (WAP) four-disulfide core domain (WFDC) genes and neighboring semenogelin (SEMG) genes encode proteins with combined roles in immunity and fertility. They display a strikingly high rate of amino acid replacement (d_N/d_S), indicative of adaptive pressures during primate evolution. In human populations, three signals of selection at the WFDC locus were described, possibly influencing the proteolytic profile and antimicrobial activities of the male reproductive tract. To evaluate the patterns of genomic variation and selection at the WFDC locus in chimpanzees, we sequenced 17 WFDC genes and 47 autosomal pseudogenes in 68 chimpanzees (15 P. t. troglodytes, 22 P. t. verus, and 31 P. t. ellioti). We found a clear differentiation of P. t. verus and estimated the divergence of P. t. troglodytes and P. t. ellioti subspecies in 0.173 Myr; further, at the WFDC locus we identified a signature of strong selective constraints common to the three subspecies in WFDC6—a recent paralog of the epididymal protease inhibitor EPPIN. Overall, chimpanzees and humans do not display similar footprints of selection across the WFDC locus, possibly due to different selective pressures between the two species related to immune response and reproductive biology.

Robustness to divergence time underestimation when inferring species trees from estimated gene trees

To infer species trees from gene trees estimated from phylogenomic data sets, tractable methods are needed that can handle dozens to hundreds of loci. We examine several computationally efficient approaches-MP-EST, STAR, STEAC, STELLS, and STEM-for inferring species trees from gene trees estimated using maximum likelihood (ML) and Bayesian approaches. Among the methods examined, we found that topology-based methods often performed better using ML gene trees and methods employing coalescent times typically performed better using Bayesian gene trees, with MP-EST, STAR, STEAC, and STELLS outperforming STEM under most conditions. We examine why the STEM tree (also called GLASS or Maximum Tree) is less accurate on estimated gene trees by comparing estimated and true coalescence times, performing species tree inference using simulations, and analyzing a great ape data set keeping track of false positive and false negative rates for inferred clades. We find that although true coalescence times are more ancient than speciation times under the multispecies coalescent model, estimated coalescence times are often more recent than speciation times. This underestimation can lead to increased bias and lack of resolution with increased sampling (either alleles or loci) when gene trees are estimated with ML. The problem appears to be less severe using Bayesian gene-tree estimates.

Different selective pressures shape the evolution of Toll-like receptors in human and African great ape populations

The study of the genetic and selective landscape of immunity genes across primates can provide insight into the existing differences in susceptibility to infection observed between human and non-human primates. Here, we explored how selection has driven the evolution of a key family of innate immunity receptors, the Toll-like receptors (TLRs), in African great ape species. We sequenced the 10 TLRs in various populations of chimpanzees and gorillas, and analysed these data jointly with a human data set. We found that purifying selection has been more pervasive in great apes than in humans. Furthermore, in chimpanzees and gorillas, purifying selection has targeted TLRs irrespectively of whether they are endosomal or cell surface, in contrast to humans where strong selective constraints are restricted to endosomal TLRs. These observations suggest important differences in the relative importance of TLR-mediated pathogen sensing, such as that of recognition of flagellated bacteria by TLR5, between humans and great apes. Lastly, we used a population genetics-phylogenetics method that jointly analyses polymorphism and divergence data to detect fine-scale variation in selection pressures at specific codons within TLR genes. We identified different codons at different TLRs as being under positive selection in each species, highlighting that functional variation at these genes has conferred a selective advantage in immunity to infection to specific primate species. Overall, this study showed that the degree of selection driving the evolution of TLRs has largely differed between human and non-human primates, increasing our knowledge on their respective biological contribution to host defence in the natural setting.

A genome-wide survey of genetic variation in gorillas using reduced representation sequencing

All non-human great apes are endangered in the wild, and it is therefore important to gain an understanding of their demography and genetic diversity. Whole genome assembly projects have provided an invaluable foundation for understanding genetics in all four genera, but to date genetic studies of multiple individuals within great ape species have largely been confined to mitochondrial DNA and a small number of other loci. Here, we present a genome-wide survey of genetic variation in gorillas using a reduced representation sequencing approach, focusing on the two lowland subspecies. We identify 3,006,670 polymorphic sites in 14 individuals: 12 western lowland gorillas (Gorilla gorilla gorilla) and 2 eastern lowland gorillas (Gorilla beringei graueri). We find that the two species are genetically distinct, based on levels of heterozygosity and patterns of allele sharing. Focusing on the western lowland population, we observe evidence for population substructure, and a deficit of rare genetic variants suggesting a recent episode of population contraction. In western lowland gorillas, there is an elevation of variation towards telomeres and centromeres on the chromosomal scale. On a finer scale, we find substantial variation in genetic diversity, including a marked reduction close to the major histocompatibility locus, perhaps indicative of recent strong selection there. These findings suggest that despite their maintaining an overall level of genetic diversity equal to or greater than that of humans, population decline, perhaps associated with disease, has been a significant factor in recent and long-term pressures on wild gorilla populations.

Inferring the evolutionary histories of divergences in Hylobates and Nomascus gibbons through multilocus sequence data

BACKGROUND

Gibbons (Hylobatidae) are the most diverse group of living apes. They exist as geographically-contiguous species which diverged more rapidly than did their close relatives, the great apes (Hominidae). Of the four extant gibbon genera, the evolutionary histories of two polyspecific genera, Hylobates and Nomascus, have been the particular focus of research but the DNA sequence data used was largely derived from the maternally inherited mitochondrial DNA (mtDNA) locus.

RESULTS

To investigate the evolutionary relationships and divergence processes of gibbon species, particularly those of the Hylobates genus, we produced and analyzed a total of 11.5 kb DNA of sequence at 14 biparentally inherited autosomal loci. We find that on average gibbon genera have a high average sequence diversity but a lower degree of genetic differentiation as compared to great ape genera. Our multilocus species tree features H. pileatus in a basal position and a grouping of the four Sundaic island species (H. agilis, H. klossii, H. moloch and H. muelleri). We conducted pairwise comparisons based on an isolation-with-migration (IM) model and detect signals of asymmetric gene flow between H. lar and H. moloch, between H. agilis and H. muelleri, and between N. leucogenys and N. siki.

CONCLUSIONS

Our multilocus analyses provide inferences of gibbon evolutionary histories complementary to those based on single gene data. The results of IM analyses suggest that the divergence processes of gibbons may be accompanied by gene flow. Future studies using analyses of multi-population model with samples of known provenance for Hylobates and Nomascus species would expand the understanding of histories of gene flow during divergences for these two gibbon genera.

Origins of HIV and the AIDS pandemic

Acquired immunodeficiency syndrome (AIDS) of humans is caused by two lentiviruses, human immunodeficiency viruses types 1 and 2 (HIV-1 and HIV-2). Here, we describe the origins and evolution of these viruses, and the circumstances that led to the AIDS pandemic. Both HIVs are the result of multiple cross-species transmissions of simian immunodeficiency viruses (SIVs) naturally infecting African primates. Most of these transfers resulted in viruses that spread in humans to only a limited extent. However, one transmission event, involving SIVcpz from chimpanzees in southeastern Cameroon, gave rise to HIV-1 group M-the principal cause of the AIDS pandemic. We discuss how host restriction factors have shaped the emergence of new SIV zoonoses by imposing adaptive hurdles to cross-species transmission and/or secondary spread. We also show that AIDS has likely afflicted chimpanzees long before the emergence of HIV. Tracing the genetic changes that occurred as SIVs crossed from monkeys to apes and from apes to humans provides a new framework to examine the requirements of successful host switches and to gauge future zoonotic risk.

Aye-aye population genomic analyses highlight an important center of endemism in northern Madagascar

We performed a population genomics study of the aye-aye, a highly specialized nocturnal lemur from Madagascar. Aye-ayes have low population densities and extensive range requirements that could make this flagship species particularly susceptible to extinction. Therefore, knowledge of genetic diversity and differentiation among aye-aye populations is critical for conservation planning. Such information may also advance our general understanding of Malagasy biogeography, as aye-ayes have the largest species distribution of any lemur. We generated and analyzed whole-genome sequence data for 12 aye-ayes from three regions of Madagascar (North, West, and East). We found that the North population is genetically distinct, with strong differentiation from other aye-ayes over relatively short geographic distances. For comparison, the average FST value between the North and East aye-aye populations--separated by only 248 km--is over 2.1-times greater than that observed between human Africans and Europeans. This finding is consistent with prior watershed- and climate-based hypotheses of a center of endemism in northern Madagascar. Taken together, these results suggest a strong and long-term biogeographical barrier to gene flow. Thus, the specific attention that should be directed toward preserving large, contiguous aye-aye habitats in northern Madagascar may also benefit the conservation of other distinct taxonomic units. To help facilitate future ecological- and conservation-motivated population genomic analyses by noncomputational biologists, the analytical toolkit used in this study is available on the Galaxy Web site.

Comparative population genomics of the ejaculate in humans and the great apes

The rapid molecular evolution of reproductive genes is nearly ubiquitous across animals, yet the selective forces and functional targets underlying this divergence remain poorly understood. Humans and closely related species of great apes show strongly divergent mating systems, providing a powerful system to investigate the influence of sperm competition on the evolution of reproductive genes. This is complemented by detailed information on male reproductive biology and unparalleled genomic resources in humans. Here, we have used custom microarrays to capture and sequence 285 genes encoding proteins present in the ejaculate as well as 101 randomly selected control genes in 21 gorillas, 20 chimpanzees, 20 bonobos, and 20 humans. In total, we have generated >25× average genomic coverage per individual for over 1 million target base pairs. Our analyses indicate high levels of evolutionary constraint across much of the ejaculate combined with more rapid evolution of genes involved in immune defense and proteolysis. We do not find evidence for appreciably more positive selection along the lineage leading to bonobos and chimpanzees, although this would be predicted given more intense sperm competition in these species. Rather, the extent of positive and negative selection depended more on the effective population sizes of the species. Thus, general patterns of male reproductive protein evolution among apes and humans depend strongly on gene function but not on inferred differences in the intensity of sperm competition among extant species.

Incomplete lineage sorting is common in extant gibbon genera

We sequenced reduced representation libraries by means of Illumina technology to generate over 1.5 Mb of orthologous sequence from a representative of each of the four extant gibbon genera (Nomascus, Hylobates, Symphalangus, and Hoolock). We used these data to assess the evolutionary relationships between the genera by evaluating the likelihoods of all possible bifurcating trees involving the four taxa. Our analyses provide weak support for a tree with Nomascus and Hylobates as sister taxa and with Hoolock and Symphalangus as sister taxa, though bootstrap resampling suggests that other phylogenetic scenarios are also possible. This uncertainty is due to short internal branch lengths and extensive incomplete lineage sorting across taxa. The true phylogenetic relationships among gibbon genera will likely require a more extensive whole-genome sequence analysis.

Eco-geographical diversification of bitter taste receptor genes (TAS2Rs) among subspecies of chimpanzees (Pan troglodytes)

Chimpanzees (Pan troglodytes) have region-specific difference in dietary repertoires from East to West across tropical Africa. Such differences may result from different genetic backgrounds in addition to cultural variations. We analyzed the sequences of all bitter taste receptor genes (cTAS2Rs) in a total of 59 chimpanzees, including 4 putative subspecies. We identified genetic variations including single-nucleotide variations (SNVs), insertions and deletions (indels), gene-conversion variations, and copy-number variations (CNVs) in cTAS2Rs. Approximately two-thirds of all cTAS2R haplotypes in the amino acid sequence were unique to each subspecies. We analyzed the evolutionary backgrounds of natural selection behind such diversification. Our previous study concluded that diversification of cTAS2Rs in western chimpanzees (P. t. verus) may have resulted from balancing selection. In contrast, the present study found that purifying selection dominates as the evolutionary form of diversification of the so-called human cluster of cTAS2Rs in eastern chimpanzees (P. t. schweinfurthii) and that the other cTAS2Rs were under no obvious selection as a whole. Such marked diversification of cTAS2Rs with different evolutionary backgrounds among subspecies of chimpanzees probably reflects their subspecies-specific dietary repertoires.

High diversity at PRDM9 in chimpanzees and bonobos

Background

The PRDM9 locus in mammals has increasingly attracted research attention due to its role in mediating chromosomal recombination and possible involvement in hybrid sterility and hence speciation processes. The aim of this study was to characterize sequence variation at the PRDM9 locus in a sample of our closest living relatives, the chimpanzees and bonobos.

Methodology/Principal Findings

PRDM9 contains a highly variable and repetitive zinc finger array. We amplified this domain using long-range PCR and determined the DNA sequences using conventional Sanger sequencing. From 17 chimpanzees representing three subspecies and five bonobos we obtained a total of 12 alleles differing at the nucleotide level. Based on a data set consisting of our data and recently published Pan PRDM9 sequences, we found that at the subspecies level, diversity levels did not differ among chimpanzee subspecies or between chimpanzee subspecies and bonobos. In contrast, the sample of chimpanzees harbors significantly more diversity at PRDM9 than samples of humans. Pan PRDM9 shows signs of rapid evolution including no alleles or ZnFs in common with humans as well as signals of positive selection in the residues responsible for DNA binding.

Conclusions and Significance

The high number of alleles specific to the genus Pan, signs of positive selection in the DNA binding residues, and reported lack of conservation of recombination hotspots between chimpanzees and humans suggest that PRDM9 could be active in hotspot recruitment in the genus Pan. Chimpanzees and bonobos are considered separate species and do not have overlapping ranges in the wild, making the presence of shared alleles at the amino acid level between the chimpanzee and bonobo species interesting in view of the hypothesis that PRDM9 plays a universal role in interspecific hybrid sterility.

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.

Genomic tools for evolution and conservation in the chimpanzee: Pan troglodytes ellioti is a genetically distinct population

In spite of its evolutionary significance and conservation importance, the population structure of the common chimpanzee, Pan troglodytes, is still poorly understood. An issue of particular controversy is whether the proposed fourth subspecies of chimpanzee, Pan troglodytes ellioti, from parts of Nigeria and Cameroon, is genetically distinct. Although modern high-throughput SNP genotyping has had a major impact on our understanding of human population structure and demographic history, its application to ecological, demographic, or conservation questions in non-human species has been extremely limited. Here we apply these tools to chimpanzee population structure, using ∼700 autosomal SNPs derived from chimpanzee genomic data and a further ∼100 SNPs from targeted re-sequencing. We demonstrate conclusively the existence of P. t. ellioti as a genetically distinct subgroup. We show that there is clear differentiation between the verus, troglodytes, and ellioti populations at the SNP and haplotype level, on a scale that is greater than that separating continental human populations. Further, we show that only a small set of SNPs (10–20) is needed to successfully assign individuals to these populations. Tellingly, use of only mitochondrial DNA variation to classify individuals is erroneous in 4 of 54 cases, reinforcing the dangers of basing demographic inference on a single locus and implying that the demographic history of the species is more complicated than that suggested analyses based solely on mtDNA. In this study we demonstrate the feasibility of developing economical and robust tests of individual chimpanzee origin as well as in-depth studies of population structure. These findings have important implications for conservation strategies and our understanding of the evolution of chimpanzees. They also act as a proof-of-principle for the use of cheap high-throughput genomic methods for ecological questions.

Chimpanzees are viewed with fondness as our closest animal relatives and are valued by scientists for the biological and evolutionary insights they provide. In spite of this, the relationships between different populations of common chimpanzees are still relatively poorly understood, a situation that potentially threatens conservation efforts. Here we have used information gathered in the Chimpanzee Genome Project to design comprehensive tests of genetic variability that show unambiguously the existence of four genetically distinct groups (or populations) of common chimpanzee. We demonstrate that previous methods based on mitochondrial DNA sequences alone are not always accurate and show the feasibility of cheap new genetic tests of individuals' origins that could play an important role in conservation.

Y-Chromosome variation in hominids: intraspecific variation is limited to the polygamous chimpanzee

Background

We have previously demonstrated that the Y-specific ampliconic fertility genes DAZ (deleted in azoospermia) and CDY (chromodomain protein Y) varied with respect to copy number and position among chimpanzees (Pan troglodytes). In comparison, seven Y-chromosomal lineages of the bonobo (Pan paniscus), the chimpanzee's closest living relative, showed no variation. We extend our earlier comparative investigation to include an analysis of the intraspecific variation of these genes in gorillas (Gorilla gorilla) and orangutans (Pongo pygmaeus), and examine the resulting patterns in the light of the species' markedly different social and mating behaviors.

Methodology/Principal Findings

Fluorescence in situ hybridization analysis (FISH) of DAZ and CDY in 12 Y-chromosomal lineages of western lowland gorilla (G. gorilla gorilla) and a single lineage of the eastern lowland gorilla (G. beringei graueri) showed no variation among lineages. Similar findings were noted for the 10 Y-chromosomal lineages examined in the Bornean orangutan (Pongo pygmaeus), and 11 Y-chromosomal lineages of the Sumatran orangutan (P. abelii). We validated the contrasting DAZ and CDY patterns using quantitative real-time polymerase chain reaction (qPCR) in chimpanzee and bonobo.

Conclusion/Significance

High intraspecific variation in copy number and position of the DAZ and CDY genes is seen only in the chimpanzee. We hypothesize that this is best explained by sperm competition that results in the variant DAZ and CDY haplotypes detected in this species. In contrast, bonobos, gorillas and orangutans—species that are not subject to sperm competition—showed no intraspecific variation in DAZ and CDY suggesting that monoandry in gorillas, and preferential female mate choice in bonobos and orangutans, probably permitted the fixation of a single Y variant in each taxon. These data support the notion that the evolutionary history of a primate Y chromosome is not simply encrypted in its DNA sequences, but is also shaped by the social and behavioral circumstances under which the specific species has evolved.

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.