High level of cryptic species diversity revealed by sympatric lineages of Southeast Asian forest frogs

Morphological Species Delimitation in The Western Pond Turtle (Actinemys): Can Machine Learning Methods Aid in Cryptic Species Identification?

As the discovery of cryptic species has increased in frequency, there has been an interest in whether geometric morphometric data can detect fine-scale patterns of variation that can be used to morphologically diagnose such species. We used a combination of geometric morphometric data and an ensemble of five supervised machine learning methods (MLMs) to investigate whether plastron shape can differentiate two putative cryptic turtle species, Actinemys marmorata and Actinemys pallida. Actinemys has been the focus of considerable research due to its biogeographic distribution and conservation status. Despite this work, reliable morphological diagnoses for its two species are still lacking. We validated our approach on two datasets, one consisting of eight morphologically disparate emydid species, the other consisting of two subspecies of Trachemys (T. scripta scripta, T. scripta elegans). The validation tests returned near-perfect classification rates, demonstrating that plastron shape is an effective means for distinguishing taxonomic groups of emydids via MLMs. In contrast, the same methods did not return high classification rates for a set of alternative phylogeographic and morphological binning schemes in Actinemys. All classification hypotheses performed poorly relative to the validation datasets and no single hypothesis was unequivocally supported for Actinemys. Two hypotheses had machine learning performance that was marginally better than our remaining hypotheses. In both cases, those hypotheses favored a two-species split between A. marmorata and A. pallida specimens, lending tentative morphological support to the hypothesis of two Actinemys species. However, the machine learning results also underscore that Actinemys as a whole has lower levels of plastral variation than other turtles within Emydidae, but the reason for this morphological conservatism is unclear.

A new odorous frog species of Odorrana (Amphibia, Anura, Ranidae) from Guizhou Province, China

The frog genus Odorrana is distributed across east and southeastern Asia. Based on morphological differences and molecular phylogenetics, a new species of the genus occurring from Leigong Mountain in Guizhou Province, China is described. Phylogenetic analyses based on DNA sequences of the mitochondrial 12S rRNA, 16S rRNA, and ND2 genes supported the new species as an independent lineage. The uncorrected genetic distances between the 12S rRNA, 16S rRNA, and ND2 genes between the new species and its closest congener were 5.0%, 4.9%, and 16.3%, respectively. The new species is distinguished from its congeners by a combination of the following characters: body size moderate (SVL 39.1-49.4 mm in males, 49.7 mm in female); head width larger than head length; tympanum distinctly visible; small rounded granules scattered all over dorsal body and limbs; dorsolateral folds absent; heels overlapping when thighs are positioned at right angles to the body; tibiotarsal articulation reaching the level between eye to nostril when leg stretched forward; vocal sacs absent in male and nuptial pads present on the base of finger I.

Amphibian diversity across three adjacent ecosystems in Área de Conservación Guanacaste, Costa Rica

Amphibians are the most threatened species-rich vertebrate group, with species extinctions and population declines occurring globally, even in protected and seemingly pristine habitats. These 'enigmatic declines' are generated by climate change and infectious diseases. However, the consequences of these declines are undocumented as no baseline ecological data exists for most affected areas. Like other neotropical countries, Costa Rica, including Área de Conservación Guanacaste (ACG) in north-western Costa Rica, experienced rapid amphibian population declines and apparent extinctions during the past three decades. To delineate amphibian diversity patterns within ACG, a large-scale comparison of multiple sites and habitats was conducted. Distance and time constrained visual encounter surveys characterised species richness at five sites-Murciélago (dry forest), Santa Rosa (dry forest), Maritza (mid-elevation dry-rain forest intersect), San Gerardo (rainforest) and Cacao (cloud forest). Furthermore, species-richness patterns for Cacao were compared with historic data from 1987-8, before amphibians declined in the area. Rainforests had the highest species richness, with triple the species of their dry forest counterparts. A decline of 45% (20 to 11 species) in amphibian species richness was encountered when comparing historic and contemporary data for Cacao. Conservation efforts sometimes focus on increasing the resilience of protected areas, by increasing their range of ecosystems. In this sense ACG is unique containing many tropical ecosystems compressed in a small geographic space, all protected and recognised as a UNESCO world heritage site. It thus provides an extraordinary platform to understand changes, past and present, and the resilience of tropical ecosystems and assemblages, or lack thereof, to climate change.

Distinguishing fanged frogs (Limnonectes) species (Amphibia: Anura: Dicroglossidae), from Thailand using high resolution melting analysis

Morphologically, species of fanged frogs (Limnonectes) are exceedingly similar, making it difficult to distinguish them within the complex. In Thailand, it has been difficult to distinguish between the sympatric species L. bannaensis and L. taylori, particularly among tadpoles, adolescents, and adult females. A precise identification contributes to a greater understanding of biodiversity, particularly for assessing distributions and population dynamics. Therefore, a novel approach is required. The objective of this study was to develop a high resolution melting analysis (HRM) for the rapid and accurate identification of six species of Limnonectes of the L. kuhlii complex found in Thailand, particularly the two sympatric fanged frogs. Here, HRM assays using 16S rRNA mitochondrial primers were designed and developed. There was as much as a 25.3% variation in the nucleotide sequence of the fragment amplified by HRM16S primers among the six species of Limnonectes. Prior to conducting an in vitro HRM, the DNA sequences were used in a simulation HRM, uMELT Quartz, to predict the melting curve for each species of Limnonectes. There were discrepancies between the predicted melting curves of each species generated by the programme. Consequently, in vitro HRM tests were conducted. The obtained melting curve and Tm values were consistent with those predicted, albeit with a slightly different Tm value and a more distinct melting curve. All evaluated species of Limnonectes could be easily distinguished from one another by comparing the melting curve shapes. The HRM assay was then used to confirm the species of 18 Limnonectes samples in comparison to the reference samples (confidence interval > 90%). In addition, the results of HRM were consistent with those of experts who used morphological analysis to identify species. The HRM was found to be useful, and therefore the method would also contribute to future ecological and systematic studies on the target species.

Species Delimitation, Phylogenomics, and Biogeography of Sulawesi Flying Lizards: A Diversification History Complicated by Ancient Hybridization, Cryptic Species, and Arrested Speciation

The biota of Sulawesi is noted for its high degree of endemism and for its substantial levels of in situ biological diversification. While the island's long period of isolation and dynamic tectonic history have been implicated as drivers of the regional diversification, this has rarely been tested in the context of an explicit geological framework. Here, we provide a tectonically informed biogeographical framework that we use to explore the diversification history of Sulawesi flying lizards (the Draco lineatus Group), a radiation that is endemic to Sulawesi and its surrounding islands. We employ a framework for inferring cryptic speciation that involves phylogeographic and genetic clustering analyses as a means of identifying potential species followed by population demographic assessment of divergence-timing and rates of bi-directional migration as means of confirming lineage independence (and thus species status). Using this approach, phylogenetic and population genetic analyses of mitochondrial sequence data obtained for 613 samples, a 50-SNP data set for 370 samples, and a 1249-locus exon-capture data set for 106 samples indicate that the current taxonomy substantially understates the true number of Sulawesi Draco species, that both cryptic and arrested speciations have taken place, and that ancient hybridization confounds phylogenetic analyses that do not explicitly account for reticulation. The Draco lineatus Group appears to comprise 15 species-9 on Sulawesi proper and 6 on peripheral islands. The common ancestor of this group colonized Sulawesi ~11 Ma when proto-Sulawesi was likely composed of two ancestral islands, and began to radiate ~6 Ma as new islands formed and were colonized via overwater dispersal. The enlargement and amalgamation of many of these proto-islands into modern Sulawesi, especially during the past 3 Ma, set in motion dynamic species interactions as once-isolated lineages came into secondary contact, some of which resulted in lineage merger, and others surviving to the present. [Genomics; Indonesia; introgression; mitochondria; phylogenetics; phylogeography; population genetics; reptiles.].

The identity of Nanorana yunnanensis (Anderson, 1879) from Thailand, with a description of its tadpole (Anura, Dicroglossidae)

Nanorana yunnanensis has been reported in Thailand for many years; however, no specimens have been verifiably described. To verify the identity of this species in Thailand, we observed Thai specimens formerly assigned to the species N. bourreti and N. yunnanensis. The 12S rRNA, 16S rRNA, and ND2 mitochondrial genes were used to clarify the taxonomic statuses of the Thai specimens. The results showed that our specimens formed a clade with N. yunnanensis, containing specimens from southern Yunnan and Vietnam. The respective mean sequence divergences of the 12S and 16S rRNA genes between populations from Thailand versus China were 0.59% and 1.09%, for populations from Thailand versus Vietnam were 0.57% and 0.61%, and for populations from China versus Vietnam were 0.59% and 1.09%. We provide descriptions of the tadpole and adult specimens of N. yunnanensis from Thailand.

Molecular phylogeny reveals distinct evolutionary lineages of the banded krait, Bungarus fasciatus (Squamata, Elapidae) in Asia

The banded krait, Bungarus fasciatus is a widespread elapid snake, likely to comprise several distinct species in different geographic regions of Asia. Therefore, based on molecular phylogenetics and comparative morphology data, we present an overview of the systematic composition of the species to delimit potential biogeographic boundaries. Our phylogenetic analyses, based on four mitochondrial genes, reveal the existence of at least three evolutionary lineages within B. fasciatus, corresponding to Indo-Myanmar, Sundaic and eastern Asian lineages. We are convinced that there are at least three taxonomic entities within the nomen B. fasciatus and restrict the distribution of B. fasciatus sensu stricto to the Indo-Myanmar region. We also provide additional natural history data of the taxon from eastern India. Finally, we advocate further studies to establish the degree of reproductive isolation among these diverging evolutionary lineages and to reassess the systematic status of this species complex especially the Sundaic and eastern Asian lineages.

Ventral and lateral spot patterns differentiation between three smooth newt species (Amphibia: Salamandridae: Lissotriton)

The location, shape and area of dark spots on the belly, throat and sides of the body of three closely related species of smooth newts (Lissotriton kosswigi, L. schmidtleri and L. vulgaris) were analysed. Differences were found between species and sexes in the location and shape of dark spots on the belly, throat and lateral surfaces of the body. Discriminant analysis of these spot characteristics allows to identificate of males of all three species with a high degree of confidence (85–91%). In females, only L. vulgaris was correctly distinguished from the other two species (accuracy 81–94%). Anatolian and Thracian populations of L. schmidtleri have very similar patterns of dark spots, which confirm their conspecificity. The differences in the location and size of dark spots make possible successful identification of the species during field research and study of museum specimens. The developed method could be useful for comparative studies of other animal species which have spotted camouflage or aposematic colouration.

Phylogeographic inference of Sumatran ranids bearing gastromyzophorous tadpoles with regard to the Pleistocene drainage systems of Sundaland

Rivers are known to act as biogeographic barriers in several strictly terrestrial taxa, while possibly serving as conduits of dispersal for freshwater-tolerant or -dependent species. However, the influence of river systems on genetic diversity depends on taxa-specific life history traits as well as other geographic factors. In amphibians, several studies have demonstrated that river systems have only minor influence on their divergence. Here, we assess the role of the paleodrainage systems of the Sunda region (with a focus on the island of Sumatra) in shaping the evolutionary history of two genera of frogs (Sumaterana and Wijayarana) whose tadpoles are highly dependent on cascading stream habitats. Our phylogenetic results show no clear association between the genetic diversification patterns of both anurans genera and the existence of paleodrainage systems. Time-calibrated phylogenies and biogeographical models suggest that these frogs colonized Sumatra and diversified on the island before the occurrence of the Pleistocene drainage systems. Both genera demonstrate phylogenetic structuring along a north–south geographic axis, the temporal dynamics of which coincide with the geological chronology of proto Sumatran and -Javan volcanic islands. Our results also highlight the chronic underestimation of Sumatran biodiversity and call for more intense sampling efforts on the island.

Molecular data confirm Triatoma pallidipennis Stål, 1872 (Hemiptera: Reduviidae: Triatominae) as a novel cryptic species complex

Triatoma pallidipennis constitues one of the most important Chagas disease vector in Mexico. Previous studies based on molecular data suggest T. pallidipennis as a complex of cryptic species. For that reason, we analyzed the phylogenetic relationships of T. pallidipennis using DNA sequences from the mitochondrial ND4 gene and the ITS-2 gene. In addition, the divergence times were estimated, and possible new taxa were delimited with three species delimitation methods. Finally, genetic distances and possible connectivity routes based on shared haplotypes were obtained among the T. pallidipennis populations. Five haplogroups (possible cryptic species) were found, based on delimitation methods and genetic distances. Haplogroup divergence began about 3 Ma, in the Pleistocene. Moreover, none of the haplogroups showed potential connectivity routes between them, evidencing lack of gene flow. Our results suggest the existence of a new cryptic species complex within what is currently recognized as a T. pallidipennis.

How to Accurately Delineate Morphologically Conserved Taxa and Diagnose Their Phenotypic Disparities: Species Delimitation in Cryptic Rhinolophidae (Chiroptera)

Systematics and taxonomy are the backbone of all components of biology and ecology, yet cryptic species present a major challenge for accurate species identification. This is especially problematic as they represent a substantial portion of undiscovered biodiversity, and have implications for not only species conservation, but even assaying potential risk of zoonotic spillover. Here, we use integrative approaches to delineate potential cryptic species in horseshoe bats (Rhinolophidae), evaluate the phenotypic disparities between cryptic species, and identify key traits for their identification. We tested the use of multispecies coalescent models (MSC) using Bayesian Phylogenetic and Phylogeography (BPP) and found that BPP was useful in delineating potential cryptic species, and consistent with acoustic traits. Our results show that around 40% of Asian rhinolophid species are potentially cryptic and have not been formally described. In order to avoid potential misidentification and allow species to be accurately identified, we identified quantitative noseleaf sella and acoustic characters as the most informative traits in delineating between potential cryptic species in Rhinolophidae. This highlights the physical differences between cryptic species that are apparent in noseleaf traits which often only qualitatively described but rarely measured. Each part of the noseleaf including the sella, lateral lappets, and lancet furrows, play roles in focusing acoustic beams and thus, provide useful characteristics to identify cryptic Rhinolophus species. Finally, species delimitation for cryptic species cannot rely on genetic data alone, but such data should be complemented by other evidence, including phenotypic, acoustic data, and geographic distributions to ensure accurate species identification and delineation.

DNA barcode reveals candidate species of Scinax and Ololygon (Anura: Hylidae) in Atlantic Forest

Molecular species delimitation methods are efficient tools to identify species, including the discovery of new taxa and cryptic organisms, thus being useful to biodiversity studies. In the present work, 16S mitochondrial sequences and cytochrome oxidase I (COI) were used to evaluate the richness of species in the genus Scinax and Ololygon from a biodiversity hotspot in Atlantic Forest. A total of 109 specimens formally belonging to eight species of Scinax and three species of Ololygon were collected in 13 localities along the state of Bahia (northeastern Brazil) and one site in Espírito Santo (southeastern Brazil). Of the Scinax species collected in this study, three were morphologically differentiated from other described species and identified as putative new species (Scinax sp.1, Scinax sp.2 and Scinax sp.3). The species delimitations were inferred using three different methods: ABGD, PTP and mPTP which allowed recognizing 11 Scinax species and five Ololygon species. Scinax sp. 1, Scinax sp. 2 and Scinax sp. 3, have been confirmed as new putative species and Ololygon argyreornata possibly contains cryptic species. We suggest additional studies, including morphological and bioacoustic data to validate these new putative species.

An insight into molecular taxonomy of bufonids, microhylids, and dicroglossid frogs: First genetic records from Pakistan

The current study was focused on documentation of amphibian assemblage in North Punjab and Islamabad Capital Territory, Pakistan, by using mitochondrial gene sequences of 16S rRNA. Our study entailed 37% of the known amphibian species of the country. We provided a phylogenetic analysis based on 74 newly generated mitochondrial 16S rRNAs from nine species of genus Microlyla, Duttaphrynus, Allopaa, Nanorana, Sphaerotheca, Minervarya, Hoplobatrachus, and Euphlyctis. We employed the maximum-likelihood inference and Bayesian analysis to assess the taxonomic status of the samples obtained from Pakistan, with respect to other congeneric species from surrounding regions. Our findings confirmed the taxonomic status of South Asian anuran species Duttaphrynus stomaticus, Duttaphrynus melanostictus, Microhyla nilphamariensis, Allopaa hazarensis, Nanorana vicina, Sphaerotheca maskeyi (synonym: S. pashchima), Minervarya pierrei, Hoplobatrachus tigerinus, and Euphlyctis kalasgramensis in Pakistan. We have reported new country records of genus Minervarya ( M. pierrei). Minervarya pierrei was previously misidentified as Fejervarya limnocharis, due to dearth of genetic information. We provided the first genetic records of our endemic species N. vicina. The results revealed the taxonomic placement of N. vicina with respect to its congeners and validated the taxonomic status of N. vicina from its type locality (Murree) for the first time. The findings of the present study also indicated the paraphyletic relationship of A.- hazarensis with Nanorana species. So, based on our phylogenetic inferences, morphological characters, and habitat preferences, validity of generic status of A. hazarensis is undecided. As our data were not enough to resolve this issue, we suggest sequencing of additional mitochondrial and nuclear genes in the future studies to get a better resolution. We recommend carrying out extensive surveys throughout the country for proper scientific documentation of amphibians of Pakistan. Many new species, some of them might be endemic to Pakistan, are expected to be discovered, and taxonomic status of other species would be resolved.

King or royal family? Testing for species boundaries in the King Cobra, Ophiophagus hannah (Cantor, 1836), using morphology and multilocus DNA analyses

In widespread species, the diverse ecological conditions in which the populations occur, and the presence of many potential geographical barriers through their range are expected to have created ample opportunities for the evolution of distinct, often cryptic lineages. In this work, we tested for species boundaries in one such widespread species, the king cobra, Ophiophagus hannah (Cantor, 1836), a largely tropical elapid snake distributed across the Oriental realm. Based on extensive geographical sampling across most of the range of the species, we initially tested for candidate species (CS) using Maximum-Likelihood analysis of mitochondrial genes. We then tested the resulting CS using both morphological data and sequences of three single-copy nuclear genes. We used snapclust to determine the optimal number of clusters in the nuclear dataset, and Bayesian Phylogenetics and Phylogeography (BPP) to test for likely species status. We used non-metric multidimensional scaling (nMDS) analysis for discerning morphological separation. We recovered four independently evolving, geographically separated lineages that we consider Confirmed Candidate Species: (1) Western Ghats lineage; (2) Indo-Chinese lineage (3) Indo-Malayan lineage; (4) Luzon Island lineage, in the Philippine Archipelago. We discuss patterns of lineage divergence, particularly in the context of low morphological divergence, and the conservation implications of recognizing several endemic king cobra lineages.

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.

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.

Genomic divergence in sympatry indicates strong reproductive barriers and cryptic species within Eucalyptus salubris

Genetic studies are increasingly detecting cryptic taxa that likely represent a significant component of global biodiversity. However, cryptic taxa are often criticized because they are typically detected serendipitously and may not receive the follow-up study required to verify their geographic or evolutionary limits. Here, we follow-up a study of Eucalyptus salubris that unexpectedly detected two divergent lineages but was not sampled sufficiently to make clear interpretations. We undertook comprehensive sampling for an independent genomic analysis (3,605 SNPs) to investigate whether the two purported lineages remain discrete genetic entities or if they intergrade throughout the species' range. We also assessed morphological and ecological traits, and sequenced chloroplast DNA. SNP results showed strong genome-wide divergence (F ST = 0.252) between two discrete lineages: one dominated the north and one the southern regions of the species' range. Within lineages, gene flow was high, with low differentiation (mean F ST = 0.056) spanning hundreds of kilometers. In the central region, the lineages were interspersed but maintained their genomic distinctiveness: an indirect demonstration of reproductive isolation. Populations of the southern lineage exhibited significantly lower specific leaf area and occurred on soils with lower phosphorus relative to the northern lineage. Finally, two major chloroplast haplotypes were associated with each lineage but were shared between lineages in the central distribution. Together, these results suggest that these lineages have non-contemporary origins and that ecotypic adaptive processes strengthened their divergence more recently. We conclude that these lineages warrant taxonomic recognition as separate species and provide fascinating insight into eucalypt speciation.

Revision shock in Pacific oysters taxonomy: the genus Magallana (formerly Crassostrea in part) is well-founded and necessary

The description of the genus Magallana provoked taxonomic ‘revision shock’ among Bayne and colleagues, who criticised: (1) the supporting evidence for this taxonomic change, (2) the procedure used and (3) its impact on the nomenclatural stability of oysters. Here, we demonstrate that the description of the genus Magallana fulfils the norms of taxonomic and indexing revisions, and is well-founded on a scientific basis. The clade named Magallana is supported by a robust phylogeny based on comprehensive taxon sampling, independent datasets and varied analytical methods. A new maximum likelihood phylogenetic analysis of oysters, based on mitochondrial gene order data (representing molecular phenotypes above the sequence level), has provided additional support for this classification. These results are fully consistent with all previously published phylogenetic studies, thus providing an unambiguous indication of the stability of the clade that meets all the currently accepted criteria for naming clades as taxa. We show that, while the criterion of ‘morphological diagnosability’ is not applicable at any level of oyster classification, several striking molecular phenotypes are diagnostic of Magallana, both at the nuclear and the mitochondrial genome level. The classification with three genera Magallana, Crassostrea and Talonostrea reflects the evolutionary diversity of Crassostreinae and is consistent with taxonomic ranking criteria adopted for other oyster subfamilies.

Does evolution of echolocation calls and morphology in Molossus result from convergence or stasis?

Although many processes of diversification have been described to explain variation of morphological traits within clades that have obvious differentiation among taxa, not much is known about these patterns in complexes of cryptic species. Molossus is a genus of bats that is mainly Neotropical, occurring from the southeastern United States to southern Argentina, including the Caribbean islands. Molossus comprises some groups of species that are morphologically similar but phylogenetically divergent, and other groups of species that are genetically similar but morphologically distinct. This contrast allows investigation of unequal trait diversification and the evolution of morphological and behavioural characters. In this study, we assessed the role of phylogenetic history in a genus of bat with three cryptic species complexes, and evaluated if morphology and behavior are evolving concertedly. The Genotype by Sequence genomic approach was used to build a species-level phylogenetic tree for Molossus and to estimate the ancestral states of morphological and echolocation call characters. We measured the correlation of phylogenetic distances to morphological and echolocation distances, and tested the relationship between morphology and behavior when the effect of phylogeny is removed. Morphology evolved via a mosaic of convergence and stasis, whereas call design was influenced exclusively through local adaptation and convergent evolution. Furthermore, the frequency of echolocation calls is negatively correlated with the size of the bat, but other characters do not seem to be evolving in concert. We hypothesize that slight variation in both morphology and behaviour among species of the genus might result from niche specialization, and that traits evolve to avoid competition for resources in similar environments.

Gene flow creates a mirage of cryptic species in a Southeast Asian spotted stream frog complex

Most new cryptic species are described using conventional tree- and distance-based species delimitation methods (SDMs), which rely on phylogenetic arrangements and measures of genetic divergence. However, although numerous factors such as population structure and gene flow are known to confound phylogenetic inference and species delimitation, the influence of these processes is not frequently evaluated. Using large numbers of exons, introns, and ultraconserved elements obtained using the FrogCap sequence-capture protocol, we compared conventional SDMs with more robust genomic analyses that assess population structure and gene flow to characterize species boundaries in a Southeast Asian frog complex (Pulchrana picturata). Our results showed that gene flow and introgression can produce phylogenetic patterns and levels of divergence that resemble distinct species (up to 10% divergence in mitochondrial DNA). Hybrid populations were inferred as independent (singleton) clades that were highly divergent from adjacent populations (7%-10%) and unusually similar (<3%) to allopatric populations. Such anomalous patterns are not uncommon in Southeast Asian amphibians, which brings into question whether the high levels of cryptic diversity observed in other amphibian groups reflect distinct cryptic species-or, instead, highly admixed and structured metapopulation lineages. Our results also provide an alternative explanation to the conundrum of divergent (sometimes nonsister) sympatric lineages-a pattern that has been celebrated as indicative of true cryptic speciation. Based on these findings, we recommend that species delimitation of continuously distributed "cryptic" groups should not rely solely on conventional SDMs, but should necessarily examine population structure and gene flow to avoid taxonomic inflation.

Determining levels of cryptic diversity within the endemic frog genera, Indirana and Walkerana, of the Western Ghats, India

A large number of species in the tropics are awaiting discovery, many due to their cryptic morphology ie. lack of discernable morphological difference. We explored the presence of cryptic lineages within the frog genera, Indirana and Walkerana, which are endemic to the Western Ghats of Peninsular India. By reconstructing a phylogeny using 5 genes and robust geographic sampling, we delimited 19 lineages along a population—species continuum, using multiple criteria including haplotype clusters, genetic distance, morphological distinctness, and geographical separation. Of these 19 lineages, 14 belonged to the genus Indirana and 5 to the genus Walkerana. Divergence dating analyses revealed that the clade comprising Indirana and Walkerana began diversifying around 71 mya and the most recent common ancestor of Indirana and Walkerana split around 43 mya. We tested for the presence of cryptic lineages by examining the relationship between genetic and morphological divergence among related pairs within a pool of 15 lineages. The pairs showed strong morphological conservatism across varying levels of genetic divergence. Our results highlight the prevalence of morphologically cryptic lineages in these ancient endemic clades of the Western Ghats. This emphasizes the significance of other axes, such as geography, in species delimitation. With increasing threats to amphibian habitats, it is imperative that cryptic lineages are identified so that appropriate conservation measures can be implemented.

Impact of the Mid-Pleistocene Revolution and Anthropogenic Factors on the Dispersion of Asian Black-Spined Toads (Duttaphrynus melanostictus)

Divergence-time estimation critically improves the understanding of biogeography processes underlying the distribution of species, especially when fossil data is not available. We hypothesise that the Asian black-spined toad, Duttaphrynus melanostictus, expanded into the Eastern Indomalaya following the Quaternary glaciations with the subsequent colonisation of new landscapes during the Last Glacial Maximum. Divergence dating inferred from 364 sequences of mitochondrial tRNAGly ND3 supported the emergence of a common ancestor to the three D. melanostictus clades around 1.85 (±0.77) Ma, matching with the Lower to Mid-Pleistocene transition. Duttaphrynus melanostictus then dispersed into Southeast Asia from the central Indo-Pacific and became isolated in the Southern Sundaic and Wallacea regions 1.43 (±0.10) Ma through vicariance as a result of sea level oscillations. The clade on the Southeast Asian mainland then colonised the peninsula from Myanmar to Vietnam and expanded towards Southeastern China at the end of the Mid-Pleistocene Revolution 0.84 (±0.32) Ma. Population dynamics further highlight an expansion of the Southeast Asian mainland population towards Taiwan, the Northeastern edge of the species' range after the last interglacial, and during the emergence of the Holocene human settlements around 7000 BP. Thus, the current divergence of D. melanostictus into three segregated clades was mostly shaped by Quaternary glaciations, followed by natural dispersion events over land bridges and accelerated by anthropogenic activities.

Integrative taxonomy reveals three new taxa within the Tylototriton asperrimus complex (Caudata, Salamandridae) from Vietnam

The Tylototritonasperrimus complex from northern Vietnam is reviewed based on morphological comparisons and analysis of the mitochondrial marker NADH dehydrogenase subunit 2 (ND2). Based on molecular divergences, which were revealed to be higher than in other congeners, in concert with morphological differences, two new species and one subspecies are described herein: Tylototritonpasmansisp. nov. differs from T.asperrimussensu stricto by 3.2 to 3.6 % genetic divergence and a combination of distinct morphological characters, such as head slightly longer than wide, distinct mid-dorsal ridge, relatively wide distance between the eyes, tips of fingers reaching the eye when foreleg is laid forward, labial and gular folds present, central belly skin with tubercles shaped like transverse wrinkles and distinct, pointy to round rib nodules. The population of T.pasmansisp. nov. consists of two subclades, the nominotypic one occurring on the eastern side of the Da River (or Black River, including Hoa Binh and Phu Tho provinces), and another occurring on the western side (including Son La and Thanh Hoa provinces). These two subclades differ by 2.5 to 3.1 % genetic divergence and distinct morphological characters. The western subclade is herein described as Tylototritonpasmansiobstissp. nov., which differs from the nominotypic form by a wider head, longer and narrower snout, shorter femur length, and an overall less granulose skin, without an increased concentration of warts on the body sides.

A second new species, Tylototritonsparreboomisp. nov. is described from Lai Chau Province. It differs from T.asperrimussensu stricto by 4.1 to 4.2 % and from Tylototritonpasmansisp. nov. by 3.6 to 4.5 % genetic divergences as well as by a combination of distinct morphological characters, such as head longer than wide, tips of fingers reaching nostril when foreleg adpressed along head, rib nodules distinct, round and relatively enlarged, and wide distance between the eyes.

Four new species of Pristimantis Jiménez de la Espada, 1870 (Anura: Craugastoridae) in the eastern Amazon

The Pristimantis genus (Anura: Craugastoridae) is the most diverse among all vertebrates with 531 described species. The highest diversity occurs in Ecuador (215 species), followed by Colombia (202), Peru (139), Venezuela (60), Brazil (30), Bolivia (17), Guyana (6) Suriname and French Guiana (5). The genus is divided into 11 species groups. Of these, the P. conspicillatus group (containing 34 species), distributed in extreme southeastern Costa Rica, Isla Taboga (Panama), northern South America (from Colombia to eastern Guyana), south Bolivia, and is the best represented in Brazil (16 species). The main characteristics of this group are the tympanic membrane and tympanic annulus distinct (except in P. johannesdei); dorsum smooth or shagreen; dorsal lateral fold present or absent; usually smooth belly, but may be weakly granular in some species; toe V slightly larger than the toe III. Most of the taxonomic inconsistencies in species of Pristimantis could be due to its much conserved morphology and the lack of comprehensive taxonomic evaluations. Thus, an ongoing challenge for taxonomists dealing with the Pristimantis genus is the ubiquitous abundance of cryptic species. In this context, accurate species delimitation should integrate evidences of morphological, molecular, bioacoustics and ecological data, among others. Based on an integrative taxonomy perspective, we utilize morphological, molecular (mtDNA) and bioacoustic evidence to describe four new species of the Pristimantis conspicillatus group from the eastern Amazon basin. Pristimantis giorgiisp. nov. is known from the Xingu/Tocantins interfluve and can be distinguished from the other Pristimantis species of the region by presenting discoidal fold, dorsolateral fold absent, vocalization composed of three to four notes and genetic distance of 7.7% (16S) and 14.8% (COI) from P. latro, the sister and sympatric species with respect P. giorgiisp. nov.. Pristimantis pictussp. nov. is known to the northern Mato Grosso state, Brazil, and can be distinguished from the other species of Pristimantis by presenting the posterior surface of the thigh with light yellow patches on a brown background, also extending to the inguinal region, vocalization consisting of four to five notes and a genetic distance of 11.6% (16S) and 19.7% (COI) from P. pluviansp. nov., which occurs in sympatry. Pristimantis pluviansp. nov. is known to the northern Mato Grosso state, Brazil, and may be distinguished from the other Pristimantis species by having a posterior surface of the thigh reddish and vocalization composed of two notes. Pristimantis moasp. nov. is known to the northern Tocantins state and southwestern Maranhão state. This species can be distinguished from the other Pristimantis species by possessing slightly perceptible canthal stripe, external thigh surface with dark yellow spots on brown background, vocalization consisting of three to five notes and genetic distance of 2.3–11.7 (16S) and 10.5–23.1 (COI) for the new Pristimantis species of this study.

Population genetics and sympatric divergence of the freshwater gudgeon, Gobiobotia filifer, in the Yangtze River inferred from mitochondrial DNA

The ecosystem and Pleistocene glaciations play important roles in population demography. The freshwater gudgeon, Gobiobotia filifer, is an endemic benthic fish in the Yangtze River and is a good model for ecological and evolutionary studies. This study aimed to decode the population structure of G. filifer in the Yangtze River and reveal whether divergence occurred before or after population radiation. A total of 292 specimens from eight locations in the upper and middle reaches of the Yangtze River were collected from 2014 to 2016 and analyzed via mitochondrial DNA Cyt b gene sequencing. A moderately high level of genetic diversity was found without structures among the population. However, phylogenetic and network topology showed two distinct haplotype groups, and each group contained a similar proportion of individuals from all sampled sites. This suggested the existence of two genetically divergent source populations in G. filifer. We deduced that a secondary contact of distinct glacial refugia was the main factor creating sympatric populations of G. filifer, and climate improvement promoted population expansion and colonization.

Species delimitation of crab-eating frogs (Fejervarya cancrivora complex) clarifies taxonomy and geographic distributions in mainland Southeast Asia

The taxonomy and geographic distributions of species of crab-eating frogs (Fejervaryacancrivora complex) in mainland Southeast Asia have been highly uncertain. Three taxonomic names are used in recent literature (F.cancrivora, F.raja, and F.moodiei) but the applications of these names to localities has been inconsistent, especially owing to the lack of available molecular data for F.raja. Morphometric and mitochondrial DNA variation was examined in these frogs, including name-bearing types and topotypes of all three species. Findings corroborate evidence for the existence of two species in coastal mainland Southeast Asia, with F.moodiei having a wide geographic distribution and F.cancrivora sensu stricto occurring only in extreme southern Thailand and peninsular Malaysia. Fejervaryaraja is shown to be only a large-bodied population of F.cancrivora sensu stricto and is synonymized with that species. Revised descriptions of F.moodiei and F.cancrivora sensu stricto are provided.

Cryptic diversity of limestone karst inhabiting land snails (Cyclophorus spp.) in northern Vietnam, their evolutionary history and the description of four new species

Limestone karsts can form terrestrial habitat islands for calcium-dependent organisms. In Vietnam, many karst habitats are threatened, while their rich biodiversity is still far from being thoroughly explored. Given that conservation of karst biota strongly relies on correct species identification, the presence of undetected cryptic species can pose severe problems. The present study focuses on cryptic diversity among karst-inhabiting land snails of the genus Cyclophorus in northern Vietnam, where specimens with a similar shell morphology have been reported from various regions. In order to examine the diversity and evolutionary history of this “widespread morphotype”, we generated a Bayesian phylogeny based on DNA sequence data. Automatic Barcode Gap Discovery (ABGD) and the Bayesian implementation of the Poisson tree processes model (bPTP) contributed to species delimitation and analyses of shell shape and size aided the morphological characterisation of individual species. We found that the examined specimens of the widespread morphotype did not form a single monophyletic group in the phylogeny but clustered into several different clades. We delimited nine different species that develop the widespread morphotype and described four of them as new. Processes of convergent evolution were probably involved in the origin of the delimited species, while their generally allopatric distribution could result from interspecific competition. Our findings indicate ongoing processes of speciation and a potential case of morphological character displacement. The high degree of morphological overlap found among the species underlines the importance of DNA sequence data for species delimitation and description in the genus Cyclophorus. Given the findings of the present study and the high potential that as yet undiscovered cryptic taxa have also evolved in other groups of karst-inhabiting organisms, we argue for a systematic and efficient detection and description of Vietnam’s karst biodiversity to provide a solid basis for future conservation planning.

Cryptic diversity in Lithobateswarszewitschii (Amphibia, Anura, Ranidae)

Lithobateswarszewitschii is a species of ranid frog distributed from southern Honduras to Panama. This species suffered severe population declines at higher elevations (above 500 m a.s.l.) from the 1980s to early 1990s, but there is more recent evidence of recovery in parts of its range. Here we advocate for the status of Lithobateswarszewitschii as a candidate cryptic species complex based on sequence data from mitochondrial genes CO1 and 16S. Using concatenated phylogenies, nucleotide diversity (K2P-π), net between group mean distance (NBGMD) (π_net) and species delimitation methods, we further elucidate cryptic diversity within this species. All phylogenies display polyphyletic lineages within Costa Rica and Panama. At both loci, observed genetic polymorphism (K2P-π) is also high within and between geographic populations, surpassing proposed species threshold values for amphibians. Additionally, patterns of phylogeographic structure are complicated for this species, and do not appear to be explained by geographic barriers or isolation by distance. These preliminary findings suggest L.warszewitschii is a wide-ranging species complex. Therefore, we propose further research within its wider range, and recommend integrative taxonomic assessment is merited to assess species status.

To split or not to split? Multilocus phylogeny and molecular species delimitation of southeast Asian toads (family: Bufonidae)

Background

Recent studies have demonstrated that Bayesian species delimitation based on the multispecies coalescent model can produce inaccurate results by misinterpreting population splits as species divergences. An approach based on the genealogical divergence index (gdi) was shown to be a viable alternative, especially for delimiting allopatric populations where gene flow is low. We implemented these analyses to assess species boundaries in Southeast Asian toads, a group that is understudied and characterized by numerous unresolved species complexes.

Results

Multilocus phylogenetic analyses showed that deep evolutionary relationships including the genera Sigalegalephrynus, Ghatophryne, Parapelophryne, Leptophryne, Pseudobufo, Rentapia, and Phrynoides remain unresolved. Comparison of genetic divergences revealed that intraspecific divergences among allopatric populations of Pelophyrne signata (Borneo vs. Peninsular Malaysia), Ingerophrynus parvus (Peninsular Malaysia vs. Myanmar), and Leptophryne borbonica (Peninsular Malaysia, Java, Borneo, and Sumatra) are consistent with interspecific divergences of other Southeast Asian bufonid taxa. Conversely, interspecific divergences between Pelophryne guentheri/P. api, Ansonia latiffi/A. leptopus, and I. gollum/I. divergens were low (< 3%) and consistent with intraspecific divergences of other closely related taxa. The BPP analysis produced variable results depending on prior settings and priors estimated from empirical data produced the best results that were also congruent with the gdi analysis.

Conclusions

This study showed that the evolutionary history of Southeast Asian toads is difficult to resolve and numerous relationships remain ambiguous. Although some results from the species delimitation analyses were inconclusive, they were nevertheless efficacious at identifying potential new species and taxonomic incompatibilities for future in-depth investigation. We also demonstrated the sensitivity of BPP to different priors and that careful selection priors based on empirical data can greatly improve the analysis. Finally, the gdi can be a robust tool to complement other species delimitation methods.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1422-3) contains supplementary material, which is available to authorized users.

Species Boundaries and Parapatric Speciation in the Complex of Alpine Shrubs, Rosa sericea (Rosaceae), Based on Population Genetics and Ecological Tolerances

Discerning species boundaries among closely related taxa is fundamental to studying evolution and biodiversity. However, species boundaries can be difficult to access in plants because ongoing divergence and speciation may leave an evolutionary footprint similar to introgression, which occurs frequently among species and genera. In this study, we sought to determine species boundaries between two closely related alpine shrubs, Rosa sericea and Rosa omeiensis, using population genetics, environmental data and ecological niche modeling, and morphological traits. We analyzed populations of R. sericea and R. omeiensis using genetic markers comprising a fragment of the single-copy nuclear gene, LEAFY, micro-satellites (EST-SSR), and plastid DNA sequences. The DNA sequence data suggested clusters of populations consistent with geography but not with previously proposed species boundaries based on morphology. Nevertheless, we found that the ecological niches of the previously proposed species only partially overlap. Thus, we suspect that these species are in the process of parapatric speciation; that is, differentiating along an ecological gradient, so that they exhibit differing morphology. Morphology has previously been the basis of recognizing the species R. sericea and R. omeiensis, which are the most widely distributed species within a broader R. sericea complex that includes several other narrow endemics. Here, we recognize R. sericea and R. omeiensis as independent species based on morphological and ecological data under the unified species concept, which emphasizes that these data types are of equal value to DNA for determining species boundaries and refining taxonomic treatments. While the DNA data did not delimit species within the R. sericea complex, we expect to develop and utilize new, robust DNA tools for understanding speciation within this group in future studies.

Assessing Genomic Admixture between Cryptic Plutella Moth Species following Secondary Contact

Cryptic species are genetically distinct taxa without obvious variation in morphology and are occasionally discovered using molecular or sequence data sets of populations previously thought to be a single species. The world-wide Brassica pest, Plutella xylostella (diamondback moth), has been a problematic insect in Australia since 1882, yet a morphologically cryptic species with apparent endemism (P. australiana) was only recognized in 2013. Plutella xylostella and P. australiana are able to hybridize under laboratory conditions, and it was unknown whether introgression of adaptive traits could occur in the field to improve fitness and potentially increase pressure on agriculture. Phylogenetic reconstruction of 29 nuclear genomes confirmed P. xylostella and P. australiana are divergent, and molecular dating with 13 mitochondrial genes estimated a common Plutella ancestor 1.96 ± 0.175 Ma. Sympatric Australian populations and allopatric Hawaiian P. xylostella populations were used to test whether neutral or adaptive introgression had occurred between the two Australian species. We used three approaches to test for genomic admixture in empirical and simulated data sets including 1) the f3 statistic at the level of the population, 2) pairwise comparisons of Nei's absolute genetic divergence (dXY) between populations, and 3) changes in phylogenetic branch lengths between individuals across 50-kb genomic windows. These complementary approaches all supported reproductive isolation of the Plutella species in Australia, despite their ability to hybridize. Finally, we highlight the most divergent genomic regions between the two cryptic Plutella species and find they contain genes involved with processes including digestion, detoxification, and DNA binding.

Hidden diversity within the depauperate genera of the snake tribe Lampropeltini (Serpentes, Colubridae)

Accurate representation of lineage diversity through complete taxon sampling is crucial to understanding the evolution of biodiversity, particularly when using molecular phylogenetics to estimate evolutionary relationships. In this interest, taxonomic diversity is often used as a proxy for lineage diversity even though the two concepts are not synonymous. We explore this within the snake tribe Lampropeltini which includes some of the most conspicuous and heavily studied snakes in North America. Both the taxonomy and hypothesized relationships within this tribe have been in flux. The number of species has increased from 23 to 51 over the last thirty years, predominately within three of the nine genera (Lampropeltis, Pantherophis, Pituophis). The remaining six depauperate genera (Arizona, Bogertophis, Cemophora, Pseudelaphe, Rhinocheilus, and Senticolis) have been poorly represented in phylogenetic studies. To estimate evolutionary relationships and determine if the dichotomy in depauperate and speciose genera within Lampropeltini is a function of taxon sampling or truly represents the lineage diversity, we estimated the phylogeny of this group using nuclear and mitochondrial loci in a concatenated and coalescent framework with the largest sampling of the six depauperate genera to date. In addition, we estimated the divergence dates among the genera to assess whether the instability of Lampropeltini phylogenetic relationships is due to an adaptive radiation. While some nodes still remain unresolved, the generic-level relationships we recovered agree with those of a recent next-generation study that used a much larger set of loci for fewer individuals. We also tested two putative species, Arizona pacata and Pseudelaphe phaescens, for the first time phylogenetically and find evidence that they are distinct lineages. Overall, we find that the taxonomic and genetic diversity are not correlated in Lampropeltini and that representing putative diversity in phylogenies will lead to a better estimate of evolutionary histories, especially in groups with complex radiations.

Cryptic diversity within grass-associated Abacarus species complex (Acariformes: Eriophyidae), with the description of a new species, Abacarus plumiger n. sp

Accurate estimation of species richness is often complex as genetic divergence is not always accompanied by appreciable morphological differentiation. In consequence, cryptic lineages or species evolve. Cryptic speciation is common especially in taxa characterized by small and simplified bodies, what makes their proper identification challenging. The cereal rust mite, Abacarus hystrix, was regarded for a long time as a species associated with a wide range of grass hosts, whereas wide host ranges are rather rare in eriophyoid mites. Therefore, the generalist status of A. hystrix was questioned. In this paper we demonstrate that the diversity within Abacarus species associated with grasses is more complex than it was previously thought. The 78 Abacarus mtDNA COI sequences used in this study formed 10 highly supported clades (bootstrap value 99%) and four more distinct genetic lineages were represented by unique sequences. The genetic distances between them ranged from 6.6 to 26.5%. Moreover, morphological study and genetic approach based on the combination of the Poisson Tree Processes model for species delimitation (PTP) and a Bayesian implementation of PTP (bPTP), and Neighbour Joining analyses led to delimitation of a new species within the Abacarus complex: Abacarus plumiger, specialized on smooth brome (Bromus inermis). Furthermore, our analyses demonstrated a pattern of host-associated differentiation within the complex. Overall, our study indicates that cryptic speciation occurs in the grass-associated Abacarus genus, and suggests the need for more extensive sampling using integrative methods.

A new minute Pristimantis (Amphibia: Anura: Strabomantidae) from the Andes of southern Ecuador

We describe a new rainfrog species (Pristimantis), from the wetland complex Oña, Nabón, Saraguro and Yacuambi, in the Andes of southern Ecuador, at altitudes ranging between 3000–3400 m a.s.l. Pristimantis tiktik sp. nov. is a small frog, displaying sexual dimorphism (the males with dorsum of various shades of gray, brown, orange or green and a whitish or pinkish yellow venter; females with brownish gray or gray dorsum and a reticulated white and black venter), with SVL ranging between 19.7–20.4 mm in females (n = 3) and 16.1–18.4 mm in males (n = 6). The skin on dorsum is tuberculated, that on venter is coarsely areolate, dorsolateral folds are absent, tympanic membrane is absent but the tympanic annulus is evident, cranial crests are absent, discs on fingers just slightly expanded, heel is lacking enlarged tubercles, inner edge of tarsus is bearing a long fold, Toe V is slightly longer than Toe III and the iris coloration is bronze with fine black reticulations. The males have a large subgular vocal sac that extends onto the chest and vocal slits but lack nuptial pads. The unique advertisement call consists of long duration series of periodically repeated clicks: “tik”. Molecular analyses place the new species in the recently resurrected P. orestes group, as the sister species of the assemblage that contains P. bambu, P. mazar, P. simonbolivari and an undescribed species.