Evolution of subterminal satellite (StSat) repeats in hominids

Subterminal satellite (StSat) repeats, consisting of 32-bp-long AT-rich units (GATATTTCCATGTT(T/C)ATACAGATAGCGGTGTA), were first found in chimpanzee and gorilla (African great apes) as one of the major components of heterochromatic regions located proximal to telomeres of chromosomes. StSat repeats have not been found in orangutan (Asian great ape) or human. This patchy distribution among species suggested that the StSat repeats were present in the common ancestor of African great apes and subsequently lost in the lineage leading to human. An alternative explanation is that the StSat repeats in chimpanzee and gorilla have different origins and the repeats did not occur in human. The purpose of the present study was quantitative evaluation of the above alternative possibilities by analyzing the nucleotide variation contained in the repeats. We collected large numbers of sequences of repeat units from genome sequence databases of chimpanzee and gorilla, and also bonobo (an African great ape phylogenetically closer to chimpanzee). We then compared the base composition of the repeat units among the 3 species, and found statistically significant similarities in the base composition. These results support the view that the StSat repeats had already formed multiple arrays in the common ancestor of African great apes. It is thus suggested that humans lost StSat repeats which had once grown to multiple arrays.     

Hotspots of Biased Nucleotide Substitutions in Human Genes

Genes that have experienced accelerated evolutionary rates on the human lineage during recent evolution are candidates for involvement in human-specific adaptations. To determine the forces that cause increased evolutionary rates in certain genes, we analyzed alignments of 10,238 human genes to their orthologues in chimpanzee and macaque. Using a likelihood ratio test, we identified protein-coding sequences with an accelerated rate of base substitutions along the human lineage. Exons evolving at a fast rate in humans have a significant tendency to contain clusters of AT-to-GC (weak-to-strong) biased substitutions. This pattern is also observed in noncoding sequence flanking rapidly evolving exons. Accelerated exons occur in regions with elevated male recombination rates and exhibit an excess of nonsynonymous substitutions relative to the genomic average. We next analyzed genes with significantly elevated ratios of nonsynonymous to synonymous rates of base substitution (d_N/d_S) along the human lineage, and those with an excess of amino acid replacement substitutions relative to human polymorphism. These genes also show evidence of clusters of weak-to-strong biased substitutions. These findings indicate that a recombination-associated process, such as biased gene conversion (BGC), is driving fixation of GC alleles in the human genome. This process can lead to accelerated evolution in coding sequences and excess amino acid replacement substitutions, thereby generating significant results for tests of positive selection.     

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

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

Genetic differences between humans and great apes

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

Human-specific subfamilies of HERV-K (HML-2) long terminal repeats: three master genes were active simultaneously during branching of hominoid lineages

Using 40 known human-specific LTR sequences, we have derived a consensus sequence for an evolutionary young HERV-K (HML-2) LTR family, which was named the HS family. In the human genome the HS family is represented by approximately 150-160 LTR sequences, 90% of them being human-specific (hs). The family can be subdivided into two subfamilies differing in five linked nucleotide substitutions: HS-a and HS-b of 5.8 and 10.3 Myr evolutionary ages, respectively. The HS-b subfamily members were transpositionally active both before the divergence of the human and chimpanzee ancestor lineages and after it in both lineages. The HS-a subfamily comprises only hs LTRs. These and other data strongly suggest that at least three "master genes" of HERV-K (HML-2) LTRs were active in the human ancestor lineage after the human-chimpanzee divergence. We also found hs HERV-K (HML-2) LTRs integrations in introns of 12 human genes and identified 13 new hs HERV-K (HML-2) LTRs.     

A second uniquely human mutation affecting sialic acid biology

Siglecs are immunoglobulin superfamily member lectins that selectively recognize different types and linkages of sialic acids, which are major components of cell surface and secreted glycoconjugates. We report here a human Siglec-like molecule (Siglec-L1) that lacks a conserved arginine residue known to be essential for optimal sialic acid recognition by previously known Siglecs. Loss of the arginine from an ancestral molecule was caused by a single nucleotide substitution that occurred after the common ancestor of humans with the great apes but before the origin of modern humans. The chimpanzee Siglec-L1 ortholog remains fully functional and preferentially recognizes N-glycolylneuraminic acid, which is a common sialic acid in great apes and other mammals. Reintroducing the ancestral arginine into the human molecule regenerates the same properties. Thus, the single base pair mutation that replaced the arginine on human Siglec-L1 is likely to be evolutionarily related to the previously reported loss of N-glycolylneuraminic acid expression in the human lineage. Siglec-L1 and its chimpanzee Siglec ortholog also have a different expression pattern from previously reported Siglecs because they are found on the lumenal edge of epithelial cell surfaces. Notably, the human genome contains several Siglec-like pseudogenes that have independent mutations that would have replaced the arginine residue required for optimal sialic acid recognition. Thus, additional changes in the biology of sialic acids may have taken place during human evolution.     

Dynamics of DNA Methylation in Recent Human and Great Ape Evolution

DNA methylation is an epigenetic modification involved in regulatory processes such as cell differentiation during development, X-chromosome inactivation, genomic imprinting and susceptibility to complex disease. However, the dynamics of DNA methylation changes between humans and their closest relatives are still poorly understood. We performed a comparative analysis of CpG methylation patterns between 9 humans and 23 primate samples including all species of great apes (chimpanzee, bonobo, gorilla and orangutan) using Illumina Methylation450 bead arrays. Our analysis identified ∼800 genes with significantly altered methylation patterns among the great apes, including ∼170 genes with a methylation pattern unique to human. Some of these are known to be involved in developmental and neurological features, suggesting that epigenetic changes have been frequent during recent human and primate evolution. We identified a significant positive relationship between the rate of coding variation and alterations of methylation at the promoter level, indicative of co-occurrence between evolution of protein sequence and gene regulation. In contrast, and supporting the idea that many phenotypic differences between humans and great apes are not due to amino acid differences, our analysis also identified 184 genes that are perfectly conserved at protein level between human and chimpanzee, yet show significant epigenetic differences between these two species. We conclude that epigenetic alterations are an important force during primate evolution and have been under-explored in evolutionary comparative genomics.     

Accelerated evolution of the pituitary adenylate cyclase-activating polypeptide precursor gene during human origin

Pituitary adenylate cyclase-activating polypeptide (PACAP) is a neuropeptide abundantly expressed in the central nervous system and involved in regulating neurogenesis and neuronal signal transduction. The amino acid sequence of PACAP is extremely conserved across vertebrate species, indicating a strong functional constraint during the course of evolution. However, through comparative sequence analysis, we demonstrated that the PACAP precursor gene underwent an accelerated evolution in the human lineage since the divergence from chimpanzees, and the amino acid substitution rate in humans is at least seven times faster than that in other mammal species resulting from strong Darwinian positive selection. Eleven human-specific amino acid changes were identified in the PACAP precursors, which are conserved from murine to African apes. Protein structural analysis suggested that a putative novel neuropeptide might have originated during human evolution and functioned in the human brain. Our data suggested that the PACAP precursor gene underwent adaptive changes during human origin and may have contributed to the formation of human cognition.     

Rh gene evolution in primates: study of intron sequences

By amplification and sequencing of RH gene intron 4 of various primates we demonstrate that an Alu-Sx-like element has been inserted in the RH gene of the common ancestor of humans, apes, Old World monkeys, and New World monkeys. The study of mouse and lemur intron 4 sequences allowed us to precisely define the insertion point of the Alu-Sx element in intron 4 of the RH gene ancestor common to Anthropoidea. Like humans, chimpanzees and gorillas possess two types of RH intron 4, characterized by the presence (human RHCE and ape RHCE-like genes) or absence (human RHD and ape RHD-like genes) of the Alu-Sx element. This led us to conclude that in the RH common ancestor of humans, chimpanzees, and gorillas, a duplication of the common ancestor gene gave rise to two genes, one differing from the other by a 654-bp deletion encompassing an Alu-Sx element. Moreover, most of chimpanzees and some gorillas posses two types of RHD-like intron 4. The introns 4 of type 1 have a length similar to that of human RHD intron 4, whereas introns 4 of type 2 display an insertion of 12 bp. The latest insertion was not found in the human genome (72 individuals tested). The study of RH intron 3 length polymorphism confirmed that, like humans, chimpanzees and gorillas possess two types of intron 3, with the RHD-type intron 3 being 289 bases shorter than the RHCE intron 3. By amplification and sequencing of regions encompassing introns 3 and 4, we demonstrated that chimpanzee and gorilla RH-like genes displayed associations of introns 3 and 4 distinct to those found in man. Altogether, the results demonstrate that, as in humans, chimpanzee and gorilla RH genes experienced intergenic exchanges.     

Cladistic relationships of extant and fossil hominoids

There is general agreement that the hominoid primates form a monophyletic group, that the extant great apes and humans form a second clade within that group with the gibbons as the sister group, and that the African apes and humans form a third clade. Although it has recently been proposed that humans and orang utans are sister taxa and also that the great apes form a clade to the exclusion of humans, our analysis, particularly of the molecular evidence, supports the existence of an African ape and human clade. The major problem in hominoid phylogeny at present is the relationships of the species within this clade: morphological data generally support the existence of an African ape clade which is the sister group to humans; some molecular data also support this conclusion, but most molecular evidence indicates the existence of a chimpanzee/human clade. We have cladistically re-analysed the DNA and protein sequence data for which apomorphic character states can be assessed. It is clear that there is a high degree of homoplasy whichever branching pattern is produced, with some characters supporting the existence of a chimpanzee/human clade and others supporting an African ape clade. When the cladistic analyses of morphological and molecular data are combined we believe that the most parsimonious interpretation of the data is that the African apes form a clade which is the sister taxon of the human (i.e., Australopithecus, Homo and Paranthropus) clade.This paper is not intended as a survey of all hominoid fossils but as a study of branching points in hominoid evolution and fossils are included which are relevant to this branching pattern. The analysis of fossil taxa in this study leads us to conclude that Proconsul is the sister taxon to the later Hominoidea. A number of middle Miocene forms such as Dryopithecus, Kenyapithecus, Heliopithecus and Afropithecus are shown to share derived characters with great apes and humans and provide evidence for the divergence of that clade from the gibbon lineage prior to 18 Ma. The position that Sivapithecus represents the sister group of the orang utan clade is supported here and shows that the orang utan lineage had diverged from the African ape and human lineage prior to 11·5 Ma. There is unfortunately no definitive fossil cvidence on branching sequences within the African ape and human clade, although a new specimen from Samburu, Kenya may be related to the gorilla.     

Detection of two distinct forms of apoC-I in great apes

ApoC-I, the smallest of the soluble apolipoproteins, associates with both TG-rich lipoproteins and HDL. Mass spectral analyses of human apoC-I previously had demonstrated that in the circulation there are two forms, either a 57 amino acid protein or a 55 amino acid protein, due to the loss of two amino acids from the N-terminus. In our analyses of the apolipoproteins of the other great apes by mass spectrometry, four forms of apoC-I were detected. Two of these showed a high degree of identity to the mature and truncated forms of human apoC-I. The other two were homologous to the virtual protein and its truncated form that are encoded by a human pseudogene. In humans, the genes for apoC-I and its pseudogene are located on chromosome 19, the pseudogene being 2.5 kb downstream from the apoC-I gene. Based on the similarity between the apoC-I gene and the pseudogene, it has been concluded that the latter arose from the former as a result of gene duplication approximately 35 Mya. Interestingly, the virtual protein encoded by the pseudogene is acidic, not basic like apoC-I. In the chimpanzee, there also are two genes for apoC-I, the one upstream encodes a basic protein and the downstream gene, rather than being a pseudogene, encodes an acidic protein (P86336). In addition to reporting on the molecular masses of great ape apoC-I, we were able to clearly demonstrate by “Top-down” sequencing that the acidic form arose from a separate gene. In our analyses, we have measured the molecular masses of apoC-I associated with the HDL of the following great apes: bonobo (Pan paniscus), chimpanzee (Pan troglodytes), and the Sumatran orangutan (Pongo abelii). Genomic variations in chromosome 19 among great apes, baboons and macaques as they relate to both genes for apoC-I and the pseudogene are compared and discussed.     

Lineage-Specific Changes in Biomarkers in Great Apes and Humans

Although human biomedical and physiological information is readily available, such information for great apes is limited. We analyzed clinical chemical biomarkers in serum samples from 277 wild- and captive-born great apes and from 312 healthy human volunteers as well as from 20 rhesus macaques. For each individual, we determined a maximum of 33 markers of heart, liver, kidney, thyroid and pancreas function, hemoglobin and lipid metabolism and one marker of inflammation. We identified biomarkers that show differences between humans and the great apes in their average level or activity. Using the rhesus macaques as an outgroup, we identified human-specific differences in the levels of bilirubin, cholinesterase and lactate dehydrogenase, and bonobo-specific differences in the level of apolipoprotein A-I. For the remaining twenty-nine biomarkers there was no evidence for lineage-specific differences. In fact, we find that many biomarkers show differences between individuals of the same species in different environments. Of the four lineage-specific biomarkers, only bilirubin showed no differences between wild- and captive-born great apes. We show that the major factor explaining the human-specific difference in bilirubin levels may be genetic. There are human-specific changes in the sequence of the promoter and the protein-coding sequence of uridine diphosphoglucuronosyltransferase 1 (UGT1A1), the enzyme that transforms bilirubin and toxic plant compounds into water-soluble, excretable metabolites. Experimental evidence that UGT1A1 is down-regulated in the human liver suggests that changes in the promoter may be responsible for the human-specific increase in bilirubin. We speculate that since cooking reduces toxic plant compounds, consumption of cooked foods, which is specific to humans, may have resulted in relaxed constraint on UGT1A1 which has in turn led to higher serum levels of bilirubin in humans.     

Molecular epidemiology of hepatitis B virus variants in nonhuman primates

We characterized hepatitis B virus (HBV) isolates from sera of 21 hepatitis B virus surface antigen-positive apes, members of the families Pongidae and Hylobatidae (19 gibbon spp., 1 chimpanzee, and 1 gorilla). Sera originate from German, French, Thai, and Vietnamese primate-keeping institutions. To estimate the phylogenetic relationships, we sequenced two genomic regions, one located within the pre-S1/pre-S2 region and one including parts of the polymerase and the X protein open reading frames. By comparison with published human and ape HBV isolates, the sequences could be classified into six genomic groups. Four of these represented new genomic groups of gibbon HBV variants. The gorilla HBV isolate was distantly related to the chimpanzee isolate described previously. To confirm these findings, the complete HBV genome from representatives of each genomic group was sequenced. The HBV isolates from gibbons living in different regions of Thailand and Vietnam could be classified into four different phylogenetically distinct genomic groups. The same genomic groups were found in animals from European zoos. Therefore, the HBV infections of these apes might have been introduced into European primate-keeping facilities by direct import of already infected animals from different regions in Thailand. Taken together, our data suggest that HBV infections are indigenous in the different apes. One event involving transmission between human and nonhuman primates in the Old World of a common ancestor of human HBV genotypes A to E and the ape HBV variants might have occurred.     

The DNA sequence of medaka chromosome LG22

We report the genomic DNA sequence of a single chromosome (linkage group 22; LG22) of the small teleost fish medaka (Oryzias latipes) as a first whole chromosome sequence from a non-mammalian vertebrate. The order and orientation of 633 protein-coding genes were deduced from 18,803,338 bp of DNA sequence, providing the opportunity to analyze chromosome evolution of vertebrate genomes by direct comparison with the human genome. The average number of genes in the "conserved gene cluster" (CGC), a strict definition of "synteny" at the sequence basis, between medaka and human was 1.6. These and other data suggest that approximately 38.8% of pair-wise gene relationships would have been broken from their common ancestor in the human and medaka lineages and further imply that approx 20,000 (15,520-23,280) breaks would have occurred from the entire genome of the common ancestor. These breaks were generated mainly by intra-chromosomal shufflings at a specific era in the vertebrate lineage. These precise comparative genomics allowed us to identify the pieces of ancient chromosomes of the common vertebrate ancestor and estimate chromosomal evolution in the vertebrate lineage.     

Structure,gene expression,and evolution of primate copper chaperone for superoxide dismutase

Copper chaperone for superoxide dismutase (CCS) is essential for transporting copper ion to Cu,Zn-superoxide dismutase (Cu,Zn-SOD). We cloned cDNAs for six primate species' CCSs. The total number of amino acid residues of primate CCSs is 274. Similarities between primates were over 96%. Important residues for the CCS function were well conserved. A phylogenetic tree of CCSs and Cu,Zn-SODs from various organisms showed that these two proteins were derived from a common ancestor, diverging very early on during eukaryote evolution. The high frequency of nonsynonymous substitutions was found in the lineage to Old World monkeys and apes. Expression of the CCS gene in various tissues of Japanese monkey was found to be high in the liver and adrenal gland, followed by the kidney and small intestine. Such expressional pattern was similar with that of Cu,Zn-SOD gene (Fukuhara et al., 2002).     

Detecting lineage-specific adaptive evolution of brain-expressed genes in human using rhesus macaque as outgroup

Comparative genetic analysis between human and chimpanzee may detect genetic divergences responsible for human-specific characteristics. Previous studies have identified a series of genes that potentially underwent Darwinian positive selection during human evolution. However, without a closely related species as outgroup, it is difficult to identify human-lineage-specific changes, which is critical in delineating the biological uniqueness of humans. In this study, we conducted phylogeny-based analyses of 2633 human brain-expressed genes using rhesus macaque as the outgroup. We identified 47 candidate genes showing strong evidence of positive selection in the human lineage. Genes with maximal expression in the brain showed a higher evolutionary rate in human than in chimpanzee. We observed that many immune-defense-related genes were under strong positive selection, and this trend was more prominent in chimpanzee than in human. We also demonstrated that rhesus macaque performed much better than mouse as an outgroup in identifying lineage-specific selection in humans.     

Evolution of eutherian cytochrome c oxidase subunit II: heterogeneous rates of protein evolution and altered interaction with cytochrome c

Cytochrome c oxidase subunit II (COII), encoded by the mitochondrial genome, exhibits one of the most heterogeneous rates of amino acid replacement among placental mammals. Moreover, it has been demonstrated that cytochrome c oxidase has undergone a structural change in higher primates which has altered its physical interaction with cytochrome c. We collected a large data set of COII sequences from several orders of mammals with emphasis on primates, rodents, and artiodactyls. Using phylogenetic hypotheses based on data independent of the COII gene, we demonstrated that an increased number of amino acid replacements are concentrated among higher primates. Incorporating approximate divergence dates derived from the fossil record, we find that most of the change occurred independently along the New World monkey lineage and in a rapid burst before apes and Old World monkeys diverged. There is some evidence that Old World monkeys have undergone a faster rate of nonsynonymous substitution than have apes. Rates of substitution at four-fold degenerate sites in primates are relatively homogeneous, indicating that the rate heterogeneity is restricted to nondegenerate sites. Excluding the rate acceleration mentioned above, primates, rodents, and artiodactyls have remarkably similar nonsynonymous replacement rates. A different pattern is observed for transversions at four-fold degenerate sites, for which rodents exhibit a higher rate of replacement than do primates and artiodactyls. Finally, we hypothesize specific amino acid replacements which may account for much of the structural difference in cytochrome c oxidase between higher primates and other mammals.