Framework of the Alu Subfamily Evolution in the Platyrrhine Three-Family Clade of Cebidae, Callithrichidae, and Aotidae

The history of Alu retroposons has been choreographed by the systematic accumulation of inherited diagnostic nucleotide substitutions to form discrete subfamilies, each having a distinct nucleotide consensus sequence. The oldest subfamily, AluJ, gave rise to AluS after the split between Strepsirrhini and what would become Catarrhini and Platyrrhini. The AluS lineage gave rise to AluY in catarrhines and to AluTa in platyrrhines. Platyrrhine Alu subfamilies Ta7, Ta10, and Ta15 were assigned names based on a standardized nomenclature. However, with the subsequent intensification of whole genome sequencing (WGS), large scale analyses to characterize Alu subfamilies using the program COSEG identified entire lineages of subfamilies simultaneously. The first platyrrhine genome with WGS, the common marmoset (Callithrix jacchus; [caljac3]), resulted in Alu subfamily names sf0 to sf94 in an arbitrary order. Although easily resolved by alignment of the consensus sequences, this naming convention can become increasingly confusing as more genomes are independently analyzed. In this study, we reported Alu subfamily characterization for the platyrrhine three-family clade of Cebidae, Callithrichidae, and Aotidae. We investigated one species/genome from each recognized family of Callithrichidae and Aotidae and of both subfamilies (Cebinae and Saimiriinae) of the family Cebidae. Furthermore, we constructed a comprehensive network of Alu subfamily evolution within the three-family clade of platyrrhines to provide a working framework for future research. Alu expansion in the three-family clade has been dominated by AluTa15 and its derivatives.

Analysis of Simian Endogenous Retrovirus (SERV) Full-Length Proviruses in Old World Monkey Genomes

Simian endogenous retrovirus, SERV, is a successful germ line invader restricted to Old World monkey (OWM) species. (1) Background: The availability of high-quality primate genomes warrants a study of the characteristics, evolution, and distribution of SERV proviruses. (2) Methods: Cercopithecinae OWM genomes from public databases were queried for the presence of full-length SERV proviruses. A dataset of 81 Cer-SERV genomes was generated and analyzed. (3) Results: Full-length Cer-SERV proviruses were mainly found in terrestrial OWM, and less so in arboreal, forest- dwelling monkeys. Phylogenetic analysis confirmed the existence of two genotypes, Cer-SERV-1 and Cer-SERV-2, with Cer-SERV-1 showing evidence of recent germ-line expansions. Long Terminal Repeat (LTR) variation indicated that most proviruses were of a similar age and were estimated to be between <0.3 and 10 million years old. Integrations shared between species were relatively rare. Sequence analysis further showed extensive CpG methylation-associated mutations, variable Primer Binding Site (PBS) use with Cer-SERV-1 using PBSlys3 and Cer-SERV-2 using PBSlys1,2, and the recent gain of LTR motifs for transcription factors active during embryogenesis in Cer-SERV-1. (4) Conclusions: sequence analysis of 81 SERV proviruses from Cercopithecinae OWM genomes provides evidence for the adaptation of this retrovirus to germ line reproduction.

Comparative Genomics Analysis Reveals High Levels of Differential Retrotransposition among Primates from the Hominidae and the Cercopithecidae Families

Mobile elements (MEs), making ∼50% of primate genomes, are known to be responsible for generating inter- and intra-species genomic variations and play important roles in genome evolution and gene function. Using a bioinformatics comparative genomics approach, we performed analyses of species-specific MEs (SS-MEs) in eight primate genomes from the families of Hominidae and Cercopithecidae, focusing on retrotransposons. We identified a total of 230,855 SS-MEs, with which we performed normalization based on evolutionary distances, and we also analyzed the most recent SS-MEs in these genomes. Comparative analysis of SS-MEs reveals striking differences in ME transposition among these primate genomes. Interesting highlights of our results include: 1) the baboon genome has the highest number of SS-MEs with a strong bias for SINEs, while the crab-eating macaque genome has a sustained extremely low transposition for all ME classes, suggesting the existence of a genome-wide mechanism suppressing ME transposition; 2) while SS-SINEs represent the dominant class in general, the orangutan genome stands out by having SS-LINEs as the dominant class; 3) the human genome stands out among the eight genomes by having the largest number of recent highly active ME subfamilies, suggesting a greater impact of ME transposition on its recent evolution; and 4) at least 33% of the SS-MEs locate to genic regions, including protein coding regions, presenting significant potentials for impacting gene function. Our study, as the first of its kind, demonstrates that mobile elements evolve quite differently among these primates, suggesting differential ME transposition as an important mechanism in primate evolution.

Convergence of human and Old World monkey gut microbiomes demonstrates the importance of human ecology over phylogeny

BACKGROUND

Comparative data from non-human primates provide insight into the processes that shaped the evolution of the human gut microbiome and highlight microbiome traits that differentiate humans from other primates. Here, in an effort to improve our understanding of the human microbiome, we compare gut microbiome composition and functional potential in 14 populations of humans from ten nations and 18 species of wild, non-human primates.

RESULTS

Contrary to expectations from host phylogenetics, we find that human gut microbiome composition and functional potential are more similar to those of cercopithecines, a subfamily of Old World monkey, particularly baboons, than to those of African apes. Additionally, our data reveal more inter-individual variation in gut microbiome functional potential within the human species than across other primate species, suggesting that the human gut microbiome may exhibit more plasticity in response to environmental variation compared to that of other primates.

CONCLUSIONS

Given similarities of ancestral human habitats and dietary strategies to those of baboons, these findings suggest that convergent ecologies shaped the gut microbiomes of both humans and cercopithecines, perhaps through environmental exposure to microbes, diet, and/or associated physiological adaptations. Increased inter-individual variation in the human microbiome may be associated with human dietary diversity or the ability of humans to inhabit novel environments. Overall, these findings show that diet, ecology, and physiological adaptations are more important than host-microbe co-diversification in shaping the human microbiome, providing a key foundation for comparative analyses of the role of the microbiome in human biology and health.

Complete nucleotide sequence characterization of DRB5 alleles reveals a homogeneous allele group that is distinct from other DRB genes

Next Generation Sequencing allows for testing and typing of entire genes of the HLA region. A better and comprehensive sequence assessment can be achieved by the inclusion of full gene sequences of all the common alleles at a given locus. The common alleles of DRB5 are under-characterized with the full exon-intron sequence of two alleles available. In the present study the DRB5 genes from 18 subjects alleles were cloned and sequenced; haplotype analysis showed that 17 of them had a single copy of DRB5 and one consanguineous subject was homozygous at all HLA loci. Methodological approaches including robust and efficient long-range PCR amplification, molecular cloning, nucleotide sequencing and de novo sequence assembly were combined to characterize DRB5 alleles. DRB5 sequences covering from 5'UTR to the end of intron 5 were obtained for DRB5*01:01, 01:02 and 02:02; partial coverage including a segment spanning exon 2 to exon 6 was obtained for DRB5*01:03, 01:08N and 02:03. Phylogenetic analysis of the generated sequences showed that the DRB5 alleles group together and have distinctive differences with other DRB loci. Novel intron variants of DRB5*01:01:01, 01:02 and 02:02 were identified. The newly characterized DRB5 intron variants of each DRB5 allele were found in subjects harboring distinct associations with alleles of DRB1, B and/or ethnicity. The new information provided by this study provides reference sequences for HLA typing methodologies. Extending sequence coverage may lead to identify the disease susceptibility factors of DRB5 containing haplotypes while the unexpected intron variations may shed light on understanding of the evolution of the DRB region.

Functional characterization of enhancer evolution in the primate lineage

BACKGROUND

Enhancers play an important role in morphological evolution and speciation by controlling the spatiotemporal expression of genes. Previous efforts to understand the evolution of enhancers in primates have typically studied many enhancers at low resolution, or single enhancers at high resolution. Although comparative genomic studies reveal large-scale turnover of enhancers, a specific understanding of the molecular steps by which mammalian or primate enhancers evolve remains elusive.

RESULTS

We identified candidate hominoid-specific liver enhancers from H3K27ac ChIP-seq data. After locating orthologs in 11 primates spanning around 40 million years, we synthesized all orthologs as well as computational reconstructions of 9 ancestral sequences for 348 active tiles of 233 putative enhancers. We concurrently tested all sequences for regulatory activity with STARR-seq in HepG2 cells. We observe groups of enhancer tiles with coherent trajectories, most of which can be potentially explained by a single gain or loss-of-activity event per tile. We quantify the correlation between the number of mutations along a branch and the magnitude of change in functional activity. Finally, we identify 84 mutations that correlate with functional changes; these are enriched for cytosine deamination events within CpGs.

CONCLUSIONS

We characterized the evolutionary-functional trajectories of hundreds of liver enhancers throughout the primate phylogeny. We observe subsets of regulatory sequences that appear to have gained or lost activity. We use these data to quantify the relationship between sequence and functional divergence, and to identify CpG deamination as a potentially important force in driving changes in enhancer activity during primate evolution.

The Targeted Sequencing of Alpha Satellite DNA in Cercopithecus pogonias Provides New Insight Into the Diversity and Dynamics of Centromeric Repeats in Old World Monkeys

Alpha satellite is the major repeated DNA element of primate centromeres. Specific evolutionary mechanisms have led to a great diversity of sequence families with peculiar genomic organization and distribution, which have till now been studied mostly in great apes. Using high throughput sequencing of alpha satellite monomers obtained by enzymatic digestion followed by computational and cytogenetic analysis, we compare here the diversity and genomic distribution of alpha satellite DNA in two related Old World monkey species, Cercopithecus pogonias and Cercopithecus solatus, which are known to have diverged about 7 Ma. Two main families of monomers, called C1 and C2, are found in both species. A detailed analysis of our data sets revealed the existence of numerous subfamilies within the centromeric C1 family. Although the most abundant subfamily is conserved between both species, our fluorescence in situ hybridization (FISH) experiments clearly show that some subfamilies are specific for each species and that their distribution is restricted to a subset of chromosomes, thereby pointing to the existence of recurrent amplification/homogenization events. The pericentromeric C2 family is very abundant on the short arm of all acrocentric chromosomes in both species, pointing to specific mechanisms that lead to this distribution. Results obtained using two different restriction enzymes are fully consistent with a predominant monomeric organization of alpha satellite DNA that coexists with higher order organization patterns in the C. pogonias genome. Our study suggests a high dynamics of alpha satellite DNA in Cercopithecini, with recurrent apparition of new sequence variants and interchromosomal sequence transfer.

Superoxide dismutase 1 is positively selected to minimize protein aggregation in great apes

Positive (adaptive) selection has recently been implied in human superoxide dismutase 1 (SOD1), a highly abundant antioxidant protein with energy signaling and antiaging functions, one of very few examples of direct selection on a human protein product (exon); the molecular drivers of this selection are unknown. We mapped 30 extant SOD1 sequences to the recently established mammalian species tree and inferred ancestors, key substitutions, and signatures of selection during the protein's evolution. We detected elevated substitution rates leading to great apes (Hominidae) at ~1 per 2 million years, significantly higher than in other primates and rodents, although these paradoxically generally evolve much faster. The high evolutionary rate was partly due to relaxation of some selection pressures and partly to distinct positive selection of SOD1 in great apes. We then show that higher stability and net charge and changes at the dimer interface were selectively introduced upon separation from old world monkeys and lesser apes (gibbons). Consequently, human, chimpanzee and gorilla SOD1s have a net charge of -6 at physiological pH, whereas the closely related gibbons and macaques have -3. These features consistently point towards selection against the malicious aggregation effects of elevated SOD1 levels in long-living great apes. The findings mirror the impact of human SOD1 mutations that reduce net charge and/or stability and cause ALS, a motor neuron disease characterized by oxidative stress and SOD1 aggregates and triggered by aging. Our study thus marks an example of direct selection for a particular chemical phenotype (high net charge and stability) in a single human protein with possible implications for the evolution of aging.

Genetic Diversity and Structure among Isolated Populations of the Endangered Gees Golden Langur in Assam, India

Gee’s golden langur (Trachypithecus geei) is an endangered colobine primate, endemic to the semi-evergreen and mixed-deciduous forests of Indo-Bhutan border. During the last few decades, extensive fragmentation has caused severe population decline and local extinction of golden langur from several fragments. However, no studies are available on the impact of habitat fragmentation and the genetic diversity of golden langur in the fragmented habitats. The present study aimed to estimate the genetic diversity in the Indian population of golden langur. We sequenced and analyzed around 500 bases of the mitochondrial DNA (mtDNA) hypervariable region-I from 59 fecal samples of wild langur collected from nine forest fragments. Overall, genetic diversity was high (h = 0.934, π = 0.0244) and comparable with other colobines. Populations in smaller fragments showed lower nucleotide diversity compared to the larger forest fragments. The median-joining network of haplotypes revealed a genetic structure that corresponded with the geographical distribution. The Aie and Champabati Rivers were found to be a barrier to gene flow between golden langur populations. In addition, it also established that T. geei is monophyletic but revealed possible hybridization with capped langur, T. pileatus, in the wild. It is hoped that these findings would result in a more scientific approach towards managing the fragmented populations of this enigmatic species.

Old world monkeys and new age science: the evolution of nonhuman primate systems virology

Nonhuman primate (NHP) biomedical models are critical to our understanding of human health and disease, yet we are still in the early stages of developing sufficient tools to support primate genomic research that allow us to better understand the basis of phenotypic traits in NHP models of disease. A mere 7 years ago, the limited NHP transcriptome profiling that was being performed was done using complementary DNA arrays based on human genome sequences, and the lack of NHP genomic information and immunologic reagents precluded the use of NHPs in functional genomic studies. Since then, significant strides have been made in developing genomics capabilities for NHP research, from the rhesus macaque genome sequencing project to the construction of the first macaque-specific high-density oligonucleotide microarray, paving the way for further resource development and additional primate sequencing projects. Complete published draft genome sequences are now available for the chimpanzee ( Chimpanzee Sequencing Analysis Consortium 2005), bonobo ( Prufer et al. 2012), gorilla ( Scally et al. 2012), and baboon ( Ensembl.org 2013), along with the recently completed draft genomes for the cynomolgus macaque and Chinese rhesus macaque. Against this backdrop of both expanding sequence data and the early application of sequence-derived DNA microarrays tools, we will contextualize the development of these community resources and their application to infectious disease research through a literature review of NHP models of acquired immune deficiency syndrome and models of respiratory virus infection. In particular, we will review the use of -omics approaches in studies of simian immunodeficiency virus and respiratory virus pathogenesis and vaccine development, emphasizing the acute and innate responses and the relationship of these to the course of disease and to the evolution of adaptive immunity.

Differential evolutionary fate of an ancestral primate endogenous retrovirus envelope gene, the EnvV syncytin, captured for a function in placentation

Syncytins are envelope genes of retroviral origin that have been co-opted for a role in placentation. They promote cell–cell fusion and are involved in the formation of a syncytium layer—the syncytiotrophoblast—at the materno-fetal interface. They were captured independently in eutherian mammals, and knockout mice demonstrated that they are absolutely required for placenta formation and embryo survival. Here we provide evidence that these “necessary” genes acquired “by chance” have a definite lifetime with diverse fates depending on the animal lineage, being both gained and lost in the course of evolution. Analysis of a retroviral envelope gene, the envV gene, present in primate genomes and belonging to the endogenous retrovirus type V (ERV-V) provirus, shows that this captured gene, which entered the primate lineage >45 million years ago, behaves as a syncytin in Old World monkeys, but lost its canonical fusogenic activity in other primate lineages, including humans. In the Old World monkeys, we show—by in situ analyses and ex vivo assays—that envV is both specifically expressed at the level of the placental syncytiotrophoblast and fusogenic, and that it further displays signs of purifying selection based on analysis of non-synonymous to synonymous substitution rates. We further show that purifying selection still operates in the primate lineages where the gene is no longer fusogenic, indicating that degeneracy of this ancestral syncytin is a slow, lineage-dependent, and multi-step process, in which the fusogenic activity would be the first canonical property of this retroviral envelope gene to be lost.

Syncytins are “new” genes encoding the envelope protein of captured endogenous retroviral elements. Their unambiguous status of “cellular gene” was recently demonstrated by knocking them out in genetically modified mice, showing their absolute requirement for placenta formation and embryo survival, via formation by cell–cell fusion of the feto-maternal syncytium interface. These genes are remarkable, as they are “necessary” for a basic function in placental mammals and yet they were acquired “by chance” on multiple occasions and independently in diverse mammalian species. We proposed that syncytins have been pivotal for the emergence of animals with a placenta from those laying eggs via the capture of a founding retroviral env gene, then subsequently replaced in the diverse mammalian lineages upon successive and independent germline infections by new retroviruses and co-optation of their env gene, each new gene providing its host with a positive selective advantage. This hypothesis would account for the diversity of the captured syncytins that can be currently found, concomitant with the diversity of placental architectures. A consequence of this paradigm is that evidence for “decaying syncytins” in eutherian mammals should exist, and this is precisely what we sought—and found—in this study.

Cytogenetic studies of small ape (Hylobatidae) chromosomes

Each genus of small apes has a highly distinctive karyotype (karyomorph) at every level of cytogenetic analysis. Early workers using classical staining and banding had problems integrating the karyolocial data with that of other primates. Chromosome painting allowed syntenic homology maps to be constructed for each of the four karyomorphs (2n = 38, 44, 50 and 52). They revealed that the great apes and Old World monkeys had strongly conserved karyotypes while those of small apes were highly rearranged. However, they provided contradictory phylogenetic results to other bio-molecular tree of small ape evolution. More recently BAC-FISH investigations using a panel of about 900 BACs defined each breakpoint by spanning or flanking BAC clones The syntenic map was refined and now includes small segments of homology which had previously gone undected, marker order (synteny block orientation) and the location of ancestral and Evolutionarily New Centromeres. However, the BAC-FISH data similar to other biomolecular methods used up to now could not resolve the phylogenetic tree of hylobatids. These difficulties may be explained by the rapid divergence of crown hylobatids, reticulate evolution and incomplete lineage sorting. The lack of significant cytogenetic landmarks at the nodes of the gibbon tree could indicate that chromosomal rearrangements did not play a primary role in hylobatid speciation.

The ABO blood group is a trans-species polymorphism in primates

The ABO histo-blood group, the critical determinant of transfusion incompatibility, was the first genetic polymorphism discovered in humans. Remarkably, ABO antigens are also polymorphic in many other primates, with the same two amino acid changes responsible for A and B specificity in all species sequenced to date. Whether this recurrence of A and B antigens is the result of an ancient polymorphism maintained across species or due to numerous, more recent instances of convergent evolution has been debated for decades, with a current consensus in support of convergent evolution. We show instead that genetic variation data in humans and gibbons as well as in Old World monkeys are inconsistent with a model of convergent evolution and support the hypothesis of an ancient, multiallelic polymorphism of which some alleles are shared by descent among species. These results demonstrate that the A and B blood groups result from a trans-species polymorphism among distantly related species and has remained under balancing selection for tens of millions of years-to date, the only such example in hominoids and Old World monkeys outside of the major histocompatibility complex.

Evolutionary history of the odd-nosed monkeys and the phylogenetic position of the newly described Myanmar snub-nosed monkey Rhinopithecus strykeri

Odd-nosed monkeys represent one of the two major groups of Asian colobines. Our knowledge about this primate group is still limited as it is highlighted by the recent discovery of a new species in Northern Myanmar. Although a common origin of the group is now widely accepted, the phylogenetic relationships among its genera and species, and the biogeographic processes leading to their current distribution are largely unknown. To address these issues, we have analyzed complete mitochondrial genomes and 12 nuclear loci, including one X chromosomal, six Y chromosomal and five autosomal loci, from all ten odd-nosed monkey species. The gene tree topologies and divergence age estimates derived from different markers were highly similar, but differed in placing various species or haplogroups within the genera Rhinopithecus and Pygathrix. Based on our data, Rhinopithecus represent the most basal lineage, and Nasalis and Simias form closely related sister taxa, suggesting a Northern origin of odd-nosed monkeys and a later invasion into Indochina and Sundaland. According to our divergence age estimates, the lineages leading to the genera Rhinopithecus, Pygathrix and Nasalis+Simias originated in the late Miocene, while differentiation events within these genera and also the split between Nasalis and Simias occurred in the Pleistocene. Observed gene tree discordances between mitochondrial and nuclear datasets, and paraphylies in the mitochondrial dataset for some species of the genera Rhinopithecus and Pygathrix suggest secondary gene flow after the taxa initially diverged. Most likely such events were triggered by dramatic changes in geology and climate within the region. Overall, our study provides the most comprehensive view on odd-nosed monkey evolution and emphasizes that data from differentially inherited markers are crucial to better understand evolutionary relationships and to trace secondary gene flow.

Evolution of the antiretroviral restriction factor TRIMCyp in Old World primates

The retroviral restriction factor TRIMCyp, which is a fusion protein derived from the TRIM5 gene, blocks replication at a post-entry step. Among Old World primates, TRIMCyp has been found in four species of Asian macaques, but not in African monkeys. To further define the evolutionary origin of Old World TRIMCyp, we examined two species of baboons (genus Papio) and three additional macaque species, including M. sylvanus, which is the only macaque species found outside Asia, and represents the earliest diverging branch of the macaque lineage. None of four P. cynocephalus anubis, one P. hamadryas, and 36 M. sylvanus had either TRIMCyp mRNA or the genetic features required for its expression. M. sylvanus genomic sequences indicated that the lack of TRIMCyp in this species was not due to genetic homogeneity among specimens studied and revealed the existence of four TRIM5α alleles, all distinct from M. mulatta and Papio counterparts. Together with existing data on macaque evolution, our findings indicate that TRIMCyp evolved in the ancestors of Asian macaques approximately 5-6 million years before present (ybp), likely as a result of a retroviral threat. TRIMCyp then became fixed in the M. nemestrina lineage after it diverged from M. nigra, approximately 2 million ybp. The macaque lineage is unique among primates studied so far due to the presence and diversity of both TRIM5 and TRIMCyp restriction factors. Studies of these antiviral proteins may provide valuable information about natural antiviral mechanisms, and give further insight into the factors that shaped the evolution of macaque species.

Noncoding sequences near duplicated genes evolve rapidly

Gene expression divergence and chromosomal rearrangements have been put forward as major contributors to phenotypic differences between closely related species. It has also been established that duplicated genes show enhanced rates of positive selection in their amino acid sequences. If functional divergence is largely due to changes in gene expression, it follows that regulatory sequences in duplicated loci should also evolve rapidly. To investigate this hypothesis, we performed likelihood ratio tests (LRTs) on all noncoding loci within 5 kb of every transcript in the human genome and identified sequences with increased substitution rates in the human lineage since divergence from Old World Monkeys. The fraction of rapidly evolving loci is significantly higher nearby genes that duplicated in the common ancestor of humans and chimps compared with nonduplicated genes. We also conducted a genome-wide scan for nucleotide substitutions predicted to affect transcription factor binding. Rates of binding site divergence are elevated in noncoding sequences of duplicated loci with accelerated substitution rates. Many of the genes associated with these fast-evolving genomic elements belong to functional categories identified in previous studies of positive selection on amino acid sequences. In addition, we find enrichment for accelerated evolution nearby genes involved in establishment and maintenance of pregnancy, processes that differ significantly between humans and monkeys. Our findings support the hypothesis that adaptive evolution of the regulation of duplicated genes has played a significant role in human evolution.

High diversity in functional properties of melanocortin 1 receptor (MC1R) in divergent primate species is more strongly associated with phylogeny than coat color

We have characterized the biochemical function of the melanocortin 1 receptor (MC1R), a critical regulator of melanin synthesis, from 9 phylogenetically diverse primate species with varying coat colors. There is substantial diversity in melanocyte-stimulating hormone (MSH) binding affinity and basal levels of activity in the cloned MC1Rs. MSH binding was lost independently in lemur and New World monkey lineages, whereas high basal levels of MC1R activity occur in lemurs and some New World monkeys and Old World monkeys. Highest levels of basal activity were found in the MC1R of ruffed lemurs, which have the E94K mutation that leads to constitutive activation in other species. In 3 species (2 lemurs and the howler monkey), we report the novel finding that binding and inhibition of MC1R by agouti signaling protein (ASIP) can occur when MSH binding has been lost, thus enabling continuing regulation of the melanin type via ASIP expression. Together, these findings can explain the previous paradox of a predominantly pheomelanic coat in the red ruffed lemur (Varecia rubra). The presence of a functional, MSH-responsive MC1R in orangutan demonstrates that the mechanism of red hair generation in this ape is different from the prevalent mechanism in European human populations. Overall, we have found unexpected diversity in MC1R function among primates and show that the evolution of the regulatory control of MC1R activity occurs by independent variation of 3 distinct mechanisms: basal MC1R activity, MSH binding and activation, and ASIP binding and inhibition. This diversity of function is broadly associated with primate phylogeny and does not have a simple relation to coat color phenotype within primate clades.

Genetic relationships of langur species using AFLP markers

Cytogenetic studies of five langur species using conventional banding pattern were investigated. All species studied have an identical number of 44 diploid chromosomes, they are assumed to have common evolutionary relationships. For in depth study, molecular markers were assessed using the Amplified Fragment Length Polymorphism (AFLP) method. With seven successful primer combinations, a total of 1043 scorable bands were generated. The percentage of polymorphic bands for each primer ranged from 48.60 to 94.12%. The resulting bands were used for dendrogram construction. From the dendrogram, the individuals of Trachypithecus species are grouped into two major clusters, T. phayrei is clustered with T. obscurus, while T. cristatus is clustered with T. francoisi. The bootstrap value between two groups is 94%. The other cluster, Presbytis femoralis is separated from the Trachypithecus species with a bootstrap value of 94%. Averages of inter-specific genetic similarity values among all langur species studied are 70.16% (between T. obscurus and P. femoralis) to 88.12% (between T. phayrei and T. obscurus). In summary T. phayrei might be a subspecies of T. obscurus. The development of specific molecular markers of a species is beneficial for genetic differentiation of this group of primates.

Mobile DNA in Old World monkeys: a glimpse through the rhesus macaque genome

The completion of the draft sequence of the rhesus macaque genome allowed us to study the genomic composition and evolution of transposable elements in this representative of the Old World monkey lineage, a group of diverse primates closely related to humans. The L1 family of long interspersed elements appears to have evolved as a single lineage, and Alu elements have evolved into four currently active lineages. We also found evidence of elevated horizontal transmissions of retroviruses and the absence of DNA transposon activity in the Old World monkey lineage. In addition, approximately 100 precursors of composite SVA (short interspersed element, variable number of tandem repeat, and Alu) elements were identified, with the majority being shared by the common ancestor of humans and rhesus macaques. Mobile elements compose roughly 50% of primate genomes, and our findings illustrate their diversity and strong influence on genome evolution between closely related species.

Cross-species analysis of gene expression in non-model mammals: reproducibility of hybridization on high density oligonucleotide microarrays

BACKGROUND

Gene expression profiles of non-model mammals may provide valuable data for biomedical and evolutionary studies. However, due to lack of sequence information of other species, DNA microarrays are currently restricted to humans and a few model species. This limitation may be overcome by using arrays developed for a given species to analyse gene expression in a related one, an approach known as "cross-species analysis". In spite of its potential usefulness, the accuracy and reproducibility of the gene expression measures obtained in this way are still open to doubt. The present study examines whether or not hybridization values from cross-species analyses are as reproducible as those from same-species analyses when using Affymetrix oligonucleotide microarrays.

RESULTS

The reproducibility of the probe data obtained hybridizing deer, Old-World primates, and human RNA samples to Affymetrix human GeneChip U133 Plus 2.0 was compared. The results show that cross-species hybridization affected neither the distribution of the hybridization reproducibility among different categories, nor the reproducibility values of the individual probes. Our analyses also show that a 0.5% of the probes analysed in the U133 plus 2.0 GeneChip are significantly associated to un-reproducible hybridizations. Such probes-called in the text un-reproducible probe sequences- do not increase in number in cross-species analyses.

CONCLUSION

Our study demonstrates that cross-species analyses do not significantly affect hybridization reproducibility of GeneChips, at least within the range of the mammal species analysed here. The differences in reproducibility between same-species and cross-species analyses observed in previous studies were probably caused by the analytical methods used to calculate the gene expression measures. Together with previous observations on the accuracy of GeneChips for cross-species analysis, our analyses demonstrate that cross-species hybridizations may provide useful gene expression data. However, the reproducibility and accuracy of these measures largely depends on the use of appropriated algorithms to derive the gene expression data from the probe data. Also, the identification of probes associated to un-reproducible hybridizations-useless for gene expression analyses- in the studied GeneChip, stress the need of a re-evaluation of the probes' performance.

Balancing selection and the evolution of functional polymorphism in Old World monkey TRIM5alpha

Retroviral restriction factor TRIM5alpha exhibits a high degree of sequence variation among primate species. It has been proposed that this diversity is the cumulative result of ancient, lineage-specific episodes of positive selection. Here, we describe the contribution of within-species variation to the evolution of TRIM5alpha. Sampling within two geographically distinct Old World monkey species revealed extensive polymorphism, including individual polymorphisms that predate speciation (shared polymorphism). In some instances, alleles were more closely related to orthologues of other species than to one another. Both silent and nonsynonymous changes clustered in two domains. Functional assays revealed consequences of polymorphism, including differential restriction of a small panel of retroviruses by very similar alleles. Together, these features indicate that the primate TRIM5alpha locus has evolved under balancing selection. Except for the MHC there are few, if any, examples of long-term balancing selection in primates. Our results suggest a complex evolutionary scenario, in which fixation of lineage-specific adaptations is superimposed on a subset of critical polymorphisms that predate speciation events and have been maintained by balancing selection for millions of years.

Independent centromere formation in a capricious, gene-free domain of chromosome 13q21 in Old World monkeys and pigs

BACKGROUND

Evolutionary centromere repositioning and human analphoid neocentromeres occurring in clinical cases are, very likely, two stages of the same phenomenon whose properties still remain substantially obscure. Chromosome 13 is the chromosome with the highest number of neocentromeres. We reconstructed the mammalian evolutionary history of this chromosome and characterized two human neocentromeres at 13q21, in search of information that could improve our understanding of the relationship between evolutionarily new centromeres, inactivated centromeres, and clinical neocentromeres.

RESULTS

Chromosome 13 evolution was studied, using FISH experiments, across several diverse superordinal phylogenetic clades spanning >100 million years of evolution. The analysis revealed exceptional conservation among primates (hominoids, Old World monkeys, and New World monkeys), Carnivora (cat), Perissodactyla (horse), and Cetartiodactyla (pig). In contrast, the centromeres in both Old World monkeys and pig have apparently repositioned independently to a central location (13q21). We compared these results to the positions of two human 13q21 neocentromeres using chromatin immunoprecipitation and genomic microarrays.

CONCLUSION

We show that a gene-desert region at 13q21 of approximately 3.9 Mb in size possesses an inherent potential to form evolutionarily new centromeres over, at least, approximately 95 million years of mammalian evolution. The striking absence of genes may represent an important property, making the region tolerant to the extensive pericentromeric reshuffling during subsequent evolution. Comparison of the pericentromeric organization of chromosome 13 in four Old World monkey species revealed many differences in sequence organization. The region contains clusters of duplicons showing peculiar features.

Heterogeneous genomic molecular clocks in primates

Using data from primates, we show that molecular clocks in sites that have been part of a CpG dinucleotide in recent past (CpG sites) and non-CpG sites are of markedly different nature, reflecting differences in their molecular origins. Notably, single nucleotide substitutions at non-CpG sites show clear generation-time dependency, indicating that most of these substitutions occur by errors during DNA replication. On the other hand, substitutions at CpG sites occur relatively constantly over time, as expected from their primary origin due to methylation. Therefore, molecular clocks are heterogeneous even within a genome. Furthermore, we propose that varying frequencies of CpG dinucleotides in different genomic regions may have contributed significantly to conflicting earlier results on rate constancy of mammalian molecular clock. Our conclusion that different regions of genomes follow different molecular clocks should be considered when inferring divergence times using molecular data and in phylogenetic analysis.

Growth hormone locus expands and diverges after the separation of New and Old World Monkeys

While most mammals including the prosimians have a single copy of the growth hormone (GH) gene, anthropoids possess a cluster of GH-related genes. Throughout the evolution of the main anthropoid groups [New World Monkeys (NWM), Old World Monkeys (OWM), and apes], two features stand out of the GH loci. The first is the appearance of chorionic somatommamotropin hormone (CSH) genes within the OWM lineage and the second is the expansion of the loci intergenic regions in the OWM and apes. In relation with this loci expansion, the NWM possess intergenic regions of homogeneous lengths (3.5 kb). In contrast, heterogeneous lengths (6 and 13 kb) have been reported for species of the OWM. At the present, none of the OWM genomic GH loci organizations have been described. Here, we report the genomic organization of the GH locus in the rhesus monkey, this locus has six GH-related genes separated by five intergenic regions. The 5' end gene (GH-1) encodes for the pituitary GH and is followed by CSH-1, GH-2, CSH-2, CSH-3 and CSH-4 genes. The five intergenic regions have heterogeneous lengths and also present more or less the same Alu distribution as the human GH locus. To analyze the events that contributed to the extension of the intergenic regions of the GH locus and the emergence of the regulatory elements, the five GH locus intergenic regions of the spider monkey (NWM) were sequenced. The results of comparing the loci from both species suggest that the long intergenic regions (13 kb) of the rhesus GH locus share a common ancestor with the 3.5 kb intergenic regions of the spider monkey. However, the observed increased length of the former is due to an insertion (approximately 8.7 kb) at their 3' end. Interestingly in this insert, we discovered a DNA element resembling the enhancer of the CSH genes of the human GH locus. On the other hand, we observed that the short intergenic regions (6 kb) increased by a different recombination event.

Multiplexed microsatellites for rapid identification and characterization of individuals and populations of Cercopithecidae

Cross-amplification of 15 human microsatellites was performed successfully in cynomolgus (Macaca fascicularis) and rhesus (M. mulatta) macaques and 11 other Cercopithecidae species of biomedical and conservation relevance. To allow for quick, efficient, and high-throughput genotyping to assess intra- and interspecific genetic variation, we performed three multiplex sets, each comprised of five markers from different parts of the genome (i.e., autosomes, the MHC region, and the X-chromosome). These multiplex sets are likely to reveal allelic divergence between taxa, which could be used for their discrimination. Population studies on three regional populations of M. fascicularis and one of M. mulatta revealed that most of the loci, with the exception of one monomorphic locus, displayed polymorphisms (the expected heterozygosities were 0.48-0.91 for M. fascicularis, and 0.61-0.93 for M. mulatta), which makes them useful for population genetics. For the multiplex set M1, including the nonlinked autosomal markers, low probabilities of identity were observed: P(ID) values ranged from 8 x 10(-7) to 3 x 10(-5). This multiplex set is reliable for forensic applications, such as individual identification, parentage testing, and kinship analysis, in wild and captive populations.

Molecular characterization of immunoglobulin D in mammals: immunoglobulin heavy constant delta genes in dogs, chimpanzees and four old world monkey species

Antibodies are adaptor molecules that neutralize pathogens and link humoral and cellular defence mechanisms. Immunoglobulin D (IgD), one of the five antibody classes present in mammals, is expressed as an antigen receptor on naïve B cells. The functional role that IgD plays in the immune response is still poorly understood, but the recent characterization of immunoglobulin heavy constant delta genes (IGHD) in a variety of species challenges the view that IgD is of minor importance and is not present in many animals. On the basis of serological studies, IgD appears to be expressed in the majority of mammalian species examined. To confirm, at the molecular level, that IgD is present in different species, we cloned and sequenced IGHD cDNA from dogs and five non-human primate species (chimpanzee, rhesus macaque, cynomolgus macaque, baboon and sooty mangabey). Our results show that in all six species, IgD heavy chains possess three immunoglobulin domains and a long hinge region encoded by two exons. Only the hinge region of non-human primates is similar to the human hinge region, with conservation of O-glycosylation sites and multiple charged residues at opposing ends. The preservation of IgD in primates, dogs and previously characterized species suggests an important functional role for IgD, possibly involving binding to a receptor. The high degree of similarity existing between the structural features of human and non-human primate IgD suggests that non-human primates are suitable for in vivo studies designed to define the role that IgD plays in the immune response.

Fast and non-invasive PCR sexing of primates: apes, Old World monkeys, New World monkeys and Strepsirrhines

Background

One of the key tools for determining the social structure of wild and endangered primates is the ability to sex DNA from small amounts of non-invasive samples that are likely to include highly degraded DNA. Traditional markers for molecular sex determination of primates are developed on the basis of the human sequence and are often non-functional in distantly related primate species. Hence, it is highly desirable to develop markers that simultaneously detect Y- and X-chromosome specific sequences and also work across many species.

Results

A novel method for sex identification in primates is described using a triple primer PCR reaction and agarose gel electrophoresis of the sex-chromosomal isoforms of the ubiquitously transcribed tetratricopeptide repeat protein gene (UTX/UTY). By comparing genomic data from several mammals we identified the UTX/UTY locus as the best candidate for a universal primate sexing marker. Using data from several species we identified a XY-conserved region, a Y conserved region and an X conserved region. This enabled the design of a triple primer PCR setup that amplifies X and Y products of different length in a single PCR reaction.

Conclusion

This simple PCR amplification of X and Y fragments is useful for sexing DNA samples from all species of primates. Furthermore, since the amplified fragments are very short the method can be applied to fragmented DNA extracted from non-invasive samples.

Pervasive adaptive evolution in primate seminal proteins

Seminal fluid proteins show striking effects on reproduction, involving manipulation of female behavior and physiology, mechanisms of sperm competition, and pathogen defense. Strong adaptive pressures are expected for such manifestations of sexual selection and host defense, but the extent of positive selection in seminal fluid proteins from divergent taxa is unknown. We identified adaptive evolution in primate seminal proteins using genomic resources in a tissue-specific study. We found extensive signatures of positive selection when comparing 161 human seminal fluid proteins and 2,858 prostate-expressed genes to those in chimpanzee. Seven of eight outstanding genes yielded statistically significant evidence of positive selection when analyzed in divergent primates. Functional clues were gained through divergent analysis, including several cases of species-specific loss of function in copulatory plug genes, and statistically significant spatial clustering of positively selected sites near the active site of kallikrein 2. This study reveals previously unidentified positive selection in seven primate seminal proteins, and when considered with findings in Drosophila, indicates that extensive positive selection is found in seminal fluid across divergent taxonomic groups.

Proteins found in seminal fluid accompanying sperm show dramatic effects on reproduction, such as manipulating female behavior. Even in primates they participate in competition between sperm of different males, and serve to protect sperm from infection by pathogens. These types of roles require the proteins to constantly adapt to stay ahead of the competition. Such adaptive pressures on proteins leave characteristic signatures in the DNA sequences that encode them. The authors used these signatures to identify adaptive evolution in primate seminal proteins and found extensive signs of adaptation when comparing thousands of seminal genes between human and chimpanzee. They further characterized outstanding genes in several primate species, including a diversity of apes and monkeys. Several of these proteins have no known function, yet by visualizing the adaptation on their three-dimensional surfaces, the authors uncovered clues to what is driving their evolution. In addition, they found several cases in which certain species lost their functional copies of these genes. Interestingly, species that showed loss of function do not participate in sperm competition. Past studies found widespread adaptation in fruit fly seminal fluid, and this study reveals extensive adaptation in primate seminal proteins. Could this be a phenomenon common among animals?

[Population genetics of Rhinopithecus bieti: a study of the mitochondrial control region]

Yunnan snub-nose monkey (Rhinopithecus bieti) is a famous endangered primate in China. So far, however, studies on its population genetics based on DNA sequences are not available. In this paper, the whole mitochondria control region of the samples from Weixi, Yunan Province as well as the whole cytochrome b gene in some individuals were sequenced. A deep divergence was observed within the Weixi population, which was confirmed after excluding the possibility of it being a nuclear pseudogene. Nonetheless, if the effects of population structure and migration are considered, the true level of polymorphism of the Weixi population may be not as high as observed.

Evolutionary scenario for acquisition of CAG repeats in human SCA1 gene

We investigated the CAG repeat sequence of the spinocerebellar ataxia type 1 (SCA1) gene in various species of primates to reveal how human has acquired the repeat structure with interruptions. Our results demonstrate no repetitive structure in the region corresponding to the human CAG repeats in prosimians and New World monkeys like in rodents, perfect (uninterrupted) CAG repeats in Old World monkeys, and interrupted CAG repeats in hominoids. Comparative analysis on the secondary structures of the primate SCA1 transcripts suggests the human prototype was built in the common ancestor of simians. We show an evolutionary scenario for acquisition of CAG repeats with interruptions in the human SCA1 gene.

The gene of retroviral origin Syncytin 1 is specific to hominoids and is inactive in Old World monkeys

Syncytin 1 is one of the best known examples of recent acquisition of a new gene from an endogenous retrovirus (HERV) in the human genome and has been implicated in placental physiology. Within primates, Syncytin 1 is conserved in all hominoids but has not been characterized in Old World monkeys (OWMs). In this study, we investigated the status of Syncytin 1 in 14 hominoid and OWM species. We show that although the HERV-W provirus responsible for the origin of this gene was present in the genome of the most recent common ancestor of hominoids and OWMs, Syncytin 1 is inactive in OWMs. In addition, we were able to determine that the evolution of Syncytin 1 in hominoids involved an accumulation of amino acid changes and showed signatures of both positive and purifying selection. Our results indicate that Syncytin 1 is indeed a hominoid-specific gene and illustrate the complex and dynamic process associated with the origin of new genes.

Human endogenous retroviruses: from infectious elements to human genes

Mammalian genomes contain a heavy load (42% in humans) of retroelements, which are mobile sequences requiring reverse transcription for their replicative transposition. A significant proportion of these elements is of retroviral origin, with thousands of sequences resembling the integrated form of infectious retroviruses, with two LTRs bordering internal regions homologous to the gag, prt, pol, and env genes. These elements, named endogenous retroviruses (ERVs), are most probably the proviral remnants of ancestral germ-line infections by active retroviruses, which have thereafter been transmitted in a Mendelian manner. The complete sequencing of the human genome now allows a comprehensive survey of human ERVs (HERVs), which can be grouped according to sequence homologies into approximately 80 distinct families, each containing a few to several hundred elements. As reviewed here, strong similarities between HERVs and present-day retroviruses can be inferred from phylogenetic analyses on the reverse transcriptase (RT) domain of the pol gene or the transmembrane subunit (TM) of the env gene, which disclose interspersion of both classes of elements and suggest a common history and shared ancestors. Similarities are also observed at the functional levels, since despite the fact that most HERVs have accumulated mutations, deletions, and/or truncations, several elements still possess some of the functions of retroviruses, with evidence for viral-like particle formation, and occurrence of envelope proteins allowing cell-cell fusion and even conferring infectivity to pseudotypes. Along this line, a genomewide screening for human retroviral genes with coding capacity has revealed 16 fully coding envelope genes. These genes are transcribed in several healthy tissues including the placenta, three of them at a very high level. Besides their impact in modelling the genome, HERVs thus appear to contain still active genes, which most probably have been subverted by the host for its benefit and should be considered as bona fide human genes. Some of their characteristic features and possible physiological roles, as well as potential pathological effects inherited from their retroviral ancestors are also reviewed.

Human SNPs reveal no evidence of frequent positive selection

We compared the single-nucleotide polymorphisms (SNPs) in humans in 182 housekeeping and 148 tissue-specific genes. SNPs were divided into rare and common polymorphisms based on their frequencies. We found that housekeeping genes tend to be less polymorphic than tissue-specific genes for both rare and common SNPs. Using mouse as a second species for computing sequence divergences, we found no evidence of positive selection: for both housekeeping and tissue-specific genes, the ratio of nonsynonymous to synonymous common SNPs per site showed no significant difference from that of divergence. Similarly, we observed no evidence of positive selection for the 289 and 149 genes that have orthologs available for divergence calculation between humans and chimpanzees and between humans and Old World monkeys, respectively. A comparison with previous SNP studies suggests that approximately 20% of the nonsynonymous SNPs in the human population are nearly neutral and that positive selection in the human genome might not be as frequent as previously thought.

Study of the mitochondrial transcription factor A (Tfam) gene in the primate Presbytis cristata

The mitochondrial transcription factor A (Tfam) is a member of the HMG-box protein family, necessary for both transcription and maintenance of mitochondrial DNA. The gene is structured in seven exons and six introns and it is estimated to span about 10 kb in mouse, human and rat. In addition to the full length mRNA of Tfam, a shorter mRNA isoform lacking exon 5 has been found to be widely distributed in human and rat tissues. Here we present the isolation and characterization of Tfam gene in the primate Presbytis cristata which belongs to the Cercopithecidae family. We have determined the complete CDS sequence, the size of all the six introns, the complete sequences of the three shorter ones (I, III, VI) and the partial sequences of the long introns (II, IV, V). The comparison with other available Tfam sequences from mammals has revealed a high degree of conservation (above 90%) both in CDS and introns. By in situ hybridization (FISH) experiments we have mapped Tfam gene on chromosome 12 which, according to other cytogenetics studies, is the homologous region of chromosome 10, where human Tfam has been mapped. Moreover we have searched for the presence of alternatively spliced isoforms through several approaches, such as RT-PCR and differential hybridization. In Presbytis cristata we have not detected the presence of any spliced isoforms lacking exons; however we have identified one isoform in which part of the intron I is retained in the mRNA. The inclusion of this portion of intron I would originate an early stop codon if translated.

Genetic evidence of a strong functional constraint of neurotrypsin during primate evolution

Neurotrypsin is one of the extra-cellular serine proteases that are predominantly expressed in the brain and involved in neuronal development and function. Mutations in humans are associated with autosomal recessive non-syndromic mental retardation (MR). We studied the molecular evolution of neurotrypsin by sequencing the coding region of neurotrypsin in 11 representative non-human primate species covering great apes, lesser apes, Old World monkeys and New World monkeys. Our results demonstrated a strong functional constraint of neurotrypsin that was caused by strong purifying selection during primate evolution, an implication of an essential functional role of neurotrypsin in primate cognition. Further analysis indicated that the purifying selection was in fact acting on the SRCR domains of neurotrypsin, which mediate the binding activity of neurotrypsin to cell surface or extra-cellular proteins. In addition, by comparing primates with three other mammalian orders, we demonstrated that the absence of the first copy of the SRCR domain (exon 2 and 3) in mouse and rat was due to the deletion of this segment in the murine lineage.

Genomic structure and paralogous regions of the inversion breakpoint occurring between human chromosome 3p12.3 and orangutan chromosome 2

Intrachromosomal duplications play a significant role in human genome pathology and evolution. To better understand the molecular basis of evolutionary chromosome rearrangements, we performed molecular cytogenetic and sequence analyses of the breakpoint region that distinguishes human chromosome 3p12.3 and orangutan chromosome 2. FISH with region-specific BAC clones demonstrated that the breakpoint-flanking sequences are duplicated intrachromosomally on orangutan 2 and human 3q21 as well as at many pericentromeric and subtelomeric sites throughout the genomes. Breakage and rearrangement of the human 3p12.3-homologous region in the orangutan lineage were associated with a partial loss of duplicated sequences in the breakpoint region. Consistent with our FISH mapping results, computational analysis of the human chromosome 3 genomic sequence revealed three 3p12.3-paralogous sequence blocks on human chromosome 3q21 and smaller blocks on the short arm end 3p26-->p25. This is consistent with the view that sequences from an ancestral site at 3q21 were duplicated at 3p12.3 in a common ancestor of orangutan and humans. Our results show that evolutionary chromosome rearrangements are associated with microduplications and microdeletions, contributing to the DNA differences between closely related species.

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.

Molecular evolution of growth hormone gene family in old world monkeys and hominoids

Growth hormone is a classic molecule in the study of the molecular clock hypothesis as it exhibits a relatively constant rate of evolution in most mammalian orders except primates and artiodactyls, where dramatically enhanced rate of evolution (25-50-fold) has been reported. The rapid evolution of primate growth hormone occurred after the divergence of tarsiers and simians, but before the separation of old world monkeys (OWM) from new world monkeys (NWM). Interestingly, this event of rapid sequence evolution coincided with multiple duplications of the growth hormone gene, suggesting gene duplication as a possible cause of the accelerated sequence evolution. Here we determined 21 different GH-like sequences from four species of OWM and hominoids. Combining with published sequences from OWM and hominoids, our analysis demonstrates that multiple gene duplications and several gene conversion events both occurred in the evolutionary history of this gene family in OWM/hominoids. The episode of recent duplications of CSH-like genes in gibbon is accompanied with rapid sequence evolution likely resulting from relaxation of purifying selection. GHN genes in both hominoids and OWM are under strong purifying selection. In contrast, CSH genes in both lineages are probably not. GHV genes in OWM and hominoids evolved at different evolutionary rates and underwent different selective constraints. Our results disclosed the complex history of the primate growth hormone gene family and raised intriguing questions on the consequences of these evolutionary events.

Ancestral population sizes and species divergence times in the primate lineage on the basis of intron and BAC end sequences

The effective sizes of ancestral populations and species divergence times of six primate species (humans, chimpanzees, gorillas, orangutans, and representatives of Old World monkeys and New World monkeys) are estimated by applying the two-species maximum likelihood (ML) method to intron sequences of 20 different loci. Examination of rate heterogeneity of nucleotide substitutions and intragenic recombination identifies five outrageous loci (ODC1, GHR, HBE, INS, and HBG). The estimated ancestral polymorphism ranges from 0.21 to 0.96% at major divergences in primate evolution. One exceptionally low polymorphism occurs when African and Asian apes diverged. However, taking into consideration the possible short generation times in primate ancestors, it is concluded that the ancestral population size in the primate lineage was no smaller than that of extant humans. Furthermore, under the assumption of 6 million years (myr) divergence between humans and chimpanzees, the divergence time of humans from gorillas, orangutans. Old World monkeys, and New World monkeys is estimated as 7.2, 18, 34, and 65 myr ago, respectively, which are generally older than traditional estimates. Beside the intron sequences, three other data sets of orthologous sequences are used between the human and the chimpanzee comparison. The ML application to these data sets including 58,156 random BAC end sequences (BES) shows that the nucleotide substitution rate is as low as 0.6-0.8 x 10(-9) per site per year and the extent of ancestral polymorphism is 0.33-0.51%. With such a low substitution rate and short generation time, the relatively high extent of polymorphism suggests a fairly large effective population size in the ancestral lineage common to humans and chimpanzees.

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.

Genomic data support the hominoid slowdown and an Early Oligocene estimate for the hominoid-cercopithecoid divergence

Several lines of indirect evidence suggest that hominoids (apes and humans) and cercopithecoids (Old World monkeys) diverged around 23-25 Mya. Importantly, although this range of dates has been used as both an initial assumption and as a confirmation of results in many molecular-clock analyses, it has not been critically assessed on its own merits. In this article we test the robusticity of the 23- to 25-Mya estimate with approximately 150,000 base pairs of orthologous DNA sequence data from two cercopithecoids and two hominoids by using quartet analysis. This method is an improvement over other estimates of the hominoid-cercopithecoid divergence because it incorporates two calibration points, one each within cercopithecoids and hominoids, and tests for a statistically appropriate model of molecular evolution. Most comparisons reject rate constancy in favor of a model incorporating two rates of evolution, supporting the "hominoid slowdown" hypothesis. By using this model of molecular evolution, the hominoid-cercopithecoid divergence is estimated to range from 29.2 to 34.5 Mya, significantly older than most previous analyses. Hominoid-cercopithecoid divergence dates of 23-25 Mya fall outside of the confidence intervals estimated, suggesting that as much as one-third of ape evolution has not been paleontologically sampled. Identifying stem cercopithecoids or hominoids from this period will be difficult because derived features that define crown catarrhines need not be present in early members of these lineages. More sites that sample primate habitats from the Oligocene of Africa are needed to better understand early ape and Old World monkey evolution.

Evolutionary conserved chromosomal segments in the human karyotype are bounded by unstable chromosome bands

In this paper an ancestral karyotype for primates, defining for the first time the ancestral chromosome morphology and the banding patterns, is proposed, and the ancestral syntenic chromosomal segments are identified in the human karyotype. The chromosomal bands that are boundaries of ancestral segments are identified. We have analyzed from data published in the literature 35 different primate species from 19 genera, using the order Scandentia, as well as other published mammalian species as out-groups, and propose an ancestral chromosome number of 2n = 54 for primates, which includes the following chromosomal forms: 1(a+c(1)), 1(b+c(2)), 2a, 2b, 3/21, 4, 5, 6, 7a, 7b, 8, 9, 10a, 10b, 11, 12a/22a, 12b/22b, 13, 14/15, 16a, 16b, 17, 18, 19a, 19b, 20 and X and Y. From this analysis, we have been able to point out the human chromosome bands more "prone" to breakage during the evolutionary pathways and/or pathology processes. We have observed that 89.09% of the human chromosome bands, which are boundaries for ancestral chromosome segments, contain common fragile sites and/or intrachromosomal telomeric-like sequences. A more in depth analysis of twelve different human chromosomes has allowed us to determine that 62.16% of the chromosomal bands implicated in inversions and 100% involved in fusions/fissions correspond to fragile sites, intrachromosomal telomeric-like sequences and/or bands significantly affected by X irradiation. In addition, 73% of the bands affected in pathological processes are co-localized in bands where fragile sites, intrachromosomal telomeric-like sequences, bands significantly affected by X irradiation and/or evolutionary chromosomal bands have been described. Our data also support the hypothesis that chromosomal breakages detected in pathological processes are not randomly distributed along the chromosomes, but rather concentrate in those important evolutionary chromosome bands which correspond to fragile sites and/or intrachromosomal telomeric-like sequences.

Fluorescence in situ hybridization to chromosomes as a tool to understand human and primate genome evolution

For the last 15 years molecular cytogenetic techniques have been extensively used to study primate evolution. Molecular probes were helpful to distinguish mammalian chromosomes and chromosome segments on the basis of their DNA content rather than solely on morphological features such as banding patterns. Various landmark rearrangements have been identified for most of the nodes in primate phylogeny while chromosome banding still provides helpful reference maps. Fluorescence in situ hybridization (FISH) techniques were used with probes of different complexity including chromosome painting probes, probes derived from chromosome sub-regions and in the size of a single gene. Since more recently, in silico techniques have been applied to trace down evolutionarily derived chromosome rearrangements by searching the human and mouse genome sequence databases. More detailed breakpoint analyses of chromosome rearrangements that occurred during higher primate evolution also gave some insights into the molecular changes in chromosome rearrangements that occurred in evolution. Hardly any "fusion genes" as known from chromosome rearrangements in cancer cells or dramatic "position effects" of genes transferred to new sites in primate genomes have been reported yet. Most breakpoint regions have been identified within gene poor areas rich in repetitive elements and/or low copy repeats (segmental duplications). The progress in various molecular and molecular-cytogenetic approaches including the recently launched chimpanzee genome project suggests that these new tools will have a significant impact on the further understanding of human genome evolution.

Evolution of the tandem repeats in thymidylate synthase enhancer region (TSER) in primates

The upstream regulatory region of the human thymidylate synthase gene (thymidylate synthase enhancer region, TSER) is length polymorphic, attributable to variable numbers of tandemly repeated copies of a 28-bp fragment. It has been found that TSER length polymorphism is correlated to malignancy risk. To further our understanding of the origin and evolution of TSER, this region was investigated among different primates, including hominoids, two subfamilies of the Old World monkeys (OWMs): colobines and cercopithecines, and two species of the New World monkeys (NWMs). In addition to humans, our results show that length polymorphism in TSER is also present in some primate populations, although it appears that this region is length monomorphic in many other primates. We identified three unique repeat motifs in TSER and defined them as R1, R2, and R3, respectively, starting from the 3' end. The same repeat motifs from different species are more similar to each other than different repeat motifs within same species are. Such a paraphyletic pattern suggests that divergence of the three repeat motifs predated divergence of the OWMs/hominoids and the NWMs. The most recent common ancestor (MRCA) of hominoids and the OWMs probably possessed triple repeats but now double and triple repeats are two dominant types in hominoids and the OWMs. In addition, our results show that each of the three repeat motifs may be lost independently. We have also found clues that recombination was involved in formation of tandem repeat polymorphism in TSER.

A probe generated by chromosome microdissection, useful for analyzing Y chromosome evolution in Old World monkeys

We isolated a DNA probe, designated MMDYZ1, using a chromosome microdissection technique from the Y chromosome of the Rhesus monkey. The probe obtained from eight whole Y chromosomes shows higher specificity for the Y short arm of the Rhesus monkey, which consists totally of constitutive heterochromatin. Two microclones (MMY#3 and MMY#4) were constructed from the Y-specific primary PCR products. Sequence analysis of these two microclones revealed that both were essentially identical to each other and the sizes were 870 and 686 bp, respectively. From alignment analysis using the Genbank database of primates, the alphoid DNA has the highest affinity with the probe. However, the total composition of this probe has extremely high homology with the Y short arm of the Rhesus monkey, as demonstrated by fluorescence in-situ hybridization (FISH). Comparative FISH-mapping disclosed that this DNA-sequence cluster was located at extremely different sites on the Y chromosome in several species of the Old World monkey. Accordingly, this probe seems to be a high-quality tool, now established for the first time, for investigating Y chromosome evolution of the Old World monkey.

Timing the origin of New World monkeys

The origin of New World monkeys (Infraorder Platyrrhini) has been an extensively debated issue. In this study, we analyzed mitochondrial genomes from Cebus (Platyrrhini), Homo, Hylobates, Pan, Pongo (Hominoids), Macaca, Papio (Cercopithecoids), and Tarsius (outgroup) to investigate this matter. Two distinct methodologies were employed on mitochondrial genes to estimate divergence times: the traditional likelihood ratio test performed in ML analyses of individual and concatenated gene sequences and the recent multigene Bayesian approach. Using the Cercopithecoid-Hominoid split as calibration point (25 MYA), our results show consistently that Platyrrhines split from Catarrhines at around 35 MYA. Although the main focus of the study is New World monkey origins, we have also estimated other primate divergence times: Homo-Pan at 5-7 MYA; Pongo-(Homo/Pan) at 13-16 MYA; Hylobates-(Pongo/Homo/Pan) at 15-19 MYA; and Macaca-Papio at 10-12 MYA. Our estimate for the origin of New World monkeys is in agreement with the hypothesis of a transatlantic journey from Africa to South America, as suggested by the fossil record.

Neocentromeres in 15q24-26 map to duplicons which flanked an ancestral centromere in 15q25

The existence of latent centromeres has been proposed as a possible explanation for the ectopic emergence of neocentromeres in humans. This hypothesis predicts an association between the position of neocentromeres and the position of ancient centromeres inactivated during karyotypic evolution. Human chromosomal region 15q24-26 is one of several hotspots where multiple cases of neocentromere emergence have been reported, and it harbors a high density of chromosome-specific duplicons, rearrangements of which have been implicated as a susceptibility factor for panic and phobic disorders with joint laxity. We investigated the evolutionary history of this region in primates and found that it contains the site of an ancestral centromere which became inactivated about 25 million years ago, after great apes/Old World monkeys diverged. This inactivation has followed a noncentromeric chromosomal fission of an ancestral chromosome which gave rise to phylogenetic chromosomes XIV and XV in human and great apes. Detailed mapping of the ancient centromere and two neocentromeres in 15q24-26 has established that the neocentromere domains map approximately 8 Mb proximal and 1.5 Mb distal of the ancestral centromeric region, but that all three map within 500 kb of duplicons, copies of which flank the centromere in Old World Monkey species. This suggests that the association between neocentromere and ancestral centromere position on this chromosome may be due to the persistence of recombinogenic duplications accrued within the ancient pericentromere, rather than the retention of "centromere-competent" sequences per se. The high frequency of neocentromere emergence in the 15q24-26 region and the high density of clinically important duplicons are, therefore, understandable in the light of the evolutionary history of this region.