Evolutionary history inferred from the de novo assembly of a nonmodel organism, the blue-eyed black lemur

Lemurs, the living primates most distantly related to humans, demonstrate incredible diversity in behaviour, life history patterns and adaptive traits. Although many lemur species are endangered within their native Madagascar, there is no high‐quality genome assembly from this taxon, limiting population and conservation genetic studies. One critically endangered lemur is the blue‐eyed black lemur Eulemur flavifrons. This species is fixed for blue irises, a convergent trait that evolved at least four times in primates and was subject to positive selection in humans, where 5′ regulatory variation of OCA2 explains most of the brown/blue eye colour differences. We built a de novo genome assembly for E. flavifrons, providing the most complete lemur genome to date, and a high confidence consensus sequence for close sister species E. macaco, the (brown‐eyed) black lemur. From diversity and divergence patterns across the genomes, we estimated a recent split time of the two species (160 Kya) and temporal fluctuations in effective population sizes that accord with known environmental changes. By looking for regions of unusually low diversity, we identified potential signals of directional selection in E. flavifrons at MITF, a melanocyte development gene that regulates OCA2 and has previously been associated with variation in human iris colour, as well as at several other genes involved in melanin biosynthesis in mammals. Our study thus illustrates how whole‐genome sequencing of a few individuals can illuminate the demographic and selection history of nonmodel species.

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

[2]  L. Duret,et al.  Comparative population genomics in animals uncovers the determinants of genetic diversity , 2014, Nature.

[3]  J. Wall Estimating ancestral population sizes and divergence times. , 2003, Genetics.

[4]  Nicholas G. Martin,et al.  Digital Quantification of Human Eye Color Highlights Genetic Association of Three New Loci , 2010, PLoS genetics.

[5]  R. Mittermeier,et al.  Lemurs of Madagascar : a strategy for their conservation 2013-2016 , 2013 .

[6]  C. Downs,et al.  A Population Estimate of Blue-Eyed Black Lemurs in Ankarafa Forest, Sahamalaza-Iles Radama National Park, Madagascar , 2011, Folia Primatologica.

[7]  Hans Eiberg,et al.  Blue eye color in humans may be caused by a perfectly associated founder mutation in a regulatory element located within the HERC2 gene inhibiting OCA2 expression , 2008, Human Genetics.

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

[9]  Liqun Luo,et al.  A molecular basis for classic blond hair color in Europeans , 2014, Nature Genetics.

[10]  Inge Jonassen,et al.  The genome sequence of Atlantic cod reveals a unique immune system , 2011, Nature.

[11]  M. Kimura Preponderance of synonymous changes as evidence for the neutral theory of molecular evolution , 1977, Nature.

[12]  Molly Przeworski,et al.  The signature of positive selection at randomly chosen loci. , 2002, Genetics.

[13]  Jason E Stajich,et al.  Likelihoods From Summary Statistics: Recent Divergence Between Species , 2005, Genetics.

[14]  J. B. Walsh,et al.  Effects of linkage on rates of molecular evolution. , 1988, Proceedings of the National Academy of Sciences of the United States of America.

[15]  D. Gudbjartsson,et al.  Recombination rate and reproductive success in humans , 2004, Nature Genetics.

[16]  Timothy B Sackton,et al.  Natural Selection Constrains Neutral Diversity across A Wide Range of Species , 2014, bioRxiv.

[17]  Rikke Bagger Jørgensen,et al.  Introgression between oilseed rape (Brassica napus L.) and its weedy relative B. rapa L. in a natural population , 2001, Genetic Resources and Crop Evolution.

[18]  K. Olsen,et al.  Phylogeographic studies in plants: problems and prospects , 1998 .

[19]  D. Fisher,et al.  MITF: master regulator of melanocyte development and melanoma oncogene. , 2006, Trends in molecular medicine.

[20]  C. Groves,et al.  Lemur Diversity in Madagascar , 2008, International Journal of Primatology.

[21]  I. Nolte,et al.  Analysis of systematic effects on congenital sensorineural deafness in German Dalmatian dogs. , 2003, Veterinary journal.

[22]  Brad T. Sherman,et al.  Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists , 2008, Nucleic acids research.

[23]  W. Oetting,et al.  Molecular basis of albinism: Mutations and polymorphisms of pigmentation genes associated with albinism , 1999, Human mutation.

[24]  K. Kidd,et al.  A global view of the OCA2-HERC2 region and pigmentation , 2011, Human Genetics.

[25]  A. Yoder The lemur revolution starts now: the genomic coming of age for a non-model organism. , 2013, Molecular phylogenetics and evolution.

[26]  D. Hartl,et al.  Principles of population genetics , 1981 .

[27]  M. Kreitman,et al.  Adaptive protein evolution at the Adh locus in Drosophila , 1991, Nature.

[28]  David R. Kelley,et al.  Quake: quality-aware detection and correction of sequencing errors , 2010, Genome Biology.

[29]  M. Przeworski,et al.  The Convergent Evolution of Blue Iris Pigmentation in Primates Took Distinct Molecular Paths , 2013, American journal of physical anthropology.

[30]  R. Woychik,et al.  Molecular characterization of the mouse agouti locus , 1992, Cell.

[31]  M. Lynch,et al.  mlRho – a program for estimating the population mutation and recombination rates from shotgun‐sequenced diploid genomes , 2010, Molecular ecology.

[32]  D. Altshuler,et al.  A map of human genome variation from population-scale sequencing , 2010, Nature.

[33]  R. Nielsen,et al.  Quantifying Population Genetic Differentiation from Next-Generation Sequencing Data , 2013, Genetics.

[34]  M. Goddard,et al.  Inferring Demography from Runs of Homozygosity in Whole-Genome Sequence, with Correction for Sequence Errors , 2013, Molecular biology and evolution.

[35]  Wei Chen,et al.  A three-single-nucleotide polymorphism haplotype in intron 1 of OCA2 explains most human eye-color variation. , 2007, American journal of human genetics.

[36]  Snæbjörn Pálsson,et al.  Genetic determinants of hair, eye and skin pigmentation in Europeans , 2007, Nature Genetics.

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

[38]  Matthew Thomas,et al.  Sequences associated with human iris pigmentation. , 2003, Genetics.

[39]  J. Matthews,et al.  Mutations in MITF and PAX3 Cause “Splashed White” and Other White Spotting Phenotypes in Horses , 2012, PLoS Genetics.

[40]  L. Gould Lemurs of Madagascar , 1996, International Journal of Primatology.

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

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

[43]  M. Picardo,et al.  The genetic determination of skin pigmentation: KITLG and the KITLG/c-Kit pathway as key players in the onset of human familial pigmentary diseases. , 2011, The Journal of investigative dermatology.

[44]  T. Korneliussen,et al.  Estimating Individual Admixture Proportions from Next Generation Sequencing Data , 2013, Genetics.

[45]  Páll Melsted,et al.  Comparative RNA sequencing reveals substantial genetic variation in endangered primates. , 2012, Genome research.

[46]  S. O’Brien,et al.  A Molecular Phylogeny of Living Primates , 2011, PLoS genetics.

[47]  A. Gnirke,et al.  High-quality draft assemblies of mammalian genomes from massively parallel sequence data , 2010, Proceedings of the National Academy of Sciences.

[48]  Keith Bradnam,et al.  Assessing the gene space in draft genomes , 2008, Nucleic acids research.

[49]  Eric S. Lander,et al.  Human genome sequence variation and the influence of gene history, mutation and recombination , 2002, Nature Genetics.

[50]  E. Myers,et al.  Basic local alignment search tool. , 1990, Journal of molecular biology.

[51]  M. Lynch,et al.  The critical effective size for a genetically secure population , 1998 .

[52]  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.

[53]  S. Jeandroz,et al.  Comparative genome mapping among Picea glauca, P. mariana × P. rubens and P. abies, and correspondence with other Pinaceae , 2006, Theoretical and Applied Genetics.

[54]  D. Stern Evolution, Development, & the Predictable Genome , 2010 .

[55]  J. Yamagiwa Some External Characters of the Japanese Monkeys (Macaca fuscata) , 1979 .

[56]  M. DePristo,et al.  A framework for variation discovery and genotyping using next-generation DNA sequencing data , 2011, Nature Genetics.

[57]  Michael G. Roth,et al.  Genome-Wide siRNA-Based Functional Genomics of Pigmentation Identifies Novel Genes and Pathways That Impact Melanogenesis in Human Cells , 2008, PLoS genetics.

[58]  J. M. Smith,et al.  The hitch-hiking effect of a favourable gene. , 1974, Genetical research.

[59]  Nicholas G Martin,et al.  A single SNP in an evolutionary conserved region within intron 86 of the HERC2 gene determines human blue-brown eye color. , 2008, American journal of human genetics.

[60]  Jing He,et al.  Peregrine and saker falcon genome sequences provide insights into evolution of a predatory lifestyle , 2013, Nature Genetics.

[61]  P. Andolfatto,et al.  Revisiting an Old Riddle: What Determines Genetic Diversity Levels within Species? , 2012, PLoS biology.

[62]  Simon H. Martin,et al.  Butterfly genome reveals promiscuous exchange of mimicry adaptations among species , 2012, Nature.

[63]  Páll Melsted,et al.  A Genome Sequence Resource for the Aye-Aye (Daubentonia madagascariensis), a Nocturnal Lemur from Madagascar , 2011, Genome biology and evolution.

[64]  Anders Albrechtsen,et al.  ANGSD: Analysis of Next Generation Sequencing Data , 2014, BMC Bioinformatics.

[65]  Shelley D. Smith,et al.  Tietz syndrome (hypopigmentation/deafness) caused by mutation ofMITF , 2000, Journal of medical genetics.

[66]  M. Fumagalli,et al.  Assessing the Effect of Sequencing Depth and Sample Size in Population Genetics Inferences , 2013, PloS one.

[67]  Brad T. Sherman,et al.  Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources , 2008, Nature Protocols.

[68]  N. Mundy,et al.  Blue eyes in lemurs and humans: same phenotype, different genetic mechanism. , 2009, American journal of physical anthropology.

[69]  Y. Rumpler,et al.  Distribution and morphological characters of intermediate forms between the black lemur (Eulemur macaco macaco) and the Sclater's lemur (E. m. flavifrons) , 1991, Primates.

[70]  Sonja J. Prohaska,et al.  Analysis of the African coelacanth genome sheds light on tetrapod evolution , 2013, Nature.

[71]  J. Lupski,et al.  Mutation of FIG4 causes neurodegeneration in the pale tremor mouse and patients with CMT4J , 2007, Nature.

[72]  P. deMenocal,et al.  Plio-Pleistocene African Climate , 1995, Science.

[73]  Robert-Jan Palstra,et al.  HERC2 rs12913832 modulates human pigmentation by attenuating chromatin-loop formation between a long-range enhancer and the OCA2 promoter. , 2012, Genome research.

[74]  Albert J. Vilella,et al.  A high-resolution map of human evolutionary constraint using 29 mammals , 2011, Nature.

[75]  M. Goossens,et al.  Review and update of mutations causing Waardenburg syndrome , 2010, Human mutation.

[76]  M. Przeworski,et al.  A new approach to estimate parameters of speciation models with application to apes. , 2007, Genome research.