Nucleotide diversity in gorillas

Comparison of the levels of nucleotide diversity in humans and apes may provide valuable information for inferring the demographic history of these species, the effect of social structure on genetic diversity, patterns of past migration, and signatures of past selection events. Previous DNA sequence data from both the mitochondrial and the nuclear genomes suggested a much higher level of nucleotide diversity in the African apes than in humans. Noting that the nuclear DNA data from the apes were very limited, we previously conducted a DNA polymorphism study in humans and another in chimpanzees and bonobos, using 50 DNA segments randomly chosen from the noncoding, nonrepetitive parts of the human genome. The data revealed that the nucleotide diversity (pi) in bonobos (0.077%) is actually lower than that in humans (0.087%) and that pi in chimpanzees (0.134%) is only 50% higher than that in humans. In the present study we sequenced the same 50 segments in 15 western lowland gorillas and estimated pi to be 0.158%. This is the highest value among the African apes but is only about two times higher than that in humans. Interestingly, available mtDNA sequence data also suggest a twofold higher nucleotide diversity in gorillas than in humans, but suggest a threefold higher nucleotide diversity in chimpanzees than in humans. The higher mtDNA diversity in chimpanzees might be due to the unique pattern in the evolution of chimpanzee mtDNA. From the nuclear DNA pi values, we estimated that the long-term effective population sizes of humans, bonobos, chimpanzees, and gorillas are, respectively, 10,400, 12,300, 21,300, and 25,200.     

Estimation of the ancestral effective population sizes of African great apes under different selection regimes

Reliable estimates of ancestral effective population sizes are necessary to unveil the population-level phenomena that shaped the phylogeny and molecular evolution of the African great apes. Although several methods have previously been applied to infer ancestral effective population sizes, an analysis of the influence of the selective regime on the estimates of ancestral demography has not been thoroughly conducted. In this study, three independent data sets under different selective regimes were used were composed to tackle this issue. The results showed that selection had a significant impact on the estimates of ancestral effective population sizes of the African great apes. The inference of the ancestral demography of African great apes was affected by the selection regime. The effects, however, were not homogeneous along the ancestral populations of great apes. The effective population size of the ancestor of humans and chimpanzees was more impacted by the selection regime when compared to the same parameter in the ancestor of humans, chimpanzees and gorillas. Because the selection regime influenced the estimates of ancestral effective population size, it is reasonable to assume that a portion of the discrepancy found in previous studies that inferred the ancestral effective population size may be attributable to the differential action of selection on the genes sampled.     

Assessing endocranial variations in great apes and humans using 3D data from virtual endocasts

Modern humans are characterized by their large, complex, and specialized brain. Human brain evolution can be addressed through direct evidence provided by fossil hominid endocasts (i.e. paleoneurology), or through indirect evidence of extant species comparative neurology. Here we use the second approach, providing an extant comparative framework for hominid paleoneurological studies. We explore endocranial size and shape differences among great apes and humans, as well as between sexes. We virtually extracted 72 endocasts, sampling all extant great ape species and modern humans, and digitized 37 landmarks on each for 3D generalized Procrustes analysis. All species can be differentiated by their endocranial shape. Among great apes, endocranial shapes vary from short (orangutans) to long (gorillas), perhaps in relation to different facial orientations. Endocranial shape differences among African apes are partly allometric. Major endocranial traits distinguishing humans from great apes are endocranial globularity, reflecting neurological reorganization, and features linked to structural responses to posture and bipedal locomotion. Human endocasts are also characterized by posterior location of foramina rotunda relative to optic canals, which could be correlated to lesser subnasal prognathism compared to living great apes. Species with larger brains (gorillas and humans) display greater sexual dimorphism in endocranial size, while sexual dimorphism in endocranial shape is restricted to gorillas, differences between males and females being at least partly due to allometry. Our study of endocranial variations in extant great apes and humans provides a new comparative dataset for studies of fossil hominid endocasts.     

Nuclear integrations of mitochondrial DNA in gorillas

Great ape systematics, particularly at the species level and below, is currently under debate, due in part to the recent influx of molecular data. The phylogenies of previously published mitochondrial control region (or D-loop) DNA sequences in gorillas show deep splits within West African gorillas (Gorilla gorilla gorilla), and very high levels of nucleotide diversity in this subspecies. Here we demonstrate that several previously reported D-loop haplotypes from West African gorillas are in all likelihood nuclear integrations of mitochondrial DNA. Revised estimates of the amount and pattern of mitochondrial DNA diversity in gorillas are provided, revealing two reciprocally monophyletic and highly divergent groups of gorillas, concurrent with their geographic distribution.     

Ancestral Population Sizes and Species Divergence Times in the Primate Lineage on the Basis of Intron and BAC End Sequences

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

Characterization of human-specific DNA sequences obtained by genome subtraction

Recent studies on molecular evolution using nucleotide sequence data to clarify phylogenetic relationships among humans and the African great apes, have revealed that humans are more closely related to chimpanzees than to gorillas. However, the genetic basis of human uniqueness remains unclear. This is because phylogenetic studies have merely evaluated the degree of similarity by calculating the accumulation of nucleotide substitutions that have occurred in neutral DNA regions commonly present in all the species examined. In contrast, the genome subtraction method recently developed by us has revealed dissimilarity even among the genomes of the most closely related species. Here we describe the characteristics of the DNA sequences obtained by genome subtraction between humans and chimpanzees.     

Estimating ancestral population sizes and divergence times

This article presents a new method for jointly estimating species divergence times and ancestral population sizes. The method improves on previous ones by explicitly incorporating intragenic recombination, by utilizing orthologous sequence data from closely related species, and by using a maximum-likelihood framework. The latter allows for efficient use of the available information and provides a way of assessing how much confidence we should place in the estimates. I apply the method to recently collected intergenic sequence data from humans and the great apes. The results suggest that the human-chimpanzee ancestral population size was four to seven times larger than the current human effective population size and that the current human effective population size is slightly >10,000. These estimates are similar to previous ones, and they appear relatively insensitive to assumptions about the recombination rates or mutation rates across loci.     

Humans and Great Apes Cohabiting the Forest Ecosystem in Central African Republic Harbour the Same Hookworms

Background

Hookworms are important pathogens of humans. To date, Necator americanus is the sole, known species of the genus Necator infecting humans. In contrast, several Necator species have been described in African great apes and other primates. It has not yet been determined whether primate-originating Necator species are also parasitic in humans.

Methodology/Principal Findings

The infective larvae of Necator spp. were developed using modified Harada-Mori filter-paper cultures from faeces of humans and great apes inhabiting Dzanga-Sangha Protected Areas, Central African Republic. The first and second internal transcribed spacers (ITS-1 and ITS-2) of nuclear ribosomal DNA and partial cytochrome c oxidase subunit 1 (cox1) gene of mtDNA obtained from the hookworm larvae were sequenced and compared. Three sequence types (I–III) were recognized in the ITS region, and 34 cox1 haplotypes represented three phylogenetic groups (A–C). The combinations determined were I-A, II-B, II-C, III-B and III-C. Combination I-A, corresponding to N. americanus, was demonstrated in humans and western lowland gorillas; II-B and II-C were observed in humans, western lowland gorillas and chimpanzees; III-B and III-C were found only in humans. Pairwise nucleotide difference in the cox1 haplotypes between the groups was more than 8%, while the difference within each group was less than 2.1%.

Conclusions/Significance

The distinctness of ITS sequence variants and high number of pairwise nucleotide differences among cox1 variants indicate the possible presence of several species of Necator in both humans and great apes. We conclude that Necator hookworms are shared by humans and great apes co-habiting the same tropical forest ecosystems.     

Using mitochondrial and nuclear DNA markers to reconstruct human evolution

Molecular genetic data have greatly improved our ability to test hypotheses about human evolution. During the past decade, a large amount of nuclear and mitochondrial data have been collected from diverse human populations. Taken together, these data indicate that modern humans are a relatively young species. African populations show the largest amount of genetic diversity, and they are the most genetically divergent population. Modern human populations expanded in size first on the African continent. These findings support a recent African origin of modern humans, but this conclusion should be tempered by the possible effects of factors such as gene flow, population size differences, and natural selection. BioEssays 20:126–136, 1998. © 1998 John Wiley & Sons, Inc.     

A HERV-K provirus in chimpanzees, bonobos and gorillas, but not humans

Evidence from DNA sequencing studies strongly indicated that humans and chimpanzees are more closely related to each other than either is to gorillas [1-4]. However, precise details of the nature of the evolutionary separation of the lineage leading to humans from those leading to the African great apes have remained uncertain. The unique insertion sites of endogenous retroviruses, like those of other transposable genetic elements, should be useful for resolving phylogenetic relationships among closely related species. We identified a human endogenous retrovirus K (HERV-K) provirus that is present at the orthologous position in the gorilla and chimpanzee genomes, but not in the human genome. Humans contain an intact preintegration site at this locus. These observations provide very strong evidence that, for some fraction of the genome, chimpanzees, bonobos, and gorillas are more closely related to each other than they are to humans. They also show that HERV-K replicated as a virus and reinfected the germline of the common ancestor of the four modern species during the period of time when the lineages were separating and demonstrate the utility of using HERV-K to trace human evolution.     

Comparison of DNA and protein polymorphisms between humans and chimpanzees

To examine the nucleotide diversity at silent (synonymous + intron + untranslated) and non-silent (nonsynonymous) sites in chimpanzees and humans, genes at six nuclear loci from two chimpanzees were sequenced. The average silent diversity was 0.19%, which was significantly higher than that in humans (0.05%). This observation suggests a significantly larger effective population size and a higher extent of neutral polymorphism in chimpanzees than in humans. On the other hand, the non-silent nucleotide diversity is similar in both species, resulting in a larger fraction of neutral mutations at non-silent sites in humans than in chimpanzees. Other types of polymorphism data were collected from the literature or databases to examine whether or not they are consistent with the nuclear DNA sequence polymorphism observed here. The nucleotide diversity at both silent and non-silent sites in mitochondrial (mt) DNA genes was compatible with that of the nuclear genes. Microsatellite loci showed a similar high extent of heterozygosity in both species, perhaps due to the combined effect of a high mutation rate and a recent population expansion in humans. At protein loci, humans are more heterozygous than chimpanzees, and the estimated fraction of neutral alleles in humans (0.84) is much larger than that in chimpanzees (0.26). These data show that the neutral fraction in non-silent changes is relatively large in the human population. This difference may be due to a relaxation of the functional constraint against proteins in the human lineage. To evaluate this possibility, it will be necessary to examine nucleotide sequences in relation to the physiological or biochemical properties of proteins.     

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.     

Evolution of a HOXB6 intergenic region within the great apes and humans

Data accumulated over the past decade from several loci suggest that nonhuman primates have a greater amount of intraspecific genetic variation relative to humans. In phylogenetic reconstructions among primates that are based on genetic data, therefore, it becomes essential to adequately sample the population(s) being analyzed. Inadequate sampling may not only underestimate variation within any given population, but such an underestimate may confound any phylogenetic or population-specific conclusions implied by the data. Here we present inter- and intraspecific data on the molecular evolution of an approximately 1.0 kb intergenic HOXB6 sequence among humans, common chimpanzees, pygmy chimpanzees, gorillas and orangutans. To date, this study represents the most comprehensive investigation of a noncoding nuclear locus among the great apes and humans that examines the nature and amount of intraspecific variation in DNA sequences. Not only do these HOXB6 data continue to support earlier findings that Homo sapiens sapiens has less genetic variation than any great ape species (Ruano et al., 1992; Deinard & Kidd, 1995), but they strongly suggest that a high level of genetic polymorphism existed within the common ancestor of the African ape clade (Homo-Pan-Gorilla). Despite detecting two nucleotide substitutions linking Pan and Homo, we caution against concluding that the HOXB6 data definitively support a Homo-Pan clade to the exclusion of Gorilla. Rather, we believe that taking into consideration the level of genetic polymorphism that is likely to have existed within the common ancestor, the most prudent conclusion that can be made from all available data, including morphological, karyotypic and genetic data, may be that speciation among Homo-Pan-Gorilla is best represented by a "trichotomy".     

Discovery of Novel Herpes Simplexviruses in Wild Gorillas,Bonobos, and Chimpanzees Supports Zoonotic Origin of HSV-2

Viruses closely related to human pathogens can reveal the origins of human infectious diseases. Human herpes simplexvirus type 1 (HSV-1) and type 2 (HSV-2) are hypothesized to have arisen via host-virus codivergence and cross-species transmission. We report the discovery of novel herpes simplexviruses during a large-scale screening of fecal samples from wild gorillas, bonobos, and chimpanzees. Phylogenetic analysis indicates that, contrary to expectation, simplexviruses from these African apes are all more closely related to HSV-2 than to HSV-1. Molecular clock-based hypothesis testing suggests the divergence between HSV-1 and the African great ape simplexviruses likely represents a codivergence event between humans and gorillas. The simplexviruses infecting African great apes subsequently experienced multiple cross-species transmission events over the past 3 My, the most recent of which occurred between humans and bonobos around 1 Ma. These findings revise our understanding of the origins of human herpes simplexviruses and suggest that HSV-2 is one of the earliest zoonotic pathogens.     

Assessing DNA sequence variations in human ESTs in a phylogenetic context using high-density oligonucleotide arrays

We have analyzed human genomic diversity in 32 individuals representing four continental populations of Homo sapiens in the context of four ape species. We used DNA resequencing chips covering 898 expressed sequence tags (ESTs), corresponding to 109 kb of sequence. Based on the intra-species data, the neutral hypothesis could not be rejected. However, the mutation rate was two times lower than typically observed in functionally unconstrained genomic segments, suggesting a certain level of selection. The worldwide diversity (297 segregating sites and nucleotide diversity of 0.054%) was partitioned among continents, with the greatest amount of variation observed in the African sample. The long-term effective population size of the human population was estimated at 13,000; a similar figure was obtained for the African sample and a 20% lower estimate was obtained for the other continents. Africans also differed in having a higher number of continental-specific polymorphisms contributing to the higher average nucleotide diversity. These results are consistent with the existence of two distinct lineages of modern humans: amalgamation of these lineages in Africa led to the higher present-day diversity on that continent, whereas colonization of other continents by one of them gave the effect of a population bottleneck.     

Likelihood and Bayes estimation of ancestral population sizes in hominoids using data from multiple loci

Polymorphisms in an ancestral population can cause conflicts between gene trees and the species tree. Such conflicts can be used to estimate ancestral population sizes when data from multiple loci are available. In this article I extend previous work for estimating ancestral population sizes to analyze sequence data from three species under a finite-site nucleotide substitution model. Both maximum-likelihood (ML) and Bayes methods are implemented for joint estimation of the two speciation dates and the two population size parameters. Both methods account for uncertainties in the gene tree due to few informative sites at each locus and make an efficient use of information in the data. The Bayes algorithm using Markov chain Monte Carlo (MCMC) enjoys a computational advantage over ML and also provides a framework for incorporating prior information about the parameters. The methods are applied to a data set of 53 nuclear noncoding contigs from human, chimpanzee, and gorilla published by Chen and Li. Estimates of the effective population size for the common ancestor of humans and chimpanzees by both ML and Bayes methods are approximately 12,000-21,000, comparable to estimates for modern humans, and do not support the notion of a dramatic size reduction in early human populations. Estimates published previously from the same data are several times larger and appear to be biased due to methodological deficiency. The divergence between humans and chimpanzees is dated at approximately 5.2 million years ago and the gorilla divergence 1.1-1.7 million years earlier. The analysis suggests that typical data sets contain useful information about the ancestral population sizes and that it is advantageous to analyze data of several species simultaneously.     

Consequences of non-intervention for infectious disease in African great apes

Infectious disease has recently joined poaching and habitat loss as a major threat to African apes. Both "naturally" occurring pathogens, such as Ebola and Simian Immunodeficiency Virus (SIV), and respiratory pathogens transmitted from humans, have been confirmed as important sources of mortality in wild gorillas and chimpanzees. While awareness of the threat has increased, interventions such as vaccination and treatment remain controversial. Here we explore both the risk of disease to African apes, and the status of potential responses. Through synthesis of published data, we summarize prior disease impact on African apes. We then use a simple demographic model to illustrate the resilience of a well-known gorilla population to disease, modeled on prior documented outbreaks. We found that the predicted recovery time for this specific gorilla population from a single outbreak ranged from 5 years for a low mortality (4%) respiratory outbreak, to 131 years for an Ebola outbreak that killed 96% of the population. This shows that mortality rates comparable to those recently reported for disease outbreaks in wild populations are not sustainable. This is particularly troubling given the rising pathogen risk created by increasing habituation of wild apes for tourism, and the growth of human populations surrounding protected areas. We assess potential future disease spillover risk in terms of vaccination rates amongst humans that may come into contact with wild apes, and the availability of vaccines against potentially threatening diseases. We discuss and evaluate non-interventionist responses such as limiting tourist access to apes, community health programs, and safety, logistic, and cost issues that constrain the potential of vaccination.     

Genetic variation in gorillas

This review summarizes what is currently known concerning genetic variation in gorillas, on both inter- and intraspecific levels. Compared to the human species, gorillas, along with the other great apes, possess greater genetic variation as a consequence of a demographic history of rather constant population size. Data and hence conclusions from analysis of mitochondrial DNA (mtDNA), the usual means of describing intraspecific patterns of genetic diversity, are limited at this time. An important task for future studies is to determine the degree of confidence with which gorilla mtDNA can be analyzed, in view of the risk that one will inadvertently analyze artifactual rather than genuine sequences. The limited information available from sequences of nuclear genomic segments does not distinguish western from eastern gorillas, and, in comparison with results from the two chimpanzee species, suggests a relatively recent common ancestry for all gorillas. In the near future, the greatest insights are likely to come from studies aimed at genetic characterization of all individual members of social groups. Such studies, addressing topics such as behavior of individuals with kin and non-kin, and the actual success of male reproductive strategies, will provide a link between behavioral and genetic studies of gorillas.     

Hypoglossal canal size in living hominoids and the evolution of human speech

The relative size of the hypoglossal canal has been proposed as a useful diagnostic tool for the identification of human-like speech capabilities in the hominid fossil record. Relatively large hypoglossal canals (standardized to oral cavity size) were observed in humans and assumed to correspond to relatively large hypoglossal nerves, the cranial nerve that controls motor function of the tongue. It was suggested that the human pattern of tongue motor innervation and associated speech potential are very different from those of African apes and australopithecines; the modern human condition apparently appeared by the time of Middle Pleistocene Homo. A broader interspecific analysis of hypoglossal canal size in primates conducted in 1999 has rejected this diagnostic and inferences based upon it. In an attempt to resolve these differences of opinion, which we believe are based in part on biased size-adjustments and/or unwarranted assumptions, a new data set was collected and analyzed from 298 extant hominoid skulls, including orangutans, gorillas, chimpanzees, bonobos, siamang, gibbons, and modern humans. Data on the absolute size of the hypoglossal nerve itself were also gathered from a small sample of humans and chimpanzee cadavers. A scale-free index of relative hypoglossal canal size (RHCS) was computed as 100 x (hypoglossal canal area(0.5)/oral cavity volume(0.333)). No significant sexual dimorphism in RHCS was discovered in any species of living hominoid, but there are significant interspecific differences in both absolute and relative sizes of the hypoglossal canal. In absolute terms, humans possess significantly larger canals than any other species except gorillas, but there is considerable overlap with chimpanzees. Humans are also characterized by large values of RHCS, but gibbons possess an even larger average mean for this index; siamang and bonobos overlap appreciably with humans in RHCS. The value of RHCS in Australopithecus afarensis is well within both human and gibbon ranges, as are the indices computed for selected representatives of fossil Homo. Furthermore, the size of the hypoglossal nerve itself, expressed as the mass of nerve per millimeter of length, does not distinguish chimpanzees from modern humans. We conclude, therefore, that the relative size of the hypoglossal canal is neither a reliable nor sufficient predictor of human-like speech capabilities, and paleoanthropology still lacks a quantifiable, morphological diagnostic for when this capability finally emerged in the human career.