Population Genomic Analysis Reveals a Rich Speciation and Demographic History of Orang-utans (Pongo pygmaeus and Pongo abelii)

18Citations
Citations of this article
83Readers
Mendeley users who have this article in their library.

Abstract

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. © 2013 Ma et al.

Figures

  • Figure 1. Phylogenetic tree of mtDNA HVRI region and allele sharing matrix of mtDNA and autosomes. (A) The phylogenetic tree among human, chimpanzee, Bornean (red) and Sumatran (blue). The nodes shaded in cyan represent the grouping of the three Sumatran individuals that has the furthest distance to the rest clustering of the Sumatran group. The 10 re-sequenced individuals from our study are colored as red (Bornean) or blue (Sumatran). The corresponding geographic origins are shown on the right with the following annotation (SU: Sumatran; TU/SL: Tuanan or Sungai Lading; DV/K: Danum Valley or Kinabatangan; SA/GP: Sabangau or Gunung Palung; EK: East Kalimantan including Kutai national park or Sangatta; SAR/DS: Semongok Wildlife Rehabilitation Centre or Danau Sentarum). (B) The IBS allele sharing matrix among the 10 individuals from our study with the same group coloring scheme. The upper and lower triangular matrix represents the IBS sharing of all mtDNA variation loci (n = 1084) and autosomal variation loci (n = 11,866,619), respectively. doi:10.1371/journal.pone.0077175.g001
  • Figure 2. SFS analysis. (A) SFS for all genomic SNPs (Genome), synonymous (S) or non-synonymous (NS) SNPs in Bornean (Red) and Sumatran (Blue) population. (B) Two-dimensional SFS (Bornean on Y axis and Sumatran on X axis). Top row: Spectra for different functional categories of SNPs. To make all the spectra directly comparable, the non-coding and synonymous spectra have been re-scaled to represent the same number of segregating SNPs as the non-synonymous spectrum. Bottom row: residuals between pairs of spectra. Bottom-left compares synonymous and noncoding spectra, and bottom-right compares non-synonymous and synonymous spectra. Red and blue entries indicate, respectively, that the first spectrum has greater or fewer SNPs in that entry than the second. doi:10.1371/journal.pone.0077175.g002
  • Table 1. Summary statistics of McDonald-Kreitman table entries for Bornean and Sumatran population using Human (hg19) as outgroup.
  • Figure 3. Distribution of the proportion SNPs within each of the PolyPhen categories based on PolyPhen-2 result. Error bars denote 95% confidence intervals on the proportion of SNPs in each category. doi:10.1371/journal.pone.0077175.g003
  • Figure 4. Inferred cumulative distribution of selection coefficients. Solid lines are for the full demographic model, dashed are for the model with no migration, and dotted are for the model with no growth. The break in the axes indicates that the scale is not continuous, because a logarithmic scale cannot span zero. doi:10.1371/journal.pone.0077175.g004
  • Table 2. Percent of mutations within various selective classes for different models of the distribution of selection coefficients and different underlying demographic models.
  • Table 3. Genes under positive selection identified by SnIPRE.

References Powered by Scopus

A Simple, Fast, and Accurate Algorithm to Estimate Large Phylogenies by Maximum Likelihood

14964Citations
N/AReaders
Get full text

BEAST: Bayesian evolutionary analysis by sampling trees

10962Citations
N/AReaders
Get full text

A method and server for predicting damaging missense mutations

10638Citations
N/AReaders
Get full text

Cited by Powered by Scopus

Inference of the distribution of selection coefficients for new nonsynonymous mutations using large samples

120Citations
N/AReaders
Get full text

Accelerated diversification of nonhuman primate malarias in Southeast Asia: Adaptive radiation or geographic speciation?

72Citations
N/AReaders
Get full text

Two to tango: Co-evolution of hominid natural killer cell receptors and MHC

39Citations
N/AReaders
Get full text

Register to see more suggestions

Mendeley helps you to discover research relevant for your work.

Already have an account?

Cite

CITATION STYLE

APA

Ma, X., Kelley, J. L., Eilertson, K., Musharoff, S., Degenhardt, J. D., Martins, A. L., … Bustamante, C. D. (2013). Population Genomic Analysis Reveals a Rich Speciation and Demographic History of Orang-utans (Pongo pygmaeus and Pongo abelii). PLoS ONE, 8(10). https://doi.org/10.1371/journal.pone.0077175

Readers over time

‘13‘14‘15‘16‘17‘18‘19‘20‘21‘22‘23‘2405101520

Readers' Seniority

Tooltip

PhD / Post grad / Masters / Doc 40

71%

Researcher 10

18%

Professor / Associate Prof. 4

7%

Lecturer / Post doc 2

4%

Readers' Discipline

Tooltip

Agricultural and Biological Sciences 49

78%

Biochemistry, Genetics and Molecular Bi... 11

17%

Mathematics 2

3%

Social Sciences 1

2%

Save time finding and organizing research with Mendeley

Sign up for free
0