Considering gene flow when using coalescent methods to delimit lineages of North American pitvipers of the genus Agkistrodon

Examining species diversity and mechanisms of speciation using coalescent models provides a framework for how regional diversity is accrued, even in well‐studied areas such as the Nearctic. It is likely, that gene flow among closely‐related species with adjacent distributions may be common. However, the absence of gene flow is a primary assumption of many phylogeographical methods that produce species trees and delimit species using Bayesian or likelihood functions in a coalescent framework. In the present study, we examine delimitation when gene flow between species is present using empirical datasets from two species of North American pitvipers of the genus Agkistrodon. We also use niche modelling to determine whether these young lineages occur in distinct environmental niches. To manage the problem of gene flow between species, we first identify admixed individuals, demonstrate that gene flow has occurred, and then identify the impact of alternative population assignments of admixed individuals on delimitation posterior probabilities. In addition, we examine the influence of mitochondrial genes relative to other loci combined in coalescent analyses that delimit species. Here, we find that the copperheads (Agkistrodon contortrix) and the cottonmouths (Agkistrodon piscivorus) are each composed of two distinct species, with each occupying different niches. Importantly, we find that species can be delimited when the amount of gene flow between lineages is low, although the methods are acutely sensitive to population assignment of individuals. © 2014 The Linnean Society of London     

Species delimitation with abc and other coalescent-based methods: a test of accuracy with simulations and an empirical example with lizards of the liolaemus darwinii complex (squamata: liolaemidae)

Species delimitation is a major research focus in evolutionary biology because accurate species boundaries are a prerequisite for the study of speciation. New species delimitation methods (SDMs) can accommodate nonmonophyletic species and gene tree discordance as a result of incomplete lineage sorting via the coalescent model, but do not explicitly accommodate gene flow after divergence. Approximate Bayesian computation (ABC) can incorporate gene flow and estimate other relevant parameters of the speciation process while testing alternative species delimitation hypotheses. We evaluated the accuracy of BPP, SpeDeSTEM, and ABC for delimiting species using simulated data and applied these methods to empirical data from lizards of the Liolaemus darwinii complex. Overall, BPP was the most accurate, ABC showed an intermediate accuracy, and SpeDeSTEM was the least accurate under most simulated conditions. All three SDMs showed lower accuracy when speciation occurred despite gene flow, as found in previous studies, but ABC was the method with the smallest decrease in accuracy. All three SDMs consistently supported the distinctness of southern and northern lineages within L. darwinii. These SDMs based on genetic data should be complemented with novel SDMs based on morphological and ecological data to achieve truly integrative and statistically robust approaches to species discovery.     

Bayesian Poisson tree processes and multispecies coalescent models shed new light on the diversification of Nawab butterflies in the Solomon Islands (Nymphalidae,Charaxinae, Polyura)

Butterflies of the genus Polyura form a widespread tropical group distributed from Pakistan to Fiji. The rare endemic Polyura epigenes Godman & Salvin, 1888 from the Solomon Islands archipelago represents a case of marked island polymorphism. We sequenced museum specimens of this species across its geographic range to study the phylogeography and genetic differentiation of populations in the archipelago. We used the Bayesian Poisson tree processes and multispecies coalescent models, to study species boundaries. We also estimated divergence times to investigate the biogeographic history of populations. Our molecular species delimitation and nuclear DNA network analyses unambiguously indicate that Malaita populations form an independent metapopulation lineage, as defined in the generalized lineage concept. This lineage, previously ranked as a subspecies, is raised to species rank under the name Polyura bicolor Turlin & Sato, 1995  stat. nov. Divergence time estimates suggest that this lineage split from its sister taxon in the late Pleistocene. At this time, the bathymetric isolation of Malaita from the rest of the archipelago probably prevented gene flow during periods of lower sea level, thereby fostering allopatric speciation. The combination of molecular species delimitation methods, morphological comparisons, and divergence time estimation is useful to study lineage diversification across intricate geographic regions.     

Miocene and Pliocene dominated diversification of the lichen-forming fungal genus Melanohalea (Parmeliaceae, Ascomycota) and Pleistocene population expansions

ABSTRACT: BACKGROUND: Factors promoting diversification in lichen symbioses remain largely unexplored. While Pleistocene events have been important for driving diversification and affecting distributions in many groups, recent estimates suggest that major radiations within some genera in the largest clade of macrolichens (Parmeliaceae, Ascomycota) vastly predate the Pleistocene. To better understand the temporal placement and sequence of diversification events in lichens, we estimated divergence times in a common lichen-forming fungal genus, Melanohalea, in the Northern Hemisphere. Divergence times were estimated using both concatenated gene tree and coalescence-based multilocus species tree approaches to assess the temporal context of major radiation events within Melanohalea. In order to complement our understanding of processes impacting genetic differentiation, we also evaluated the effects of Pleistocene glacial cycles on population demographics of distinct Melanohalea lineages, differing in reproductive strategies. RESULTS: We found that divergence estimates, from both concatenated gene tree and coalescence-based multilocus species tree approaches, suggest that diversification within Melanohalea occurred predominantly during the Miocene and Pliocene, although estimated of divergence times differed by up to 8.3 million years between the two methods. These results indicate that, in some cases, taxonomically diagnostic characters may be maintained among divergent lineages for millions of years. In other cases, similar phenotypic characters among non-sister taxa, including reproductive strategies, suggest the potential for convergent evolution due to similar selective pressures among distinct lineages. Our analyses provide evidence of population expansions predating the last glacial maximum in the sampled lineages. These results suggest that Pleistocene glaciations were not inherently unfavorable or restrictive for some Melanohalea species, albeit with apparently different demographic histories between sexually and vegetatively reproducing lineages. CONCLUSIONS: Our results contribute to the understanding of how major changes during the Miocene and Pliocene have been important in promoting diversification within common lichen-forming fungi in the northern Hemisphere. Additionally, we provide evidence that glacial oscillations have influenced current population structure of broadly distributed lichenized fungal species throughout the Holarctic.     

Late Pleistocene climate change promoted divergence between Picea asperata and P. crassifolia on the Qinghai–Tibet Plateau through recent bottlenecks

Divergence during the early stage of speciation can be driven by a population bottleneck via reduced gene flow and enhanced lineage sorting. In this study, we aimed to examine whether such bottlenecks occurred during the initial speciation of two closely related spruce species Picea asperata and P. crassifolia occurring on the Qinghai–Tibet Plateau (QTP). We analyzed sequences of three chloroplast, two mitochondrial DNA fragments and a further 13 nuclear loci from 216 individuals of the two species. Both species showed a low level of genetic diversity in contrast to other congeners occurring in the QTP and adjacent regions. The estimated population sizes of P. asperata and P. crassifolia are less than the ancestral population size before splitting. These results together with multiple statistical tests (Tajima's D, Fu and Li's D* and F*) suggest that these two species underwent recent bottlenecks. Based on approximate Bayesian computation (ABC), we also determined that the period of the population shrinkage was consistent with the interspecific divergence during the late Pleistocene. The reduced population sizes and the divergent selection may together have triggered the initial divergence under high gene flow between these two species. Our results therefore highlight the importance of climatic oscillations during the late Pleistocene in promoting speciation through changing demographic sizes of the ancestral species on the QTP and in adjacent regions.     

Testing species boundaries in an ancient species complex with deep phylogeographic history: genus Xantusia (Squamata: Xantusiidae)

Identification of species in natural populations has recently received increased attention with a number of investigators proposing rigorous methods for species delimitation. Morphologically conservative species (or species complexes) with deep phylogenetic histories (and limited gene flow) are likely to pose particular problems when attempting to delimit species, yet this is crucial to comparative studies of the geography of speciation. We apply two methods of species delimitation to an ancient group of lizards (genus Xantusia) that occur throughout southwestern North America. Mitochondrial cytochrome b and nicotinamide adenine dinucleotide dehydrogenase subunit 4 gene sequences were generated from samples taken throughout the geographic range of Xantusia. Maximum likelihood, Bayesian, and nested cladogram analyses were used to estimate relationships among haplotypes and to infer evolutionary processes. We found multiple well-supported independent lineages within Xantusia, for which there is considerable discordance with the currently recognized taxonomy. High levels of sequence divergence (21.3%) suggest that the pattern in Xantusia may predate the vicariant events usually hypothesized for the fauna of the Baja California peninsula, and the existence of deeply divergent clades (18.8%-26.9%) elsewhere in the complex indicates the occurrence of ancient sundering events whose genetic signatures were not erased by the late Wisconsin vegetation changes. We present a revised taxonomic arrangement for this genus consistent with the distinct mtDNA lineages and discuss the phylogeographic history of this genus as a model system for studies of speciation in North American deserts.     

DNA-based species delimitation in algae

Given the problems of species delimitation in algae using morphology or sexual compatibility, molecular data are becoming the standard for delimiting species and testing their traditional boundaries. The idea that species are separately evolving metapopulation lineages, along with theoretical progress in phylogenetic and population genetic analyses, has led to the development of new methods of species delimitation. We review these recent developments in DNA-based species delimitation methods, and discuss how they have changed and continue to change our understanding of algal species boundaries. Although single-locus approaches have proven effective for a first rapid and large-scale assessment of species diversity, species delimitation based on single gene trees falls short due to gene tree–species tree incongruence, caused by confounding processes like incomplete lineage sorting, trans-species polymorphism, hybridization and introgression. Data from unlinked loci and multi-species coalescent methods, which combine principles from phylogenetics and population genetics, may now be able to account for these complicating factors. Several of these methods also provide statistical support regarding species boundaries, which is important because speciation is a process and therefore uncertainty about precise species boundaries is inevitable in recently diverged lineages.     

Sky island diversification meets the multispecies coalescent – divergence in the spruce‐fir moss spider (Microhexura montivaga,Araneae, Mygalomorphae) on the highest peaks of southern Appalachia

Microhexura montivaga is a miniature tarantula‐like spider endemic to the highest peaks of the southern Appalachian mountains and is known only from six allopatric, highly disjunct montane populations. Because of severe declines in spruce‐fir forest in the late 20th century, M. montivaga was formally listed as a US federally endangered species in 1995. Using DNA sequence data from one mitochondrial and seven nuclear genes, patterns of multigenic genetic divergence were assessed for six montane populations. Independent mitochondrial and nuclear discovery analyses reveal obvious genetic fragmentation both within and among montane populations, with five to seven primary genetic lineages recovered. Multispecies coalescent validation analyses [guide tree and unguided Bayesian Phylogenetics and Phylogeography (BPP), Bayes factor delimitation (BFD)] using nuclear‐only data congruently recover six or seven distinct lineages; BFD analyses using combined nuclear plus mitochondrial data favour seven or eight lineages. In stark contrast to this clear genetic fragmentation, a survey of secondary sexual features for available males indicates morphological conservatism across montane populations. While it is certainly possible that morphologically cryptic speciation has occurred in this taxon, this system may alternatively represent a case where extreme population genetic structuring (but not speciation) leads to an oversplitting of lineage diversity by multispecies coalescent methods. Our results have clear conservation implications for this federally endangered taxon and illustrate a methodological issue expected to become more common as genomic‐scale data sets are gathered for taxa found in naturally fragmented habitats.     

Incorporating the speciation process into species delimitation

The “multispecies” coalescent (MSC) model that underlies many genomic species-delimitation approaches is problematic because it does not distinguish between genetic structure associated with species versus that of populations within species. Consequently, as both the genomic and spatial resolution of data increases, a proliferation of artifactual species results as within-species population lineages, detected due to restrictions in gene flow, are identified as distinct species. The toll of this extends beyond systematic studies, getting magnified across the many disciplines that rely upon an accurate framework of identified species. Here we present the first of a new class of approaches that addresses this issue by incorporating an extended speciation process for species delimitation. We model the formation of population lineages and their subsequent development into independent species as separate processes and provide for a way to incorporate current understanding of the species boundaries in the system through specification of species identities of a subset of population lineages. As a result, species boundaries and within-species lineages boundaries can be discriminated across the entire system, and species identities can be assigned to the remaining lineages of unknown affinities with quantified probabilities. In addition to the identification of species units in nature, the primary goal of species delimitation, the incorporation of a speciation model also allows us insights into the links between population and species-level processes. By explicitly accounting for restrictions in gene flow not only between, but also within, species, we also address the limits of genetic data for delimiting species. Specifically, while genetic data alone is not sufficient for accurate delimitation, when considered in conjunction with other information we are able to not only learn about species boundaries, but also about the tempo of the speciation process itself.     

Evolutionary History of the Morocco lizard-Fingered Geckos of the <Emphasis Type="Italic">Saurodactylus brosseti</Emphasis> Complex

Studies of biodiversity in the Maghreb have revealed high genetic diversity and divergent genetic lineages among many taxa including squamates. Geographic barriers such as the Atlas Mountains are one of the key factors promoting genetic differentiation and the high levels of endemism. The lizard-fingered gecko Saurodactylus brosseti is endemic to Morocco. Its range includes both sides of the Atlas Mountains, and although high diversity was previously identified within the species, much of the range was unsampled. To understand the evolutionary and biogeographical history of this species, we used mitochondrial and nuclear DNA sequence data from 64 populations sampled across the entire species range. We employed phylogenetic methods based on gene trees and species trees as well as a time calibrated Bayesian genealogy and coalescent species delimitation approaches. We uncovered four highly divergent and allopatric mitochondrial lineages that did not share haplotypes at variable nuclear loci, suggesting the four groups have been evolving independently since the Miocene, according to our molecular dating estimates. These results coupled with the geographic pattern of genetic diversity suggest a possible role of the Atlas Mountains for the divergence observed between the four lineages of S. brosseti, while each lineage probably later underwent several episodes of fragmentation followed by (re-) expansion during Pleistocene climatic oscillations. Bayesian species delimitation analysis indicates that the four lineages may well be distinct species but we suggest that detailed morphological analyses are needed prior to taxonomic changes. The four lineages represent ancient independent evolutionary units, and deserve conservation management as distinct entities.     

Multilocus phylogeographic assessment of the California Mountain Kingsnake (Lampropeltis zonata) suggests alternative patterns of diversification for the California Floristic Province

Phylogeographic inference can determine the timing of population divergence, historical demographic processes, patterns of migration, and when extended to multiple species, the history of communities. Single‐locus analyses can mislead interpretations of the evolutionary history of taxa and comparative analyses. It is therefore important to revisit previous single‐locus phylogeographic studies, particularly those that have been used to propose general patterns for regional biotas and the processes responsible for generating inferred patterns. Here, we employ a multilocus statistical approach to re‐examine the phylogeography of Lampropeltis zonata. Using nonparametic and Bayesian species delimitation, we determined that there are two well‐supported species within L. zonata. Ecological niche modelling supports the delimitation of these taxa, suggesting that the two species inhabit distinct climatic environments. Gene flow between the two taxa is low and appears to occur unidirectionally. Further, our data suggest that gene flow was mediated by females, a rare pattern in snakes. In contrast to previous analyses, we determined that the divergence between the two lineages occurred in the late Pliocene (c. 2.07 Ma). Spatially and temporally, the divergence of these lineages is associated with the inundation of central California by the Monterey Bay. The effective population sizes of the two species appear to have been unaffected by Pleistocene glaciation. Our increased sampling of loci for L. zonata, combined with previously published multilocus analyses of other sympatric species, suggests that previous conclusions reached by comparative phylogeographic studies conducted within the California Floristic Province should be reassessed.     

Specialist and generalist symbionts show counterintuitive levels of genetic diversity and discordant demographic histories along the Florida Reef Tract

Specialist and generalist life histories are expected to result in contrasting levels of genetic diversity at the population level, and symbioses are expected to lead to patterns that reflect a shared biogeographic history and co-diversification. We test these assumptions using mtDNA sequencing and a comparative phylogeographic approach for six co-occurring crustacean species that are symbiotic with sea anemones on western Atlantic coral reefs, yet vary in their host specificities: four are host specialists and two are host generalists. We first conducted species discovery analyses to delimit cryptic lineages, followed by classic population genetic diversity analyses for each delimited taxon, and then reconstructed the demographic history for each taxon using traditional summary statistics, Bayesian skyline plots, and approximate Bayesian computation to test for signatures of recent and concerted population expansion. The genetic diversity values recovered here contravene the expectations of the specialist–generalist variation hypothesis and classic population genetics theory; all specialist lineages had greater genetic diversity than generalists. Demography suggests recent population expansions in all taxa, although Bayesian skyline plots and approximate Bayesian computation suggest the timing and magnitude of these events were idiosyncratic. These results do not meet the a priori expectation of concordance among symbiotic taxa and suggest that intrinsic aspects of species biology may contribute more to phylogeographic history than extrinsic forces that shape whole communities. The recovery of two cryptic specialist lineages adds an additional layer of biodiversity to this symbiosis and contributes to an emerging pattern of cryptic speciation in the specialist taxa. Our results underscore the differences in the evolutionary processes acting on marine systems from the terrestrial processes that often drive theory. Finally, we continue to highlight the Florida Reef Tract as an important biodiversity hotspot.

     

Mitochondrial DNA phylogeography and mating compatibility reveal marked genetic structuring and speciation in the NE Atlantic bryozoan Celleporella hyalina

The marine bryozoan Celleporella hyalina is a species complex composed of many highly divergent and mostly allopatric genetic lineages that are reproductively isolated but share a remarkably similar morphology. One such lineage commonly encrusts macroalgae throughout the NE Atlantic coast. To explore the processes leading to geographical diversification, reproductive isolation and speciation in this taxon, we (i) investigated NE Atlantic C. hyalina mitochondrial DNA phylogeography, and (ii) used breeding trials between geographical isolates to ascertain reproductive isolation. We find that haplotype diversity is geographically variable and there is a strong population structure, with significant isolation by distance. NE Atlantic C. hyalina is structured into two main parapatric lineages that appear to have had independent Pleistocene histories. Range expansions have resulted in two contact zones in Spain and W Ireland. Lineage 1 is found from Ireland to Spain and has low haplotype diversity, with closely related haplotypes, suggesting a recent population expansion into the Irish Sea, S Ireland, S England and Spain. Lineage 2 is found from Iceland to Spain and has high haplotype diversity. Complete reproductive isolation was found between some geographical isolates representing both lineages, whereas it was incomplete or asymmetric between others, suggesting these latter phylogeographical groups probably represent incipient species. The phylogeographical distribution of NE Atlantic C. hyalina does not fall easily into a pattern of southern refugia, and we discuss likely differences between terrestrial and marine system responses to Pleistocene glacial cycles.     

The species concept as an emergent property of population biology

Resurgent interest in the genetics of population divergence and speciation coincides with recent critical evaluation of species concepts and proposals for species delimitation. An important result of these parallel trends is a slight but important conceptual shift in focus away from species diagnoses based on prior species concepts or definitions, and toward analyses of the processes acting on lineages of metapopulations that eventually lead to differences recognizable as species taxa. An advantage of this approach is that it identifies quantitative metapopulation differences in continuous variables, rather than discrete entities that do or do not conform to a prior species concept, and species taxa are recognized as an emergent property of population-level processes. The tension between species concepts and diagnosis versus emergent recognition of species taxa is at least as old as Darwin, and is unlikely to be resolved soon in favor of either view, because the products of both approaches (discrete utilitarian taxon names for species, process-based understanding of the origins of differentiated metapopulations) continue to have important applications.     

Divergence,gene flow,and speciation in eight lineages of trans‐Beringian birds

Determining how genetic diversity is structured between populations that span the divergence continuum from populations to biological species is key to understanding the generation and maintenance of biodiversity. We investigated genetic divergence and gene flow in eight lineages of birds with a trans‐Beringian distribution, where Asian and North American populations have likely been split and reunited through multiple Pleistocene glacial cycles. Our study transects the speciation process, including eight pairwise comparisons in three orders (ducks, shorebirds and passerines) at population, subspecies and species levels. Using ultraconserved elements (UCEs), we found that these lineages represent conditions from slightly differentiated populations to full biological species. Although allopatric speciation is considered the predominant mode of divergence in birds, all of our best divergence models included gene flow, supporting speciation with gene flow as the predominant mode in Beringia. In our eight lineages, three were best described by a split‐migration model (divergence with gene flow), three best fit a secondary contact scenario (isolation followed by gene flow), and two showed support for both models. The lineages were not evenly distributed across a divergence space defined by gene flow (M) and differentiation (F_ST), instead forming two discontinuous groups: one with relatively shallow divergence, no fixed single nucleotide polymorphisms (SNPs), and high rates of gene flow between populations; and the second with relatively deeply divergent lineages, multiple fixed SNPs, and low gene flow. Our results highlight the important role that gene flow plays in avian divergence in Beringia.     

Methods for delimiting species via population genetics and phylogenetics using genotype data

In species delimitation, a formidable goal in the discipline of systematic biology, we identify and describe species morphologically and ecologically based on phenotypic data. Efficient genotyping technologies produce genetic and genomic data with relative ease, which promotes species discovery and validation using genotype data. For the last two decades, we have seen the development of species delimitation methods based on genetic distances and phylogenetic trees using genotype data. However, speciation processes via evolutionary relationship among species were mostly divorced from species delimitation. Recent approaches to drawing species boundaries use multi-locus sequence data to account for evolutionary processes including speciation and gene flow. They allow us to learn of jointly speciation and species delimitation, leveraging computational and statistical techniques developed in population genetics and phylogenetics. Here, we review the recent progress in the development of species delimitation using genotype data and discuss the future outlook for the research of developing species delimitation methods.     

Dealing with incomplete taxon sampling and diversification of a large clade of mushroom-forming fungi

The absence of an adequate fossil record can hinder understanding the process of diversification that underlies the evolutionary history of a given group. In such cases, investigators have used ultrametric trees derived from molecular data from extant taxa to gain insights into processes of speciation and extinction over time. Inadequate taxon sampling, however, impairs such inferences. In this study, we use simulations to investigate the effect of incomplete taxon sampling on the accumulation of lineages through time for a clade of mushroom-forming fungi, the Hebelomateae. To achieve complete taxon sampling, we use a new Bayesian approach that incorporates substitute lineages to estimate diversification rates. Unlike many studies of animals and plants, we find no evidence of a slowdown in speciation. This indicates the Hebelomateae has not undergone an adaptive radiation. Rather, these fungi have evolved under a relatively constant rate of diversification since their most recent common ancestor, which we date back to the Eocene. The estimated net diversification rate (0.08-0.19 spp./lineage/Ma) is comparable with that of many plants and animals. We suggest that continuous diversification in the Hebelomateae has been facilitated by climatic and vegetation changes throughout the Cenozoic. We also caution against modeling multiple genes as a single partition when performing phylogenetic dating analyses.     

Evolutionary history of bulldog bats (genus Noctilio): recent diversification and the role of the Caribbean in Neotropical biogeography

We present phylogenetic relationships and phylogeographic patterns of the two species of bulldog bats, genus Noctilio. Using a comprehensive sampling of 118 individuals throughout the species distribution, we investigated the distribution of molecular variation in one nuclear and two mitochondrial markers. Phylogenetic trees do not recover Noctilio albiventris as a monophyletic group and point to three similar‐age intraspecific genetic lineages, suggesting cryptic diversity in this taxon. These lineages correspond to the subspecies previously proposed, and are strongly associated with major river basins in South America. Analyses also suggest a very recent origin for the fishing bat Noctilio leporinus, which probably originates from N. albiventris, with a population expansion corresponding to its invasion in South America. Based on our analysis, the speciation event was dated in the Pleistocene epoch and seems to be associated with the variation of the sea level in the Caribbean islands. The present work indicates how phylogeographic studies support the identification of independent evolutionary lineages, driving new systematic/taxonomic investigations, while at the same time shed light on the role of the Caribbean in shaping Neotropical bat fauna diversity. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ?? , ??–??.     

Testing hypotheses of speciation timing in Dicamptodon copei and Dicamptodon aterrimus (Caudata: Dicamptodontidae)

Giant salamanders of the genus Dicamptodon are members of the mesic forest ecosystem that occurs in the Pacific Northwest of North America. We estimate the phylogeny of the genus to test several hypotheses concerning speciation and the origin of current species distributions. Specifically, we test competing a priori hypotheses of dispersal and vicariance to explain the disjunct inland distribution of the Idaho giant salamander (D. aterrimus) and to test the hypothesis of Pleistocene speciation of Cope's giant salamander (D. copei) using Bayesian hypothesis testing. We determined that available outgroups were too divergent to root the phylogeny effectively, and we calculated Bayesian posterior probabilities for each of the 15 possible root placements for this four-taxon group. This analysis placed the root on the branch leading to D. aterrimus, indicating that current distribution and speciation of D. aterrimus fit the ancient vicariance hypothesis and are attributable to the orogeny of the Cascade Mountains rather than recent inland dispersal. Furthermore, test results indicate that D. copei is distantly related to other coastal lineages and likely originated much earlier than the Pleistocene. These results suggest that speciation within the genus is attributable to ancient geologic events, while more recent Pleistocene glaciation has shaped genetic variation and distributions within the extant species.     

Phylogeography of the heavily poached African common pangolin (Pholidota,Manis tricuspis) reveals six cryptic lineages as traceable signatures of Pleistocene diversification

Knowledge on faunal diversification in African rainforests remains scarce. We used phylogeography to assess (i) the role of Pleistocene climatic oscillations in the diversification of the African common pangolin (Manis tricuspis) and (ii) the utility of our multilocus approach for taxonomic delineation and trade tracing of this heavily poached species. We sequenced 101 individuals for two mitochondrial DNA (mtDNA), two nuclear DNA and one Y‐borne gene fragments (totalizing 2602 bp). We used a time‐calibrated, Bayesian inference phylogenetic framework and conducted character‐based, genetic and phylogenetic delineation of species hypotheses within African common pangolins. We identified six geographic lineages partitioned into western Africa, Ghana, the Dahomey Gap, western central Africa, Gabon and central Africa, all diverging during the Middle to Late Pleistocene. MtDNA (cytochrome b + control region) was the sole locus to provide diagnostic characters for each of the six lineages. Tree‐based Bayesian delimitation methods using single‐ and multilocus approaches gave high support for ‘species’ level recognition of the six African common pangolin lineages. Although the diversification of African common pangolins occurred during Pleistocene cyclical glaciations, causative correlation with traditional rainforest refugia and riverine barriers in Africa was not straightforward. We conclude on the existence of six cryptic lineages within African common pangolins, which might be of major relevance for future conservation strategies. The high discriminative power of the mtDNA markers used in this study should allow an efficient molecular tracing of the regional origin of African common pangolin seizures.