Genomic relationships and speciation times of human, chimpanzee, and gorilla inferred from a coalescent hidden Markov model

The genealogical relationship of human, chimpanzee, and gorilla varies along the genome. We develop a hidden Markov model (HMM) that incorporates this variation and relate the model parameters to population genetics quantities such as speciation times and ancestral population sizes. Our HMM is an analytically tractable approximation to the coalescent process with recombination, and in simulations we see no apparent bias in the HMM estimates. We apply the HMM to four autosomal contiguous human–chimp–gorilla–orangutan alignments comprising a total of 1.9 million base pairs. We find a very recent speciation time of human–chimp (4.1 ± 0.4 million years), and fairly large ancestral effective population sizes (65,000 ± 30,000 for the human–chimp ancestor and 45,000 ± 10,000 for the human–chimp–gorilla ancestor). Furthermore, around 50% of the human genome coalesces with chimpanzee after speciation with gorilla. We also consider 250,000 base pairs of X-chromosome alignments and find an effective population size much smaller than 75% of the autosomal effective population sizes. Finally, we find that the rate of transitions between different genealogies correlates well with the region-wide present-day human recombination rate, but does not correlate with the fine-scale recombination rates and recombination hot spots, suggesting that the latter are evolutionarily transient.     

Evolution of X-Degenerate Y Chromosome Genes in Greater Apes: Conservation of Gene Content in Human and Gorilla,But Not Chimpanzee

Compared with the X chromosome, the mammalian Y chromosome is considerably diminished in size and has lost most of its ancestral genes during evolution. Interestingly, for the X-degenerate region on the Y chromosome, human has retained all 16 genes, while chimpanzee has lost 4 of the 16 genes since the divergence of the two species. To uncover the evolutionary forces governing ape Y chromosome degeneration, we determined the complete sequences of the coding exons and splice sites for 16 gorilla Y chromosome genes of the X-degenerate region. We discovered that all studied reading frames and splice sites were intact, and thus, this genomic region experienced no gene loss in the gorilla lineage. Higher nucleotide divergence was observed in the chimpanzee than the human lineage, particularly for genes with disruptive mutations, suggesting a lack of functional constraints for these genes in chimpanzee. Surprisingly, our results indicate that the human and gorilla orthologues of the genes disrupted in chimpanzee evolve under relaxed functional constraints and might not be essential. Taking mating patterns and effective population sizes of ape species into account, we conclude that genetic hitchhiking associated with positive selection due to sperm competition might explain the rapid decline in the Y chromosome gene number in chimpanzee. As we found no evidence of positive selection acting on the X-degenerate genes, such selection likely targets other genes on the chimpanzee Y chromosome. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Examination of hominid evolution by DNA sequence homology

Hominoid phylogeny was investigated in terms of unique DNA sequence homologies. In comparisons from the human standpoint the ΔTe₅₀ DNA values were Man 0, chimpanzee 0·7, gorilla 1·4, gibbon 2·7, orangutan 2·9, and African green monkey 5·7. In comparisons from the orangutan standpoint the ΔTe₅₀ DNA values were orangutan 0, chimpanzee 1·8, Man 1·9, gorilla 2·3, gibbon 2·4 and African green monkey 4·3. These results indicate that chimpanzee and gorilla are cladistically closer to Man than to orangutan and other primates, and that gorilla DNA may have diverged slightly more from the ancestral state than chimpanzee or human DNA. Comparisons from chimpanzee and gorilla DNA standpoints are needed to achieve a more definitive picture of hominoid phylogeny.     

人与大猩猩,黑猩猩和猩猩亲缘关系的探讨

有关人锆超科的系统发育仍然存在刍议。争论焦点在与大猩猩和黑猩猩哪 个关系更近一点。酪氨酸酶是黑色素合成中的关键酶,酪氨酶基因的突变将导致白化病。测定了人猿科中大猩猩,黑猩猩、猩猩和长臂锆产基因全部５个外显子的ＤＮＡ序列。     

Primate genes for glycophorins carrying MN blood group antigens

Glycophorin A, B, and E genes were derived from a common ancestral gene and this gene family appeared during primate evolution, probably between orangutan and gorilla divergences. Based on the study of genomic structures of these human glycophorins and the genetic and immunological study of primate glycophorins, we hypothesize that chimpanzee and gorilla glycophorin B could possess a longer extracellular region and carry a stronger N blood group antigenicity compared with that of the human.     

Concerted evolution of the primate immunoglobulin alpha-gene through gene conversion.

We determined four nucleotide sequences of the hominoid immunoglobulin alpha (C alpha) genes (chimpanzee C alpha 2, gorilla C alpha 2, and gibbon C alpha 1 and C alpha 2 genes), which made possible the examination of gene conversions in all hominoid C alpha genes. The following three methods were used to detect gene conversions: 1) phenetic tree construction; 2) detection of a DNA segment with extremely low variability between duplicated C alpha genes; and 3) a site by site search of shared nucleotide changes between duplicated C alpha genes. Results obtained from method 1 indicated a concerted evolution of the duplicated C alpha genes in the human, chimpanzee, gorilla, and gibbon lineages, while results obtained from method 2 suggested gene conversions in the human, gorilla, and gibbon C alpha genes. With method 3 we identified clusters of shared nucleotide changes between duplicated C alpha genes in human, chimpanzee, gorilla, and gibbon lineages, and in their hypothetical ancestors. In the present study converted regions were identified over the entire C alpha gene region excluding a few sites in the coding region which have escaped from gene conversion. This indicates that gene conversion is a general phenomenon in evolution, that can be clearly observed in non-functional regions.     

Synteny comparison between apes and human using fine-mapping of the genome

Comparing the genomes of the great apes and human should provide novel information concerning the origins of humankind. Relative to the great apes, the human karyotype has one fewer chromosome pair, as human chromosome 2 derived from the telomeric fusion of two ancestral primate chromosomes. To identify the genomic rearrangements that accompanied human speciation, we initiated a comparative study between human, chimpanzee, and gorilla. Using the HAPPY mapping method, an acellular adaptation of the radiation hybrid method, we mapped a few hundred markers on the human, chimpanzee, and gorilla genomes. This allowed us to identify several chromosome rearrangements, in particular a pericentric inversion and a translocation. We precisely localized the synteny breakpoint that led to the formation of human chromosome 2. This breakpoint was confirmed by FISH mapping.     

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.     

Characterization of a Novel Class of Interspersed LTR Elements in Primate Genomes: Structure, Genomic Distribution, and Evolution

Retrovirus-like sequences and their solitary (solo) long terminal repeats (LTRs) are common repetitive elements in eukaryotic genomes. We reported previously that the tandemly arrayed genes encoding U2 snRNA (the RNU2 locus) in humans and apes contain a solo LTR (U2-LTR) which was presumably generated by homologous recombination between the two LTRs of an ancestral provirus that is retained in the orthologous baboon RNU2 locus. We have now sequenced the orthologous U2-LTRs in human, chimpanzee, gorilla, orangutan, and baboon and examined numerous homologs of the U2-LTR that are dispersed throughout the human genome. Although these U2-LTR homologs have been collectively referred to as LTR13 in the literature, they do not display sequence similarity to any known retroviral LTRs; however, the structure of LTR13 closely resembles that of other retroviral LTRs with a putative promoter, polyadenylation signal, and a tandemly repeated 53-bp enhancer-like element. Genomic blotting indicates that LTR13 is primate-specific; based on sequence analysis, we estimate there are about 2,500 LTR13 elements in the human genome. Comparison of the primate U2-LTR sequences suggests that the homologous recombination event that gave rise to the solo U2-LTR occurred soon after insertion of the ancestral provirus into the ancestral U2 tandem array. Phylogenetic analysis of the LTR13 family confirms that it is diverse, but the orthologous U2-LTRs form a coherent group in which chimpanzee is closest to the humans; orangutan is a clear outgroup of human, chimpanzee, and gorilla; and baboon is a distant relative of human, chimpanzee, gorilla, and orangutan. We compare the LTR13 family with other known LTRs and consider whether these LTRs might play a role in concerted evolution of the primate RNU2 locus. Received: 29 September 1997 / Accepted: 16 January 1998     

A Y-chromosomal DNA fragment is conserved in human and chimpanzee.

A human male-specific Y-chromosomal DNA fragment (lambda YH2D6) has been isolated. By deletion-mapping analysis, 2D6 has been localized to the euchromatic portion of the long arm (Yq11) of the human Y chromosome. Among great apes, this fragment was found to be conserved in male chimpanzee but was lacking in male gorilla and male orangutan. No homologous fragments were detected in females of orangutan, gorilla, chimpanzee, or human. Nucleotide sequence analysis indicated the presence of partial-Alu-elements and of sequences similar to the GATA repeats of the snake Bkm sequence.     

Fluorescent (F) bodies in the spermatozoa of man and the great apes

Mature spermatozoa of the chimpanzee (Pan troglodytes), the gorilla (Gorilla gorilla), and the orangutan (Pongo pygmaeus) were stained with quinacrine dihydrochloride. Fluorescent (F) bodies were visualized in the spermatozoa of the chimpanzee and gorilla but were absent in the orangutan, in which there is no brilliant fluorescence in any chromosome. The F bodies appeared to be randomly located in the sperm heads of these two species, as they usually are in human spermatozoa. However, the proportion of sperm showing one or more F bodies in the chimpanzee and gorilla was not comparable to what is usually found in man. The F bodies in the chimpanzee presumably represent brilliant regions in the autosomes, since the Y chromosome has no brilliant fluorescence in this species. This is contrary to man, in which the F body is an useful indicator of the Y chromosome. In the gorilla, the F bodies probably correspond to both the Y chromosome and to some brilliant regions in the autosomes.     

The sequence of the gorilla fetal globin genes: evidence for multiple gene conversions in human evolution

Two fetal globin genes (G gamma and A gamma) from one chromosome of a lowland gorilla (Gorilla gorilla gorilla) have been sequenced and compared to three human loci (a G gamma-gene and two A gamma-alleles). A comparison of regions of local homology among these five sequences indicates that long after the duplication that produced the two nonallelic gamma-globin loci of catarrhine primates, about 35 million years (Myr) ago, at least one gene conversion event occurred between these loci. This conversion occurred not long before the ancestral divergence (about 6 Myr ago) of Homo and Gorilla. After this ancestral divergence, a minimum of three more gene conversion events occurred in the human lineage. Each human A gamma-allele shares specific sequence features with the gorilla A gamma-gene; one such distinctive allelic feature involves the simple repeated sequence in IVS 2. This suggests that early in the human lineage the A gamma-genes may have undergone a crossing-over event mediated by this simple repeated sequence. The DNA sequences from coding regions of both G gamma- and A gamma-loci, a comparison of 292 codons in the corresponding gorilla and human genes, show an unusually low evolutionary rate, with only two nonsilent differences and, surprisingly, not even one silent substitution. The two nonsynonymous substitutions observed predict a glycine at codon 73 and an arginine at codon 104 in the gorilla A gamma-sequence rather than aspartic acid and lysine, respectively, in human A gamma. Because only arginine has been found at position 104 in gamma-chains of Old World monkeys, it may represent the ancestral residue lost in gorilla and human G gamma-chains and in the human A gamma-chain. Possibly the arginine codon (AGG) was replaced by the lysine codon (AAG) in the G gamma-gene of a common ancestor of Homo and Gorilla and then was transferred to the A gamma-gene by subsequent conversions in the human lineage. DNA sequence conversions, similar to that attributed to the fetal gamma-globin genes, appear to be relatively frequent phenomena and, if widespread throughout the genome, may have profound evolutionary consequences.      

Evolution of Genomic Structures on Mammalian Sex Chromosomes

Throughout mammalian evolution, recombination between the two sex chromosomes was suppressed in a stepwise manner. It is thought that the suppression of recombination led to an accumulation of deleterious mutations and frequent genomic rearrangements on the Y chromosome. In this article, we review three evolutionary aspects related to genomic rearrangements and structures, such as inverted repeats (IRs) and palindromes (PDs), on the mammalian sex chromosomes. First, we describe the stepwise manner in which recombination between the X and Y chromosomes was suppressed in placental mammals and discuss a genomic rearrangement that might have led to the formation of present pseudoautosomal boundaries (PAB). Second, we describe ectopic gene conversion between the X and Y chromosomes, and propose possible molecular causes. Third, we focus on the evolutionary mode and timing of PD formation on the X and Y chromosomes. The sequence of the chimpanzee Y chromosome was recently published by two groups. Both groups suggest that rapid evolution of genomic structure occurred on the Y chromosome. Our re-analysis of the sequences confirmed the species-specific mode of human and chimpanzee Y chromosomal evolution. Finally, we present a general outlook regarding the rapid evolution of mammalian sex chromosomes.     

A genomic region evolving toward different GC contents in humans and chimpanzees indicates a recent and regionally limited shift in the mutation pattern

DNA sequences evolving differently in the human and chimpanzee genomes signal recent and regionally limited changes in the process of DNA sequence evolution. Here we present the comparison of 90 kb from the nonrecombining part of the human Y chromosome to the corresponding part of the chimpanzee genome using gorilla as out-group. Our results reveal a significant difference in the region-specific substitution process among the human and chimpanzee lineages. As a consequence, this region experiences a change in its GC content on the human lineage while it resides in compositional equilibrium on the chimpanzee lineage. Based on our analysis, we suggest a recent and species-specific shift in the region's mutation pattern as the cause of its differing evolution in humans and chimpanzees.     

Meiotic Gene Conversion Tract Length Distribution within the Rosy Locus of Drosophila Melanogaster

Employing extensive co-conversion data for selected and unselected sites of known molecular location in the rosy locus of Drosophila melanogaster, we determine the parameters of meiotic gene conversion tract length distribution. The tract length distribution for gene conversion events can be approximated by the equation P(L >/= n) = (n) where P is the probability that tract length (L) is greater than or equal to a specified number of nucleotides (n). From the co-conversion data, a maximum likelihood estimate with standard error for is 0.99717 +/- 0.00026, corresponding to a mean conversion tract length of 352 base pairs. (Thus, gene conversion tract lengths are sufficiently small to allow for extensive shuffling of DNA sequence polymorphisms within a gene.) For selected site conversions there is a bias towards recovery of longer tracts. The distribution of conversion tract lengths associated with selected sites can be approximated by the equation P(L >/= n| selected = (n)(1 - n + n/), where P is now the probability that a selected site tract length (L) is greater than or equal to a specified number of nucleotides (n). For the optimal value of determined from the co-conversion analysis, the mean conversion tract length for selected sites is 706 base pairs. We discuss, in the light of this and other studies, the relationship between meiotic gene conversion and P element excision induced gap repair and determine that they are distinct processes defined by different parameters and, possibly, mechanisms.     

Nucleotide sequence of the beta-globin genes in gorilla and macaque: The origin of nucleotide polymorphisms in human

Summary Part of the beta-globin genes ofMacaca cynomolgus andGorilla gorilla has been cloned and sequenced. Ten putatively neutral nucleotide polymorphisms have been described at the beta-globin locus in humans. They are associated in seven combinations, which define seven different haplotypes of the beta-globin gene: four major frameworks—1, 2, 3, and 3^*—and three minor frameworks, which we term KI1, KA1, and OR1. The nucleotide sequences of these frameworks are compared with those of homologous sequences in chimpanzee, colobus, macaque, and gorilla. This comparison provides strong evidence that framework 2 was the earliest framework in the human lineage. From framework 2, a rooted parsimonious tree for the six other frameworks is constructed. This phylogenetic tree is discussed in terms of the evolution of nucleotide polymorphisms as well as in terms of genetic affinities between human populations.For each position at which there is base difference in comparing human, gorilla, and chimpanzee beta-globin genes, the phyletic lineage where the corresponding substitution occurred has been identified using the maximum parsimony procedure. The data provide evidence that polymorphisms may represent a significant component of differences between closely related species. If so, nucleotide polymorphisms may strongly bias estimates of small evolutionary distances.     

Phylogenetic isolation of a human alu flounder gene: Drift to new subfamily identity

A severe bottleneck in the size of the PV Alu subfamily in the common ancestor of human and gorilla has been used to isolate an Alu source gene. The human PV Alu subfamily consists of about one thousand members which are absent in gorilla and chimpanzee DNA. Exhaustive library screening shows that there are as few as two PV Alus in the gorilla genome. One is gorilla-specific, i.e., absent in the orthologous loci in both human and chimpanzee, suggesting the independent retrotranspositional activity of the PV subfamily in the gorilla lineage. The second of these two gorilla PV Alus is present in both human and chimpanzee DNAs and is the single PV Alu known to precede the radiation of these three species. The orthologous Alu in gibbon DNA resembles the next older Alu subfamily. Thus, this Alu locus is originally templated by a non-PV source gene and acquired characteristic PV sequence variants by mutational drift in situ, consequently becoming the first member and presumptive founder of this PV subfamily. Correspondence to: C.W. Schmid