Integrative taxonomy of the primitively segmented spider genus Ganthela (Araneae: Mesothelae: Liphistiidae): DNA barcoding gap agrees with morphology

Species delimitation is difficult for taxa in which the morphological characters are poorly known because of the rarity of adult morphs or sexes, and in cryptic species. In primitively segmented spiders, family Liphistiidae, males are often unknown, and female genital morphology – usually species‐specific in spiders – exhibits considerable intraspecific variation. Here, we report on an integrative taxonomic study of the liphistiid genus Ganthela Xu & Kuntner, 2015, endemic to south‐east China, where males are only available for two of the seven morphological species (two known and five undescribed). We obtained DNA barcodes (cytochrome c oxidase subunit I gene, COI) for 51 newly collected specimens of six morphological species and analysed them using five species‐delimitation methods: DNA barcoding gap, species delimitation plugin [P ID(Liberal)], automatic barcode gap discovery (ABGD), generalized mixed Yule‐coalescent model (GMYC), and statistical parsimony (SP). Whereas the first three agreed with the morphology, GMYC and SP indicate several additional species. We used the consensus results to delimit and diagnose six Ganthela species, which in addition to the type species Ganthela yundingensis Xu, 2015, completes the revision of the genus. Although multi‐locus phylogenetic approaches may be needed for complex taxonomic delimitations, our results indicate that even single‐locus analyses based on the COI barcodes, if integrated with morphological and geographical data, may provide sufficiently reliable species delimitation. © 2015 The Linnean Society of London     

Species delimitation with gene flow: A methodological comparison and population genomics approach to elucidate cryptic species boundaries in Malaysian Torrent Frogs

Accurately delimiting species boundaries is a nontrivial undertaking that can have significant effects on downstream inferences. We compared the efficacy of commonly used species delimitation methods (SDMs) and a population genomics approach based on genomewide single‐nucleotide polymorphisms (SNPs) to assess lineage separation in the Malaysian Torrent Frog Complex currently recognized as a single species (Amolops larutensis). First, we used morphological, mitochondrial DNA and genomewide SNPs to identify putative species boundaries by implementing noncoalescent and coalescent‐based SDMs (mPTP, iBPP, BFD*). We then tested the validity of putative boundaries by estimating spatiotemporal gene flow (fastsimcoal2 , ABBA‐BABA) to assess the extent of genetic isolation among putative species. Our results show that the A. larutensis complex runs the gamut of the speciation continuum from highly divergent, genetically isolated lineages (mean F_st = 0.9) to differentiating populations involving recent gene flow (mean F_st = 0.05; N_m > 5). As expected, SDMs were effective at delimiting divergent lineages in the absence of gene flow but overestimated species in the presence of marked population structure and gene flow. However, using a population genomics approach and the concept of species as separately evolving metapopulation lineages as the only necessary property of a species, we were able to objectively elucidate cryptic species boundaries in the presence of past and present gene flow. This study does not discount the utility of SDMs but highlights the danger of violating model assumptions and the importance of carefully considering methods that appropriately fit the diversification history of a particular system.     

Estimating the temporal and spatial extent of gene flow among sympatric lizard populations (genus Sceloporus) in the southern Mexican highlands

Interspecific gene flow is pervasive throughout the tree of life. Although detecting gene flow between populations has been facilitated by new analytical approaches, determining the timing and geography of hybridization has remained difficult, particularly for historical gene flow. A geographically explicit phylogenetic approach is needed to determine the overlap of ancestral populations. In this study, we performed population genetic analyses, species delimitation, simulations and a recently developed approach of species tree diffusion to infer the phylogeographic history, timing and geographic extent of gene flow in lizards of the Sceloporus spinosus group. The two species in this group, S. spinosus and S. horridus, are distributed in eastern and western portions of Mexico, respectively, but populations of these species are sympatric in the southern Mexican highlands. We generated data consisting of three mitochondrial genes and eight nuclear loci for 148 and 68 individuals, respectively. We delimited six lineages in this group, but found strong evidence of mito‐nuclear discordance in sympatric populations of S. spinosus and S. horridus owing to mitochondrial introgression. We used coalescent simulations to differentiate ancestral gene flow from secondary contact, but found mixed support for these two models. Bayesian phylogeography indicated more than 60% range overlap between ancestral S. spinosus and S. horridus populations since the time of their divergence. Isolation–migration analyses, however, revealed near‐zero levels of gene flow between these ancestral populations. Interpreting results from both simulations and empirical data indicate that despite a long history of sympatry among these two species, gene flow in this group has only recently occurred.     

Caribbean golden orbweaving spiders maintain gene flow with North America

The Caribbean archipelago offers one of the best natural arenas for testing biogeographic hypotheses. The intermediate dispersal model of biogeography (IDM) predicts variation in species richness among lineages on islands to relate to their dispersal potential. To test this model, one would need background knowledge of dispersal potential of lineages and their biogeographic patterns, which has been problematic as evidenced by our prior work on the Caribbean tetragnathid spiders. In order to investigate the biogeographic imprint of an excellent disperser, we study Trichonephila in the Americas. Trichonephila is a nephilid genus that contains globally distributed species known to overcome long, overwater distances. The results of our phylogenetic and population genetic analyses on T. clavipes suggest that populations over the Caribbean and North America maintain a lively gene flow. However, the single species status of T. clavipes over the entire New World is challenged by our species delimitation analyses. Combined with prior evidence from spider genera of different dispersal ability, these patterns coming from an excellent disperser (Trichonephila) that is species-poor and of a relatively homogenous genetic structure, support the IDM predictions.     

Digest: Species delimitation in the face of demographic processes*

Existing approaches to species delimitation use the extent of divergence between taxa. However, processes, such as gene flow during divergence or secondary contact, as well as population expansion and migration, complicate this task. Smith and Carstens introduce the R package delimitR, which uses machine learning to integrate gene flow into species delimitation inference.     

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     

GENE FLOW AND SPECIES DELIMITATION: A CASE STUDY OF TWO PINE SPECIES WITH OVERLAPPING DISTRIBUTIONS IN SOUTHEAST CHINA

Species delimitation detected by molecular markers is complicated by introgression and incomplete lineage sorting between species. Recent modeling suggests that fixed genetic differences between species are highly related to rates of intraspecific gene flow. However, it remains unclear whether such differences are due to high levels of intraspecific gene flow overriding the spread of introgressed alleles or favoring rapid lineage sorting between species. In pines, chloroplast (cp) and mitochondrial (mt) DNAs are normally paternally and maternally inherited, respectively, and thus their relative rates of intraspecific gene flow are expected to be high and low, respectively. In this study, we used two pine species with overlapping geographical distributions in southeast China, P. massoniana and P. hwangshanensis, as a model system to examine the association between organelle gene flow and variation within and between species. We found that cpDNA variation across these two pine species is more species specific than mtDNA variation and almost delimits taxonomic boundaries. The shared mt/cp DNA genetic variation between species shows no bias in regard to parapatric versus allopatric species’ distributions. Our results therefore support the hypothesis that high intraspecific gene flow has accelerated cpDNA lineage sorting between these two pine species.     

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.     

When taxonomy meets genomics: lessons from a common songbird

Taxonomy is being increasingly informed by genomics. Traditionally, taxonomy has relied extensively on phenotypic traits for the identification and delimitation of species, though with a growing influence from molecular phylogenetics in recent decades. Now, genomics opens up new and more powerful tools for analysing the evolutionary history and relatedness among species, as well as understanding the genetic basis for phenotypic traits and their role in reproductive isolation. New insights gained from genomics will therefore have major effects on taxonomic classifications and species delimitation. How a genomics approach can inform a flawed taxonomy is nicely exemplified by Mason & Taylor ( 2015 ) in this issue of Molecular Ecology. They studied redpolls, which comprise a genus (Acanthis) of fringillid finches with a wide distribution in the Holarctic region, and whose species taxonomy has been a matter of much controversy for decades (Fig.  1 ). Current authoritative checklists classify them into one, two or three species, and five or six subspecies, based largely on geographical differences in phenotypic traits. Previous studies, including a recent one of the subspecies on Iceland (Amouret et al. 2015 ), have found no evidence of differentiation between these taxa in conventional molecular markers. The lack of genetic structure has been interpreted as incomplete lineage sorting among rapidly evolving lineages. Now Mason & Taylor ( 2015 ), using a large data set of genomewide SNPs, verify that they all belong to a single gene pool with a common evolutionary history, and with little or no geographical structuring. They also show that phenotypic traits used in taxonomic classifications (plumage and bill morphology) are closely associated with polygenic patterns of gene expression, presumably driven by ecological selection on a few regulatory genes. Several lessons can be learned from this study. Perhaps the most important one for taxonomy is the risk of taxonomic inflation resulting from overemphasizing phenotypic traits under local adaptation and ignoring a lack of phylogenetic signal in molecular markers.     

Non-equilibrium estimates of gene flow inferred from nuclear genealogies suggest that Iberian and North African wall lizards (<Emphasis Type="Italic">Podarcis</Emphasis> spp.) are an assemblage of incipient species

Background  

The study of recently-diverged species offers significant challenges both in the definition of evolutionary entities and in the estimation of gene flow among them. Iberian and North African wall lizards (Podarcis) constitute a cryptic species complex for which previous assessments of mitochondrial DNA (mtDNA) and allozyme variation are concordant in describing the existence of several highly differentiated evolutionary units. However, these studies report important differences suggesting the occurrence of gene flow among forms. Here we study sequence variation in two nuclear introns, β-fibint7 and 6-Pgdint7, to further investigate overall evolutionary dynamics and test hypotheses related to species delimitation within this complex.     

An assessment of species limits of the South American mouse genus Oligoryzomys (Rodentia,Cricetidae) using unilocus delimitation methods

Oligoryzomys, as currently understood is formed by 25 living species, is the most diverse genus of the tribe Oryzomyini of the New World subfamily Sigmodontinae of cricetid rodents. Nonetheless, the species richness of Oligoryzomys seems to be an underestimate, given some species complex has been proposed in previous studies, at the time that large geographic areas remain to be sampled, and several taxonomic forms have not been assessed with contemporary approaches. In this study, we present a new assessment of the species diversity of Oligoryzomys based on multiple unilocus species delimitation methods (ABGD, BPP, PTP, GMYC and b GMYC), using 665 cytb gene sequences as evidence (532 gathered from Genbank and 133 obtained in this study). We sampled representatives of almost all currently known species of Oligoryzomys, at the time that extending the geographic coverage to the Central Andes, a large area that was largely unrepresented in previous studies. Phylogenetic relationships, based on a non‐redundant alignment, were inferred via maximum likelihood and Bayesian inference; an ultrametric tree, used in species delimitation analyses, was obtained using multiple secondary calibration points. Results of species delimitation methods are discussed at the light of previous knowledge (e.g., taxonomic history and geographic provenance of samples in relation to type localities) and the morphological assessments of some specimens. Results of the distinct delimitation methods are mostly congruent, being BPP and PTP the most sensible to estimate species delimitation, allowing us to suggest that Oligoryzomys is composed of 30 lineages of species level. Of these, 22 correspond to forms currently considered species; some of these include in their synonymy some forms currently considered valid species (e.g., yatesi would be a synonym of longicaudatus). The remaining eight lineages are candidate species that need to be further evaluated. This study, by advancing taxonomic hypothesis that should be further tested in future studies, constitutes a stepping‐stone for upcoming taxonomic and biogeographic studies centred on Oligoryzomys.     

Coalescence-based species delimitation using genome-wide data reveals hidden diversity in a cosmopolitan group of lichens

Although there is an abundance of species delimitation methods on the market, most approaches depend on predefined assignment of specimens to species or populations. Assignment-free methods, which can simultaneously infer boundaries and relationships among species, are of high importance in cases, when correct pre-assignment is difficult or not at all possible. In this study, we use assignment-free multispecies coalescent-based species delimitation (STACEY, tr2-delimitation, and BP&P), phylogenetic methods, and clustering algorithms to investigate the inter- and infraspecific relationships within a common and widespread group of lichens with contentious species boundaries. The Cetraria aculeata group presents a good example of extreme morphological variability and unclear species delimitation in lichens. Based on DNA-sequence data from 26 fungal loci and 10 microsatellite loci, as well as morphological and chemical data, our results provide evidence for the occurrence of five different taxa within the group and highlight the difficulties of morphologically distinguishing these species. We discovered a separate lineage (clade C) within C. aculeata s. str., which does not fully coincide with any of the a priori identified species C. aculeata, C. crespoae, or C. steppae and conclude that this clade constitutes a semi-cryptic, genetically isolated lineage within C. aculeata. We recognize this lineage at subspecific rank as C. aculeata subsp. steppae and synonymize Cetraria crespoae with C. aculeata subsp. aculeata. Epitypes are designated for all involved names to stabilize their usage. The PKS8 gene locus is recommended as a barcode for the separation of C. aculeata subsp. aculeata and subsp. steppae. We demonstrate the potential use of microsatellite data for species delimitation in lichens that might offer an alternative insight or be used to test species delimitation hypotheses, when dealing with closely related or potentially cryptic species. Our results also confirm the presence of an undescribed sister lineage to C. odontella previously misidentified as C. muricata and extend the known range of this lineage to Central Asia (Altay Mts.) and the Central European Alps (France, Switzerland), which calls for a critical reappraisal of records of C. aculeata and C. muricata from these mountain ranges.

     

Integrative analyses of Nervilia (Orchidaceae) section Linervia reveal further undescribed cryptic diversity in Thailand

The delimitation of cryptic species is necessary to accurately classify and appropriately conserve biodiversity. Integrative analyses can be incisive in detecting and circumscribing cryptic diversity, especially in species complexes whose members are delineated by minor or overlapping morphological variation. We adopt an integrative approach to assess species relationships and resolve species boundaries in the taxonomically difficult Nervilia adolphi/punctata species alliance of N. sect. Linervia, an Old World complex of reduced, one-flowered terrestrial orchids that is both species-rich and poorly known in tropical and warm temperate Asia. We sampled 12 of the 27 known species of the alliance in Asia, including all four species reported from Thailand and a further 20 plants collected in that country that could not be satisfactorily identified using morphology alone. Phylogenetic analyses using one nuclear (ITS) and two plastid (matK and trnL-F) markers confirmed both N. sect. Linervia and the alliance itself as monophyletic, and corroborated 11 of the 12 sampled species; N. punctata proved polyphyletic, with the Thai samples referred to this Indonesian species falling sister to the Himalayan N. mackinnonii. The 20 unidentified Thai samples formed three distinct, strongly supported clades. STACEY, a Bayesian coalescence approach to species delimitation, resolved the same three clusters, but provided evidence suggesting that one comprised two distinct sub-clades. Building on this genetic evidence, we identify subtle morphological differences and invoke a diagnosable species concept to circumscribe three previously unrecognized cryptic species from Thailand. This objective approach to species delimitation validates ostensibly minor morphological differences as a basis for differentiating species within the alliance, paving the way for a global analysis of species boundaries throughout the genus as a whole.     

Take up the challenge! Opportunities for evolution research from resolving conflict in integrative taxonomy

What's in a species? The multiple connotations of the question tend to lack simple answers, and not surprisingly so. For example, speciation is a gradual process. Can we say when exactly a child has become an adult? We have precocious youngsters and late bloomers, and often, adults are in some ways childish. There are many triggers for and routes to adolescence. All this holds for speciation, and delimiting species can therefore be a tricky task. Recently, the field of integrative taxonomy has emerged—species delimitation based on multiple sources of evidence. Given that we expect species to exhibit peculiarities in at least one or a few aspects, might it be their alleles of a gene, their morphology, chemistry, behaviour, ecology, reproductive compatibility, or whatever, investigating not just one but several of these aspects makes it more likely that we capture such peculiarities. If the same pattern is found multiply, we talk about agreement among disciplines, and species delimitation is easy. But what if different disciplines tell different stories? Such disagreement makes species delimitation more difficult but is also an opportunity for evolutionary biology (Schlick‐Steiner et al. 2010). In this issue of Molecular Ecology, Andújar et al. (2014) present a comprehensive integrative‐taxonomic case study of Mesocarabus ground beetles including nomenclatural consequences. They resolve extensive disagreement among disciplines by invoking evolutionary explanations, and the process of conflict resolution thus advances knowledge on species boundaries and evolutionary processes simultaneously.     

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.     

Species delimitation methods reveal cryptic diversity in the Hypnea cornuta complex (Cystocloniaceae,Rhodophyta)

The simple and convergent morphologies of many red algae make these species difficult to identify using traditional morphological characters. Many cryptic species have been described in recent years based on molecular datasets, and this has led to the application of an integrative taxonomy approach in species delimitation. Here, we performed several species delimitation methods (mBGD, ABGD, SPN, PTP, GMYCs and GMYCm) based on two different loci (COI-5P and rbcL) in species of the Hypnea cornuta complex. These methods were combined with morphological and phylogenetic data, extensive sampling, analysis of topotype material, and historically relevant herbarium samples. Our findings demonstrate that the groups morphologically assigned to H. cornuta and H. stellulifera consist of five different cryptic species. H. cornuta is a polyphyletic taxon composed of three well-separated lineages, thus requiring sequencing of type or topotype specimens to determine which one is Hypnea cornuta sensu stricto. We have revealed that the distribution of H. stellulifera is limited to Asia, while the Brazilian specimens initially assigned to this species were clarified as a new endemic species: Hypnea cryptica sp. nov. Our results indicated that only an integrative approach combining several lines of evidence, including morphology, nomenclature history, molecular data, biogeography and ecology can correctly solve the taxonomic status of widely distributed cryptic species.     

Integrative taxonomy uncovers hidden species diversity in the rheophilic genus Potamometra (Hemiptera: Gerridae)

The genus Potamometra Bianchi, 1896 represents big rheophilic semi-aquatic bugs that typically inhabit middle-altitude mountainous streams. Here, we integrated molecular and morphological data for delimiting species boundaries and understanding the evolutionary history of the genus Potamometra. Twenty-seven complete mitochondrial genomes of Potamometra were sequenced, with samples representing most of the known geographically distributed locations around the Sichuan Basin. The results of different species delimitation methods (ABGD, bPTP, GMYC and BPP) based on the monolocus or multilocus data strongly supported the existence of two cryptic new species (Potamometra anderseni Zheng, Ye & Bu, sp. nov. and Potamometra zhengi Zheng, Ye & Bu, sp. nov.) although more entities were found in the tree-based delimitation methods. The two new species were successfully validated using morphological characters within a detailed taxonomic framework. Phylogenetic analyses supported the reciprocal monophyly of the seven highly node-supported clades, which were matched with the five known species and two new taxa. A novel gene arrangement pattern that two trnF (trnF1 and trnF2) genes separated by an intergenic spacer (IGS) were found in all the species except the sister group of Potamometra berezowskii Bianchi, 1896 and Potamometra linnavuorii Chen, Nieser & Bu, 2016. This gene rearrangement event could be explained by the tandem duplication and random loss (TDRL) model. Our study emphasized that the combination of molecular sequence data, morphological characters and mitochondrial structural information could improve the accuracy of species delimitation.     

The influence of locus number and information content on species delimitation: an empirical test case in an endangered Mexican salamander

Perhaps the most important recent advance in species delimitation has been the development of model‐based approaches to objectively diagnose species diversity from genetic data. Additionally, the growing accessibility of next‐generation sequence data sets provides powerful insights into genome‐wide patterns of divergence during speciation. However, applying complex models to large data sets is time‐consuming and computationally costly, requiring careful consideration of the influence of both individual and population sampling, as well as the number and informativeness of loci on species delimitation conclusions. Here, we investigated how locus number and information content affect species delimitation results for an endangered Mexican salamander species, Ambystoma ordinarium. We compared results for an eight‐locus, 137‐individual data set and an 89‐locus, seven‐individual data set. For both data sets, we used species discovery methods to define delimitation models and species validation methods to rigorously test these hypotheses. We also used integrated demographic model selection tools to choose among delimitation models, while accounting for gene flow. Our results indicate that while cryptic lineages may be delimited with relatively few loci, sampling larger numbers of loci may be required to ensure that enough informative loci are available to accurately identify and validate shallow‐scale divergences. These analyses highlight the importance of striking a balance between dense sampling of loci and individuals, particularly in shallowly diverged lineages. They also suggest the presence of a currently unrecognized, endangered species in the western part of A. ordinarium's range.