Sooty mangabey genome sequence provides insight into AIDS resistance in a natural SIV host.

In contrast to infections with human immunodeficiency virus (HIV) in humans and simian immunodeficiency virus (SIV) in macaques, SIV infection of a natural host, sooty mangabeys (Cercocebus atys), is non-pathogenic despite high viraemia. Here we sequenced and assembled the genome of a captive sooty mangabey. We conducted genome-wide comparative analyses of transcript assemblies from C. atys and AIDS-susceptible species, such as humans and macaques, to identify candidates for host genetic factors that influence susceptibility. We identified several immune-related genes in the genome of C. atys that show substantial sequence divergence from macaques or humans. One of these sequence divergences, a C-terminal frameshift in the toll-like receptor-4 (TLR4) gene of C. atys, is associated with a blunted in vitro response to TLR-4 ligands. In addition, we found a major structural change in exons 3-4 of the immune-regulatory protein intercellular adhesion molecule 2 (ICAM-2); expression of this variant leads to reduced cell surface expression of ICAM-2. These data provide a resource for comparative genomic studies of HIV and/or SIV pathogenesis and may help to elucidate the mechanisms by which SIV-infected sooty mangabeys avoid AIDS.

Evidence from Cameroon reveals differences in the genetic structure and histories of chimpanzee populations.

The history of the genus Pan is a topic of enduring interest. Chimpanzees (Pan troglodytes) are often divided into subspecies, but the population structure and genetic history of chimpanzees across Africa remain unclear. Some population genetics studies have led to speculation that, until recently, this species constituted a single population with ongoing gene flow across its range, which resulted in a continuous gradient of allele frequencies. Chimpanzees, designated here as P. t. ellioti, occupy the Gulf of Guinea region that spans southern Nigeria and western Cameroon at the center of the distribution of this species. Remarkably, few studies have included individuals from this region, hindering the examination of chimpanzee population structure across Africa. Here, we analyzed microsatellite genotypes of 94 chimpanzees, including 32 designated as P. t. ellioti. We find that chimpanzees fall into three major populations: (i) Upper Guinea in western Africa (P. t. verus); (ii) the Gulf of Guinea region (P. t. ellioti); and (iii) equatorial Africa (P. t. troglodytes and P. t. schweinfurthii). Importantly, the Gulf of Guinea population is significantly different genetically from the others, sharing a last common ancestor with the populations in Upper Guinea ~0.46 million years ago (mya) and equatorial Africa ~0.32 mya. Equatorial chimpanzees are subdivided into up to three populations occupying southern Cameroon, central Africa, and eastern Africa, which may have constituted a single population until ~0.10-0.11 mya. Finally, occasional hybridization may be occurring between the Gulf of Guinea and southern Cameroon populations.

Successive radiations, not stasis, in the South American primate fauna.

The earliest Neotropical primate fossils complete enough for taxonomic assessment, Dolichocebus, Tremacebus, and Chilecebus, date to approximately 20 Ma. These have been interpreted as either closely related to extant forms or as extinct stem lineages. The former hypothesis of morphological stasis requires most living platyrrhine genera to have diverged before 20 Ma. To test this hypothesis, we collected new complete mitochondrial genomes from Aotus lemurinus, Saimiri sciureus, Saguinus oedipus, Ateles belzebuth, and Callicebus donacophilus. We combined these with published sequences from Cebus albifrons and other primates to infer the mitochondrial phylogeny. We found support for a cebid/atelid clade to the exclusion of the pitheciids. Then, using Bayesian methods and well-supported fossil calibration constraints, we estimated that the platyrrhine most recent common ancestor (MRCA) dates to 19.5 Ma, with all major lineages diverging by 14.3 Ma. Next, we estimated catarrhine divergence dates on the basis of platyrrhine divergence scenarios and found that only a platyrrhine MRCA less than 21 Ma is concordant with the catarrhine fossil record. Finally, we calculated that 33% more change in the rate of evolution is required for platyrrhine divergences consistent with the morphologic stasis hypothesis than for a more recent radiation. We conclude that Dolichocebus, Tremacebus, and Chilecebus are likely too old to be crown platyrrhines, suggesting they were part of an extinct early radiation. We note that the crown platyrrhine radiation was concomitant with the radiation of 2 South American xenarthran lineages and follows a global temperature peak and tectonic activity in the Andes.

Evolution of the mammalian lysozyme gene family.

BACKGROUND: Lysozyme c (chicken-type lysozyme) has an important role in host defense, and has been extensively studied as a model in molecular biology, enzymology, protein chemistry, and crystallography. Traditionally, lysozyme c has been considered to be part of a small family that includes genes for two other proteins, lactalbumin, which is found only in mammals, and calcium-binding lysozyme, which is found in only a few species of birds and mammals. More recently, additional testes-expressed members of this family have been identified in human and mouse, suggesting that the mammalian lysozyme gene family is larger than previously known. RESULTS: Here we characterize the extent and diversity of the lysozyme gene family in the genomes of phylogenetically diverse mammals, and show that this family contains at least eight different genes that likely duplicated prior to the diversification of extant mammals. These duplicated genes have largely been maintained, both in intron-exon structure and in genomic context, throughout mammalian evolution. CONCLUSIONS: The mammalian lysozyme gene family is much larger than previously appreciated and consists of at least eight distinct genes scattered around the genome. Since the lysozyme c and lactalbumin proteins have acquired very different functions during evolution, it is likely that many of the other members of the lysozyme-like family will also have diverse and unexpected biological properties.

Mitochondrial and nuclear markers suggest Hanuman langur (Primates: Colobinae) polyphyly: implications for their species status.

Recent molecular studies on langurs of the Indian subcontinent suggest that the widely-distributed and morphologically variable Hanuman langurs (Semnopithecus entellus) are polyphyletic with respect to Nilgiri and purple-faced langurs. To further investigate this scenario, we have analyzed additional sequences of mitochondrial cytochrome b as well as nuclear protamine P1 genes from these species. The results confirm Hanuman langur polyphyly in the mitochondrial tree and the nuclear markers suggest that the Hanuman langurs share protamine P1 alleles with Nilgiri and purple-faced langurs. We recommend provisional splitting of the so-called Hanuman langurs into three species such that the taxonomy is consistent with their evolutionary relationships.

Molecular phylogeny and biogeography of langurs and leaf monkeys of South Asia (Primates: Colobinae).

The two recently proposed taxonomies of the langurs and leaf monkeys (Subfamily Colobinae) provide different implications to our understanding of the evolution of Nilgiri and purple-faced langurs. Groves (2001) [Groves, C.P., 2001. Primate Taxonomy. Smithsonian Institute Press, Washington], placed Nilgiri and purple-faced langurs in the genus Trachypithecus, thereby suggesting disjunct distribution of the genus Trachypithecus. [Brandon-Jones, D., Eudey, A.A., Geissmann, T., Groves, C.P., Melnick, D.J., Morales, J.C., Shekelle, M., Stewart, C.-B., 2003. Asian primate classification. Int. J. Primatol. 25, 97-162] placed these langurs in the genus Semnopithecus, which suggests convergence of morphological characters in Nilgiri and purple-faced langurs with Trachypithecus. To test these scenarios, we sequenced and analyzed the mitochondrial cytochrome b gene and two nuclear DNA-encoded genes, lysozyme and protamine P1, from a variety of colobine species. All three markers support the clustering of Nilgiri and purple-faced langurs with Hanuman langur (Semnopithecus), while leaf monkeys of Southeast Asian (Trachypithecus) form a distinct clade. The phylogenetic position of capped and golden leaf monkeys is still unresolved. It is likely that this species group might have evolved due to past hybridization between Semnopithecus and Trachypithecus clades.

Mitochondrial data support an odd-nosed colobine clade.

To obtain a more complete understanding of the evolutionary history of the leaf-eating monkeys we have examined the mitochondrial genome sequence of two African and six Asian colobines. Although taxonomists have proposed grouping the "odd-nosed" colobines (proboscis monkey, douc langur, and the snub-nosed monkey) together, phylogenetic support for such a clade has not been tested using molecular data. Phylogenetic analyses using parsimony, maximum likelihood, and Bayesian methods support a monophyletic clade of odd-nosed colobines consisting of Nasalis, Pygathrix, and Rhinopithecus, with tentative support for Nasalis occupying a basal position within this clade. The African and Asian colobine lineages are inferred to have diverged by 10.8 million years ago (mya or Ma). Within the Asian colobines the odd-nosed clade began to diversify by 6.7 Ma. These results augment our understanding of colobine evolution, particularly the nature and timing of the colobine expansion into Asia. This phylogenetic information will aid those developing conservation strategies for these highly endangered, diverse, and unique primates.

Duplication and divergence of 2 distinct pancreatic ribonuclease genes in leaf-eating African and Asian colobine monkeys.

Unique among primates, the colobine monkeys have adapted to a predominantly leaf-eating diet by evolving a foregut that utilizes bacterial fermentation to breakdown and absorb nutrients from such a food source. It has been hypothesized that pancreatic ribonuclease (pRNase) has been recruited to perform a role as a digestive enzyme in foregut fermenters, such as artiodactyl ruminants and the colobines. We present molecular analyses of 23 pRNase gene sequences generated from 8 primate taxa, including 2 African and 2 Asian colobine species. The pRNase gene is single copy in all noncolobine primate species assayed but has duplicated more than once in both the African and Asian colobine monkeys. Phylogenetic reconstructions show that the pRNase-coding and noncoding regions are under different evolutionary constraints, with high levels of concerted evolution among gene duplicates occurring predominantly in the noncoding regions. Our data suggest that 2 functionally distinct pRNases have been selected for in the colobine monkeys, with one group adapting to the role of a digestive enzyme by evolving at an increased rate with loss of positive charge, namely arginine residues. Conclusions relating our data to general hypotheses of evolution following gene duplication are discussed.

Catarrhine primate divergence dates estimated from complete mitochondrial genomes: concordance with fossil and nuclear DNA evidence.

Accurate divergence date estimates improve scenarios of primate evolutionary history and aid in interpretation of the natural history of disease-causing agents. While molecule-based estimates of divergence dates of taxa within the superfamily Hominoidea (apes and humans) are common in the literature, few such estimates are available for the Cercopithecoidea (Old World monkeys), the sister taxon of the hominoids in the primate infraorder Catarrhini. To help fill this gap, we have sequenced the entire mitochondrial DNA (mtDNA) genomes from a representative of three cercopithecoid tribes, Cercopithecini (Chlorocebus aethiops), Colobini (Colobus guereza), and Presbytini (Trachypithecus obscurus), and analyzed these new data together with other catarrhine mtDNA genomes available in public databases. Molecular divergence date estimates are dependent on calibration points gleaned from the paleontological record. We defined criteria for the selection of good calibration points and identified three points meeting these criteria: Homo-Pan, 6.0 Ma; Pongo-hominines, 14.0 Ma; hominoid/cercopithecoid, 23.0 Ma. Because a uniform molecular clock does not fit the catarrhine mtDNA data, we estimated divergence dates using a penalized likelihood and a Bayesian method, both of which take into account the effects of rate differences on lineages, phylogenetic tree structure, and multiple calibration points. The penalized likelihood method applied to the coding regions of the mtDNA genome yielded the following divergence date estimates, with approximate 95% confidence intervals: cercopithecine-colobine, 16.2 (14.4-17.9) Ma; colobin-presbytin, 10.9 (9.6-12.3) Ma; cercopithecin-papionin, 11.6 (10.3-12.9) Ma; and Macaca-Papio, 9.8 (8.6-10.9) Ma. Within the hominoids, the following dates were inferred: hylobatid-hominid, 16.8 (15.0-18.5) Ma; Gorilla-Homo+Pan, 8.1 (7.1-9.0) Ma; Pongo pygmaeus pygmaeus-P. p. abelii, 4.1 (3.5-4.7) Ma; and Pan troglodytes-P. paniscus, 2.4 (2.0-2.7) Ma. These dates were similar to those found using penalized likelihood on other subsets of the data, but slightly younger than several of the Bayesian estimates.

Evolution of base-substitution gradients in primate mitochondrial genomes.

Inferences of phylogenies and dates of divergence rely on accurate modeling of evolutionary processes; they may be confounded by variation in substitution rates among sites and changes in evolutionary processes over time. In vertebrate mitochondrial genomes, substitution rates are affected by a gradient along the genome of the time spent being single-stranded during replication, and different types of substitutions respond differently to this gradient. The gradient is controlled by biological factors including the rate of replication and functionality of repair mechanisms; little is known, however, about the consistency of the gradient over evolutionary time, or about how evolution of this gradient might affect phylogenetic analysis. Here, we evaluate the evolution of response to this gradient in complete primate mitochondrial genomes, focusing particularly on A-->G substitutions, which increase linearly with the gradient. We developed a methodology to evaluate the posterior probability densities of the response parameter space, and used likelihood ratio tests and mixture models with different numbers of classes to determine whether groups of genomes have evolved in a similar fashion. Substitution gradients usually evolve slowly in primates, but there have been at least two large evolutionary jumps: on the lineage leading to the great apes, and a convergent change on the lineage leading to baboons (Papio). There have also been possible convergences at deeper taxonomic levels, and different types of substitutions appear to evolve independently. The placements of the tarsier and the tree shrew within and in relation to primates may be incorrect because of convergence in these factors.

A mobile element based phylogeny of Old World monkeys.

SINEs (Short INterspersed Elements) are a class of non-autonomous mobile elements that are <500 bp in length and have no open reading frames. Individual SINE elements are essentially homoplasy free with known ancestral states, making them useful genetic systems for phylogenetic studies. Alu elements are the most successful SINE in primate genomes and have been utilized for resolving primate phylogenetic relationships and human population genetics. However, no Alu based phylogenetic analysis has yet been performed to resolve relationships among Old World monkeys. Using both a computational approach and polymerase chain reaction display methodology, we identified 285 new Alu insertions from sixteen Old World monkey taxa that were informative at various levels of catarrhine phylogeny. We have utilized these elements along with 12 previously reported loci to construct a phylogenetic tree of the selected taxa. Relationships among all major clades are in general agreement with other molecular and morphological data sets but have stronger statistical support.

Ancestral sequence reconstruction in primate mitochondrial DNA: compositional bias and effect on functional inference.

Reconstruction of ancestral DNA and amino acid sequences is an important means of inferring information about past evolutionary events. Such reconstructions suggest changes in molecular function and evolutionary processes over the course of evolution and are used to infer adaptation and convergence. Maximum likelihood (ML) is generally thought to provide relatively accurate reconstructed sequences compared to parsimony, but both methods lead to the inference of multiple directional changes in nucleotide frequencies in primate mitochondrial DNA (mtDNA). To better understand this surprising result, as well as to better understand how parsimony and ML differ, we constructed a series of computationally simple "conditional pathway" methods that differed in the number of substitutions allowed per site along each branch, and we also evaluated the entire Bayesian posterior frequency distribution of reconstructed ancestral states. We analyzed primate mitochondrial cytochrome b (Cyt-b) and cytochrome oxidase subunit I (COI) genes and found that ML reconstructs ancestral frequencies that are often more different from tip sequences than are parsimony reconstructions. In contrast, frequency reconstructions based on the posterior ensemble more closely resemble extant nucleotide frequencies. Simulations indicate that these differences in ancestral sequence inference are probably due to deterministic bias caused by high uncertainty in the optimization-based ancestral reconstruction methods (parsimony, ML, Bayesian maximum a posteriori). In contrast, ancestral nucleotide frequencies based on an average of the Bayesian set of credible ancestral sequences are much less biased. The methods involving simpler conditional pathway calculations have slightly reduced likelihood values compared to full likelihood calculations, but they can provide fairly unbiased nucleotide reconstructions and may be useful in more complex phylogenetic analyses than considered here due to their speed and flexibility. To determine whether biased reconstructions using optimization methods might affect inferences of functional properties, ancestral primate mitochondrial tRNA sequences were inferred and helix-forming propensities for conserved pairs were evaluated in silico. For ambiguously reconstructed nucleotides at sites with high base composition variability, ancestral tRNA sequences from Bayesian analyses were more compatible with canonical base pairing than were those inferred by other methods. Thus, nucleotide bias in reconstructed sequences apparently can lead to serious bias and inaccuracies in functional predictions.

Comprehensive analysis of two Alu Yd subfamilies.

Alu elements have inserted in the human genome throughout primate evolution. A small number of Alu insertions have occurred after the divergence of humans from nonhuman primates and therefore should not be present in nonhuman primate genomes. Most of these recently integrated Alu elements are contained with a series of discrete Alu subfamilies that are related to each other based upon diagnostic nucleotide substitutions. We have extracted members of the Alu Yd subfamily that are derivatives of the Alu Y subfamily that share a common 12-bp deletion that defines the Yd lineage from the draft sequence of the human genome. Analysis of the Yd Alu elements resulted in the recovery of two new Alu subfamilies, Yd3 and Yd6, which contain a total of 295 members (198 Yd3 and 97 Yd6). DNA sequence analysis of each of the Alu Yd subfamilies yielded age estimates of 8.02 and 1.20 million years old for the Alu Yd3 and Yd6 subfamilies, respectively. Two hundred Alu Yd3 and Yd6 loci were screened using polymerase chain reaction (PCR) assays to determine their phylogenetic origin and associated levels of human genomic diversity. The Alu Yd3 subfamily appears to have started amplifying relatively early in primate evolution and continued propagating albeit at a low level as many of its members are found in a variety of hominoid (humans, greater and lesser ape) genomes. Only two of the elements are polymorphic in the human genome and absent from the genomes of nonhuman primates. By contrast all of the members of the Alu Yd6 subfamily are restricted to the human genome, with 12% of the elements representing insertion polymorphisms in human populations. A single Alu Yd6 locus contained an independent parallel forward insertion of a paralogous Alu Sq sequence in the owl monkey. These Alu subfamilies are a source of genomic fossil relics for the study of primate phylogenetics and human population genetics.