Nucleotide sequences of immunoglobulin-epsilon pseudogenes in man and apes and their phylogenetic relationships

To understand the phylogenetic relationships between hominoids, the nucleotide sequences of immunoglobulin-epsilon processed pseudogenes from chimpanzee, gorilla and orangutan were determined. The basic structures of these processed pseudogenes agreed with their human counterpart. Although the degrees of nucleotide differences between man and the African apes had no statistical significance, all the analytical data examined supported the theory that chimpanzee is the closest relative of man. This result was consistent with that deduced by our recent qualitative study. Studies on the nucleotide sequences of globin genes have suggested that the molecular clock runs more slowly in hominoids than in non-hominoid primates. According to the present data, however, further retardation of the evolutionary rate was not observed in the human lineage. Assuming that orangutan diverged 14 million years ago and that the evolutionary rate between the orangutan lineage and the lineage leading to the other three species is constant, the divergence dates of chimpanzee and gorilla were estimated to be 4.9(+/- 0.9) and 5.9(+/- 0.9) million years ago, respectively.     

Human and ape molecular clocks and constraints on paleontological hypotheses

Although the relationships of the living hominoid primates (humans and apes) are well known, the relationships of the fossil species, times of divergence of both living and fossil species, and the biogeographic history of hominoids are not well established. Divergence times of living species, estimated from molecular clocks, have the potential to constrain hypotheses of the relationships of fossil species. In this study, new DNA sequences from nine protein-coding nuclear genes in great apes are added to existing datasets to increase the precision of molecular time estimates bearing on the evolutionary history of apes and humans. The divergence of Old World monkeys and hominoids at the Oligocene-Miocene boundary (approximately 23 million years ago) provides the best primate calibration point and yields a time and 95% confidence interval of 5.4 +/- 1.1 million years ago (36 nuclear genes) for the human-chimpanzee divergence. Older splitting events are estimated as 6.4 +/- 1.5 million years ago (gorilla, 31 genes), 11.3 +/- 1.3 million years ago (orangutan, 33 genes), and 14.9 +/- 2.0 million years ago (gibbon, 27 genes). Based on these molecular constraints, we find that several proposed phylogenies of fossil hominoid taxa are unlikely to be correct.     

Gorilla and orangutan c-myc nucleotide sequences: Inference on hominoid phylogeny

The nucleotide sequences of the gorilla and orangutan myc loci have been determined by the dideoxy nucleotide method. As previously observed in the human and chimpanzee sequences, an open reading frame (ORF) of 188 codons overlapping exon 1 could be deduced from the gorilla sequence. However, no such ORF appeared in the orangutan sequence.The two sequences were aligned with those of human and chimpanzee as hominoids and of gibbon and marmoset as outgroups of hominoids. The branching order in the evolution of primates was inferred from these data by different methods: maximum parsimony and neighborjoining.Our results support the view that the gorilla lineage branched off before the human and chimpanzee diverged and strengthen the hypothesis that chimpanzee and gorilla are more related to human than is orangutan. Correspondence to: F. Galibert     

Reexamination of the African hominoid trichotomy with additional sequences from the primate beta-globin gene cluster.

Additional DNA sequence information from a range of primates, including 13.7 kb from pygmy chimpanzee (Pan paniscus), was added to data sets of beta-globin gene cluster sequence alignments that span the gamma 1, gamma 2, and psi eta loci and their flanking and intergenic regions. This enlarged body of data was used to address the issue of whether the ancestral separations of gorilla, chimpanzee, and human lineages resulted from only one trichotomous branching or from two dichotomous branching events. The degree of divergence, corrected for superimposed substitutions, seen in the beta-globin gene cluster between human alleles is about a third to a half that observed between two species of chimpanzee and about a fourth that between human and chimpanzee. The divergence either between chimpanzee and gorilla or between human and gorilla is slightly greater than that between human and chimpanzee, suggesting that the ancestral separations resulted from two closely spaced dichotomous branchings. Maximum parsimony analysis further strengthened the evidence that humans and chimpanzees share the longest common ancestry. Support for this human-chimpanzee clade is statistically significant at P = 0.002 over a human-gorilla clade or a chimpanzee-gorilla clade. An analysis of expected and observed homoplasy revealed that the number of sequence changes uniquely shared by human and chimpanzee lineages is too large to be attributed to homoplasy. Molecular clock calculations that accommodated lineage variations in rates of molecular evolution yielded hominoid branching times that ranged from 17-19 million years ago (MYA) for the separation of gibbon from the other hominoids to 5-7 MYA for the separation of chimpanzees from humans. Based on the relatively late dates and mounting corroborative evidence from unlinked nuclear genes and mitochondrial DNA for the close sister grouping of humans and chimpanzees, a cladistic classification would place all apes and humans in the same family. Within this family, gibbons would be placed in one subfamily and all other extant hominoids in another subfamily. The later subfamily would be divided into a tribe for orangutans and another tribe for gorillas, chimpanzees, and humans. Finally, gorillas would be placed in one subtribe with chimpanzees and humans in another, although this last division is not as strongly supported as the other divisions.     

A comparison of TSPY genes from Y-chromosomal DNA of the great apes and humans: Sequence,evolution, and phylogeny

The genes for testis-specific protein Y (TSPY) were sequenced from chimpanzee (Pan troglodytes), gorilla (Gorilla gorilla), orangutan (Pongo pygmaeus), and baboon (Papio hamadryas). The sequences were compared with each other and with the published human sequence. Substitutions were detected at 144 of the 755 nucleotide positions compared. In overviewing five sequences, one deletion in human, four successive nucleotide insertions in orangutan, and seven deletions/insertions in baboon sequence were noted. The present sequences differed from that of human by 1.9% (chimpanzee), 4.0% (gorilla), 8.2% (orangutan), and 16.8% (baboon), respectively. The phylogenetic tree constructed by the neighbor-joining method suggests that human and chimpanzee are more closely related to each other than either of them is to gorilla, and this result is also supported by maximum likelihood and strict consensus maximum parsimony trees. The number of nucleotide substitutions per site between human and chimpanzee, gorilla, and orangutan for TSPY intron were 0.024, 0.048, and 0.094, respectively. The rates of nucleotide substitutions per site per year were higher in the TSPY intron than in the TSPY exon, and higher in the TSPY intron than in the ZFY (Zinc Finger Y) intron in human and apes. © 1996 Wiley-Liss, Inc.     

Estimation of hominoid ancestral population sizes under bayesian coalescent models incorporating mutation rate variation and sequencing errors

Estimation of population parameters for the common ancestors of humans and the great apes is important in understanding our evolutionary history. In particular, inference of population size for the human-chimpanzee common ancestor may shed light on the process by which the 2 species separated and on whether the human population experienced a severe size reduction in its early evolutionary history. In this study, the Bayesian method of ancestral inference of Rannala and Yang (2003. Bayes estimation of species divergence times and ancestral population sizes using DNA sequences from multiple loci. Genetics. 164:1645-1656) was extended to accommodate variable mutation rates among loci and random species-specific sequencing errors. The model was applied to analyze a genome-wide data set of approximately 15,000 neutral loci (7.4 Mb) aligned for human, chimpanzee, gorilla, orangutan, and macaque. We obtained robust and precise estimates for effective population sizes along the hominoid lineage extending back approximately 30 Myr to the cercopithecoid divergence. The results showed that ancestral populations were 5-10 times larger than modern humans along the entire hominoid lineage. The estimates were robust to the priors used and to model assumptions about recombination. The unusually low X chromosome divergence between human and chimpanzee could not be explained by variation in the male mutation bias or by current models of hybridization and introgression. Instead, our parameter estimates were consistent with a simple instantaneous process for human-chimpanzee speciation but showed a major reduction in X chromosome effective population size peculiar to the human-chimpanzee common ancestor, possibly due to selective sweeps on the X prior to separation of the 2 species.     

DNA hybridization evidence of hominoid phylogeny: Results from an expanded data set

Summary The living hominoids are human, the two species of chimpanzees, gorilla, orangutan, and nine species of gibbons. The cercopithecoids (Old World monkeys) are the sister group of the hominoids. A consensus about the phylogeny of the hominoids has been reached for the branching order of the gibbons (earliest) and the orangutan (next earliest), but the branching order among gorilla, chimpanzees, and human remains in contention. In 1984 we presented DNA-DNA hybridization data, based on 183 DNA hybrids, that we interpreted as evidence that the branching order, from oldest to most recent, was gibbons, orangutan, gorilla, chimpanzees, and human. In the present paper we report on an expanded data set totaling 514 DNA hybrids, which supports the branching order given above. The ranges for the datings of divergence nodes are Old World monkeys, 25–34 million years (Myr) ago; gibbons, 16.4–23 Myr ago; orangutan, 12.2–17 Myr ago; gorilla, 7.7–11 Myr ago; chimpanzees-human, 5.5–7.7 Myr ago. The possible effects of differences in age at first breeding are discussed, and some speculations about average genomic rates of evolution are presented.     

Man's place in hominoidea revealed by mitochondrial DNA genealogy

Summary Molecular biology has resurrected C. Darwin and T.H. Huxley's question about the origin of humans, but the precise branching pattern and dating remain controversial. To settle this issue, a large amount of sequence information is required. We determined mitochondrial (mt) DNA sequences for five hominoids; pygmy and common chimpanzees, gorilla, orangutan, and siamang. The common region compared with the known human sequence is 4759 by long, encompassing genes for 11 transfer RNAs and 6 proteins. Because of the high substitution rates in mammalian mtDNA and an unprecedentedly large region compared, the sequence differences clearly indicate that the closest relatives to human are chimpanzees rather than gorilla. For dating the divergences of human, chimpanzee, and gorilla, we used only unsaturated parts of sequence differences in which the mtDNA genealogy is not obscured by multiple substitutions. The result suggests that gorilla branched off 7.7 ± 0.7 million years (Myr) ago and human 4.7 ± 0.5 Myr ago; the time difference between these divergences being as long as 3 Myr.Offprint requests to: S. Horai     

Alu recombination-mediated structural deletions in the chimpanzee genome

With more than 1.2 million copies, Alu elements are one of the most important sources of structural variation in primate genomes. Here, we compare the chimpanzee and human genomes to determine the extent of Alu recombination-mediated deletion (ARMD) in the chimpanzee genome since the divergence of the chimpanzee and human lineages (~6 million y ago). Combining computational data analysis and experimental verification, we have identified 663 chimpanzee lineage-specific deletions (involving a total of ~771 kb of genomic sequence) attributable to this process. The ARMD events essentially counteract the genomic expansion caused by chimpanzee-specific Alu inserts. The RefSeq databases indicate that 13 exons in six genes, annotated as either demonstrably or putatively functional in the human genome, and 299 intronic regions have been deleted through ARMDs in the chimpanzee lineage. Therefore, our data suggest that this process may contribute to the genomic and phenotypic diversity between chimpanzees and humans. In addition, we found four independent ARMD events at orthologous loci in the gorilla or orangutan genomes. This suggests that human orthologs of loci at which ARMD events have already occurred in other nonhuman primate genomes may be “at-risk” motifs for future deletions, which may subsequently contribute to human lineage-specific genetic rearrangements and disorders.     

Molecular phylogeny of the family of apes and humans

The morphological picture of primate phylogeny has not unambiguously identified the nearest outgroup of Anthropoidea and has not resolved the branching pattern within Hominoidea. The molecular picture provides more resolution and clarifies the systematics of Hominoidea. Protein and DNA evidence divides Hominoidea into Hylobatidae (gibbons) and Hominidae, family Hominidae into Ponginae (orangutan) and Homininae, and subfamily Homininae into two tribes, one for Gorilla, and the other for Pan (chimpanzee) and Homo. Parsimony and maximum likelihood analyses, carried out on orthologous noncoding nucleotide sequences from primate beta-globin gene clusters, provide significant evidence for the human-chimpanzee tribe and overwhelming evidence for the human-chimpanzee-gorilla clade. These analyses also indicate that the rate of molecular evolution became slower in hominoids than in other primates and mammals.     

Maximum-likelihood models for combined analyses of multiple sequence data

Models of nucleotide substitution were constructed for combined analyses of heterogeneous sequence data (such as those of multiple genes) from the same set of species. The models account for different aspects of the heterogeneity in the evolutionary process of different genes, such as differences in nucleotide frequencies, in substitution rate bias (for example, the transition/transversion rate bias), and in the extent of rate variation across sites. Model parameters were estimated by maximum likelihood and the likelihood ratio test was used to test hypotheses concerning sequence evolution, such as rate constancy among lineages (the assumption of a molecular clock) and proportionality of branch lengths for different genes. The example data from a segment of the mitochondrial genome of six hominoid species (human, common and pygmy chimpanzees, gorilla, orangutan, and siamang) were analyzed. Nucleotides at the three codon positions in the protein-coding regions and from the tRNA-coding regions were considered heterogeneous data sets. Statistical tests showed that the amount of evolution in the sequence data reflected in the estimated branch lengths can be explained by the codon-position effect and lineage effect of substitution rates. The assumption of a molecular clock could not be rejected when the data were analyzed separately or when the rate variation among sites was ignored. However, significant differences in substitution rate among lineages were found when the data sets were combined and when the rate variation among sites was accounted for in the models. Under the assumption that the orangutan and African apes diverged 13 million years ago, the combined analysis of the sequence data estimated the times for the human-chimpanzee separation and for the separation of the gorilla as 4.3 and 6.8 million years ago, respectively.     

Evolutionary conservation of 5' upstream sequence of nine genes between human and great apes

Nucleotide sequences of nine 5' upstream gene regions for human, chimpanzee, gorilla, and orangutan were determined. We estimated nucleotide differences (d) for each region between human and great apes. The overall d was 0.027 (ranged from 0.004 to 0.052). Rates of nucleotide substitution were estimated by using d and divergence times of human, chimpanzee, gorilla, and orangutan. The overall rate of nucleotide substitution between human and other hominoids was estimated to be 0.52-0.85 x 10(-9). This rate in 5' upstream regions was lower than that of synonymous sites, suggesting that 5' upstream regions have evolved under some functional constraints. Because lower rates have been reported for coding sequences in primates compared to rodents, we also estimated the rate (1.17-1.76 x 10(-9)) of nucleotide substitutions for the corresponding 5' upstream regions in rodents (mouse/rat comparison). Thus the primate rate was lower than rodent rate also for the 5' upstream regions.     

Molecular phylogeny of the hominoids: inferences from multiple independent DNA sequence data sets

Consensus on the evolutionary relationships of humans, chimpanzees, and gorillas has not been reached, despite the existence of a number of DNA sequence data sets relating to the phylogeny, partly because not all gene trees from these data sets agree. However, given the well-known phenomenon of gene tree-species tree mismatch, agreement among gene trees is not expected. A majority of gene trees from available DNA sequence data support one hypothesis, but is this evidence sufficient for statistical confidence in the majority hypothesis? All available DNA sequence data sets showing phylogenetic resolution among the hominoids are grouped according to genetic linkage of their corresponding genes to form independent data sets. Of the 14 independent data sets defined in this way, 11 support a human- chimpanzee clade, 2 support a chimpanzee-gorilla clade, and one supports a human-gorilla clade. The hypothesis of a trichotomous speciation event leading to Homo; Pan, and Gorilla can be firmly rejected on the basis of this data set distribution. The multiple-locus test (Wu 1991), which evaluates hypotheses using gene tree-species tree mismatch probabilities in a likelihood ratio test, favors the phylogeny with a Homo-Pan clade and rejects the other alternatives with a P value of 0.002. When the probabilities are modified to reflect effective population size differences among different types of genetic loci, the observed data set distribution is even more likely under the Homo-Pan clade hypothesis. Maximum-likelihood estimates for the time between successive hominoid divergences are in the range of 300,000-2,800,000 years, based on a reasonable range of estimates for long-term hominoid effective population size and for generation time. The implication of the multiple-locus test is that existing DNA sequence data sets provide overwhelming and sufficient support for a human-chimpanzee clade: no additional DNA data sets need to be generated for the purpose of estimating hominoid phylogeny. Because DNA hybridization evidence (Caccone and Powell 1989) also supports a Homo-Pan clade, the problem of hominoid phylogeny can be confidently considered solved.      

人猿超科和旧大陆猴7种高等灵长类FKN全基因序列的测定及进化分析

测定人猿超科(人、黑猩猩、大猩猩、红毛猩猩和长臂猿)和旧大陆猴(猕猴和叶猴)7种高等灵长类FKN全基因序列, 探讨其系统进化分析。用简并引物PCR(Degenerated PCR)法分别扩增FKN的3个外显子, 其产物经琼脂糖凝胶回收、纯化后测序, 然后用BioEdit软件剪切拼接FKN基因全序列, 用DNAStar比对后比较基因和氨基酸序列同源性, Mega软件重构FKN基因进化树, 应用Datamonkey分析FKN的负选择位点。序列分析发现人猿超科较旧大陆猴FKN基因除了有散在的点突变外, 还有一明显的30 bp的核苷酸缺失突变; 人FKN基因序列与黑猩猩、大猩猩、红毛猩猩、长臂猿、猕猴和叶猴的同源性分别是99.2%、98.4%、98.1%、96.5%、95.9%和93.8%, 由此推导的氨基酸序列同源性分别是98.5%、98.0%、97.7%、94.7%、93.7%和90.5%; FKN基因进化树表明人与黑猩猩关系更近, FKN基因进化和通常认为的物种进化一致; Datamonkey分析结果显示FKN存在3个负选择位点53Q、84D、239N。成功获得人、黑猩猩、大猩猩、红毛猩猩、长臂猿、猕猴和叶猴7种高等灵长类物种FKN全基因序列, 为后续探讨FKN在高等灵长类物种进化过程中免疫学功能演变及其结构与功能的关系奠定基础。     

Slow molecular clocks in Old World monkeys,apes, and humans

Two longstanding issues on the molecular clock hypothesis are studied in this article. First, is there a global molecular clock in mammals? Although many authors have observed unequal rates of nucleotide substitution among mammalian lineages, some authors have proposed a global clock for all eutherians, i.e., a single global rate of 2.2 x 10(-9) substitutions per nucleotide site per year. We reexamine this issue using noncoding, nonrepetitive DNA from Old World monkeys (OWMs), chimpanzee, and human. First, using the minimal date of 6 MYA for the human-chimpanzee divergence and more than 2.5 million base pairs of genomic sequences from human and chimpanzee, we estimate a maximal rate of 0.99 x 10(-9) for noncoding, nonrepetitive genomic regions for these two species. This estimate is less than half of the proposed global rate and much smaller than the commonly used rate (3.5 x 10(-9)) for eutherians. Further, using a minimal date of 23 MYA for the human-OWM divergence, we estimate a maximal rate of 1.5 x 10(-9) for both introns and fourfold degenerate sites in humans and OWMs. In addition, with the New World monkey (NWM) lineage as an outgroup, we estimate that the rate of substitution in introns is 30% higher in the OWM lineage than in the human lineage. Clearly, there is no global molecular clock in eutherians. Second, although many studies have indicated considerable variation in the mutation rate among regions of the mammalian genome, a recent study proposed a uniform rate. Using new and existing intron sequence data from higher primates, we find significant rate variation among genomic regions and a positive correlation between the rate of substitution and the GC content, refuting the claim of a uniform rate.     

Molecular phylogeny and evolution of the human endogenous retrovirus HERV-W LTR family in hominoid primates

Long terminal repeats (LTRs) of human endogenous retrovirus (HERV) have contributed to the structural change or genetic variation of primate genome that are connected to speciation and evolution. Using genomic DNAs that were derived from hominoid primates (chimpanzee, gorilla, orangutan, and gibbon), we performed PCR amplification and identified thirty HERV-W LTR elements. These LTR elements showed a 82-98% sequence similarity with HERV-W LTR (AF072500). Specifically, additional sequences (GCCACCACCACTGTTT in the gorilla and TGCTGCTGACTCCCATCC in the gibbon) were noticed. Clone OR3 from the orangutan and clone GI2 from the gibbon showed a 100% sequence similarity, although they are different species. This indicates that both LTR elements were proliferated during the last 2 to 5 million years from the integration of the original LTR element. A phylogenetic tree that was obtained by the neighbor-joining method revealed a wide overlap of the LTR elements across species, suggesting that the HERV-W LTR family evolved independently during the hominoid evolution.     

Identification of human specific gene duplications relative to other primates by array CGH and quantitative PCR

In order to identify human lineage specific (HLS) copy number differences (CNDs) compared to other primates, we performed pair wise comparisons (human vs. chimpanzee, gorilla and orangutan) by using cDNA array comparative genomic hybridization (CGH). A set of 23 genes with HLS duplications were identified, as well as other lineage differences in gene copy number specific of chimpanzee, gorilla and orangutan. Each species has gained more copies of specific genes rather than losing gene copies. Eleven of the 23 genes have only been observed to have undergone HLS duplication in Fortna et al. (2004) and in the present study. Then, seven of these 11 genes were analyzed by quantitative PCR in chimpanzee, gorilla and orangutan, as well as in other six primate species (Hylobates lar, Cercopithecus aethiops, Papio hamadryas, Macaca mulatta, Lagothrix lagothricha, and Saimiri sciureus). Six genes confirmed array CGH data, and four of them appeared to have bona fide HLS duplications (ABCB10, E2F6, CDH12, and TDG genes). We propose that these gene duplications have a potential to contribute to specific human phenotypes.     

A young Alu subfamily amplified independently in human and African great apes lineages.

A variety of Alu subfamilies amplified in primate genomes at different evolutionary time periods. Alu Sb2 belongs to a group of young subfamilies with a characteristic two-nucleotide deletion at positions 65/66. It consists of repeats having a 7-nucleotide duplication of a sequence segment involving positions 246 through 252. The presence of Sb2 inserts was examined in five genomic loci in 120 human DNA samples as well as in DNAs of higher primates. The lack of the insertional polymorphism seen at four human loci and the absence of orthologous inserts in apes indicated that the examined repeats retroposed early in the human lineage, but following the divergence of great apes. On the other hand, similar analysis of the fifth locus (butyrylcholinesterase gene) suggested contemporary retropositional activity of this subfamily. By a semi-quantitative PCR, using a primer pair specific for Sb2 repeats, we estimated their copy number at about 1500 per human haploid genome; the corresponding numbers in chimpanzee and gorilla were two orders of magnitude lower, while in orangutan and gibbon the presence of Sb2 Alu was hardly detectable. Sequence analysis of PCR-amplified Sb2 repeats from human and African great apes is consistent with the model in which the founding of Sb2 subfamily variants occurred independently in chimpanzee, gorilla and human lineages.     

Paleo-demography of the Drosophila melanogaster subgroup: application of the maximum likelihood method

The species divergence times and demographic histories of Drosophila melanogaster and its three sibling species, D. mauritiana, D. simulans, and D. yakuba, were investigated using a maximum likelihood (ML) method. Thirty-nine orthologous loci for these four species were retrieved from DDBJ/EMBL/GenBank database. Both autosomal and X-linked loci were used in this study. A significant degree of rate heterogeneity across loci was observed for each pair of species. Most loci have the GC content greater than 50% at the third codon position. The codon usage bias in Drosophila loci is considered to result in the high GC content and the heterogenous rates across loci. The chi-square, G, and Fisher's exact tests indicated that data sets with 11, 23, and 9 pairs of DNA sequences for the comparison of D. melanogaster with D. mauritiana, D. simulans, and D. yakuba, respectively, retain homogeneous rates across loci. We applied the ML method to these data sets to estimate the DNA sequence divergences before and after speciation of each species pair along with their standard deviations. Using 1.6 x 10(-8) as the rate of nucleotide substitutions per silent site per year, our results indicate that the D. melanogaster lineage split from D. yakuba approximately 5.1 +/- 0.8 million years ago (mya), D. mauritiana 2.7 +/- 0.4 mya, and D. simulans 2.3 +/- 0.3 mya. It implies that D. melanogaster became distinct from D. mauritiana and D. simulans at approximately the same time and from D. yakuba no earlier than 10 mya. The effective ancestral population size of D. melanogaster appears to be stable over evolutionary time. Assuming 10 generations per year for Drosophila, the effective population size in the ancestral lineage immediately prior to the time of species divergence is approximately 3 x 10(6), which is close to that estimated for the extant D. melanogaster population. The D. melanogaster did not encounter any obvious bottleneck during the past 10 million years.     

Methods for computing the standard errors of branching points in an evolutionary tree and their application to molecular data from humans and apes

Statistical methods for computing the standard errors of the branching points of an evolutionary tree are developed. These methods are for the unweighted pair-group method-determined (UPGMA) trees reconstructed from molecular data such as amino acid sequences, nucleotide sequences, restriction-sites data, and electrophoretic distances. They were applied to data for the human, chimpanzee, gorilla, orangutan, and gibbon species. Among the four different sets of data used, DNA sequences for an 895-nucleotide segment of mitochondrial DNA (Brown et al. 1982) gave the most reliable tree, whereas electrophoretic data (Bruce and Ayala 1979) gave the least reliable one. The DNA sequence data suggested that the chimpanzee is the closest and that the gorilla is the next closest to the human species. The orangutan and gibbon are more distantly related to man than is the gorilla. This topology of the tree is in agreement with that for the tree obtained from chromosomal studies and DNA-hybridization experiments. However, the difference between the branching point for the human and the chimpanzee species and that for the gorilla species and the human-chimpanzee group is not statistically significant. In addition to this analysis, various factors that affect the accuracy of an estimated tree are discussed.