DNA archives and our nearest relative: the trichotomy problem revisited

Knuckle-walking anteater: a convergence test of adaptation for purported knuckle-walking features of African Hominidae

Appeals to synapomorphic features of the wrist and hand in African apes, early hominins, and modern humans as evidence of knuckle-walking ancestry for the hominin lineage rely on accurate interpretations of those features as adaptations to knuckle-walking locomotion. Because Gorilla, Pan, and Homo share a relatively close common ancestor, the interpretation of such features is confounded somewhat by phylogeny. The study presented here examines the evolution of a similar locomotor regime in New World anteaters (order Xenarthra, family Myrmecophagidae) and uses the terrestrial giant anteater (Myrmecophaga tridactyla) as a convergence test of adaptation for purported knuckle-walking features of the Hominidae. During the stance phase of locomotion, Myrmecophaga transmits loads through flexed digits and a vertical manus, with hyperextension occurring at the metacarpophalangeal joints of the weight-bearing rays. This differs from the locomotion of smaller, arboreal anteaters of outgroup genera Tamandua and Cyclopes that employ extended wrist postures during above-branch quadrupedality. A number of features shared by Myrmecophaga and Pan and Gorilla facilitate load transmission or limit extension, thereby stabilizing the wrist and hand during knuckle-walking, and distinguish these taxa from their respective outgroups. These traits are a distally extended dorsal ridge of the distal radius, proximal expansion of the nonarticular surface of the dorsal capitate, a pronounced articular ridge on the dorsal aspects of the load-bearing metacarpal heads, and metacarpal heads that are wider dorsally than volarly. Only the proximal expansion of the nonarticular area of the dorsal capitate distinguishes knuckle-walkers from digitigrade cercopithecids, but features shared with digitigrade primates might be adaptive to the use of a vertical manus of some sort in the stance phase of terrestrial locomotion. The appearance of capitate nonarticular expansion and the dorsal ridge of the distal radius in the hominin lineage might be indicative of a knuckle-walking ancestry for bipedal hominins if interpreted within the biomechanical and phylogenetic context of hominid locomotor evolution.

More genes or more taxa? The relative contribution of gene number and taxon number to phylogenetic accuracy

The relative contribution of taxon number and gene number to accuracy in phylogenetic inference is a major issue in phylogenetics and of central importance to the choice of experimental strategies for the successful reconstruction of a broad sketch of the tree of life. Maximization of the number of taxa sampled is the strategy favored by most phylogeneticists, although its necessity remains the subject of debate. Vast increases in gene number are now possible due to advances in genomics, but large numbers of genes will be available for only modest numbers of taxa, raising the question of whether such genome-scale phylogenies will be robust to the addition of taxa. To examine the relative benefit of increasing taxon number or gene number to phylogenetic accuracy, we have developed an assay that utilizes the symmetric difference tree distance as a measure of phylogenetic accuracy. We have applied this assay to a genome-scale data matrix containing 106 genes from 14 yeast species. Our results show that increasing taxon number correlates with a slight decrease in phylogenetic accuracy. In contrast, increasing gene number has a significant positive effect on phylogenetic accuracy. Analyses of an additional taxon-rich data matrix from the same yeast clade show that taxon number does not have a significant effect on phylogenetic accuracy. The positive effect of gene number and the lack of effect of taxon number on phylogenetic accuracy are also corroborated by analyses of two data matrices from mammals and angiosperm plants, respectively. We conclude that, for typical data sets, the number of genes utilized may be a more important determinant of phylogenetic accuracy than taxon number.

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

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

Evolution of hominoids and the search for a genetic basis for creating humanness

The phylogenetic relationship of human and apes are reviewed. The history of molecular phylogenetic studies in this field is then discussed, as is the role of natural selection at the molecular level. It is argued that approximately 10,000 genetic changes are responsible for creating human specific phenotypes. A genome-wide comparison is necessary to decipher those changes.

Morphometrics and hominoid phylogeny: Support for a chimpanzee-human clade and differentiation among great ape subspecies

Taxonomic and phylogenetic analyses of great apes and humans have identified two potential areas of conflict between molecular and morphological data: phylogenetic relationships among living species and differentiation of great ape subspecies. Here we address these problems by using morphometric data. Three-dimensional landmark data from the hominoid temporal bone effectively quantify the shape of a complex element of the skull. Phylogenetic analysis using distance-based methods corroborates the molecular consensus on African ape and human phylogeny, strongly supporting a Pan-Homo clade. Phenetic differentiation of great ape subspecies is pronounced, as suggested previously by mitochondrial DNA and some morphological studies. These results show that the hominoid temporal bone contains a strong phylogenetic signal and reveal the potential for geometric morphometric analysis to shed light on phylogenetic relationships.

Human-specific amino acid changes found in 103 protein-coding genes

We humans have many characteristics that are different from those of the great apes. These human-specific characters must have arisen through mutations accumulated in the genome of our direct ancestor after the divergence of the last common ancestor with chimpanzee. Gene trees of human and great apes are necessary for extracting these human-specific genetic changes. We conducted a systematic analysis of 103 protein-coding genes for human, chimpanzee, gorilla, and orangutan. Nucleotide sequences for 18 genes were newly determined for this study, and those for the remaining genes were retrieved from the DDBJ/EMBL/GenBank database. The total number of amino acid changes in the human lineage was 147 for 26,199 codons (0.56%). The total number of amino acid changes in the human genome was, thus, estimated to be about 80,000. We applied the acceleration index test and Fisher's synonymous/nonsynonymous exact test for each gene tree to detect any human-specific enhancement of amino acid changes compared with ape branches. Six and two genes were shown to have significantly higher nonsynonymous changes at the human lineage from the acceleration index and exact tests, respectively. We also compared the distribution of the differences of the nonsynonymous substitutions on the human lineage and those on the great ape lineage. Two genes were more conserved in the ape lineage, whereas one gene was more conserved in the human lineage. These results suggest that a small proportion of protein-coding genes started to evolve differently in the human lineage after it diverged from the ape lineage.

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.

Inferring nonneutral evolution from human-chimp-mouse orthologous gene trios

Even though human and chimpanzee gene sequences are nearly 99% identical, sequence comparisons can nevertheless be highly informative in identifying biologically important changes that have occurred since our ancestral lineages diverged. We analyzed alignments of 7645 chimpanzee gene sequences to their human and mouse orthologs. These three-species sequence alignments allowed us to identify genes undergoing natural selection along the human and chimp lineage by fitting models that include parameters specifying rates of synonymous and nonsynonymous nucleotide substitution. This evolutionary approach revealed an informative set of genes with significantly different patterns of substitution on the human lineage compared with the chimpanzee and mouse lineages. Partitions of genes into inferred biological classes identified accelerated evolution in several functional classes, including olfaction and nuclear transport. In addition to suggesting adaptive physiological differences between chimps and humans, human-accelerated genes are significantly more likely to underlie major known Mendelian disorders.

Evolution of the Human ASPM Gene, a Major Determinant of Brain Size

The size of human brain tripled over a period of ∼2 million years (MY) that ended 0.2–0.4 MY ago. This evolutionary expansion is believed to be important to the emergence of human language and other high-order cognitive functions, yet its genetic basis remains unknown. An evolutionary analysis of genes controlling brain development may shed light on it. ASPM (abnormal spindle-like microcephaly associated) is one of such genes, as nonsense mutations lead to primary microcephaly, a human disease characterized by a 70% reduction in brain size. Here I provide evidence suggesting that human ASPM went through an episode of accelerated sequence evolution by positive Darwinian selection after the split of humans and chimpanzees but before the separation of modern non-Africans from Africans. Because positive selection acts on a gene only when the gene function is altered and the organismal fitness is increased, my results suggest that adaptive functional modifications occurred in human ASPM and that it may be a major genetic component underlying the evolution of the human brain.

Genome-scale approaches to resolving incongruence in molecular phylogenies

One of the most pervasive challenges in molecular phylogenetics is the incongruence between phylogenies obtained using different data sets, such as individual genes. To systematically investigate the degree of incongruence, and potential methods for resolving it, we screened the genome sequences of eight yeast species and selected 106 widely distributed orthologous genes for phylogenetic analyses, singly and by concatenation. Our results suggest that data sets consisting of single or a small number of concatenated genes have a significant probability of supporting conflicting topologies. By contrast, analyses of the entire data set of concatenated genes yielded a single, fully resolved species tree with maximum support. Comparable results were obtained with a concatenation of a minimum of 20 genes; substantially more genes than commonly used but a small fraction of any genome. These results have important implications for resolving branches of the tree of life.

Alu elements and hominid phylogenetics

Alu elements have inserted in primate genomes throughout the evolution of the order. One particular Alu lineage (Ye) began amplifying relatively early in hominid evolution and continued propagating at a low level as many of its members are found in a variety of hominid genomes. This study represents the first conclusive application of short interspersed elements, which are considered nearly homoplasy-free, to elucidate the phylogeny of hominids. Phylogenetic analysis of Alu Ye5 elements and elements from several other subfamilies reveals high levels of support for monophyly of Hominidae, tribe Hominini and subtribe Hominina. Here we present the strongest evidence reported to date for a sister relationship between humans and chimpanzees while clearly distinguishing the chimpanzee and human lineages.

Low genetic divergence obscures phylogeny among populations of Sphenodon, remnant of an ancient reptile lineage

Tuatara (two species of Sphenodon) are the last representatives of a branch of an ancient reptilian lineage, Sphenodontia, that have been isolated on the New Zealand landmass for 82 million years. We present analyses of geographic variation in allozymes, mitochondrial DNA, nuclear DNA sequences, and one-way albumin immunological comparisons. These all confirm a surprisingly low level of genetic diversity within Sphenodon for such an ancient lineage. We hypothesise a recent extended population bottleneck, probably during the Pliocene/Pleistocene glaciation cycles, to explain the current paucity of variation. All data sets reveal clear genetic differentiation between the northern populations and those in Cook Strait, but offer conflicting views of the history and taxonomic relationships of the Cook Strait population on North Brother Island, currently recognised as Sphenodon guntheri. Allozymes show this population to be the most divergent of all tuatara populations, but preliminary mitochondrial DNA data indicate few differences between S. guntheri and Cook Strait Sphenodon punctatus. Interpretation of the trees is confounded by the lack of a suitable outgroup. As in other cases of conflicting nuclear and mitochondrial data sets, the different data sets likely reveal different aspects of the animals' evolutionary history, and introgression is not uncommon between species pairs.

Human-specific regulation of alpha 2-6-linked sialic acids

Many microbial pathogens and toxins recognize animal cells via cell surface sialic acids (Sias) that are alpha 2-3- or alpha 2-8-linked to the underlying glycan chain. Human influenza A/B viruses are unusual in preferring alpha 2-6-linked Sias, undergoing a switch from alpha 2-3 linkage preference during adaptation from animals to humans. This correlates with the expression of alpha 2-6-linked Sias on ciliated human airway epithelial target cells and of alpha 2-3-linked Sias on secreted soluble airway mucins, which are unable to inhibit virus binding. Given several known differences in Sia biology between humans and apes, we asked whether this pattern of airway epithelial Sia linkages is also human-specific. Indeed, we show that since the last common ancestor with apes, humans underwent a concerted bidirectional switch in alpha 2-6-linked Sia expression between airway epithelial cell surfaces and secreted mucins. This can explain why the chimpanzee appears relatively resistant to experimental infection with human Influenza viruses. Other tissues showed additional examples of human-specific increases or decreases in alpha 2-6-linked Sia expression and only one example of a change specific to certain great apes. Furthermore, while human and great ape leukocytes both express alpha 2-6-linked Sias, only human erythrocytes have markedly up-regulated expression. These cell type-specific changes in alpha 2-6-Sia expression during human evolution represent another example of a human-specific change in Sia biology. Because the data set involves multiple great apes, we can also conclude that Sia linkage expression patterns can be conserved during millions of years of evolution within some vertebrate taxa while undergoing sudden major changes in other closely related ones.

Molecular systematics of armadillos (Xenarthra, Dasypodidae): contribution of maximum likelihood and Bayesian analyses of mitochondrial and nuclear genes

The 30 living species of armadillos, anteaters, and sloths (Mammalia: Xenarthra) represent one of the three major clades of placentals. Armadillos (Cingulata: Dasypodidae) are the earliest and most speciose xenarthran lineage with 21 described species. The question of their tricky phylogeny was here studied by adding two mitochondrial genes (NADH dehydrogenase subunit 1 [ND1] and 12S ribosomal RNA [12S rRNA]) to the three protein-coding nuclear genes (alpha2B adrenergic receptor [ADRA2B], breast cancer susceptibility exon 11 [BRCA1], and von Willebrand factor exon 28 [VWF]) yielding a total of 6869 aligned nucleotide sites for thirteen xenarthran species. The two mitochondrial genes were characterized by marked excesses of transitions over transversions-with a strong bias toward CT transitions for the 12S rRNA-and exhibited two- to fivefold faster evolutionary rates than the fastest nuclear gene (ADRA2B). Maximum likelihood and Bayesian phylogenetic analyses supported the monophyly of Dasypodinae, Tolypeutinae, and Euphractinae, with the latter two armadillo subfamilies strongly clustering together. Conflicting branching points between individual genes involved relationships within the subfamilies Tolypeutinae and Euphractinae. Owing to a greater number of informative sites, the overall concatenation favored the mitochondrial topology with the classical grouping of Cabassous and Priodontes within Tolypeutinae, and a close relationship between Euphractus and Chaetophractus within Euphractinae. However, low statistical support values associated with almost equal distributions of apomorphies among alternatives suggested that two parallel events of rapid speciation occurred within these two armadillo subfamilies.

Divergence of the genes on human chromosome 21 between human and other hominoids and variation of substitution rates among transcription units

The study of genomic divergence between humans and primates may provide insight into the origins of human beings and the genetic basis of unique human traits and diseases. Chromosome 21 is the smallest chromosome in the human genome, and some of its regions have been implicated in mental retardation and other diseases. In this study, we sequenced the coding and regulatory regions of 127 known genes on human chromosome 21 in DNA samples from human and chimpanzees and a part of the corresponding genes from orangutan, gorilla, and macaque. Overall, 3,003 nucleotide differences between human and chimpanzee were identified over approximately 400 kb. The differences in coding, promoter, and exon-intron junction regions were 0.51 +/- 0.02%, 0.88 +/- 0.03%, and 0.85 +/- 0.02%, respectively, much lower than the previously reported 1.23% in genomic regions, which suggests the presence of purifying selection. Significant variation in substitution rate among genes was observed by comparing the divergence between human and chimpanzee. Furthermore, by implementing a bioinformatics-based approach, we showed that the identification of genetic variants specific to the human lineage might lead to an understanding of the mechanisms that are attributable to the phenotypes that unique to humans, by changing the structure and/or dosage of the proteins expressed. A phylogenetic analysis unambiguously confirms the conclusion that chimpanzees were our closest relatives to the exclusion of other primates and the relative divergence of the Homo-Pan and that of (Homo-Pan)-Gorilla are 4.93 million years and 7.26 million years, respectively.

The application of molecular genetic approaches to the study of human evolution

The past decade of advances in molecular genetic technology has heralded a new era for all evolutionary studies, but especially the science of human evolution. Data on various kinds of DNA variation in human populations have rapidly accumulated. There is increasing recognition of the importance of this variation for medicine and developmental biology and for understanding the history of our species. Haploid markers from mitochondrial DNA and the Y chromosome have proven invaluable for generating a standard model for evolution of modern humans. Conclusions from earlier research on protein polymorphisms have been generally supported by more sophisticated DNA analysis. Co-evolution of genes with language and some slowly evolving cultural traits, together with the genetic evolution of commensals and parasites that have accompanied modern humans in their expansion from Africa to the other continents, supports and supplements the standard model of genetic evolution. The advances in our understanding of the evolutionary history of humans attests to the advantages of multidisciplinary research.

Large differences between LINE-1 amplification rates in the human and chimpanzee lineages

The genomic evolution and causes of phenotypic variation among humans and great apes remain largely unknown, although the phylogenetic relationships among them have been extensively explored. Previous studies that focus on differences at the amino acid and nucleotide sequence levels have revealed a high degree of similarity between humans and chimpanzees, suggesting that other types of genomic change may have contributed to the relatively large phenotypic differences between them. For example, the activity of long interspersed element 1 (LINE-1) retrotransposons may impose significant changes on genomic structure and function and, consequently, on phenotype. Here we investigate the relative rates of LINE-1 amplification in the lineages leading to humans, bonobos (Pan paniscus), and chimpanzees (P. troglodytes). Our data indicate that LINE-1 insertions have accumulated at significantly greater rates in bonobos and chimpanzees than in humans, provide insights into the timing of major LINE-1 amplification events during great ape evolution, and identify a Pan-specific LINE-1 subfamily.

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.

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.

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.

Contribution of homoplasy and of ancestral polymorphism to the evolution of genes in anthropoid primates

Molecular phylogenies of lineages that split from one another in short succession are often difficult to resolve because different loci and different sites within the same locus yield incongruent relationships. The incongruity is commonly attributed to two causes: differential assortment of ancestral polymorphisms and homoplasy. To assess the relative contribution of these two causes, sequences of 57 segments from 51 loci in six primate lineages (human, chimpanzee, gorilla, orangutan, macaque, and tamarin, abbreviated as H, C, G, O, M, and T, respectively) were subjected to "partitioning" analysis, in which phylogenetically informative sites were identified in all 15 pairwise comparisons of each of the 57 segments and tallied for their support or lack thereof for each of the theoretically possible phylogenies. The six lineages include one of the best known cases of a difficult-to-resolve phylogeny: the trichotomy (H, C, G), in which the three lineages may have diverged from each other within a short period of time. In this period many of the ancestral polymorphisms apparently persisted and yielded phylogenetically incongruent signals. By contrast, no ancestral polymorphism is expected to have survived during the interval separating the divergences of the O, M, and T lineages from the ancestor of the (H, C, G) group. Any phylogenetic incompatibilities at sites in the O, M, and T lineages relative to the (H, C, G) group are therefore presumably the result of homoplasy. The frequency of homoplasy estimated in this manner is unexpectedly high: 12% for the (H, C, G) clade and 19% for the (H, C, G, O) clade. At least three-quarters of the 48% incompatibility observed in the (H, C) clade is attributable to the sorting out of ancestral polymorphisms coupled with intragenic recombination. Possible reasons for this high level of homoplasy in the O, M, and T lineages are discussed, and a computer simulation has been carried out to produce a model explaining the observed data.

Direct evidence for the Homo-Pan clade

For a long time, the evolutionary relationship between human and African apes, the 'trichotomy problem', has been debated with strong differences in opinion and interpretation. Statistical analyses of different molecular DNA data sets have been carried out and have primarily supported a Homo-Pan clade. An alternative way to address this question is by the comparison of evolutionarily relevant chromosomal breakpoints. Here, we made use of a P1-derived artificial chromosome (PAC)/bacterial artificial chromosome (BAC) contig spanning approximately 2.8 Mb on the long arm of the human Y chromosome, to comparatively map individual PAC clones to chromosomes from great apes, gibbons, and two species of Old World monkeys by fluorescence in-situ hybridization. During our search for evolutionary breakpoints on the Y chromosome, it transpired that a transposition of an approximately 100-kb DNA fragment from chromosome 1 onto the Y chromosome must have occurred in a common ancestor of human, chimpanzee and bonobo. Only the Y chromosomes of these three species contain the chromosome-1-derived fragment; it could not be detected on the Y chromosomes of gorillas or the other primates examined. Thus, this shared derived (synapomorphic) trait provides clear evidence for a Homo-Pan clade independent of DNA sequence analysis.

Loss of urate oxidase activity in hominoids and its evolutionary implications

We have determined and compared the promoter, coding, and intronic sequences of the urate oxidase (Uox) gene of various primate species. Although we confirm the previous observation that the inactivation of the gene in the clade of the human and the great apes results from a single CGA to TGA nonsense mutation in exon 2, we find that the inactivation in the gibbon lineage results from an independent nonsense mutation at a different CGA codon in exon 2 or from either one-base deletion in exon 3 or one-base insertion in exon 5, contrary to the previous claim that the cause is a 13-bp deletion in exon 2. We also find that compared with other organisms, the primate functional Uox gene is exceptional in terms of usage of CGA codons which are prone to TGA nonsense mutations. Nevertheless, we demonstrate rather strong selective constraint against nonsynonymous sites of the functional Uox gene and argue that this observation is consistent with the fact that the Uox gene is unique in the genome and evolutionarily conserved not only among animals but also among eukaryotes. Another finding that there are a few substitutions in the cis-acting element or CAAT-box (or both) of primate functional Uox genes may explain the lowered transcriptional activity. We suggest that although the inactivation of the hominoid Uox gene was caused by independent nonsense or frameshift mutations, the gene has taken a two-step deterioration process, first in the promoter and second in the coding region during primate evolution. It is also argued that the high concentration of uric acid in the blood of humans and nonhuman primates has developed molecular coevolution with the xanthine oxidoreductase in purine metabolism. However, it remains to be answered whether loss of Uox activity in hominoids is related to protection from oxidative damage and the prolonged life span.

The effects of nucleotide substitution model assumptions on estimates of nonparametric bootstrap support

The use of parameter-rich substitution models in molecular phylogenetics has been criticized on the basis that these models can cause a reduction both in accuracy and in the ability to discriminate among competing topologies. We have explored the relationship between nucleotide substitution model complexity and nonparametric bootstrap support under maximum likelihood (ML) for six data sets for which the true relationships are known with a high degree of certainty. We also performed equally weighted maximum parsimony analyses in order to assess the effects of ignoring branch length information during tree selection. We observed that maximum parsimony gave the lowest mean estimate of bootstrap support for the correct set of nodes relative to the ML models for every data set except one. For several data sets, we established that the exact distribution used to model among-site rate variation was critical for a successful phylogenetic analysis. Site-specific rate models were shown to perform very poorly relative to gamma and invariable sites models for several of the data sets most likely because of the gross underestimation of branch lengths. The invariable sites model also performed poorly for several data sets where this model had a poor fit to the data, suggesting that addition of the gamma distribution can be critical. Estimates of bootstrap support for the correct nodes often increased under gamma and invariable sites models relative to equal rates models. Our observations are contrary to the prediction that such models cause reduced confidence in phylogenetic hypotheses. Our results raise several issues regarding the process of model selection, and we briefly discuss model selection uncertainty and the role of sensitivity analyses in molecular phylogenetics.

The probability of topological concordance of gene trees and species trees

The concordance of gene trees and species trees is reconsidered in detail, allowing for samples of arbitrary size to be taken from the species. A sense of concordance for gene tree and species tree topologies is clarified, such that if the "collapsed gene tree" produced by a gene tree has the same topology as the species tree, the gene tree is said to be topologically concordant with the species tree. The term speciodendric is introduced to refer to genes whose trees are topologically concordant with species trees. For a given three-species topology, probabilities of each of the three possible collapsed gene tree topologies are given, as are probabilities of monophyletic concordance and concordance in the sense of N. Takahata (1989), Genetics 122, 957-966. Increasing the sample size is found to increase the probability of topological concordance, but a limit exists on how much the topological concordance probability can be increased. Suggested sample sizes beyond which this probability can be increased only minimally are given. The results are discussed in terms of implications for molecular studies of phylogenetics and speciation.

Genome of the apes

The Human Genome Project has generated both the information and technological infrastructure needed to accelerate genetic comparisons between humans and the African great apes (chimpanzees and gorillas). Sequence and chromosomal organization differences between these highly related genomes will provide clues to the genetic basis for recently evolved, specifically human traits such as bipedal gait and advanced cognitive function. Recent studies comparing the primate genomes have the potential to affect many aspects of human biomedical research and could benefit primate conservation efforts.

Alu-mediated inactivation of the human CMP- N-acetylneuraminic acid hydroxylase gene

Inactivation of the CMP-N-acetylneuraminic acid hydroxylase gene has provided an example of human-specific genomic mutation that results in a widespread biochemical difference between human and nonhuman primates. We have found that, although a region containing a 92-bp exon and an AluSq element in the hydroxylase gene is intact in all nonhuman primates examined, the same region in the human genome is replaced by an AluY element that was disseminated at least one million years ago. We propose a mechanistic model for this Alu-mediated replacement event, which deleted the 92-bp exon and thus inactivated the human hydroxylase gene. It is suggested that Alu elements have played potentially important roles in genotypic and phenotypic evolution in the hominid lineage.

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.

Mutations of the 3' untranslated region of the SDF1 gene in apes and monkeys: potential impact on sensitivity to AIDS induced by lentiviruses

The comparison of the stromal cell-derived factor-1 (SDF1) gene 3' untranslated region (3'UTR) of four great ape and four monkey species with their human counterparts shows that the human SDF1-3'A mutation is present in primate species that are the most susceptible to lentivirus-induced AIDS and is absent in species that are particularly resistant to lentivirus-induced AIDS. The results enlighten the possible relationship between SDF1-3'UTR polymorphism and sensitivity to AIDS.

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.

Cytochrome b phylogeny and the taxonomy of great apes and mammals

In the Linnaean system of classification, the generic status of a species is part of its binomial name, and it is therefore important that the classification at the level of genus is consistent at least in related groups of organisms. Using maximum-likelihood phylogenetic trees constructed from a large number of complete or nearly complete mammalian cytochrome b sequences, I show that the distributions of intrageneric and intergeneric distances derived from these trees are clearly separated, which allows the limits for a more rational generic classification of mammals to be established. The analysis of genetic distances among hominids in this context provides strong support for the inclusion of humans and chimpanzees in the same genus. It is also of interest to decipher the main reasons for the possible biases in the mammalian classification. I found by correlation analysis that the classification of mammals of large body size tends to be oversplit, whereas that of small mammals has an excess of lumping, which may be a manifestation of the larger difficulty in finding diagnostic characters in the classification of small animals. In addition, and contrary to some previous observations, there is no correlation between body size and rate of cytochrome b evolution in mammals, which excludes the difference in evolutionary rates as the cause of the observed body size taxonomic bias.

Genomic divergences between humans and other hominoids and the effective population size of the common ancestor of humans and chimpanzees

To study the genomic divergences among hominoids and to estimate the effective population size of the common ancestor of humans and chimpanzees, we selected 53 autosomal intergenic nonrepetitive DNA segments from the human genome and sequenced them in a human, a chimpanzee, a gorilla, and an orangutan. The average sequence divergence was only 1.24% +/- 0.07% for the human-chimpanzee pair, 1.62% +/- 0.08% for the human-gorilla pair, and 1.63% +/- 0.08% for the chimpanzee-gorilla pair. These estimates, which were confirmed by additional data from GenBank, are substantially lower than previous ones, which included repetitive sequences and might have been based on less-accurate sequence data. The average sequence divergences between orangutans and humans, chimpanzees, and gorillas were 3.08% +/- 0.11%, 3.12% +/- 0.11%, and 3.09% +/- 0.11%, respectively, which also are substantially lower than previous estimates. The sequence divergences in other regions between hominoids were estimated from extensive data in GenBank and the literature, and Alus showed the highest divergence, followed in order by Y-linked noncoding regions, pseudogenes, autosomal intergenic regions, X-linked noncoding regions, synonymous sites, introns, and nonsynonymous sites. The neighbor-joining tree derived from the concatenated sequence of the 53 segments--24,234 bp in length--supports the Homo-Pan clade with a 100% bootstrap value. However, when each segment is analyzed separately, 22 of the 53 segments (approximately 42%) give a tree that is incongruent with the species tree, suggesting a large effective population size (N(e)) of the common ancestor of Homo and Pan. Indeed, a parsimony analysis of the 53 segments and 37 protein-coding genes leads to an estimate of N(e) = 52,000 to 96,000. As this estimate is 5 to 9 times larger than the long-term effective population size of humans (approximately 10,000) estimated from various genetic polymorphism data, the human lineage apparently had experienced a large reduction in effective population size after its separation from the chimpanzee lineage. Our analysis assumes a molecular clock, which is in fact supported by the sequence data used. Taking the orangutan speciation date as 12 to 16 million years ago, we obtain an estimate of 4.6 to 6.2 million years for the Homo-Pan divergence and an estimate of 6.2 to 8.4 million years for the gorilla speciation date, suggesting that the gorilla lineage branched off 1.6 to 2.2 million years earlier than did the human-chimpanzee divergence.

Genetic differences between humans and great apes

The remarkable similarity among the genomes of humans and the African great apes could warrant their classification together as a single genus. However, whereas there are many similarities in the biology, life history, and behavior of humans and great apes, there are also many striking differences that need to be explained. The complete sequencing of the human genome creates an opportunity to ask which genes are involved in those differences. A logical approach would be to use the chimpanzee genome for comparison and the other great ape genomes for confirmation. Until such a great ape genome project can become reality, the next best approach must be educated guesses of where the genetic differences may lie and a careful analysis of differences that we do know about. Our group recently discovered a human-specific inactivating mutation in the CMP-sialic acid hydroxylase gene, which results in the loss of expression of a common mammalian cell-surface sugar throughout all cells in the human body. We are currently investigating the implications of this difference for a variety of issues relevant to humans, ranging from pathogen susceptibility to brain development. Evaluating the uniqueness of this finding has also led us to explore the existing literature on the broader issue of genetic differences between humans and great apes. The aim of this brief review is to consider a listing of currently known genetic differences between humans and great apes and to suggest avenues for future research. The differences reported between human and great ape genomes include cytogenetic differences, differences in the type and number of repetitive genomic DNA and transposable elements, abundance and distribution of endogenous retroviruses, the presence and extent of allelic polymorphisms, specific gene inactivation events, gene sequence differences, gene duplications, single nucleotide polymorphisms, gene expression differences, and messenger RNA splicing variations. Evaluation of the reported findings in all these categories indicates that the CMP-sialic hydroxylase mutation is the only one that has so far been shown to result in a global biochemical and structural difference between humans and great apes. Several of the other known genetic dissimilarities deserve more exploration at the functional level. Among the areas of focus for the future should be genes affecting development, mental maturation, reproductive biology, and other aspects of life history. The approaches taken should include both going from the genome up to the adaptive potential of the organisms and going from novel adaptive regimes down to the relevant repercussions in the genome. Also, as much as we desire a simple genetic explanation for the human phenomenon, it is much more probable that our evolution occurred in multiple genetic steps, many of which must have left detectable footprints in our genomes. Ultimately, we need to know the exact number of genetic steps, the order in which they occurred, and the temporal, spatial, environmental, and cultural contexts that determined their impact on human evolution.

Catarrhine phylogeny: noncoding DNA evidence for a diphyletic origin of the mangabeys and for a human-chimpanzee clade

Maximum-parsimony and maximum-likelihood analyses of two of the serum albumin gene's intron sequences from 24 catarrhines (17 cercopithecid and 7 hominid) and 3 platyrrhines (an outgroup to the catarrhines) yielded results on catarrhine phylogeny that are congruent with those obtained with noncoding sequences of the gamma(1)-gamma(2) globin gene genomic region, using only those flanking and intergenic gamma sequences that in their history were not involved in gene conversion. A data set that combined in a tandem alignment these two sets of noncoding DNA orthologues from the two unlinked nuclear genomic loci yielded the following confirmatory results both on the course of cladistic branchings (the divisions in a cladistic classification of higher ranking taxa into subordinate taxa) and on the ages of the taxa (each taxon representing a clade). The cercopithecid branch of catarrhines, at approximately 14 Ma (mega annum) divided into Colobini (the leaf-eating Old World monkeys) and Cercopithecini (the cheek-pouched Old World monkeys). At approximately 10-9 Ma, Colobini divided into an African clade, Colobina, and an Asian clade, Presbytina; similarly at this time level, Cercopithecini divided into Cercopithecina (the guenons, patas, and green monkeys) and Papionina. At approximately 7 Ma, Papionina divided into Macaca, Cercocebus, and Papio. At approximately 5 Ma, Cercocebus divided subgenerically into C. (Cercocebus) for terrestrial mangabeys and C. (Mandrillus) for drills and mandrills, while at approximately 4 Ma Papio divided subgenerically into P. (Locophocebus) for arboreal mangabeys, P. (Theropithecus) for gelada baboons, and P. (Papio) for hamadryas baboons. In turn, the hominid branch of catarrhines at approximately 18 Ma divided into Hylobatini (gibbons and siamangs) and Hominini; at approximately 14 Ma, Hominini divided into Pongina (orangutans) and Hominina; at approximately 7 Ma, Hominina divided into Gorilla and Homo; and at approximately 6-5 Ma, Homo divided subgenerically into H. (Homo) for humans and H. (Pan) for common and bonobo chimpanzees. Rates of noncoding DNA evolution were assessed using a data set of noncoding gamma sequence orthologues that represented 18 catarrhines, 16 platyrrhines, 3 non-anthropoid primates (2 tarsiers and 1 strepsirhine), and rabbit (as outgroup to the primates). Results obtained with this data set revealed a faster rate of nucleotide substitutions in the early primate lineage to the anthropoid (platyrrhine/catarrhine) ancestor than from that ancestor to the present. Rates were slower in catarrhines than in platyrrhines, slower in the cheek-pouched than in the leaf-eating cercopithecids, and slower yet in the hominids. On relating these results to data on brain sizes and life spans, it was suggested that life-history strategies that favor intelligence and longer life spans also select for decreases in de novo mutation rates.