Amphibians are the focus of a recent debate and public attention owing to the global decline in their populations worldwide. Amphibians are one of the most threatened and poorly known groups of vertebrates in several geographic areas, even though they play a central role in their own ecosystems. At different levels, amphibians make their contribution to proper ecosystem functioning. They act as regulators of the food web and nutrient cycling, and they also provide several valuable ecosystem services, e.g., as a food source and as animal models for lab research. In this sense, it seems clear that the maintenance of amphibian diversity should be one of the major goals for the several countries where their population decline is observed. However, we are still struggling with the very first step of this process, i.e., the correct identification of the amphibian species diversity. Over the past few decades, research on molecular identification of amphibians using DNA barcoding has encountered some difficulties related to high variability in the mitochondrial genome of amphibians, and a research gap is noticeable in the literature. We herein evaluated both COI and 16S rRNA mitochondrial genes for the molecular identification of frogs and tadpoles in a large fragment of the South American Atlantic Rainforest in Rio de Janeiro, Brazil. Our results suggest that both COI and 16S rRNA are informative markers for the molecular identification of the amphibian specimens with all specimens unambiguously identified at the species level. We also made publicly available 12 new sequences of Atlantic Rainforest amphibian species for the first time, and we discussed some conservation issues related to amphibians within the Atlantic Rainforest domains in the state of Rio de Janeiro, Brazil.
Citation: Amaral CRL, Chaves ACS, Borges Júnior VNT, Pereira F, Silva BM, Silva DA, et al. (2019) Amphibians on the hotspot: Molecular biology and conservation in the South American Atlantic Rainforest. PLoS ONE 14(10): e0224320. https://doi.org/10.1371/journal.pone.0224320
Editor: Pedro Abellán, Universidad de Sevilla, SPAIN
Received: May 3, 2019; Accepted: October 11, 2019; Published: October 23, 2019
Copyright: © 2019 Amaral et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the manuscript and its Supporting Information files.
Funding: The authors benefited from grants provided to CFDR (304791/2010-5; 470265/2010-8 and 302974/2015-6) from Conselho Nacional do Desenvolvimento Científico e Tecnológico (CNPq) and through “Cientistas do Nosso Estado” Program from FAPERJ to CFDR (process No. E-26/102.765.2012 E-26/202.920.2015). We would like to thank C. Haddad (Universidade Estadual Paulista) and M. L. Lyra (Universidade Estadual Paulista) for some of the specimens and the primers used in the present manuscript and L. A. Fusinatto, who kindly revised the manuscript offering helpful suggestions. FP was supported by the Portuguese Foundation for Science and Technology (FCT) [IF/01356/2012] and Northern Regional Operational Programme (NORTE2020) through the European Regional Development Fund (ERDF) [MARINFO NORTE-01-0145-FEDER-000031]. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Amphibians are the focus of a recent debate and public attention owing to the global decline in their populations worldwide. The Class Amphibia has an estimated 8,004 species  with Brazil home to the greatest known species richness, having about 1080 species, of which 1039 are anurans [2–3]. Of these, 543 inhabit the South American Atlantic Rainforest biome in Brazil .
One of the world’s most threatened biomes, the South American Atlantic Rainforest originally spread from northeastern Argentina and Paraguay to northeastern Brazil, nowadays has only about 12% of its original cover  and is considered one of the world’s biodiversity hotspots . Among all Brazilian states, about 20% of Rio de Janeiro is covered by Atlantic Rainforest remnants . Located between 20–24° S and 45–41°W, it features a heterogeneous landscape which favors the maintenance of high rates of biodiversity and endemism for several groups [8–9] such as Amphibians [9, 10–13].
In the state of Rio de Janeiro, as well as Brazil as a whole, intense industrial and agroindustrial activities are putting exponentially growing pressure on biodiversity . Oil- and gas-related activities, the expansion of agriculture and livestock areas, the changes caused by urban development, and, surely, the effects of climate change are seriously changing the patterns of land use and deeply influencing biotic communities all around the world [15–17].
Among these communities, amphibians are one of the most threatened and poorly known groups of vertebrates in several geographic areas [18–19], even though they play a central role in their own ecosystems. At different levels, amphibians make their contribution to proper ecosystem functioning. They act as regulators of the food web and nutrient cycling [20–21], and they also provide several valuable ecosystem services, e.g., as a food source and as animal models for lab research [21–23]. In this sense, it seems clear that the maintenance of amphibian diversity should be one of the major goals for the several countries where their population decline is observed. However, we are still struggling with the very first step of this process, i.e., the correct identification of the amphibian species diversity.
Most amphibian species have a complex life cycle, with aquatic larval and terrestrial adult stages, being susceptible to disturbances in both environments [24–26]. Therefore, amphibians are commonly used as a biological indicator of habitat degradation and environmental health status in both aquatic and terrestrial habitats [18,27–31].
Several studies on anuran have explored the variation and characteristics of the postmetamorphic stages, with tadpole morphological characters and ecology comparatively less studied [32–33]. However, the tadpole is the stage most frequently found and easiest registered in nature for several anuran species, by remaining in the breeding spot for longer periods than adults [34–36]. Although the number of studies on tadpoles has increased [33,37–38], the tadpoles of the majority of tropical anuran species remain unknown . Amphibians have been declining from the planet at an alarming rate, and related to the spread of pathogens into wild (including in protected areas) or due to unknown causes [18,33,40–42]. In this sense, conservation and management actions regarding anurans must involve information on the different life stages to ensure efficient and comprehensive conservation plans for the species .
A major drawback for the use of tadpole data in conservation programs is the difficulties of appropriated species identification in such early developmental stages. Molecular methods offer a straightforward alternative for the identification of several organisms such as anurans, including its eggs and larval stages. The DNA barcoding methodology uses a short, standardized gene fragment for species identification [43–44] being a rapid and low-cost method . This method has been already used in the identification of many taxa, such as amphibians [31,42,46], birds [47–48], bivalves , butterflies [50–51], spiders , fishes [53–54], and mammals . In the case of amphibians, a useful application of DNA barcoding is the identification of larvae of some species, especially of those species whose tadpoles are considerably distinct morphologically from their adult phase [31,35,42]. However, sequence variation at the COI priming sites leads to problems in the use of universal primers in amphibian barcoding [31,42,56]. A mixture of several primers is needed to reliably amplify the COI gene from all amphibian species. The use of alternative markers with more conserved priming sites, such as 16S rRNA is often used for identifying amphibian species [31,57]. This calls for the collection of field data and species inventories to guide management and conservation actions and policies all around the world, including Brazil . However, despite the increasing number of amphibian and/or reptile studies for various localities in the South American Atlantic Rainforest, species lists for various localities are still needed in order to provide a comprehensive characterization of their herpetofaunas .
Here, we present a comprehensive identification of adult frogs and tadpoles in the REGUA area. Aimed at precise identification of REGUA amphibian fauna, we lay the groundwork for proper management strategies for this area and the remaining Atlantic Forest fragments. To ensure this precision, we have used DNA barcoding with partial COI and 16S rRNA sequences of adult frogs and tadpoles.
Material and methods
Tissue samples from 88 specimens, including adults (n = 50 specimens from 23 species) and tadpoles (n = 38 specimens from 12 species), were collected in an Atlantic Rainforest fragment in Cachoeiras de Macacu (REGUA), Rio de Janeiro, Brazil (Fig 1). A complete list of all species grouped into its eight distinct families is presented in Table 1. The Reserva Ecológica de Guapiaçu (REGUA—22°24’S, 42°44’W) comprises an area of about 7,600 ha. REGUA is part of a much larger continuous and protected forest area in Rio de Janeiro, Brazil. Together with other conservation units, such as the Parque Estadual dos Três Picos (ca. 46,000 ha), Parque Nacional da Serra dos Órgãos (ca.11,800 ha), and Estação Ecológica Estadual do Paraíso (ca.5,000 ha), REGUA is part of a large contiguous protected area of wet Atlantic Forest, mostly represented by montane and low montane rainforest.
Map of Rio de Janeiro State, Brazil, showing the Reserva Ecológica de Guapiaçú (REGUA) from which specimens were sampled. Atlantic Rainforest remaining fragments are represented in dark and light green; continental waters in dark blue; urbanized area in gray. The map was made using the QGIS software with shapefiles from the IBGE Database (www.ibge.gov.br).
The adults were morphologically identified and classified into 23 species. The tadpoles were cleared and stained for morphological observation and classification (Fig 2). All tissue samples (muscle) were extracted from dead specimens collected with permission under the approved authorization number 18684 issued by ICMBio (Instituto Chico Mendes de Conservação da Biodiversidade). All techniques used to capture, sample, and euthanize were performed to minimize animal suffering and followed the guidelines provided by the Herpetological Animal Care and Use Committee (HACC) of the American Society of Ichthyologists and Herpetologists. The collected individuals were euthanized with an anesthetic application (5% Lidocaine) over the skin. The voucher specimens were fixed in formalin and deposited on the reference collection of the Ecology Department. The samples were preserved in alcohol 70%.
Lateral view (a) and oral disk (b) of 1. Aplastodiscus arildae; 2. Bokermannohyla circumdata; 3. Scinax albicans; 4. Scinax flavoguttatus; 5. Crossodactylus aeneus; 6. Hylodes asper; 7. Hylodes charadranaetes; 8. Hylodes lateristrigatus; 9. Hylodes pipilans; 10. Physalaemus signifer; 11. Proceratophrys appendiculata; and 12. Proceratophrys boiei.
For molecular identification, small fragments of liver and leg muscle were removed from adult amphibians. In the case of tadpoles, fragments of tail and muscle were collected. All samples were submitted to the standard protocol for DNA extraction and purification from the Qiagen QIAmp® DNA FFPE Tissue Kit. Amplification and sequencing of 16S rDNA gene fragments used the following primers: 16SrA-L (5’-CGC CTG TTT ATC AAA AAC AT-3’) and 16SrB-H (5’ -CCG GTC TGA ACT CAG ATC ACG T-3’) . PCR amplifications were performed in a 20 μL volume reaction containing 1 μL (2 ng) of genomic DNA, 12 μL of ddH2O, 0.5 μM of 4X Platinum buffer (Thermo Fisher), 10 μM of each primer and 0.2 μL (1U) of Taq Platinum. The PCR program consisted of 2 minutes at 94°C, 35 cycles of 30 seconds at 94°, 40 seconds at 52°C, and 1 minute at 72°C, with a final extension step for 10 minutes at 72°C. The COI 5’ end gene region was amplified using the following primers: AnF1 (5’-HAA YCA YAA AGA YAT YGG-3’) and AnR1 (5’-CCR AAR AAT CAR AAD ARR TGT TG-3’)  Frog samples for Adenomera marmoratus, Phyllomedusa burmeisteri and Thoropa miliaris were provided by donation by M. Lyra from the Universidade Estadual Paulista (UNESP) and individual information about the specimens could be checked on . PCR amplifications were performed in a 25μL volume reaction containing 2 μL (4 ng) of genomic DNA, 16.5 μL of ddH2O, 0.2 μM of 4X Platinum buffer (Thermo Fisher), 10 μM of each primer and 0.2 μL (1U) of Taq Platinum. The PCR profiles consisted of 3 minutes at 95°C, 5 cycles of 30 seconds at 95°, 30 seconds at 48°C, and 1 minute at 60°C; 30 cycles of 30 seconds at 95°, 30 seconds at 50°C, and 1 minute at 60°C with a final extension step for 5 minutes at 60°C, as previously described for the 16S. The sequencing of the COI gene used M13 extensions tails M13F (5’-TGT AAA ACG ACG GCC AGT-3’) and M13R (5’CAG GAA ACA GCT ATG AC-3’) to maximize the sequencing reliability on the first bases of the COI gene.
The sequencing reactions were performed with BigDye® Terminator v3.1 Cycle Sequencing kit (Thermo Fisher) with 10 seconds at 95°C, 25 cycles of 5 minutes at 50°C, and 4 minutes at 60°C. Sequencing products were processed in an ABI 3500 capillary system (Thermo Fisher). To ensure the absence of saturation in both datasets, both databases (COI and 16S rRNA) were evaluated with and without the 3rd codon position. The neighbor-joining (NJ) and maximum-likelihood (ML) trees were constructed using MEGA 7.0 software . The NJ sequence divergences were calculated based on the Kimura 2-parameter (K2P) distance model with bootstrap support analysis (1,000 replicates), while the maximum-likelihood (ML) result was also generated with bootstrap support analysis (1,000 replicates) with the GTR+G+I nucleotide substitution model chosen as the best fit model based on the AKAYKE (AIC) criterion.
The 50 adult frogs sampled for this study represent 23 distinct species, which is about 31% of the 73 amphibian species registered for the area of REGUA . Tadpoles were represented by 38 individuals. The degenerate PCR primers successfully amplified and sequenced the COI gene region for all 88 analyzed specimens, suggesting that it is a uniquely useful primer pair for the molecular identification of amphibians. No saturation signals on the constructed datasets were observed after evaluation of resultant topologies with and without the 3rd codon position. GenBank accession numbers of the retrieved COI and 16S rRNA reference sequences are available in Table 2. Field and voucher data are available in the S1 File.
New publicly available amphibian sequences
Among all the retrieved sequences of the 23 species found within the REGUA area in the present manuscript, 12 represent new sequences from species that were, for the first time, deposited and made publicly available. They are the sequences for the following species: Crossodactylus aeneus (KU234706), Cycloramphus brasiliensis (KU234692), Dendropsophus bipunctatus (KU234698), Dendropsophus meridianus (KU234699), Hylodes charadranaetes (KU234708), Hylodes lateristrigatus (KU234709), Hypsiboas secedens (KU234700), Phyllomedusa burmeisteri (KU234702), Physalaemus signifier (KU234712), Scinax albicans (KU234703), Scinax flavoguttatus (KU234704), and Scinax humilis (KU234705). The sequences for Hylodes pipilans (KU234710) and Megaelosia goeldii (KU234711) are only the second available molecular record for these two species. The databases, including all sequences used in the present manuscript, are available in the S2 File.
The collected data corresponds to eight amphibian families, as follows: Bufonidae, Cycloramphidae, Hylidae, Hylodidae, Leiuperidae, Leptodactylidae, Odontophrynidae, and Phyllomedusidae. All exhibit both adult and tadpole specimens except for the Cycloramphidae, the Leiuperidae, and the Leptodactylidae. All Neighbor-Joining (K2P) trees can be seen in the S3 File. All divergence values in this section were based on the COI K2P distances, and the datasets used for distance calculations are also available upon request. Only nodes presenting over 50% of bootstrap support values were presented and discussed in the following text.
The Bufonidae was represented by one unique species, Rhinella ornata. The obtained sequences were grouped with previously available sequences of R. ornata and R. crucifer in a clade exhibiting 0.8% of mean distance and with 9.3% divergence of the closest neighbor, R. icterica.
The Cycloramphidae was represented by two species, Cycloramphus brasiliensis and Thoropa miliaris. C. brasiliensis was clearly recovered as a monophyletic clade and diverged 20.2% from its closest neighbor, C. bandeirensis. Thoropa miliaris was also recovered and grouped as a monophyletic clade although it exhibits 8.3% of divergence within the species and 16.9% of divergence from its closest neighbor, T. taophora.
The Hylidae was by far the most represented family within our dataset. It was represented by 10 species, including Aplastodiscus arildae, Dendropsophus berthalutzae, D. bipunctatus, D. meridianus, Hypsiboas secedens, H. semilineatus, Scinax albicans, S. flavoguttatus, and S. humilis. Initially regarded as Aplastodiscus eugenioi based on morphological identification, the REGUA specimens of Aplastodiscus were recovered and grouped with A. arildae, exhibiting 12.4% of divergence from previously published sequences of A. arildae. The closest neighbors were A. albofrenatus with 12.9% of divergence and A. eugenioi, which exhibited 13.9% of divergence. The genus Dendropsophus was represented by three species. D. berthalutzae grouped with the unique sequence publicly available for the species although with a divergence of 18.9%. The remaining two species were D. meridianus and D. bipunctatus, both newly available sequences for the species. D. meridianus was recovered in a node together with D. vraemi in half of the analyses, presenting a divergence of 20.2%, while D. bipunctatus was recovered in a node with low bootstrap support value. The genus Hypsiboas was represented by two species, H. secedens and H. semilineatus. The obtained sequence for H. secedens was the first one available for this species. It was also recovered in a low support node; however, H. semilineatus was clearly recovered together with previously available sequences of the same species with 5.3% of divergence The group formed by both collected and previously available sequences of H. semilineatus was recovered in a group with the species H. boans, diverging 13.5%, and H. wavrini, which presented 13.4% of divergence. The last Hylidae genus in our dataset, Scinax, was represented by three species, including Scinax albicans, S. flavoguttatus, and S. humilis. Although recovered clearly grouped, both S. flavoguttatus and S. humilis are newly sequenced species, and reference sequences are not available for comparison. Scinax albicans was recovered as monophyletic and grouped to its closest neighbor, S. obtriangulatus, exhibiting a distance of 17.3%.
The Hylodidae was represented by five species grouped into two genera, Crossodactylus and Hylodes. Crossodactylus aeneus sequences are new, and the unique available sequences for comparison came from C. caramaschii, which diverged from C. aeneus by 19.1%. C. aeneus also exhibited within divergence of 3.7%. The genus Hylodes was respresented by four species, including Hylodes asper and H. pipilans with available sequences for comparison and the new H. charadranaetes and H. lateristrigatus. Hylodes asper grouped with available H. asper sequences although they exhibited 9.2% of divergence from them, while H. pipilans also grouped with available H. pipilans sequences, but presented only 6.2% of intraspecific divergence.
The Leiuperidae was represented uniquely by Physalaemus signifier which was recovered in a group that includes the basal P. nattereri, and, sequentially, P. signifier, P. bokermanni, P. atlanticus, and P. moreirae. P. signifier exhibited divergence of 15.5% from P. bokermanni, its closest neighbor within the group.
The Leptodactylidae was represented by two species, Adenomera marmorata and Megaelosia goeldii. Adenomera marmorata was recovered and grouped with available sequences of A. marmorata and presented A. ajurauna as the closest neighbor, exhibiting 13.3% of divergence from it and 1.3% of intraspecific divergence, while A. ajurauna exhibited 5.9% of divergence within the species. Megaelosia goeldi was also recovered and grouped with previously available M. goeldi sequences although it exhibited divergence of 7.2% from the unique previously published sequence.
The family Odontophrynidae was represented by Proceratophrys boiei and P. appendiculata. Both species were successfully recovered and grouped with previously available sequences for the species, diverging 7.1% and 2.3% from them, respectively.
The unique Phyllomedusidae Phyllomedusa burmeisteri is newly sequenced, and no previous sequences are available for comparison. It was recovered together with P. distincta, exhibiting 8.9% of divergence.
Finally, within distances among the several species herein surveyed ranged from quite small distances, such as 0.1% found within Scinax albicans sequences, to 3.2% − 3.7% found for Proceratophrys boiei and Crossodactylus aeneus and 7.5% − 8.3% found within Aplastodiscus arildae and Thoropa miliaris, respectively.
The topologies observed for both COI and 16S rRNA were comparable and supported by high bootstrap values. All tadpoles were unambiguously identified at species level based on the reference database built with the adult frogs. The COI and 16S NJ trees, including both tadpoles and specimens from the reference database (adult frogs), show that all species were recovered as monophyletic groups with high bootstrap support (Figs 3 and 4). The overall COI K2P mean distance is 0.245, ranging from 0.187 between Proceratophrys boiei and Proceratophrys appendiculata to 0.340 between Scinax flavoguttatus and Rhinella ornata.
Bootstrap probabilities near each internal node. Tadpoles are marked with a T after the species epithet. All remaining specimens comprise adults from the reference database.
Bootstrap probabilities near each internal node. Tadpoles are marked with a T after the species epithet. All remaining specimens comprise adults from the reference database.
As an alternative molecular marker, 16S rRNA barcoding correctly identified all species in the present manuscript. The neighbor-joining (NJ) and maximum-likelihood (ML) resultant trees were presented in Fig 4A and 4B, and they show that all individuals were recovered as monophyletic group with high bootstrap support. For both NJ and ML trees, specimen vnbj22 Aplatodiscus arildae was recovered outside of the group formed by the remaining individuals of the same species. The 16S rRNA interspecific distances were also calculated for all species. The values ranged from 0.039 between Proceratophrys boiei and P. appendiculata to 0.189 between Hylodes asper and Aplastodiscus arildae.
Discussion and conclusions
The state of Rio de Janeiro is home to about 40% of the South American Atlantic Rainforest amphibian species. About 19% is regarded as data deficient, meaning that the conservation status of these species has not yet been evaluated and that the number of threatened species within the South American Atlantic Rainforest could be higher.
Together with the anthropogenic pressures related to habitat fragmentation and destruction, the presence of exotic invasive species also threatens survivability of amphibian fauna within the South American Atlantic Rainforest. Van Sluys et al.  reported the expansion of the exotic American bullfrog (Lithobates catesbeianus) within the state of Rio de Janeiro, pointing to the need for population monitoring and management programs for removal of individuals from natural habitats and, thereby, favoring the maintenance of native species and populations.
Recently,  presented the last checklist of amphibian species from the Atlantic Rainforest in Rio de Janeiro, indicating that the amphibian richness presented is still underestimated, tending to increase in the following years as a result of new taxonomic studies and local herpetofaunal surveys. The authors also pointed out the importance of Rio de Janeiro state as a key spot for amphibian diversity, not only for the Atlantic Rainforest, but for Brazilian amphibian fauna as a whole.
Amphibians are commonly regarded as a challenging group for molecular identification approaches and the usefulness of these methods within amphibians has already been discussed in several manuscripts over the last decades [31,64–66]. Difficulties related to the high variability of their mitochondrial genome and the genetic structure of their recognized species make identification of amphibian species based on the well- accepted DNA barcoding methodology a kind of outlier among herpetologists. The DNA barcoding sequence, or the first 650bp of the Cytochrome C oxidase subunit I gene, as proposed by , is highly polymorphic among amphibians, and its use for species identification depends on a mixture of various primers to correctly amplify all species, mainly from the high level of polymorphisms observed on the common COI priming sites. Difficulties could be also observed when assessing the population level of amphibian diversity where the analysis and comparison of mitochondrial and nuclear datasets is rarely concordant, and events of occasional hybridization and introgression commonly affect amphibian species and their population substructure. In this sense, it is clear that species delimitation among amphibians does not seem to be an easy task, especially using only the COI barcoding approach. It is well accepted that the use of alternative markers, such as the more preserved 16S rRNA, should be taken into account for amphibians, the species of which are relatively old entities and exhibit a high amount of mutations .
Our results for both COI and 16S rRNA (Figs 3 and 4) showed that at least within the limits of the obtained databases, all species were clearly recovered as monophyletic units, therefore making the DNA barcoding approach herein used a successful methodology for species diversity assessment within the analyzed fragment of the South American Atlantic Rainforest. As expected, the interspecific divergence was consistent, ranging from 0.187 to 0.340, while the intraspecific divergence was tiny for almost all the species, ranging from 0.002 to 0.007, except for Hylodes charadranaetes, which presented an intraspecific divergence of 0.042. However, independent of the result observed for H. charadranaetes, the divergence observed was sufficient to ensure a useful barcoding gap among all the analyzed species.
The recent study of  also reached high levels of success using COI barcoding on the molecular identification of amphibian species from west Central Africa, as well as those already obtained by [46,61].
Vences et al.  also pointed out the tadpole / adult correlation as one of the major goals for using molecular identification methodologies. The similarity observed for several species included in the present survey (Fig 2) is remarkable, even with stained sections, such as the buccal apparatus, a morphological complex commonly used for species determination on tadpoles. The identification of tadpoles is surely one of the major challenges in studying amphibians by virtue of being the most fragile life stage of anurans [67–68]. The patterns of convergence and parallelism among relatively closely related taxa constitute a particular challenge to the study of tadpoles in the field . The tadpoles of many anuran species have not been morphologically described [39,70]. Furthermore, tadpoles of related species can be externally similar in morphology. In such cases, genetic analyses have been considered a successful tool to identify tadpoles at the species level [31,70]. This success was surely reached herein since all the tadpoles were correctly assigned to their species. Tadpole molecular monitoring is a positive step for determining the amphibian species richness of a given area. They are the most common and easily found form for several amphibian species within the South American Atlantic Rainforest. They are the most persistent stage, staying in the same area/pond for several days and probably weeks, making them an easy target for sampling.
Finally, we do believe, as already observed by [46,61], that the DNA barcoding methodology could be clearly used for amphibian species identification and richness estimation, therefore guiding the development of appropriate conservation and management strategies. We also agree, as already mentioned by , that it is far more advantageous to include amphibians in a DNA barcoding initiative rather than excluding them. The BOLD database is a huge achievement for several zoological groups all around the world, but amphibians were left by the wayside in recent years because of a supposed inability to use the technology for species identification. However, it seems obvious that correct taxonomic identification depends on the comprehensiveness of the dataset built for comparative purposes, and this concern is surely not restricted to amphibians. In this sense and regarding all threats amphibians are facing nowadays, the construction of reliable databases is urgent, and we should focus on the construction of publicly available datasets of reference sequences that include geographical and ecological information, as well as any other kind of data that support proper species identification, in the effort to support ecologists and conservation biologists as they struggle to develop strategies of amphibian fauna management and conservation within the South American Atlantic Rainforest.
S2 File. Databases used in the present manuscript.
S3 File. NJ trees.
NJ K2P trees for each studied taxon.
S4 File. COI distances.
COI K2P distances for the analyzed groups.
This study is a portion of the results of the “Programa de Pesquisas em Biodiversidade da Mata Atlântica (PPBio Mata Atlântica Program)" of Ministério do Meio Ambiente do Brasil (MMA) and was supported by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) (Process No. 457458/212-7).
- 1. Amphibiaweb (2019) Amphibiaweb: information on amphibian biology and conservation. Berkeley, California. Amphibiaweb—http://amphibiaweb.org/. (Accessed: April 19, 2019).
- 2. Segalla MV, Caramaschi U, Cruz CAG, Grant T, Haddad CFB, Garcia PCA, Berneck BVM, Langone JA (2016) Brazilian amphibians: list of species. Herpetologia Brasileira, 5, 34–46.
- 3. SBH (2010) Brazilian amphibians–list of species. http://www.sbherpetologia.org.br. (accessed: April 19, 2019).
- 4. Haddad CFB, Toledo LF, Prado CPA, Loebmann D, Gasparini L, Sazima I (2013) Guia dos anfíbios da Mata Atlantica–Diversidade e Biologia. São Paulo, Anolis Books.
- 5. Ribeiro MC, Metzger JP, Martensen AC, Ponzoni FJ, Hirota MM (2009) The Brazilian Atlantic Forest: How much is left, and how is the remaining forest distributed? Implication for conservation. Biological Conservation, 142, 1144–1156.
- 6. Mittermeier RA, Gil RP, Hoffman M, Pilgrim J, Brooks T, Mittermeier CG, Lamoreux J, Fonseca GAB (2005) Hotspots revisited: earth’s biologically richest and most endangered terrestrial ecoregions. Chicago, University of Chicago Press.
- 7. Bergallo HG, Uzêda MC, Fidalgo ECC, Alves MAS, Rocha CFD, van Sluys M, Correa e Castro TC, Bebianno MC, Santos MA & Cozzolino ACR (2009) Conservação da biodiversidade da Mata Atlântica no estado do Rio de Janeiro: uma nova abordagem. In: Bergallo HG, Uzêda MC, Fidalgo ECC, Alves MAS, Rocha CFD, van Sluys M, Correa e Castro TC, Bebianno MC, Santos MA & Cozzolino ACR (eds.) Estratégias e ações para a conservação da biodiversidade no Estado do Rio de Janeiro. Rio de Janeiro, Instituto Biomas, 23–32.
- 8. Jenkins CN, Pimm SL (2006) Definindo prioridades de conservação em um hotspot de biodiversidade global. In: Rocha CFD, Bergallo HG, van Sluys M, Alves MAS (eds.) Biologia da conservação: essências. São Carlos, RiMa. 41–52.
- 9. Nascimento JL, Campos IB (2011) Atlas da fauna brasileira ameaçada de extinção em unidades de conservação federais. Brasília, ICMBIO.
- 10. Rocha CFD, Bergallo HG, Alves MAS, van Sluys M (2005) Endemic and threatened tetrapods in the restingas of the biodiversity corridors of Serra do Mar and of the central of Mata Atlântica in eastern Brazil. Brazilian Journal of Biology, 65, 159–168.
- 11. Cruz CAG & Feio RN (2007) Endemismos em anfíbios em áreas de altitude na Mata Atlântica no sudeste do Brasil. In: Nascimento LB & Oliveira ME (eds.) Herpetologia no Brasil II. Belo Horizonte, Sociedade Brasileira de Herpetologia, 117–126.
- 12. van Sluys M, Cruz CAG, Vrcibradic D, Silva HR, Almeida-Gomes M, Rocha CFD (2009) Anfíbios nos remanescentes florestais de Mata Atlântica no Estado do Rio de Janeiro. In: Bergallo HG, Fidalgo ECC, Rocha CFD, Uzêda MC, Bebianno MC, Alves MAS, van Sluys M, Santos MA, Correa e Castro TC, Cozzolino ACR (eds.) Estratégias e ações para a conservação da biodiversidade no Estado do Rio de Janeiro. Rio de Janeiro, Instituto Biomas, 175–182.
- 13. Dorigo TA, Vrcbradic D, Rocha CFD (2018) The amphibians of the state of Rio de Janeiro, Brazil: an updated and commented list. Papéis Avulsos de Zoologia, 58, e20185805.
- 14. Wilkie D, Shaw E, Rotberg F, Morelli G, Auzel P (2000) Roads, development, and conservation in the Congo Basin. Conservation Biology, 14, 1614–1622.
- 15. Gibbs HK, Ruesch AS, Achard F, Clayton MK, Holmgren P, Ramankutty N, et al. (2010) Tropical forests were the primary sources of new agricultural land in the 1980s and 1990s. Proceedings of the National Academy of Sciences of the Unites States of America, 107, 16732–16737. pmid:20807750
- 16. Osti M, Coad L, Fisher J, Bomhard B, Hutton J (2011) Oil and gas development in the World Heritage and wider protected area network in sub-Saharan Africa. Biodiversity and Conservation, 20, 1863–1877.
- 17. Hof C, Araujo MB, Jetz W, Rahbek C (2011) Additive threats from pathogens, climate and land-use change for global amphibian diversity. Nature, 480, 516–519. pmid:22089134
- 18. Stuart SN, Chanson JS, Cox NA, Young BE, Rodrigues ASL, Fischman DL, Waller RW (2004) Status and trends of amphibian declines and extinctions worldwide. Science, 306, 1783–1786. pmid:15486254
- 19. IUCN (2015) The IUCN Red List of Threatened Species <http://www.iucnredlist.org>. Downloaded on 19 November2015.
- 20. Cortés-Gomez AM, Ruiz-Agudelo CA, Valencia-Aguilar A, Ladle RJ (2015) Ecological functions of neotropical amphibians and reptiles: a review. Universitas Scientiarum, 20, 229–245.
- 21. Hocking DJ, Babbitt KJ (2014) Amphibian contributions to ecosystem services. Herpetological Conservation and Biology, 9, 1–17.
- 22. Warkentin IG, Bickford D, Sodhi NS, Bradshaw CJA (2009) Eating frogs to extinction. Conservation Biology, 23, 1056–1059. pmid:19210303
- 23. Deichmann JL, Mulcahy DG, Vanthomme H, Tobi E, Wynn AH, Zimkus BM, et al. (2017) How many species and under what names? Using DNA barcoding and GenBank data for west Central African amphibian conservation. PLoSONE, 12, e0187283. pmid:29131846
- 24. Duellman WE, Trueb L (1994) Biology of Amphibians. Baltimore and London. McGraw-Hill.
- 25. Becker CG, Fonseca CR, Haddad CFB, Batista RF, Prado PI (2007) Habitat split and the global decline of amphibians. Science, 318, 1775–1777. pmid:18079402
- 26. Almeida-Gomes M, Rocha CFD (2014) Landscape connectivity may explain anuran species distribution in an Atlantic forest fragmented area. Landscape Ecology, 29, 29–40.
- 27. Welsh HH, Ollivier LM (1998) Stream amphibians as indicators of ecosystem stress: a case study from California’s redwoods. Ecological Applications, 8, 1118–1132.
- 28. Alford RA, Richards SJ (1999) Global amphibian declines: a problem in applied ecology. Annual Review of Ecology and Systematics. 30, 133–165.
- 29. Welsh HH, Droege S (2001) A case for using plethodontid salamanders for monitoring biodiversity and ecosystem integrity of North American forests. Conservation Biology, 15, 558–569.
- 30. Davic RD, Welsh HH (2004) On the ecological role of salamanders. Annual Review of Ecology and Systematics, 35, 405–434.
- 31. Vences M, Van der Meijden A, Chiari Y, Vieites DR (2005a) Comparative performance of the 16S rRNA gene in DNA barcoding of amphibians. Frontiers in Zoology, 2.
- 32. Larson PM, de Sá RO (1998) Chondrocranial morphology of Leptodactylus larvae (Leptodactylidae: Leptodactylinae): its utility in phylogenetic reconstruction. Journal of Morpholology, 238, 287–305.
- 33. Provete DB, Garey MV, Da Silva FR, Jordani MX (2012) Knowledge gaps and bibliographical revision about descriptions of free-swimming anuran larvae from Brazil. North-Western Journal of Zoology, 8, 283–286.
- 34. Lips KR, Savage JM (1996) Key to known Tadpoles (Amphibia: Anura) of Costa Rica. Studies on Neotropical Fauna and Environments, 31, 17–26.
- 35. Altig R, McDiarmid RW (1999) Diversity: Familial and Generic Characterizations. In Tadpoles. In: McDiarmid RW & Altig R (eds.) The Biology of Anuran Larvae. University of Chicago Press, Chicago and London, 295–337.
- 36. Rossa-Feres DC, Nomura F (2006) Characterization and taxonomic key for tadpoles (Amphibia: Anura) from the northwestern region of São Paulo State, Brazil. Biota Neotropica, 5.
- 37. McDiarmid RW, Altig R (1999) Tadpoles: The biology of Anuran Larvae. University of Chicago Press.
- 38. Natale GS, Alcalde L, Herrera R, et al. (2011) Underwater acoustic communication in the macrophagic carnivorous larvae of Ceratophrys ornata (Anura: Ceratophryidae). Acta Zoologica, 92, 46–53.
- 39. Laia RC, Rocha CFD(2012). Adults and tadpoles of species of Hylodidae (Anura): History and taxonomy perspectives. Zoologia, 29, 89–94.
- 40. Berger L, Speare R, Daszak P, Green DE, Cunningham AA, Goggin CL, Slocombe R, Ragan MA, Hyatt AD, McDonald KR, Hines HB, Lips KR, Marantelli G, Parkes H. (1998) Chytridiomycosis causes amphibian mortality associated with population declines in the rainforests of Australia and Central America. Proceedings of the National Academy of Sciences of the United States of America, 95, 9031–9036. pmid:9671799
- 41. Daszak P, Cunningham AA, Hyatt AD (2003) Infectious disease and amphibian population declines. Diversity and Distributions, 9, 141–150.
- 42. Vences M, Thomas M, Bonett RM, Vieites DR (2005b) Deciphering amphibian diversity through DNA barcoding: Chances and challenges. Philosophical Transactions of the Royal Society B: Biological Sciences, 360, 1859–1868.
- 43. Hebert PDN, Cywinska A, Ball SL, deWaard JR (2003) Biological identifications through DNA barcodes. Proceedings of the Royal Society B Biological Sciences, 270, 313–321.
- 44. Hebert PDN, Gregory TR (2005) The promise of DNA barcoding for taxonomy. Systematic Biology, 54, 852–859. pmid:16243770
- 45. Golding RDH, Hebert PDN (2009) Preface. Molecular Ecology Resources, 9(Suppl.1), iv–vi.
- 46. Smith MA, Poyarkov NA Jr, Hebert PD (2007) DNA barcoding amphibians: take the chance, meet the challenge. Molecular Ecology Notes, 8, 235–46.
- 47. Hebert PDN, Stoeckle MY, Zemlak TS, Francis CM (2004) Identification of birds through DNA barcodes. PLoS Biology, 2, 1657–1663.
- 48. Kerr KCR, Lijtmaer DA, Barreira AS, Hebert PDN, Tubaro PL (2009) Probing evolutionary patterns in Neotropical birds through DNA barcodes. PloSONE 4, e4379.
- 49. Jämergren J, Schander C, Sneli JA, Ronningen V, Young CM (2007) Four genes, morphology and ecology: distinguishing a new species of Acesta (Mollusca; Bivalvia) from the Gulf of Mexico. Marine Biology, 152, 43–55.
- 50. Janzen DH, Hajibabaei M, Burns JM (2005) Wedding biodiversity inventory of a large and complex Lepidoptera fauna with DNA barcoding. Philosophical Transactions of the Royal Society B: Biological Sciences, 270, 96–99.
- 51. Silva-Brandão KL, Lyra ML, Freitas AVL (2009) Barcoding lepidoptera: current situation and perspectives on the usefulness of a contentious technique. Neotropical Entomology, 38, 441–451. pmid:19768260
- 52. Barret RDH, Hebert PDN (2005) Identifying spiders through DNA barcodes. Canadian Journal of Zoology, 83, 481–491.
- 53. Ward RD, Zemlak TS, Innes BH, Last PR, Hebert PDN (2005) DNA barcoding Australia’s fish species. Philosophical Transactions of the Royal Society B: Biological Sciences, 360, 1847–1857.
- 54. Amaral CR, Brito PM, Silva DA, Carvalho EF (2013) A new cryptic species of South American freshwater pufferfish of the genus Colomesus (Tetraodontidae), based on both morphology and DNA data. PloS One, 8, e74397. pmid:24040239
- 55. Borisenko AB, Lim BK, Ivanova NV, Hanner RH, Hebert PDN (2008) DNA barcoding in surveys of small mammal communities: a field study in Suriname. Molecular Ecology Resources, 8, 471–479. pmid:21585824
- 56. Vieites DR, Wollenbergb KC, Andreone F, Kohler J, Glaw F, Vences M (2009) Vast underestimation of Madagascar’s biodiversity evidenced by an integrative amphibian inventory. Proceedings of the National Academy of Sciences of the Unites States of America, 106, 8267–8272.
- 57. Vences M, Gehara M, Köhler J, Glaw F (2012) Description of a new Malagasy treefrog (Boophis) occurring syntopically with its sister species, and a plea for studies on non-allopatric speciation in tropical amphibians. Amphibia-Reptilia, 33, 503–520.
- 58. Verdade VK, Valdujo PH, Carnaval AC, Schiesari L, Toledo LF, et al. (2012) A leap further: the Brazilian amphibian conservation action plan. Alytes, 29, 8–43.
- 59. Almeida-Gomes M, Siqueira CC, Borges-Junior VNT, Vrcibradic D, Fusinatto LA, Rocha CFD (2014) Herpetofauna of the Reserva Ecológica de Guapiaçu (REGUA) and its surrounding areas, in the state of Rio de Janeiro, Brazil. Biota Neotropica, 14(3), e20130078.
- 60. Palumbi SR (1996) Nucleic acids II: the polymerase chain reaction. Molecular Systematics, 205–247.
- 61. Lyra ML, Hadadd CFB, Azeredo-Espin AML (2017) Meeting the challenge of DNA barcoding Neotropical amphibians: polymerase chain reaction optimization and new COI primers. Molecular Ecology, 17, 966–980. pmid:28029226
- 62. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S (2011) MEGA: Molecular evolutionary genetics analysis using Maximum Likelihood, evolutionary distance, and Maximum Parsimony methods. Molecular Biology and Evolution, 10, 1093.
- 63. van Sluys M, Cruz CAG, Vrcibradic D, Silva HR, Almeida-Gomes M, Rocha CFD (2009) Anfíbios nos remanescentes florestais de Mata Atlântica no Estado do Rio de Janeiro. In: Bergallo HG, Fidalgo ECC, Rocha CFD, Uzêda MC, Bebianno MC, Alves MAS, van Sluys M, Santos MA, Correa e Castro TC, Cozzolino ACR (eds.) Estratégias e ações para a conservação da biodiversidade no Estado do Rio de Janeiro. Rio de Janeiro, Instituto Biomas, 175–182.
- 64. Arntzen JW (1989) The identification of eggs of some species of palaearctic newts (genus Triturus) by starch gel electrophoresis. Isozyme Bulletin. 22, 69.
- 65. Veith M, Kosuch J, Feldmann R, Martens H, Seitz A (2000) A test for correct species declaration of frog legs imports from Indonesia into the European Union. Biodiversity Conservation, 9, 333–341.
- 66. Wilson AC, Maxson LR, Sarich VM (1974) Two types of molecular evolution. Evidence from studies of interspecific hybridization. Proceedings of the National Academy of Sciences of the Unites States of America, 71, 2843–2847.
- 67. Heyer WR, McDiarmid RW, Weigmann DL (1975) Tadpoles, predation, and pond habitats in the tropics. Biotropica, 7, 100–111.
- 68. Fatorelli PC & Rocha CFD (2008) O que molda a distribuição das guildas de girinos tropicais? Quarenta anos de busca por padrões. Oecologia brasiliensis, 12, 733–742.
- 69. Wiens JJ, & Graham CH (2005) Niche conservatism: integrating evolution, ecology, and conservation biology. Annual Review of Ecology, Evolution, and Systematics, 36, 519–539.
- 70. Thomas M, Raharivololoniaina L, Glaw F, Vences M, Vieites D.R (2005) Montane tadpoles in Madagascar: molecular identification and description of the larval stages of Mantidactylus elegans, Mantidactylus madecassus, and Boophis laurenti from the Andringitra Massif. Copeia, 1, 74–183.