Revealing hidden species diversity in closely related species using nuclear SNPs, SSRs and DNA sequences – a case study in the tree genus Milicia

BackgroundSpecies delimitation in closely related plant taxa can be challenging because (i) reproductive barriers are not always congruent with morphological differentiation, (ii) use of plastid sequences might lead to misinterpretation, (iii) rare species might not be sampled. We revisited molecular-based species delimitation in the African genus Milicia, currently divided into M. regia (West Africa) and M. excelsa (from West to East Africa). We used 435 samples collected in West, Central and East Africa. We genotyped SNP and SSR loci to identify genetic clusters, and sequenced two plastid regions (psbA-trnH, trnC-ycf6) and a nuclear gene (At103) to confirm species’ divergence and compare species delimitation methods. We also examined whether ecological niche differentiation was congruent with sampled genetic structure.ResultsWest African M. regia, West African and East African M. excelsa samples constituted three well distinct genetic clusters according to SNPs and SSRs. In Central Africa, two genetic clusters were consistently inferred by both types of markers, while a few scattered samples, sympatric with the preceding clusters but exhibiting leaf traits of M. regia, were grouped with the West African M. regia cluster based on SNPs or formed a distinct cluster based on SSRs. SSR results were confirmed by sequence data from the nuclear region At103 which revealed three distinct ‘Fields For Recombination’ corresponding to (i) West African M. regia, (ii) Central African samples with leaf traits of M. regia, and (iii) all M. excelsa samples. None of the plastid sequences provide indication of distinct clades of the three species-like units. Niche modelling techniques yielded a significant correlation between niche overlap and genetic distance.ConclusionsOur genetic data suggest that three species of Milicia could be recognized. It is surprising that the occurrence of two species in Central Africa was not reported for this well-known timber tree. Globally, our work highlights the importance of collecting samples in a systematic way and the need for combining different nuclear markers when dealing with species complexes. Recognizing cryptic species is particularly crucial for economically exploited species because some hidden taxa might actually be endangered as they are merged with more abundant species.

[1]  A. Couloux,et al.  Haplowebs as a graphical tool for delimiting species: a revival of Doyle's "field for recombination" approach and its application to the coral genus Pocillopora in Clipperton , 2010, BMC Evolutionary Biology.

[2]  O. Hardy,et al.  spagedi: a versatile computer program to analyse spatial genetic structure at the individual or population levels , 2002 .

[3]  S. Cavers,et al.  Plio-Pleistocene history and phylogeography of Acacia senegal in dry woodlands and savannahs of sub-Saharan tropical Africa: evidence of early colonisation and recent range expansion , 2012, Heredity.

[4]  Bryan C. Carstens,et al.  How to fail at species delimitation , 2013, Molecular ecology.

[5]  S. Cherny,et al.  Utility of the trnH–psbA Intergenic Spacer Region and Its Combinations as Plant DNA Barcodes: A Meta-Analysis , 2012, PloS one.

[6]  J. Doebley,et al.  Rate and pattern of mutation at microsatellite loci in maize. , 2002, Molecular biology and evolution.

[7]  M. Feldman,et al.  Phylogenetic resolution power of microsatellites and various single-nucleotide polymorphism types assessed in 10 divergent chicken populations. , 2014, Animal genetics.

[8]  D. Pollock,et al.  Launching microsatellites: a review of mutation processes and methods of phylogenetic interference. , 1997, The Journal of heredity.

[9]  R. Ennos,et al.  Next-generation hybridization and introgression , 2011, Heredity.

[10]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[11]  Olivier François,et al.  Bayesian clustering algorithms ascertaining spatial population structure: a new computer program and a comparison study , 2007 .

[12]  Jean‐François Flot,et al.  Delimiting Species-Poor Data Sets using Single Molecular Markers: A Study of Barcode Gaps, Haplowebs and GMYC. , 2015, Systematic biology.

[13]  L. Katz,et al.  Diversity of diversity: conceptual and methodological differences in biodiversity estimates of eukaryotic microbes as compared to bacteria. , 2014, Trends in microbiology.

[14]  G. White.  The Place of Morphological * Studies in the Investigation of AnopheZes Species Complexes , 1977 .

[15]  H. Reeve,et al.  Estimating effective paternity number in social insects and the effective number of alleles in a population , 2003, Molecular ecology.

[16]  Brett R. Scheffers,et al.  What we know and don't know about Earth's missing biodiversity. , 2012, Trends in ecology & evolution.

[17]  Vivienne L. Williams,et al.  Phylogeography of the Afromontane Prunus africana reveals a former migration corridor between East and West African highlands , 2011, Molecular ecology.

[18]  Michael Balke,et al.  The Effect of Geographical Scale of Sampling on DNA Barcoding , 2012, Systematic biology.

[19]  B. Kempenaers,et al.  Heterozygosity–fitness correlations in zebra finches: microsatellite markers can be better than their reputation , 2012, Molecular ecology.

[20]  R. Petit,et al.  Blind population genetics survey of tropical rainforest trees , 2006, Molecular ecology.

[21]  Hans Peter Linder,et al.  Plant diversity and endemism in sub‐Saharan tropical Africa , 2001 .

[22]  D. McKey,et al.  Genetic diversity of Tamarindus indica populations: Any clues on the origin from its current distribution? , 2007 .

[23]  J. Mullikin,et al.  Microsatellites are molecular clocks that support accurate inferences about history. , 2009, Molecular biology and evolution.

[24]  R. Henry,et al.  Microsatellites retain phylogenetic signals across genera in eucalypts (Myrtaceae) , 2007 .

[25]  Development of SNP markers for Ayous (Triplochiton scleroxylon K. Schum) an economically important tree species from tropical West and Central Africa , 2016, Conservation Genetics Resources.

[26]  M. Nei,et al.  POPTREE2: Software for constructing population trees from allele frequency data and computing other population statistics with Windows interface. , 2010, Molecular biology and evolution.

[27]  D. Khasa,et al.  Novel microsatellite DNA markers for the threatened African endemic tree species, Milicia excelsa (Moraceae), and cross‐species amplification in Milicia regia , 2006 .

[28]  T. Schoener.  The Anolis Lizards of Bimini: Resource Partitioning in a Complex Fauna , 1968 .

[29]  A. Furman,et al.  Cryptic diversity of the bent-wing bat, Miniopterus schreibersii (Chiroptera: Vespertilionidae), in Asia Minor , 2010, BMC Evolutionary Biology.

[30]  H. Linder,et al.  Species delimitation and relationships: The dance of the seven veils , 2015 .

[31]  Steven J. Phillips,et al.  Sample selection bias and presence-only distribution models: implications for background and pseudo-absence data. , 2009, Ecological applications : a publication of the Ecological Society of America.

[32]  D. Janzen,et al.  Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[33]  L. Kruckenhauser,et al.  Species inflation and taxonomic artefacts—A critical comment on recent trends in mammalian classification , 2013 .

[34]  B. Husband,et al.  RAPID ADAPTIVE DIVERGENCE IN NEW WORLD ACHILLEA, AN AUTOPOLYPLOID COMPLEX OF ECOLOGICAL RACES , 2008, Evolution; international journal of organic evolution.

[35]  David J. Lohman,et al.  Cryptic species as a window on diversity and conservation. , 2007, Trends in ecology & evolution.

[36]  D. J. Funk,et al.  Species-Level Paraphyly and Polyphyly: Frequency, Causes, and Consequences, with Insights from Animal Mitochondrial DNA , 2003 .

[37]  Q. Wheeler,et al.  The perils of DNA barcoding and the need for integrative taxonomy. , 2005, Systematic biology.

[38]  H. Saedler,et al.  Development of COS genes as universally amplifiable markers for phylogenetic reconstructions of closely related plant species , 2008 .

[39]  L. Rieseberg,et al.  Phylogenetic consequences of cytoplasmic gene flow in plants. , 1991 .

[40]  M. Möller,et al.  DNA barcoding for the discrimination of Eurasian yews (Taxus L., Taxaceae) and the discovery of cryptic species , 2011, Molecular ecology resources.

[41]  H. Bandelt,et al.  Median-joining networks for inferring intraspecific phylogenies. , 1999, Molecular biology and evolution.

[42]  L. Excoffier,et al.  Gene flow and species delimitation. , 2009, Trends in ecology & evolution.

[43]  Noah A. Rosenberg,et al.  CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure , 2007, Bioinform..

[44]  G. Dauby,et al.  Comparative phylogeography of eight herbs and lianas (Marantaceae) in central African rainforests , 2014, Front. Genet..

[45]  A. Drummond,et al.  Bayesian Inference of Species Trees from Multilocus Data , 2009, Molecular biology and evolution.

[46]  A. Lambert,et al.  ABGD, Automatic Barcode Gap Discovery for primary species delimitation , 2012, Molecular ecology.

[47]  J. L. Parra,et al.  Very high resolution interpolated climate surfaces for global land areas , 2005 .

[48]  A. Yoder,et al.  Molecular Evidence of Reproductive Isolation in Sympatric Sibling Species of Mouse Lemurs , 2002, International Journal of Primatology.

[49]  C. Cruaud,et al.  How Effective Are DNA Barcodes in the Identification of African Rainforest Trees? , 2013, PloS one.

[50]  C. Hutter,et al.  Microsatellite variation in Drosophila melanogaster and Drosophila simulans: a reciprocal test of the ascertainment bias hypothesis. , 1998, Molecular biology and evolution.

[51]  B. Hansson,et al.  Do microsatellites reflect genome‐wide genetic diversity in natural populations? A comment on Väli et al. (2008) , 2010, Molecular ecology.

[52]  D. J. Funk,et al.  Divergent selection and heterogeneous genomic divergence , 2009, Molecular ecology.

[53]  J. Merilä,et al.  Characterizing genic and nongenic molecular markers: comparison of microsatellites and SNPs , 2013, Molecular ecology resources.

[54]  Martin Lascoux,et al.  Cryptic population genetic structure: the number of inferred clusters depends on sample size , 2010, Molecular ecology resources.

[55]  M W Feldman,et al.  An evaluation of genetic distances for use with microsatellite loci. , 1994, Genetics.

[56]  A. M. Noreen,et al.  Assessing hidden species diversity in the coral Pocillopora damicornis from Eastern Australia , 2012, Coral Reefs.

[57]  G. Mahy,et al.  Identité et écologie des espèces forestières commerciales d'Afrique Centrale: le cas de Milicia spp. , 2012 .

[58]  Jean‐François Flot.  champuru 1.0: a computer software for unraveling mixtures of two DNA sequences of unequal lengths , 2007 .

[59]  J. Doucet,et al.  Phenological patterns in a natural population of a tropical timber tree species, Milicia excelsa (Moraceae): Evidence of isolation by time and its interaction with feeding strategies of dispersers. , 2012, American journal of botany.

[60]  R. Petit,et al.  Current trends in microsatellite genotyping , 2011, Molecular ecology resources.

[61]  J. McKay,et al.  Adaptive population divergence: markers, QTL and traits , 2002 .

[62]  M. Sykes,et al.  Methods and uncertainties in bioclimatic envelope modelling under climate change , 2006 .

[63]  S. Tishkoff,et al.  SNP ascertainment bias in population genetic analyses: Why it is important, and how to correct it , 2013, BioEssays : news and reviews in molecular, cellular and developmental biology.

[64]  J. Doucet,et al.  Spatial genetic structure in Milicia excelsa (Moraceae) indicates extensive gene dispersal in a low‐density wind‐pollinated tropical tree , 2009, Molecular ecology.

[65]  E. Nevo,et al.  Microsatellites: genomic distribution, putative functions and mutational mechanisms: a review , 2002, Molecular ecology.

[66]  Pablo Librado,et al.  DnaSP v5: a software for comprehensive analysis of DNA polymorphism data , 2009, Bioinform..

[67]  J. Doyle.  The Irrelevance of Allele Tree Topologies for Species Delimitation, and a Non-Topological Alternative , 1995 .

[68]  Robert P. Anderson,et al.  Maximum entropy modeling of species geographic distributions , 2006 .

[69]  G. Luikart,et al.  SNPs in ecology, evolution and conservation , 2004 .

[70]  S. Narum,et al.  Comparison of SNPs and microsatellites for fine‐scale application of genetic stock identification of Chinook salmon in the Columbia River Basin , 2011, Molecular ecology resources.

[71]  Damon P. Little,et al.  Choosing and Using a Plant DNA Barcode , 2011, PloS one.

[72]  M. Fortin,et al.  Measuring ecological niche overlap from occurrence and spatial environmental data , 2012 .

[73]  J. Doucet,et al.  Speciation slowing down in widespread and long-living tree taxa: insights from the tropical timber tree genus Milicia (Moraceae) , 2014, Heredity.

[74]  Aurélien Miralles,et al.  The integrative future of taxonomy , 2010, Frontiers in Zoology.

[75]  D. Janzen,et al.  Use of DNA barcodes to identify flowering plants. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[76]  Ryan J. Haasl,et al.  Multi-locus inference of population structure: a comparison between single nucleotide polymorphisms and microsatellites , 2011, Heredity.

[77]  D. Goldstein,et al.  Launching microsatelites: a review of mutation processes and methods of phylogenetic inference , 1997 .

[78]  J. Doucet,et al.  Forest refugia revisited: nSSRs and cpDNA sequences support historical isolation in a wide‐spread African tree with high colonization capacity, Milicia excelsa (Moraceae) , 2010, Molecular ecology.

[79]  Ragupathy Subramanyam,et al.  Testing plant barcoding in a sister species complex of pantropical Acacia (Mimosoideae, Fabaceae) , 2009, Molecular ecology resources.