Primate origins: evolutionary change in digital ray patterning and segmentation

This study presents evidence that the first primates share with extant lemurs, tarsiers, and anthropoids hand proportions unlike those of their close relatives, the tree shrews (Scandentia), colugos (Dermoptera), and plesiadapiforms. Specifically, early primates as well as modern strepsirhines and haplorhines have relatively short metacarpals and long proximal phalanges giving them a grasping, prehensile hand. Limb development was studied in the primate Microcebus murinus and a comparative sample of rodents, artiodactyls, and marsupials to investigate the role of embryonic patterning in the morphogenesis and evolution of primate hand proportions. Comparative analysis shows that the derived finger proportions of primates are generated during the early phases of digital ray patterning and segmentation, when the interzone cells marking the presumptive metacarpo- and interphalangeal joints first appear. Interspecific variation in relative digit and metapodial proportions therefore has high developmental penetrance; that is, adult differences are observed at early ontogenetic stages. The paleontological, comparative, and developmental data are therefore consistent with the hypothesis that the early Cenozoic origin of primates involved an evolutionary change in digital ray pattern formation ultimately yielding a grasping, prehensile hand.     

Molecular evolution of cytochrome c oxidase subunit I in primates: is there coevolution between mitochondrial and nuclear genomes?

Phylogenetic analyses carried out on cytochrome c oxidase (COX) subunit I mitochondrial genes from 14 primates representing the major branches of the order and four outgroup nonprimate eutherians revealed that transversions and amino acid replacements (i.e., the more slowly occurring sequence changes) contained lower levels of homoplasy and thus provided more accurate information on cladistic relationships than transitions (i.e., the more rapidly occurring sequence changes). Several amino acids, each with a high likelihood of functionality involving the binding of cytochrome c or interaction with COX VIII, have changed in Anthropoidea, the primate suborder grouping New World monkey, Old World monkey, ape, and human lineages. They are conserved in other mammalian lineages and in nonanthropoid primates. Maximum-likelihood ancestral COX I nucleotide sequences were determined utilizing a near most parsimonious branching arrangement for the primate sequences that was consistent with previously hypothesized primate cladistic relationships based on larger and more diverse data sets. Relative rate tests of COX I mitochondrial sequences showed an elevated nonsynonymous (N) substitution rate for anthropoid-nonanthropoid comparisons. This finding for the largest mitochondrial (mt) DNA-encoded subunit is consistent with previous observations of elevated nonsynonymous substitution/synonymous substitution (S) rates in primates for mt-encoded COX II and for the nuclear-encoded COX IV and COX VIIa-H. Other COX-related proteins, including cytochrome c and cytochrome b, also show elevated amino acid replacement rates or N/S during similar time frames, suggesting that this group of interacting genes is likely to have coevolved during primate evolution.     

Molecular phylogeny and evolution of prosimians based on complete sequences of mitochondrial DNAs

Prosimians (tarsiers and strepsirrhini) represent the basal lineages in primates and have a close bearing on the origin of primates. Although major lineages among anthropoidea (humans, apes and monkeys) are well represented by complete mitochondrial DNA (mtDNA) sequence data, only one complete mtDNA sequence from a representative of each of the infraorders in prosimians has been described until quite recently, and therefore we newly determined complete mtDNA sequences from 5 lemurs, 4 lorises, one tarsier and one platyrrhini. These sequences were provided to phylogenetic analyses in combination with the sequences from the 15 primates species reported to the database. The position of tarsiers among primates could not be resolved by the maximum likelihood (ML) and neighbor-joining (NJ) analyses with several data sets. As to the position of tarsiers, any of the three alternative topologies (monophyly of haplorhini, monophyly of prosimians, and tarsiers being basal in primates) was not rejected at the significance level of 5%, neither at the nucleotide nor at the amino acid level. In addition, the significant variations of C and T compositions were observed across primates species. Furthermore, we used AGY data sets for phylogenetic analyses in order to remove the effect of different C/T composition bias across species. The analyses of AGY data sets provided a medium support for the monophyly of haplorhini, which might have been screened by the variation in base composition of mtDNA across species. To estimates the speciation dates within primates, we analyzed the amino acid sequences of mt-proteins with a Bayesian method of Thorne and Kishino. Divergence dates were estimated as follows for the crown groups: about 35.4 million years ago (mya) for lorisiformes, 55.3 mya for lemuriformes, 64.5 mya for strepsirrhini, 70.1 mya for haplorhini and 76.0 mya for primates. Furthermore, we reexamined the biogeographic scenarios which have been proposed for the origin of strepsirrhini (lemuriformes and lorisiformes) and for the dispersal of the lemuriformes and lorisiformes.     

Allelic lineages of the ficolin genes (FCNs) are passed from ancestral to descendant primates

The ficolins recognize carbohydrates and acetylated compounds on microorganisms and dying host cells and are able to activate the lectin pathway of the complement system. In humans, three ficolin genes have been identified: FCN1, FCN2 and FCN3, which encode ficolin-1, ficolin-2 and ficolin-3, respectively. Rodents have only two ficolins designated ficolin-A and ficolin-B that are closely related to human ficolin-1, while the rodent FCN3 orthologue is a pseudogene. Ficolin-2 and ficolin-3 have so far only been observed in humans. Thus, we performed a systematic investigation of the FCN genes in non-human primates. The exons and intron-exon boundaries of the FCN1-3 genes were sequenced in the following primate species: chimpanzee, gorilla, orangutan, rhesus macaque, cynomolgus macaque, baboon and common marmoset. We found that the exon organisation of the FCN genes was very similar between all the non-human primates and the human FCN genes. Several variations in the FCN genes were found in more than one primate specie suggesting that they were carried from one species to another including humans. The amino acid diversity of the ficolins among human and non-human primate species was estimated by calculating the Shannon entropy revealing that all three proteins are generally highly conserved. Ficolin-1 and ficolin-2 showed the highest diversity, whereas ficolin-3 was more conserved. Ficolin-2 and ficolin-3 were present in non-human primate sera with the same characteristic oligomeric structures as seen in human serum. Taken together all the FCN genes show the same characteristics in lower and higher primates. The existence of trans-species polymorphisms suggests that different FCN allelic lineages may be passed from ancestral to descendant species.     

Parasite-mediated evolution of the functional part of the MHC in primates

The major histocompatibility complex (MHC) is a key model of genetic polymorphism, but the mechanisms underlying its extreme variability are debated. Most hypotheses for MHC diversity focus on pathogen-driven selection and predict that MHC polymorphism evolves under the pressure of a diverse parasite fauna. Several studies reported that certain alleles offer protection against certain parasites, yet it remains unclear whether variation in parasite pressure more generally covaries with allelic diversity and rates of molecular evolution of MHC across species. We tested this prediction in a comparative study of 41 primate species. We characterized polymorphism of the exon 2 of DRB region of the MHC class II. Our phylogenetic analyses controlled for the potential effects of neutral mutation rate, population size, geographic origin and body mass and revealed that nematode species richness associates positively with nonsynonymous nucleotide substitution rate at the functional part of the molecule. We failed to find evidence for allelic diversity being strongly related to parasite species richness. Continental distribution was a strong predictor of both allelic diversity and substitution rate, with higher values in Malagasy and Neotropical primates. These results indicate that parasite pressure can influence the different estimates of MHC polymorphism, whereas geography plays an independent role in the natural history of MHC.     

Different selective pressures shape the molecular evolution of color vision in chimpanzee and human populations

A population genetic analysis of the long-wavelength opsin (OPN1LW, "red") color vision gene in a global sample of 236 human nucleotide sequences had previously discovered nine amino acid replacement single nucleotide polymorphisms, which were found at high frequencies in both African and non-African populations and associated with an unusual haplotype diversity. Although this pattern of nucleotide diversity is consistent with balancing selection, it has been argued that a recombination "hot spot" or gene conversion within and between X-linked color vision genes alone may explain these patterns. The current analysis investigates a closely related primate with trichromatism to determine whether color vision gene amino acid polymorphism and signatures of adaptive evolution are characteristic of humans alone. Our population sample of 56 chimpanzee (Pan troglodytes) OPN1LW sequences shows three singleton amino acid polymorphisms and no unusual recombination or linkage disequilibrium patterns across the approximately 5.5-kb region analyzed. Our comparative population genetic approach shows that the patterns of OPN1LW variation in humans and chimpanzees are consistent with positive and purifying selection within the two lineages, respectively. Although the complex role of color vision has been greatly documented in primate evolution in general, it is surprising that trichromatism has followed very different selective trajectories even between humans and our closest relatives.     

Control Region Length Dynamics Potentially Drives Amino Acid Evolution in Tarsier Mitochondrial Genomes

Patterns and processes of molecular evolution critically influence inferences in phylogeny and phylogeography. Within primates, a shift in evolutionary rates has been identified as the rationale for contrasting findings from mitochondrial and nuclear DNA studies as to the position of Tarsius. While the latter now seems settled, we sequenced complete mitochondrial genomes of three Sulawesi tarsiers (Tarsius dentatus, T. lariang, and T. wallacei) and analyzed substitution rates among tarsiers and other primates to infer driving processes of molecular evolution. We found substantial length polymorphism of the D-loop within tarsier individuals, but little variation of predominant lengths among them, regardless of species. Length variation was due to repetitive elements in the CSB domain—minisatellite motifs of 35 bp length and microsatellite motifs of 6 bp length. Amino acid evolutionary rates were second highest among major primate taxa relative to nucleotide substitution rates. We observed many radical possibly function-altering amino acid changes that were rarely driven by positive selection and thus potentially slightly deleterious or neutral. We hypothesize that the observed pattern of an increased amino acid evolutionary rate in tarsier mitochondrial genomes may be caused by hitchhiking of slightly deleterious mutations with favored D-loop length variants selected for maximizing replication success within the cell or the mitochondrion.     

Comparative genomic analysis of the false killer whale (Pseudorca crassidens) LMBR1 locus

The sequencing and comparative genomic analysis of LMBR1 loci in mammals or other species, including human, would be very important in understanding evolutionary genetic changes underlying the evolution of limb development. In this regard, comparative genomic annotation of the false killer whale LMBR1 locus could shed new light on the evolution of limb development. We sequenced two false killer whale BAC clones, corresponding to 156 kb and 144 kb, respectively, harboring the tightly linked RNF32, LMBR1, and NOM1 genes. Our annotation of the false killer whale LMBR1 gene showed that it consists of 17 exons (1473 bp), in contrast to 18 exons (1596 bp) in human, and it displays 93.1% and 95.6% nucleotide and amino acid sequence similarity, respectively, compared with the human gene. In particular, we discovered that exon 10, deleted in the false killer whale LMBR1 gene, is present only in primates, and this fact strongly implies that exon 10 might be crucial in determining primate-specific limb development. ZRS and TFBS sequences have been well conserved across 11 species, suggesting that these regions could be involved in an important function of limb development and limb patterning. The neighboring gene RNF32 showed several lineage-conserved exons, such as exons 2 through 9 conserved in eutherian mammals, exons 3 through 9 conserved in mammals, and exons 5 through 9 conserved in vertebrates. The other neighboring gene, NOM1, had undergone a substitution (ATG→GTA) at the start codon, giving rise to a 36 bp shorter N-terminal sequence compared with the human sequence. Our comparative analysis of the false killer whale LMBR1 genomic locus provides important clues regarding the genetic regions that may play crucial roles in limb development and patterning.     

Dynamics of DNA Methylation in Recent Human and Great Ape Evolution

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

Phylogenetic Analysis and Selection Pressures of 5-HT Receptors in Human and Non-human Primates: Receptor of an Ancient Neurotransmitter

Abstract

Neurotransmitter serotonin (5-hydroxytryptamine, 5-HT) an ancient neurotransmitter, involved in several neurophysiological and behavioral functions, acts by interacting with multiple receptors (5-HT₁-5-HT₇). Alterations in serotonergic signalling have also been implicated in various psychiatric disorders. The availability of the genome data of nonhuman primates permits comparative analysis of human 5-HT receptors with sequences of non-human primates to understand evolutionary divergence. We compared and analyzed serotonergic receptor sequences from human and non-human primates. Phylogenetic analysis by Maximum Likelihood (ML) method classified human and primate 5-HT receptors into six unique clusters. There was considerable conservation of 5-HT receptor sequences between human and non-human primates; however, a greater diversity at the sub-group level was observed. Compared to the other subgroups, larger multiplicity and expansion was seen within the 5-HT₄ receptor subtype in both human and non-human primates. Analysis of non-synonymous and synonymous substitution ratios (Ka/Ks ratio) using the Nei-Gojobori method suggests that 5-HT receptor sequences have undergone negative (purifying) selection over the course of evolution in human, chimpanzee and rhesus monkey. Abnormal human and non-human primate psychopathalogy and behavior, in the context of these variations is discussed. Analysis of these 5-HT receptors in other species will help understand the molecular evolution of 5-HT receptors, and its possible influence on complex behaviors, and psychiatric disorders.     

Genomic data reject the hypothesis of a prosimian primate clade

The phylogenetic position of tarsiers within the primates has been a controversial subject for over a century. Despite numerous morphological and molecular studies, there has been weak support for grouping tarsiers with either strepsirrhine primates in a prosimian clade or with anthropoids in a haplorrhine clade. Here, we take advantage of the recently released whole genome assembly of the Philippine tarsier, Tarsius syrichta, in order to infer the phylogenetic relationship of Tarsius within the order Primates. We also present estimates of divergence times within the primates. Using a 1.26 million base pair multiple sequence alignment derived from 1078 orthologous genes, we provide overwhelming statistical support for the presence of a haplorrhine clade. We also present divergence date estimates using local relaxed molecular clock methods. The estimated time of the most recent common ancestor of extant Primates ranged from 64.9 Ma to 72.6 Ma, and haplorrhines were estimated to have a most recent common ancestor between 58.9 Ma and 68.6 Ma. Examination of rates of nucleotide substitution in the three major extant primate clades show that anthropoids have a slower substitution rate than either strepsirrhines or tarsiers. Our results provide the framework on which primate morphological, reproductive, and genomic features can be reconstructed in the broader context of mammalian phylogeny.     

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.     

Brain reorganization,not relative brain size,primarily characterizes anthropoid brain evolution

Comparative analyses of primate brain evolution have highlighted changes in size and internal organization as key factors underlying species diversity. It remains, however, unclear (i) how much variation in mosaic brain reorganization versus variation in relative brain size contributes to explaining the structural neural diversity observed across species, (ii) which mosaic changes contribute most to explaining diversity, and (iii) what the temporal origin, rates and processes are that underlie evolutionary shifts in mosaic reorganization for individual branches of the primate tree of life. We address these questions by combining novel comparative methods that allow assessing the temporal origin, rate and process of evolutionary changes on individual branches of the tree of life, with newly available data on volumes of key brain structures (prefrontal cortex, frontal motor areas and cerebrocerebellum) for a sample of 17 species (including humans). We identify patterns of mosaic change in brain evolution that mirror brain systems previously identified by electrophysiological and anatomical tract-tracing studies in non-human primates and functional connectivity MRI studies in humans. Across more than 40 Myr of anthropoid primate evolution, mosaic changes contribute more to explaining neural diversity than changes in relative brain size, and different mosaic patterns are differentially selected for when brains increase or decrease in size. We identify lineage-specific evolutionary specializations for all branches of the tree of life covered by our sample and demonstrate deep evolutionary roots for mosaic patterns associated with motor control and learning.     

Primate home range and GRIN2A, a receptor gene involved in neuronal plasticity: implications for the evolution of spatial memory

N -methyl- d -aspartate (NMDA) glutamate receptors play crucial roles in neuronal synaptic plasticity, learning and memory. However, as to whether different NMDA subunits are implicated in specific forms of memory is unclear. Moreover, nothing is known about the interspecific genetic variability of the GRIN2A subunit and how this variation can potentially explain evolutionary changes in behavioral phenotypes. Here, we used 28 primate GRIN2A sequences and various proxies of memory across primates to investigate the role of GRIN2A . Codon-specific sequence analysis on these sequences showed that GRIN2A in primates coevolved with a likely ecological proxy of spatial memory (relative home-range size) but not with other indices of non-spatial learning and memory such as social memory and social learning. Models based on gene averages failed to detect positive selection in primate branches with major changes in relative home-range size. This implies that accelerated evolution is concentrated in specific parts of the protein expressed by GRIN2A . Overall, our molecular evolution study, the first on GRIN2A , supports the notion that different NMDA subunits may play a role in specific forms of memory and that phenotypic diversity along with genetic evolution can be used to investigate the link between genes and behavior across evolutionary time.     

Molecular evolution of a primate-specific microRNA family

Lineage-specific microRNA (miRNA) families may contribute to developmental novelties during evolution. However, little is known about the origin and evolution of new miRNA families. We report evidence of an Alu-mediated rapid expansion of miRNA genes in a previously identified primate-specific miRNA family, drawn from sequencing and comparative analysis of 9 diverse primate species. Evolutionary analysis reveals similar divergence among miRNA copies whether they are within or between species, lineage-specific gain and loss of miRNAs, and gene pseudolization in multiple species. These observations support a birth-and-death process of miRNA genes in this family, implicating functional diversification during primate evolution. In addition, both secondary structure conservation and reduced single nucleotide polymorphisms density attest to functional constraint of this family in primates. Finally, we observed preferential expression of miRNAs in human placenta and fetal brain, suggesting a functional importance of this family for primate development.     

Characterization of Novel MSX1 Mutations Identified in Japanese Patients with Nonsyndromic Tooth Agenesis

Since MSX1 and PAX9 are linked to the pathogenesis of nonsyndromic tooth agenesis, we performed detailed mutational analysis of these two genes sampled from Japanese patients. We identified two novel MSX1 variants with an amino acid substitution within the homeodomain; Thr174Ile (T174I) from a sporadic hypodontia case and Leu205Arg (L205R) from a familial oligodontia case. Both the Thr174 and Leu205 residues in the MSX1 homeodomain are highly conserved among different species. To define possible roles of mutations at these amino acids in the pathogenesis of nonsyndromic tooth agenesis, we performed several functional analyses. It has been demonstrated that MSX1 plays a pivotal role in hard tissue development as a suppressor for mesenchymal cell differentiation. To evaluate the suppression activity of the variants in mesenchymal cells, we used the myoD-promoter, which is one of convenient reporter assay system for MSX1. Although the gene products of these MSX1 variants are stable and capable of normal nuclear localization, they do not suppress myoD-promoter activity in differentiated C2C12 cells. To clarify the molecular mechanisms underlying our results, we performed further analyses including electrophoretic mobility shift assays, and co-immunoprecipitation assays to survey the molecular interactions between the mutant MSX1 proteins and the oligonucleotide DNA with MSX1 consensus binding motif or EZH2 methyltransferase. Since EZH2 is reported to interact with MSX1 and regulate MSX1 mediated gene suppression, we hypothesized that the T174I and L205R substitutions would impair this interaction. We conclude from the results of our experiments that the DNA binding ability of MSX1 is abolished by these two amino acid substitutions. This illustrates a causative role of the T174I and L205R MSX1 homeodomain mutations in tooth agenesis, and suggests that they may influence cell proliferation and differentiation resulting in lesser tooth germ formation in vivo.