PCR primers and amplification methods for 12S ribosomal DNA, the control region, cytochrome oxidase I, and cytochrome b in bufonids and other frogs, and an overview of PCR primers which have amplified DNA in amphibians successfully

Deep divergences among inconspicuously colored clades of Epipedobates poison frogs

Poison frogs (Dendrobatidae) are famous for their aposematic species, having a combination of diverse color patterns and defensive skin toxins, yet most species in this family are inconspicuously colored and considered non-aposematic. Epipedobates is among the youngest genus-level clades of Dendrobatidae that includes both aposematic and inconspicuous species. Using Sanger-sequenced mitochondrial and nuclear markers, we demonstrate deep genetic divergences among inconspicuous species of Epipedobates but relatively shallow genetic divergences among conspicuous species. Our phylogenetic analysis includes broad geographic sampling of the inconspicuous lineages typically identified as E. boulengeri and E. espinosai, which reveals two putative new species, one in west-central Colombia (E. sp. 1) and the other in north-central Ecuador (E. aff. espinosai). We conclude that E. darwinwallacei is a junior subjective synonym of E. espinosai. We also clarify the geographic distributions of inconspicuous Epipedobates species including the widespread E. boulengeri. We provide a qualitative assessment of the phenotypic diversity in each nominal species, with a focus on the color and pattern of inconspicuous species. We conclude that Epipedobates contains eight known valid species, six of which are inconspicuous. A relaxed molecular clock analysis suggests that the most recent common ancestor of Epipedobates is ∼11.1 million years old, which nearly doubles previous estimates. Last, genetic information points to a center of species diversity in the Chocó at the southwestern border of Colombia with Ecuador. A Spanish translation of this text is available in the supplementary materials.

Easternmost distribution of Bufo bufo (Linnaeus, 1758) in Türkiye: implications for the putative contact zone between B. bufo and B. verrucosissimus

The geographic range of a species is crucial for obtaining information on the exact distribution of the species. The geographic data are important for delimiting distinct species or exploring the degree of differentiation among different populations of a species. The local details of species boundaries facilitate the study of the importance of phylogeographic background, secondary contacts, and hybrid zones, along with the relations between the species and its extrinsic environmental factors. In the present study, the range boundaries of Bufo bufo and Bufo verrucosissimus in the north-eastern region of Türkiye were delineated using an integrative taxonomic approach that utilized a combination of molecular and morphological data. According to the mtDNA results of the present study, B. bufo inhabits a single distribution from İyidere town to Çayeli town in Rize, while B. verrucosissimus is distributed from Şavşat town of Artvin to Ardeşen town in Rize. In addition, the two species coexist in Pazar, Hemşin, and Çamlıhemşin towns in Rize. The demographic analyses indicated a distinct population expansion for the B. verrucosissimus species after the Last Glacial Maximum, while the same did not occur for B. bufo. The univariate and multivariate statistical analyses conducted for the morphological data of the two species corroborated the presence of a putative contact zone between B. bufo and B. verrucosissimus. In summary, the present study resolved the non-distinct geographic boundaries between B. bufo and B. verrucosissimus species and also revealed the easternmost distribution of B. bufo in Türkiye. In addition, important evidence on the putative contact zone between the two species was indicated using an integrative taxonomic approach.

Evaluation of partial 12S rRNA, 16S rRNA, COI and Cytb gene sequence datasets for potential single DNA barcode for hylids (Anura: Hylidae)

We evaluated the extent of intraspecific and interspecific genetic distances and the effectiveness of predefined threshold values using the main genes for estimates of biodiversity and specimens' identification in anurans. Partial sequences of the mitochondrial genes for small (12S) and large (16S) ribosomal subunits, cytochrome c oxidase subunit I (COI) and Cytochrome b (Cytb) of the family Hylidae were downloaded from GenBank and curated for length, coverage, and potential contaminations. We performed analyses for all sequences of each gene and the same species present in these datasets by distance and tree (monophyly)-based evaluations. We also evaluated the ability to identify specimens using these datasets applying "nearest neighbor" (NN), "best close match" (BCM) and "BOLD ID" tests. Genetic distance thresholds were generated by the function 'threshVal' and "localMinima" from SPIDER package and traditional threshold values (1%, 3%, 6% and 10%) were also evaluated. Coding genes, especially COI, had a better identification capacity than non-coding genes on barcoding gap and monophyly analysis and NN, BCM, BOLD ID tests. Considering the multiple factors involved in global DNA barcoding evaluations, we present a critical assessment of the use of these genes for biodiversity estimation and specimens' identification in anurans (e.g. hylids).

Microfluidic Continuous Flow DNA Fragmentation based on a Vibrating Sharp-tip

Fragmentation of DNA into short fragments is of great importance for detecting and studying DNAs. Current microfluidic methods of DNA fragmentation are either inefficient for generating small fragments or rely on microbubbles. Here, we report a DNA fragmentation method in a 3D-printed microfluidic device, which allows efficient continuous flow fragmentation of genomic DNAs without the need for microbubbles. This method is enabled by localized acoustic streaming induced by a single vibrating sharp-tip. Genomic DNAs were fragmented into 700 to 3000 bp fragments with a low power consumption of ~140 mW. The system demonstrated successful fragmentation under a wide range of flow rates from 1 to 50 μL/min without the need for air bubbles. Finally, the utility of the continuous DNA fragmentation method was demonstrated to accelerate the DNA hybridization process for biosensing. Due to the small footprint, continuous flow and bubble-free operation, and high fragmentation efficiency, this method demonstrated great potential for coupling with other functional microfluidic units to achieve an integrated DNA analysis platform.

Six new species of Pristimantis (Anura: Strabomantidae) from Llanganates National Park and Sangay National Park in Amazonian cloud forests of Ecuador

We describe six new species of rainfrogs of the genus Pristimantis (Strabomantidae) from Amazonian cloud forests in Ecuador. We also present a phylogeny showing the relationships of the new species. The phylogeny is based on mitochondrial genes 16S rRNA (16S), 12 rRNA (12S), NADH-ubiquinone oxidoreductase chain 1 (ND1) and the nuclear gene recombination-activating 1 (RAG1). We also describe the osteology of two of the new species using high-resolution x-ray computed tomography. The new species belong to two clades. The first clade is sister to the subgenus Huicundomantis and includes P. tamia sp. nov., P. miktos, and P. mallii. Pristimantis tamia sp. nov. is morphologically similar to P. miktos, P. mallii, P. martiae, and P. incomptus, but differs from them by lacking vocal slits and tympanic membrane and by having light greenish blue iris. Based in our results we expand the subgenus Huicundomantis to include the P. miktos species group. The second clade is remarkable by being highly divergent and consisting exclusively of new species: P. anaiae sp. nov., P. glendae sp. nov., P. kunam sp. nov., P. resistencia sp. nov., and P. venegasi sp. nov. The new species resemble P. roni, P. yanezi, P. llanganati, P. katoptroides, P. verecundus, and P. mutabilis but can be distinguished from them by lacking vocal slits and tympanic membrane and by having large dark round areas with thin clear borders in the sacral region. All six new species occur in the eastern slopes of the Ecuadorian Andes and are known from a single locality in Llanganates or Sangay National Park. We recommend assigning all of them to the Data Deficient (DD) Red List category. Based in our high-resolution x-ray tomographies, we report the presence of structures that appear to be intercalary elements. This would be the first report of such structures in Terrarana.

Genetic structure of American bullfrog populations in Brazil

Non-native species are a major problem affecting numerous biomes around the globe. Information on their population genetics is crucial for understanding their invasion history and dynamics. We evaluated the population structure of the non-native American bullfrog, Aquarana catesbeiana, in Brazil on the basis of 324 samples collected from feral and captive groups at 38 sites in seven of the nine states where feral populations occur. We genotyped all samples using previously developed, highly polymorphic microsatellite loci and performed a discriminant analysis of principal components together with Jost’s D index to quantify pairwise differentiation between populations. We then amplified 1,047 base pairs of the mitochondrial cytochrome b (cytb) gene from the most divergent samples from each genetic population and calculated their pairwise differences. Both the microsatellite and cytb data indicated that bullfrogs comprise two populations. Population grouping 1 is widespread and possesses two cytb haplotypes. Population grouping 2 is restricted to only one state and possesses only one of the haplotypes from Population grouping 1. We show that there were two imports of bullfrogs to Brazil and that there is low genetic exchange between population groupings. Also, we find that there is no genetic divergence among feral and captive populations suggesting continuous releases. The limited genetic variability present in the country is associated to the small number of introductions and founders. Feral bullfrogs are highly associated to leaks from farms, and control measures should focus on preventing escapes using other resources than genetics, as feral and captive populations do not differ.

Widely used, short 16S rRNA mitochondrial gene fragments yield poor and erratic results in phylogenetic estimation and species delimitation of amphibians

Background

The 16S mitochondrial rRNA gene is the most widely sequenced molecular marker in amphibian systematic studies, making it comparable to the universal CO1 barcode that is more commonly used in other animal groups. However, studies employ different primer combinations that target different lengths/regions of the 16S gene ranging from complete gene sequences (~ 1500 bp) to short fragments (~ 500 bp), the latter of which is the most ubiquitously used. Sequences of different lengths are often concatenated, compared, and/or jointly analyzed to infer phylogenetic relationships, estimate genetic divergence (p-distances), and justify the recognition of new species (species delimitation), making the 16S gene region, by far, the most influential molecular marker in amphibian systematics. Despite their ubiquitous and multifarious use, no studies have ever been conducted to evaluate the congruence and performance among the different fragment lengths.

Results

Using empirical data derived from both Sanger-based and genomic approaches, we show that full-length 16S sequences recover the most accurate phylogenetic relationships, highest branch support, lowest variation in genetic distances (pairwise p-distances), and best-scoring species delimitation partitions. In contrast, widely used short fragments produce inaccurate phylogenetic reconstructions, lower and more variable branch support, erratic genetic distances, and low-scoring species delimitation partitions, the numbers of which are vastly overestimated. The relatively poor performance of short 16S fragments is likely due to insufficient phylogenetic information content.

Conclusions

Taken together, our results demonstrate that short 16S fragments are unable to match the efficacy achieved by full-length sequences in terms of topological accuracy, heuristic branch support, genetic divergences, and species delimitation partitions, and thus, phylogenetic and taxonomic inferences that are predicated on short 16S fragments should be interpreted with caution. However, short 16S fragments can still be useful for species identification, rapid assessments, or definitively coupling complex life stages in natural history studies and faunal inventories. While the full 16S sequence performs best, it requires the use of several primer pairs that increases cost, time, and effort. As a compromise, our results demonstrate that practitioners should utilize medium-length primers in favor of the short-fragment primers because they have the potential to markedly improve phylogenetic inference and species delimitation without additional cost.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-01994-y.

Peripheral Auditory System Divergence Does Not Explain Species Differences in Call Preference

Receiver sensory systems have long been cited as an important source of variation in mate preferences that could lead to signal diversification and behavioral isolation between lineages, with a general assumption that animals prefer the most conspicuous signals. The matched filter hypothesis posits that tuning of the frog peripheral auditory system matches dominant frequencies in advertisement calls used to attract mates. However, little work has characterized species with frequency modulation in their calls. In this study, we extend prior work characterizing the lack of correlated evolution between auditory tuning and spectral properties of male calls in Engystomops (=Physalaemus) frogs. We analyze auditory sensitivity of three cryptic species that differ consistently in female mate preferences for calls of different frequencies. The audiograms of these species differ, but the frequency at which the frog is maximally sensitive is not the most relevant difference in tuning of the auditory periphery. Rather, we identify species differences in overall sensitivity within specific frequency ranges, and we model the effects of these sensitivity differences on neural responses to natural calls. We find a general mismatch between auditory brainstem responses and behavioral preferences of these taxa and rule out the matched filter hypothesis as explaining species differences in male calls and mate preferences in this group.

A new cryptic species of the Pristimantis lacrimosus group (Anura, Strabomantidae) from the eastern slopes of the Ecuadorian Andes

With 566 species, the neotropical genus Pristimantis is the most speciose vertebrate genus. As a result of its striking diversity, taxonomic reviews remain a challenge. Herein, we present an updated phylogeny of the Pristimantis lacrimosus group and describe a new species from Llanganates and Sangay National Parks. We also report, for the first time, the phylogenetic position of Pristimantis degener, P. eugeniae, P. katoptroides, and P. petersi. Based on our phylogeny, we add two species to the Pristimantis lacrimosus group. Through the integration of molecular and bioacoustic evidence, we describe a new species which was hidden under “Pristimantis petersi”. Pristimantis petersioidessp. nov. is most closely related to Pristimantis petersi and an undescribed species from Peru. It can be distinguished from P. petersi by its advertisement call and large genetic differences (uncorrected p-genetic distances 7.9% to 8.4% for gene 16S). Moreover, the new species and P. petersi are not sister species. We suggest assigning the new species to the Endangered Red List category because it has a small distribution range with deforestation as result of agriculture and other anthropogenic influences.

Revisiting the phylogenetic predicament of the genus Huia (Amphibia: Ranidae) using molecular data and tadpole morphology

Despite a considerable amount of research, the systematics of the ranid genus Huia have remained unresolved, mostly owing to insufficient sampling and morphological similarities. As currently circumscribed, Huia consists of five species, but multiple studies have consistently demonstrated that it is not a monophyletic genus. However, no study has approached the problem with adequate data and provided a systematically sound solution, leaving the genus to languish in a classification that is phylogenetically incoherent. We generated the most comprehensive sampling of Huia to date, based on extensive fieldwork in Java and Sumatra. Using an integrative taxonomy framework, we analysed four mitochondrial and two nuclear markers and, in conjunction with tadpole morphology, investigated the phylogenetics of Huia and its congeners. Corroborating previous studies, Huia is recovered as a paraphyletic group. Huia cavitympanum emerges as the sister taxon to Meristogenys. The remaining members of Huia form a monophyletic group, sister to the H. cavitympanum + Meristogenys clade. Our extensive geographical sampling in Sumatra and Java reveals multiple highly divergent lineages that potentially represent undescribed species. Using our expanded molecular and morphological dataset, we resolve the paraphyly of Huia by restricting the genus to its type species and propose a new genus to accommodate the strongly supported clade of Sumatran and Javan populations previously belonging to Huia.

Cytonuclear discordance, reticulation and cryptic diversity in one of North America's most common frogs

Complicated phylogenetic histories benefit from diverse sources of inference. Pseudacris crucifer (spring peeper) spans most of eastern North America and comprises six mtDNA lineages that form multiple contact zones. The putative Miocene or early Pliocene origins of the oldest lineages within Pseudacris crucifer imply sufficient time for species-level divergence. To understand why this species appears unified while congeners have radiated, we analyze and compare male advertisement calls, mitochondrial, and nuclear markers and speak to the complex processes that have potentially influenced its contemporary patterns. We find extensive geographic and topological mitonuclear discordance, with three nuclear lineages containing 6 more-structured mtDNA lineages, and nuclear introgression at some contact zones. Male advertisement call differentiation is incongruent with the genetic structure as only one lineage appears differentiated. Occupying the Interior Highlands of the central United States, this Western lineage also has the most concordant mitochondrial and nuclear geographic patterns. Based on our findings we suggest that the antiquity of common ancestors was not as important as the maintenance of allopatry in the divergence in P. crucifer genetic lineages. We use multiple lines of evidence to generate hypotheses of isolation, reticulation, and discordance within this species and to expand our understanding of the early stages of speciation.

Three new species of frogs of the genus Pristimantis (Anura, Strabomantidae) with a redefinition of the P. lacrimosus species group

A new phylogeny for the Pristimantis lacrimosus species group is presented, its species content reviewed, and three new species described from the eastern slopes of the Ecuadorian Andes. Our phylogeny includes, for the first time, samples of P. aureolineatus, P. bromeliaceus, and P. lacrimosus. The morphology of hyperdistal subarticular tubercles is also assessed among 21 species of Pristimantis. The P. lacrimosus species group is composed of 36 species distributed in the Chocó, Guiana, and Amazon regions of tropical South America with a single species reaching Central America. Ancestral area reconstruction indicates that, despite its high diversity in the Amazon region, the P. lacrimosus group originated in the Pacific basin, Chocó region of Ecuador and Colombia. Pristimantis amaguanae sp. nov. is most closely related to P. bromeliaceus. It differs from P. bromeliaceus by being smaller, having transversal dark bands in the hindlimbs (absent or faint in P. bromeliaceus) and the absence of discoidal fold (present in P. bromeliaceus). Pristimantis nankints sp. nov. and P. romeroae sp. nov. are part of a clade of predominantly light-green frogs that includes P. acuminatus, P. enigmaticus, P. limoncochensis, and P. omeviridis. Pristimantis nankints sp. nov. and P. romeroae sp. nov. can be distinguished from all of them by the presence of a dark dorsolateral stripe that borders a light green band on a green background. Hyperdistal tubercles are present in all examined species of the P. lacrimosus species group and its sister clade. Species with hyperdistal tubercles are characterized by having relatively long terminal phalanges and narrow T-shaped expansion at the end of the terminal phalange. We discuss the phylogenetic distribution of these characters and their potential diagnostic significance.

Skull diversity and evolution in miniaturized amphibians, genus Brachycephalus (Anura: Brachycephalidae)

Miniaturized amphibians of the genus Brachycephalus are phenotypically diverse. The species of Brachycephalus have bufoniform or leptodactyliform Baupläne and any of three skeletal states: nonhyperossified, hyperossified without dorsal shield, and hyperossified with dorsal shield. We integrate high-resolution microcomputed tomography, geometric morphometrics, and an estimate of molecular phylogenetic relationships to investigate skull diversity in shape and size-shape space in selected species of Brachycephalus. Skull diversity amongst species of Brachycephalus can be partitioned into shape and size-shape space according to the four conditions of skeletal states-Baupläne, namely, nonhyperossified leptodactyliform, nonhyperossified bufoniform, hyperossified bufoniform without dorsal shield, and hyperossified bufoniform with dorsal shield. Skull diversity in shape and size-shape space in nonhyperossified leptodactyliform species of Brachycephalus is markedly larger, when compared to skull diversity in species of the three other conditions of skeletal states-Baupläne. Variation in skull shape scales with size across Brachycephalus and, therefore, can be explained by allometry. Skull diversity, Baupläne, and skeletal states covary to a large extent with monophyletic lineages of Brachycephalus, as revealed by a mitochondrial DNA species tree. Nonhyperossified bufoniform species and hyperossified bufoniform species with or without dorsal shield are monophyletic lineages, as inferred from a mitochondrial DNA species tree. Nonhyperossified leptodactyliform species of Brachycephalus do not share, however, a most recent common ancestor. The nonhyperossified leptodactyliform species of Brachycephalus, due to their marked skull diversity and lack of monophyly, emerge as evolutionarily complex. Therefore, further sampling of the nonhyperossified leptodactyliform condition of skeletal states-Baupläne will be necessary to further understand the evolutionary history of Brachycephalus.

Phylogeographic analysis and species distribution modelling of the wood frog Batrachyla leptopus (Batrachylidae) reveal interglacial diversification in south western Patagonia

Background

The evolutionary history of southern South American organisms has been strongly influenced by Pleistocene climate oscillations. Amphibians are good models to evaluate hypotheses about the influence of these climate cycles on population structure and diversification of the biota, because they are sensitive to environmental changes and have restricted dispersal capabilities. We test hypotheses regarding putative forest refugia and expansion events associated with past climatic changes in the wood frog Batrachyla leptopus distributed along ∼1,000 km of length including glaciated and non-glaciated areas in southwestern Patagonia.

Methods

Using three mitochondrial regions (D-loop, cyt b, and coI) and two nuclear loci (pomc and crybA1), we conducted multilocus phylogeographic analyses and species distribution modelling to gain insights of the evolutionary history of this species. Intraspecific genealogy was explored with maximum likelihood, Bayesian, and phylogenetic network approaches. Diversification time was assessed using molecular clock models in a Bayesian framework, and demographic scenarios were evaluated using approximate Bayesian computation (ABC) and extended Bayesian skyline plot (EBSP). Species distribution models (SDM) were reconstructed using climatic and geographic data.

Results

Population structure and genealogical analyses support the existence of four lineages distributed north to south, with moderate to high phylogenetic support (Bootstrap > 70%; BPP > 0.92). The diversification time of B. leptopus’ populations began at ∼0.107 mya. The divergence between A and B lineages would have occurred by the late Pleistocene, approximately 0.068 mya, and divergence between C and D lineages was approximately 0.065 mya. The ABC simulations indicate that lineages coalesced at two different time periods, suggesting the presence of at least two glacial refugia and a postglacial colonization route that may have generated two southern lineages (p = 0.93, type I error: <0.094, type II error: 0.134). EBSP, mismatch distribution and neutrality indexes suggest sudden population expansion at ∼0.02 mya for all lineages. SDM infers fragmented distributions of B. leptopus associated with Pleistocene glaciations. Although the present populations of B. leptopus are found in zones affected by the last glacial maximum (∼0.023 mya), our analyses recover an older history of interglacial diversification (0.107–0.019 mya). In addition, we hypothesize two glacial refugia and three interglacial colonization routes, one of which gave rise to two expanding lineages in the south.

Lost and found: Frogs in a biodiversity hotspot rediscovered with environmental DNA

Declines and extinctions are increasing globally and challenge conservationists to keep pace with biodiversity monitoring. Organisms leave DNA traces in the environment, e.g., in soil, water, and air. These DNA traces are referred to as environmental DNA (eDNA). The analysis of eDNA is a highly sensitive method with the potential to rapidly assess local diversity and the status of threatened species. We searched for DNA traces of 30 target amphibian species of conservation concern, at different levels of threat, using an environmental DNA metabarcoding approach, together with an extensive sequence reference database to analyse water samples from six montane sites in the Atlantic Coastal Forest and adjacent Cerrado grasslands of Brazil. We successfully detected DNA traces of four declined species (Hylodes ornatus, Hylodes regius, Crossodactylus timbuhy, and Vitreorana eurygnatha); two locally disappeared (Phasmahyla exilis and Phasmahyla guttata); and one species that has not been seen since 1968 (putatively assigned to Megaelosia bocainensis). We confirm the presence of species undetected by traditional methods, underscoring the efficacy of eDNA metabarcoding for biodiversity monitoring at low population densities, especially in megadiverse tropical sites. Our results support the potential application of eDNA in conservation biology, to evaluate persistence and distribution of threatened species in surveyed habitats or sites, and improve accuracy of red lists, especially for species undetected over long periods.

Rapid mitochondrial genome sequencing based on Oxford Nanopore Sequencing and a proxy for vertebrate species identification

Molecular information is crucial for species identification when facing challenging morphology-based specimen identifications. The use of DNA barcodes partially solves this problem, but in some cases when PCR is not an option (i.e., primers are not available, problems in reaction standardization), amplification-free approaches could be an optimal alternative. Recent advances in DNA sequencing, like the MinION device from Oxford Nanopore Technologies (ONT), allow to obtain genomic data with low laboratory and technical requirements, and at a relatively low cost. In this study, we explore ONT sequencing for molecular species identification from a total DNA sample obtained from a neotropical rodent and we also test the technology for complete mitochondrial genome reconstruction via genome skimming. We were able to obtain "de novo" the complete mitogenome of a specimen from the genus Melanomys (Cricetidae: Sigmodontinae) with average depth coverage of 78X using ONT-only data and by combining multiple assembly routines. Our pipeline for an automated species identification was able to identify the sample using unassembled sequence data (raw) in a reasonable computing time, which was substantially reduced when a priori information related to the organism identity was known. Our findings suggest ONT sequencing as a suitable candidate to solve species identification problems in metazoan nonmodel organisms and generate complete mtDNA datasets.

Development of a quantitative polymerase chain reaction assay and environmental DNA sampling methods for Giant Gartersnake (Thamnophis gigas)

The Giant Gartersnake (Thamnophis gigas) is a low density visually evasive species with a low detection probability based on standard field survey methods (e.g., traps, visual census). Habitat loss has resulted in extirpations or serious declines for T. gigas populations throughout the southern two thirds of its historic range. Uncertainty regarding its current distribution and occupancy present management challenges for the species. Enhancing survey sensitivity through development of environmental DNA sampling (eDNA) methods would improve compliance monitoring under the Endangered Species Act, recovery planning for T. gigas, and evaluation of California’s Central Valley tule marsh habitat on which this species depends. To address these needs, we designed and validated diagnostic quantitative Polymerase Chain Reaction (qPCR) assays for identifying portions of the Cytochrome B (CytB) and the Nicotinamide adenine dinucleotide (NADH) dehydrogenase subunit 4 (ND4) genes of the T. gigas mitochondrial genome. The designed ND4 qPCR assay was not specific to T. gigas DNA and amplified DNA from a closely related and spatially co-occurring Thamnophis species (T.s. fitchi). The CytB T. gigas qPCR assay proved specific to a species level with a sensitivity that reliably detected T. gigas DNA at a concentration of 2.0x10^-5 ng μL^-1. To assess detection range, coordinated field sampling was conducted at aquatic sites with an observed and documented population of T. gigas. The T. gigas qPCR assay reliably detected DNA from samples taken 300m downstream from the known source. We then used environmental eDNA sampling and qPCR analysis to augment unsuccessful trap surveys in the southern range of T. gigas and detected DNA in 28 of the 52 locations sampled, confirming that T. gigas was still present at some sites where physical trapping failed to identify presence. QPCR-based DNA detection coupled with eDNA sampling methods provides an effective means to obtain critical population metrics from this otherwise cryptic, federally protected and hard to study organism, offering great promise for elucidating patterns of occupancy with greater efficiency and at far less cost than trapping methods, particularly where detection probabilities are low.

Identification of leopard frogs (Anura: Ranidae: Lithobates) distributed in some localities of the Southern Mexican Plateau using mitochondrial DNA sequences

Species of the genus Lithobates are difficult to identify, especially on the 'Rana pipiens' complex, because the morphological differences among the species are often subtle. In this study, we used 12S ribosomal RNA gene partial sequences to identify specimens of leopard frogs from some new localities (not sampled on previous studies) of the Southern Mexican Plateau, to carry out a phylogenetic analysis, and also a morphometric analysis of some morphological features were analyzed to evaluate their morphological variation. A phylogenetic analysis using partial sequences of 12S rDNA mitochondrial gene for 31 samples was carried out using Bayesian Inference, Maximum Parsimony, and Maximum Likelihood. In addition, 20 morphological lineal measurements from 97 specimens were analyzed by principal component analysis (PCA) to evaluate if the frogs studied are conspecific. Partial sequences of the 12S rDNA obtained from frogs distributed in the Southern Mexican Plateau show two haplotypes with a divergence of 0.4%. Phylogenetic hypothesis shows an exclusive group with a previous sequence of Lithobates montezumae. The PCA indicates that variables are not linearly correlated and specimens belong to a single group. Evidence found, let us consider that specimens studied belong to one species of the Lithobates montezumae subgroup, distributed principally in the Southern Mexican Plateau.

A new treefrog from Cordillera del Cóndor with comments on the biogeographic affinity between Cordillera del Cóndor and the Guianan Tepuis (Anura, Hylidae, Hyloscirtus)

The Hyloscirtuslarinopygion group is a clade of 16 species of large hylids that inhabit cascading Andean streams. They have brown coloration that, in most species, contrasts with bright marks. Herein morphological and genetic evidence is used to describe a new species of the group from Cordillera del Cóndor, a sub-Andean mountain chain that has phytogeographic affinities with the Guianan Tepuis. The new species is characterized by dark-brown coloration with contrasting bright orange flecks and by the presence of an enlarged and curved prepollex protruding as a spine. The new species is closely related to H.tapichalaca and an undescribed species from the southern Andes of Ecuador. The genetic distance between H.hillisi sp. n. and its closest relative, H.tapichalaca, is 2.9% (gene 16S mtDNA). Our phylogeny and a review of recently published phylogenies show that amphibians from Cordillera del Cóndor have close relationships with either Andean or Amazonian species. Amphibians do not show the Condor-Guianan Tepuis biogeographic link that has been documented in plants.

Speciation in a biodiversity hotspot: Phylogenetic relationships, species delimitation, and divergence times of Patagonian ground frogs from the Eupsophus roseus group (Alsodidae)

The alsodid ground frogs of the Eupsophus genus are divided into two groups, the roseus (2n = 30) and vertebralis (2n = 28), which are distributed throughout the temperate Nothofagus forests of South America. Currently, the roseus group is composed by four species, while the vertebralis group consists of two. Phylogenetic relationships and species delimitation within each group are controversial. In fact, previous analyses considered that the roseus group was composed of between four to nine species. In this work, we evaluated phylogenetic relationships, diversification times, and species delimitation within the roseus group using a multi-locus dataset. For this purpose, mitochondrial (D-loop, Cyt b, and COI) and nuclear (POMC and CRYBA1) partial sequences from 164 individuals were amplified, representing all species. Maximum Likelihood (ML) and Bayesian approaches were used to reconstruct phylogenetic relationships. Species tree was estimated using BEAST and singular value decomposition scores for species quartets (SVDquartets). Species limits were evaluated with six coalescent approaches. Diversification times were estimated using mitochondrial and nuclear rates with LogNormal relaxed clock in BEAST. Nine well-supported monophyletic lineages were recovered in Bayesian, ML, and SVDquartets, including eight named species and a lineage composed by specimens from the Villarrica population (Bootstrap:>70, PP:> 0.99). Single-locus species delimitation analyses overestimated the species number in E. migueli, E. calcaratus, and E. roseus lineages, while multi-locus analyses recovered as species the nine lineages observed in phylogenetic analyses (Ctax = 0.69). It is hypothesized that Eupsophus diversification occurred during Mid-Pleistocene (0.42-0.14 Mya), with most species having originated after the Last Southern Patagonian Glaciation (0.18 Mya). Our results revitalize the hypothesis that the E. roseus group is composed of eight species and support the Villarrica lineage as a new putative species.

Identification of Arenin, a Novel Kunitz-Like Polypeptide from the Skin Secretions of Dryophytes arenicolor

Amphibian skin secretions are enriched with complex cocktails of bioactive molecules such as proteins, peptides, biogenic amines, alkaloids guanidine derivatives, steroids and other minor components spanning a wide spectrum of pharmacological actions exploited for centuries in folk medicine. This study presents evidence on the protein profile of the skin secretions of the canyon tree frog, Dryophytes arenicolor. At the same time, it presents the reverse-phase liquid chromatography isolation, mass spectrometry characterization and identification at mRNA level of a novel 58 amino acids Kunitz-like polypeptide from the skin secretions of Dryophytes arenicolor, arenin. Cell viability assays performed on HDFa, CaCo2 and MCF7 cells cultured with different concentrations of arenin showed a discrete effect at low concentrations (2, 4, 8 and 16 µg/mL) suggesting a multi-target interaction in a hormetic-like dose-response. Further work is required to investigate the mechanisms underlying the variable effect on cell viability produced by different concentrations of arenin.

Advertisement calls and DNA sequences reveal a new species of Scinax (Anura: Hylidae) on the Pacific lowlands of Ecuador

Scinax is a speciose genus of Neotropical hylid frogs. We describe a new species from western Ecuador (elevations between 0 and 1207 m) using morphology, vocalizations, and DNA sequences. We also present a new phylogeny for Scinax based on mitochondrial DNA genes 12S rRNA, Cytochrome Oxidase sub-unit I, Cytochrome B, 16S rRNA, NADH dehydrogenase subunit 1, and adjacent tRNAs. The new species, Scinax tsachila sp. nov. was previously confused with S. quinquefasciatus, a morphologically similar sympatric species. They differ by having markedly different advisement calls, distinct skin texture in the dorsum, and different bone coloration. The new species is sister to S. elaeochroa, a species that differs in advertisement call and color pattern. We provide an updated species account for Scinax quinquefasciatus and a redescription of its holotype.

Mitochondrial genomes of the South American frogs Eupsophus vertebralis and E. emiliopugini (Neobatrachia: Alsodidae) and their phylogenetic relationships

We report the sequencing and compare the mitochondrial genomes of the South American ground frogs Eupsophus vertebralis and E. emiliopugini and reconstruct phylogenetic relationships among Eupsophus species. These genomes consist of 16,156 and 16,711 bp in length, respectively and contain 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes (tRNA), and partial non-coding D-loop region. Both genomes share 94.5% identity with 879 variable sites. A phylogenetic analysis with other available mitogenomes recovered both species as the sister clade of Alsodes gargola. Sequences from D-loop, CO1, and Cyt b, amplified and sequenced with primers developed from the mitochondrial genomes, allowed us to reconstruct phylogenetic relationships among Eupsophus species. Since our report represents the first mitogenomes for the genus Eupsophus, we expect these data will be valuable for further studies on conservation genetics and on the evolution of Patagonian amphibians.

Phylogeography of Aegean green toads (Bufo viridis subgroup): continental hybrid swarm vs. insular diversification with discovery of a new island endemic

Background

Debated aspects in speciation research concern the amount of gene flow between incipient species under secondary contact and the modes by which post-zygotic isolation accumulates. Secondary contact zones of allopatric lineages, involving varying levels of divergence, provide natural settings for comparative studies, for which the Aegean (Eastern Mediterranean) geography offers unique scenarios. In Palearctic green toads (Bufo viridis subgroup or Bufotes), Plio-Pleistocene (~ 2.6 Mya) diverged species show a sharp transition without contemporary gene flow, while younger lineages, diverged in the Lower-Pleistocene (~ 1.9 Mya), admix over tens of kilometers. Here, we conducted a fine-scale multilocus phylogeographic analysis of continental and insular green toads from the Aegean, where a third pair of taxa, involving Mid-Pleistocene diverged (~ 1.5 Mya) mitochondrial lineages, earlier tentatively named viridis and variabilis, (co-)occurs.

Results

We discovered a new lineage, endemic to Naxos (Central Cyclades), while coastal islands and Crete feature weak genetic differentiation from the continent. In continental Greece, both lineages, viridis and variabilis, form a hybrid swarm, involving massive mitochondrial and nuclear admixture over hundreds of kilometers, without obvious selection against hybrids.

Conclusions

The genetic signatures of insular Aegean toads appear governed by bathymetry and Quaternary sea level changes, resulting in long-term isolation (Central Cyclades: Naxos) and recent land-bridges (coastal islands). Conversely, Crete has been isolated since the end of the Messinian salinity crisis (5.3 My) and Cretan populations thus likely result from human-mediated colonization, at least since Antiquity, from Peloponnese and Anatolia. Comparisons of green toad hybrid zones support the idea that post-zygotic hybrid incompatibilities accumulate gradually over the genome. In this radiation, only one million years of divergence separate a scenario of complete reproductive isolation, from a secondary contact resulting in near panmixia.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1179-0) contains supplementary material, which is available to authorized users.

Molecular phylogenetic analysis of a taxonomically unstable ranid from Sumatra, Indonesia, reveals a new genus with gastromyzophorous tadpoles and two new species

The presence of an adhesive abdominal sucker (gastromyzophory) allows tadpoles of certain species of anurans to live in fast-flowing streams. Gastromyzophorous tadpoles are rare among anurans, known only in certain American bufonids and Asian ranids. To date,Huia sumatrana, which inhabits cascading streams, has been the only Sumatran ranid known to possess gastromyzophorous tadpoles. In the absence of thorough sampling and molecular barcoding of adults and larvae, it has remained to be confirmed whether other Sumatran ranid species living in similar habitats, i.e.,Chalcorana crassiovis, possesses this larval type. Moreover, the taxonomic status of this species has long been uncertain and its taxonomic position within the Ranidae, previously based exclusively on morphological characters, has remained unresolved. To study the diversity and relationships of these frogs and to establish the identity of newly collected gastromyzophorous tadpoles from Sumatra, we compared genetic sequences ofC. crassiovis-like taxa from a wide range of sites on Sumatra. We conducted bayesian and maximum likelihood phylogenetic analyses on a concatenated dataset of mitochondrial (12S rRNA, 16S rRNA, and tRNAval) and nuclear (RAG1 and TYR) gene fragments. Our analyses recoveredC. crassiovisto be related toClinotarsus,Huia, andMeristogenys. The DNA barcodes of the gastromyzophorous tadpoles matched adults from the same sites. Herein, we provide a re-description of adultC. crassiovisand propose “C. kampeni” as a synonym of this species. The molecular evidence, morphological features, and distribution suggest the presence of two related new species. The two new species andC. crassiovistogether represent a distinct phylogenetic clade possessing unique molecular and morphological synapomorphies, thus warranting a new genus.

Profound genetic divergence and asymmetric parental genome contributions as hallmarks of hybrid speciation in polyploid toads

The evolutionary causes and consequences of allopolyploidization, an exceptional pathway to instant hybrid speciation, are poorly investigated in animals. In particular, when and why hybrid polyploids versus diploids are produced, and constraints on sources of paternal and maternal ancestors, remain underexplored. Using the Palearctic green toad radiation (including bisexually reproducing species of three ploidy levels) as model, we generate a range-wide multi-locus phylogeny of 15 taxa and present four new insights: (i) at least five (up to seven) distinct allotriploid and allotetraploid taxa have evolved in the Pleistocene; (ii) all maternal and paternal ancestors of hybrid polyploids stem from two deeply diverged nuclear clades (6 Mya, 3.1-9.6 Mya), with distinctly greater divergence than the parental species of diploid hybrids found at secondary contact zones; (iii) allotriploid taxa possess two conspecific genomes and a deeply diverged allospecific one, suggesting that genomic imbalance and divergence are causal for their partly clonal reproductive mode; (iv) maternal versus paternal genome contributions exhibit asymmetry, with the maternal nuclear (and mitochondrial) genome of polyploids always coming from the same clade, and the paternal genome from the other. We compare our findings with similar patterns in diploid/polyploid vertebrates, and suggest deep ancestral divergence as a precondition for successful allopolyploidization.

Living quarters of a living fossil-Uncovering the current distribution pattern of the rediscovered Hula painted frog (Latonia nigriventer) using environmental DNA

One of the greatest challenges of effective conservation measures is the correct identification of sites where rare and elusive organisms reside. The recently rediscovered Hula painted frog (Latonia nigriventer) has not been seen for many decades and was therefore categorized extinct. Since its rediscovery in 2011, individuals from the critically endangered species have been found, with great effort, only in four restricted sites. We applied the environmental DNA (eDNA) approach to search for new populations of the Hula painted frog in suitable aquatic habitats. We further used the eDNA data to classify the landscape factors associated with the species distribution and to predict its suitable habitats. We sampled 52 aquatic sites in the Hula Valley during the spring of 2015 and 2016 and amplified the samples with a species-specific qPCR assay. DNA of the Hula painted frog was detected in 22 of the sites, all of which clustered within three main areas. A boosting classification model showed that soil type, vegetation cover and the current and former habitats are all key predictors of the frog's current distribution. Intriguingly, the habitat suitability models reveal a high affinity of the species to its long-lost habitat of the historical wetlands. Our findings encourage a series of informed searches for new populations of this threatened frog and provide guidance for future conservation management programmes. In the era of global conservation crisis of amphibians, developing the eDNA approach, a reliable detection method for many critically endangered and elusive amphibians, is of particular importance.

Evolutionarily significant units of the critically endangered leaf frog Pithecopus ayeaye (Anura, Phyllomedusidae) are not effectively preserved by the Brazilian protected areas network

Protected areas (PAs) are essential for biodiversity conservation, but their coverage is considered inefficient for the preservation of all species. Many species are subdivided into evolutionarily significant units (ESUs) and the effectiveness of PAs in protecting them needs to be investigated. We evaluated the usefulness of the Brazilian PAs network in protecting ESUs of the critically endangered Pithecopus ayeaye through ongoing climate change. This species occurs in a threatened mountaintop ecosystem known as campos rupestres. We used multilocus DNA sequences to delimit geographic clusters, which were further validated as ESUs with a coalescent approach. Ecological niche modeling was used to estimate spatial changes in ESUs' potential distributions, and a gap analysis was carried out to evaluate the effectiveness of the Brazilian PAs network to protect P. ayeaye in the face of climate changes. We tested the niche overlap between ESUs to gain insights for potential management alternatives for the species. Pithecopus ayeaye contains at least three ESUs isolated in distinct mountain regions, and one of them is not protected by any PA. There are no climatic niche differences between the units, and only 4% of the suitable potential area of the species is protected in present and future projections. The current PAs are not effective in preserving the intraspecific diversity of P. ayeaye in its present and future range distributions. The genetic structure of P. ayeaye could represent a typical pattern in campos rupestres endemics, which should be considered for evaluating its conservation status.