Whole-genome sequencing of the snub-nosed monkey provides insights into folivory and evolutionary history

Colobines are a unique group of Old World monkeys that principally eat leaves and seeds rather than fruits and insects. We report the sequencing at 146× coverage, de novo assembly and analyses of the genome of a male golden snub-nosed monkey (Rhinopithecus roxellana) and resequencing at 30× coverage of three related species (Rhinopithecus bieti, Rhinopithecus brelichi and Rhinopithecus strykeri). Comparative analyses showed that Asian colobines have an enhanced ability to derive energy from fatty acids and to degrade xenobiotics. We found evidence for functional evolution in the colobine RNASE1 gene, encoding a key secretory RNase that digests the high concentrations of bacterial RNA derived from symbiotic microflora. Demographic reconstructions indicated that the profile of ancient effective population sizes for R. roxellana more closely resembles that of giant panda rather than its congeners. These findings offer new insights into the dietary adaptations and evolutionary history of colobine primates.

[1]  Joel Dudley,et al.  MEGA: A biologist-centric software for evolutionary analysis of DNA and protein sequences , 2008, Briefings Bioinform..

[2]  M. Hicks,et al.  Expression of type 2 cystatin genes CST1-CST5 in adult human tissues and the developing submandibular gland. , 2002, DNA and cell biology.

[3]  W. Miller,et al.  Polar and brown bear genomes reveal ancient admixture and demographic footprints of past climate change , 2012, Proceedings of the National Academy of Sciences.

[4]  Z. An,et al.  Late Pliocene-Pleistocene changes in mass accumulation rates of eolian deposits on the central Chine , 2005 .

[5]  P. Bork,et al.  iPath: interactive exploration of biochemical pathways and networks. , 2008, Trends in biochemical sciences.

[6]  William A. Walters,et al.  QIIME allows analysis of high-throughput community sequencing data , 2010, Nature Methods.

[7]  A. Hobolth,et al.  Estimating Divergence Time and Ancestral Effective Population Size of Bornean and Sumatran Orangutan Subspecies Using a Coalescent Hidden Markov Model , 2011, PLoS genetics.

[8]  J. Rougemont,et al.  A rapid bootstrap algorithm for the RAxML Web servers. , 2008, Systematic biology.

[9]  Ziheng Yang PAML 4: phylogenetic analysis by maximum likelihood. , 2007, Molecular biology and evolution.

[10]  J. Oates,et al.  Colobine Monkeys: Their Ecology, Behaviour and Evolution , 1995 .

[11]  Alexandros Stamatakis,et al.  RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models , 2006, Bioinform..

[12]  Hong-chun Li,et al.  High-resolution climate variability of southwest China during 57–70 ka reflected in a stalagmite δ18O record from Xinya Cave , 2007 .

[13]  R. Durbin,et al.  Inference of human population history from individual whole-genome sequences. , 2011, Nature.

[14]  Z. Yang,et al.  Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. , 1998, Molecular biology and evolution.

[15]  Damian Szklarczyk,et al.  eggNOG v2.0: extending the evolutionary genealogy of genes with enhanced non-supervised orthologous groups, species and functional annotations , 2009, Nucleic Acids Res..

[16]  Robert C. Edgar,et al.  UPARSE: highly accurate OTU sequences from microbial amplicon reads , 2013, Nature Methods.

[17]  M. Nei,et al.  Accuracies of ancestral amino acid sequences inferred by the parsimony, likelihood, and distance methods , 2009, Journal of Molecular Evolution.

[18]  Baoguo Li,et al.  Extinction of Snub-Nosed Monkeys in China During the Past 400 Years , 2002, International Journal of Primatology.

[19]  F. Wei,et al.  Current status and conservation of the gray snub-nosed monkey Rhinopithecus brelichi (Colobinae) in Guizhou, China. , 2009 .

[20]  B. Zheng,et al.  The relationship between climate change and Quaternary glacial cycles on the Qinghai–Tibetan Plateau: review and speculation , 2002 .

[21]  M. Nei,et al.  Extensive Gains and Losses of Olfactory Receptor Genes in Mammalian Evolution , 2007, PloS one.

[22]  Steven Henikoff,et al.  SIFT: predicting amino acid changes that affect protein function , 2003, Nucleic Acids Res..

[23]  P. Bork,et al.  A method and server for predicting damaging missense mutations , 2010, Nature Methods.

[24]  N. Jablonski,et al.  Evolutionary History of the Odd-Nosed Monkeys and the Phylogenetic Position of the Newly Described Myanmar Snub-Nosed Monkey Rhinopithecus strykeri , 2012, PloS one.

[25]  J.-F. Cheng,et al.  Adaptation to herbivory by the Tammar wallaby includes bacterial and glycoside hydrolase profiles different from other herbivores , 2010, Proceedings of the National Academy of Sciences.

[26]  K. Milton Physiological Ecology of Howlers (Alouatta): Energetic and Digestive Considerations and Comparison with the Colobinae , 1998, International Journal of Primatology.

[27]  Steven Salzberg,et al.  BIOINFORMATICS ORIGINAL PAPER , 2004 .

[28]  Bernard B. Suh,et al.  Reconstructing contiguous regions of an ancestral genome. , 2006, Genome research.

[29]  U. Francke,et al.  Molecular cloning of cDNAs encoding rat and human medium-chain acyl-CoA dehydrogenase and assignment of the gene to human chromosome 1. , 1986, Proceedings of the National Academy of Sciences of the United States of America.

[30]  Chankyu Park,et al.  The complete swine olfactory subgenome: expansion of the olfactory gene repertoire in the pig genome , 2012, BMC Genomics.

[31]  Huanming Yang,et al.  Whole-genome sequencing of giant pandas provides insights into demographic history and local adaptation , 2012, Nature Genetics.

[32]  M. Harris,et al.  Food and evolution : toward a theory of human food habits , 1987 .

[33]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[34]  S. Pääbo,et al.  Loss of Olfactory Receptor Genes Coincides with the Acquisition of Full Trichromatic Vision in Primates , 2004, PLoS biology.

[35]  Manuel A. R. Ferreira,et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. , 2007, American journal of human genetics.

[36]  Thierry Heidmann,et al.  Human LINE retrotransposons generate processed pseudogenes , 2000, Nature Genetics.

[37]  R. S. W. van de Wal,et al.  North American ice-sheet dynamics and the onset of 100,000-year glacial cycles , 2008, Nature.

[38]  R. Nielsen,et al.  Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular level. , 2005, Molecular biology and evolution.

[39]  Robert C. Edgar,et al.  MUSCLE: multiple sequence alignment with high accuracy and high throughput. , 2004, Nucleic acids research.

[40]  Jianzhi Zhang Parallel adaptive origins of digestive RNases in Asian and African leaf monkeys , 2006, Nature Genetics.

[41]  R. Durbin,et al.  GeneWise and Genomewise. , 2004, Genome research.

[42]  K. Dyer,et al.  Eosinophil Cationic Protein and Eosinophil-derived Neurotoxin , 1995, Journal of Biological Chemistry.

[43]  Arcadi Navarro,et al.  Great ape genetic diversity and population history , 2013, Nature.

[44]  C. Wadelius,et al.  Spectrum of mutations and sequence variants in the FALDH gene in patients with Sjögren‐Larsson syndrome , 1998, Human mutation.

[45]  Doron Lancet,et al.  Prediction of the odorant binding site of olfactory receptor proteins by human–mouse comparisons , 2004, Protein science : a publication of the Protein Society.

[46]  W. F. Fricke,et al.  Immune status, antibiotic medication and pH are associated with changes in the stomach fluid microbiota , 2013, The ISME Journal.

[47]  Ya-ping Zhang,et al.  Adaptive evolution of digestive RNASE1 genes in leaf-eating monkeys revisited: new insights from ten additional colobines. , 2010, Molecular biology and evolution.

[48]  Vanessa E. Gray,et al.  Evolutionary diagnosis method for variants in personal exomes , 2012, Nature Methods.

[49]  K. Südekum,et al.  First identification of tannin‐binding proteins in saliva of Papio hamadryas using MS/MS mass spectrometry , 2011, American journal of primatology.

[50]  Albert J. Vilella,et al.  Insights into hominid evolution from the gorilla genome sequence , 2012, Nature.

[51]  Koichiro Tamura,et al.  Estimating divergence times in large molecular phylogenies , 2012, Proceedings of the National Academy of Sciences.

[52]  P. Deininger Alu elements: know the SINEs , 2011, Genome Biology.

[53]  A. Hobolth,et al.  Genomic Relationships and Speciation Times of Human, Chimpanzee, and Gorilla Inferred from a Coalescent Hidden Markov Model , 2006, PLoS genetics.

[54]  Christina Freytag,et al.  Colobine Monkeys Their Ecology Behaviour And Evolution , 2016 .

[55]  J. Dudley,et al.  Phylomedicine: an evolutionary telescope to explore and diagnose the universe of disease mutations. , 2011, Trends in genetics : TIG.

[56]  M. Raymo,et al.  A Pliocene‐Pleistocene stack of 57 globally distributed benthic δ18O records , 2005 .

[57]  Jianzhi Zhang,et al.  Adaptive evolution of a duplicated pancreatic ribonuclease gene in a leaf-eating monkey , 2002, Nature Genetics.

[58]  Qi Wu,et al.  Evidence of cellulose metabolism by the giant panda gut microbiome , 2011, Proceedings of the National Academy of Sciences.

[59]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[60]  M. Borodovsky,et al.  Ab initio gene identification in metagenomic sequences , 2010, Nucleic acids research.

[61]  Huanming Yang,et al.  De novo assembly of human genomes with massively parallel short read sequencing. , 2010, Genome research.

[62]  Natalia N. Ivanova,et al.  Metagenomic and functional analysis of hindgut microbiota of a wood-feeding higher termite , 2007, Nature.

[63]  Yang Zhong,et al.  CRDB: Database of Chemosensory Receptor Gene Families in Vertebrate , 2012, PloS one.

[64]  T Hashimoto,et al.  Novel fatty acid beta-oxidation enzymes in rat liver mitochondria. II. Purification and properties of enoyl-coenzyme A (CoA) hydratase/3-hydroxyacyl-CoA dehydrogenase/3-ketoacyl-CoA thiolase trifunctional protein. , 1992, The Journal of biological chemistry.

[65]  M. Ruggero,et al.  Similarity of Traveling-Wave Delays in the Hearing Organs of Humans and Other Tetrapods , 2007, Journal for the Association for Research in Otolaryngology.

[66]  R. Holt,et al.  Duplication and divergence of 2 distinct pancreatic ribonuclease genes in leaf-eating African and Asian colobine monkeys. , 2006, Molecular biology and evolution.