Phylogenetic estimation error can decrease the accuracy of species delimitation: a Bayesian implementation of the general mixed Yule-coalescent model

Taxonomic inflation due to inadequate sampling: are girdled lizards (Cordylus minor species complex) from the Great Karoo one and the same?

The Great Karoo and Namaqualand of South Africa are home to a species complex of morphologically conserved lizards that occur in allopatry (Karoo: Cordylus aridus, Cordylus cloetei, Cordylus minor; Namaqualand: Cordylus imkeae). However, there are negligible morphological differences and a lack of obvious physical or climatic barriers, particularly among the three Karoo species. We hypothesized that poor geographic coverage in previous studies and lack of an explicit species concept has caused taxonomic inflation. We therefore tested species boundaries by examining multiple criteria: multi-gene phylogenetics, niche distribution modelling and re-examination of diagnostic morphological features with a larger sample size. We found that C. aridus, C. cloetei and C. minor lack diagnosable differences for both genetics and morphology. Distribution modelling, ranging from present day to the last interglacial period, show connectivity has been maintained especially during cooler periods. Conversely, C. imkeae is morphologically diagnosable, genetically distinct and lacks connectivity with the other taxa. By evaluating multiple operational criteria, we conclude that the C. minor species complex comprises only two species, C. minor (with C. aridus and C. cloetei as junior synonyms) and C. imkeae, demonstrating that species defined from inadequate data and lack of an explicit species concept can lead to taxonomic inflation.

Recent speciation and hybridization in Icelandic deep-sea isopods: An integrative approach using genomics and proteomics

The crustacean marine isopod species Haploniscus bicuspis (Sars, 1877) shows circum-Icelandic distribution in a wide range of environmental conditions and along well-known geographic barriers, such as the Greenland-Iceland-Faroe (GIF) Ridge. We wanted to explore population genetics, phylogeography and cryptic speciation as well as investigate whether previously described, but unaccepted subspecies have any merit. Using the same set of specimens, we combined mitochondrial COI sequences, thousands of nuclear loci (ddRAD), and proteomic profiles, plus selected morphological characters using confocal laser scanning microscopy (CLSM). Five divergent genetic lineages were identified by COI and ddRAD, two south and three north of the GIF Ridge. Assignment of populations to the three northern lineages varied and detailed analyses revealed hybridization and gene flow between them, suggesting a single northern species with a complex phylogeographic history. No apparent hybridization was observed among lineages south of the GIF Ridge, inferring the existence of two more species. Differences in proteomic profiles between the three putative species were minimal, implying an ongoing or recent speciation process. Population differentiation was high, even among closely associated populations, and higher in mitochondrial COI than nuclear ddRAD loci. Gene flow is apparently male-biased, leading to hybrid zones and instances of complete exchange of the local nuclear genome through immigrating males. This study did not confirm the existence of subspecies defined by male characters, which probably instead refer to different male developmental stages.

Molecular species delimitation of shrub frogs of the genus Pseudophilautus (Anura, Rhacophoridae)

Sri Lanka is an amphibian hotspot of global significance. Its anuran fauna is dominated by the shrub frogs of the genus Pseudophilautus. Except for one small clade of four species in Peninsular India, these cool-wet adapted frogs, numbering some 59 extant species, are distributed mainly across the montane and lowland rain forests of the island. With species described primarily by morphological means, the diversification has never yet been subjected to a molecular species delimitation analysis, a procedure now routinely applied in taxonomy. Here we test the species boundaries of Pseudophilautus in the context of the phylogenetic species concept (PSC). We use all the putative species for which credible molecular data are available (nDNA-Rag-1; mt-DNA- 12S rRNA, 16S rRNA) to build a well resolved phylogeny, which is subjected to species delimitation analyses. The ABGD, bPTP, mPTP and bGMYC species delimitation methods applied to the 16S rRNA frog barcoding gene (for all species), 12S rRNA and Rag-1 nDNA grouped P. procax and P. abundus; P. hallidayi and P. fergusonianus; P. reticulatus and P. pappilosus; P. pleurotaenia and P. hoipolloi; P. hoffmani and P. asankai; P. silvaticus and P. limbus; P. dilmah and P. hankeni; P. fulvus and P. silus.. Surprisingly, all analyses recovered 14 unidentified potential new species as well. The geophylogeny affirms a distribution across the island's aseasonal 'wet zone' and its three principal hill ranges, suggestive of allopatric speciation playing a dominant role, especially between mountain masses. Among the species that are merged by the delimitation analyses, a pattern leading towards a model of parapatric speciation emerges-ongoing speciation in the presence of gene flow. This delimitation analysis reinforces the species hypotheses, paving the way to a reasonable understanding of Sri Lankan Pseudophilautus, enabling both deeper analyses and conservation efforts of this remarkable diversification. http://zoobank.org/urn:lsid:zoobank.org:pub:DA869B6B-870A-4ED3-BF5D-5AA3F69DDD27.

Ancient volcanos as species pumps: A case study of freshwater amphipods in Northeast Asia

Volcano-tectonic processes have been viewed as primary drivers in the formation of present-day diversity. Volcanos associated with mountain uplifts drive allopatric speciation through vicariance and may impact the surrounding areas like species pump or species attractor. However, the application of these hypotheses to aquatic fauna has rarely been tested explicitly. We tested these hypotheses in the Changbai Mountains (Mts), which are one of the most typical, active volcanic ranges in Northeast (NE) Asia with a long and turbulent geological history. The Gammarus nekkensis species complex of amphipod crustaceans, widely distributed throughout NE Asia with poor dispersal abilities and a long evolutionary history, is a suitable model for testing hypotheses of species pump or species attractor. Phylogenetic and ancestral range reconstructions demonstrated that the studied amphipod originated from the Changbai Mts ~27 Ma and diverged into eastern (Clade I) and western (Clade II) clades, which corresponds well with the initial volcanic eruption of the Changbai Mts in the Late Oligocene. The subsequent diversifications of subclades CI-3, CII-1a and CII-2a were probably driven by second and third eruptions of the Changbai Mts during the Miocene. In particular, the Changbai lineages had spread to the Russian Far East multiple times since the Early Miocene, and widely colonized the region during the Pleistocene. Our discoveries suggest that the ancient volcanos of the Changbai Mts act as species pumps in NE Asia, resulted in burst of diversification around the Changbai Mts and subsequent dispersals into adjacent regions.

Hidden diversity revealed in the freshwater snails, Bythinella and Pseudamnicola, in the Island of Crete

Crete with its complex geomorphological history is the island with the highest number of endemism observed in animal and plant taxa throughout the Aegean archipelago. While other groups of organisms within Crete are well-studied, the freshwater gastropod fauna still remains poorly investigated. Bythinella and Pseudamnicola, are 2 genera of freshwater springsnails, both present on the island, inhabiting springs and other freshwater habitats. Here, we conduct a comprehensive study on the distribution of the different genetic lineages of the 2 gastropod genera in order to assess the mode of their differentiation on the island and infer the actual number of species present in the island. Towards these aims, sequence data from the mitochondrial gene were used and analyzed within a phylogenetic framework. For Bythinella, our results strongly support at least 5 delineated Bythinella spp. inhabiting Crete, which correspond to the already described species from previous studies with the addition of a new one. Bythinella analyses reveal an old time-frame of differentiation with vicariant phenomena being more likely the main drivers shaping the present-day distribution of the genus' genetic lineages. For Pseudamnicola, our data indicate the presence of at least 2 delineated Pseudamnicola spp. with a differentiation more consistent to an isolation-by-distance pattern of a relatively recent origin. Dispersion processes followed by isolation of the populations and/or recent speciation, seem to be the underlying process for the current distribution of Pseudamnicola lineages.

Diversity and taxonomic revision of tribes Rhipileae and Rhipiliopsideae (Halimedaceae, Chlorophyta) based on molecular and morphological data

Genera and species of the tribes Rhipileae and Rhipiliopsideae are abundant in most coral reef ecosystems worldwide. However, the group has been largely overlooked, and very little genetic data is available to accurately assess its diversity, phylogenetic relationships, and geographic distribution. Our study provided an in-depth reassessment of tribes Rhipileae and Rhipiliopsideae based on a species-rich dataset and the combination of molecular species delimitation, multilocus phylogenetic analyses (tufA, rbcL, and 18S rDNA), and morpho-anatomic observations. Our results revealed an unexpected diversity of 38 morphologically validated species hypotheses, including 20 new species, two of which are described in this paper and one resurrected species (Rhipilia diaphana). Based on our phylogenetic results we proposed to redefine the genera Rhipilia and Rhipiliopsis and described two new genera, Kraftalia gen. nov. (Rhipileae) and Rhipiliospina gen. nov. (Rhipiliopsideae). Finally, we validated Rhipiliella Kraft and included it in the tribe Rhipileae. Although Rhipilia and Rhipiliopsis have a pantropical distribution, none of the species studied here appeared cosmopolitan; instead, they have restricted distributions.

Blue mussels of the Mytilus edulis species complex from South America: The application of species delimitation models to DNA sequence variation

Smooth-shelled blue mussels, Mytilus spp., have a worldwide antitropical distribution and are ecologically and economically important. Mussels of the Mytilus edulis species complex have been the focus of numerous taxonomic and biogeographical studies, in particular in the Northern hemisphere, but the taxonomic classification of mussels from South America remains unclear. The present study analysed 348 mussels from 20 sites in Argentina, Chile, Uruguay and the Falkland Islands on the Atlantic and Pacific coasts of South America. We sequenced two mitochondrial locus, Cytochrome c Oxidase subunit I (625 bp) and 16S rDNA (443 bp), and one nuclear gene, ribosomal 18S rDNA (1770 bp). Mitochondrial and nuclear loci were analysed separately and in combination using maximum likelihood and Bayesian inference methods to identify the combination of the most informative dataset and model. Species delimitation using five different models (GMYC single, bGMYC, PTP, bPTP and BPP) revealed that the Mytilus edulis complex in South America is represented by three species: native M. chilensis, M. edulis, and introduced Northern Hemisphere M. galloprovincialis. However, all models failed to delimit the putative species Mytilus platensis. In contrast, however, broad spatial scale genetic structure in South America using Geneland software to analyse COI sequence variation revealed a group of native mussels (putatively M. platensis) in central Argentina and the Falkland Islands. We discuss the scope of species delimitation methods and the use of nuclear and mitochondrial genetic data to the recognition of species within the Mytilus edulis complex at regional and global scales.

Subtle environmental variation affects phenotypic differentiation of shallow divergent treefrog lineages in Amazonia

Amazonia harbours a vast biotic and ecological diversity, enabling investigation of the effects of microevolutionary processes and environmental variation on species diversification. Integrative approaches combining phenotypic and genetic variation can improve our knowledge on diversification processes in megadiverse regions. Here, we investigate the influence of environmental and geographic variation on the genetic and morphological differentiation in the Amazonian Boana calcarata-fasciata (Anura: Hylidae) species complex. We analysed the variation of one mtDNA gene from individuals of different forest environments, and assessed their phylogenetic relationships and species limits to define the lineages to perform a phenotypic-environmental approach. We collected morphological data (head shape and size) using 3D models and investigated the phylogenetic signal, evolutionary model and influence of environmental variables on morphology. We verified associations between environmental and geographical distances with morphological and genetic variation using distance-based redundancy analyses and Mantel tests. We found an even higher cryptic diversity than already recognized within the species complex. Body size and head shape varied among specimens, but did not present phylogenetic signal, diverging under a selective evolutionary model. Our results show that diverse factors have influenced morphological and genetic variation, but environmental conditions such as vegetation cover, precipitation and climate change velocity influenced morphological diversification. Possible population-level mechanisms such as parallel morphological evolution or plastic responses to similar environments could account for such patterns in these typical Amazonian treefrogs.

Phylogeography and molecular species delimitation reveal cryptic diversity in Potamolithus (Caenogastropoda: Tateidae) of the southwest basin of the Andes

The species of the genus Potamolithus inhabiting the southwestern basin of the Andes are difficult to distinguish due to small size and similar shell morphology. Only Potamolithus australis and Potamolithus santiagensis have been traditionally recognized in this region, but the occurrence of several morphologically similar undescribed populations could increase the regional richness. Here we delimit described and potentially undescribed cryptic species of the genus using partial sequences of the mitochondrial cytochrome c oxidase subunit I (COI) gene. Network analysis and diversity indices inferred six highly differentiated haplogroups, many of them sympatric and widespread in the study area. Phylogeographic analyses suggest a scenario of recent diversification and the occurrence of multiple refuges during the successive Pleistocene glaciations. Phylogenetic analysis also recovered six major clades that showed no relationship with physiography. Species delimitation analyses consistently recognized three or four candidate species apart from P. australis and P. santiagensis. Divergence times indicate that speciation of Chilean Potamolithus began at the end of the Pliocene, probably driven by climatic rather than geographic events. Considering the high inter- and intra-basin genetic diversity, conservation efforts should be focused on protecting sympatric taxa in the basins with the highest species richness.

Excluding spatial sampling bias does not eliminate oversplitting in DNA-based species delimitation analyses

DNA barcoding and DNA-based species delimitation are major tools in DNA taxonomy. Sampling has been a central debate in this context, because the geographical composition of samples affects the accuracy and performance of DNA barcoding. Performance of complex DNA-based species delimitation is to be tested under simpler conditions in absence of geographic sampling bias. Here, we present an empirical dataset sampled from a single locality in a Southeast-Asian biodiversity hotspot (Laos: Phou Pan mountain). We investigate the performance of various species delimitation approaches on a megadiverse assemblage of herbivorous chafer beetles (Coleoptera: Scarabaeidae) to infer whether species delimitation suffers in the same way from exaggerate infraspecific variation despite the lack of geographic genetic variation that led to inconsistencies between entities from DNA-based and morphology-based species inference in previous studies. For this purpose, a 658 bp fragment of the mitochondrial cytochrome c oxidase subunit 1 (cox1) was analyzed for a total of 186 individuals of 56 morphospecies. Tree-based and distance-based species delimitation methods were used. All approaches showed a rather limited match ratio (max. 77%) with morphospecies. Poisson tree process (PTP) and statistical parsimony network analysis (TCS) prevailingly over-splitted morphospecies, while 3% clustering and Automatic Barcode Gap Discovery (ABGD) also lumped several species into one entity. ABGD revealed the highest congruence between molecular operational taxonomic units (MOTUs) and morphospecies. Disagreements between morphospecies and MOTUs have to be explained by historically acquired geographic genetic differentiation, incomplete lineage sorting, and hybridization. The study once again highlights how important morphology still is in order to correctly interpret the results of molecular species delimitation.

The hidden diversity of the genus Lyurella Derzhavin, 1939 (Crustacea: Amphipoda: Crangonyctidae): four new species from the subterranean habitats of the northwestern Caucasus, Russia

Four new species of the Palaearctic crangonyctid amphipod genus Lyurella Derzhavin, 1939 (Crustacea: Amphipoda: Crangonyctidae), L. mikhailovi sp. n., L. fanagorica sp. n., L. fontinalis sp. n. and L. asheensis sp. n., are described based on an integrative approach from the subterranean habitats of the southwestern foothills of the Greater Caucasian Ridge (the north-eastern Black Sea coast). Despite the relative proximity of the habitats, the interspecific genetic divergence (by COI mtDNA gene marker) between the newly outlined Caucasian species of the genus varied from 11 to 21%, demonstrating a long-term isolation and lack of gene flow for at least 37Mya, starting from the Pliocene. The lowest genetic divergence between L. shepsiensis Sidorov, 2015 and L. asheensis sp. n., estimated as 4%, is also considered species-specific due to the presence of distinct morphological differences. We discuss the phylogeny, morphology, and distribution and provide a key for all known species of Lyurella. DNA barcoding data for all species, including the type species of the genus, Lyurella hyrcana Derzhavin, 1939, are presented for the first time.

Multiple evolutionary lineages for the main vector of Leishmania guyanensis, Lutzomyia umbratilis (Diptera: Psychodidae), in the Brazilian Amazon

Lutzomyia umbratilis is the main vector of Leishmania guyanensis in the Brazilian Amazon and in neighboring countries. Previous biological and molecular investigations have revealed significant differences between L. umbratilis populations from the central Brazilian Amazon region. Here, a phylogeographic survey of L. umbratilis populations collected from nine localities in the Brazilian Amazon was conducted using two mitochondrial genes. Statistical analyses focused on population genetics, phylogenetic relationships and species delimitations. COI genetic diversity was very high, whereas Cytb diversity was moderate. COI genealogical haplotypes, population structure and phylogenetic analyses identified a deep genetic differentiation and three main genetic groups. Cytb showed a shallower genetic structure, two main haplogroups and poorly resolved phylogenetic trees. These findings, allied to absence of isolation by distance, support the hypothesis that the Amazon and Negro Rivers and interfluves are the main evolutionary forces driving L. umbratilis diversification. The main three genetic groups observed represent three evolutionary lineages, possibly species. The first lineage occurs north of the Amazon River and east of Negro River, where Le. guyanensis transmission is intense, implying that L. umbratilis is an important vector there. The second lineage is in the interfluve between north of Amazon River and west of Negro River, an area reported to be free of Le. guyanensis transmission. The third lineage, first recorded in this study, is in the interfluve between south of Amazonas River and west of Madeira River, and its involvement in the transmission of this parasite remains to be elucidated.

Cryptic diversity and diversification processes in three cis-Andean Rhamdia species (Siluriformes: Heptapteridae) revealed by DNA barcoding

The wide distribution of the Neotropical freshwater catfish Rhamdia offers an excellent opportunity to investigate the historical processes responsible for modeling South America’s hydrogeological structure. We used sequences from cis-Andean and Mesoamerican Rhamdia species to reconstruct and estimate divergence times among cis-Andean lineages, correlating the results with known geological events. Species delimitation methods based on distance (DNA barcoding and BIN) and coalescence (GMYC) approaches identified nine well-supported lineages from the cis-Andean region from sequences available in the BOLD dataset. The cis-Andean Rhamdia lineages diversification process began in Eocene and represented the split between cis-Andean and Mesoamerican clades. The cis-Andean clade contains two principal groups: Northwest clade (MOTUs from Amazon, Essequibo, Paraguay, and Itapecuru basins) and Southeast clade (Eastern Brazilian shield basins (Paraná, Uruguay, Iguaçu, and São Francisco) plus eastern coastal basins). The diversification of the cis-Andean Rhamdia lineages results from vicariance and geodispersion events, which played a key role in the current intricate distribution pattern of the Rhamdia lineages. The wide geographical distribution and large size of the specimens make it attractive to cultivate in different countries of the Neotropical region. The lineages delimitation minimizes identification mistakes, unintentional crossings by aquaculture, and reduces natural stocks contamination.

A taxonomist's nightmare - Cryptic diversity in Caribbean intertidal arthropods (Arachnida, Acari, Oribatida)

There has been a long controversy about what defines a species and how to delimitate them which resulted in the existence of more than two dozen different species concepts. Recent research on so-called "cryptic species" heated up this debate as some scientists argue that these cryptic species are only a result of incompatible species concepts. While this may be true, we should keep in mind that all concepts are nothing more than human constructs and that the phenomenon of high phenotypic similarity despite reproductive isolation is real. To investigate and understand this phenomenon it is important to classify and name cryptic species as it allows to communicate them with other fields of science that use Linnaean binomials. To provide a common framework for the description of cryptic species, we propose a possible protocol of how to formally name and describe these taxa in practice. The most important point of this protocol is to explain which species concept was used to delimitate the cryptic taxon. As a model, we present the case of the allegedly widespread Caribbean intertidal mite Thalassozetes barbara, which in fact consists of seven phenotypically very similar but genetically distinct species. All species are island or short-range endemics with poor dispersal abilities that have evolved in geographic isolation. Stabilizing selection caused by the extreme conditions of the intertidal environment is suggested to be responsible for the morphological stasis of this cryptic species complex.

Divergence time and species delimitation of microbivalves in the Southern Ocean: the case of Kidderia species

The systematics of Subantarctic and Antarctic near-shore marine benthic invertebrates requires major revision and highlights the necessity to incorporate additional sources of information in the specimen identification chart in the Southern Ocean (SO). In this study, we aim to improve our understanding of the biodiversity of Kidderia (Dall 1876) through molecular and morphological comparisons of Antarctic and Subantarctic taxa. The microbivalves of the genus Kidderia are small brooding organisms that inhabit intertidal and shallow subtidal rocky ecosystems. This genus represents an interesting model to test the vicariance and dispersal hypothesis in the biogeography of the SO. However, the description of Kidderia species relies on a few morphological characters and biogeographic records that raise questions about the true diversity in the group. Here we will define the specimens collected with genetic tools, delimiting their respective boundaries across provinces of the SO, validating the presence of two species of Kidderia. Through the revision of taxonomic issues and species delimitation, it was possible to report that the Antarctic species is Kidderia subquadrata and the species recorded in the Subantarctic islands Diego Ramirez, South Georgia and the Kerguelen Archipelago is Kidderia minuta. The divergence time estimation suggests the origin and diversification of Kidderia lineages are related to historical vicariant processes probably associated with the separation of the continental landmasses close to the late Eocene.

Re-examination of species limits in Aspergillus section Flavipedes using advanced species delimitation methods and description of four new species

Since the last revision in 2015, the taxonomy of section Flavipedes evolved rapidly along with the availability of new species delimitation techniques. This study aims to re-evaluate the species boundaries of section Flavipedes members using modern delimitation methods applied to an extended set of strains (n = 90) collected from various environments. The analysis used DNA sequences of three house-keeping genes (benA, CaM, RPB2) and consisted of two steps: application of several single-locus (GMYC, bGMYC, PTP, bPTP) and multi-locus (STACEY) species delimitation methods to sort the isolates into putative species, which were subsequently validated using DELINEATE software that was applied for the first time in fungal taxonomy. As a result, four new species are introduced, i.e. A. alboluteus, A. alboviridis, A. inusitatus and A. lanuginosus, and A. capensis is synonymized with A. iizukae. Phenotypic analyses were performed for the new species and their relatives, and the results showed that the growth parameters at different temperatures and colonies characteristics were useful for differentiation of these taxa. The revised section harbors 18 species, most of them are known from soil. However, the most common species from the section are ecologically diverse, occurring in the indoor environment (six species), clinical samples (five species), food and feed (four species), droppings (four species) and other less common substrates/environments. Due to the occurrence of section Flavipedes species in the clinical material/hospital environment, we also evaluated the susceptibility of 67 strains to six antifungals (amphotericin B, itraconazole, posaconazole, voriconazole, isavuconazole, terbinafine) using the reference EUCAST method. These results showed some potentially clinically relevant differences in susceptibility between species. For example, MICs higher than those observed for A. fumigatus wild-type were found for both triazoles and amphotericin B for A. ardalensis, A. iizukae, and A. spelaeus whereas A. lanuginosus, A. luppiae, A. movilensis, A. neoflavipes, A. olivimuriae and A. suttoniae were comparable to or more susceptible as A. fumigatus. Finally, terbinafine was in vitro active against all species except A. alboviridis.

Molecular Phylogeny of Cryptonanus (Didelphidae: Thylamyini): Evidence for a recent and complex diversification in South American open biomes

Systematic revisions of South American marsupials have contributed to our knowledge about genus and species diversity in the last decades, including studies of the most recently described genus Cryptonanus (Didelphidae), currently comprising four recognized species. Herein we provide the first phylogeny for these mouse opossums based on comprehensive sampling, including representatives from all nominal taxa, encompassing most of the geographic distribution of the genus while also extending its known range. The taxonomic status of Cryptonanus species was explored by analyses of multiple mitochondrial and nuclear DNA markers to assess phylogenetic relationships and to provide divergence time estimates, species delimitations and biogeographical hypotheses. Cryptonanus monophyly remained highly supported despite the inclusion of abundant new data from more than a hundred specimens, comprising 10 independent evolutionary lineages. Species-complexes within valid nominal taxa reveal higher species richness in the genus. Based on divergence estimates from a dated phylogeny, we suggest that Cryptonanus diversified along the Quaternary, with speciation events occurring well into the Pleistocene. The best supported biogeographical hypothesis endorses speciation by vicariance and subset speciation across open formations in shaping the evolutionary history of this didelphid genus, strongly associated with dry tropical landscapes of South America.

Unexpectedly high levels of lineage diversity in Sundaland puddle frogs (Dicroglossidae: Occidozyga Kuhl and van Hasselt, 1822)

One of the most urgent contemporary tasks for taxonomists and evolutionary biologists is to estimate the number of species on earth. Recording alpha diversity is crucial for protecting biodiversity, especially in areas of elevated species richness, which coincide geographically with increased anthropogenic environmental pressures - the world's so-called biodiversity hotspots. Although the distribution of Puddle frogs of the genus Occidozyga in South and Southeast Asia includes five biodiversity hotspots, the available data on phylogeny, species diversity, and biogeography are surprisingly patchy. Samples analyzed in this study were collected throughout Southeast Asia, with a primary focus on Sundaland and the Philippines. A mitochondrial gene region comprising ~ 2000 bp of 12S and 16S rRNA with intervening tRNA Valine and three nuclear loci (BDNF, NTF3, POMC) were analyzed to obtain a robust, time-calibrated phylogenetic hypothesis. We found a surprisingly high level of genetic diversity within Occidozyga, based on uncorrected p-distance values corroborated by species delimitation analyses. This extensive genetic diversity revealed 29 evolutionary lineages, defined by the > 5% uncorrected p-distance criterion for the 16S rRNA gene, suggesting that species diversity in this clade of phenotypically homogeneous forms probably has been underestimated. The comparison with results of other anuran groups leads to the assumption that anuran species diversity could still be substantially underestimated in Southeast Asia in general. Many genetically divergent lineages of frogs are phenotypically similar, indicating a tendency towards extensive morphological conservatism. We present a biogeographic reconstruction of the colonization of Sundaland and nearby islands which, together with our temporal framework, suggests that lineage diversification centered on the landmasses of the northern Sunda Shelf. This remarkably genetically structured group of amphibians could represent an exceptional case for future studies of geographical structure and diversification in a widespread anuran clade spanning some of the most pronounced geographical barriers on the planet (e.g., Wallace's Line). Studies considering gene flow, morphology, ecological and bioacoustic data are needed to answer these questions and to test whether observed diversity of Puddle frog lineages warrants taxonomic recognition.

DNA barcoding of Austrian snow scorpionflies (Mecoptera, Boreidae) reveals potential cryptic diversity in Boreus westwoodi

Background

Snow scorpionflies (genus Boreus) belong to a family of Mecoptera, Boreidae, that has been vastly neglected by entomological researchers due to their shift in seasonality to the winter months. Their activity during this time is regarded as a strategy for predator avoidance and regular sightings on snow fields suggest that this also facilitates dispersal. However, many aspects about snow scorpionflies, especially systematics, taxonomy, distribution of species, phylogenetics and phylogeography have remained fairly unexplored until today. In this study, we fill some of these gaps by generating a reference DNA barcode database for Austrian snow scorpionflies in the frame of the Austrian Barcode of Life initiative and by characterising morphological diversity in the study region.

Methods

Initial species assignment of all 67 specimens was based on male morphological characters previously reported to differ between Boreus species and, for females, the shape of the ovipositor. DNA barcoding of the mitochondrial cytochrome c oxidase subunit 1 (COI) gene was carried out for all 67 samples and served as a basis for BIN assignment, genetic distance calculations, as well as alternative species delimitation analyses (ABGD, GMYC, bGMYC, bPTP) and a statistical parsimony network to infer phylogenetic relationships among individual samples/sampling sites.

Results

Morphological investigations suggested the presence of both Boreus hyemalis and Boreus westwoodi in Austria. DNA barcoding also separated the two species, but resulted in several divergent clades, the paraphyly of B. westwoodi in Austria, and high levels of phylogeographic structure on a small geographic scale. Even though the different molecular species delimitation methods disagreed on the exact number of species, they unequivocally suggested the presence of more than the traditionally recognized two Boreus species in Austria, thus indicating potential cryptic species within the genus Boreus in general and especially in B. westwoodi.

Bosminopsis deitersi (Crustacea: Cladocera) as an ancient species group: a revision

Water fleas (Crustacea: Cladocera) of the Family Bosminidae have been studied since the founding of paleolimnology and freshwater ecology. However, one species, Bosminopsis deitersi, stands out for its exceptional multicontinental range and broad ecological requirements. Here we use an integrated morphological and multilocus genetic approach to address the species problem in B. deitersi. We analyzed 32 populations of B. deitersi s. lat. Two nuclear and two mitochondrial loci were used to carry out the bGMYC, mPTP and STACEY algorithms for species delimitation. Detailed morphological study was also carried out across continents. The evidence indicated a widely distributed cryptic species in the Old World (Bosminopsis zernowi) that is genetically divergent from B. deitersi s.str. We revised the taxonomy and redescribed the species in this complex. Our sampling indicated that B. zernowi had weak genetic differentiation across its range. A molecular clock and biogeographic analysis with fossil calibrations suggested a Mesozoic origin for the Bosminopsis deitersi group. Our evidence rejects the single species hypothesis for B. deitersi and is consistent with an ancient species group (potentially Mesozoic) that shows marked morphological conservation. The family Bosminidae, then, has examples of both rapid morphological evolution (Holocene Bosmina), and morphological stasis (Bosminopsis).

Molecular biogeography of the Mediterranean Buthus species complex (Scorpiones: Buthidae) at its southern Palaearctic margin

Neogene orogenesis and climatic cycles strongly influenced inter- and intraspecific differentiation and variability of taxa. In this study, we focused on the southern margin of the western Palaearctic, known to be a geographically complex region. We performed mitochondrial DNA analyses of Buthus scorpions from the Moroccan Atlas Mountains, from the Hoggar Mountains in Algeria and from Tunisia, Iberia and Israel. Molecular species delimitation suggests the existence of ≥ 24 molecular operational taxonomic units. The data confirm complex differentiation patterns across the Atlas Mountains of Morocco, but structures in Iberia, Algeria and Tunisia have considerably lower complexity. This identifies the Atlas Mountain region as the most important differentiation centre of Buthus scorpions. Samples from the Hoggar Mountains (southern Algeria) cluster with those from the southernmost parts of Morocco in the middle and upper parts of the Draa Valley. This reinforces a recent connection of these regions. Samples from Israel are genetically similar to individuals from eastern Algeria and Tunisia. This suggests a widespread group across major parts of North Africa. Divergence time estimates indicate that differentiation in the genus began during the late Miocene, a period characterized by strong tectonic activities in this region. Further differentiation could be linked to subsequent climatic changes that have occurred since the end of the Miocene, with an increasing aridification of the Moroccan area. This also produced many microrefugia in the mountains of the area during the Pleistocene climatic fluctuations.

Multi-drug resistance of Staphylococcus aureus Strains in Baqiyatallah hospital: a Primary Step Towards Digital Health Biomonitoring Systems

The aim of the study was to evaluate the drug-resistance patterns of Staphylococcus aureus infections in Baqiyatallah hospital within 2010–2019 and to present a novel monitoring and detection system making use of molecular laboratory methods teamed with molecular delimitation analyses. This in turn is a primary step to establishment of a digital health system within Baqiyatallah hospital as a perfect pilot instance for other hospitals to follow upon. Totally, 100 patients of Baqiyatallah hospital suspicious of Staphylococcus aureus infections were sampled. Bacterial identity confirmations were done using routine biochemical test. Antibiograms were made for all the patients in this study. Consequently, bacterial total DNA was extracted and 16S rDNA gene amplified and sequenced for all patients. To uncover any cryptic strain grouping within the samples, a molecular delimitation method, i.e. automated barcode gap discovery (ABGD), was done. Our results showed Ceftaroline to be the most and Erythromycin and Oxacillin the least effective drugs. Delimitation uncovered 19 groups out of which group 19 seemed to have location-specific genetic signals in regards to susceptibility of Erythromycin and Oxacillin. Our results indicate the importance of genetic identification of bacteria with respect to their genetic patterns before antibiotic administration in order to both reduce unnecessary medicine use and to biomonitor the bacterial patterns in respect to their behavior towards general antibiotics.

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.

Quaternary climatic fluctuations influence the demographic history of two species of sky-island endemic amphibians in the Neotropics

We evaluated the role of Quaternary climatic fluctuations on the demographic history and population structure of amphibian species endemic to the 'campo rupestre' in the Neotropics, evaluating their distributional shifts, demographic changes, and lineage formation from the end of Pleistocene to present. We chose two anurans endemic to the high-elevation 'campo rupestre' in the Espinhaço Range (ER) in northeastern and southeastern Brazil (Bokermannohyla alvarengai and Bokermannohyla oxente), as models to test the role of Quaternary climatic fluctuations over their distribution range in this region. We collected tissue samples throughout their distribution range and used statistical phylogeography to examine processes of divergence and population demography. We generated spatial-temporal reconstructions using Bayesian inference in a coalescent framework in combination with hind-cast projections of species distribution models (SDMs). We also used the results and literature information to test alternative diversification scenarios via approximate Bayesian computation (ABC). Our results show that Quaternary climatic fluctuations influenced the geographic ranges of both species showing population expansion during the last glacial maximum (LGM) and range contraction during interglacial periods, as inferred from selected ABC models and from past projections of SDMs. We recovered Pleistocene diversification for both species occuring in distinctly unique periods for each taxon. An older and range-restricted lineage was recovered in a geographically isolated geological massif, deserving conservation and further taxonomic study. The diversification and distribution of these amphibian species endemic to the Neotropical 'campo rupestre' were influenced by Quaternary climatic fluctuations. The expansion of cold adapted species restricted to higher elevations during glacial periods and their concomitant retraction during interglacial periods may have been crucial for producing patterns of species richness and endemism along elevation gradients in tropical and subtropical domains. Such processes may influence the evolution of the biota distributed in heterogeneous landscapes with varied topography.

The number of species of degus (genus Octodon) is currently underestimated: An appraisal of species limits and their phylogenetic relationships (Rodentia: Hystricomorpha: Octodontidae)

As currently understood, the genus Octodon contains five species degus, lunatus, bridgesii, pacificus, and ricardojeda. Previous phylogenetic studies suggest that genus specific diversity is underestimated. In order to evaluate the taxonomic diversity of Octodon, we implemented unilocus (cytochrome-b) and multilocus (cytochrome-b + 4 nuclear genes) species delimitation methods. Octodon degus was recovered as a sister of the other species of the genus. The unilocus bGMYC and mPTP methods, based on cytochrome-b sequences, delimits 11 and 7 candidate species respectively, and both methods fail to recognize O. pacificus from O. ricardojeda. Results of the multilocus analysis (BPP) vary as a function of the dataset used. When the five genes are used 11 species are delimited, while eight species are delimited when only the nuclear genes are used. Octodon bridgesii is shown as comprising at least two species (one on the Pacific coast and the typical form found on the Andean slopes), while O. ricardojeda may comprise two species (one on the Chilean side of the Andes and the other in Argentina). Likewise, both multilocus matrices recover O. pacificus as a distinct species. This shows that species diversity of Octodon is underestimated. Remarkably, many of the delimited species based on genetic data are morphologically differentiated in cranio-dental characteristics. However, a pair of species has not achieved morphological differentiation, being cryptic species. Finally, the incongruence between mitochondrial and nuclear phylogenies suggests that processes such as incomplete lineage sorting and/or introgression have been present during the radiation of the genus.

Unexpected mitochondrial lineage diversity within the genus Alonella Sars, 1862 (Crustacea: Cladocera) across the Northern Hemisphere

Representatives of the genus Alonella Sars (Crustacea: Cladocera: Chydorinae) belong to the smallest known water fleas. Although species of Alonella are widely distributed and often abundant in acidic and mountain water bodies, their diversity is poorly studied. Morphological and genetic approaches have been complicated by the minute size of these microcrustaceans. As a result, taxonomists have avoided revising these species. Here, we present genetic data on Alonella species diversity across the Northern Hemisphere with particular attention to the A. excisa species complex. We analyzed 82 16S rRNA sequences (all newly obtained), and 78 COI sequences (39 were newly obtained). The results revealed at least twelve divergent phylogenetic lineages, possible cryptic species, of Alonella, with different distribution patterns. As expected, the potential species diversity of this genus is significantly higher than traditionally accepted. The A. excisa complex is represented by nine divergent clades in the Northern Hemisphere, some of them have relatively broad distribution ranges and others are more locally distributed. Our results provide a genetic background for subsequent morphological analyses, formal descriptions of Alonella species and detailed phylogeographical studies.

Cryptic Clitellata: Molecular Species Delimitation of Clitellate Worms (Annelida): An Overview

Methods for species delimitation using molecular data have developed greatly and have become a staple in systematic studies of clitellate worms. Here we give a historical overview of the data and methods used to delimit clitellates from the mid-1970s to today. We also discuss the taxonomical treatment of the cryptic species, including the recommendation that cryptic species, as far as possible, should be described and named. Finally, we discuss the prospects and further development of the field.

An integrative approach to address species limits in the southernmost members of the Liolaemus kingii group (Squamata: Liolaemini)

Recent conceptual and methodological advances have enabled an increasing number of studies to address the problem of species delimitation in a comprehensive manner. This is of particular interest in cases of species whose divergence times are recent and/or effective population sizes are large, where the conclusions obtained from a single source of evidence may lead to erroneous estimations of true species numbers or incorrect assignment of individuals to species. Iguanian lizards of the Liolaemus kingii group (13 species) comprise an important component of the endemic fauna of Patagonia. The southernmost species of this group (namely L. baguali, L. escarchadosi, L. sarmientoi, and L. tari) show widely overlapping distributions across southern Patagonia, also, their phylogenetic relationships are ambiguous and species boundaries have not been explicitly tested. Here we use a comprehensive approach to assess species limits through the use of molecular and morphological information (mitochondrial cytb, nuclear sequences collected by ddRADseq, and linear, meristic and landmark-based morphometrics). We found support for the current taxonomy given that the different analyses recognized the nominal species (4 entities), also a candidate species was supported by mitochondrial and morphological data. In addition, we detected signs of admixture between some of the species. Our results indicate that the L. kingii group can serve as a model system in studies of diversification accompanied by hybridization in nature, which in turn might have been promoted by past climatic oscillations and generalist morphologies. We emphasize the importance of using multiple lines of evidence in order to solve evolutionary stories, and minimizing potential erroneous results that may arise when relying on a single source of information.

ASAP: assemble species by automatic partitioning

Here, we describe Assemble Species by Automatic Partitioning (ASAP), a new method to build species partitions from single locus sequence alignments (i.e., barcode data sets). ASAP is efficient enough to split data sets as large 10⁴ sequences into putative species in several minutes. Although grounded in evolutionary theory, ASAP is the implementation of a hierarchical clustering algorithm that only uses pairwise genetic distances, avoiding the computational burden of phylogenetic reconstruction. Importantly, ASAP proposes species partitions ranked by a new scoring system that uses no biological prior insight of intraspecific diversity. ASAP is a stand-alone program that can be used either through a graphical web-interface or that can be downloaded and compiled for local usage. We have assessed its power along with three others programs (ABGD, PTP and GMYC) on 10 real COI barcode data sets representing various degrees of challenge (from small and easy cases to large and complicated data sets). We also used Monte-Carlo simulations of a multispecies coalescent framework to assess the strengths and weaknesses of ASAP and the other programs. Through these analyses, we demonstrate that ASAP has the potential to become a major tool for taxonomists as it proposes rapidly in a full graphical exploratory interface relevant species hypothesis as a first step of the integrative taxonomy process.

Miocene geologic dynamics of the Australian Sahul Shelf determined the biogeographic patterns of freshwater planorbid snails (Miratestinae) in the Indo-Australian Archipelago

The complex geological and climatic processes that have shaped the Indo-Australian Archipelago since the Cenozoic likely also gave rise to its species-rich biota. Strictly freshwater organisms might be particularly suitable for understanding the influence of these abiotic factors on their biogeography in such a insular setting as their distribution may reflect past abiotic events at large and small geographical scales. We here investigate the historical biogeography of the Miratestinae, a subfamily of Planorbidae. These freshwater gastropods are widely distributed in the eastern IAA from Australia, New Guinea, the Moluccas, and Sulawesi to the Philippines. The first comprehensive molecular phylogeny of the Miratestinae was inferred based on two mitochondrial and two nuclear genetic markers using maximum likelihood and Bayesian inference. Four species delimitation methods were applied to identify molecular operational taxonomic units (MOTUs). Divergence times were inferred using an uncorrelated lognormal relaxed-clock model by applying a taxon- and marker-specific substitution rate. Ancestral geographic ranges were estimated based on the dated phylogeny using BioGeoBEARS. The species delimitation revealed a total of 23 MOTUs, 16 of which might represent species new to science. The BioGeoBEARS analyses suggest an Australian origin for the Miratestinae at c. 22 Ma and identified jump dispersal to be the main process of colonization. The first colonization events from Australia to the IAA occurred in the Middle-Late Miocene (12-13 Ma), whereas intra-island diversification took mainly place since the Late Miocene-Pliocene. Colonization and diversification events remarkably coincide with major geologic events that shaped the geography of the region. The increasing availability of landmasses along the Sahul Shelf likely promoted stepping-stone dispersal to New Guinea, Sulawesi and the Philippines as early as the islands emerged. Major geological and climatic events such as the amalgamation of the island Sulawesi, the regional aridification in Australia or the uplift of massive mountain ranges in New Guinea likely played a considerable role for intra-island diversification.

Multiple species delimitation approaches applied to the avian lark genus Alaudala

Species delimitation has advanced from a purely phenotypic exercise to a branch of science that integrates multiple sources of data to identify independently evolving lineages that can be treated as species. We here test species limits in the avian Lesser Short-toed Lark Alaudala rufesens-Sand Lark A. raytal complex, which has an intricate taxonomic history, ranging from a single to three recognised species, with different inclusiveness in different treatments. Our integrative taxonomic approach is based on a combination of DNA sequences, plumage, biometrics, songs, song-flights, geographical distributions, habitat, and bioclimatic data, and using various methods including a species delimitation program (STACEY) based on the multispecies coalescent model. We propose that four species should be recognised: Lesser Short-toed Lark A. rufescens (sensu stricto), Heine's Short-toed Lark A. heinei, Asian Short-toed Lark A. cheleensis and Sand Lark A. raytal. There is also some evidence suggesting lineage separation within A. cheleensis and A. raytal, but additional data are required to evaluate this. The species delimitation based on STACEY agrees well with the non-genetic data. Although computer-based species delimitation programs can be useful in identifying independently evolving lineages, we stress that whenever possible, species hypotheses proposed by these programs should be tested by independent, non-genetic data. Our results highlight the difficulty and subjectivity of delimiting lineages and species, especially at early stages in the speciation process.

Cryptic Diversity in the Monotypic Neotropical Micromoth Genus Angelabella (Lepidoptera: Gracillariidae) in the Peru-Chile Desert

Simple Summary

The Neotropical Region harbors a highly diverse and poorly known fauna of leaf miners of the micromoth family Gracillariidae (Lepidoptera). Angelabella is a genus of Gracillariidae whose geographic range is restricted to a few valleys of the arid environments of the Peru-Chile desert. Only one species is currently included in this genus. The aims of this study were to explore the geographic range, determine the spatial distribution of mitochondrial lineages, and test lineage conspecificity hypotheses in Angelabella. The spatial distribution of genetic diversity indicated four spatial clusters, three of which are north of the previously known geographic range. These groups were defined as different species by four species delimitation methods. These results suggest that Angelabella harbors at least four morphologically cryptic species with restricted, not overlapping geographic ranges. This study shows that adequate single locus sequence analysis can be useful to discover surprising biodiversity patterns in underexplored environments, providing the base to plan further studies involving little-known organisms.

Abstract

Angelabella (Lepidoptera: Gracillariidae: Oecophyllembiinae) is considered a monotypic Neotropical genus of leaf miner micromoths known only from a few valleys of the arid environments of the Peru-Chile desert, particularly the southernmost part of Peru and northernmost part of Chile (type locality), where natural populations of its primary host plant occur. The geographic distribution of potential host plants provides a scenario for a wider range for this micromoth genus. The aims of this study were to explore the geographic range of Angelabella, determine the spatial distribution of mitochondrial lineages, and test lineage conspecificity hypotheses. The spatial distribution of genetic diversity indicated the presence of four spatial clusters, three of which are north of the previously known geographic range. Genetic distances were 0.2–0.8% and 3.6–8.3% (K2P) between haplotypes of the same and different spatial clusters, respectively. Phylogenetic relationships indicated reciprocal monophyly among the four spatial clusters, suggesting that allopatric differentiation processes have governed the recent history of Angelabella in these arid environments. These groups were defined as different species by four species delimitation methods, suggesting that Angelabella is not a monotypic genus, but harbors at least four morphologically cryptic allopatric species with restricted geographic ranges, including the type species and three candidate species.

Systematics, biogeography and evolution of the Saharo-Arabian naked-toed geckos genus Tropiocolotes

Plate tectonics constitute one of the main mechanisms of biological diversification on Earth, often being associated with cladogenetic events at different phylogenetic levels, as well as with exchange of faunas and floras across previously isolated biogeographic regions. North Africa and Arabia share a complex geological history that dates back to the break-up of the Arabian plate from the African plate ~30-25 Mya, followed by various geological events, such as the formation of the Red Sea or the connection between the African, Arabian and Eurasian plates. Species with Saharo-Arabian distributions have shown a close association between their evolutionary history and these geological events. In this study, we investigate the systematics, biogeography and evolution of the genus Tropiocolotes, a group of small ground-dwelling geckos, comprised by 12 species distributed from the Atlantic coast of North Africa to southwestern Iran. Species delimitation analyses uncovered the existence of high levels of undescribed diversity, with forms here considered at the species level including Tropiocolotes tripolitanus (Mauritania and southern Morocco), T. nattereri (southern Israel) and T. scorteccii (Yemen and Oman). Phylogenetic and biogeographic analyses recovered two main clades, an exclusively African clade and a Saharo-Arabian clade, that split ~25 Mya following the vicariant event mediated by the separation of the Arabian and African plates. The complex geological activity around the Red Sea is associated with the diversification within the Saharo-Arabian clade, including the colonization of North Africa from a second Tropiocolotes group. Results also provide new insights into the geographic distribution of Tropiocolotes nubicus, previously considered as exclusively associated to the Nile River valley, extending its known distribution further west, up to the Central Mountains of the Sahara. Accordingly, the Nile River seems to act as a major biogeographic barrier, separating Tropiocolotes nubicus and T. steudneri in their western and eastern margins, respectively.

Ultraconserved yet informative for species delimitation: Ultraconserved elements resolve long-standing systematic enigma in Central European bees

Accurate and testable species hypotheses are essential for measuring, surveying and managing biodiversity. Taxonomists often rely on mitochondrial DNA barcoding to complement morphological species delimitations. Although COI-barcoding has largely proven successful in assisting identifications for most animal taxa, there are nevertheless numerous cases where mitochondrial barcodes do not reflect species hypotheses. For instance, what is regarded as a single species can be associated with two distinct DNA barcodes, which can point either to cryptic diversity or to within-species mitochondrial divergences without reproductive isolation. In contrast, two or more species can share barcodes, for instance due to mitochondrial introgression. These intrinsic limitations of DNA barcoding are commonly addressed with nuclear genomic markers, which are expensive, may have low repeatability and often require high-quality DNA. To overcome these limitations, we examined the use of ultraconserved elements (UCEs) as a quick and robust genomic approach to address such problematic cases of species delimitation in bees. This genomic method was assessed using six different species complexes suspected to harbour cryptic diversity, mitochondrial introgression or mitochondrial paraphyly. The sequencing of UCEs recovered between 686 and 1,860 homologous nuclear loci and provided explicit species delimitation in all investigated species complexes. These results provide strong evidence for the suitability of UCEs as a fast method for species delimitation even in recently diverged lineages. Furthermore, we provide the first evidence for both mitochondrial introgression among distinct bee species, and mitochondrial paraphyly within a single bee species.