Search

From the recent genome assembly NHGRI_mPonAbe1-v2.0_NCBI (GCF_028885655.2) of orangutan chromosome 13, we computed the precise alpha satellite higher-order repeat (HOR) structure using the novel high-precision GRM2023 algorithm with Global Repeat Map (GRM) and Monomer Distance (MD) diagrams. This study rigorously identified alpha satellite HORs in the centromere of orangutan chromosome 13, discovering a novel 59mer HOR-the longest HOR unit identified in any primate to date. Additionally, it revealed the first intertwined sequence of three HORs, 18mer/27mer/45mer HORs, with a common aligned "backbone" across all HOR copies. The major 7mer HOR exhibits a Willard's-type canonical copy, although some segments of the array display significant irregularities. In contrast, the 14mer HOR forms a regular Willard's-type HOR array. Surprisingly, the GRM2023 high-precision analysis of chromosome 13 of human genome assembly T2T-CHM13v2.0 reveals the presence of only a 7mer HOR, despite both the orangutan and human genome assemblies being derived from whole genome shotgun sequences.

We identified five distinct full-length human mineralocorticoid receptor (MR) genes containing either 984 amino acids (MR-984) or 988 amino acids (MR-988), which can be distinguished by the presence or absence of Lys, Cys, Ser, and Trp (KCSW) in their DNA-binding domain (DBD) and mutations at codons 180 and 241 in their amino-terminal domain (NTD). Two human MR-KCSW genes contain either (Val-180, Val-241) or (Ile-180, Val-241) in their NTD, and three human MR-984 genes contain either (Ile-180, Ala-241), (Val-180, Val-241), or (Ile-180, Val-241). Human MR-KCSW with (Ile-180, Ala-241) has not been cloned. In contrast, chimpanzees contain four MRs: two MR-988s with KCSW in their DBD, or two MR-984s without KCSW in their DBD. Chimpanzee MRs only contain (Ile180, Val-241) in their NTD. A chimpanzee MR with either (Val-180, Val-241) or (Ile-180, Ala-241) in the NTD has not been cloned. Gorillas and orangutans each contain one MR-988 with KCSW in the DBD and one MR-984 without KCSW, and these MRs only contain (Ile-180, Val-241) in their NTD. A gorilla MR or orangutan MR with either (Val-180, Val-241) or (Ile-180, Ala-241) in the NTD has not been cloned. Together, these data suggest that human MRs with (Val-180, Val-241) or (Ile-180, Ala-241) in the NTD evolved after humans and chimpanzees diverged from their common ancestor. Considering the multiple functions in human development of the MR in kidney, brain, heart, skin, and lungs, as well as MR activity in interaction with the glucocorticoid receptor, we suggest that the evolution of human MRs that are absent in chimpanzees may have been important in the evolution of humans from chimpanzees. Investigation of the physiological responses to corticosteroids mediated by the MR in humans, chimpanzees, gorillas, and orangutans may provide insights into the evolution of humans and their closest relatives.

Human centromeres have been traditionally very difficult to sequence and assemble owing to their repetitive nature and large size 1 . As a result, patterns of human centromeric variation and models for their evolution and function remain incomplete, despite centromeres being among the most rapidly mutating regions 2,3 . Here, using long-read sequencing, we completely sequenced and assembled all centromeres from a second human genome and compared it to the finished reference genome 4,5 . We find that the two sets of centromeres show at least a 4.1-fold increase in single-nucleotide variation when compared with their unique flanks and vary up to 3-fold in size. Moreover, we find that 45.8% of centromeric sequence cannot be reliably aligned using standard methods owing to the emergence of new α-satellite higher-order repeats (HORs). DNA methylation and CENP-A chromatin immunoprecipitation experiments show that 26% of the centromeres differ in their kinetochore position by >500 kb. To understand evolutionary change, we selected six chromosomes and sequenced and assembled 31 orthologous centromeres from the common chimpanzee, orangutan and macaque genomes. Comparative analyses reveal a nearly complete turnover of α-satellite HORs, with characteristic idiosyncratic changes in α-satellite HORs for each species. Phylogenetic reconstruction of human haplotypes supports limited to no recombination between the short (p) and long (q) arms across centromeres and reveals that novel α-satellite HORs share a monophyletic origin, providing a strategy to estimate the rate of saltatory amplification and mutation of human centromeric DNA., (© 2024. The Author(s).)

Cysts and trophozoites of vestibuliferid ciliates and larvae of Strongyloides were found in fecal samples from captive orangutans Pongo pygmaeus and P. abelii from Czech and Slovak zoological gardens. As comparative material, ciliates from semi-captive mandrills Mandrillus sphinx from Gabon were included in the study. Phylogenetic analysis of the detected vestibuliferid ciliates using ITS1-5.8s-rRNA-ITS2 and partial 18S ribosomal deoxyribonucleic acid (rDNA) revealed that the ciliates from orangutans are conspecific with Balantioides coli lineage A, while the ciliates from mandrills clustered with Buxtonella-like ciliates from other primates. Morphological examination of the cysts and trophozoites using light microscopy did not reveal differences robust enough to identify the genera of the ciliates. Phylogenetic analysis of detected L1 larvae of Strongyloides using partial cox1 revealed Strongyloides stercoralis clustering within the cox1 lineage A infecting dogs, humans, and other primates. The sequences of 18S rDNA support these results. As both B. coli and S. stercoralis are zoonotic parasites and the conditions in captive and semi-captive settings may facilitate transmission to humans, prophylactic measures should reflect the findings., (© 2023 The Authors. American Journal of Primatology published by Wiley Periodicals LLC.)

The growth hormone (GH) locus has experienced a dramatic evolution in primates, becoming multigenic and diverse in anthropoids. Despite sequence information from a vast number of primate species, it has remained unclear how the multigene family was favored. We compared the structure and composition of apes' GH loci as a prerequisite to understanding their origin and possible evolutionary role. These thorough analyses of the GH loci of the chimpanzee, gorilla, and orangutan were done by resorting to previously sequenced bacterial artificial chromosomes (BACs) harboring them, as well as to their respective genome projects data available in GenBank. The GH loci of modern man, Neanderthal, gibbon, and wild boar were retrieved from GenBank. Coding regions, regulatory elements, and repetitive sequences were identified and compared among species. The GH loci of all the analyzed species are flanked by the genes CD79B (5') and ICAM-1 (3'). In man, Neanderthal, and chimpanzee, the loci were integrated by five almost indistinguishable genes; however, in the former two, they rendered three different hormones, and in the latter, four different proteins were derived. Gorilla exhibited six genes, gibbon seven, and orangutan four. The sequences of the proximal promoters, enhancers, P-elements, and a locus control region (LCR) were highly conserved. The locus evolution might have implicated duplications of the ancestral pituitary gene ( GH-N ) and subsequent diversification of the copies, leading to the placental single GH-V gene and the multiple CSH genes.

A syndrome of chronic upper and lower airway disease leading to increased morbidity and mortality occurs primarily in captive orangutans. Similarities in symptoms to the inherited human respiratory disease, cystic fibrosis, led us to hypothesize that orangutan respiratory disease is a result of variants in the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) gene. We identified the nonsense variant, c.484A>T (p.Lys162X), in heterozygosity in an unaffected orangutan. Analysis of the pedigree of this orangutan confirmed that both his sire and deceased fetus also harbored the c.484A>T allele. An expression minigene harboring c.484A>T produced no full-length CFTR protein in HEK293 cells. Finally, the c.484A>T CFTR messenger RNA abundance was severely reduced in primary nasal epithelial cells of the orangutan indicating that c.484A>T (p.Lys162X) is potentially lethal. Genetic screening of the captive orangutan population could be used to prevent transmission of this potentially lethal variant, and thus aid in the conservation of this critically endangered species., (© 2020 Wiley Periodicals, Inc.)

The divergence of chimpanzee and bonobo provides one of the few examples of recent hominid speciation 1,2 . Here we describe a fully annotated, high-quality bonobo genome assembly, which was constructed without guidance from reference genomes by applying a multiplatform genomics approach. We generate a bonobo genome assembly in which more than 98% of genes are completely annotated and 99% of the gaps are closed, including the resolution of about half of the segmental duplications and almost all of the full-length mobile elements. We compare the bonobo genome to those of other great apes 1,3-5 and identify more than 5,569 fixed structural variants that specifically distinguish the bonobo and chimpanzee lineages. We focus on genes that have been lost, changed in structure or expanded in the last few million years of bonobo evolution. We produce a high-resolution map of incomplete lineage sorting and estimate that around 5.1% of the human genome is genetically closer to chimpanzee or bonobo and that more than 36.5% of the genome shows incomplete lineage sorting if we consider a deeper phylogeny including gorilla and orangutan. We also show that 26% of the segments of incomplete lineage sorting between human and chimpanzee or human and bonobo are non-randomly distributed and that genes within these clustered segments show significant excess of amino acid replacement compared to the rest of the genome.

Background: Orang-utans comprise three critically endangered species endemic to the islands of Borneo and Sumatra. Though whole-genome sequencing has recently accelerated our understanding of their evolutionary history, the costs of implementing routine genome screening and diagnostics remain prohibitive. Capitalizing on a tri-fold locus discovery approach, combining data from published whole-genome sequences, novel whole-exome sequencing, and microarray-derived genotype data, we aimed to develop a highly informative gene-focused panel of targets that can be used to address a broad range of research questions., Results: We identified and present genomic co-ordinates for 175,186 SNPs and 2315 Y-chromosomal targets, plus 185 genes either known or presumed to be pathogenic in cardiovascular (N = 109) or respiratory (N = 43) diseases in humans - the primary and secondary causes of captive orang-utan mortality - or a majority of other human diseases (N = 33). As proof of concept, we designed and synthesized 'SeqCap' hybrid capture probes for these targets, demonstrating cost-effective target enrichment and reduced-representation sequencing., Conclusions: Our targets are of broad utility in studies of orang-utan ancestry, admixture and disease susceptibility and aetiology, and thus are of value in addressing questions key to the survival of these species. To facilitate comparative analyses, these targets could now be standardized for future orang-utan population genomic studies. The targets are broadly compatible with commercial target enrichment platforms and can be utilized as published here to synthesize applicable probes.

The mammalian male-specific Y chromosome plays a critical role in sex determination and male fertility. However, because of its repetitive and haploid nature, it is frequently absent from genome assemblies and remains enigmatic. The Y chromosomes of great apes represent a particular puzzle: their gene content is more similar between human and gorilla than between human and chimpanzee, even though human and chimpanzee share a more recent common ancestor. To solve this puzzle, here we constructed a dataset including Ys from all extant great ape genera. We generated assemblies of bonobo and orangutan Ys from short and long sequencing reads and aligned them with the publicly available human, chimpanzee, and gorilla Y assemblies. Analyzing this dataset, we found that the genus Pan , which includes chimpanzee and bonobo, experienced accelerated substitution rates. Pan also exhibited elevated gene death rates. These observations are consistent with high levels of sperm competition in Pan Furthermore, we inferred that the great ape common ancestor already possessed multicopy sequences homologous to most human and chimpanzee palindromes. Nonetheless, each species also acquired distinct ampliconic sequences. We also detected increased chromatin contacts between and within palindromes (from Hi-C data), likely facilitating gene conversion and structural rearrangements. Our results highlight the dynamic mode of Y chromosome evolution and open avenues for studies of male-specific dispersal in endangered great ape species., Competing Interests: The authors declare no competing interest., (Copyright © 2020 the Author(s). Published by PNAS.)

Exploration is essential for skill acquisition and strongly facilitates cognitive performance. In humans, it is widely known that exploration and later cognitive performance are highly dependent on early social inputs. Here, we aim to shed light on the evolutionary roots of this process by studying the effects of variation in opportunities for social learning on the exploratory tendency of immature orangutans ( Pongo spp.) in nature. We based our analyses on mixed cross-sectional, longitudinal data of exploration by immatures and their mothers. Current exploration rates were correlated with levels of past experienced sociability, but not with current food abundance or with maternal condition, and only partly with genetic similarity. We conclude that the dependence of cognitive development on socially triggered exploration, which underlies the construction of cognitive skills such as intelligence, existed before the emergence of the human lineage., (Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).)

Changes in transcriptional regulation are thought to be a major contributor to the evolution of phenotypic traits, but the contribution of changes in chromatin accessibility to the evolution of gene expression remains almost entirely unknown. To address this important gap in knowledge, we developed a new method to identify DNase I Hypersensitive (DHS) sites with differential chromatin accessibility between species using a joint modeling approach. Our method overcomes several limitations inherent to conventional threshold-based pairwise comparisons that become increasingly apparent as the number of species analyzed rises. Our approach employs a single quantitative test which is more sensitive than existing pairwise methods. To illustrate, we applied our joint approach to DHS sites in fibroblast cells from five primates (human, chimpanzee, gorilla, orangutan, and rhesus macaque). We identified 89,744 DHS sites, of which 41% are identified as differential between species using the joint model compared with 33% using the conventional pairwise approach. The joint model provides a principled approach to distinguishing single from multiple chromatin accessibility changes among species. We found that nondifferential DHS sites are enriched for nucleotide conservation. Differential DHS sites with decreased chromatin accessibility relative to rhesus macaque occur more commonly near transcription start sites (TSS), while those with increased chromatin accessibility occur more commonly distal to TSS. Further, differential DHS sites near TSS are less cell type-specific than more distal regulatory elements. Taken together, these results point to distinct classes of DHS sites, each with distinct characteristics of selection, genomic location, and cell type specificity., (© The Author(s) 2019. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Background: Integrating demography and adaptive evolution is pivotal to understanding the evolutionary history and conservation of great apes. However, little is known about the adaptive evolution of our closest relatives, in particular if and to what extent adaptions to environmental differences have occurred. Here, we used whole-genome sequencing data from critically endangered orangutans from North Sumatra (Pongo abelii) and Borneo (P. pygmaeus) to investigate adaptive responses of each species to environmental differences during the Pleistocene., Results: Taking into account the markedly disparate demographic histories of each species after their split ~ 1 Ma ago, we show that persistent environmental differences on each island had a strong impact on the adaptive evolution of the genus Pongo. Across a range of tests for positive selection, we find a consistent pattern of between-island and species differences. In the more productive Sumatran environment, the most notable signals of positive selection involve genes linked to brain and neuronal development, learning, and glucose metabolism. On Borneo, however, positive selection comprised genes involved in lipid metabolism, as well as cardiac and muscle activities., Conclusions: We find strikingly different sets of genes appearing to have evolved under strong positive selection in each species. In Sumatran orangutans, selection patterns were congruent with well-documented cognitive and behavioral differences between the species, such as a larger and more complex cultural repertoire and higher degrees of sociality. However, in Bornean orangutans, selective responses to fluctuating environmental conditions appear to have produced physiological adaptations to generally lower and temporally more unpredictable food supplies.

We introduce ABLE (Approximate Blockwise Likelihood Estimation), a novel simulation-based composite likelihood method that uses the blockwise site frequency spectrum to jointly infer past demography and recombination. ABLE is explicitly designed for a wide variety of data from unphased diploid genomes to genome-wide multi-locus data (for example, RADSeq) and can also accommodate arbitrarily large samples. We use simulations to demonstrate the accuracy of this method to infer complex histories of divergence and gene flow and reanalyze whole genome data from two species of orangutan. ABLE is available for download at https://github.com/champost/ABLE .

Gene conversion is the copying of a genetic sequence from a "donor" region to an "acceptor." In nonallelic gene conversion (NAGC), the donor and the acceptor are at distinct genetic loci. Despite the role NAGC plays in various genetic diseases and the concerted evolution of gene families, the parameters that govern NAGC are not well characterized. Here, we survey duplicate gene families and identify converted tracts in 46% of them. These conversions reflect a large GC bias of NAGC. We develop a sequence evolution model that leverages substantially more information in duplicate sequences than used by previous methods and use it to estimate the parameters that govern NAGC in humans: a mean converted tract length of 250 bp and a probability of [Formula: see text] per generation for a nucleotide to be converted (an order of magnitude higher than the point mutation rate). Despite this high baseline rate, we show that NAGC slows down as duplicate sequences diverge-until an eventual "escape" of the sequences from its influence. As a result, NAGC has a small average effect on the sequence divergence of duplicates. This work improves our understanding of the NAGC mechanism and the role that it plays in the evolution of gene duplicates., Competing Interests: The authors declare no conflict of interest.

Six extant species of non-human great apes are currently recognized: Sumatran and Bornean orangutans, eastern and western gorillas, and chimpanzees and bonobos [1]. However, large gaps remain in our knowledge of fine-scale variation in hominoid morphology, behavior, and genetics, and aspects of great ape taxonomy remain in flux. This is particularly true for orangutans (genus: Pongo), the only Asian great apes and phylogenetically our most distant relatives among extant hominids [1]. Designation of Bornean and Sumatran orangutans, P. pygmaeus (Linnaeus 1760) and P. abelii (Lesson 1827), as distinct species occurred in 2001 [1, 2]. Here, we show that an isolated population from Batang Toru, at the southernmost range limit of extant Sumatran orangutans south of Lake Toba, is distinct from other northern Sumatran and Bornean populations. By comparing cranio-mandibular and dental characters of an orangutan killed in a human-animal conflict to those of 33 adult male orangutans of a similar developmental stage, we found consistent differences between the Batang Toru individual and other extant Ponginae. Our analyses of 37 orangutan genomes provided a second line of evidence. Model-based approaches revealed that the deepest split in the evolutionary history of extant orangutans occurred ∼3.38 mya between the Batang Toru population and those to the north of Lake Toba, whereas both currently recognized species separated much later, about 674 kya. Our combined analyses support a new classification of orangutans into three extant species. The new species, Pongo tapanuliensis, encompasses the Batang Toru population, of which fewer than 800 individuals survive. VIDEO ABSTRACT., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

The characterization of mutational processes that generate sequence diversity in the human genome is of paramount importance both to medical genetics and to evolutionary studies. To understand how the age and sex of transmitting parents affect de novo mutations, here we sequence 1,548 Icelanders, their parents, and, for a subset of 225, at least one child, to 35× genome-wide coverage. We find 108,778 de novo mutations, both single nucleotide polymorphisms and indels, and determine the parent of origin of 42,961. The number of de novo mutations from mothers increases by 0.37 per year of age (95% CI 0.32-0.43), a quarter of the 1.51 per year from fathers (95% CI 1.45-1.57). The number of clustered mutations increases faster with the mother's age than with the father's, and the genomic span of maternal de novo mutation clusters is greater than that of paternal ones. The types of de novo mutation from mothers change substantially with age, with a 0.26% (95% CI 0.19-0.33%) decrease in cytosine-phosphate-guanine to thymine-phosphate-guanine (CpG>TpG) de novo mutations and a 0.33% (95% CI 0.28-0.38%) increase in C>G de novo mutations per year, respectively. Remarkably, these age-related changes are not distributed uniformly across the genome. A striking example is a 20 megabase region on chromosome 8p, with a maternal C>G mutation rate that is up to 50-fold greater than the rest of the genome. The age-related accumulation of maternal non-crossover gene conversions also mostly occurs within these regions. Increased sequence diversity and linkage disequilibrium of C>G variants within regions affected by excess maternal mutations indicate that the underlying mutational process has persisted in humans for thousands of years. Moreover, the regional excess of C>G variation in humans is largely shared by chimpanzees, less by gorillas, and is almost absent from orangutans. This demonstrates that sequence diversity in humans results from evolving interactions between age, sex, mutation type, and genomic location.

Simian foamy viruses (SFVs) infect most nonhuman primate species and appears to co-evolve with its hosts. This co-evolutionary signal is particularly strong among great apes, including orangutans (genus Pongo). Previous studies have identified three distinct orangutan SFV clades. The first of these three clades is composed of SFV from P. abelii from Sumatra, the second consists of SFV from P. pygmaeus from Borneo, while the third clade is mixed, comprising an SFV strain found in both species of orangutan. The existence of the mixed clade has been attributed to an expansion of P. pygmaeus into Sumatra following the Mount Toba super-volcanic eruption about 73,000years ago. Divergence dating, however, has yet to be performed to establish a temporal association with the Toba eruption. Here, we use a Bayesian framework and a relaxed molecular clock model with fossil calibrations to test the Toba hypothesis and to gain a more complete understanding of the evolutionary history of orangutan SFV. As with previous studies, our results show a similar three-clade orangutan SFV phylogeny, along with strong statistical support for SFV-host co-evolution in orangutans. Using Bayesian inference, we date the origin of orangutan SFV to >4.7 million years ago (mya), while the mixed species clade dates to approximately 1.7mya, >1.6 million years older than the Toba super-eruption. These results, combined with fossil and paleogeographic evidence, suggest that the origin of SFV in Sumatran and Bornean orangutans, including the mixed species clade, likely occurred on the mainland of Indo-China during the Late Pliocene and Calabrian stage of the Pleistocene, respectively., (Copyright © 2017 Elsevier B.V. All rights reserved.)

The immune and reproductive functions of human NK cells are regulated by interactions of the C1 and C2 epitopes of HLA-C with C1-specific and C2-specific lineage III killer cell Ig-like receptors (KIR). This rapidly evolving and diverse system of ligands and receptors is restricted to humans and great apes. In this context, the orangutan has particular relevance because it represents an evolutionary intermediate, one having the C1 epitope and corresponding KIR but lacking the C2 epitope. Through a combination of direct sequencing, KIR genotyping, and data mining from the Great Ape Genome Project, we characterized the KIR alleles and haplotypes for panels of 10 Bornean orangutans and 19 Sumatran orangutans. The orangutan KIR haplotypes have between 5 and 10 KIR genes. The seven orangutan lineage III KIR genes all locate to the centromeric region of the KIR locus, whereas their human counterparts also populate the telomeric region. One lineage III KIR gene is Bornean specific, one is Sumatran specific, and five are shared. Of 12 KIR gene-content haplotypes, 5 are Bornean specific, 5 are Sumatran specific, and 2 are shared. The haplotypes have different combinations of genes encoding activating and inhibitory C1 receptors that can be of higher or lower affinity. All haplotypes encode an inhibitory C1 receptor, but only some haplotypes encode an activating C1 receptor. Of 130 KIR alleles, 55 are Bornean specific, 65 are Sumatran specific, and 10 are shared., (Copyright © 2017 by The American Association of Immunologists, Inc.)

The site frequency spectrum (SFS) has long been used to study demographic history and natural selection. Here, we extend this summary by examining the SFS conditional on the alleles found at the same site in other species. We refer to this extension as the "phylogenetically-conditioned SFS" or cSFS. Using recent large-sample data from the Exome Aggregation Consortium (ExAC), combined with primate genome sequences, we find that human variants that occurred independently in closely related primate lineages are at higher frequencies in humans than variants with parallel substitutions in more distant primates. We show that this effect is largely due to sites with elevated mutation rates causing significant departures from the widely-used infinite sites mutation model. Our analysis also suggests substantial variation in mutation rates even among mutations involving the same nucleotide changes. In summary, we show that variable mutation rates are key determinants of the SFS in humans., Competing Interests: The authors have declared that no competing interests exist.

The great apes are the closest living relatives of humans. Chimpanzees and bonobos group together with humans, while gorillas and orangutans are more divergent from humans. Here, we review insights into their evolution pertaining to the topology of species and subspecies and the reconstruction of their demography based on genome-wide variation. These advances have only become possible recently through next-generation sequencing technologies. Given the close relationship to humans, they provide an important evolutionary context for human genetics., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

In primates, the unigenic growth hormone (GH) locus of prosimians expressed primarily in the anterior pituitary, evolved by gene duplications, independently in New World Monkeys (NWM) and Old World Monkeys (OWMs)/apes, to give complex clusters of genes expressed in the pituitary and placenta. In human and chimpanzee, the GH locus comprises five genes, GH-N being expressed as pituitary GH, whereas GH-V (placental GH) and CSHs (chorionic somatomammotropins) are expressed (in human and probably chimpanzee) in the placenta; the CSHs comprise CSH-A, CSH-B and the aberrant CSH-L (possibly a pseudogene) in human, and CSH-A1, CSH-A2 and CSH-B in chimpanzee. Here, the GH locus in two additional great apes, gorilla (Gorilla gorilla gorilla) and orangutan (Pongo abelii), is shown to contain six and four GH-like genes, respectively. The gorilla locus possesses six potentially expressed genes, gGH-N, gGH-V and four gCSHs, whereas the orangutan locus has just three functional genes, oGH-N, oGH-V and oCSH-B, plus a pseudogene, oCSH-L. Analysis of regulatory sequences, including promoter, enhancer and P-elements, shows significant variation; in particular the proximal Pit-1 element of GH-V genes differs markedly from that of other genes in the cluster. Phylogenetic analysis shows that the initial gene duplication led to distinct GH-like and CSH-like genes and that a second duplication provided separate GH-N and GH-V. However, evolution of the CSH-like genes remains unclear. Rapid adaptive evolution gave rise to the distinct CSHs, after the first duplication, and to GH-V after the second duplication. Analysis of transcriptomic databases derived from gorilla tissues establishes that the gGH-N, gGH-V and several gCSH genes are expressed, but the significance of the many CSH genes in gorilla remains unclear.

Genetic differences that specify unique aspects of human evolution have typically been identified by comparative analyses between the genomes of humans and closely related primates, including more recently the genomes of archaic hominins. Not all regions of the genome, however, are equally amenable to such study. Recurrent copy number variation (CNV) at chromosome 16p11.2 accounts for approximately 1% of cases of autism and is mediated by a complex set of segmental duplications, many of which arose recently during human evolution. Here we reconstruct the evolutionary history of the locus and identify bolA family member 2 (BOLA2) as a gene duplicated exclusively in Homo sapiens. We estimate that a 95-kilobase-pair segment containing BOLA2 duplicated across the critical region approximately 282 thousand years ago (ka), one of the latest among a series of genomic changes that dramatically restructured the locus during hominid evolution. All humans examined carried one or more copies of the duplication, which nearly fixed early in the human lineage--a pattern unlikely to have arisen so rapidly in the absence of selection (P < 0.0097). We show that the duplication of BOLA2 led to a novel, human-specific in-frame fusion transcript and that BOLA2 copy number correlates with both RNA expression (r = 0.36) and protein level (r = 0.65), with the greatest expression difference between human and chimpanzee in experimentally derived stem cells. Analyses of 152 patients carrying a chromosome 16p11. rearrangement show that more than 96% of breakpoints occur within the H. sapiens-specific duplication. In summary, the duplicative transposition of BOLA2 at the root of the H. sapiens lineage about 282 ka simultaneously increased copy number of a gene associated with iron homeostasis and predisposed our species to recurrent rearrangements associated with disease.

Background: The rise of the primate lineage is accompanied by an outstanding emergence of microRNAs, small non-coding RNAs with a prominent role in gene regulation. In spite of their biological importance little is known about the way in which natural selection has influenced microRNAs in the human lineage. To study the recent evolutionary history of human microRNAs and to analyze the signatures of natural selection in genomic regions harbouring microRNAs we have investigated the nucleotide substitution rates of 1,872 human microRNAs in the human and chimpanzee lineages., Results: We produced a depurated set of microRNA alignments of human, chimpanzee and orang-utan orthologs combining BLAT and liftOver and selected 1,214 microRNA precursors presenting optimal secondary structures. We classified microRNAs in categories depending on their genomic organization, duplication status and conservation along evolution. We compared substitution rates of the aligned microRNAs between human and chimpanzee using Tajima's Relative Rate Test taking orang-utan as out-group and found several microRNAs with particularly high substitution rates in either the human or chimpanzee branches. We fitted different models of natural selection on these orthologous microRNA alignments and compared them using a likelihood ratio test that uses ancestral repeats and microRNA flanking regions as neutral sequences. We found that although a large fraction of human microRNAs is highly conserved among the three species studied, significant differences in rates of molecular evolution exist among microRNA categories. Particularly, primate-specific microRNAs, which are enriched in isolated and single copy microRNAs, more than doubled substitution rates of those belonging to older, non primate-specific microRNA families., Conclusions: Our results corroborate the remarkable conservation of microRNAs, a proxy of their functional relevance, and indicate that a subset of human microRNAs undergo nucleotide substitutions at higher rates, which may be suggestive of the action of positive selection.

Vocal fold control was critical to the evolution of spoken language, much as it today allows us to learn vowel systems. It has, however, never been demonstrated directly in a non-human primate, leading to the suggestion that it evolved in the human lineage after divergence from great apes. Here, we provide the first evidence for real-time, dynamic and interactive vocal fold control in a great ape during an imitation "do-as-I-do" game with a human demonstrator. Notably, the orang-utan subject skilfully produced "wookies" - an idiosyncratic vocalization exhibiting a unique spectral profile among the orang-utan vocal repertoire. The subject instantaneously matched human-produced wookies as they were randomly modulated in pitch, adjusting his voice frequency up or down when the human demonstrator did so, readily generating distinct low vs. high frequency sub-variants. These sub-variants were significantly different from spontaneous ones (not produced in matching trials). Results indicate a latent capacity for vocal fold exercise in a great ape (i) in real-time, (ii) up and down the frequency spectrum, (iii) across a register range beyond the species-repertoire and, (iv) in a co-operative turn-taking social setup. Such ancestral capacity likely provided the neuro-behavioural basis of the more fine-tuned vocal fold control that is a human hallmark.

De novo protein-coding gene origination is increasingly recognized as an important evolutionary mechanism. However, there remains a large amount of uncertainty regarding the frequency of these events and the mechanisms and speed of gene establishment. Here, we describe a rigorous search for cases of de novo gene origination in the great apes. We analyzed annotated proteomes as well as full genomic DNA and transcriptional and translational evidence. It is notable that results vary between database updates due to the fluctuating annotation of these genes. Nonetheless we identified 35 de novo genes: 16 human-specific; 5 human and chimpanzee specific; and 14 that originated prior to the divergence of human, chimpanzee, and gorilla and are found in all three genomes. The taxonomically restricted distribution of these genes cannot be explained by loss in other lineages. Each gene is supported by an open reading frame-creating mutation that occurred within the primate lineage, and which is not polymorphic in any species. Similarly to previous studies we find that the de novo genes identified are short and frequently located near pre-existing genes. Also, they may be associated with Alu elements and prior transcription and RNA-splicing at the locus. Additionally, we report the first case of apparent independent lineage sorting of a de novo gene. The gene is present in human and gorilla, whereas chimpanzee has the ancestral noncoding sequence. This indicates a long period of polymorphism prior to fixation and thus supports a model where de novo genes may, at least initially, have a neutral effect on fitness., (© The Author(s) 2016. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

microRNAs are crucial post-transcriptional regulators of gene expression involved in a wide range of biological processes. Although microRNAs are highly conserved among species, the functional implications of existing lineage-specific changes and their role in determining differences between humans and other great apes have not been specifically addressed. We analyzed the recent evolutionary history of 1,595 human microRNAs by looking at their intra- and inter-species variation in great apes using high-coverage sequenced genomes of 82 individuals including gorillas, orangutans, bonobos, chimpanzees and humans. We explored the strength of purifying selection among microRNA regions and found that the seed and mature regions are under similar and stronger constraint than the precursor region. We further constructed a comprehensive catalogue of microRNA species-specific nucleotide substitutions among great apes and, for the first time, investigated the biological relevance that human-specific changes in microRNAs may have had in great ape evolution. Expression and functional analyses of four microRNAs (miR-299-3p, miR-503-3p, miR-508-3p and miR-541-3p) revealed that lineage-specific nucleotide substitutions and changes in the length of these microRNAs alter their expression as well as the repertoires of target genes and regulatory networks. We suggest that the studied molecular changes could have modified crucial microRNA functions shaping phenotypes that, ultimately, became human-specific. Our work provides a frame to study the impact that regulatory changes may have in the recent evolution of our species.

Background: Orangutans are an endangered species whose natural habitats are restricted to the Southeast Asian islands of Borneo and Sumatra. Along with the African great apes, orangutans are among the closest living relatives to humans. For potential species conservation and functional genomics studies, we derived induced pluripotent stem cells (iPSCs) from cryopreserved somatic cells obtained from captive orangutans., Results: Primary skin fibroblasts from two Sumatran orangutans were transduced with retroviral vectors expressing the human OCT4, SOX2, KLF4, and c-MYC factors. Candidate orangutan iPSCs were characterized by global gene expression and DNA copy number analysis. All were consistent with pluripotency and provided no evidence of large genomic insertions or deletions. In addition, orangutan iPSCs were capable of producing cells derived from all three germ layers in vitro through embryoid body differentiation assays and in vivo through teratoma formation in immune-compromised mice., Conclusions: We demonstrate that orangutan skin fibroblasts are capable of being reprogrammed into iPSCs with hallmark molecular signatures and differentiation potential. We suggest that reprogramming orangutan somatic cells in genome resource banks could provide new opportunities for advancing assisted reproductive technologies relevant for species conservation efforts. Furthermore, orangutan iPSCs could have applications for investigating the phenotypic relevance of genomic changes that occurred in the human, African great ape, and/or orangutan lineages. This provides opportunities for orangutan cell culture models that would otherwise be impossible to develop from living donors due to the invasive nature of the procedures required for obtaining primary cells.

Androgen receptor (AR) mediated signalling is necessary for normal development of the prostate gland and also drives prostate cancer (PCa) cell growth and survival, with many studies showing a correlation between increased receptor levels and therapy resistance with progression to fatal castrate recurrent PCa (CRPC). Although it has been held for some time that the transcription factor Sp1 is the main stimulator of AR gene transcription, comprehensive knowledge of the regulation of the AR gene remains incomplete. Here we describe and characterise in detail two novel active regulatory elements in the 5'UTR of the human AR gene. Both of these elements contain overlapping binding sites for the positive transcription factor Sp1 and the repressor protein pur-α. Aberrant cell signalling is characteristic of PCa and the transcriptional activity of the AR promoter in PCa cells is dependent upon the relative amounts of the two transcription factors. Together with our corroboration of the dominant role of Sp1, the findings support the rationale of targeting this transcription factor to inhibit tumour progression. This should be of particular therapeutic relevance in CRPC where the levels of the repressor pur-α are reduced.

Small noncoding regulatory RNA exist in wide spectrum of organisms ranging from prokaryote bacteria to humans. In human, a systematic search for noncoding RNA is mainly limited to the nuclear and cytosolic compartments. To investigate whether endogenous small regulatory RNA are present in cell organelles, human mitochondrial genome was also explored for prediction of precursor microRNA (pre-miRNA) and mature miRNA (miRNA) sequences. Six novel miRNA were predicted from the organelle genome by bioinformatics analysis. The structures are conserved in other five mammals including chimp, orangutan, mouse, rat, and rhesus genome. Experimentally, six human miRNA are well accumulated or deposited in human mitochondria. Three of them are expressed less prominently in Northern analysis. To ascertain their presence in human skeletal muscles, total RNA was extracted from enriched mitochondria by an immunomagnetic method. The expression of six novel pre-miRNA and miRNA was confirmed by Northern blot analysis; however, low level of remaining miRNA was found by sensitive Northern analysis. Their presence is further confirmed by real time RT-PCR. The six miRNA find their multiple targets throughout the human genome in three different types of software. The luciferase assay was used to confirm that MT-RNR2 gene was the potential target of hsa-miR-mit3 and hsa-miR-mit4.

Understanding the demographic history of populations and species is a central issue in evolutionary biology and molecular ecology. In this work, we develop a maximum-likelihood method for the inference of past changes in population size from microsatellite allelic data. Our method is based on importance sampling of gene genealogies, extended for new mutation models, notably the generalized stepwise mutation model (GSM). Using simulations, we test its performance to detect and characterize past reductions in population size. First, we test the estimation precision and confidence intervals coverage properties under ideal conditions, then we compare the accuracy of the estimation with another available method (MSVAR) and we finally test its robustness to misspecification of the mutational model and population structure. We show that our method is very competitive compared with alternative ones. Moreover, our implementation of a GSM allows more accurate analysis of microsatellite data, as we show that the violations of a single step mutation assumption induce very high bias toward false contraction detection rates. However, our simulation tests also showed some limits, which most importantly are large computation times for strong disequilibrium scenarios and a strong influence of some form of unaccounted population structure. This inference method is available in the latest implementation of the MIGRAINE software package., (© The Author 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

We claim that the apparently enhanced CpG transversions in the form CpG to CpC/GpG or to ApG/CpT are caused by the hypermutable CpG to CpA/TpG transition. The nucleotide replacement counts obtained from the human/chimpanzee/gorilla/orangutan sequence alignments representing the replacements due to the evolutionary species divergence and the results of 1000 genomes project that provide us with the differences due to the intraspecies diversification were analyzed to estimate the ratio of CpG versus non-CpG transversion probabilities. The trinucleotide replacement counts were extracted from the regions that are free of functional constraints. The CpG transversion probabilities based upon the genomic comparisons were found to exceed more than twice the non-CpG transversions. The diversity data emerging from 14 population groups were partitioned in five classes as a function of the parameter quantifying the spread of the polymorphic allele among the group of individuals. The results based upon the human polymorphism exhibit a trend where CpG over non-CpG transversion probability ratio is less and less exceeding unity as the values of the derived allele frequency (DAF) of snps are diminishing. A computer simulation of a simplified model indicates that the phenomenon of the apparent enhancement of CpG transversions can have its source in the interference of the entropic effects with the maximum likelihood methodologies.

Background: Hidden Markov models are widely used for genome analysis as they combine ease of modelling with efficient analysis algorithms. Calculating the likelihood of a model using the forward algorithm has worst case time complexity linear in the length of the sequence and quadratic in the number of states in the model. For genome analysis, however, the length runs to millions or billions of observations, and when maximising the likelihood hundreds of evaluations are often needed. A time efficient forward algorithm is therefore a key ingredient in an efficient hidden Markov model library., Results: We have built a software library for efficiently computing the likelihood of a hidden Markov model. The library exploits commonly occurring substrings in the input to reuse computations in the forward algorithm. In a pre-processing step our library identifies common substrings and builds a structure over the computations in the forward algorithm which can be reused. This analysis can be saved between uses of the library and is independent of concrete hidden Markov models so one preprocessing can be used to run a number of different models.Using this library, we achieve up to 78 times shorter wall-clock time for realistic whole-genome analyses with a real and reasonably complex hidden Markov model. In one particular case the analysis was performed in less than 8 minutes compared to 9.6 hours for the previously fastest library., Conclusions: We have implemented the preprocessing procedure and forward algorithm as a C++ library, zipHMM, with Python bindings for use in scripts. The library is available at http://birc.au.dk/software/ziphmm/.

To gain insights into evolutionary forces that have shaped the history of Bornean and Sumatran populations of orang-utans, we compare patterns of variation across more than 11 million single nucleotide polymorphisms found by previous mitochondrial and autosomal genome sequencing of 10 wild-caught orang-utans. Our analysis of the mitochondrial data yields a far more ancient split time between the two populations (~3.4 million years ago) than estimates based on autosomal data (0.4 million years ago), suggesting a complex speciation process with moderate levels of primarily male migration. We find that the distribution of selection coefficients consistent with the observed frequency spectrum of autosomal non-synonymous polymorphisms in orang-utans is similar to the distribution in humans. Our analysis indicates that 35% of genes have evolved under detectable negative selection. Overall, our findings suggest that purifying natural selection, genetic drift, and a complex demographic history are the dominant drivers of genome evolution for the two orang-utan populations.

The genomes of related species contain valuable information on the history of the considered taxa. Great apes in particular exhibit variation of evolutionary patterns along their genomes. However, the great ape data also bring new challenges, such as the presence of incomplete lineage sorting and ancestral shared polymorphisms. Previous methods for genome-scale analysis are restricted to very few individuals or cannot disentangle the contribution of mutation rates and fixation biases. This represents a limitation both for the understanding of these forces as well as for the detection of regions affected by selection. Here, we present a new model designed to estimate mutation rates and fixation biases from genetic variation within and between species. We relax the assumption of instantaneous substitutions, modeling substitutions as mutational events followed by a gradual fixation. Hence, we straightforwardly account for shared ancestral polymorphisms and incomplete lineage sorting. We analyze genome-wide synonymous site alignments of human, chimpanzee, and two orangutan species. From each taxon, we include data from several individuals. We estimate mutation rates and GC-biased gene conversion intensity. We find that both mutation rates and biased gene conversion vary with GC content. We also find lineage-specific differences, with weaker fixation biases in orangutan species, suggesting a reduced historical effective population size. Finally, our results are consistent with directional selection acting on coding sequences in relation to exonic splicing enhancers.

It is a challenge to classify protein-coding or non-coding transcripts, especially those re-constructed from high-throughput sequencing data of poorly annotated species. This study developed and evaluated a powerful signature tool, Coding-Non-Coding Index (CNCI), by profiling adjoining nucleotide triplets to effectively distinguish protein-coding and non-coding sequences independent of known annotations. CNCI is effective for classifying incomplete transcripts and sense-antisense pairs. The implementation of CNCI offered highly accurate classification of transcripts assembled from whole-transcriptome sequencing data in a cross-species manner, that demonstrated gene evolutionary divergence between vertebrates, and invertebrates, or between plants, and provided a long non-coding RNA catalog of orangutan. CNCI software is available at http://www.bioinfo.org/software/cnci.

A long-standing question in evolutionary biology concerns the effect of recombination in shaping the genomic architecture of organisms and, in particular, how this impacts the speciation process. Despite efforts employed in the last decade, the role of chromosomal reorganizations in the human-chimpanzee speciation process remains unresolved. Through whole-genome comparisons, we have analyzed the genome-wide impact of genomic shuffling in the distribution of human recombination rates during the human-chimpanzee speciation process. We have constructed a highly refined map of the reorganizations and evolutionary breakpoint regions in the human and chimpanzee genomes based on orthologous genes and genome sequence alignments. The analysis of the most recent human and chimpanzee recombination maps inferred from genome-wide single-nucleotide polymorphism data revealed that the standardized recombination rate was significantly lower in rearranged than in collinear chromosomes. In fact, rearranged chromosomes presented significantly lower recombination rates than chromosomes that have been maintained since the ancestor of great apes, and this was related with the lineage in which they become fixed. Importantly, inverted regions had lower recombination rates than collinear and noninverted regions, independently of the effect of centromeres. Our observations have implications for the chromosomal speciation theory, providing new evidences for the contribution of inversions in suppressing recombination in mammals.

Introduction: Recent evidence from chimpanzees and gorillas has raised doubts that preeclampsia is a uniquely human disease. The deep extravillous trophoblast (EVT) invasion and spiral artery remodeling that characterizes our placenta (and is abnormal in preeclampsia) is shared within great apes, setting Homininae apart from Hylobatidae and Old World Monkeys, which show much shallower trophoblast invasion and limited spiral artery remodeling. We hypothesize that the evolution of a more invasive placenta in the lineage ancestral to the great apes involved positive selection on genes crucial to EVT invasion and spiral artery remodeling. Furthermore, identification of placentally-expressed genes under selection in this lineage may identify novel genes involved in placental development., Methods: We tested for positive selection in approximately 18,000 genes using the ratio of non-synonymous to synonymous amino acid substitution for protein-coding DNA. DAVID Bioinformatics Resources identified biological processes enriched in positively selected genes, including processes related to EVT invasion and spiral artery remodeling., Results: Analyses revealed 295 and 264 genes under significant positive selection on the branches ancestral to Hominidae (Human, Chimp, Gorilla, Orangutan) and Homininae (Human, Chimp, Gorilla), respectively. Gene ontology analysis of these gene sets demonstrated significant enrichments for several functional gene clusters relevant to preeclampsia risk, and sets of placentally-expressed genes that have been linked with preeclampsia and/or trophoblast invasion in other studies., Conclusion: Our study represents a novel approach to the identification of candidate genes and amino acid residues involved in placental pathologies by implicating them in the evolution of highly-invasive placenta., (Copyright © 2012 Elsevier Ltd. All rights reserved.)

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., Competing Interests: The authors have declared that no competing interests exist.

Unbiased readings of fossils are well known to contradict some of the popular molecular groupings among primates, particularly with regard to great apes and tarsiers. The molecular methodologies today are however flawed as they are based on a mistaken theoretical interpretation of the genetic equidistance phenomenon that originally started the field. An improved molecular method the 'slow clock' was here developed based on the Maximum Genetic Diversity hypothesis, a more complete account of the unified changes in genotypes and phenotypes. The method makes use of only slow evolving sequences and requires no uncertain assumptions or mathematical corrections and hence is able to give definitive results. The findings indicate that humans are genetically more distant to orangutans than African apes are and separated from the pongid clade ∼17.6 million years ago. Also, tarsiers are genetically closer to lorises than simian primates are. Finally, the fossil times for the radiation of mammals at the K/T boundary and for the Eutheria-Metatheria split in the Early Cretaceous were independently confirmed from molecular dating calibrated using the fossil split times of gorilla-orangutan, mouse-rat, and opossum-kangaroo. Therefore, the re-established primate phylogeny indicates a remarkable unity between molecules and fossils.

Mating systems are thought to be an important determinant of dispersal strategies in most animals, including the great apes. As the most basal taxon of all great apes, orang-utans can provide information about the evolution of mating systems and their consequences for population structure in this Family. To assess the sex-specific population structure in orang-utans, we used a combination of paternally transmitted Y-chromosomal genetic markers and maternally transmitted mitochondrial DNA sequences. Markers transmitted through the more philopatric sex are expected to show stronger differentiation among populations than the ones transmitted through the dispersing sex. We studied these patterns using 70 genetic samples from wild orang-utans from seven Bornean and two Sumatran populations. We found pronounced population structure in haplotype networks of mitochondrial sequence data, but much less so for male-specific markers. Similarly, mitochondrial genetic differentiation was twice as high among populations compared to Y-chromosomal variation. We also found that genetic distance increased faster with geographic distance for mitochondrial than for Y-linked markers, leading to estimates of male dispersal distances that are several-fold higher than those of females. These findings provide evidence for strong male-biased dispersal in orang-utans. The transition to predominantly female-biased dispersal in the great ape lineage appears to be correlated with life in multimale groups and may reflect the associated fitness benefits of reliable male coalitions with relatives or known partners, a feature that is absent in orang-utans., (© 2012 Blackwell Publishing Ltd.)

The heritability of human personality is well-established. Recent research indicates that nonadditive genetic effects, such as dominance and epistasis, play a large role in personality variation. One possible explanation for the latter finding is that there has been recent selection on human personality. To test this possibility, we estimated additive and nonadditive genetic variance in personality and subjective well-being of zoo-housed orangutans. More than half of the genetic variance in these traits could be attributed to nonadditive genetic effects, modeled as dominance. Subjective well-being had genetic overlap with personality, though less so than has been found in humans or chimpanzees. Since a large portion of nonadditive genetic variance in personality is not unique to humans, the nonadditivity of human personality is not sufficient evidence for recent selection of personality in humans. Nonadditive genetic variance may be a general feature of the genetic structure of personality in primates and other animals.

Loss of gene function is implicated in the emergence of novel phenotypes during organism evolution. Here, we report the inactivation of the MSLNL gene encoding mesothelin-like protein in African great ape evolution. Human MSLNL has a nonsense mutation in exon 10 and two polymorphic mutations: a frameshift in exon 3 and a nonsense codon in exon 8. The gorilla gene also shows multiple deleterious mutations, including a premature stop codon, a deletion, and a splice site mutation. Molecular evolutionary analysis indicated relaxed selection pressure on MSLNL in African great ape lineages, which suggested that MSLNL might have become inactivated before the divergence of human, chimpanzee and gorilla. The mouse Mslnl gene is highly expressed in olfactory epithelium and moderately expressed in several other tissues. We propose that the loss of MSLNL may be associated with the evolution of the olfactory system in African great apes including human., (Copyright © 2012 Elsevier B.V. All rights reserved.)

Gorillas are humans' closest living relatives after chimpanzees, and are of comparable importance for the study of human origins and evolution. Here we present the assembly and analysis of a genome sequence for the western lowland gorilla, and compare the whole genomes of all extant great ape genera. We propose a synthesis of genetic and fossil evidence consistent with placing the human-chimpanzee and human-chimpanzee-gorilla speciation events at approximately 6 and 10 million years ago. In 30% of the genome, gorilla is closer to human or chimpanzee than the latter are to each other; this is rarer around coding genes, indicating pervasive selection throughout great ape evolution, and has functional consequences in gene expression. A comparison of protein coding genes reveals approximately 500 genes showing accelerated evolution on each of the gorilla, human and chimpanzee lineages, and evidence for parallel acceleration, particularly of genes involved in hearing. We also compare the western and eastern gorilla species, estimating an average sequence divergence time 1.75 million years ago, but with evidence for more recent genetic exchange and a population bottleneck in the eastern species. The use of the genome sequence in these and future analyses will promote a deeper understanding of great ape biology and evolution.

We present a hidden Markov model (HMM) for inferring gradual isolation between two populations during speciation, modelled as a time interval with restricted gene flow. The HMM describes the history of adjacent nucleotides in two genomic sequences, such that the nucleotides can be separated by recombination, can migrate between populations, or can coalesce at variable time points, all dependent on the parameters of the model, which are the effective population sizes, splitting times, recombination rate, and migration rate. We show by extensive simulations that the HMM can accurately infer all parameters except the recombination rate, which is biased downwards. Inference is robust to variation in the mutation rate and the recombination rate over the sequence and also robust to unknown phase of genomes unless they are very closely related. We provide a test for whether divergence is gradual or instantaneous, and we apply the model to three key divergence processes in great apes: (a) the bonobo and common chimpanzee, (b) the eastern and western gorilla, and (c) the Sumatran and Bornean orang-utan. We find that the bonobo and chimpanzee appear to have undergone a clear split, whereas the divergence processes of the gorilla and orang-utan species occurred over several hundred thousands years with gene flow stopping quite recently. We also apply the model to the Homo/Pan speciation event and find that the most likely scenario involves an extended period of gene flow during speciation., Competing Interests: The authors have declared that no competing interests exist.

Background: Several studies suggested great ape cultures, arguing that human cumulative culture presumably evolved from such a foundation. These focused on conspicuous behaviours, and showed rich geographic variation, which could not be attributed to known ecological or genetic differences. Although geographic variation within call types (accents) has previously been reported for orang-utans and other primate species, we examine geographic variation in the presence/absence of discrete call types (dialects). Because orang-utans have been shown to have geographic variation that is not completely explicable by genetic or ecological factors we hypothesized that this will be similar in the call domain and predict that discrete call type variation between populations will be found., Methodology/principal Findings: We examined long-term behavioural data from five orang-utan populations and collected fecal samples for genetic analyses. We show that there is geographic variation in the presence of discrete types of calls. In exactly the same behavioural context (nest building and infant retrieval), individuals in different wild populations customarily emit either qualitatively different calls or calls in some but not in others. By comparing patterns in call-type and genetic similarity, we suggest that the observed variation is not likely to be explained by genetic or ecological differences., Conclusion/significance: These results are consistent with the potential presence of 'call cultures' and suggest that wild orang-utans possess the ability to invent arbitrary calls, which spread through social learning. These findings differ substantially from those that have been reported for primates before. First, the results reported here are on dialect and not on accent. Second, this study presents cases of production learning whereas most primate studies on vocal learning were cases of contextual learning. We conclude with speculating on how these findings might assist in bridging the gap between vocal communication in non-human primates and human speech.

Although geographic variation in an organism's traits is often seen as a consequence of selection on locally adaptive genotypes accompanied by canalized development [1], developmental plasticity may also play a role [2, 3], especially in behavior [4]. Behavioral plasticity includes both individual learning and social learning of local innovations ("culture"). Cultural plasticity is the undisputed and dominant explanation for geographic variation in human behavior. It has recently also been suggested to hold for various primates and birds [5], but this proposition has been met with widespread skepticism [6-8]. Here, we analyze parallel long-term studies documenting extensive geographic variation in behavioral ecology, social organization, and putative culture of orangutans [9] (genus Pongo). We show that genetic differences among orangutan populations explain only very little of the geographic variation in behavior, whereas environmental differences explain much more, highlighting the importance of developmental plasticity. Moreover, variation in putative cultural variants is explained by neither genetic nor environmental differences, corroborating the cultural interpretation. Thus, individual and cultural plasticity provide a plausible pathway toward local adaptation in long-lived organisms such as great apes and formed the evolutionary foundation upon which human culture was built., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

A human mitochondrial DNA (mtDNA) transition, m.1555A>G, in the 12S rRNA gene causes non-syndromic hearing loss. However, this pathological mutation is the wild-type allele in orangutan mtDNA. Here we rule out different genetic factors as the reason for its fixation in orangutans and show that aminoglycosides negatively affect the oxidative phosphorylation function by decreasing the synthesis of mtDNA-encoded proteins and the amount and activity of respiratory complex IV. These drugs also diminish the growth rate of orangutan cells. The m.1555G nucleotide is also the wild-type allele in other mammal species and they might be at risk of suffering a mitochondrial disorder if treated with aminoglycosides. Therefore, pharmacogenomic approaches should be used to confirm this possibility. These observations are important for human health. Due to the fact that old age and high frequency are criteria widely used in mitochondrial medicine to rule out a genetic change as being a pathological mutation, our results prevent against simplistic genetic approaches that do not consider the potential effect of environmental conditions. Hence, these results suggest that some ancient and highly frequent human population polymorphisms, such as those defining mtDNA haplogroups, in mitochondrial rRNA genes can be deleterious in association with new environmental conditions. Therefore, as the discovery of ribosomal antibiotics has allowed to fight infectious diseases and this breakthrough can be considered an important scientific advance or 'progress', our results suggest that 'progress' can also have a negative counterpart and render detrimental many of these mtDNA genotypes.

Copy number variants (CNVs) are increasingly acknowledged as an important source of evolutionary novelties in the human lineage. However, our understanding of their significance is still hindered by the lack of primate CNV data. We performed intraspecific comparative genomic hybridizations to identify loci harboring copy number variants in each of the four great apes: bonobos, chimpanzees, gorillas, and orangutans. For the first time, we could analyze differences in CNV location and frequency in these four species, and compare them with human CNVs and primate segmental duplication (SD) maps. In addition, for bonobo and gorilla, patterns of CNV and nucleotide diversity were studied in the same individuals. We show that CNVs have been subject to different selective pressures in different lineages. Evidence for purifying selection is stronger in gorilla CNVs overlapping genes, while positive selection appears to have driven the fixation of structural variants in the orangutan lineage. In contrast, chimpanzees and bonobos present high levels of common structural polymorphism, which is indicative of relaxed purifying selection together with the higher mutation rates induced by the known burst of segmental duplication in the ancestor of the African apes. Indeed, the impact of the duplication burst is noticeable by the fact that bonobo and chimpanzee share more CNVs with gorilla than expected. Finally, we identified a number of interesting genomic regions that present high-frequency CNVs in all great apes, while containing only very rare or even pathogenic structural variants in humans.

Background: The substitution rates within different nucleotide contexts are subject to varying levels of bias. The most well known example of such bias is the excess of C to T (C > T) mutations in CpG (CG) dinucleotides. The molecular mechanisms underlying this bias are important factors in human genome evolution and cancer development. The discovery of other nucleotide contexts that have profound effects on substitution rates can improve our understanding of how mutations are acquired, and why mutation hotspots exist., Results: We compared rates of inherited mutations in 1-4 bp nucleotide contexts using reconstructed ancestral states of human single nucleotide polymorphisms (SNPs) from intergenic regions. Chimp and orangutan genomic sequences were used as outgroups. We uncovered 3.5 and 3.3-fold excesses of T > C mutations in the second position of ATTG and ATAG words, respectively, and a 3.4-fold excess of A > C mutations in the first position of the ACAA word., Conclusions: Although all the observed biases are less pronounced than the 5.1-fold excess of C > T mutations in CG dinucleotides, the three 4 bp mutation contexts mentioned above (and their complementary contexts) are well distinguished from all other mutation contexts. This provides a challenge to discover the underlying mechanisms responsible for the observed excesses of mutations.