The complete sequence and comparative analysis of ape sex chromosomes.

Apes possess two sex chromosomes-the male-specific Y chromosome and the X chromosome, which is present in both males and females. The Y chromosome is crucial for male reproduction, with deletions being linked to infertility1. The X chromosome is vital for reproduction and cognition2. Variation in mating patterns and brain function among apes suggests corresponding differences in their sex chromosomes. However, owing to their repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the methodology developed for the telomere-to-telomere (T2T) human genome, we produced gapless assemblies of the X and Y chromosomes for five great apes (bonobo (Pan paniscus), chimpanzee (Pan troglodytes), western lowland gorilla (Gorilla gorilla gorilla), Bornean orangutan (Pongo pygmaeus) and Sumatran orangutan (Pongo abelii)) and a lesser ape (the siamang gibbon (Symphalangus syndactylus)), and untangled the intricacies of their evolution. Compared with the X chromosomes, the ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements-owing to the accumulation of lineage-specific ampliconic regions, palindromes, transposable elements and satellites. Many Y chromosome genes expand in multi-copy families and some evolve under purifying selection. Thus, the Y chromosome exhibits dynamic evolution, whereas the X chromosome is more stable. Mapping short-read sequencing data to these assemblies revealed diversity and selection patterns on sex chromosomes of more than 100 individual great apes. These reference assemblies are expected to inform human evolution and conservation genetics of non-human apes, all of which are endangered species.

Understanding Language Evolution: Beyond Pan-Centrism.

Language does not fossilize but this does not mean that the language's evolutionary timeline is lost forever. Great apes provide a window back in time on our last prelinguistic ancestor's communication and cognition. Phylogeny and cladistics implicitly conjure Pan (chimpanzees, bonobos) as a superior (often the only) model for language evolution compared with earlier diverging lineages, Gorilla and Pongo (orangutans). Here, in reviewing the literature, it is shown that Pan do not surpass other great apes along genetic, cognitive, ecologic, or vocal traits that are putatively paramount for language onset and evolution. Instead, revived herein is the idea that only by abandoning single-species models and learning about the variation among great apes, there might be a chance to retrieve lost fragments of the evolutionary timeline of language.

Evolutionary history of orangutan plasmodia revealed by phylogenetic analysis of complete mtDNA genomes and new biogeographical divergence dating calibration models.

During the past 15 years, researchers have shown a renewed interest in the study of the Plasmodium parasites that infect orangutans. Most recently, studies examined the phylogenetic relationships and divergence dates of these parasites in orangutans using complete mitochondrial DNA genomes. Questions regarding the dating of these parasites, however, remain. In the present study, we provide a new calibration model for dating the origins of Plasmodium parasites in orangutans using a modified date range for the origin of macaques in Asia. Our Bayesian phylogenetic analyses of complete Plasmodium sp. mitochondrial DNA genomes inferred two clades of plasmodia in orangutans (Pongo 1 and Pongo 2), and that these clades likely represent the previously identified species Plasmodium pitheci and Plasmodium silvaticum. However, we cannot identify which Pongo clade is representative of the morphologically described species. The most recent common ancestor of both Pongo sp. plasmodia, Plasmodium. hylobati, and Plasmodium. inui dates to 3-3.16 million years ago (mya) (95% highest posterior density [HPD]: 2.09-4.08 mya). The Pongo 1 parasite diversified 0.33-0.36 mya (95% HPD: 0.12-0.63), while the Pongo 2 parasite diversified 1.15-1.22 mya (95% HPD: 0.63-1.82 mya). It now seems likely that the monkey Plasmodium (P. inui) is the result of a host switch event from the Pongo 2 parasite to sympatric monkeys, or P. hylobati. Our new estimates for the divergence of orangutan malaria parasites, and subsequent diversification, are all several hundred thousand years later than previous Bayesian estimates.

Nomenclature for the KIR of non-human species.

The increasing number of Killer Immunoglobulin-like Receptor (KIR) sequences available for non-human primate species and cattle has prompted development of a centralized database, guidelines for a standardized nomenclature, and minimum requirements for database submission. The guidelines and nomenclature are based on those used for human KIR and incorporate modifications made for inclusion of non-human species in the companion IPD-NHKIR database. Included in this first release are the rhesus macaque (Macaca mulatta), chimpanzee (Pan troglodytes), orangutan (Pongo abelii and Pongo pygmaeus), and cattle (Bos taurus).

Complex MHC Class I Gene Transcription Profiles and Their Functional Impact in Orangutans.

MHC haplotypes of humans and the African great ape species have one copy of the MHC-A, -B, and -C genes. In contrast, MHC haplotypes of orangutans, the Asian great ape species, exhibit variation in the number of gene copies. An in-depth analysis of the MHC class I gene repertoire in the two orangutan species, Pongo abelii and Pongo pygmaeus, is presented in this article. This analysis involved Sanger and next-generation sequencing methodologies, revealing diverse and complicated transcription profiles for orangutan MHC-A, -B, and -C. Thirty-five previously unreported MHC class I alleles are described. The data demonstrate that each orangutan MHC haplotype has one copy of the MHC-A gene, and that the MHC-B region has been subject to duplication, giving rise to at least three MHC-B genes. The MHC-B*03 and -B*08 lineages of alleles each account for a separate MHC-B gene. All MHC-B*08 allotypes have the C1-epitope motif recognized by killer cell Ig-like receptor. At least one other MHC-B gene is present, pointing to MHC-B alleles that are not B*03 or B*08. The MHC-C gene is present only on some haplotypes, and each MHC-C allotype has the C1-epitope. The transcription profiles demonstrate that MHC-A alleles are highly transcribed, whereas MHC-C alleles, when present, are transcribed at very low levels. The MHC-B alleles are transcribed to a variable extent and over a wide range. For those orangutan MHC class I allotypes that are detected by human monoclonal anti-HLA class I Abs, the level of cell-surface expression of proteins correlates with the level of transcription of the allele.

Comparative and demographic analysis of orang-utan genomes.

'Orang-utan' is derived from a Malay term meaning 'man of the forest' and aptly describes the southeast Asian great apes native to Sumatra and Borneo. The orang-utan species, Pongo abelii (Sumatran) and Pongo pygmaeus (Bornean), are the most phylogenetically distant great apes from humans, thereby providing an informative perspective on hominid evolution. Here we present a Sumatran orang-utan draft genome assembly and short read sequence data from five Sumatran and five Bornean orang-utan genomes. Our analyses reveal that, compared to other primates, the orang-utan genome has many unique features. Structural evolution of the orang-utan genome has proceeded much more slowly than other great apes, evidenced by fewer rearrangements, less segmental duplication, a lower rate of gene family turnover and surprisingly quiescent Alu repeats, which have played a major role in restructuring other primate genomes. We also describe a primate polymorphic neocentromere, found in both Pongo species, emphasizing the gradual evolution of orang-utan genome structure. Orang-utans have extremely low energy usage for a eutherian mammal, far lower than their hominid relatives. Adding their genome to the repertoire of sequenced primates illuminates new signals of positive selection in several pathways including glycolipid metabolism. From the population perspective, both Pongo species are deeply diverse; however, Sumatran individuals possess greater diversity than their Bornean counterparts, and more species-specific variation. Our estimate of Bornean/Sumatran speciation time, 400,000 years ago, is more recent than most previous studies and underscores the complexity of the orang-utan speciation process. Despite a smaller modern census population size, the Sumatran effective population size (N(e)) expanded exponentially relative to the ancestral N(e) after the split, while Bornean N(e) declined over the same period. Overall, the resources and analyses presented here offer new opportunities in evolutionary genomics, insights into hominid biology, and an extensive database of variation for conservation efforts.

Evolutionary toggling of the MAPT 17q21.31 inversion region.

Using comparative sequencing approaches, we investigated the evolutionary history of the European-enriched 17q21.31 MAPT inversion polymorphism. We present a detailed, BAC-based sequence assembly of the inverted human H2 haplotype and compare it to the sequence structure and genetic variation of the corresponding 1.5-Mb region for the noninverted H1 human haplotype and that of chimpanzee and orangutan. We found that inversion of the MAPT region is similarly polymorphic in other great ape species, and we present evidence that the inversions occurred independently in chimpanzees and humans. In humans, the inversion breakpoints correspond to core duplications with the LRRC37 gene family. Our analysis favors the H2 configuration and sequence haplotype as the likely great ape and human ancestral state, with inversion recurrences during primate evolution. We show that the H2 architecture has evolved more extensive sequence homology, perhaps explaining its tendency to undergo microdeletion associated with mental retardation in European populations.

Reconstructing the demographic history of orang-utans using Approximate Bayesian Computation.

Investigating how different evolutionary forces have shaped patterns of DNA variation within and among species requires detailed knowledge of their demographic history. Orang-utans, whose distribution is currently restricted to the South-East Asian islands of Borneo (Pongo pygmaeus) and Sumatra (Pongo abelii), have likely experienced a complex demographic history, influenced by recurrent changes in climate and sea levels, volcanic activities and anthropogenic pressures. Using the most extensive sample set of wild orang-utans to date, we employed an Approximate Bayesian Computation (ABC) approach to test the fit of 12 different demographic scenarios to the observed patterns of variation in autosomal, X-chromosomal, mitochondrial and Y-chromosomal markers. In the best-fitting model, Sumatran orang-utans exhibit a deep split of populations north and south of Lake Toba, probably caused by multiple eruptions of the Toba volcano. In addition, we found signals for a strong decline in all Sumatran populations ~24 ka, probably associated with hunting by human colonizers. In contrast, Bornean orang-utans experienced a severe bottleneck ~135 ka, followed by a population expansion and substructuring starting ~82 ka, which we link to an expansion from a glacial refugium. We showed that orang-utans went through drastic changes in population size and connectedness, caused by recurrent contraction and expansion of rainforest habitat during Pleistocene glaciations and probably hunting by early humans. Our findings emphasize the fact that important aspects of the evolutionary past of species with complex demographic histories might remain obscured when applying overly simplified models.

Estimating divergence time and ancestral effective population size of Bornean and Sumatran orangutan subspecies using a coalescent hidden Markov model.

Due to genetic variation in the ancestor of two populations or two species, the divergence time for DNA sequences from two populations is variable along the genome. Within genomic segments all bases will share the same divergence-because they share a most recent common ancestor-when no recombination event has occurred to split them apart. The size of these segments of constant divergence depends on the recombination rate, but also on the speciation time, the effective population size of the ancestral population, as well as demographic effects and selection. Thus, inference of these parameters may be possible if we can decode the divergence times along a genomic alignment. Here, we present a new hidden Markov model that infers the changing divergence (coalescence) times along the genome alignment using a coalescent framework, in order to estimate the speciation time, the recombination rate, and the ancestral effective population size. The model is efficient enough to allow inference on whole-genome data sets. We first investigate the power and consistency of the model with coalescent simulations and then apply it to the whole-genome sequences of the two orangutan sub-species, Bornean (P. p. pygmaeus) and Sumatran (P. p. abelii) orangutans from the Orangutan Genome Project. We estimate the speciation time between the two sub-species to be thousand years ago and the effective population size of the ancestral orangutan species to be , consistent with recent results based on smaller data sets. We also report a negative correlation between chromosome size and ancestral effective population size, which we interpret as a signature of recombination increasing the efficacy of selection.

The interplay between genome organization and nuclear architecture of primate evolutionary neo-centromeres.

An Evolutionary Neo-Centromere (ENC) is a centromere that emerged in an ectopic region of a chromosome during evolution. It is thought that the old centromere must be inactivated because dicentric chromosomes are not viable. The aim of the present study was to investigate whether 3D arrangement in the interphase nucleus of the novel and old centromeric domains was affected by the repositioning event. The data we present here strongly indicate that the ENC phenomenon does not affect the 3D location of either novel or old centromeres. Very likely, other features, such as gene density, rather than the newly acquired or lost functions, define positioning in the nucleus.

Impact of replication timing on non-CpG and CpG substitution rates in mammalian genomes.

Neutral nucleotide substitutions occur at varying rates along genomes, and it remains a major issue to unravel the mechanisms that cause these variations and to analyze their evolutionary consequences. Here, we study the role of replication in the neutral substitution pattern. We obtained a high-resolution replication timing profile of the whole human genome by massively parallel sequencing of nascent BrdU-labeled replicating DNA. These data were compared to the neutral substitution rates along the human genome, obtained by aligning human and chimpanzee genomes using macaque and orangutan as outgroups. All substitution rates increase monotonously with replication timing even after controlling for local or regional nucleotide composition, crossover rate, distance to telomeres, and chromatin compaction. The increase in non-CpG substitution rates might result from several mechanisms including the increase in mutation-prone activities or the decrease in efficiency of DNA repair during the S phase. In contrast, the rate of C --> T transitions in CpG dinucleotides increases in later-replicating regions due to increasing DNA methylation level that reflects a negative correlation between timing and gene expression. Similar results are observed in the mouse, which indicates that replication timing is a main factor affecting nucleotide substitution dynamics at non-CpG sites and constitutes a major neutral process driving mammalian genome evolution.

Genome assembly quality: assessment and improvement using the neutral indel model.

We describe a statistical and comparative-genomic approach for quantifying error rates of genome sequence assemblies. The method exploits not substitutions but the pattern of insertions and deletions (indels) in genome-scale alignments for closely related species. Using two- or three-way alignments, the approach estimates the amount of aligned sequence containing clusters of nucleotides that were wrongly inserted or deleted during sequencing or assembly. Thus, the method is well-suited to assessing fine-scale sequence quality within single assemblies, between different assemblies of a single set of reads, and between genome assemblies for different species. When applying this approach to four primate genome assemblies, we found that average gap error rates per base varied considerably, by up to sixfold. As expected, bacterial artificial chromosome (BAC) sequences contained lower, but still substantial, predicted numbers of errors, arguing for caution in regarding BACs as the epitome of genome fidelity. We then mapped short reads, at approximately 10-fold statistical coverage, from a Bornean orangutan onto the Sumatran orangutan genome assembly originally constructed from capillary reads. This resulted in a reduced gap error rate and a separation of error-prone from high-fidelity sequence. Over 5000 predicted indel errors in protein-coding sequence were corrected in a hybrid assembly. Our approach contributes a new fine-scale quality metric for assemblies that should facilitate development of improved genome sequencing and assembly strategies.

Effects of Pleistocene glaciations and rivers on the population structure of Bornean orangutans (Pongo pygmaeus).

Sundaland, a tropical hotspot of biodiversity comprising Borneo and Sumatra among other islands, the Malay Peninsula, and a shallow sea, has been subject to dramatic environmental processes. Thus, it presents an ideal opportunity to investigate the role of environmental mechanisms in shaping species distribution and diversity. We investigated the population structure and underlying mechanisms of an insular endemic, the Bornean orangutan (Pongo pygmaeus). Phylogenetic reconstructions based on mtDNA sequences from 211 wild orangutans covering the entire range of the species indicate an unexpectedly recent common ancestor of Bornean orangutans 176 ka (95% highest posterior density, 72-322 ka), pointing to a Pleistocene refugium. High mtDNA differentiation among populations and rare haplotype sharing is consistent with a pattern of strong female philopatry. This is corroborated by isolation by distance tests, which show a significant correlation between mtDNA divergence and distance and a strong effect of rivers as barriers for female movement. Both frequency-based and Bayesian clustering analyses using as many as 25 nuclear microsatellite loci revealed a significant separation among all populations, as well as a small degree of male-mediated gene flow. This study highlights the unique effects of environmental and biological features on the evolutionary history of Bornean orangutans, a highly endangered species particularly vulnerable to future climate and anthropogenic change as an insular endemic.

Analysis of the beta-globin gene in DNA of suspected thalassemic great apes.

DNA was recovered from teeth of 2 great ape skeletons, Pan troglodytes (Ptr) and Pongo pygmaeus (Ppy), belonging to a 19th-century zoological collection. The skeletons presented morphological alterations possibly associated with ß-thalassemia: Ptr had deformation of the calvaria and oro-maxillo-facial bones with porotic hyperostosis and extended osteoporotic lesions of the skeleton, while Ppy showed a general marked widening of the calvarial diploe but moderate osteoporotic signs on the post-cranial skeleton. We screened Ptr and Ppy for mutations in the ß-globin gene (exons 1, 2, and 3) because we suspected thalassemia. Ptr ß-globin sequences showed the highest degree of similarity with the human ones (99.8%), while those of Ppy were slightly different (98.2%). The results were consistent with the phylogenetic relationships between their β-globin gene sequences. We did not find any mutation in the ß-globin gene of Ptr and Ppy; therefore, we conclude that, in spite of skeletal alterations, the 2 subjects analyzed were not affected by ß-thalassemia.

A subterminal satellite located adjacent to telomeres in chimpanzees is absent from the human genome.

One of the significant unresolved differences between the karyotypes of humans and African apes is the presence of positively staining G-bands at the ends of many chromosome arms in the chimpanzee and gorilla but absent from human chromosomes. Using a telomere anchored PCR strategy, we have isolated DNA from a subterminal satellite, composed of a 32 basepair A-T rich repeat, from the chimpanzee genome that hybridizes to all the additional terminal bands and at two interstitial sites. The satellite is more abundant in gorillas and is not detected in humans or orangutans. Furthermore, there is no similarity between other chimpanzee telomere-junction clones and human subterminal sequences, and therefore the organization of sequences adjacent to telomeres is very different between these closely related primates.

Great ape DNA sequences reveal a reduced diversity and an expansion in humans.

The extent of DNA sequence variation of chimpanzees is several-fold greater than that of humans. It is unclear, however, if humans or chimpanzees are exceptional among primates in having low and high amounts of DNA sequence diversity, respectively. To address this, we have determined approximately 10,000 bp of noncoding DNA sequences at Xq13.3 (which has been extensively studied in both humans and chimpanzees) from 10 western lowland gorillas (Gorilla gorilla gorilla) and 1 mountain gorilla (Gorilla gorilla beringei; that is, from 2 of the 3 currently recognized gorilla subspecies), as well as 8 Bornean (Pongo pygmaeus pygmaeus) and 6 Sumatran (Pongo pygmaeus abelii) orang-utans, representing both currently recognized orang-utan subspecies. We show that humans differ from the great apes in having a low level of genetic variation and a signal of population expansion.

Sex-biased dispersal and volcanic activities shaped phylogeographic patterns of extant Orangutans (genus: Pongo).

The Southeast Asian Sunda archipelago harbors a rich biodiversity with a substantial proportion of endemic species. The evolutionary history of these species has been drastically influenced by environmental forces, such as fluctuating sea levels, climatic changes, and severe volcanic activities. Orangutans (genus: Pongo), the only Asian great apes, are well suited to study the relative impact of these forces due to their well-documented behavioral ecology, strict habitat requirements, and exceptionally slow life history. We investigated the phylogeographic patterns and evolutionary history of orangutans in the light of the complex geological and climatic history of the Sunda archipelago. Our study is based on the most extensive genetic sampling to date, covering the entire range of extant orangutan populations. Using data from three mitochondrial DNA (mtDNA) genes from 112 wild orangutans, we show that Sumatran orangutans, Pongo abelii, are paraphyletic with respect to Bornean orangutans (P. pygmaeus), the only other currently recognized species within this genus. The deepest split in the mtDNA phylogeny of orangutans occurs across the Toba caldera in northern Sumatra and, not as expected, between both islands. Until the recent past, the Toba region has experienced extensive volcanic activity, which has shaped the current phylogeographic patterns. Like their Bornean counterparts, Sumatran orangutans exhibit a strong, yet previously undocumented structuring into four geographical clusters. However, with 3.50 Ma, the Sumatran haplotypes have a much older coalescence than their Bornean counterparts (178 kya). In sharp contrast to the mtDNA data, 18 Y-chromosomal polymorphisms show a much more recent coalescence within Sumatra compared with Borneo. Moreover, the deep geographic structure evident in mtDNA is not reflected in the male population history, strongly suggesting male-biased dispersal. We conclude that volcanic activities have played an important role in the evolutionary history of orangutans and potentially of many other forest-dwelling Sundaland species. Furthermore, we demonstrate that a strong sex bias in dispersal can lead to conflicting patterns in uniparentally inherited markers even at a genus-wide scale, highlighting the need for a combined usage of maternally and paternally inherited marker systems in phylogenetic studies.

The adaptive significance of unproductive alternative splicing in primates.

Alternative gene splicing is pervasive in metazoa, particularly in humans, where the majority of genes generate splice variant transcripts. Characterizing the biological significance of alternative transcripts is methodologically difficult since it is impractical to assess thousands of splice variants as to whether they actually encode proteins, whether these proteins are functional, or whether transcripts have a function independent of protein synthesis. Consequently, to elucidate the functional significance of splice variants and to investigate mechanisms underlying the fidelity of mRNA splicing, we used an indirect approach based on analyzing the evolutionary conservation of splice variants among species. Using DNA polymerase ß as an indicator locus, we cloned and characterized the types and frequencies of transcripts generated in primary cell lines of five primate species. Overall, we found that in addition to the canonical DNA polymerase ß transcript, there were 25 alternative transcripts generated, most containing premature terminating codons. We used a statistical method borrowed from community ecology to show that there is significant diversity and little conservation in alternative splicing patterns among species, despite high sequence similarity in the underlying genomic (exonic) sequences. However, the frequency of alternative splicing at this locus correlates well with life history parameters such as the maximal longevity of each species, indicating that the alternative splicing of unproductive splice variants may have adaptive significance, even if the specific RNA transcripts themselves have no function. These results demonstrate the validity of the phylogenetic conservation approach in elucidating the biological significance of alternative splicing.

Parentage-based pedigree reconstruction reveals female matrilineal clusters and male-biased dispersal in nongregarious Asian great apes, the Bornean orang-utans (Pongo pygmaeus).

Philopatry and sex-biased dispersal have a strong influence on population genetic structure, so the study of species dispersal patterns and evolutionary mechanisms shaping them are of great interest. Particularly nongregarious mammalian species present an underexplored field of study: despite their lower levels of sociality compared to group-living species, interactions among individuals do occur, providing opportunities for cryptic kin selection. Among the least gregarious primates are orang-utans (genus: Pongo), in which preferential associations among females have nevertheless been observed, but for which the presence of kin structures was so far unresolved because of the equivocal results of previous genetic studies. To clarify relatedness and dispersal patterns in orang-utans, we examined the largest longitudinal set of individuals with combined genetic, spatial and behavioural data. We found that males had significantly higher mitochondrial DNA (mtDNA) variation and more unique haplotypes, thus underscoring their different maternal ancestries compared to females. Moreover, pedigree reconstruction based on 24 highly polymorphic microsatellite markers and mtDNA haplotypes demonstrated the presence of three matrilineal clusters of generally highly related females with substantially overlapping ranges. In orang-utans and possibly other nongregarious species, comparing average biparental relatedness (r) of males and females to infer sex-biased dispersal is extremely problematic. This is because the opportunistic sampling regime frequently employed in nongregarious species, combined with overlapping space use of distinct matrilineal clusters, leads to a strong downward bias when mtDNA lineage membership is ignored. Thus, in nongregarious species, correct inferences of dispersal can only be achieved by combining several genetic approaches with detailed spatial information.

Expression profiling in primates reveals a rapid evolution of human transcription factors.

Although it has been hypothesized for thirty years that many human adaptations are likely to be due to changes in gene regulation, almost nothing is known about the modes of natural selection acting on regulation in primates. Here we identify a set of genes for which expression is evolving under natural selection. We use a new multi-species complementary DNA array to compare steady-state messenger RNA levels in liver tissues within and between humans, chimpanzees, orangutans and rhesus macaques. Using estimates from a linear mixed model, we identify a set of genes for which expression levels have remained constant across the entire phylogeny (approximately 70 million years), and are therefore likely to be under stabilizing selection. Among the top candidates are five genes with expression levels that have previously been shown to be altered in liver carcinoma. We also find a number of genes with similar expression levels among non-human primates but significantly elevated or reduced expression in the human lineage, features that point to the action of directional selection. Among the gene set with a human-specific increase in expression, there is an excess of transcription factors; the same is not true for genes with increased expression in chimpanzee.