Conservation and loss of the ERV3 open reading frame in primates

The human endogenous retrovirus ERV3 possesses an open reading frame for a truncated envelope, which is expressed as mRNA and protein. Here we examine the env sequence in primates for evidence of evolutionary conservation. ERV3 sequences were amplified by PCR from genomic DNA of great ape and Old World primates but not from New World primates or gorilla, suggesting an integration event more than 30 million years ago with a subsequent loss in one species. In the chimpanzee, the protein sequence of Env is 98.18% identical to that of human. In other species the identity falls (93.71% in rhesus macaque) in proportion to the separation from the human lineage. Start and stop codons and domains of functional significance in the envelope protein are conserved. The evolutionary conservation of the ERV3 envelope suggests a beneficial function, though the loss from gorilla shows that it is not essential for survival or reproduction.     

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.     

Evolution of human cytoplasmic actin gene sequences: chromosome mapping and structural characterizations of three cytoplasmic actin-like pseudogenes including one on the Y chromosome

Eight recombinant phage clones containing cytoplasmic actin-like gene sequences have been isolated from a human genomic library for structural characterization. Kpn I family repeat sequences flank six of these actin genes isolated, and Alu family repeats are scattered throughout the DNA inserts of all eight phage clones. Three of these genes are γ actin-like, and the other five are β actin-like. The complete nucleotide sequence analysis of one β and one γ actin-like genes and their flanking regions demonstrates that they both are processed pseudogenes. Using unique DNA sequences flanking these two pseudogenes as hybridization probes for human-mouse somatic cell hybrid DNAs, we have mapped the two actin pseudogenes on human chromosomes 8 and 3, respectively. We have also determined the DNA sequence of a human Y chromosome-linked, processed actin pseudogene. The different values of sequence divergence of these processed pseudogenes and their functional counterparts allow us to estimate the time of generation of the pseudogenes. The results suggest that the cDNA insertion events generating the human cytoplasmic actin-like pseudogenes have occurred at significantly different times during the evolution of primates, after their separation from other mammalian species.     

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 type I hair keratin pseudogene ϕ hHaA has functional orthologs in the chimpanzee and gorilla: evidence for recent inactivation of the human gene after the Pan-Homo divergence

In addition to nine functional genes, the human type I hair keratin gene cluster contains a pseudogene, phihHaA (KRTHAP1), which is thought to have been inactivated by a single base-pair substitution that introduced a premature TGA termination codon into exon 4. Large-scale genotyping of human, chimpanzee, and gorilla DNAs revealed the homozygous presence of the phihHaA nonsense mutation in humans of different ethnic backgrounds, but its absence in the functional orthologous chimpanzee (cHaA) and gorilla (gHaA) genes. Expression analyses of the encoded cHaA and gHaA hair keratins served to highlight dramatic differences between the hair keratin phenotypes of contemporary humans and the great apes. The relative numbers of synonymous and non-synonymous substitutions in the phihHaA and cHaA genes, as inferred by using the gHaA gene as an outgroup, suggest that the human hHaA gene was inactivated only recently, viz., less than 240,000 years ago. This implies that the hair keratin phenotype of hominids prior to this date, and after the Pan-Homo divergence some 5.5 million years ago, could have been identical to that of the great apes. In addition, the homozygous presence of the phihHaA exon 4 nonsense mutation in some of the earliest branching lineages among extant human populations lends strong support to the "single African origin" hypothesis of modern humans.     

Structure and chromosome localization of the human eosinophil-derived neurotoxin and eosinophil cationic protein genes: Evidence for intronless coding sequences in the ribonuclease gene superfamily

Human genomic DNAs for the eosinophil granule proteins, eosinophil-derived neurotoxin (EDN) and eosinophil cationic protein (ECP), were isolated from genomic libraries. Alignment of EDN (RNS2) and ECP (RNS3) gene sequences demonstrated remarkable nucleotide similarities in noncoding sequences, introns, and flanking regions, as well as in the previously known coding regions. Detailed examination of the 5'-noncoding regions yielded putative TATA and CAAT boxes, as well as similarities to promoter motifs from unrelated genes. A single intron of 230 bases was found in the 5' untranslated region and we suggest that a single intron in this region and an intronless coding region are features common to many members of the RNase gene superfamily. The RNS2 and RNS3 genes were localized to the q24-q31 region of human chromosome 14. It is likely that these two genes arose as a consequence of a gene duplication event that took place approximately 25-40 million years ago and that a subset of anthropoid primates possess both of these genes or closely related genes.     

Accelerated Evolution of the ASPM Gene Controlling Brain Size Begins Prior to Human Brain Expansion

Primary microcephaly (MCPH) is a neurodevelopmental disorder characterized by global reduction in cerebral cortical volume. The microcephalic brain has a volume comparable to that of early hominids, raising the possibility that some MCPH genes may have been evolutionary targets in the expansion of the cerebral cortex in mammals and especially primates. Mutations in ASPM, which encodes the human homologue of a fly protein essential for spindle function, are the most common known cause of MCPH. Here we have isolated large genomic clones containing the complete ASPM gene, including promoter regions and introns, from chimpanzee, gorilla, orangutan, and rhesus macaque by transformation-associated recombination cloning in yeast. We have sequenced these clones and show that whereas much of the sequence of ASPM is substantially conserved among primates, specific segments are subject to high Ka/Ks ratios (nonsynonymous/synonymous DNA changes) consistent with strong positive selection for evolutionary change. The ASPM gene sequence shows accelerated evolution in the African hominoid clade, and this precedes hominid brain expansion by several million years. Gorilla and human lineages show particularly accelerated evolution in the IQ domain of ASPM. Moreover, ASPM regions under positive selection in primates are also the most highly diverged regions between primates and nonprimate mammals. We report the first direct application of TAR cloning technology to the study of human evolution. Our data suggest that evolutionary selection of specific segments of the ASPM sequence strongly relates to differences in cerebral cortical size.     

Accelerated evolution of the ASPM gene controlling brain size begins prior to human brain expansion

Primary microcephaly (MCPH) is a neurodevelopmental disorder characterized by global reduction in cerebral cortical volume. The microcephalic brain has a volume comparable to that of early hominids, raising the possibility that some MCPH genes may have been evolutionary targets in the expansion of the cerebral cortex in mammals and especially primates. Mutations in ASPM, which encodes the human homologue of a fly protein essential for spindle function, are the most common known cause of MCPH. Here we have isolated large genomic clones containing the complete ASPM gene, including promoter regions and introns, from chimpanzee, gorilla, orangutan, and rhesus macaque by transformation-associated recombination cloning in yeast. We have sequenced these clones and show that whereas much of the sequence of ASPM is substantially conserved among primates, specific segments are subject to high Ka/Ks ratios (nonsynonymous/synonymous DNA changes) consistent with strong positive selection for evolutionary change. The ASPM gene sequence shows accelerated evolution in the African hominoid clade, and this precedes hominid brain expansion by several million years. Gorilla and human lineages show particularly accelerated evolution in the IQ domain of ASPM. Moreover, ASPM regions under positive selection in primates are also the most highly diverged regions between primates and nonprimate mammals. We report the first direct application of TAR cloning technology to the study of human evolution. Our data suggest that evolutionary selection of specific segments of the ASPM sequence strongly relates to differences in cerebral cortical size.     

Composite human VK genes and a model of their evolution.

A phage library and two cosmid libraries were screened for human VK genes. Two recombinant phage and four cosmid clones were analysed in detail by restriction mapping and sequencing. Each one contained a single VKI sequence. Two of these six sequences are potentially functional VK genes and four are pseudogenes. Two pseudogenes derived from different genomic DNAs are highly homologous and are therefore either allelic variants or the products of a recent duplication event. Comparisons of our sequences with all fully determined human VKI amino acid and DNA sequences reveal identical segments which at first sight appear like minigenes. But these segments do not coincide with the subregions and some of the segments include both, framework and complementarity determining regions (FR, CDR, ref. 2). The findings may be explained by an evolutionary model generating composite genes by gene conversion and selection.     

Comparative genomics search for losses of long-established genes on the human lineage

Taking advantage of the complete genome sequences of several mammals, we developed a novel method to detect losses of well-established genes in the human genome through syntenic mapping of gene structures between the human, mouse, and dog genomes. Unlike most previous genomic methods for pseudogene identification, this analysis is able to differentiate losses of well-established genes from pseudogenes formed shortly after segmental duplication or generated via retrotransposition. Therefore, it enables us to find genes that were inactivated long after their birth, which were likely to have evolved nonredundant biological functions before being inactivated. The method was used to look for gene losses along the human lineage during the approximately 75 million years (My) since the common ancestor of primates and rodents (the euarchontoglire crown group). We identified 26 losses of well-established genes in the human genome that were all lost at least 50 My after their birth. Many of them were previously characterized pseudogenes in the human genome, such as GULO and UOX. Our methodology is highly effective at identifying losses of single-copy genes of ancient origin, allowing us to find a few well-known pseudogenes in the human genome missed by previous high-throughput genome-wide studies. In addition to confirming previously known gene losses, we identified 16 previously uncharacterized human pseudogenes that are definitive losses of long-established genes. Among them is ACYL3, an ancient enzyme present in archaea, bacteria, and eukaryotes, but lost approximately 6 to 8 Mya in the ancestor of humans and chimps. Although losses of well-established genes do not equate to adaptive gene losses, they are a useful proxy to use when searching for such genetic changes. This is especially true for adaptive losses that occurred more than 250,000 years ago, since any genetic evidence of the selective sweep indicative of such an event has been erased.     

Comparative and evolutionary analysis of the rhesus macaque extended MHC class II region

The sequence-based map of a part of the rhesus macaque major histocompatibility complex (MHC) extended class II region is presented. The sequenced region encompasses 67,401 bp and contains the SACM2L, RING1, FABGL and KE4 genes, as well as the HTATSF1-like and ZNF-like pseudogenes. Similar to human, but different from rat and mouse, no class I genes are found in the SACM2L- RING1 interval. The rhesus macaque extended MHC class II region shows a high degree of conservation of exonic as well as intronic and intergenic sequences compared with the respective human region. It is concluded that this particular genomic organization of the extended class II region-i.e., the absence of class I genes and the presence of the HTATSF1-like and ZNF-like pseudogenes-can be traced back to a common ancestor of humans and rhesus macaques about 23 million years ago.     

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.      

Baboon immunoglobulin constant region heavy chains: identification of four <Emphasis Type="Italic">IGHG</Emphasis> genes

The increasing use of nonhuman primate models in biomedical research and especially in vaccine development requires the characterization of their immunoglobulin genes and corresponding products. Therefore, we sequenced, cloned and characterized the four immunoglobulin gamma chain constant region genes ( IGHG) present in baboons. These four genes were designated IGHG1, IGHG2, IGHG3 and IGHG4 on the basis of sequence similarities with the four human genes encoding the IgG1, IgG2, IgG3 and IgG4 subclasses and the three known rhesus macaque IGHG genes. Specifically, the baboon IgG1, IgG2, IgG3 and IgG4 sequences exhibit 90.3%, 88.3%, 86.7% and 89.6% amino acid identity to their human counterpart. The percent of amino acid identity of baboon IgG1, IgG2 and IgG3 to the corresponding rhesus macaque sequences is 98.5, 93.1 and 94.4, respectively. Therefore, baboon and rhesus macaque IGHG genes are highly homologous to each other. The majority of differences existing between baboon and human sequences are clustered in the hinge region, with the upper hinge being the most diverse and containing several proline residues. Similar to rhesus macaques, the hinge regions of all baboon IGHG genes consist of a single exon, whereas in humans the IgG3 molecule is encoded by multiple exons. These results confirm the evolutionary instability of the hinge region and indicate that functional properties associated with the hinge regions of baboon and human IgG molecules might differ between the two species.     

Evolution of the β-globin gene cluster in man and the primates

The arrangement of primate β-related globin genes has been determined by restriction endonuclease mapping of genomic DNA from species ranging from prosimians to man. The arrangement of the entire ?γγδβ-globin gene cluster in the gorilla and the yellow baboon is indistinguishable from that of man. Restriction site differences between these species are consistent with a surprisingly low overall rate of intergenic DNA sequence divergence of approximately 1% in 5 million years. A new world monkey (owl monkey) has a single γ-globin gene, suggesting that the ^Gγ-^Aγ-globin gene duplication in man is ancient, and occurred about 20 to 40 million years ago. The β-globin gene cluster in the brown lemur, a prosimian, is remarkably short (about 20,000 base-pairs) and contains single ?-, γ- and β-globin genes. The γ- and β-globin genes in this animal are separated by a curious gene containing the 3′ end of a β-globin gene preceded by sequences related to the 5′ end of the ?-globin gene.     

Origin and phylogenetic distribution of Alu DNA repeats: irreversible events in the evolution of primates.

Over the past 60 million years, or so, approximately one million copies of Alu DNA repeats have accumulated in the genome of primates, in what appears to be an ongoing process. We determined the phylogenetic distribution of specific Alu (and other) DNA repeats in the genome of several primates: human, chimpanzee, gorilla, orangutan, baboon, rhesus, and macaque. At the population level studied, the majority of the repeats was found to be fixed in the primate species. Our data suggest that new Alu elements arise in unique, irreversible events, in a mechanism that seems to preclude precise excision and loss. The same insertions did not arise independently in two species. Once inserted and genetically fixed, the DNA elements are retained in all descendant lineages. The irreversible expansion of Alu s introduces a vector of time into the evolutionary process, and provides realistic (rather than statistical) answers to questions on phylogenies. In contrast to point mutations, the present distribution of individual Alu s is congruent with just one phylogeny. We submit that only irreversible and taxonomically relevant events are at the molecular basis of evolution. Most point mutations do not belong to this category.     

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.     

An anthropoid-specific segmental duplication on human chromosome 1q22

Segmental duplications (SDs) play a key role in genome evolution by providing material for gene diversification and creation of variant or novel functions. They also mediate recombinations, resulting in microdeletions, which have occasionally been associated with human genetic diseases. Here, we present a detailed analysis of a large genomic region (about 240 kb), located on human chromosome 1q22, that contains a tandem SD, SD1q22. This duplication occurred about 37 million years ago in a lineage leading to anthropoid primates, after their separation from prosimians but before the Old and New World monkey split. We reconstructed the hypothetical unduplicated ancestral locus and compared it with the extant SD1q22 region. Our data demonstrate that, as a consequence of the duplication, new anthropoid-specific genetic material has evolved in the resulting paralogous segments. We describe the emergence of two new genes, whose new functions could contribute to the speciation of anthropoid primates. Moreover, we provide detailed information regarding structure and evolution of the SD1q22 region that is a prerequisite for future studies of its anthropoid-specific functions and possible linkage to human genetic disorders.     

Molecular characterization of the HERV-W env gene in humans and primates: expression, FISH, phylogeny, and evolution

We characterized the human endogenous retrovirus (HERV-W) family in humans and primates. In silico expression data indicated that 22 complete HERV-W families from human chromosomes 1-3, 5-8, 10-12, 15, 19, and X are randomly expressed in various tissues. Quantitative real-time RT-PCR analysis of the HERV-W env gene derived from human chromosome 7q21.2 indicated predominant expression in the human placenta. Several copies of repeat sequences (SINE, LINE, LTR, simple repeat) were detected within the complete or processed pseudo HERV-W of the human, chimpanzee, and rhesus monkey. Compared to other regions (5'LTR, Gag, Gag-Pol, Env, 3'LTR), the repeat family has been mainly integrated into the region spanning the 5'LTRs of Gag (1398 bp) and Pol (3242 bp). FISH detected the HERV-W probe (fosWE1) derived from a gorilla fosmid library in the metaphase chromosomes of all primates (five hominoids, three Old World monkeys, two New World monkeys, and one prosimian), but not in Tupaia. This finding was supported by molecular clock and phylogeny data using the divergence values of the complete HERV-W LTR elements. The data suggested that the HERV-W family was integrated into the primate genome approximately 63 million years (Myr) ago, and evolved independently during the course of primate radiation.