Time of the deepest root for polymorphism in human mitochondrial DNA

Summary A molecular clock analysis was carried out on the nucleotide sequences of parts of the major noncoding region of mitochondrial DNA (mtDNA) from the major geographic populations of humans. Dates of branchings in the mtDNA tree among humans were estimated with an improved maximum likelihood method. Two species of chimpanzees were used as an outgroup, and the mtDNA clock was calibrated by assuming that the chimpanzee/human split occurred 4 million years ago, following our earlier works. A model of homogeneous evolution among sites does not fit well with the data even within hypervariable segments, and hence an additional parameter that represents a proportion of variable sites was introduced. Taking account of this heterogeneity among sites, the date for the deepest root of the mtDNA tree among humans was estimated to be 280,000±50,000 years old (±1 SE), although there remains uncertainty about the constancy of the evolutionary rate among lineages. The evolutionary rate of the most rapidly evolving sites in mtDNA was estimated to be more than 100 times greater than that of a nuclear pseudogene.     

Evaluation of the maximum likelihood estimate of the evolutionary tree topologies from DNA sequence data,and the branching order in hominoidea

Summary A maximum likelihood method for inferring evolutionary trees from DNA sequence data was developed by Felsenstein (1981). In evaluating the extent to which the maximum likelihood tree is a significantly better representation of the true tree, it is important to estimate the variance of the difference between log likelihood of different tree topologies. Bootstrap resampling can be used for this purpose (Hasegawa et al. 1988; Hasegawa and Kishino 1989), but it imposes a great computation burden. To overcome this difficulty, we developed a new method for estimating the variance by expressing it explicitly.The method was applied to DNA sequence data from primates in order to evaluate the maximum likelihood branching order among Hominoidea. It was shown that, although the orangutan is convincingly placed as an outgroup of a human and African apes clade, the branching order among human, chimpanzee, and gorilla cannot be determined confidently from the DNA sequence data presently available when the evolutionary rate constancy is not assumed.     

Probability distribution of molecular evolutionary trees: A new method of phylogenetic inference

A new method is presented for inferring evolutionary trees using nucleotide sequence data. The birth-death process is used as a model of speciation and extinction to specify the prior distribution of phylogenies and branching times. Nucleotide substitution is modeled by a continuous-time Markov process. Parameters of the branching model and the substitution model are estimated by maximum likelihood. The posterior probabilities of different phylogenies are calculated and the phylogeny with the highest posterior probability is chosen as the best estimate of the evolutionary relationship among species. We refer to this as the maximum posterior probability (MAP) tree. The posterior probability provides a natural measure of the reliability of the estimated phylogeny. Two example data sets are analyzed to infer the phylogenetic relationship of human, chimpanzee, gorilla, and orangutan. The best trees estimated by the new method are the same as those from the maximum likelihood analysis of separate topologies, but the posterior probabilities are quite different from the bootstrap proportions. The results of the method are found to be insensitive to changes in the rate parameter of the branching process. Correspondence to: Z. Yang     

Mitochondrial DNA sequences of primates: Tempo and mode of evolution

Summary We cloned and sequenced a segment of mitochondrial DNA from human, chimpanzee, gorilla, orangutan, and gibbon. This segment is 896 bp in length, contains the genes for three transfer RNAs and parts of two proteins, and is homologous in all 5 primates. The 5 sequences differ from one another by base substitutions at 283 positions and by a deletion of one base pair. The sequence differences range from 9 to 19% among species, in agreement with estimates from cleavage map comparisons, thus confirming that the rate of mtDNA evolution in primates is 5 to 10 times higher than in nuclear DNA. The most striking new finding to emerge from these comparisons is that transitions greatly outnumber transversions. Ninety-two percent of the differences among the most closely related species (human, chimpanzee, and gorilla) are transitions. For pairs of species with longer divergence times, the observed percentage of transitions falls until, in the case of comparisons between primates and non-primates, it reaches a value of 45. The time dependence is probably due to obliteration of the record of transitions by multiple substitutions at the same nucleotide site. This finding illustrates the importance of choosing closely related species for analysis of the evolutionary process. The remarkable bias toward transitions in mtDNA evolution necessitates the revision of equations that correct for multiple substitutions at the same site. With revised equations, we calculated the incidence of silent and replacement substitutions in the two protein-coding genes. The silent substitution rate is 4 to 6 times higher than the replacement rate, indicating strong functional constraints at replacement sites. Moreover, the silent rate for these two genes is about 10% per million years, a value 10 times higher than the silent rate for the nuclear genes studied so far. In addition, the mean substitution rate in the three mitochondrial tRNA genes is at least 100 times higher than in nuclear tRNA genes. Finally, genealogical analysis of the sequence differences supports the view that the human lineage branched off only slightly before the gorilla and chimpanzee lineages diverged and strengthens the hypothesis that humans are more related to gorillas and chimpanzees than is the orangutan.Abbreviations mtDNA mitochondrial DNA - bp base pair - URF unidentified reading frame     

Forty Million Years of Independent Evolution: A Mitochondrial Gene and Its Corresponding Nuclear Pseudogene in Primates

Sequences from nuclear mitochondrial pseudogenes (numts) that originated by transfer of genetic information from mitochondria to the nucleus offer a unique opportunity to compare different regimes of molecular evolution. Analyzing a 1621-nt-long numt of the rRNA specifying mitochondrial DNA residing on human chromosome 3 and its corresponding mitochondrial gene in 18 anthropoid primates, we were able to retrace about 40 MY of primate rDNA evolutionary history. The results illustrate strengths and weaknesses of mtDNA data sets in reconstructing and dating the phylogenetic history of primates. We were able to show the following. In contrast to numt-DNA, (1) the nucleotide composition of mtDNA changed dramatically in the different primate lineages. This is assumed to lead to significant misinterpretations of the mitochondrial evolutionary history. (2) Due to the nucleotide compositional plasticity of primate mtDNA, the phylogenetic reconstruction combining mitochondrial and nuclear sequences is unlikely to yield reliable information for either tree topologies or branch lengths. This is because a major part of the underlying sequence evolution model — the nucleotide composition — is undergoing dramatic change in different mitochondrial lineages. We propose that this problem is also expressed in the occasional unexpected long branches leading to the “common ancestor” of orthologous numt sequences of different primate taxa. (3) The heterogeneous and lineage-specific evolution of mitochondrial sequences in primates renders molecular dating based on primate mtDNA problematic, whereas the numt sequences provide a much more reliable base for dating.[Reviewing Editor: Dr. Rafael Zardoya]     

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.      

Two independent mutational events in the loss of urate oxidase during hominoid evolution

Summary Urate oxidase was lost in hominoids during primate evolution. The mechanism and biological reason for this loss remain unknown. In an attempt to address these questions, we analyzed the sequence of urate oxidase genes from four species of hominoids: human (Homo sapiens), chimpanzee (Pan troglodytes), orangutan (Pongo pygmaeus), and gibbon (Hylobates). Two nonsense mutations at codon positions 33 and 187 and an aberrant splice site were found in the human gene. These three deleterious mutations were also identified in the chimpanzee. The nonsense mutation at codon 33 was observed in the orangutan urate oxidase gene. None of the three mutations was present in the gibbon; in contrast, a 13-bp deletion was identified that disrupted the gibbon urate oxidase reading frame. These results suggest that the loss of urate oxidase during the evolution of hominoids could be caused by two independent events after the divergence of the gibbon lineage; the nonsense mutation at codon position 33 resulted in the loss of urate oxidase activity in the human, chimpanzee, and orangutan, whereas the 13-bp deletion was responsible for the urate oxidase deficiency in the gibbon. Because the disruption of a functional gene by independent events in two different evolutionary lineages is unlikely to occur on a chance basis, our data favor the hypothesis that the loss of urate oxidase may have evolutionary advantages. Offprint requests to: C.T. Caskey     

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     

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.     

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.     

Karyotype of the gibbons hylobates lar and h. moloch inversion in chromosome 7.

A karyotype of the gibbon, Hylobates, has been prepared based on the chromosome banding patterns produced by quinacrine, trypsin-Giemsa, and centromeric heterochromatin stains. The banding patterns of H. lar and H. moloch are virtually identical. No brilliant quinacrine-fluorescent areas are present. The banding pattern of most of the gibbon chromosomes show less resemblance to those of the human, chimpanzee, gorilla, or orangutan than the chromosomes of the higher primates do to each other, suggesting a relatively large evolutionary separation of the gibbon from the higher primates. A pericentric inversion of chromosome 7 is present in one gibbon.     

Comparative immunochemical studies of primate hemoglobins

The antigenic properties of a number of chromatographically purified primate hemoglobins were compared to those of normal human hemoglobin using a sensitive radioimmunochemical procedure. The degree of inhibition of the antigen-antibody reaction with heterologous hemoglobins appeared to be related to the structural similarity of these proteins to the normal human hemoglobin immunogen. With the exception of the baboon hemoglobin, the antigenicity of the hemoglobins paralleled the phylogeny of the primates. The gorilla and chimpanzee hemoglobins were antigenically identical to normal human hemoglobin, whereas the gibbon and orangutan hemoglobins were substantially more variable. Of the Old World monkey hemoglobins examined, the baboon produced lower inhibition values, suggesting a greater degree of structural dissimilarity than other Cercopithecoidea hemoglobins, which is compatible with a greater rate of evolutionary change occurring in this protein. Using the known amino acid sequences of human and other primate hemoglobins, we have attempted to identify antigenic determinant areas of the proteins.     

Evolution of glycophorin A in the hominoid primates studied with monoclonal antibodies, and description of a sialoglycoprotein analogous to human glycophorin B in chimpanzee

Comparison of human and primate erythrocyte membrane sialoglycoproteins showed that common chimpanzee, dwarf chimpanzee, gorilla, orangutan, and gibbon have major periodic acid Schiff-positive proteins resembling human glycophorin A (GPA) monomer and dimer in electrophoretic mobility on sodium dodecyl sulfate-polyacrylamide gels. Immunoperoxidase staining of Western blots with monoclonal antibodies to human GPA showed that these primate bands express some GPA antigenic determinants. A new sialoglycoprotein analogous to human glycophorin B (GPB) was detected in common chimpanzee. Although human MN blood group phenotype results from an amino acid polymorphism of GPA, Western blots showed that in chimpanzee sialoglycoprotein (GPAch) always expresses the M blood group, whereas chimpanzee sialoglycoprotein (GPBch) expresses either the N blood group or a null phenotype. This result explains the detection of M and MN, but not of N, blood group phenotypes in chimpanzee. GPBch has higher apparent m.w. than human GPB, is present in the erythrocyte membrane in greater quantity than human GPB, and contains trypsin cleavage site(s) and the 10F7 determinant (both found on human GPA but not GPB). Expression of human GPA antigenic determinants was consistent with the phylogeny of the hominoid primates; common and dwarf chimpanzee expressed most of the determinants tested, gorilla and orangutan an intermediate number, and gibbon and siamang the least. Of the GPA antigenic determinants examined, the MN blood group determinants were most consistently expressed during evolution of the hominoid primates. The results suggested that variability in expression of GPA antigenic determinants between species was due to both differences in amino acid sequence and glycosylation.     

Silencing, positive selection and parallel evolution: busy history of primate cytochromes C

Cytochrome c (cyt c) participates in two crucial cellular processes, energy production and apoptosis, and unsurprisingly is a highly conserved protein. However, previous studies have reported for the primate lineage (i) loss of the paralogous testis isoform, (ii) an acceleration and then a deceleration of the amino acid replacement rate of the cyt c somatic isoform, and (iii) atypical biochemical behavior of human cyt c. To gain insight into the cause of these major evolutionary events, we have retraced the history of cyt c loci among primates. For testis cyt c, all primate sequences examined carry the same nonsense mutation, which suggests that silencing occurred before the primates diversified. For somatic cyt c, maximum parsimony, maximum likelihood, and Bayesian phylogenetic analyses yielded the same tree topology. The evolutionary analyses show that a fast accumulation of non-synonymous mutations (suggesting positive selection) occurred specifically on the anthropoid lineage root and then continued in parallel on the early catarrhini and platyrrhini stems. Analysis of evolutionary changes using the 3D structure suggests they are focused on the respiratory chain rather than on apoptosis or other cyt c functions. In agreement with previous biochemical studies, our results suggest that silencing of the cyt c testis isoform could be linked with the decrease of primate reproduction rate. Finally, the evolution of cyt c in the two sister anthropoid groups leads us to propose that somatic cyt c evolution may be related both to COX evolution and to the convergent brain and body mass enlargement in these two anthropoid clades.     

Primate evolution at the DNA level and a classification of hominoids

Summary The genetic distances among primate lineages estimated from orthologous noncoding nucleotide sequences of -type globin loci and their flanking and intergenic DNA agree closely with the distances (delta T₅₀H values) estimated by cross hybridization of total genomic single-copy DNAs. These DNA distances and the maximum parsimony tree constructed for the nucleotide sequence orthologues depict a branching pattern of primate lineages that is essentially congruent with the picture from phylogenetic analyses of morphological characters. The molecular evidence, however, resolves ambiguities in the morphological picture and provides an objective view of the cladistic position of humans among the primates. The molecular data group humans with chimpanzees in subtribe Hominina, with gorillas in tribe Hominini, orangutans in subfamily Homininae, gibbons in family Hominidae, Old World monkeys in infraorder Catarrhini, New World monkeys in semisuborder Anthropoidea, tarsiers in suborder Haplorhini, and strepsirhines (lemuriforms and lorisiforms) in order Primates. A seeming incongruency between organismal and molecular levels of evolution, namely that morphological evolution appears to have speeded up in higher primates, especially in the lineage to humans, while molecular evolution has slowed down, may have the trivial explanation that relatively small genetic changes may sometimes result in marked phenotypic changes.     

Evolution of DNA on Y-chromosome in hominoid primates as examined by PCR

Every species of non-human primates, especially those of hominoids, has a variety of reproductive structures and accompanying male traits, such as sexual dimorphism and relative size of testis to body weight, which may be at least partly triggered by DNA on the Y-chromosome. Recently, a panel of PCR (Polymerase Chain Reaction) primer sets were designed to amplify various DNA segments spread over the human Y-chromosome. We applied these primer sets for amplification of DNA segments on the Y-chromosome of hominoid species: chimpanzee, bonobo (Pygmy chimpanzee), gorilla, orangutan, whitehanded gibbon, agile gibbon, and Japanese monkey as an out group. The DNA segments including SRY, testis determining factor, and ZFX/ZFY could be amplified clearly in males of all species examined. These highly conserved genes may serve important biological functions. However, as the phylogenic distance from humans increased, some of the DNA segments could not be amplified. For example, DYZ1 (SY160) could be amplified only using human DNA as a template, and DYF60S1 (SY61), DYZ217 (SY126) and DYS233 (SY148) could be amplified only using human and African great ape DNA. It is interesting to note that locus DYS250 (SY17) could not be amplified in chimpanzee and bonobo but amplified in gorilla and orangutan. Locus DYS251 (SY18) was amplified in all species except the white-handed gibbon. These results indicate that a variety of evolutionary events including mutation, deletion, insertion, and rearrangement occurred in Y-chromosome DNA during primate evolution.     

灵长类比较基因组学的研究进展

随着人类和黑猩猩全基因组测序工作宣布完成,以及其他灵长类基因组测序工作的逐步开展,目前已经积累了大量的灵长类基因组数据,一个崭新的研究领域——灵长类比较基因组学应运而生。该文主要通过对人类和其他非人灵长类系统关系和基因组结构的比较,从系统进化、基因组结构和基因表达调控等方面评述该领域的研究进展,阐述人类、黑猩猩与其他非人灵长类之间的主要生物学差异,揭示人类进化的生物学机制。     

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.     

A molecular phylogeny of living primates

Comparative genomic analyses of primates offer considerable potential to define and understand the processes that mold, shape, and transform the human genome. However, primate taxonomy is both complex and controversial, with marginal unifying consensus of the evolutionary hierarchy of extant primate species. Here we provide new genomic sequence (~8 Mb) from 186 primates representing 61 (~90%) of the described genera, and we include outgroup species from Dermoptera, Scandentia, and Lagomorpha. The resultant phylogeny is exceptionally robust and illuminates events in primate evolution from ancient to recent, clarifying numerous taxonomic controversies and providing new data on human evolution. Ongoing speciation, reticulate evolution, ancient relic lineages, unequal rates of evolution, and disparate distributions of insertions/deletions among the reconstructed primate lineages are uncovered. Our resolution of the primate phylogeny provides an essential evolutionary framework with far-reaching applications including: human selection and adaptation, global emergence of zoonotic diseases, mammalian comparative genomics, primate taxonomy, and conservation of endangered species.