Large-scale population genomics of Malayan pangolins reveals deep diversification and a new species

Background Archipelagos and oceanic islands often present high percentage of endemism due to rapid speciation. The Malayan pangolin is a species distributing at both mainland (southern Yunnan, China) and oceanic islands via Malayan peninsula, which may result in deep differentiation among populations. In-depth investigation of population structure and genetic consequences for such species is of vital importance for their protection and conservation, practically for the critically endangered Malayan pangolin that is suffering from poaching, illegal trade, and habitat loss. Results Here we carried out a large-scale population genomic analysis for Malayan pangolins, and revealed three highly distinct genetic populations in this species, two of which are now being reported for the first time. Based on multiple lines of genomic and morphological evidence, we postulate the existence of a new pangolin species (Manis_1). Genetic diversity and recent inbreeding were both at a moderate level for both Malayan pangolins and Manis_1, but mainland Malayan pangolins presented relatively lower genetic diversity, higher inbreeding and fitness cost than island populations. Conclusions We found extremely deep and graded differentiation in Malayan pangolins, with two newly discovered genetic populations and a new pangolin species that is closely related to the Philippine pangolin than the typical Malayan pangolin, but a distant relative of the Indian pangolin. Anthropogenic factors did not significantly weaken the basis of genetic sustainability for Malayan pangolins, but mainland Malayan pangolins should be paid more attention for conservation due to higher genetic risks than island populations.

[1]  Sunil Kumar Sahu,et al.  Whole-genome resequencing of Chinese pangolins reveals a population structure and provides insights into their conservation , 2022, Communications Biology.

[2]  Sunil Kumar Sahu,et al.  Chromosome‐scale genomes reveal genomic consequences of inbreeding in the South China tiger: A comparative study with the Amur tiger , 2022, Molecular ecology resources.

[3]  S. Choo,et al.  Successful captive breeding of a Malayan pangolin population to the third filial generation , 2021, Communications Biology.

[4]  P. Gaubert,et al.  A timely systematic review on pangolin research, commercialization, and popularization to identify knowledge gaps and produce conservation guidelines , 2021 .

[5]  A. Albrechtsen,et al.  High genetic diversity and low differentiation reflect the ecological versatility of the African leopard , 2021, Current Biology.

[6]  Sunil Kumar Sahu,et al.  Chromosome-scale genomes provide new insights into subspecies divergence and evolutionary characteristics of the giant panda. , 2021, Science bulletin.

[7]  Lifeng Zhu,et al.  Genome-wide signatures of mammalian skin covering evolution , 2021, Science China Life Sciences.

[8]  B. Deák,et al.  Linking environmental heterogeneity and plant diversity: The ecological role of small natural features in homogeneous landscapes. , 2020, The Science of the total environment.

[9]  Thomas M. Keane,et al.  Twelve years of SAMtools and BCFtools , 2020, GigaScience.

[10]  C. Roos,et al.  Molecular Genetics Supports a Potential Fifth Asian Pangolin Species (Mammalia, Pholidota, Manis) , 2020, Zoological Science.

[11]  D. Hong Gen‐morph species concept—A new and integrative species concept for outbreeding organisms , 2020 .

[12]  Fei Gao,et al.  CNGBdb: China National GeneBank DataBase. , 2020, Yi chuan = Hereditas.

[13]  Ashutosh Kumar Singh,et al.  Pangolin Indexing System: implications in forensic surveillance of large seizures , 2020, International Journal of Legal Medicine.

[14]  Ya-ping Zhang,et al.  Genomic consequences of population decline in critically endangered pangolins and their demographic histories , 2020, National science review.

[15]  Jirimutu,et al.  Whole-genome sequencing of 128 camels across Asia reveals origin and migration of domestic Bactrian camels , 2020, Communications Biology.

[16]  Fei Gao,et al.  CNSA: a data repository for archiving omics data , 2020, bioRxiv.

[17]  F. Delsuc,et al.  Skull shape variation in extant pangolins (Pholidota: Manidae): allometric patterns and systematic implications , 2019, Zoological Journal of the Linnean Society.

[18]  Ryan W Taylor,et al.  Recent Evolutionary History of Tigers Highlights Contrasting Roles of Genetic Drift and Selection , 2019, bioRxiv.

[19]  Martin Kuhlwilm,et al.  The Genomic Footprints of the Fall and Recovery of the Crested Ibis , 2019, Current Biology.

[20]  N. Schmidt,et al.  The Muskox Lost a Substantial Part of Its Genetic Diversity on Its Long Road to Greenland , 2018, Current Biology.

[21]  J. Flowers,et al.  Origins and geographic diversification of African rice (Oryza glaberrima) , 2018, bioRxiv.

[22]  S. Choo,et al.  Conservation genomics reveals possible illegal trade routes and admixture across pangolin lineages in Southeast Asia , 2018, Conservation Genetics.

[23]  Sudhir Kumar,et al.  MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. , 2018, Molecular biology and evolution.

[24]  A. Antunes,et al.  The Complete Phylogeny of Pangolins: Scaling Up Resources for the Molecular Tracing of the Most Trafficked Mammals on Earth , 2018, The Journal of heredity.

[25]  S. Pfeifer,et al.  The demographic and adaptive history of the African green monkey , 2017, bioRxiv.

[26]  Robert Lanfear,et al.  PartitionFinder 2: New Methods for Selecting Partitioned Models of Evolution for Molecular and Morphological Phylogenetic Analyses. , 2016, Molecular biology and evolution.

[27]  J. Ross,et al.  Where did all the pangolins go? International CITES trade in pangolin species , 2016 .

[28]  Guat Jah Wong,et al.  Pangolin genomes and the evolution of mammalian scales and immunity , 2016, Genome research.

[29]  Diego Ortega-Del Vecchyo,et al.  Genomic Flatlining in the Endangered Island Fox , 2016, Current Biology.

[30]  M. Stephens,et al.  New evidence for hybrid zones of forest and savanna elephants in Central and West Africa , 2015, Molecular Ecology.

[31]  David Bryant,et al.  popart: full‐feature software for haplotype network construction , 2015 .

[32]  Shaoshan Li,et al.  Reproductive parameters of the Sunda pangolin, Manis javanica , 2015, Folia Zoologica.

[33]  D. Reich,et al.  Genetic evidence for two founding populations of the Americas , 2015, Nature.

[34]  P. Gaubert,et al.  Molecular tracing of confiscated pangolin scales for conservation and illegal trade monitoring in Southeast Asia , 2015 .

[35]  F. Hailer Introgressive hybridization: brown bears as vectors for polar bear alleles , 2015, Molecular ecology.

[36]  I. Stirling,et al.  Genomic evidence of geographically widespread effect of gene flow from polar bears into brown bears , 2015, Molecular ecology.

[37]  A. von Haeseler,et al.  IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies , 2014, Molecular biology and evolution.

[38]  Q. Weng,et al.  Fauna in decline: plight of the pangolin. , 2014, Science.

[39]  H. Burbano,et al.  The Genomic Architecture of Population Divergence between Subspecies of the European Rabbit , 2014, PLoS genetics.

[40]  S. Harris,et al.  High Risks of Losing Genetic Diversity in an Endemic Mauritian Gecko: Implications for Conservation , 2014, PloS one.

[41]  R. Durbin,et al.  Inferring human population size and separation history from multiple genome sequences , 2014, Nature Genetics.

[42]  D. Macdonald,et al.  Scaling up pangolin protection in China , 2014 .

[43]  B. Browning,et al.  Improving the Accuracy and Efficiency of Identity-by-Descent Detection in Population Data , 2013, Genetics.

[44]  Swapan Mallick,et al.  Ancient Admixture in Human History , 2012, Genetics.

[45]  J. Losos,et al.  The island–mainland species turnover relationship , 2012, Proceedings of the Royal Society B: Biological Sciences.

[46]  Maxim Teslenko,et al.  MrBayes 3.2: Efficient Bayesian Phylogenetic Inference and Model Choice Across a Large Model Space , 2012, Systematic biology.

[47]  Gonçalo R. Abecasis,et al.  The variant call format and VCFtools , 2011, Bioinform..

[48]  James Rosindell,et al.  A unified model of island biogeography sheds light on the zone of radiation. , 2011, Ecology letters.

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

[50]  C. Klingenberg MorphoJ: an integrated software package for geometric morphometrics , 2011, Molecular ecology resources.

[51]  P. Visscher,et al.  GCTA: a tool for genome-wide complex trait analysis. , 2011, American journal of human genetics.

[52]  Harry Campbell,et al.  Genomic Runs of Homozygosity Record Population History and Consanguinity , 2010, PloS one.

[53]  Josyf Mychaleckyj,et al.  Robust relationship inference in genome-wide association studies , 2010, Bioinform..

[54]  M. Pärtel,et al.  Environmental heterogeneity, species diversity and co‐existence at different spatial scales , 2010 .

[55]  H. Hakonarson,et al.  ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data , 2010, Nucleic acids research.

[56]  Richard Durbin,et al.  Fast and accurate long-read alignment with Burrows–Wheeler transform , 2010, Bioinform..

[57]  David H. Alexander,et al.  Fast model-based estimation of ancestry in unrelated individuals. , 2009, Genome research.

[58]  T. Gaudin,et al.  The Phylogeny of Living and Extinct Pangolins (Mammalia, Pholidota) and Associated Taxa: A Morphology Based Analysis , 2009, Journal of Mammalian Evolution.

[59]  V. Lowe,et al.  Hybridization between Red deer (Cervus elaphus) and Sika deer (Cervus nippon) with particular reference to stocks in N.W. England , 2009 .

[60]  W. Barthlott,et al.  A global assessment of endemism and species richness across island and mainland regions , 2009, Proceedings of the National Academy of Sciences.

[61]  D. Strayer Implications for Conservation , 2008 .

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

[63]  A. Antunes,et al.  ASSESSING THE TAXONOMIC STATUS OF THE PALAWAN PANGOLIN MANIS CULIONENSIS (PHOLIDOTA) USING DISCRETE MORPHOLOGICAL CHARACTERS , 2005 .

[64]  Q. Wan,et al.  A NEW SUBSPECIES OF GIANT PANDA (Ailuropoda melanoleuca) FROM SHAANXI, CHINA , 2005 .

[65]  R. Frankham,et al.  Most species are not driven to extinction before genetic factors impact them. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[66]  A. Meyer,et al.  SHAPE ANALYSIS OF SYMMETRIC STRUCTURES: QUANTIFYING VARIATION AMONG INDIVIDUALS AND ASYMMETRY , 2002, Evolution; international journal of organic evolution.

[67]  L. Keller,et al.  Inbreeding effects in wild populations. , 2002 .

[68]  D. Schluter,et al.  Analysis of an evolutionary species–area relationship , 2000, Nature.

[69]  R. Mittermeier,et al.  Biodiversity hotspots for conservation priorities , 2000, Nature.

[70]  D. Roff,et al.  A COMPARISON OF INBREEDING DEPRESSION IN LIFE‐HISTORY AND MORPHOLOGICAL TRAITS IN ANIMALS , 1999, Evolution; international journal of organic evolution.

[71]  N. Barton,et al.  Introgression through rare hybridization: A genetic study of a hybrid zone between red and sika deer (genus Cervus) in Argyll, Scotland. , 1999, Genetics.

[72]  K. S. Kim,et al.  The complete nucleotide sequence of the domestic dog (Canis familiaris) mitochondrial genome. , 1998, Molecular phylogenetics and evolution.

[73]  Christian Peter Klingenberg,et al.  GEOMETRIC MORPHOMETRICS OF DEVELOPMENTAL INSTABILITY: ANALYZING PATTERNS OF FLUCTUATING ASYMMETRY WITH PROCRUSTES METHODS , 1998, Evolution; international journal of organic evolution.

[74]  J. V. López,et al.  Complete nucleotide sequences of the domestic cat (Felis catus) mitochondrial genome and a transposed mtDNA tandem repeat (Numt) in the nuclear genome. , 1996, Genomics.

[75]  M. Donoghue,et al.  PHYLOGENETIC SYSTEMATICS AND SPECIES REVISITED , 1990 .

[76]  F. Rohlf,et al.  Extensions of the Procrustes Method for the Optimal Superimposition of Landmarks , 1990 .

[77]  T. Bonebrake,et al.  Meat and medicine: historic and contemporary use in Asia , 2020 .

[78]  J. Jimerson,et al.  The rescue, rehabilitation and release of pangolins , 2020 .

[79]  Zhigang Jiang,et al.  China’s mammal diversity (2nd edition) , 2017 .

[80]  Wen-Hsiung Li,et al.  Nonrandomness of point mutation as reflected in nucleotide substitutions in pseudogenes and its evolutionary implications , 2005, Journal of Molecular Evolution.

[81]  Nai-fa Liu,et al.  Measurement and comparison of skull variables in Chinese pangolin and Malayan pangolin , 2004 .