Toward a phylogenetic classification of Primates based on DNA evidence complemented by fossil evidence

Mitogenomic analyses of eutherian relationships

Reasonably correct phylogenies are fundamental to the testing of evolutionary hypotheses. Here, we present phylogenetic findings based on analyses of 67 complete mammalian mitochondrial (mt) genomes. The analyses, irrespective of whether they were performed at the amino acid (aa) level or on nucleotides (nt) of first and second codon positions, placed Erinaceomorpha (hedgehogs and their kin) as the sister group of remaining eutherians. Thus, the analyses separated Erinaceomorpha from other traditional lipotyphlans (e.g., tenrecs, moles, and shrews), making traditional Lipotyphla polyphyletic. Both the aa and nt data sets identified the two order-rich eutherian clades, the Cetferungulata (comprising Pholidota, Carnivora, Perissodactyla, Artiodactyla, and Cetacea) and the African clade (Tenrecomorpha, Macroscelidea, Tubulidentata, Hyracoidea, Proboscidea, and Sirenia). The study corroborated recent findings that have identified a sister-group relationship between Anthropoidea and Dermoptera (flying lemurs), thereby making our own order, Primates, a paraphyletic assembly. Molecular estimates using paleontologically well-established calibration points, placed the origin of most eutherian orders in Cretaceous times, 70-100 million years before present (MYBP). The same estimates place all primate divergences much earlier than traditionally believed. For example, the divergence between Homo and Pan is estimated to have taken place approximately 10 MYBP, a dating consistent with recent findings in primate paleontology.

Natural selection shaped regional mtDNA variation in humans

Human mtDNA shows striking regional variation, traditionally attributed to genetic drift. However, it is not easy to account for the fact that only two mtDNA lineages (M and N) left Africa to colonize Eurasia and that lineages A, C, D, and G show a 5-fold enrichment from central Asia to Siberia. As an alternative to drift, natural selection might have enriched for certain mtDNA lineages as people migrated north into colder climates. To test this hypothesis we analyzed 104 complete mtDNA sequences from all global regions and lineages. African mtDNA variation did not significantly deviate from the standard neutral model, but European, Asian, and Siberian plus Native American variations did. Analysis of amino acid substitution mutations (nonsynonymous, Ka) versus neutral mutations (synonymous, Ks) (kaks) for all 13 mtDNA protein-coding genes revealed that the ATP6 gene had the highest amino acid sequence variation of any human mtDNA gene, even though ATP6 is one of the more conserved mtDNA proteins. Comparison of the kaks ratios for each mtDNA gene from the tropical, temperate, and arctic zones revealed that ATP6 was highly variable in the mtDNAs from the arctic zone, cytochrome b was particularly variable in the temperate zone, and cytochrome oxidase I was notably more variable in the tropics. Moreover, multiple amino acid changes found in ATP6, cytochrome b, and cytochrome oxidase I appeared to be functionally significant. From these analyses we conclude that selection may have played a role in shaping human regional mtDNA variation and that one of the selective influences was climate.

Orangutan cultures and the evolution of material culture

Geographic variation in some aspects of chimpanzee behavior has been interpreted as evidence for culture. Here we document similar geographic variation in orangutan behaviors. Moreover, as expected under a cultural interpretation, we find a correlation between geographic distance and cultural difference, a correlation between the abundance of opportunities for social learning and the size of the local cultural repertoire, and no effect of habitat on the content of culture. Hence, great-ape cultures exist, and may have done so for at least 14 million years.

Historical contingency in the evolution of primate color vision

Primates are unique among eutherian mammals for possessing three types of retinal cone. Curiously, catarrhines, platyrrhines, and strepsirhines share this anatomy to different extents, and no hypothesis has hitherto accounted for this variability. Here we propose that the historical biogeography of figs and arborescent palms accounts for the global variation in primate color vision. Specifically, we suggest that primates invaded Paleogene forests characterized by figs and palms, the fruits of which played a keystone function. Primates not only relied on such resources, but also provided high-quality seed dispersal. In turn, figs and palms lost or simply did not evolve conspicuous coloration, as this conferred little advantage for attracting mammals. We suggest that the abundance and coloration of figs and palms offered a selective advantage to foraging groups with mixed capabilities for chromatic distinction. Climatic cooling at the end of the Eocene and into the Neogene resulted in widespread regional extinction or decimation of palms and (probably) figs. In regions where figs and palms became scarce, we suggest primates evolved routine trichromatic vision in order to exploit proteinaceous young leaves as a replacement resource. A survey of the hue and biogeography of extant figs and palms provides some empirical support. Where these resources are infrequent, primates are routinely trichromatic and consume young leaves during seasonal periods of fruit dearth. These results imply a link between the differential evolution of primate color vision and climatic changes during the Eocene-Oligocene transition.

Phylogenetic analyses of Zostera species based on rbcL and matK nucleotide sequences: Implications for the origin and diversification of seagrasses in Japanese waters

Seagrasses are composed of four families belonging to angiosperms and they are thought to become adaptive to aquatic life independently. Zosteraceae is one such family and because of the relatively high species diversity around Japan and Korea coast areas, the family might have arisen therefrom. To elucidate the origin and evolution of Zosteraceae which consists of three genera, Phyllospadix, Zostera, and Heterozostera, 2.8 kb nucleotide sequences of rbcL and matK genes in the chloroplast genome were examined for various species, including cosmopolitan Z. marina and endemic Z. caulescens. The phylogenetic analysis reveals the following three features. First, based on the synonymous nucleotide substitution rate of the rice chloroplast genome, we estimated the divergence times between Zosteraceae and its closest relative, Potamogetonaceae, and between different genera, Zostera and Phyllospadix, as approximately 100 million years (myr) and 36 myr, respectively, suggesting that Zosteraceae emerged somewhere in the period from 36 myr ago to 100 myr ago. Second, two subgenera of Zostera, Zostera and Zosterella, exhibit their reciprocal monophyly and appear to have differentiated from each other approximately 33 myr ago. However, the third genus Heterozostera branched off only 5 myr ago from the stem lineage leading to Zosterella and this seems too recent in comparison with the ancient divergence of the two subgenera. Third, we estimated the most recent common ancestor of subgenus Zostera as 6 myr. In Z. marina four haplotypes were found in the sample and have diversified in the past 1.5 myr. One haplotype is shared by both sides of the Japan Archipelago and its closely related haplotypes occur also in eastern Pacific Ocean. Based on these phylogeographic analyses, we propose a provisional age related classification of Zosteraceae to argue the origin and evolution.

Molecular systematics and biogeography of the Neotropical monkey genus, Alouatta

We take advantage of the broad distribution of howler monkeys from Mexico to Argentina to provide a historical biogeographical analysis on a regional scale that encompasses the entire Neotropics. The phylogenetic relationships among 9 of the 10 recognized Alouatta species were inferred using three mitochondrial and two nuclear genes. The nuclear gene regions provided no phylogenetic resolution among howler monkey species, and were characterized by very low levels of sequence divergence between Alouatta and the Ateles outgroup. The mtDNA genes, on the other hand, produced a well-resolved phylogeny, which indicated that the earliest split among howler monkeys separated cis- and trans-Andean clades. Eight monophyletic mtDNA haplotype clades were identified, representing six named species in South America, including Alouatta seniculus, Alouatta sara, Alouatta macconelli, Alouatta caraya, Alouatta belzebul, and Alouatta guariba, and two in Mesoamerica, Alouatta pigra and Alouatta palliata. Molecular clock-based estimates of branching times indicated that contemporary howler monkey species originated in the late Miocene and Pliocene, not the Pleistocene. The causes of Alouatta diversification were more difficult to pin down, although we posit that the initial cis-, trans-Andean split in the genus was caused by the late Miocene completion of the northern Andes. Riverine barriers to dispersal and putative forest refuges can neither be discounted nor distinguished as causes of speciation in many cases, and one, the other or both have likely played a role in the diversification of South American howler monkeys. Finally, we estimated the separation of Mesoamerican A. pigra and A. palliata at 3Ma, which corresponds to the completion date of the Panama Isthmus promoting a role for this earth history event in the speciation of Central American howler monkeys.

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

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

Quantifying the slightly deleterious mutation model of molecular evolution

We have attempted to quantify the frequency and effects of slightly deleterious mutations (SDMs), those that have selective effects close to the reciprocal of the effective population size of a species, by comparing the level of selective constraint in protein-coding genes of related species that have different present-day effective population sizes. In our two comparisons, the species with the smaller effective population size showed lower constraint, implying that SDMs had become fixed. The fixation of SDMs was supported by the observation of a higher fraction of radical to conservative amino acid substitutions in species with smaller effective population sizes. The fraction of strongly deleterious mutations (which rarely become fixed) is >70% in most species. Only approximately 10% or fewer of mutations seem to behave as SDMs, but SDMs could comprise a substantial fraction of mutations in protein-coding genes that have a chance of becoming fixed between species.

Adaptive loss of ultraviolet-sensitive/violet-sensitive (UVS/VS) cone opsin in the blind mole rat (Spalax ehrenbergi)

In previous studies, fully functional rod and long-wavelength-sensitive (LWS) cone photopigments have been isolated from the eye of the subterranean blind mole rat (Spalax ehrenbergi superspecies). Spalax possesses subcutaneous atrophied eyes and lacks any ability to respond to visual images. By contrast this animal retains the ability to entrain circadian rhythms of locomotor behaviour to environmental light cues. As this is the only known function of the eye, the rod and LWS photopigments are thought to mediate this response. Most mammals are dichromats possessing, in addition to a single rod photopigment, two classes of cone photopigment, LWS and ultraviolet-sensitive/violet-sensitive (UVS/VS) with differing spectral sensitivities which mediate colour vision. In this paper we explore whether Spalax is a dichromat and has the potential to use colour discrimination for photoentrainment. Using immunocytochemistry and molecular approaches we demonstrate that Spalax is a LWS monochromat. Spalax lacks a functional UVS/VS cone photopigment due to the accumulation of several deleterious mutational changes that have rendered the gene nonfunctional. Using phylogenetic analysis we show that the loss of this class of photoreceptor is likely to have arisen from the visual ecology of this species, and is not an artefact of having an ancestor which lacked a functional UVS/VS cone photopigment. We conclude that colour discrimination is not a prerequisite for photoentrainment in this species.

Identification and analysis of over 2000 ribosomal protein pseudogenes in the human genome

Mammals have 79 ribosomal proteins (RP). Using a systematic procedure based on sequence-homology, we have comprehensively identified pseudogenes of these proteins in the human genome. Our assignments are available at http://www.pseudogene.org or http://bioinfo.mbb.yale.edu/genome/pseudogene. In total, we found 2090 processed pseudogenes and 16 duplications of RP genes. In relation to the matching parent protein, each of the processed pseudogenes has an average relative sequence length of 97% and an average sequence identity of 76%. A small number (258) of them do not contain obvious disablements (stop codons or frameshifts) and, therefore, could be mistaken as functional genes, and 178 are disrupted by one or more repetitive elements. On average, processed pseudogenes have a longer truncation at the 5' end than the 3' end, consistent with the target-primed-reverse-transcription (TPRT) mechanism. Interestingly, on chromosome 16, an RPL26 processed pseudogene was found in the intron region of a functional RPS2 gene. The large-scale distribution of RP pseudogenes throughout the genome appears to result, chiefly, from random insertions with the numbers on each chromosome, consequently, proportional to its size. In contrast to RP genes, the RP pseudogenes have the highest density in GC-intermediate regions (41%-46%) of the genome, with the density pattern being between that of LINEs and Alus. This can be explained by a negative selection theory as we observed that GC-rich RP pseudogenes decay faster in GC-poor regions. Also, we observed a correlation between the number of processed pseudogenes and the GC content of the associated functional gene, i.e., relatively GC-poor RPs have more processed pseudogenes. This ranges from 145 pseudogenes for RPL21 down to 3 pseudogenes for RPL14. We were able to date the RP pseudogenes based on their sequence divergence from present-day RP genes, finding an age distribution similar to that for Alus. The distribution is consistent with a decline in retrotransposition activity in the hominid lineage during the last 40 Myr. We discuss the implications for retrotransposon stability and genome dynamics based on these new findings.

Structure, molecular evolution, and gene expression of primate superoxide dismutases

Mn- and Cu,Zn-superoxide dismutase (SOD) cDNAs of eight primate species, Pan troglodytes, Pongo pygmaeus, Hylobates lar, Macaca fuscata, Macaca fascicularis, Macaca mulatta, Cebus apella, and Callithrix jacchus, were cloned. The whole protein-coding sequences were covered, comparing 198 and 153 (or 154) amino acids, for Mn- and Cu,Zn-SODs, respectively. Residues forming metal ligands were completely conserved in the two primate SODs and nucleotide/amino acid substitutions were more frequent in Cu,Zn-SODs than in Mn-SODs. Molecular evolutionary analyses showed Mn-SOD to have evolved at a constant rate and its phylogenetic tree well reflected primate phylogeny. Cu,Zn-SOD was shown to have evolved differently between primate lineages. The significant high ratio of a non-synonymous/synonymous rate was found in the lineage leading to great apes and humans, showing that this lineage underwent positive Darwinian selection. Southern hybridization suggested that the genes for primate Mn- and Cu,Zn-SOD exist as single copies. Northern analysis in various Japanese monkey tissues showed Mn- and Cu,Zn-SOD expression to be high in the liver, kidneys, and adrenal glands.

Direct evidence for the Homo-Pan clade

For a long time, the evolutionary relationship between human and African apes, the 'trichotomy problem', has been debated with strong differences in opinion and interpretation. Statistical analyses of different molecular DNA data sets have been carried out and have primarily supported a Homo-Pan clade. An alternative way to address this question is by the comparison of evolutionarily relevant chromosomal breakpoints. Here, we made use of a P1-derived artificial chromosome (PAC)/bacterial artificial chromosome (BAC) contig spanning approximately 2.8 Mb on the long arm of the human Y chromosome, to comparatively map individual PAC clones to chromosomes from great apes, gibbons, and two species of Old World monkeys by fluorescence in-situ hybridization. During our search for evolutionary breakpoints on the Y chromosome, it transpired that a transposition of an approximately 100-kb DNA fragment from chromosome 1 onto the Y chromosome must have occurred in a common ancestor of human, chimpanzee and bonobo. Only the Y chromosomes of these three species contain the chromosome-1-derived fragment; it could not be detected on the Y chromosomes of gorillas or the other primates examined. Thus, this shared derived (synapomorphic) trait provides clear evidence for a Homo-Pan clade independent of DNA sequence analysis.

Chromosomal homologies between Cebus and Ateles (primates) based on ZOO-FISH and G-banding comparisons

ZOO-FISH (Fluorescent "in vitro" hybridization) was used to establish the chromosomal homology between humans (HSA) and Cebus nigrivitatus (CNI) and Ateles belzebuth hybridus (ABH). These two species belong to different New World monkey families (Cebidae and Atelidae, respectively) which differ greatly in chromosome number and in chromosome morphology. The molecular results were followed by a detailed banding analysis. The ancestral karyotype of Cebus was then determined by a comparison of in situ hybridization results, as well as chromosomal morphology and banding in other Platyrrhini species. The karyotypes of the four species belonging to the genus Cebus differ from each other by three inversions and one fusion as well as in the location and amounts of heterochromatin. Results obtained by ZOO-FISH in ABH are in general agreement with previous gene-mapping and in situ hybridization data in Ateles, which show that spider monkeys have highly derived genomes. The chromosomal rearrangements detected between HSA and ABH on a band-to-band basis were 27 fusions/fissions, 12 centromeric shifts, and six pericentric inversions. The ancestral karyotype of Cebus was then compared with that of Ateles. The rearrangements detected were 20 fusions/fissions, nine centromeric shifts, and five inversions. Atelidae species are linked by a fragmentation of chromosome 4 into three segments forming an association of 4/15, while Ateles species are linked by 13 derived associations. The results also helped clarify the content of the ancestral platyrrhine karyotype and the mode of chromosomal evolution in these primates. In particular, associations 2/16 and 5/7 should be included in the ancestral karyotype of New World monkeys.

Cloning and characterization of human Siglec-11. A recently evolved signaling molecule that can interact with SHP-1 and SHP-2 and is expressed by tissue macrophages, including brain microglia

Siglecs are sialic acid-recognizing animal lectins of the immunoglobulin superfamily. We have cloned and characterized a novel human molecule, Siglec-11, that belongs to the subgroup of CD33/Siglec-3-related Siglecs. As with others in this subgroup, the cytosolic domain of Siglec-11 is phosphorylated at tyrosine residue(s) upon pervanadate treatment of cells and then recruits the protein-tyrosine phosphatases SHP-1 and SHP-2. However, Siglec-11 has several novel features relative to the other CD33/Siglec-3-related Siglecs. First, it binds specifically to alpha2-8-linked sialic acids. Second, unlike other CD33/Siglec-3-related Siglecs, Siglec-11 was not found on peripheral blood leukocytes. Instead, we observed its expression on macrophages in various tissues, such as liver Kupffer cells. Third, it was also expressed on brain microglia, thus becoming the second Siglec to be found in the nervous system. Fourth, whereas the Siglec-11 gene is on human chromosome 19, it lies outside the previously described CD33/Siglec-3-related Siglec cluster on this chromosome. Fifth, analyses of genome data bases indicate that Siglec-11 has no mouse ortholog and that it is likely to be the last canonical human Siglec to be reported. Finally, although Siglec-11 shows marked sequence similarity to human Siglec-10 in its extracellular domain, the cytosolic tail appears only distantly related. Analysis of genomic regions surrounding the Siglec-11 gene suggests that it is actually a chimeric molecule that arose from relatively recent gene duplication and recombination events, involving the extracellular domain of a closely related ancestral Siglec gene (which subsequently became a pseudogene) and a transmembrane and cytosolic tail derived from another ancestral Siglec.

Functional promiscuity of squirrel monkey growth hormone receptor toward both primate and nonprimate growth hormones

Primate growth hormone (GH) has evolved rapidly, having undergone approximately 30% amino acid substitutions from the inferred ancestral eutherian sequence. Nevertheless, human growth hormone (hGH) is physiologically effective when administered to nonprimate mammals. In contrast, its functional counterpart, the human growth hormone receptor (hGHR), has evolved species specificity so that it responds only to Old World primate GHs. It has been proposed that this species specificity of the hGHR is largely caused by the Leu --> Arg change at position 43 after a prior His --> Asp change at position 171 of the GH. Sequence analyses supported this hypothesis and revealed that the transitional phase in the GH:GHR coevolution still persists in New World monkeys. For example, although the GH of the squirrel monkey has the His --> Asp substitution at position 171, residue 43 of its GHR is a Leu, the nonprimate residue. If the squirrel monkey truly represents an intermediate stage of GH:GHR coevolution, its GHR should respond to both hGH and nonprimate GH. Also, if the emergence of species specificity was a result of the selection for a more efficient GH:GHR interaction, then changing residue 43 of the squirrel monkey growth hormone receptor (smGHR) to Arg should increase its binding affinity toward higher primate GH. To test these hypotheses, we performed protein-binding assays between the smGHR and both human and rat GHs, using the surface plasmon resonance methodology. Furthermore, the effects of reciprocal mutations at position 43 of human and squirrel monkey GHRs are measured for their binding affinities toward human and squirrel monkey GHs. The results from the binding kinetic assays clearly demonstrate that the smGHR is in the intermediate state of the evolution of species specificity. Interestingly, the altered residue Arg at position 43 of the smGHR does not lead to an increased binding affinity. The implications of these results on the evolution of the GH:GHR interaction and on functional evolution are discussed.

A map of the common chimpanzee genome

The completion of the chimpanzee genome will greatly help us determine which genetic changes are unique to humanity. Chimpanzees are our closest living relative, and a recent study has made considerable progress towards decoding the genome of our sister taxon.1 Over 75,000 common chimpanzee (Pan troglodytes) bacterial artificial chromosome end sequences were aligned and mapped to the human genome. This study shows the remarkable genetic similarity (98.77%) between humans and chimpanzees, while highlighting intriguing areas of potential difference. If we wish to understand the genetic basis of humankind, the completion of the chimpanzee genome deserves high priority.

Search for genes positively selected during primate evolution by 5'-end-sequence screening of cynomolgus monkey cDNAs

It is possible to assess positive selection by using the ratio of K(a) (nonsynonymous substitutions per plausible nonsynonymous sites) to K(s) (synonymous substitutions per plausible synonymous sites). We have searched candidate genes positively selected during primate evolution by using 5'-end sequences of 21,302 clones derived from cynomolgus monkey (Macaca fascicularis) brain cDNA libraries. Among these candidates, 10 genes that had not been shown by previous studies to undergo positive selection exhibited a K(a)/K(s) ratio > 1. Of the 10 candidate genes we found, 5 were included in the mitochondrial respiratory enzyme complexes, suggesting that these nuclear-encoded genes coevolved with mitochondrial-encoded genes, which have high mutation rates. The products of other candidate genes consisted of a cell-surface protein, a member of the lipocalin family, a nuclear transcription factor, and hypothetical proteins.

Biomedical applications and studies of molecular evolution: a proposal for a primate genomic library resource

The anticipated completion of two of the most biomedically relevant genomes, mouse and human, within the next three years provides an unparalleled opportunity for the large-scale exploration of genome evolution. Targeted sequencing of genomic regions in a panel of primate species and comparison to reference genomes will provide critical insight into the nature of single-base pair variation, mechanisms of chromosomal rearrangement, patterns of selection, and species adaptation. Although not recognized as model "genetic organisms" because of their longevity and low fecundity, 30 of the approximately 300 primate species are targets of biomedical research. The existence of a human reference sequence and genomic primate BAC libraries greatly facilitates the recovery of genes/genomic regions of high biological interest because of an estimated maximum neutral nucleotide sequence divergence of 25%. Primate species, therefore, may be regarded as the ideal model "genomic organisms". Based on existing BAC library resources, we propose the construction of a panel of primate BAC libraries from phylogenetic anchor species for the purpose of comparative medicine as well as studies of genome evolution.

Evolution of the phosphoglycerate mutase processed gene in human and chimpanzee revealing the origin of a new primate gene

Processed genes are created by retroposition from messenger RNA of expressed genes. The estimated amount of processed copies of genes in the human genome is 10,000-14,000. Some of these might be pseudogenes with the expected pattern for nonfunctional sequences, but some others might be an important source of new genes. We have studied the evolution of a Phosphoglycerate mutase processed gene (PGAM3) described in humans and believed to be a pseudogene. We sequenced PGAM3 in chimpanzee and macaque and obtained polymorphism data for human coding region. We found evidence that PGAM3 likely produces a functional protein, as an example of addressing functionality for human processed pseudogenes. First, the open reading frame was intact despite many deletions that occurred in the 3' untranslated region. Second, it appears that the gene is expressed. Finally, interspecies and intraspecies variation for PGAM3 was not consistent with the neutral model proposed for pseudogenes, suggesting that a new functional primate gene has originated. Amino acid divergence was significantly higher than synonymous divergence in PGAM3 lineage, supporting positive selection acting in this gene. This role of selection was further supported by the excess of rare alleles in a population genetic analysis. PGAM3 is located in a region of very low recombination; therefore, it is conceivable that the rapid fixation events in this newly arising gene may have contributed to a selective sweep of variation in the region.

Complementary advantageous substitutions in the evolution of an antiviral RNase of higher primates

An improved understanding of the evolution of gene function at the molecular level may provide significant insights into the origin of biological novelty and adaptation. With the approach of ancestral protein reconstruction, we here address the question of how a dramatically enhanced ribonucleolytic activity and the related antiviral activity evolved in a recently duplicated ribonuclease (eosinophil-derived neurotoxin) gene of higher primates. We show that the mother gene of the duplicated genes had already possessed a weak antiviral activity before duplication. After duplication, substitutions at two interacting sites (Arg-64-->Ser and Thr-132-->Arg) resulted in a 13-fold enhancement of the ribonucleolytic activity of eosinophil-derived neurotoxin. These substitutions are also necessary for the potent antiviral activity, with contributions from additional amino acid changes at interacting sites. Our observation that a change in eosinophil-derived neurotoxin function occurs only when both interacting sites are altered indicates the importance of complementary substitutions in protein evolution. Thus, neutral substitutions are not simply "noises" in protein evolution, as many have thought. They may play constructive roles by setting the intramolecular microenvironment for further complementary advantageous substitutions, which can lead to improved or altered function. Overall, our study illustrates the power of the "paleomolecular biochemistry" approach in delineating the complex interplays of amino acid substitutions in evolution and in identifying the molecular basis of biological innovation.

The complete mitochondrial sequence of Tarsius bancanus: evidence for an extensive nucleotide compositional plasticity of primate mitochondrial DNA

Inconsistencies between phylogenetic interpretations obtained from independent sources of molecular data occasionally hamper the recovery of the true evolutionary history of certain taxa. One prominent example concerns the primate infraordinal relationships. Phylogenetic analyses based on nuclear DNA sequences traditionally represent Tarsius as a sister group to anthropoids. In contrast, mitochondrial DNA (mtDNA) data only marginally support this affiliation or even exclude Tarsius from primates. Two possible scenarios might cause this conflict: a period of adaptive molecular evolution or a shift in the nucleotide composition of higher primate mtDNAs through directional mutation pressure. To test these options, the entire mt genome of Tarsius bancanus was sequenced and compared with mtDNA of representatives of all major primate groups and mammals. Phylogenetic reconstructions at both the amino acid (AA) and DNA level of the protein-coding genes led to faulty tree topologies depending on the algorithms used for reconstruction. We propose that these artifactual affiliations rather reflect the nucleotide compositional similarity than phylogenetic relatedness and favor the directional mutation pressure hypothesis because: (1) the overall nucleotide composition changes dramatically on the lineage leading to higher primates at both silent and nonsilent sites, and (2) a highly significant correlation exists between codon usage and the nucleotide composition at the third, silent codon position. Comparisons of mt genes with mt pseudogenes that presumably transferred to the nucleus before the directional mutation pressure took place indicate that the ancestral DNA composition is retained in the relatively fossilized mtDNA-like sequences, and that the directed acceleration of the substitution rate in higher primates is restricted to mtDNA.

Adaptive evolution of a duplicated pancreatic ribonuclease gene in a leaf-eating monkey

Although the complete genome sequences of over 50 representative species have revealed the many duplicated genes in all three domains of life, the roles of gene duplication in organismal adaptation and biodiversity are poorly understood. In addition, the evolutionary forces behind the functional divergence of duplicated genes are often unknown, leading to disagreement on the relative importance of positive Darwinian selection versus relaxation of functional constraints in this process. The methodology of earlier studies relied largely on DNA sequence analysis but lacked functional assays of duplicated genes, frequently generating contentious results. Here we use both computational and experimental approaches to address these questions in a study of the pancreatic ribonuclease gene (RNASE1) and its duplicate gene (RNASE1B) in a leaf-eating colobine monkey, douc langur. We show that RNASE1B has evolved rapidly under positive selection for enhanced ribonucleolytic activity in an altered microenvironment, a response to increased demands for the enzyme for digesting bacterial RNA. At the same time, the ability to degrade double-stranded RNA, a non-digestive activity characteristic of primate RNASE1, has been lost in RNASE1B, indicating functional specialization and relaxation of purifying selection. Our findings demonstrate the contribution of gene duplication to organismal adaptation and show the power of combining sequence analysis and functional assays in delineating the molecular basis of adaptive evolution.

The effects of nucleotide substitution model assumptions on estimates of nonparametric bootstrap support

The use of parameter-rich substitution models in molecular phylogenetics has been criticized on the basis that these models can cause a reduction both in accuracy and in the ability to discriminate among competing topologies. We have explored the relationship between nucleotide substitution model complexity and nonparametric bootstrap support under maximum likelihood (ML) for six data sets for which the true relationships are known with a high degree of certainty. We also performed equally weighted maximum parsimony analyses in order to assess the effects of ignoring branch length information during tree selection. We observed that maximum parsimony gave the lowest mean estimate of bootstrap support for the correct set of nodes relative to the ML models for every data set except one. For several data sets, we established that the exact distribution used to model among-site rate variation was critical for a successful phylogenetic analysis. Site-specific rate models were shown to perform very poorly relative to gamma and invariable sites models for several of the data sets most likely because of the gross underestimation of branch lengths. The invariable sites model also performed poorly for several data sets where this model had a poor fit to the data, suggesting that addition of the gamma distribution can be critical. Estimates of bootstrap support for the correct nodes often increased under gamma and invariable sites models relative to equal rates models. Our observations are contrary to the prediction that such models cause reduced confidence in phylogenetic hypotheses. Our results raise several issues regarding the process of model selection, and we briefly discuss model selection uncertainty and the role of sensitivity analyses in molecular phylogenetics.

The solitary long terminal repeats of ERV-9 endogenous retrovirus are conserved during primate evolution and possess enhancer activities in embryonic and hematopoietic cells

The solitary long terminal repeats (LTRs) of ERV-9 endogenous retrovirus contain the U3, R, and U5 regions but no internal viral genes. They are middle repetitive DNAs present at 2,000 to 4,000 copies in primate genomes. Sequence analyses of the 5" boundary area of the erythroid beta-globin locus control region (beta-LCR) and the intron of the embryonic axin gene show that a solitary ERV-9 LTR has been stably integrated in the respective loci for at least 15 million years in the higher primates from orangutan to human. Functional studies utilizing the green fluorescent protein (GFP) gene as the reporter in transfection experiments show that the U3 region of the LTRs possesses strong enhancer activity in embryonic cells of widely different tissue origins and in adult cells of blood lineages. In both the genomic LTRs of embryonic placental cells and erythroid K562 cells and transfected LTRs of recombinant GFP plasmids in K562 cells, the U3 enhancer activates synthesis of RNAs that are initiated from a specific site 25 bases downstream of the AATAAA (TATA) motif in the U3 promoter. A second AATAAA motif in the R region does not serve as the TATA box or as the polyadenylation signal. The LTR-initiated RNAs extend through the R and U5 regions into the downstream genomic DNA. The results suggest that the ERV-9 LTR-initiated transcription process may modulate transcription of the associated gene loci in embryonic and hematopoietic cells.

Vertebrate genome sequencing: building a backbone for comparative genomics

The human genome sequence provides a reference point from which we can compare ourselves with other organisms. Interspecies comparison is a powerful tool for inferring function from genomic sequence and could ultimately lead to the discovery of what makes humans unique. To date, most comparative sequencing has focused on pair-wise comparisons between human and a limited number of other vertebrates, such as mouse. Targeted approaches now exist for mapping and sequencing vertebrate bacterial artificial chromosomes (BACs) from numerous species, allowing rapid and detailed molecular and phylogenetic investigation of multi-megabase loci. Such targeted sequencing is complementary to current whole-genome sequencing projects, and would benefit greatly from the creation of BAC libraries from a diverse range of vertebrates.

Statistical genetic comparison of two techniques for assessing molar crown size in pedigreed baboons

Dental anthropologists and paleoanthropologists commonly use an estimated molar crown area (mesiodistal length multiplied by buccolingual width) to describe and compare individuals, populations, and species. Advances in digital imaging now allow researchers to measure the actual crown area of a molar in an occlusal two-dimensional plane. Because error is reduced by this more accurate measurement, actual crown area is thought to be a better representation of the mechanisms that determine tooth crown size, meriting the additional time required to collect it. We tested this assumption by estimating the heritability of both these measurements for the second left mandibular molar from a sample of individuals (n = 332) from a captive breeding colony of baboons. Heritability estimates of both the actual and estimated crown areas of molars are approximately 0.83. Therefore, both measurements are informative as population descriptors, with no significant difference between the accuracy of either to reflect additive genetic contributions to molar crown size. This is fortunate, because genetic studies and inference can be based on estimated areas rather than actual crown area. The heritability estimates for mesiodistal length and buccolingual width are both substantial but lower: approximately 0.67 and approximately 0.73, respectively. The best fitting models in these analyses show that sex, body size, and subspecific affinity differentially affect molar length and width. We interpret these results to suggest that potentially some of the genetics underlying these covariates also underlie tooth size. As such, measurements designed to describe molar crown size are useful for general descriptive purposes, but do not conform to the assumption of independence inherent in phylogenetic analyses, such as cladistics (Hennig [1966] Phylogenetic Systematics. Urbana: University of Illinois Press). Therefore, if variables like actual crown area and estimated crown area are to be used in phylogenetic parsimony analyses, we suggest that researchers account for the effects of covariates such as sex and body size in their analyses.

Construction and analysis of a human-chimpanzee comparative clone map

The recently released human genome sequences provide us with reference data to conduct comparative genomic research on primates, which will be important to understand what genetic information makes us human. Here we present a first-generation human-chimpanzee comparative genome map and its initial analysis. The map was constructed through paired alignment of 77,461 chimpanzee bacterial artificial chromosome end sequences with publicly available human genome sequences. We detected candidate positions, including two clusters on human chromosome 21 that suggest large, nonrandom regions of difference between the two genomes.

Prediction of unidentified human genes on the basis of sequence similarity to novel cDNAs from cynomolgus monkey brain

BACKGROUND

The complete assignment of the protein-coding regions of the human genome is a major challenge for genome biology today. We have already isolated many hitherto unknown full-length cDNAs as orthologs of unidentified human genes from cDNA libraries of the cynomolgus monkey (Macaca fascicularis) brain (parietal lobe and cerebellum). In this study, we used cDNA libraries of three other parts of the brain (frontal lobe, temporal lobe and medulla oblongata) to isolate novel full-length cDNAs.

RESULTS

The entire sequences of novel cDNAs of the cynomolgus monkey were determined, and the orthologous human cDNA sequences were predicted from the human genome sequence. We predicted 29 novel human genes with putative coding regions sharing an open reading frame with the cynomolgus monkey, and we confirmed the expression of 21 pairs of genes by the reverse transcription-coupled polymerase chain reaction method. The hypothetical proteins were also functionally annotated by computer analysis.

CONCLUSIONS

The 29 new genes had not been discovered in recent explorations for novel genes in humans, and the ab initio method failed to predict all exons. Thus, monkey cDNA is a valuable resource for the preparation of a complete human gene catalog, which will facilitate post-genomic studies.

Using secondary structure to identify ribosomal numts: cautionary examples from the human genome

The identification of inadvertently sequenced mitochondrial pseudogenes (numts) is critical to any study employing mitochondrial DNA sequence data. Failure to discriminate numts correctly can confound phylogenetic reconstruction and studies of molecular evolution. This is especially problematic for ribosomal mtDNA genes. Unlike protein-coding loci, whose pseudogenes tend to accumulate diagnostic frameshift or premature stop mutations, functional ribosomal genes are not constrained to maintain a reading frame and can accumulate insertion-deletion events of varying length, particularly in nonpairing regions. Several authors have advocated using structural features of the transcribed rRNA molecule to differentiate functional mitochondrial rRNA genes from their nuclear paralogs. We explored this approach using the mitochondrial 12S rRNA gene and three known 12S numts from the human genome in the context of anthropoid phylogeny and the inferred secondary structure of primate 12S rRNA. Contrary to expectation, each of the three human numts exhibits striking concordance with secondary structure models, with little, if any, indication of their pseudogene status, and would likely escape detection based on structural criteria alone. Furthermore, we show that the unwitting inclusion of a particularly ancient (18-25 Myr old) and surprisingly cryptic human numt in a phylogenetic analysis would yield a well-supported but dramatically incorrect conclusion regarding anthropoid relationships. Though we endorse the use of secondary structure models for inferring positional homology wholeheartedly, we caution against reliance on structural criteria for the discrimination of rRNA numts, given the potential fallibility of this approach.

Astragalar morphology of late Eocene anthropoids from the Fayum Depression (Egypt) and the origin of catarrhine primates

The phylogenetic relationships of the late Eocene anthropoids Catopithecus browni and Proteopithecus sylviae are currently a matter of debate, with opinion divided as to whether these taxa are stem or crown anthropoids. The phylogenetic position of Catopithecus is of particular interest, for, unlike the highly generalized genus Proteopithecus, this taxon shares apomorphic dental and postcranial features with more derived undoubted catarrhines that appear in the same region 1-2 Ma later. If these apomorphies are homologous and Catopithecus is a stem catarrhine, the unique combination of plesiomorphic and apomorphic features preserved in this anthropoid would have important implications for our understanding of the crown anthropoid morphotype and the pattern of morphological character transformations that occurred during the early phases of stem catarrhine evolution.Well-preserved astragali referrable to Proteopithecus, Catopithecus, and the undoubted early Oligocene stem catarrhine Aegyptopithecus have provided additional morphological evidence that allows us to further evaluate competing hypotheses of interrelationships among Eocene-Oligocene Afro-Arabian anthropoids. Qualitative observations and multivariate morphometric analyses reveal that the astragalar morphology of Proteopithecus is very similar to that of early Oligocene parapithecids and living and extinct small-bodied platyrrhines, and strengthens the hypothesis that the morphological pattern shared by these taxa is primitive within crown Anthropoidea. In contrast, Catopithecus departs markedly from the predicted crown anthropoid astragalar morphotype and shares a number of apomorphic features (e.g., deep cotylar fossa, laterally projecting fibular facet, trochlear asymmetry, mediolaterally wide astragalar head) with Aegyptopithecus and Miocene-Recent catarrhines. The evidence from the astragalus complements other independent data from the dentition, humerus and femur of Catopithecus that support this taxon's stem catarrhine status, and we continue to maintain that oligopithecines are stem catarrhines that constitute the sister group of a clade containing propliopithecines and Miocene-Recent catarrhines.

The influence of selection on the evolutionary distance estimated from the base changes observed between homologous nucleotide sequences

In most studies of molecular evolution, the nucleotide base at a site is assumed to change with the apparent rate under functional constraint, and the comparison of base changes between homologous genes is thought to yield the evolutionary distance corresponding to the site-average change rate multiplied by the divergence time. However, this view is not sufficiently successful in estimating the divergence time of species, but mostly results in the construction of tree topology without a time-scale. In the present paper, this problem is investigated theoretically by considering that observed base changes are the results of comparing the survivals through selection of mutated bases. In the case of weak selection, the time course of base changes due to mutation and selection can be obtained analytically, leading to a theoretical equation showing how the selection has influence on the evolutionary distance estimated from the enumeration of base changes. This result provides a new method for estimating the divergence time more accurately from the observed base changes by evaluating both the strength of selection and the mutation rate. The validity of this method is verified by analysing the base changes observed at the third codon positions of amino acid residues with four-fold codon degeneracy in the protein genes of mammalian mitochondria; i.e. the ratios of estimated divergence times are fairly well consistent with a series of fossil records of mammals. Throughout this analysis, it is also suggested that the mutation rates in mitochondrial genomes are almost the same in different lineages of mammals and that the lineage-specific base-change rates indicated previously are due to the selection probably arising from the preference of transfer RNAs to codons.

Dynamics of microsatellite divergence under stepwise mutation and proportional slippage/point mutation models

Recently Kruglyak, Durrett, Schug, and Aquadro showed that microsatellite equilibrium distributions can result from a balance between polymerase slippage and point mutations. Here, we introduce an elaboration of their model that keeps track of all parts of a perfect repeat and a simplification that ignores point mutations. We develop a detailed mathematical theory for these models that exhibits properties of microsatellite distributions, such as positive skewness of allele lengths, that are consistent with data but are inconsistent with the predictions of the stepwise mutation model. We use our theoretical results to analyze the successes and failures of the genetic distances (delta(mu))(2) and D(SW) when used to date four divergences: African vs. non-African human populations, humans vs. chimpanzees, Drosophila melanogaster vs. D. simulans, and sheep vs. cattle. The influence of point mutations explains some of the problems with the last two examples, as does the fact that these genetic distances have large stochastic variance. However, we find that these two features are not enough to explain the problems of dating the human-chimpanzee split. One possible explanation of this phenomenon is that long microsatellites have a mutational bias that favors contractions over expansions.

Drosophila euchromatic LTR retrotransposons are much younger than the host species in which they reside

The recent release of the complete euchromatic genome sequence of Drosophila melanogaster offers a unique opportunity to explore the evolutionary history of transposable elements (TEs) within the genome of a higher eukaryote. In this report, we describe the annotation and phylogenetic comparison of 178 full-length long terminal repeat (LTR) retrotransposons from the sequenced component of the D. melanogaster genome. We report the characterization of 17 LTR retrotransposon families described previously and five newly discovered element families. Phylogenetically, these families can be divided into three distinct lineages that consist of members from the canonical Copia and Gypsy groups as well as a newly discovered third group containing BEL, mazi, and roo elements. Each family consists of members with average pairwise identities > or =99% at the nucleotide level, indicating they may be the products of recent transposition events. Consistent with the recent transposition hypothesis, we found that 70% (125/178) of the elements (across all families) have identical intra-element LTRs. Using the synonymous substitution rate that has been calculated previously for Drosophila (.016 substitutions per site per million years) and the intra-element LTR divergence calculated here, the average age of the remaining 30% (53/178) of the elements was found to be 137,000 +/-89,000 yr. Collectively, these results indicate that many full-length LTR retrotransposons present in the D. melanogaster genome have transposed well after this species diverged from its closest relative Drosophila simulans, 2.3 +/-.3 million years ago.

Evolution of natural killer cell receptors: coexistence of functional Ly49 and KIR genes in baboons

Natural killer (NK) cells represent an important first line of defense against viruses and malignancy [1]. NK cells express a variety of inhibitory and activating receptors that interact with classical major histocompatibility complex (MHC) class I molecules on potential target cells and determine the NK cell response [2-4]. Mouse NK receptors are encoded by the C-type lectin multigene family Ly49. However, in humans, a completely different family of receptors, the immunoglobulin-like killer inhibitory receptors (KIRs), performs the same function [2-4]. One Ly49-like gene, Ly49L, exists in humans but is incorrectly spliced and assumed to be nonfunctional [5, 6]. Mouse KIR-like genes have not been found, and evidence suggests that the primate KIRs amplified after rodents and primates diverged [7, 8]. Thus, two structurally dissimilar families, Ly49 and KIR, have evolved to play similar roles in mouse and human NK cells. This apparent example of functional convergent evolution raises several questions. It is unknown, for example, when the Ly49L gene became nonfunctional and if this event affected the functional evolution of the KIRs. The distribution of these gene families in different mammals is unstudied, and it is not known if any species uses both types of receptors. Here, we demonstrate that the Ly49L gene shows evidence of conservation in other mammals and that the human gene likely became nonfunctional 6-10 million years ago. Furthermore, we show that baboon lymphocytes express both full-length Ly49L transcripts and multiple KIR genes.

Excess of rare amino acid polymorphisms in the Toll-like receptor 4 in humans

The Toll-like receptor 4 protein acts as the transducing subunit of the lipopolysaccharide receptor complex and assists in the detection of Gram-negative pathogens within the mammalian host. Several lines of evidence support the view that variation at the TLR4 locus may alter host susceptibility to Gram-negative infection or the outcome of infection. Here, we surveyed TLR4 sequence variation in the complete coding region (2.4 kb) in 348 individuals from several population samples; in addition, a subset of the individuals was surveyed at 1.1 kb of intronic sequence. More than 90% of the chromosomes examined encoded the same structural isoform of TLR4, while the rest harbored 12 rare amino acid variants. Conversely, the variants at silent sites (intronic and synonymous positions) occur at both low and high frequencies and are consistent with a neutral model of mutation and random drift. The spectrum of allele frequencies for amino acid variants shows a significant skew toward lower frequencies relative to both the neutral model and the pattern observed at linked silent sites. This is consistent with the hypothesis that weak purifying selection acted on TLR4 and that most mutations affecting TLR4 protein structure have at least mildly deleterious phenotypic effects. These results may imply that genetic variants contributing to disease susceptibility occur at low frequencies in the population and suggest strategies for optimizing the design of disease-mapping studies.

Molecular estimates of primate divergences and new hypotheses for primate dispersal and the origin of modern humans

The concept of recent hominoid divergences has been a mainstay in molecular primatology since the 1970's. However, the ages allocated to the calibration points used to establish these divergence times and the estimates resulting from their application, notably the commonly accepted divergence between Pan (chimpanzees) and Homo 5 million years before present (MYBP), are now palaeontologically refutable. Here we estimate the ages of various primate divergences using three references with a more detailed fossil record than any of the traditional primate calibration points. Our findings suggest that the latter yield datings that are too recent by a factor of about two. For example, our estimates place the divergence between Pan and Homo 10.5-13 MYBP. The revised estimates of primate divergence times suggest a new hypothesis for primate evolution and dispersal: that the divergence between strepsirhines (lorises, lemurs) and anthropoids was contemporary with the break-up of Southern continents about 90 MYBP, with strepsirhines becoming isolated on Madagascar and later dispersing to Africa (and Asia) and anthropoids evolving in South America and subsequently colonizing Africa (and Asia), or possibly North America. In addition we present a new hypothesis, which accommodates the strikingly similar coalescence times for human mitochondrial DNA and the Y-chromosome. This hypothesis posits a common mitochondrial and Y-chromosome bottleneck about 400,000 years ago, associated with the origination of the human 2n = 46 karyotype, obstructing genetic exchange with the 2n = 48 Homo contemporaries.

A HERV-K provirus in chimpanzees, bonobos and gorillas, but not humans

Evidence from DNA sequencing studies strongly indicated that humans and chimpanzees are more closely related to each other than either is to gorillas [1-4]. However, precise details of the nature of the evolutionary separation of the lineage leading to humans from those leading to the African great apes have remained uncertain. The unique insertion sites of endogenous retroviruses, like those of other transposable genetic elements, should be useful for resolving phylogenetic relationships among closely related species. We identified a human endogenous retrovirus K (HERV-K) provirus that is present at the orthologous position in the gorilla and chimpanzee genomes, but not in the human genome. Humans contain an intact preintegration site at this locus. These observations provide very strong evidence that, for some fraction of the genome, chimpanzees, bonobos, and gorillas are more closely related to each other than they are to humans. They also show that HERV-K replicated as a virus and reinfected the germline of the common ancestor of the four modern species during the period of time when the lineages were separating and demonstrate the utility of using HERV-K to trace human evolution.

Human DNA sequence variation in a 6.6-kb region containing the melanocortin 1 receptor promoter

An approximately 6.6-kb region located upstream from the melanocortin 1 receptor (MC1R) gene and containing its promoter was sequenced in 54 humans (18 Africans, 18 Asians, and 18 Europeans) and in one chimpanzee, gorilla, and orangutan. Seventy-six polymorphic sites were found among the human sequences and the average nucleotide diversity (pi) was 0.141%, one of the highest among all studies of nuclear sequence variation in humans. Opposite to the pattern observed in the MC1R coding region, in the present region pi is highest in Africans (0.136%) compared to Asians (0.116%) and Europeans (0.122%). The distributions of pi, theta, and Fu and Li's F-statistic are nonuniform along the sequence and among continents. The pattern of genetic variation is consistent with a population expansion in Africans. We also suggest a possible phase of population size reduction in non-Africans and purifying selection acting in the middle subregion and parts of the 5' subregion in Africans. We hypothesize diversifying selection acting on some sites in the 5' and 3' subregions or in the MC1R coding region in Asians and Europeans, though we cannot reject the possibility of relaxation of functional constraints in the MC1R gene in Asians and Europeans. The mutation rate in the sequenced region is 1.65 x 10(-9) per site per year. The age of the most recent common ancestor for this region is similar to that for the other long noncoding regions studied to date, providing evidence for ancient gene genealogies. Our population screening and phylogenetic footprinting suggest potentially important sites for the MC1R promoter function.

The evolution of body armor in mammals: plantigrade constraints of large body size

Predictions associated with opposing selection generating minimum variance in basal metabolic rate (BMR) in mammals at a constrained body mass (CBM; 358 g) were tested. The CBM is presumed to be associated with energetic constraints linked to predation and variable resources at intermediate sizes on a logarithmic mass scale. Opposing selection is thought to occur in response to energetic constraints associated with predation and unpredictable resources. As body size approaches and exceeds the CBM, mammals face increasing risks of predation and daily energy requirements. Fast running speeds may require high BMRs, but unpredictable and low resources may select for low BMRs, which also reduce foraging time and distances and thus predation risks. If these two selection forces oppose each other persistently, minimum BMR variance may result. However, extreme BMR outliers at and close to the CBM should be indicative of unbalanced selection and predator avoidance alternatives (escapers vs. defenders), and may therefore provide indirect support for opposing selection. It was confirmed that body armor in defenders evolves at and above the CBM, and armored mammals had significantly lower BMRs than their nonarmored counterparts. However, analyses comparing the BMR of escapers--the fastest nonarmored runners (Lagomorpha)--with similar-sized counterparts were inconclusive and were confounded by limb morphology associated with speed optimization. These analyses suggest that the risks and costs of predation and the speed limitations of the plantigrade foot may constrain the evolution of large body sizes in plantigrade mammals.

Birth of a gene: locus of neuronal BC200 snmRNA in three prosimians and human BC200 pseudogenes as archives of change in the Anthropoidea lineage

The gene encoding brain-specific dendritic BC200 small non-messenger RNA is limited to the primate order and arose from a monomeric Alu element. It is present and neuronally expressed in all Anthropoidea examined. By comparing the human sequence of about 13.2 kb with each of the prosimian (lemur 14.6 kb, galago 12 kb, and tarsier 13.8 kb) orthologous loci, we could establish that the BC200 RNA gene is absent from the prosimian lineages. In Strepsirhini (lemurs and lorises), a dimeric AluJ-like element integrated very close to the BC200 insertion point, while the corresponding tarsier region is devoid of any repetitive element. Consequently, insertion of the Alu monomer that gave rise to the BC200 RNA gene must have occurred after the anthropoid lineage diverged from the prosimian lineage(s). Shared insertions of other repetitive elements favor proximity of simians and tarsiers in support of their grouping into Haplorhini and the omomyid hypothesis. On the other hand, the nucleotide sequences in the segment that is available for comparison in all four species reveal less exchanges between Strepsirhini (lemur and galago) and human than between tarsier and human. Our data imply that the early activity of dimeric Alu sequences must have been concurrent with the activity of monomeric Alu elements that persisted longer than is usually thought. As BC200 RNA gave rise to more than 200 pseudogenes, we used their consensus sequence variations as a molecular archive recording the BC200 RNA sequence changes in the anthropoid lineage leading to Homo sapiens and timed these alterations over the past 35-55 million years.

Episodic Evolution of Growth Hormone in Primates and Emergence of the Species Specificity of Human Growth Hormone Receptor

Growth hormone (GH) evolution is very conservative among mammals, except for primates and ruminant artiodactyls. In fact, most known mammalian GH sequences differ from the inferred ancestral mammalian sequence by only a few amino acids. In contrast, the human GH sequence differs from the inferred ancestral sequence by 59 amino acids. However, it is not known when this rapid evolution of GH occurred during primate evolution or whether it was due to positive selection. Also, human growth hormone receptor (GHR) displays species specificity; i.e., it can interact only with human (or rhesus monkey) GH, not with nonprimate GHS: The species specificity of human GHR is largely due to the Leu-->Arg change at position 43, and it has been hypothesized that this change must have been preceded by the His-->Asp change at position 171 of GH. Is this hypothesis true? And when did these changes occur? To address the above issues, we sequenced GH and GHR genes in prosimians and simians. Our data supported the above hypothesis and revealed that the species specificity of human GHR actually emerged in the common ancestor of Old World primates, but the transitional phase still persists in New World monkeys. Our data showed that the rapid evolution of primate GH occurred during a relatively short period (in the common ancestor of higher primates) and that the rate of change was especially high at functionally important sites, suggesting positive selection. However, the nonsynonymous rate/synonymous rate ratio at these sites was <1, so relaxation of purifying selection might have played a role in the rapid evolution of the GH gene in simians, possibly as a result of multiple gene duplications. Similar to GH, GHR displayed an accelerated rate of evolution in primates. Our data revealed proportionally more amino acid replacements at the functionally important sites in both GH and GHR in simians but, surprisingly, showed few coincidental replacements of amino acids forming the same intermolecular contacts between the two proteins.

Molecular phylogeny of the major hylobatid divisions

We describe DNA sequences for the mitochondrial control region and phenylalanine-tRNA from the four extant gibbon subgenera. In contrast to earlier studies on gibbon phylogeny that used other parts of the mtDNA, the control region depicts the crested gibbons (Nomascus) as the most basal group of the Hylobatidae, followed by Symphalangus, with Bunopithecus and Hylobates as the last to diverge. Our data show that the molecular distances among the four gibbon subgenera are in the same range as those between Homo and Pan, or even higher. As a consequence of these findings, we propose to raise all four gibbon subgenera to genus rank.

TTY2: a multicopy Y-linked gene family

Genes involved in human male sex determination and spermatogenesis are likely to be located on the Y chromosome. In an effort to identify Y-linked, testis-expressed genes, a cDNA selection library was generated by selecting testis cDNA with Y-cosmid clones. Resultant clones containing repetitive or vector material were eliminated, and 79 of the remaining clones were sequenced. Nineteen cDNAs showed homology with the TTY2 gene, and indicated that TTY2 is part of a large gene family. Screening of a panel of Y-linked cosmids revealed that the TTY2 gene family includes at least 26 members organized in 14 subfamilies. Further investigation revealed that TTY2 genes are arranged in tandemly arrayed clusters on both arms of the Y chromosome, and each gene comprises a series of tandemly arranged repeats. RT-PCR studies for two of these genes revealed that they are expressed in adult and fetal testis, as well as in the adult kidney. None of the genes investigated in detail contain an open reading frame. We conclude that the TTY2 gene family is composed of multiple copies, some of which may function as noncoding RNA transcripts and some may be pseudogenes.

Proteomic comparison of human and great ape blood plasma reveals conserved glycosylation and differences in thyroid hormone metabolism

Most blood plasma proteins are glycosylated. These glycoproteins typically carry sialic acid-bearing sugar chains, which can modify the observed molecular weights and isoelectric points of those proteins during electrophoretic analyses. To explore changes in protein expression and glycosylation that occurred during great ape and human evolution, we subjected multiple blood plasma samples from all these species to high-resolution proteomic analysis. We found very few species-specific differences, indicating a remarkable degree of conservation of plasma protein expression and glycosylation during approximately 12 million years of evolution. A few lineage-specific differences in protein migration were noted among the great apes. The only obvious differences between humans and all great apes were an apparent decrease in transthyretin (prealbumin) and a change in haptoglobin isoforms (the latter was predictable from prior genetic studies). Quantitative studies of transthyretin in samples of blood plasma (synthesized primarily by the liver) and of cerebrospinal fluid (synthesized locally by the choroid plexus of the brain) confirmed approximately 2-fold higher levels in chimpanzees compared to humans. Since transthyretin binds thyroid hormones, we next compared plasma thyroid hormone parameters between humans and chimpanzees. The results indicate significant differences in the status of thyroid hormone metabolism, which represent the first known endocrine difference between these species. Notably, thyroid hormones are known to play major roles in the development, differentiation, and metabolism of many organs and tissues, including the brain and the cranium. Also, transthyretin is known to be the major carrier of thyroid hormone in the cerebrospinal fluid, likely regulating delivery of this hormone to the brain. A potential secondary difference in retinoid (vitamin A) metabolism is also noted. The implications of these findings for explaining unique features of human evolution are discussed.

Cytochrome b phylogeny and the taxonomy of great apes and mammals

In the Linnaean system of classification, the generic status of a species is part of its binomial name, and it is therefore important that the classification at the level of genus is consistent at least in related groups of organisms. Using maximum-likelihood phylogenetic trees constructed from a large number of complete or nearly complete mammalian cytochrome b sequences, I show that the distributions of intrageneric and intergeneric distances derived from these trees are clearly separated, which allows the limits for a more rational generic classification of mammals to be established. The analysis of genetic distances among hominids in this context provides strong support for the inclusion of humans and chimpanzees in the same genus. It is also of interest to decipher the main reasons for the possible biases in the mammalian classification. I found by correlation analysis that the classification of mammals of large body size tends to be oversplit, whereas that of small mammals has an excess of lumping, which may be a manifestation of the larger difficulty in finding diagnostic characters in the classification of small animals. In addition, and contrary to some previous observations, there is no correlation between body size and rate of cytochrome b evolution in mammals, which excludes the difference in evolutionary rates as the cause of the observed body size taxonomic bias.

Episodic evolution of pyrin in primates: human mutations recapitulate ancestral amino acid states

Familial Mediterranean fever (FMF; MIM 249100) is an autosomal recessive disease characterized by recurrent attacks of fever with synovial, pleural or peritoneal inflammation. The disease is caused by mutations in the gene encoding the pyrin protein. Human population studies have revealed extremely high allele frequencies for several different pyrin mutations, leading to the conclusion that the mutant alleles confer a selective advantage. Here we examine the ret finger protein (rfp) domain (which contains most of the disease-causing mutations) of pyrin during primate evolution. Amino acids that cause human disease are often present as wild type in other species. This is true at positions 653 (a novel mutation), 680, 681, 726, 744 and 761. For several of these human mutations, the mutant represents the reappearance of an ancestral amino acid state. Examination of lineage-specific dN/dS ratios revealed a pattern consistent with the signature of episodic positive selection. Our data, together with previous human population studies, indicate that selective pressures may have caused functional evolution of pyrin in humans and other primates.