Delimiting Species with Single-Locus DNA Sequences.

DNA sequences are increasingly used for large-scale biodiversity inventories. Because these genetic data avoid the time-consuming initial sorting of specimens based on their phenotypic attributes, they have been recently incorporated into taxonomic workflows for overlooked and diverse taxa. Major statistical developments have accompanied this new practice, and several models have been proposed to delimit species with single-locus DNA sequences. However, proposed approaches to date make different assumptions regarding taxon lineage history, leading to strong discordance whenever comparisons are made among methods. Distance-based methods, such as Automatic Barcode Gap Discovery (ABGD) and Assemble Species by Automatic Partitioning (ASAP), rely on the detection of a barcode gap (i.e., the lack of overlap in the distributions of intraspecific and interspecific genetic distances) and the associated threshold in genetic distances. Network-based methods, as exemplified by the REfined Single Linkage (RESL) algorithm for the generation of Barcode Index Numbers (BINs), use connectivity statistics to hierarchically cluster-related haplotypes into molecular operational taxonomic units (MOTUs) which serve as species proxies. Tree-based methods, including Poisson Tree Processes (PTP) and the General Mixed Yule Coalescent (GMYC), fit statistical models to phylogenetic trees by maximum likelihood or Bayesian frameworks.Multiple webservers and stand-alone versions of these methods are now available, complicating decision-making regarding the most appropriate approach to use for a given taxon of interest. For instance, tree-based methods require an initial phylogenetic reconstruction, and multiple options are now available for this purpose such as RAxML and BEAST. Across all examined species delimitation methods, judicious parameter setting is paramount, as different model parameterizations can lead to differing conclusions. The objective of this chapter is to guide users step-by-step through all the procedures involved for each of these methods, while aggregating all necessary information required to conduct these analyses. The "Materials" section details how to prepare and format input files, including options to align sequences and conduct tree reconstruction with Maximum Likelihood and Bayesian inference. The Methods section presents the procedure and options available to conduct species delimitation analyses, including distance-, network-, and tree-based models. Finally, limits and future developments are discussed in the Notes section. Most importantly, species delimitation methods discussed herein are categorized based on five indicators: reliability, availability, scalability, understandability, and usability, all of which are fundamental properties needed for any approach to gain unanimous adoption within the DNA barcoding community moving forward.

The barcode trap-Description of a new species of Microglanis, with a review of the status of Microglanis cibelae (Siluriformes: Pseudopimelodidae).

In a recent study based on the generalized mixed Yule coalescent method for delimiting species, a threshold of 2% genetic distance using cytochrome c oxidase subunit I sequences was used to delimit the species of Microglanis. That action resulted in assembling several populations of Microglanis from Atlantic coastal rivers between Rio Grande do Sul and São Paulo states as a single species, Microglanis cottoides, including Microglanis cibelae as a junior synonym. We reexamined these populations and found three species diagnosed by their morphology and that constitute separate mtDNA lineages, including a new species. The synonym of M. cibelae and M. cottoides is reviewed and refuted based on morphological and molecular evidence. M. cibelae and the new species are sympatric and occasionally syntopic in the Tramandaí, Mampituba, and Araranguá river basins. The new species is distinguished from M. cibelae and M. cottoides by the anterior margin of the posttemporosupracleitrum narrow articulated with the epioccipital, the short mental and maxillary barbels, and depressed head and body.

SPdel: A pipeline to compare and visualize species delimitation methods for single-locus datasets.

An accurate species delimitation is critical for biological studies. In this context, the use of molecular techniques along with species delimitation methods would help to a rapid and accurate biodiversity assessment. The species delimitation methods cluster data sets of orthologous sequences in molecular operational taxonomic units (MOTU). In particular, the methods based on a single gene are easily integrated with the widely used DNA barcoding approach. We developed SPdel a user-friendly pipeline to integrate different single-gene species delimitation methods. SPdel is designed to calculate and compare MOTUs obtained by different species delimitation approaches. SPdel also outputs diverse ready-to-publish quality figures, that facilitate the interpretation of results. SPdel aims to help researchers use species delimitation methods that would improve biodiversity studies.

Comparison between p-distance and single-locus species delimitation models for delineating reproductively tested strains of pennate diatoms (Bacillariophyceae) using cox1, rbcL and ITS.

Several automated molecular methods have emerged for distinguishing eukaryote species based on DNA sequence data. However, there are knowledge gaps around which of these single-locus methods is more accurate for the identification of microalgal species, such as the highly diverse and ecologically relevant diatoms. We applied genetic divergence, Automatic Barcode Gap Discovery for primary species delimitation (ABGD), Assemble Species by Automatic Partitioning (ASAP), Statistical Parsimony Network Analysis (SPNA), Generalized Mixed Yule Coalescent (GMYC) and Poisson Tree Processes (PTP) using partial cox1, rbcL, 5.8S + ITS2, ITS1 + 5.8S + ITS2 markers to delineate species and compare to published polyphasic identification data (morphological features, phylogeny and sexual reproductive isolation) to test the resolution of these methods. ASAP, ABGD, SPNA and PTP models resolved species of Eunotia, Seminavis, Nitzschia, Sellaphora and Pseudo-nitzschia corresponding to previous polyphasic identification, including reproductive isolation studies. In most cases, these models identified diatom species in similar ways, regardless of sequence fragment length. GMYC model presented smallest number of results that agreed with previous published identification. Following the recommendations for proper use of each model presented in the present study, these models can be useful tools to identify cryptic or closely related species of diatoms, even when the datasets have relatively few sequences.

Comparing the Efficiency of Single-Locus Species Delimitation Methods within Trochoidea (Gastropoda: Vetigastropoda).

In the context of diminishing global biodiversity, the validity and practicality of species delimitation methods for the identification of many neglected and undescribed biodiverse species have been paid increasing attention. DNA sequence-based species delimitation methods are mainly classified into two categories, namely, distance-based and tree-based methods, and have been widely adopted in many studies. In the present study, we performed three distance-based (ad hoc threshold, ABGD, and ASAP) and four tree-based (sGMYC, mGMYC, PTP, and mPTP) analyses based on Trochoidea COI data and analyzed the discordance between them. Moreover, we also observed the performance of these methods at different taxonomic ranks (the genus, subfamily, and family ranks). The results suggested that the distance-based approach is generally superior to the tree-based approach, with the ASAP method being the most efficient. In terms of phylogenetic methods, the single threshold version performed better than the multiple threshold version of GMYC, and PTP showed higher efficiency than mPTP in delimiting species. Additionally, GMYC was found to be significantly influenced by taxonomic rank, showing poorer efficiency in datasets at the genus level than at higher levels. Finally, our results highlighted that cryptic diversity within Trochoidea (Mollusca: Vetigastropoda) might be underestimated, which provides quantitative evidence for excavating the cryptic lineages of these species.

A survey of aquatic macroinvertebrates in a river from the dry corridor of Nicaragua using biological indices and DNA barcoding.

Aquatic macroinvertebrates are widely used as indicators for water quality assessment around the world. Modern strategies for environmental assessment implement molecular analysis to delimitate species of aquatic macroinvertebrates. Delimitation methods have been established to determine boundaries between species units using sequencing data from DNA barcodes and serve as first exploratory tools for taxonomic revisions. This is useful in regions such as the neotropics where aquatic macroinvertebrate habitats are threatened by human interference and DNA databases remain understudied. We asked whether the biodiversity of aquatic macroinvertebrates in a stream in Nicaragua, within the Central American Dry Corridor, could be characterized with biological indices and DNA barcoding. In this study, we combined regional biological indices (BMWP-CR, IBF-SV-2010) along with distance-based (ASAP, BIN) and tree-based (GMYC, bPTP) delimitation methods, as well as nucleotide BLAST in public barcode databases. We collected samples from the upper, middle, and low reaches of the Petaquilla river. The three sites presented excellent water quality with the BMWP-CR index, but evidence of high organic pollution was found in the middle reach with the IBF-SV-2010 index. We report a total of 219 COI sequences successfully generated from 18 families and 8 orders. Operational taxonomic units (OTUs) designation ranged from 69 to 73 using the four methods, with a congruency of 92% for barcode assignation. Nucleotide BLAST identified 14 species (27.4% of barcodes) and 33 genera (39.3% of barcodes) from query sequences in GenBank and BOLD system databases. This small number of identified OTUs may be explained by the paucity of molecular data from the Neotropical region. Our study provides valuable information about the characterization of macroinvertebrate families that are important biological indicators for the assessment of water quality in Nicaragua. The application of molecular approaches will allow the study of local diversity and further improve the application of molecular techniques for biomonitoring.

Coalescent-based species delimitation in North American pinyon pines using low-copy nuclear genes and plastomes.

PREMISE: Accurate species delimitation is essential for evolutionary biology, conservation, and biodiversity management. We studied species delimitation in North American pinyon pines, Pinus subsection Cembroides, a natural group with high levels of incomplete lineage sorting. METHODS: We used coalescent-based methods and multivariate analyses of low-copy number nuclear genes and nearly complete high-copy number plastomes generated with the Hyb-Seq method. The three coalescent-based species delimitation methods evaluated were the Generalized Mixed Yule Coalescent (GMYC), Poisson Tree Process (PTP), and Trinomial Distribution of Triplets (Tr2). We also measured admixture in populations with possible introgression. RESULTS: Our results show inconsistencies among GMYC, PTP, and Tr2. The single-locus based GMYC analysis of plastid DNA recovered a higher number of species (up to 24 entities, including singleton lineages and clusters) than PTP and the multi-locus coalescent approach. The PTP analysis identified 10 species whereas Tr2 recovered 13, which agreed closely with taxonomic treatments. CONCLUSIONS: We found that PTP and GMYC identified species with low levels of ILS and high morphological divergence (P. maximartinezii, P. pinceana, and P. rzedowskii). However, GMYC method oversplit species by identification of more divergent samples as singletons. Moreover, both PTP and GMYC were incapable of identifying some species that are readily identified morphologically. We suggest that the divergence times between lineages within North American pinyon pines are so disparate that GMYC results are unreliable. Results of the Tr2 method coincided well with previous delimitations based on morphology, DNA, geography, and secondary chemistry.

Species delimitation and microalgal cryptic diversity analysis of the genus Micractinium (Chlorophyta).

In this article, the system of the green microalgal genus Micractinium, based on morphological, physiological, ecological and molecular data, is considered. The main diagnostic species characteristics and the taxonomic placement of some taxa are also discussed. Phylogenetic analysis showed that the genus Micractinium is characterized by high cryptic diversity. The algorithms used for species delimitation had different results on the number of potentially species-level clusters allocated. The ABGD method was less "sensitive". The tree-based approaches GMYC and PTP showed a more feasible taxonomy of the genus Micractinium, being an effective additional tool for distinguishing species. The clustering obtained by the latter two methods is in good congruence with morphological (cell size and shape, ability to form colonies, production of bristles, chloroplast type), physiological (vitamin requirements, reaction to high and low temperatures), molecular (presence of introns, level of genetic differences, presence of CBCs or special features of the secondary structure in ITS1 and ITS2) and ecological characteristics (habitat). The polyphyly of the holotype of the genus M. pusillum as well as M. belenophorum is shown. The intron was effective as an additional tool for distinguishing species, and the results of the intron analysis should be taken into account together with other characteristics. The CBC approach, based on the search for compensatory base changes in conservative ITS2 regions, was successful only for distinguishing cryptic species from "true" members of M. pusillum. Therefore, to distinguish species, it is more effective to take into account all the CBC in ITS1 and ITS2 and analyze characteristic structural differences (molecular signatures) in the secondary structure of internal transcribed spacers. The genetic distances analysis of 18S-ITS1-5.8S-ITS2 nucleotide sequences showed that intraspecific differences in the genus ranged from 0 to 0.5 % and interspecific differences, from 0.6 to 4.7 %. Due to the polyphasic approach, it was possible to characterize 29 clusters and phylogenetic lines at the species level within the genus Micractinium and to make assumptions about the species.

Revalidation of Mazama rufa (Illiger 1815) (Artiodactyla: Cervidae) as a Distinct Species out of the Complex Mazama americana (Erxleben 1777).

The red brocket deer Mazama americana Erxleben, 1777 is considered a polyphyletic complex of cryptic species with wide chromosomal divergence. Evidence indicates that the observed chromosomal divergences result in reproductive isolation. The description of a neotype for M. americana allowed its genetic characterization and represented a comparative basis to resolve the taxonomic uncertainties of the group. Thus, we designated a neotype for the synonym Mazama rufa Illiger, 1815 and tested its recognition as a distinct species from the M. americana complex with the analysis of morphological, cytogenetic and molecular data. We also evaluated its distribution by sampling fecal DNA in the wild. Morphological data from craniometry and body biometry indicated an overlap of quantitative measurements between M. rufa and the entire M. americana complex. The phylogenetic hypothesis obtained through mtDNA confirmed the reciprocal monophyly relationship between M. americana and M. rufa, and both were identified as distinct molecular operational taxonomic units by the General Mixed Yule Coalescent species delimitation analysis. Finally, classic cytogenetic data and fluorescence in situ hybridization with whole chromosome painting probes showed M. rufa with a karyotype of 2n = 52, FN = 56. Comparative analysis indicate that at least fifteen rearrangements separate M. rufa and M. americana (sensu stricto) karyotypes, which confirmed their substantial chromosomal divergence. This divergence should represent an important reproductive barrier and allow its characterization as a distinct and valid species. Genetic analysis of fecal samples demonstrated a wide distribution of M. rufa in the South American continent through the Atlantic Forest, Cerrado and south region of Amazon. Thus, we conclude for the revalidation of M. rufa as a distinct species under the concept of biological isolation, with its karyotype as the main diagnostic character. The present work serves as a basis for the taxonomic review of the M. americana complex, which should be mainly based on cytogenetic characterization and directed towards a better sampling of the Amazon region, the evaluation of available names in the species synonymy and a multi-locus phylogenetic analysis.

Diaporthe amygdali, a species complex or a complex species?

Delimitation of species boundaries within the fungal genus Diaporthe has been challenging, but the analyses of combined multilocus DNA sequences has become an important tool to infer phylogenetic relationships and to circumscribe species. However, analyses of congruence between individual gene genealogies and the application of the genealogical concordance principle have been somehow overlooked. We noted that a group of species including D. amygdali, D. garethjonesii, D. sterilis, D. kadsurae, D. ternstroemia, D. ovoicicola, D. fusicola, D. chongqingensis and D. mediterranea, commonly known as D. amygdali complex, occupy a monophyletic clade in Diaporthe phylogenies but the limits of all species within the complex are not entirely clear. To assess the boundaries of species within this complex we employed the Genealogical Concordance Phylogenetic Species Recognition principle (GCPSR) and coalescence-based models: General Mixed Yule-Coalescent (GMYC) and Poisson Tree Processes (PTP). The incongruence detected between individual gene phylogenies, as well as the results of coalescent methods do not support the recognition of lineages within the complex as distinct species. Moreover, results support the absence of reproductive isolation and barriers to gene flow in this complex, thus providing further evidence that the D. amygdali species complex constitutes a single species. This study highlights the relevance of the application of the GCPSR principle, showing that concatenation analysis of multilocus DNA sequences, although being a powerful tool, might lead to an erroneous definition of species limits. Additionally, it further shows that coalescent methods are useful tools to assist in a more robust delimitation of species boundaries in the genus Diaporthe.

Species delimitation of Margattea cockroaches from China, with seven new species (Blattodea, Ectobiidae, Pseudophyllodromiinae).

Nearly 450 Margattea specimens were collected from 27 locations in China and their morphology was examined. Then 68 Margattea COI sequences were obtained and used to carry out phylogenetic analyses as well as species delimitation analyses using General Mixed Yule Coalescent (GMYC), Automatic Barcode Gap Discovery (ABGD), and Poisson Tree Processes (bPTP). GMYC analysis resulted in 21 molecular operational taxonomic units (MOTUs) (confidence interval: 20-22), which was completely consistent with the result of the bPTP. There were 15 MOTUs using the ABGD method. The number of MOTUs was slightly different from the assigned morphospecies (16). As to the incongruence between molecular and morphological results, we checked the specimens again and made sure that most morphological differences were determined to be intraspecific differences (except the difference between M. angusta and M. mckittrickae), although a large genetic distance existed. Finally, 16 Margattea species from China were defined in this study, of which, seven new species are established, i.e. Margattea deltodonta J-J He & Z-Q Wang, sp. nov., Margattea cuspidata J-J He & Z-Q Wang, sp. nov., Margattea caudata J-J He & Z-Q Wang, sp. nov., Margattea paratransversa J-J He & Z-Q Wang, sp. nov., Margattea disparilis J-J He & Z-Q Wang, sp. nov., Margattea transversa J-J He & Z-Q Wang, sp. nov., and Margattea bicruris J-J He & Z-Q Wang, sp. nov.

Factors affecting the efficiency of molecular species delimitation in a species-rich insect family.

In the context of global biodiversity loss, molecular species delimitation approaches can be very useful for accelerating species discovery through DNA taxonomy and inventory through DNA metabarcoding. In this study, the effect of some intrinsic factors on the efficiency of various single-marker species delimitation methods (fixed and variable nucleotide distance thresholds, ABGD, ASAP, GMYC, mPTP) was tested on more than 90 empirical data sets, derived from a set of 7,237 COI sequences attributed to 542 leaf beetles species (Coleoptera: Chrysomelidae). The considered factors were: (i) the number of haplotypes per species (as a proxy for genetic diversity), (ii) the geographic distance among conspecific collection localities (as a proxy of sampling width), (iii) the difficulty related to morphological identification of species, and (iv) the taxonomic rank. Distance-based methods, with on average more than 70% of match with morphological identification, outperformed those relying on phylogenetic trees, with less than 59%. A high number of haplotypes per species was found to have a negative effect on delimitation efficiency, whereas large geographic distances within species had a positive effect. All methods delimitations (except for GMYC) were significantly affected by the presence of species that are difficult to be identified, decreasing their efficiency. Finally, the only method influenced by the taxonomic rank of the data set was GMYC, showing lower efficiency in data sets at the genus than at higher levels. The observed biases we highlighted affecting efficiency could be accounted for when developing input data sets for species delimitation analyses to obtain a more reliable representation of biological diversity.

Use of species delimitation approaches to tackle the cryptic diversity of an assemblage of high Andean butterflies (Lepidoptera: Papilionoidea).

Cryptic biological diversity has generated ambiguity in taxonomic and evolutionary studies. Single-locus methods and other approaches for species delimitation are useful for addressing this challenge, enabling the practical processing of large numbers of samples for identification and inventory purposes. This study analyzed an assemblage of high Andean butterflies using DNA barcoding and compared the identifications based on the current morphological taxonomy with three methods of species delimitation (automatic barcode gap discovery, generalized mixed Yule coalescent model, and Poisson tree processes). Sixteen potential cryptic species were recognized using these three methods, representing a net richness increase of 11.3% in the assemblage. A well-studied taxon of the genus Vanessa, which has a wide geographical distribution, appeared with the potential cryptic species that had a higher genetic differentiation at the local level than at the continental level. The analyses were useful for identifying the potential cryptic species in Pedaliodes and Forsterinaria complexes, which also show differentiation along altitudinal and latitudinal gradients. This genetic assessment of an entire assemblage of high Andean butterflies (Papilionoidea) provides baseline information for future research in a region characterized by high rates of endemism and population isolation.

Two ways to be endemic. Alps and Apennines are different functional refugia during climatic cycles.

Endemics co-occur because they evolved in situ and persist regionally or because they evolved ex situ and later dispersed to shared habitats, generating evolutionary or ecological endemicity centres, respectively. We investigate whether different endemicity centres can intertwine in the region ranging from Alps to Sicily, by studying their butterfly fauna. We gathered an extensive occurrence data set for butterflies of the study area (27,123 records, 269 species, in cells of 0.5 × 0.5 degrees of latitude-longitude). We applied molecular-based delimitation methods (GMYC model) to 26,557 cytochrome c oxidase subunit 1 (COI) sequences of Western Palearctic butterflies. We identified entities based on molecular delimitations and/or the checklist of European butterflies and objectively attributed occurrences to their most probable entity. We obtained a zoogeographic regionalisation based on the 69 endemics of the area. Using phylogenetic ANOVA we tested if endemics from different centres differ from each other and from nonendemics for key ecological traits and divergence time. Endemicity showed high incidence in the Alps and Southern Italy. The regionalisation separated the Alps from the Italian Peninsula and Sicily. The endemics of different centres showed a high turnover and differed in phylogenetic distances, phenology and distribution traits. Endemics are on average younger than nonendemics and the Peninsula-Sicily endemics also have lower variance in divergence than those from the Alps. The observed variation identifies Alpine endemics as paleoendemics, now occupying an ecological centre, and the Peninsula-Sicily ones as neoendemics, that diverged in the region since the Pleistocene. The results challenge the common view of the Alpine-Apennine area as a single "Italian refugium".

Molecular diversity and species delimitation in the family Gasteruptiidae (Hymenoptera: Evanioidea).

Gasteruptiidae Ashmead is an easily recognised family of wasps with ∼589 described species worldwide. Although well characterised by traditional taxonomy, multiple authors have commented on the extreme morphological uniformity of the group, making species-level identification difficult. This problem is enhanced by the lack of molecular data and molecular phylogenetic research for the group. We used 187 cytochrome c oxidase subunit I (COI) barcodes to explore the efficiency of sequence data to delimitate species in Gasteruptiidae. We undertook a graphical and discussion-based comparison of six methods for species delimitation, with the success of methods judged based on known species boundaries and morphology. Both distance-based (ABGD and jMOTU threshold analysis) and tree-based (GMYC and PTP) methods compared across multiple parameters recovered variable molecular operational taxonomic units (MOTUs), ranging from 55 to 123 MOTUs. Tree-based methods tended to split known morphological species less than distance-based methods, with the single-threshold GMYC method the most concordant with known morphospecies. Our results suggest that the incorporation of molecular species delimitation techniques provides a powerful tool to assist in the interpretation of species and help direct informed decisions with taxonomic uncertainty in the family.

Cryptic speciation in populations of the genus Aphyocharax (Characiformes: Characidae) from eastern Amazon coastal river drainages and surroundings revealed by single locus species delimitation methods

Recent studies in eastern Amazon coastal drainages and their surroundings have revealed new fish species that sometimes exhibit little morphological differentiation (cryptic species). Thus, we used a DNA-based species delimitation approach to test if populations showing the morphotype and typical character states of the Aphyocharax avary holotype correspond either to A. avary or A. brevicaudatus, two known species from the region, or if they form independent lineages, indicating cryptic speciation. WP and GMYC analyses recovered five lineages (species) in the ingroup, while a bPTP analysis delimited three lineages. ABGD analyses produced two possible results: one corroborating the WP and GMYC methods and another corroborating the bPTP method. All methods indicate undescribed cryptic species in the region and show variation from at least 1 to 4 species in the ingroup, depending on the approach, corroborating previous studies, and revealing this region as a possible hotspot for discovering undescribed fish species.(AU)

Estudos recentes nas drenagens costeiras da Amazônia oriental e seus arredores revelaram novas espécies de peixes que às vezes exibem pouca diferenciação morfológica (espécies crípticas). Assim, usamos uma abordagem de delimitação de espécies baseada em DNA para testar se as populações que apresentam o morfotipo e os estados de caráter típicos do holótipo Aphyocharax avary correspondem a A. avary ou A. brevicaudatus, duas espécies conhecidas da região, ou se formam linhagens independentes, indicando especiação críptica. As análises de WP e GMYC recuperaram cinco linhagens (espécies) no grupo interno, enquanto uma análise de bPTP delimitou três linhagens. As análises ABGD produziram dois resultados possíveis: um corroborando os métodos WP e GMYC e outro corroborando o método bPTP. Todos os métodos indicam espécies crípticas não descritas na região e apresentam variação de pelo menos uma a quatro espécies no grupo interno, dependendo da abordagem, corroborando estudos anteriores, e revelando esta região como um possível "hotspot" para descoberta de espécies de peixes não descritas.(AU)

Three new species of Rugitermes (Isoptera, Kalotermitidae) from Peru and Bolivia.

The soldier of Rugitermes aridus sp. nov. is described from a xeric, termite-depauperate region of central Peru. Rugitermes rufus sp. nov. and R. volcanensis sp. nov. are described from soldiers and dealated imagos collected in a mesic forest of Amboró National Park in western Bolivia. The imago of R. rufus is unique among all described Rugitermes species in that the head capsule is reddish orange and the pronotum is brown. The imago head and pronotum are both brown in R. volcanensis. A phylogenetic and GMYC barcode analyses were performed with the COI gene. These analyses confirmed the three new species and revealed a high undescribed diversity of Rugitermes in the New World.

Cryptic diversity and species boundaries within the Paragalago zanzibaricus species complex.

The recently described genus Paragalago is a complex of several nocturnal and morphologically cryptic species distributed in the forests of eastern Africa. Species diversity within this genus has been mainly described using species-specific differences in their loud calls. However, molecular data are still lacking for this group and species boundaries remain unclear. In this study, we explore species diversity within the zanzibaricus-complex using a combination of mitochondrial and nuclear data and comparing multiple species delimitation methods. Our results consistently support the existence of three independent lineages, P. cocos, P. zanzibaricus, and P. granti, confirming previous hypotheses based on vocal data. We conclude that these three lineages represent valid cryptic species and we hypothesize that speciation within this complex was characterized by cycles of forest expansion and contraction in the Plio-Pleistocene.

A molecular assessment of species diversity in Tympanopleura and Ageneiosus catfishes (Auchenipteridae: Siluriformes).

In order to test the congruence of genetic data to the morphologically defined Neotropical catfish genera Tympanopleura and Ageneiosus and explore species diversity, we generated 17 DNA barcodes from five of six species of Tympanopleura and 12 of 13 species of Ageneiosus. To discriminate limits between species, an automatic barcode gap discovery (ABGD), a generalised mixed yule-coalescent model (GYMC) and fixed distance thresholds Kimura two-parameter (K2P; 3%) were used to discriminate putative species limits from the DNA barcodes. The ABGD, GMYC and K2P methods agreed by each generating 13 clusters: six in Tympanopleura (five nominal plus one undescribed species) and seven in Ageneiosus. These clusters corresponded broadly to the described species, except in the case of the Ageneiosus ucayalensis group (A. akamai, A. dentatus, A. intrusus, A. ucayalensis, A. uranophthalmus and A. vittatus). Haplotype sharing and low divergences may have prevented molecular methods from distinguishing these species. We hypothesise that this is the result of a recent radiation of a sympatric species group distributed throughout the Amazon Basin. One putative new species of Tympanopleura was also supported by the molecular data. These results taken together highlight the utility of molecular methods such as DNA barcoding in understanding patterns of diversification across large geographic areas and in recognising overlooked diversity.

Cryptic species delineation in freshwater planarians of the genus Dugesia (Platyhelminthes, Tricladida): Extreme intraindividual genetic diversity, morphological stasis, and karyological variability.

The keystone of planarian taxonomy traditionally has been the anatomy of the copulatory apparatus. However, many planarian species comprise asexual fissiparous populations, with the fissiparous animals not developing a copulatory apparatus, thus precluding their morphological identification. Incorporation of molecular data into planarian systematics has been of great value, not only in the identification of fissiparous individuals but also as an additional source of information for determining species boundaries. Nevertheless, the discrepancy between morphological and molecular data has highlighted the need for extra sources of taxonomic information. Moreover, a recent study has pointed out that fissiparous reproduction may lead to high levels of intraindividual genetic diversity in planarians, which may mislead molecular analyses. In the present study we aim to test a new up-to-date integrative taxonomic procedure for planarians, including intraindividual genetic data and additional sources of taxonomic information, besides morphology and DNA, using Dugesia subtentaculata sensu lato as a model organism, a species with an intricate taxonomic history. First, we used three different methods for molecular species delimitation on single locus datasets, both with and without intraindividual information, for formulating Primary Species Hypotheses (PSHs). Subsequently, Secondary Species Hypotheses (SSHs) were formulated on the basis of three types of information: (1) a coalescent-based species delimitation method applied to multilocus data, (2) morphology of the copulatory apparatus, and (3) karyological metrics. This resulted in the delimitation of four morphologically cryptic species within the nominal species D. subtentaculata. Our results provide evidence that the analysis of intraindividual genetic data is essential for properly developing PSHs in planarians. Our study reveals also that karyological differentiation, rather than morphological differentiation, may play an important role in speciation processes in planarians, thus suggesting that the currently known diversity of the group could be highly underestimated.