Next Article in Journal
Intramuscular Fatty Acids in Meat Could Predict Enteric Methane Production by Fattening Lambs
Next Article in Special Issue
Amphibian and Reptile Road Mortality in Special Nature Reserve Obedska Bara, Serbia
Previous Article in Journal
SARS-CoV-2 Clinical Outcome in Domestic and Wild Cats: A Systematic Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Update on Distribution and Conservation Status of Amphibians in the Democratic People’s Republic of Korea: Conclusions Based on Field Surveys, Environmental Modelling, Molecular Analyses and Call Properties

by
Amaël Borzée
1,2,*,
Spartak N. Litvinchuk
3,4,
Kyongsim Ri
5,
Desiree Andersen
6,
Tu Yong Nam
7,
Gwang Hyok Jon
8,
Ho Song Man
9,
Jong Sik Choe
9,
Sera Kwon
6,
Siti N. Othman
6,
Kevin Messenger
10,
Yoonhyuk Bae
1,6,
Yucheol Shin
1,
Ajoung Kim
6,
Irina Maslova
11,
Jennifer Luedtke
2,12,
Louise Hobin
2,
Nial Moores
13,
Bernhard Seliger
14,
Felix Glenk
14 and
Yikweon Jang
15
add Show full author list remove Hide full author list
1
Laboratory of Animal Behaviour and Conservation, College of Biology and the Environment, Nanjing Forestry University, 159 Longpan Rd, Nanjing 210037, China
2
Amphibian Specialist Group, IUCN Species Survival Commission, Toronto, ON L5A, Canada
3
Institute of Cytology, Russian Academy of Sciences, Tikhoretsky pr. 4, 194064 St. Petersburg, Russia
4
Department of Zoology and Physiology, Dagestan State University, Gadzhiyev str. 43-a, Makhachkala, 3367000 Dagestan, Russia
5
Department of International Economic Cooperation, Ministry of Land and Environment Protection, Pyongyang, Democratic People’s Republic of Korea
6
Interdisciplinary Program of Eco Creative, Ewha Womans University, Seoul 03760, Korea
7
Institute of Zoology, State Academy of Science, Daesong-dong, Daesong District, Pyongyang, Democratic People’s Republic of Korea
8
Department of Ecology, State Academy of Science, Daesong-dong, Daesong District, Pyongyang, Democratic People’s Republic of Korea
9
Department of Ecology, Life Science College, Kim Il Sung University, Ryongnam-dong, Daesong-dong, Daesong District, Pyongyang, Democratic People’s Republic of Korea
10
Herpetology and Applied Conservation Lab, College of Biology and the Environment, Nanjing Forestry University, 159 Longpan Rd, Nanjing 210037, China
11
Federal Scientific Center of the East Asia Terrestrial Biodiversity, Far Eastern Branch of Russian Academy of Sciences, 690022 Vladivostok, Russia
12
Re:wild, Austin, TX 78746, USA
13
Birds Korea, 101-1902, Hyundai I Park, Busan 48559, Korea
14
Hanns Seidel Foundation, Seoul 04419, Korea
15
Department of Life Sciences and Division of EcoScience, Ewha Womans University, Seoul 03760, Korea
*
Author to whom correspondence should be addressed.
Animals 2021, 11(7), 2057; https://doi.org/10.3390/ani11072057
Submission received: 26 April 2021 / Revised: 25 May 2021 / Accepted: 1 July 2021 / Published: 9 July 2021
(This article belongs to the Special Issue Conserving Reptiles and Amphibians)

Abstract

:

Simple Summary

In this study, we used field and literature surveys, call recordings, photographs, landscape models and molecular tools to estimate the presence, range and status of amphibians in the Democratic People’s Republic of Korea. We found 18 native species and the suspected presence of Karsenia koreana and two Onychodactylus species. We also determined northern range boundaries for Rana uenoi and Dryophytes japonicus with molecular tools. Based on distribution and modelling, we can expect the contact zone between species within the Rana and Onychodactylus genera to be located along the Changbai Massif, a high mountain range. The species richness was higher in the lowlands and at lower latitudes, with up to 11 species, while species richness in northern areas was half that value. Following the categories and criteria of The IUCN Red List of species, ecological models and known threats, we recommended ten species to be listed as threatened at the national level. The ecology of species in the DPR Korea is understudied, although species reliant on agricultural wetlands (e.g., rice paddies) are not as threatened as species living in forested areas due to the enduring presence of extensive agricultural landscapes.

Abstract

Determining the range, status, ecology and behaviour of species from areas where surveys and samplings are uncommon or difficult to conduct is a challenge, such as in the Democratic People’s Republic of Korea (DPR Korea). Here, we used genetic samples, field surveys, call recordings, photographic identification and a literature review to estimate the presence, range and status of amphibians in the DPR Korea. From our combined results and based on the IUCN Red List categories and criteria, we were able to estimate the national threat levels for most species. Our results demonstrated the presence of 18 native species and the suspected presence of Karsenia koreana and two Onychodactylus species. We reported the first record for Rana uenoi in the vicinity of Pyongyang using molecular tools and similarly confirmed the presence of Dryophytes japonicus at the same location. Based on distribution and modelling, we can expect the contact zone between species within the Rana and Onychodactylus genera to be located along the Changbai Massif, a mountain range that marks a shift in ecoregions and acts as a barrier to dispersion. The species richness was higher in the lowlands and at lower latitudes, with such areas populated by up to 11 species, while more northern regions were characterised by species richness of about half of that value. The combination of ecological models and known threats resulted in the recommendation of ten species as threatened at the national level following the IUCN Red List categories and criteria. This high number of threatened species was anticipated based on the high threat level to amphibians in bordering nations and globally. While the ecology of species in the DPR Korea is still understudied, we argue that species relying on agricultural wetlands such as rice paddies are not under imminent threat due to the enduring presence of extensive agricultural landscapes with low rates of chemical use and mechanisation. The maintenance of such landscapes is a clear benefit to amphibian species, in contrast to more industrialised agricultural landscapes in neighbouring nations. In comparison, the status of species dependent on forested habitats is unclear and threat levels are likely to be higher because of deforestation, as in neighbouring nations.

1. Introduction

Delineating the boundaries of species distributions is a preliminary requisite to determining potential presence, taxonomic research and threat levels that may potentially require conservation actions. The distributions of large vertebrates such as mammals and birds are now coarsely defined for most species, but details are missing for a large number of species from other taxonomic groups [1]. For instance, as of March 2019, 21.35% of described amphibians (1443 species) were still classified as data-deficient on the IUCN Red List of threatened species, a reasonable proxy for missing distribution data when in conjunction with a handful of other variables. Many of these amphibian species are from areas where it is difficult to collect data. This is the case for the Democratic People’s Republic of Korea (hereafter, DPR Korea), where recent scientific exploration and research has not been as intense as in many other nations. Some exceptions do exist, however, especially for birds [2,3,4]. Presence points for some amphibian species in DPR Korea are available from non-focal references, such as wetland inventories, but these do not constitute descriptions of distribution ranges by themselves [5,6,7,8].
Fifteen species are listed as present in DPR Korea in the update by Kim and Han [9], and a sixteenth invasive species, Lithobates catesbeianus, is also listed [10]. While most of the species are valid [11], phylogenetic developments have shown that Rana ornativentris [12,13], R. temporaria [14] and R. chensinensis [15,16,17] do not occur in DPR Korea. In addition, Bufo sambangensis has not been studied since its description in 1996, and the status of the species as well as its potential synonymy with B. stejnegeri [18] have not been confirmed through molecular tools [19]. Similarly, Dryophytes suweonensis (following [20]; previously Hyla suweonensis [21]) is missing from the list [22,23], while other genera such as Onychodactylus and Karsenia are expected to be present, although their presence in the nation has not yet been confirmed [24,25,26].
In addition, the absence of data on ranges is especially true for amphibians, as they are generally harder to detect and are sometimes considered to be of lower interest than other groups of species [27,28]. Combined with DPR Korea’s geography, i.e., being flanked by sea both to the east and the west, this means that the distribution of species in the nation is especially difficult to estimate. Species present both to the south and to the north of DPR Korea are expected to be present throughout the nation; however, when a species is present only to the north or the south, the boundaries of such ranges are unknown. For example, the range of Pelophylax chosenicus [29] was expected to extend into DPR Korea, as it is known to occur very close to its southern border [30,31,32]. Accordingly, a recent survey found the species substantially further north than originally expected [11,33].
In addition, the fact that a species occurs both to the north and to the south of DPR Korea does not necessarily make its distribution easier to predict. For instance, D. japonicus [34] is expected to be present in most of the nation [9] as it occurs both to the north and to the south [35] and because the species’ ecological requirements can be met throughout the majority of the area [36]; however, D. japonicus is divided into several clades [37] and the southern distribution limit of individuals assigned to the clade, referred to as the “ussuriensis clade” by some authors [38,39], remains unknown.
A geographically similar situation involving several species arises for the Rana genus where the range of multiple species borders the DPR Korea. While some species are present only to the south (e.g., R. uenoi [40]) or to the north (e.g., R. dybowskii [41,42], R. amurensis [43,44] and R. chensinensis [45]), others such as R. huanrenensis [46,47,48] and R. coreana [45,49] are present both to the north and south. Consequently, some of the species present in the Primorsky Krai in the Russian Federation (hereafter Russia), in the north-eastern People’s Republic of China (hereafter PR China) and along the northern border of the Republic of Korea (hereafter R Korea), may be present in DPR Korea. In this case, the species are likely in contact and potentially hybridising. In theory, and with the understanding that habitat segregation between Rana species is a potential factor for speciation, clarifying the presence of species through molecular tools should be straightforward if genetic samples are available (e.g., [9,50,51,52]).
Here, we combine the use of molecular tools on old museum collections and our field surveys to determine the distributions of all species occurring in the nation and to assess their conservation status where possible. This project relies on several types of data and analyses—first mitochondrial genetic data from historical specimens, followed by photographic identification of field and museum animals, aural and visual field surveys and acoustic analyses. Finally, we assess species diversity, model species ranges and provided a preliminary national conservation status for all amphibian species in DPR Korea. The last steps are conducted in order to prioritise conservation targets, both for species and habitat types, which would benefit a large range of species. We expect several species to be under higher levels of threat, in line with the current amphibian crisis [27,53,54], while other species are expected to be widespread and found in large numbers.

2. Materials and Methods

2.1. Taxonomic Sampling

The genetic samples used in this project were collected at an unspecified location in the vicinity of Pyongyang, DPR Korea (Figure 1). The samples used here are muscle tissues from Rana sp. (n = 13) and Dryophytes sp. (n = 12). The Rana sp. collection was conducted by the North Korea Academy of Sciences (Chosun Academy of Science) and samples (voucher IDs in Table 1) were sent to the ZISP (Zoological Institute of the Russian Academy of Sciences, St. Petersburg, Russia) during “Soviet times”, prior to the 1990’s. The Dryophytes sp. collection comes from two sources. Eight adult specimens have the same origin as the Rana sp. samples. An adult and two juvenile specimens were collected in July 1947 and housed in the collection of the S.M. Kirov Military Medical Academy (St. Petersburg, Russia) until 2018, after which they were transferred by I.V. Doronin to the ZISP collection. Three of the Dryophytes samples were stored in formalin, while all other Rana and Dryophytes samples were stored in alcohol. Although it is unknown how the samples were fixed, we suspect fixation to have been conducted in formalin based on the sampling period. Consequently, fixation and dehydration were unlikely to have been conducted following current recommendations, such as fixation times of 14–24 h and thorough dehydration prior to embedding in order to avoid severe DNA fragmentation.

2.2. Genetic Identification

Genetic identification of individuals was needed when target species did not present discrete morphologies, such as between R. uenoi and R. dybowskii. Misidentifications are especially problematic when estimating distributions and generate inaccurate results [56], potentially overestimating or underestimating the threat levels to species.
Advancements in technologies have made it possible to extract ancient DNA (e.g., [57]) or DNA from biological remnants such as insect exuviae [58]. Similarly, it is possible to retrieve DNA from samples fixed and stored in formalin [59] or other fixation chemicals [60]; however, this task is made increasingly difficult when samples are fixed in an unknown chemical and stored for decades in another chemical with unknown intermediate treatment.
For this reason, we adapted the Qiagen Formalin DNA extraction protocol (QIAamp® DNA FFPE Tissue; Venlo, The Netherlands). Here, not all samples were entirely lysed within 3 h, already 2 h longer than described in step 11 of the manufacturer standard protocol [61]. We avoided longer incubation times or higher incubation temperatures as these may have resulted in more fragmented DNA; thus, we removed the pellet 3 h after the beginning of incubation as described in step 11 and restarted the procedure from step 3 with the pellet extracted, i.e., for the remaining tissue. This multiplication of steps was conducted so that tissues already lysed were not submitted to further degradation from heat. All resulting samples were treated independently.
Once the genomic DNA was extracted and eluted, we used it for barcoding with the primers and protocols described in the literature [42,62,63]. All PCRs were run in duplicate and as gradient PCR, with a two-degree increment framing the recommended annealing temperature in two steps, such as −4 °C, −2 °C, Ta, +2 °C and +4 °C. For all samples, we amplified a fragment of the 12S ribosomal small subunit. For Dryophytes sp. samples, the primer pair and protocol used followed [62], while for Rana sp., the primer pair and protocols were the ones used in [42]. PCR reactions were performed in a SimpliAmp™ Thermal Cycler (Applied Biosystems, Waltham, MA, USA). Products were visualised on 1.5% agarose gel loaded with three microliters of PCR products, run on an Agaro-Power™ System (A-7020, Bioneer, Daejeon, R Korea) and visualised with a Nucleic Acid Bioimaging Instrument Blue Illuminator (S, NeoScience, Suwon, R Korea) using TopGreen Nucleic Acid 6× Loading Dye (GenomicBase, Seoul, R Korea).
Most PCRs did not yield visible bands, but all samples with a single band visible were sent for sequencing. When multiple bands were visible, the band with the target product size only was extracted (FavorPrep Gel/PCR purification mini kit; Favorgen, Taiwan) and sent for sequencing. All samples were sent for purification and both forward and reverse sequencing by Cosmogenetech (Cosmogenetech Co., Ltd., Seoul, R Korea).
The resulting sequences were trimmed and analysed in Geneious v 11.0.2 (Biomatters Limited, Auckland, New Zealand). Only two sequences for Dryophytes and eight sequences for Rana were good enough to be retained. To assess which clades the samples belonged to, we imported georeferenced sequences from GenBank (Table 2), based on the closest match when blasted against the NCBI database and matching sequences for all species listed [9]. The missing georeferences for GenBank samples were extracted from the literature [64,65,66,67,68] for Dryophytes sp. and from [42,48,69,70] for Rana sp.
The sequences were aligned independently for each genera using MUSCLE (v3.8.31; [71]), implemented through the Geneious plug-in, with a maximum of 10 iterations following default parameters and further revised manually when needed. We then reconstructed a maximum likelihood (ML) tree for reach genera, inferred from the 12S fragment and 349 bp long for Dryophytes sp. and 293 bp long for Rana sp. We also added Hyla orientalis as the outgroup for the Dryophytes sp. dataset and R. kukunoris as the outgroup clade for the Rana sp. dataset (GenBank accession numbers shown in Table 2). All analyses were performed using the PHYML plugin in Geneious [72]. This method implements a fast and accurate heuristic estimate of maximum likelihood phylogenies, with GTR selected as the substitution model for flexibility reasons. We performed the run with 50,000 bootstraps under default variables.

2.3. Pictures Identification

Among the individuals that could not be identified by molecular markers, some could be identified through morphological cues. Rana coreana [73] and R. amurensis have warts on their lateral side, which are usually black [44] but sometimes red in R. amurensis during the breeding season [74,75]. In contrast, these warts are absent on Rana samples not pertaining to the R. amurensis group [76], specifically R. chensinensis, R. pirica [69], R. huanrenensis, R. dybowskii and incidentally R. uenoi here (see [40,45]).
For hylids, D. suweonensis is smaller and more slender than D. japonicus, and the angle between the eyes and ipsilateral nostrils can be used for species identification [77]. Although morphologically differentiating traits with D. immaculatus and D. flaviventris are present [78], the species is not present in the area [79] and D. suweonensis is known to be present further north than Pyongyang [22,23].

2.4. Field Surveys

Field surveys to detect species presence were conducted in 2016, 2018 and 2019. The surveys conducted in Mundeok (39.5489° N, 125.4390° E; datum = WGS84) in March 2017 and May 2018 were mostly opportunistic, as they were originally focused on avifauna, but the area was surveyed again in June 2019 with a specific focus on amphibians. Opportunistic surveys were also conducted in Kumya, South Hamgyong (39.5386° N, 127.2206° E) in late May 2018. Finally, amphibian-focused surveys were conducted between 24 and 26 March and between 2 and 6 June 2018 in the area of Rason, within the localities of Rajin and Sonbong, and limited north by the Tumen River, marking the border between DPR Korea and Russia (42.3446° N, 130.4712° E; Figure 1). This coastal area is characterised by high hills in the south and wetlands in the north, including a Ramsar site, the Rason Migratory Bird Reserve. The forested habitat is substantially, degraded except for isolated patches; however, agricultural wetlands are numerous and most of the agricultural work is conducted by hand, providing adequate substitute wetlands for amphibians.
During amphibian-focused surveys, each time a habitat potentially hosting amphibians was encountered along a predetermined route, the site was investigated through aural and visual surveys for amphibian presence. Surveys were conducted through aural and visual detection, following established protocols whenever possible [36,80,81]; however, surveys were required to avoid encounters with local people and protocols were adapted following this and other local requirements. Daytime surveys were conducted under general search-and-encounter protocols, including substrate flipping when allowed. Night-time surveys followed spotlight and call count protocols [36,82]. We also took ecological measurements, whereby we measured the water quality with a PCSTestr 35 multimetre (Oakton Instruments; Vernon Hills, IL, USA) for salinity (ppm), pH, temperature (°C), conductivity (μS) and total dissolved solids (tds; ppm). Several types of landscapes were investigated, including natural wetlands, agricultural wetlands, grassy hills and forested areas.
In addition, snout–vent–length (SVL) measurements were performed down to the nearest 0.1 mm for ten individual D. japonicus samples caught in a rice paddy close to Rajin (42.304538° N, 130.390762° E; callipers 1108-150 W, Insize; Suzhou, China). The SVL values for these 10 individuals were compared through a t-test (SPSS, v21, IBM SPSS Statistics Inc., Chicago, IL, USA) to the SVL values of 308 individuals measured in R Korea, PR China, Mongolia, Russia and Japan (data for individual frogs outside of DPR Korea not provided here) to characterise the variation within the species.

2.5. Call Properties

Five acoustic files between 31 and 63 s in length were obtained from Mundeok (39.5489° N, 125.4390° E) in May 2017. The acoustic files were recorded in MP3 format, which was converted into WAV format using VLC media player v. 3.0.4. The files all had high levels of background noise, which was reduced using Audacity® (v2.2.0). Data were then extracted using R [83] and the acoustic packages WarbleR [84], seewave [85] and tuneR [86]. We then described call properties based on oscillograms and spectrograms and measured call properties such as fundamental and dominant frequencies, before comparing them to calls from other parts of the species range [87].

2.6. Distribution Estimates and Species Richness

To model distribution estimates for all species, we used the datapoints presented above and investigated all accessible literature to assess the presence of all amphibian species expected to be extant in DPR Korea. To date, the most extensive and accurate studies on the subject were the ones conducted by [88,89], which were mostly transcribed by [51] and completed by [9]. The presence points are reported on Figure 1 and GPS coordinates are reported in Supplementary Materials Table S2. All studies citing amphibians from DPR Korea available were analysed, although some references only reported the potential presence of the species or did not provide clear georeferenced datapoints associated with specific species [6,18,29,90,91,92,93,94,95,96]; thus, only clear and georeferenced presence points are listed below. Several presence points were acknowledged by [11] and used as test data for the accuracy of the data collected. In addition, datapoints available from the North Korean Human Geography website (http://www.cybernk.net (accessed on 1 July 2019)) were directly extracted from, meaning the reference was not added to the list below. Bufo sambangi [19], also spelled B. sambangensis [9,10], was not included in this analysis due to the paucity of information supporting the existence of the species and morphological similarities with B. stejnegeri (presented by [9]; Figure 2). Additional datapoints were provided from past field observations during a workshop in Pyongyang in June 2019, but only 11 were retained as they were within 20 km of a confirmed presence point.
Furthermore, all georeferenced mentions of R. uenoi, R. dybowskii, R. temporaria, R. chensinensis, R. ornativentris or any of the combination above using either species name as species or subspecies were assigned to the corresponding clade based on the recent nomenclature when available, or were not used in models when no consensus could be reached (although the GPS data are presented in Supplementary Table S2 under the label Rana sp.). Morphological cues to discriminate between R. uenoi and R. dybowskii have not yet been described and the range of the two species is expected to be continuous, while the location of the contact zone between the two species has still not been determined. Rana huanrenensis is not included in any studies from DPR Korea, but data from the R Korea and PR China were used to model the range of the species. Finally, we considered Onychodactylus to be present in DPR Korea based on [24], and considered O. zhaoermii as the species with the largest range in northern regions and O. zhangyapingi as potentially present; however, neither of the species has been recorded under those names in the nation yet. The datapoints used in the models are listed below.

2.6.1. Anura

Bufo gargarizans [97]: [9,50,76,89,91,98,99]; Bufo stejnegeri: [5,100,101]; Strauchbufo raddei [102]: [7,51,89]; Dryophytes japonicus: [5,9,50,51,91,98,103,104]; Dryophytes suweonensis: [22,23]; Bombina orientalis [105]: [9,50,51,52,73,91,98,100,106,107,108,109,110]; Kaloula borealis [111]: [9,50]; Rana amurensis: ZISP collection, [73]; Rana coreana and R. amurensis: [9,50,73,107,112,113,114,115]; Rana uenoi and R. dybowskii: [50,51,73,91,107] (using identification cues from [116] and integrating the work from [88,89]) [5,9,117]; Glandirana emeljanovi [118]: [9,50,51,73,91]; Pelophylax nigromaculatus [119]: [50,51,73,98,99,100,107,120]); Pelophylax chosenicus: [9,33].

2.6.2. Caudata

Hynobius leechii: [9,50,99,100,121,122]; Onychodactylus koreanus: [50,51] (originating from [89]) [5,7]; Salamandrella tridactyla: [51] (originating from [89]) [9,112,123].
It should be noted that three references of potential interest, as cited in other work, could not be traced, namely [124,125]. These were not included in our models, despite the potential presence of valuable datapoints. For our models, we also used the datapoints from GBIF.org (doi.org/10.15468/dl.qetyb9 (accessed on 26 February 2019)), filtered for amphibians from DPR Korea, R Korea, PR China and Russia. The non-focal species of this selection were filtered out, and in total we used 24,193 data points for all species confounded (Supplementary Table S1 for all datapoints for all species; Supplementary Table S2 for datapoints within DPR Korea). Historical datapoints from [126] were not used for the models. Additional fine tuning for specific species was also conducted.
Rana coreana: The northernmost point in Ryanggang [107] was excluded as the morphological cues to discriminate between R. coreana and R. amurensis had not yet been determined at that time.
Rana uenoi and R. dybowskii: The cut-off point used for the range of these two species was the latitude matching with the locality closest to Pyongyang, based on the genetic analysis conducted in this study.
Rana huanrenensis: We added georeferenced datapoints from the literature [48,127,128] to avoid north–south bias in the number of points.
Onychodactylus sp.: Due to the unknown range boundaries of all four species potentially present in DPR Korea and the absence of data, a single model was created, including all records for the genus (including additional datapoints from [24,25]). Only O. fischeri has been listed from DPR Korea so far, despite its presence being unlikely or limited to the extreme north-east of the nation. As O. koreanus is the only species for which sufficient confirmed datapoints were available for an independent species model, we ran two models, one for Onychodactylus sp. and one for O. koreanus.
Salamandrella tridactyla: Additional datapoints from [25,129] were used for the models.
Karsenia koreana: This species has not been found in the nation yet and was only recently described [130]; it may be present in the south-east of DPR Korea [26]. Due to the paucity of information, we did not include this species in the list of modelled species.
Lithobates catesbeianus: The presence of the species has been confirmed [10] in DPR Korea, but we did not include models for the species as it is invasive.
We used maximum entropy modelling (MaxEnt; [131]) of 19 bioclimatic factors (Hijmans et al., 2005: percent contribution, permutation importance and training gain for each bioclimatic variables in Supplementary Materials Tables S3–S5, respectively) to delineate potential ranges of each amphibian species. We decided to include all variables despite the potential correlations, as ecological requirements are not known for the species selected here and they are bound to differ greatly. A subselection of variables would result in less accurate models for some of the species, as well as an inability to compare variable responses among species; therefore, to avoid the exclusion of relevant variables [132], but also to create a baseline for future studies, we decided to include all bioclimatic variables and to ensure the absence of preconceived bias on the relevant ecological variables. For each species and model, MaxEnt was trained for ten bootstrap replicates with 20% random test percentage and then projected to 30 arc second (0.00833 dd) rasters of the Korean Peninsula. To balance the sampling bias generated by more numerous datapoints in R Korea compared to DPR Korea, we constructed two bias layers, namely 1 (DPR Korea):10 (R Korea) and 1 (DPR Korea):5 (R Korea)—equivalent to the distribution of collected unique occurrence points. Additionally, duplicate presence records (presences within the same grid cells of the environmental layers) were removed. For each species, the model selected was based on the highest true skill statistic (TSS) of the 10-percentile training omission threshold (Table 3). We modified models further for R. amurensis and R. huanrenensis; the models were trained with 0.02 dd rasters of Northeast Asia (without bias layer) and projected to 0.00833 dd rasters of the Korean Peninsula to remedy the insufficient presence of data points in DPR Korea. Predicted presence was determined with the following thresholds: minimum presence, 95% presence, 10-percentile training omission. The 10-percentile training omission threshold was further used for TSS reporting [133]. The total surface areas of the minimum presence and 10-percentile training omission polygons were calculated as each species’ expected land surface area. We also compiled a species richness map for the all species merged together for DPR Korea, using a sum of all presence probabilities for species distribution models, where richness values were classified as 0 (not suitable/present), 0.25 (minimum presence), 0.5 (95%) and 1.0 (10-percentile training omission).

2.7. Extinction Risk Assessment

One of the most common and most reliable ways to conduct assessments of extinction risk is through the use of the IUCN Red List categories and criteria [134]. These have been used to inform and catalyse the conservation of numerous species [135,136]. The data on the IUCN Red List are used in conjunction with other tools to determine conservation priorities (e.g., [137,138]), the identification and delineation of protected areas (e.g., [139]) and as an indicator of the health of the planet’s biodiversity [140]. The IUCN Red List categories and criteria are applied at the global scale [141], but they have also been adapted for national [142] and regional scales [143,144].
IUCN requires that a species be evaluated against the quantitative thresholds of the five criteria that determine the risk of extinction. When a species meets one or more of the criteria, it is assigned to one of the three threatened categories based on the threshold that is met: vulnerable, endangered or critically endangered. The criteria are: A, population size reduction; B, geographic range size; C, small population size and decline; D, very small population or restricted distribution; E, quantitative analysis of extinction risk. A species that does not meet the threshold for any of these criteria is placed into one of the other IUCN categories; however, as very few animal populations are monitored regularly and in their totality [134], protocols have been established for small datasets, which can encompass a broad type of data and account for varying degrees of uncertainty in the event that no robust data are available [145,146].
Although some authors have provided estimates of population sizes (Table 4), and others have provided data on the type of threat to be considered for specific species, population dynamics in DPR Korea are mostly unknown. We conducted a workshop in Pyongyang in June 2019, during which 12 herpetologists with knowledge on species within DPR Korea were asked to provide population trends for the species they knew about (Table 4). From our own surveys, no exact numbers can be reported due to restrictions in field surveys, but the assessment is a comparison with known population sizes in PR China, Russia and R Korea.
Dryophytes japonicus: Population surveys in Rason, Kumya and Mundeok (Figure 1) recorded large population sizes at all sites.
Dryophytes suweonensis: Population surveys in Mundeok (Figure 1) recorded a large population size.
Pelophylax nigromaculatus: Population surveys in Rason, Kumya and Mundeok (Figure 1) recorded large population sizes at all sites.
Pelophylax chosenicus: Population surveys in Mundeok (Figure 1) recorded a large population size.
When population size estimates were not available, extinction risk evaluations were conducted according to estimates of geographic range size (criterion B). The IUCN Red List categories and criteria include two different measures applied to the geographic range. One is the extent of occurrence (EOO), defined as the area contained within the shortest continuous imaginary boundary that can encompass all the known, inferred or projected sites of present occurrence of a taxon, excluding cases of vagrancy, and is a measure of the spread of extinction risk across a species’ range. This requires the delineation of a species’ range, which we modelled in the section above. The second is the area of occurrence (AOO), defined as the area within the EOO that is occupied by the species, which cannot be used here as it requires fairly comprehensive field surveys to generate sufficient occurrence data, which are then assigned to a 2 × 2 km grid cell. Process-based modelling using small-scale environmental variables (e.g., microclimate) can be applied to refine the species distribution models and better determine the possible geographic ranges of species. As an example, here the range of D. suweonensis can be modelled based on datapoints available both within and outside of DPR Korea (see species section) and the distribution model can be further restricted based on the maximal elevation known for the species [36]; therefore, we used the result of the MaxEnt models to assess the extent of occurrence (EOO) and assess threat levels following IUCN Red List thresholds and classification schemes, following the guidelines and examples from [134] and threat definitions [147]. The preliminary extinction estimates presented here are guidelines and include all species present and likely to be present in DPR Korea, set against the Red List Regional Guidelines.

3. Results

3.1. Genetic Identification

Out of the 12 Dryophytes sp. samples from which DNA was extracted, lysis was incomplete within the original 3 h for all but four samples. We transferred the four samples to a new lysis buffer and they were totally lysed on the second round. This additional step resulted in a total of 16 samples for which DNA was extracted. The DNA concentration was <0.2 ng/µL for all resulting samples, but we used all of the samples as templates for subsequent PCRs. PCRs and sequencing of adequate quality for the analysis conducted here were successful for only two samples (Supplementary Table S6) out of the 160 PCRs conducted (16 samples ran in duplicate on gradient PCR with five temperature increments). The DNA concentration obtained was slightly higher for Rana sp. in comparison with Dryophytes sp. Out of 13 samples, seven were duplicated for a second round of lysis, resulting in 20 samples with DNA concentrations <0.5 ng/µL, leading to 200 PCR reactions. Counting samples from both genera, seven samples from the initial lysis were successfully sequenced, while three were successfully sequenced from the secondary lysis (Supplementary Table S6). Successful sequencing arose from a single PCR at the recommended annealing temperature, while five samples resulted in successful sequencing for recommended Ta + 2 °C and four samples resulted in successful sequencing for recommended Ta + 4 °C (Supplementary Table S6). All sequences presented here have unresolved calls in sequencing that do not prevent species identification but would not be adequate for analyses further than the one presented. Consequently, sequences were not submitted to GenBank due to further need for protein annotation but are presented in Supplementary Table S6.
The reconstructed ML tree showed that the two Dryophytes samples clustered with D. japonicus (ZISP.14002 and 14004; Figure 3). The analysis for Rana sp. showed that all samples collected were clustered with R. uenoi (ZISP.13968–13981; Figure 4) and are the first report of the species in the vicinity of Pyongyang and in DPR Korea. For this analysis (see Figure 4), the species identification for Rana dybowskii individuals from R Korea was arbitrarily reassigned to R. uenoi based on [40]. The analysis presented here is adequate in that it depicts R. dybowskii and R. pirica as sister clades, and highlights the fact that none of the sample for which DNA was extracted belong to the species named in [9].

3.2. Picture Identification

From the collection at the ZISP (number 14036) and based on available morphological cues (Figure 5), we identified one individual as belonging to either R. amurensis or R. coreana, the first record in the vicinity of Pyongyang for either species. In addition to the presence of warts on the lateral sides, the examined individual had a non-interrupted white upper lip (Figure 5A), unlike other Northeast Asian Rana sp. We also identified individuals as D. suweonensis (ZISP.14019–14021; Figure 5B), Kaloula borealis (ZISP.14038–14044) and Pelophylax nigromaculatus (ZISP.14022–14033), the first records for these species in the vicinity of Pyongyang.

3.3. Field Surveys

Aural field surveys in Mundeok in 2017 resulted in the detection of D. japonicus, D. suweonensis, P. nigromaculatus, P. chosenicus and Bombina orientalis. The same species, with the exception of B. orientalis, were detected in 2018 and 2019. Surveys in Kumya in 2018 only detected P. nigromaculatus and D. japonicus.
The surveys conducted in Rason in March 2018 resulted in the detection of a single Rana sp. egg clutch, which could not be closely investigated and was, thus, assigned to Rana sp. The absence of detection of other species does not indicate their absence and may be related to other variables such as weather. During the surveys, the average air temperature was 13.42 ± 2.42 °C (mean ± SD), the average air relative humidity was 43.26 ± 6.01% and the average water temperature was 8.48 ± 4.07 °C (min = 2.10 °C). Based on surveys conducted in R Korea, P.R. China and Russia, the air and water temperatures and air humidity were adequate for spawning by early-breeding species such as Rana sp., Hynobius sp., Salamandrella sp. and Onychodactylus sp. The temperatures were slightly lower than that at which B. gargarizans and Strauchbufo raddei usually breed and too low for other amphibian species potentially present in the area to be active. The other environmental variables collected, namely salinity, pH, conductivity and total dissolved solids, were also within the values accepted by the species suspected to be present [148,149,150]. The only variable with slightly lower values than otherwise expected was total dissolved solids (120.87 ± 74.01 ppm), likely due to recently melted snow.
The surveys in Rason in June 2018 resulted in the detection of numerous species, including all anuran species expected to be present in the area based on models described below, with the exception of caudata, in a wide variety of wetlands (Figure 6). The water quality was not generally different to that used by the same species at similar latitudes elsewhere. At the coastal site where S. raddei, P. nigromaculatus and B. orientalis were found, the variables were (14:30 hh:mm) air temperature = 30.6 °C, relative humidity = 32.9%, air pressure = 1007.8 hPa, water temperature = 29.2 °C, conductivity = 125.4 µS, pH = 8.04, salinity = 171 ppm and tds = 343 ppm. At one of the agricultural wetlands where D. japonicus and P. nigromaculatus were recorded, the variables were (20:20) air temperature = 17.5 °C, relative humidity = 71.9%, air pressure = 1008.3 hPa, water temperature = 22.1 °C, conductivity = 203 µS, pH = 8.42, salinity = 145 ppm, and tds = 102 ppm. On the grassy hills where B. orientalis was found to breed in small water holes, the variables were (17:00) air temperature = 18.4 °C, relative humidity = 69.3%, air pressure = 1009.4 hPa, water temperature = 31.4 °C, conductivity = 146 µS, pH = 8.91, salinity = 43.4 ppm and tds = 69.3 ppm. In the natural wetland where B. orientalis was found (11:00), the variables were air temperature = 27.4 °C, relative humidity = 71.0%, air pressure = 1007.8 hPa, water temperature = 29.2 °C, conductivity = 85.54 µS, pH = 8.89, salinity = 49.3 ppm and tds = 62.7 ppm. Finally, at the natural wetlands where D. japonicus tadpoles were found, the environmental variables were (15:40) air temperature = 24.2 °C; relative humidity = 71%, air pressure = 1007.8 hPa, water temperature = 29.2 °C, conductivity = 85.5 µS, pH = 8.89, salinity = 49.3 ppm and tds = 62.7 ppm. GPS coordinates for all sites are provided in Supplementary Table S2 and are visualised in Figure 6.
Dryophytes japonicus was the most commonly detected species at night via mating calls. Despite some instances of individuals calling before sunset and as early as midday, the large majority of individuals called after sunset and choruses were characterised by a calling index equal to three for all populations surveyed after sunset [36]. Tadpoles were also detected in a variety of water bodies, although not at a Gosner stage >25 [104,151], likely due to the early season. The species was widespread in all types of wetlands and present in large numbers in agricultural wetlands. When comparing SVL values of individuals caught in Rason (n = 10) to those from neighbouring countries (n = 308), the differences were not significant (t-test; p = 0.063). The SVL value for D. japonicus in Rason was 34.62 ± 4.73 cm (mean ± SD).
Pelophylax nigromaculatus was the second most abundant species based on call surveys and the most abundant species based on daytime encounter surveys. Tadpoles were, however, not present in large numbers, indicating a delayed breeding season compared to D. japonicus. This species was also widespread in all types of wetlands and was present in large numbers in agricultural wetlands.
Bombina orientalis was detected at three different sites in densities comparable to those of populations found in R Korea and Russia (n < 10). The species was detected through it spawn at one site, best described as grassy hills with tire tracks on a non-asphalted road, which were filled with water. The species was present in small and shallow water bodies but was not found in agricultural wetlands.
Strauchbufo raddei tadpoles were detected in a shallow wetland in a sandy landscape close to the sea. Bufo gargarizans was detected once only, in a mountainous area during a short survey, and consequently expected to be locally abundant in this habitat. Rana dybowskii was also detected once in the same area as B. gargarizans, and the species is consequently equally expected to be locally abundant in this habitat.

3.4. Call Properties

The manual inspection of the recording files from Mundeok revealed the presence of D. suweonensis, P. nigromaculatus and P. chosenicus. Analyses through R managed to extract data for D. suweonensis using automated detection, resulting in 89 separate call extractions from the three files the species was present in. Frequency analysis showed a mean frequency of 3.27 kHz (Figure 7). Automated and manual detection in R failed for the other two species. The call properties of the species were not different from those in the southern parts of its range [87], although the variation could not be statistically assessed due to the low sample size.

3.5. Distribution Estimates and Species Richness

The results of the MaxEnt modelling highlighted a higher amphibian biodiversity and species richness in the southwest of the nation and along the northern section of the coast of the Yellow Sea (number of species ≥10; Figure 8). The lowest species richness was found in the northeast, close to the border with PR China and the Changbai and Baekdu mountains (n ≤ 5), while it was higher along the east coast (5 ≤ n ≤ 8). Generally, lower altitude and latitude were associated with higher species richness and vice versa.
Individual species distributions generally followed the same pattern (Figure 9); however, some species were more widely distributed than others, such as B. orientalis, B. gargarizans, D. japonicus, G. emeljanovi, P. nigromaculatus and R. huanrenensis. These species are also the ones with the largest described range in Northeast Asia. Other species were limited to the southwestern areas, also matching with lowlands (D. suweonensis, H. leechii, K. borealis and R. coreana). Bufo stejnegeri, R. uenoi and O. koreanus were remarkable due to their modelled presence in the southeast of the nation only. Salamandrella tridactyla was modelled as present in the northeastern region only or in high-elevation habitats. Rana dybowskii followed the same pattern, although reaching lower latitudes. The model for Onychodactylus sp. showed that the species are likely abundant in DPR Korea, although species assignment is so far unknown. The distribution models for P. chosenicus and S. raddei provided larger ranges than expected when compared to their distributions in adjacent nations [152], likely due to missing data, while the model for R. amurensis was very likely imperfect for the same reasons, as well as other unknown ecological variables, such as potential exclusive competition with R. coreana.
Based on the distribution models, we estimated the surface area based on the minimum suitable area and the 10-percentile omission area (Table 5). When comparing the top five species with the broadest estimated range for the two estimates, only B. orientalis, D. japonicus and P. nigromaculatus were overlapping.

3.6. Extinction Risk Assessment

As a general rule, biodiversity on the Korean Peninsula has declined over the past 200 years [153] and conservation actions are urgently required [154]. As of 2013, 6.3% of terrestrial land area is protected in DPR Korea, a high number when compared to some neighbouring countries, but a figure lower than the worldwide average of 12% [155]. General threats to ecosystems in DPR Korea have been shown to be mostly linked to land use [153,156], with 21% of land used for agricultural purposes and 46% preserved as forests [155]. Some areas were, however, still largely deforested at the beginning of the millennium [157], despite strong conservation policies being recently implemented for reforestation [156], maintenance of good water quality [157] and conservation and mitigation of environmental impacts [158]. On the other hand, threats because of invasive species may be lower due to trading policies [153]. It is also important to note that numerous protected areas have been designated, some under international conventions such as Ramsar and the UNESCO Man and the Biosphere Programme [153].
Threat levels to amphibians are higher than for other animals in DPR Korea [153], in line with the global pattern [27,53,54,159]. While no species is known to be endemic to DPR Korea, among the species present in the nation, D. suweonensis, P. chosenicus, O. koreanus, K. koreana and R. uenoi are so far described as endemic to the Korean Peninsula. Among the species listed in this study, only D. suweonensis and P. chosenicus are listed in one of the threatened categories by the IUCN, as endangered [160] and vulnerable [152], respectively. Here, we list all threats impacting the species independently of their severity. The global conservation status categories are the ones currently available and are listed for reference only. The results of our analyses for conservation status within DPR Korea are discussed below and shown in Table 5.

3.6.1. Bombina orientalis

This species is expected to be among the most abundant amphibian species. The species is widespread both throughout the nation and mainland Northeast Asia, with apparently large and abundant populations in a large variety of habitats. In addition, no parasites were found on the seven individuals tested [98]. The principal threat is habitat modification. Recommendation for IUCN National Red List assessment (national level): least concern. Current IUCN Global Red List assessment: least concern [161].

3.6.2. Bufo gargarizans

Despite once thought to have been introduced [29], the presence of the species in Northeast Asia demonstrates the native status of the species. The species is widespread both throughout the nation and mainland East Asia, with apparently large and abundant populations. Parasites have been found for this species but are not assessed as threatening, including cestode or nematode parasites on six individuals [98]. The principal threat is habitat modification. Recommendation for IUCN National Red List assessment (national level): least concern. Current IUCN Global Red List assessment: least concern [162].

3.6.3. Bufo stejnegeri

The species is vulnerable under criterion B1 on the national red list, as its EOO within the nation meets the IUCN threshold, occurring in fewer than 10 threat-defined locations, with estimated ongoing decline in its: (i) extent of occurrence; (ii) area of occupancy; (iii) area, extent or quality of habitat, potentially resulting in local extirpations. The species is restricted to high elevations and is generally found in small populations. The principal threats are habitat modification, environmental pollution and climate change. Recommendation for IUCN National Red List assessment (national level): vulnerable B1ab(i,ii,iii). Current IUCN Global Red List assessment: least concern [163].

3.6.4. Dryophytes japonicus

The species is expected to be one of the most abundant amphibian species nationwide. The species is widespread, with apparently large and abundant populations in a large variety of habitats. In addition, no parasites were found on the only individual tested [98]. The principal threats are habitat modification and environmental pollution, but climate change is known not to be a threat [164]. Recommendation for IUCN National Red List assessment (national level): least concern. Current IUCN Global Red List assessment: least concern [35].

3.6.5. Dryophytes suweonensis

The species is Endangered under criterion B1, as it is present within four independent locations and its EOO meets the IUCN threshold due to estimated ongoing decline in its: (i) extent of occurrence; (ii) area of occupancy; (iii) area, extent or quality of habitat, potentially resulting in local extirpations. In addition, as the species is restricted to lowlands <120 m a.s.l, its range is severely fragmented. The principal threats are habitat modification [80], environmental pollution [30,165] and hybridisation [166]. There is also an increased extinction risk due to the presence of the invasive Lithobates catesbeianus [10], which is abundant in habitats similar to those used by the species [167] and which plays a substantial role in pathogen transmission [168]. Recommendation for IUCN National Red List assessment (national level): endangered B1ab(i,ii,ii). Current IUCN Global Red List assessment: endangered [169].

3.6.6. Glandirana emeljanovi

The species is widespread, with apparently large and abundant populations in a large variety of habitats. The principal threat is habitat modification. Recommendation for IUCN National Red List assessment (national level): least concern. Current IUCN Global Red List assessment: least concern [170].

3.6.7. Hynobius leechii

The species is widespread, with apparently large and abundant populations restricted to mid-elevation habitats. The principal threat is habitat modification. Recommendation for IUCN National Red List assessment (national level): least concern. Current IUCN Global Red List assessment: least concern [171].

3.6.8. Kaloula borealis

The species is relatively widespread across the nation and in mainland Northeast Asia, although to a lesser extent than other species in the nation. The species is also susceptible to habitat degradation, which could potentially cause range contractions in the near future bringing the national EOO within the vulnerable threshold under B1 and could meet the other threshold subcriteria. The principal threats are habitat modification, chemical pollution and exploitation. Recommendation for IUCN National Red List assessment (national level): near threatened. Current IUCN Global Red List assessment: least concern [172].

3.6.9. Onychodactylus sp.

It is not possible to assess the conservation status of species in this genus as it represents a group of several species within unknown boundaries; however, they seem relatively widespread, although restricted to high elevation habitats, which are under threat of habitat modification. Consequently, the principal threats are habitat modification, environmental pollution and climate change. Three of the species are absent from Table 5 as their presence is not yet confirmed in DPR Korea.

3.6.10. Onychodactylus fischeri

Unconfirmed presence and predicted restricted distribution in the nation if present. Threats are likely to include habitat loss and climate change. Recommendation for IUCN National Red List assessment (national level): data-deficient. Current IUCN Global Red List assessment: least concern [173].

3.6.11. Onychodactylus zhaoermii

Threats are likely to include habitat loss and climate change. Recommendation for IUCN National Red List assessment (national level): data-deficient. Current IUCN Global Red List assessment: not evaluated.

3.6.12. Onychodactylus zhangyapingi

Threats are likely to include habitat loss and climate change. Recommendation for IUCN National Red List assessment (national level): data-deficient. Current IUCN Global Red List assessment: not evaluated.

3.6.13. Onychodactylus koreanus

The species is vulnerable under criterion B1 on the national red list as its national EOO meets the IUCN threshold, it occurs in fewer than 10 threat-defined locations and due to estimated ongoing decline in its: (i) extent of occurrence; (ii) area of occupancy; (iii) area, extent or quality of habitat, potentially resulting in local extirpations. The principal threats are habitat modification, environmental pollution and climate change; however, population dynamics are not known for the species and the species could be assessed as EN under the A criterion if further research were to demonstrate that the population is declining, in view of the extremely limited expected range of the species (Table 5). Recommendation for IUCN National Red List assessment (national level): vulnerable B1ab(i,ii,ii). Current IUCN Global Red List assessment: not evaluated.

3.6.14. Pelophylax chosenicus

The species is Vulnerable under criterion B1 on the national red list as its EOO within the nation meets the IUCN threshold, it occurs in fewer than 10 threat-defined locations and due to estimated ongoing decline in its: (i) extent of occurrence; (ii) area of occupancy; (iii) area, extent or quality of habitat, potentially resulting in local extirpations. The species is present in western areas, and while it seems to be thriving at some localities, it is restricted to lowlands and to a single type of habitat: wetlands. This makes the species vulnerable to habitat modification and environmental pollution. In addition, several subpopulations have reportedly declined and become extirpated as a result of droughts, likely linked to climate change. Recommendation for IUCN National Red List assessment (national level): vulnerable B1ab(i,ii,ii). Current IUCN Global Red List assessment: vulnerable [152].

3.6.15. Pelophylax nigromaculatus

This species is expected to be among the most abundant amphibian species in the nation. Parasites were found on one tested individual but were not assessed as threatening (nematode) [98]. The species is widespread in DPR Korea and in Northeast Asia, with apparently large and abundant populations in a large variety of habitats. The principal threats are habitat modification and environmental pollution. Recommendation for IUCN National Red List assessment (national level): near threatened. Current IUCN Global Red List assessment: near threatened [174].

3.6.16. Strauchbufo raddei

The species is vulnerable under criterion B1 on the national red list, as its EOO within the nation meets the IUCN threshold, it occurs in fewer than 10 threat-defined locations and due to estimated ongoing decline in its: (i) extent of occurrence; (ii) area of occupancy; (iii) area, extent or quality of habitat, potentially resulting in local extirpations; however, it is widespread throughout Northeast Asia [175] and there is no evidence of severe population decline or range contractions, meeting the IUCN thresholds at the global level. Recommendation for IUCN National Red List assessment (national level): vulnerable B1ab(i,ii,ii). Current IUCN Global Red List assessment: least concern [176].

3.6.17. Rana amurensis

The species is limited to northern regions but is widespread in continental northern Asia, with apparently large and abundant populations in a large variety of habitats. The principal threats are habitat modification, chemical pollution, trade and exploitation. Recommendation for IUCN National Red List assessment (national level): least concern. Current IUCN Global Red List assessment: least concern [177].

3.6.18. Rana coreana

The species is restricted to southern regions and susceptible to habitat degradation, which is projected to continue and potentially bring the EOO of the species in the nation within the vulnerable category in the near future. In addition, several subpopulations have reportedly declined and become extirpated as a result of droughts, likely linked to climate change. The principal threats are habitat modification, chemical pollution, trade and exploitation. Recommendation for IUCN National Red List assessment (national level): near threatened. Current IUCN Global Red List assessment: least concern [178].

3.6.19. Rana dybowskii

This species is expected to be among the most abundant amphibian species in the nation. The species is limited to northern regions but is widespread in continental northern Asia, with apparently large and abundant populations in a large variety of habitats. The principal threats are habitat modification, chemical pollution, trade and exploitation. Recommendation for IUCN National Red List assessment (national level): least concern. Current IUCN Global Red List assessment: least concern [179].

3.6.20. Rana huanrenensis

This species is expected to be among the most abundant amphibian species in the nation. The species is widespread in the nation and in continental Northeast Asia, with apparently large and abundant populations in a large variety of habitats. The principal threats are habitat modification, chemical pollution, trade and exploitation. Recommendation for IUCN National Red List assessment (national level): least concern. Current IUCN Global Red List assessment: least concern [180].

3.6.21. Rana uenoi

The species is vulnerable under criterion B1 on the national red list, as its EOO within the nation meets the IUCN threshold, it occurs in fewer than 10 threat-defined locations and due to estimated ongoing decline in its: (i) extent of occurrence; (ii) area of occupancy; (iii) area, extent or quality of habitat, potentially resulting in local extirpations. The principal threats are habitat modification, chemical pollution, trade and exploitation. Recommendation for IUCN National Red List assessment (national level): vulnerable B1ab(i,ii,ii). Current IUCN Global Red List assessment: not evaluated.

3.6.22. Salamandrella tridactyla

The species is restricted to high elevation habitats and is susceptible to habitat degradation, which is projected to continue and potentially bring the EOO of the species in the nation within the vulnerable category in the near future. The principal threats are habitat modification and climate change. Recommendation for IUCN National Red List assessment (national level): near threatened. Current IUCN Global Red List assessment: not evaluated.

4. Discussion

Our work represents major progress in the assimilation of knowledge about amphibians in DPR Korea. Through museum collections, molecular tools and field surveys, we were able to determine the presence of species at localities where species identification had previously been unclear. The combination of field surveys, literature review and molecular analyses also provided the first dataset, which was robust enough to model the distribution of the majority of amphibian species in DPR Korea. In addition, we determined the continuity in call properties between local populations and conspecific populations in other regions. Finally, the combination of all these variables, other risks and the analysis of the literature allowed us to determine the threat levels for most species.
The species richness estimate in this paper is the first that combines all described species, giving higher values than those that had previously been published. This is due to recent descriptions (e.g., [24]), additional data from the field [10] and modelling of species that are expected to be present but have not been found yet, such as Rana huanrenensis. In addition, some species are likely to be found in the future, such as Karsenia koreana [26], although this will require additional surveys and field exploration. Anecdotally, the presence of Karsenia koreana was mentioned during the workshop in Pyongyang in June 2019, as a researcher had been shown a picture of a very similar looking salamander by a citizen from Geumgang mountain in the southeast of the DPR Korea, within the area where the species could potentially be present. It seems important to add that R. chensinensis was not found to be potentially present in DPR Korea and is unlikely to be present [181]. This point of contention likely arises from shifts in taxonomy in the region. In 1972, two species were described in the Russian Far East, R. temporaria (currently R. dybowskii) and R. cruenta (currently R. amurensis; [182]), with R. temporaria subsequently listed as R. chensinensis [183,184]. While R. chensinensis is more closely related to R. huanrenensis than R. dybowskii [17], records were previously attributed to the wrong species and likely still are [117].
Despite some of the target species not being found during field surveys, it is not possible to confirm their absence, as demonstrated by the presence of Bufo gargarizans and Rana dybowskii in Rason in May 2018, while they were not detected in March of the same year. When species are not present in the habitat where they are expected to occur, the adequacy of the environment for their presence can be tested. This includes landscape adequacy [185], absence of urban areas [186], absence of predators [187] or adequate water quality [188]. among others. Despite the difficulty of conducting additional surveys, obtaining water quality samples from wetlands along the west coast where Dryophytes suweonensis and Pelophylax chosenicus are expected to occur would allow testing for nitrate and phosphate concentrations, as these chemicals are known to have negative impacts on the presence of these species elsewhere in their range [30].
A single Rana sp. egg clutch was detected in Rason in March 2018, an unexpectedly low number, as the genus was breeding during the same week in R Korea and in Russia, and the environmental variables measured were within the range of preference for the genus. Furthermore, a local farmer confirmed that “small brown frogs” are generally common in the area at the end of winter; thus, not detecting early breeding species during our surveys in Rason area in opposition to the presence of late breeding species [52] may have been related to differences in ecological requirements and behaviour. The absence of large forested landscapes in the area may be a negative point for the presence of amphibians requiring shelters to hibernate, such as R. dybowskii, as opposed to species that hibernate under water, such as Pelophylax sp. or R. coreana [189]. This, however, does not match with the large presence of D. japonicus in the area, a species that migrates to forested hills for brumation and hibernations [190].
Despite the age of samples used for molecular analyses, we were able to assess the presence of R. uenoi as far north as Pyongyang (Figure 1). This places the contact zone between R. uenoi and R. dybowskii somewhere between Pyongyang and the border with PR China. This pattern is relatively common for the species studied here; for instance, R. coreana is present in R Korea, but R. amurensis and R. kukunoris are present further north in PR China and Russia. The pattern is the same for R. huanrenensis, Onychodactylus koreanus and O. zhaoermii. This contact zone between sister clades likely results from the shift in ecological properties in the region ([55]; Figure 1): the density and presence of deciduous forests decreases along with wetness and vegetation on a northward gradient, matching with a decrease in species richness (Figure 8) and increases in Manchurian and Changbai mixed forests. We expect the species to be locally adapted to the habitat type, such as commonly seen in amphibians [191,192,193,194]. This pattern can also be generalised for other fauna, as it is also observed in mammals [195] and birds [196].
This contact zone may also have been shaped by climatic variations during ice ages, as geological and paleoclimatic variations had significant impacts on speciation [197,198,199,200]. When these variations result in the loss of connectivity between populations, such as during isolation in refugia, populations can split and diverge into segregated clades [201,202]. In the Northern Hemisphere, glacial maxima pushed most species further south, including populations in Northeast Asia [203,204], and consequently increased geographic distances between metapopulations [205,206]. These older boundaries may form the current contact zone between species, such as seen in Hyla orientalis in the Balkans [207], or result in population expansion from refugia, as hypothesised for raccoon dogs in north-east Asia [204]. Here, we hypothesise that this pattern of contact zones arose from climatic variations for the Rana complex, which diverged as a result of the orogenesis of the Changbai Mountain Range during the Miocene [42]. In addition, this orogenesis event may also have had an impact on some other amphibian species, as seen through temporally matching divergences in other amphibian genera [37,208,209].
Interestingly, the lower number of species in the northern regions of DPR Korea contrasts with species distributed around the Yellow Sea, where weaker patterns may result from the absence of elevational barriers [208,209,210,211,212,213,214,215]. Species bordering the Yellow Sea are generally distributed in lowlands and show lower interspecific variation than those present to the north and south of DPR Korea. This is demonstrated by D. suweonensis and D. immaculatus [79], which are distributed on either sides of the Yellow Sea and recently diverged [37]. Other species, such as Kaloula borealis [216] and Pelophylax sp. [208,217,218], exhibit only low genetic divergence along the same geographic gradient, a potential result of the Yellow Sea being drained of salt water repeatedly during geological ages, enabling gene exchange between the Asian mainland, the Korean Peninsula and the Japanese archipelago [219,220,221,222].
Bufo gargarizans is especially interesting, in that it is distributed both south and north of the Changbai Mountain Range, east and west of the Yellow Sea [162] and as far south as Vietnam [223]. It is a monophyletic species throughout its range [209,224,225,226,227,228]; however, two clades are present within B. gargarizans [209], and further genetic analyses may provide a different picture for this species. One of the reasons why B. gargarizans is not segregated into several clades like Rana sp. is that the recent colonisation of more northern latitudes may have followed on from a drop in sea level, connecting refugia on the Korean Peninsula and the Chinese mainland south of the permafrost line [226,229,230]. These movements would have been facilitated by the drainage basins of the Han, Amur, Liao and Yellow rivers into a single water body flowing into the ocean south of the current Yellow Sea [231,232] and by the development of the monsoon system [233,234].
In terms of species richness, the pattern found through models was in line with expected results, with fewer species at higher latitude and altitude [235,236,237,238,239,240], which is a consistent pattern in amphibians [236,241,242,243,244]. The patterns found here could be slightly different if additional data were available for Caudata, as their elevation distribution patterns are different from that of Anura on the Korean Peninsula, which are found in higher numbers at intermediate elevations, and also the species for which we have the least data available in the analyses presented here. We recommend numerous additional surveys with a specific focus on species boundaries for Onychodactylus sp. and the presence of K. koreana. Models will, however, have to take into account omission and commission errors generally associated with field surveys [134,245].
Based on available foreign reports on the ecology of DPR Korea, low population size and diversity would have been expected [153,156]; however, when considering the populations sizes of all amphibian species encountered during surveys in Mundeok, Kumya and Rason, the populations appeared to be comparatively healthy, in agreement with older [51] and relevant literature [5,9]. For instance, audio recordings of D. suweonensis in Mundeok (Figure 1) and D. japonicus in Kumya would indicate population sizes comparable to the largest ones in R Korea [36,80,81,165], while population sizes of P. nigromaculatus, D. japonicus and B. orientalis in Rason would likewise be comparable. This may be due to several factors, primarily the presence of large and uninterrupted agricultural wetlands, increasing the surface of available habitat matching with the ecological requirements of lowland amphibian species [246,247,248,249,250] and improving population connectivity [251,252,253]. Other reasons for large populations may be the adequacy of agricultural practices with the life cycle and behaviour of the species, where flooding matches with the breeding activities of the species and harvesting occurs after the dispersion of juveniles [254,255,256,257,258,259,260,261,262]. The comparison in population sizes is also supported by the similar morphology, as it relates to similar growth pressures, following the methodology of [263]. Furthermore, several natural wetlands connected by agricultural wetlands are now protected under the Ramsar convention, which should benefit species in the area.
In addition, the conservation status of species present at higher altitudes, such as B. stejnegeri and Onychodactylus sp., could have been improved by the designation of protected areas in DPR Korea, including biosphere reserves recognised by the UNESCO Man and the Biosphere Programme [153]. These are mostly located around Mount Baekdu and the Changbai Mountain, connected with the Biosphere Reserve in China, protecting the highest mountain on mainland Northeast Asia and also extending south along the Baekdu Mountain Range all the way to R Korea, where much of the mountain range is also included in protected area systems [153].
In term of specific threats to species and from the data gathered here, it would seem that P. nigromaculatus, D. japonicus, K. borealis, B. orientalis, G. emeljanovi, B. gargarizans and most Rana sp. are not overly threatened in DPR Korea, although small EOOs and the presence of ongoing threats do result in a threatened status for some species. The situation for D. suweonensis, S. tridactyla, B. stejnegeri and S. raddei would seem the most critical because of the small EOOs, likely even smaller AOOs and a combination of other factors, including but are not limited to habitat degradation, especially deforestation, climate change and pollution. In the case of Onychodactylus sp., it is currently not possible to clarify the conservation status, as there are too little data available regarding the identity of records included in the current concept.
The ecological data available from the DPR Korea are slowly increasing through bird and wetland resurveys, many of which have been conducted in recent years as part of national accession to the Ramsar Convention and formal participation in the East Asian–Australasian Flyway Partnership. Although understanding of the distribution and abundance of better-studied organisms such as birds is improving in DPR Korea, the scope and depth of such research still remains very limited compared with many other regions, and data on other species groups such as amphibians and fish are even more limited. Further field surveys combined with molecular and modelling tools have the potential to refine current the understanding of species boundaries and conservation status and to detect additional species. Such research would also be very helpful in supporting national and regional efforts to conserve biodiversity on the Korean Peninsula and throughout Northeast Asia.

5. Conclusions

Our results showed that species richness is higher at lower altitudes and latitudes, with up to 11 species being present. In opposition, the species richness in northern regions was half that of the lower areas. Among the amphibians species found in DPR Korea, we recommend ten species as threatened at the national level following the IUCN Red List categories and criteria. This high number is in line with threat levels to amphibians in bordering nations and globally. Species relying on agricultural wetlands such as rice paddies seem more abundant and are not under imminent threat due to the enduring presence of extensive agricultural landscapes providing adequate habitats for the species. The maintenance of such landscapes is a clear benefit to the conservation of amphibian species.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ani11072057/s1: Table S1: Dataset containing GPS coordinates. Table S2: GPS coordinates for all geolocalised amphibian presence in DPR Korea. Table S3: Percent contribution of bioclimate variables. Table S4: Permutation importance of bioclimate variables. Table S5: Training gain with individual bioclimate variables. Table S6: List of samples and matching sequences used in this study.

Author Contributions

Conceptualisation, A.B. and S.N.L.; methodology, A.B., S.N.L., D.A., S.K., S.N.O., J.L. and L.H.; software, A.B. and D.A.; validation, A.B., S.N.L., K.R., D.A., T.Y.N., G.H.J., H.S.M., J.S.C., S.K., S.N.O., K.M., Y.B., Y.S., I.M., J.L., L.H., N.M., B.S. and F.G.; formal analysis, A.B., D.A., S.N.O. and A.K.; investigation, A.B., S.K., S.N.O., K.M., N.M., B.S. and F.G.; resources, A.B., S.N.L., K.R., T.Y.N., B.S. and F.G.; data curation, A.B., S.N.L., D.A. and A.K.; writing—original draft preparation, A.B.; writing—review and editing, A.B., S.N.L., K.R., D.A., T.Y.N., G.H.J., H.S.M., J.S.C., S.K., S.N.O., K.M., Y.B., Y.S., A.K., I.M., J.L., L.H., N.M., B.S., F.G. and Y.J.; visualisation, A.B., S.N.L., D.A., T.Y.N., G.H.J., H.S.M., J.S.C., S.K., S.N.O., K.M., Y.B., Y.S., A.K., I.M., J.L., L.H., N.M., B.S., F.G. and Y.J.; supervision, A.B. and N.M.; project administration, A.B., B.S., F.G. and K.R.; funding acquisition, A.B. All authors have read and agreed to the published version of the manuscript. A.B., S.N.L., D.A., S.K., S.N.O., K.M., Y.B., Y.S., A.K., I.M., J.L., L.H., N.M., B.S., F.G. and Y.J. agreed in writing to all versions on this manuscript. K.R., T.Y.N., G.H.J., H.S.M. and J.S.C. had access to the latest version of the manuscript but did not manage to provide a formal agreement in writing due to the restrictions related to the COVID-19 pandemic. The version agreed on in writing can be found in the revised version of this manuscript.

Funding

This project was financially supported by the Hanns Seidel Foundation. The APC was funded by Nanjing Forestry University.

Institutional Review Board Statement

Not applicable. Surveys were conducted under the supervision of the Ministry of Land and Environment Protection of the Democratic People’s Republic of Korea.

Data Availability Statement

The data for this study are provided in the supplementary files and are available from the Mendeley Data Repository http://dx.doi.org/10.17632/z7kgyy8chp.1 (published on 28 April 2021).

Acknowledgments

Despite not being able to trace back to a specific name, we would like to acknowledge everyone from the Chosun Academy of Science involved in data collection, without whom the molecular part of this work would not have been possible. We are grateful to I.V. Doronin, K.D. Milto, N.B. Ananjeva, V.I. Kazakov, E.A. Kazakova, V.F. Orlova and employees of the S.M. Kirov Military Medical Academy who allowed us to study the museum collections. Finally, we are grateful to Catharina Karlsson for helping with the acoustic analyses. With the exception of museum vouchers, no specimens were handled and field work was conducted under the supervision of the Ministry of Land and Environment Protection of the Democratic People’s Republic of Korea.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. IUCN. The IUCN Red List of Threatened Species; IUCN, Ed.; Version 2020-1; IUCN: Gland, Switzerland, 2020; Available online: https://www.iucnredlist.org (accessed on 19 March 2020).
  2. Tomek, T. The birds of North Korea. Non-Passeriformes. Acta Zool. Crac. 1999, 42, 1–21. [Google Scholar]
  3. Tomek, T. The birds of North Korea. Passeriformes. Acta Zool. Crac. 2002, 45, 1–235. [Google Scholar]
  4. Moores, N. Birds and their conservation in Rason Special Economic Zone, Democratic People’s Republic of Korea. Forktail 2018, 33, 124–133. [Google Scholar]
  5. MAB National Committee of DPR Korea. Biodiversity of Mt. Kumgang; UNESCO-MAB Programme: Pyongang, Democratic People’s Republic of Korea, 2003; Volume 92. [Google Scholar]
  6. Li, M. Research on the organic components of toads. Biol. Sci. 2014, 103, 57. [Google Scholar]
  7. MAB National Committee of DPR Korea. Red Data Book of DPRK; National Committee of DPR Korea, Ed.; Science and Technology Publisher: Pyongyang, Democratic People’s Republic of Korea, 2016; Volume Animal. [Google Scholar]
  8. Chong, J.-R.; Hong, Y.J.; Hong, Y.G.; Ie, D.S. A Directory of Wetlands of D. P. R. Korea; Wild Bird Society of Japan and Korea University: Pyongyang, Democratic People’s Republic of Korea, 2016. [Google Scholar]
  9. Kim, L.; Han, G. Chosun Animal Encyclopedia, Herpetology Volume; Roh, H., Ed.; Science and Technology Publisher: Pyongyang, Democratic People’s Republic of Korea, 2009; p. 138. [Google Scholar]
  10. Ri, J. Systematics of the frogs (Anura) of our country. Biol. Sci. 2018, 107, 64. [Google Scholar]
  11. Song, J.-Y. Current status and distribution of Amphibian in North Korea. J. Natl. Park Res. 2016, 7, 211–217. [Google Scholar]
  12. Werner, F. Über Reptilien und Batrachier aus Guatemala und China in der zoologischen Staats-Sammlung in München nebst einem Abhang über seltene Formen aus Anderer gebieten. Abh. Der Math. Phys. Cl. Königlich Bayer. Akad. Wiss. 1903, 22, 342–384. [Google Scholar]
  13. Kaneko, Y.; Matsui, M. The IUCN Red List of Threatened Species; IUCN, Ed.; IUCN: Gland, Switzerland, 2004; Volume e.T58688A11811732. [Google Scholar]
  14. Kuzmin, S.; Ishchenko, V.; Tuniyev, B.; Beebee, T.; Andreone, F.; Nyström, P.; Anthony, B.P.; Schmidt, B.; Ogrodowczyk, A.; Ogielska, M.; et al. The IUCN Red List of Threatened Species; IUCN, Ed.; IUCN: Gland, Switzerland, 2009; Volume e.T58734A86470817. [Google Scholar]
  15. Xie, F.; Ye, C.; Fei, L.; Jiang, J.; Zeng, X.; Matsui, M. Taxonomical studies on Brown frogs (Rana) from Noertheastern China (amphibian: Ranidae). Acta Zootaxonomica Sin. 1999, 24, 224–231. [Google Scholar]
  16. Lu, Y.-Y.; Li, P.-P. A brief review on advance of Wood frogs research in China and discussing the importance of studying the Wood frog biodiversity around Bohai. Sichuan J. Zool. 2004, 24, 271–275. [Google Scholar]
  17. Yang, J.-X.; Zhou, W.-W.; Rao, D.-Q.; Poyarkov, A.N.; Kuzmin, S.L.; Che, J. Validity and systematic position of Rana altaica (Rana: Ranidae): Results of a phylogenetic analysis. Zool. Res. 2010, 31, 353–360. [Google Scholar]
  18. Schmidt, K.P. A new toad from Korea. Copeia 1931, 1931, 93–94. [Google Scholar] [CrossRef]
  19. Lim, C.H. Abour a new speice Bufo sambangi sp. n. Sci. Bull. 1996, 258, 39–41. [Google Scholar]
  20. Duellman, W.E.; Marion, A.B.; Hedges, S.B. Phylogenetics, classification, and biogeography of the treefrogs (Amphibia: Anura: Arboranae). Zootaxa 2016, 4104, 1–109. [Google Scholar] [CrossRef] [Green Version]
  21. Kuramoto, M. Mating calls of treefrogs (genus Hyla) in the far east, with description of a new species from Korea. Copeia 1980, 1, 100–108. [Google Scholar] [CrossRef]
  22. Chun, S.; Chung, E.; Voloshina, I.; Chong, J.R.; Lee, H.; Min, M.-S. Genetic diversity of korean tree frogs (Hyla suweonensis and Hyla japonica): Assessed by mitochondrial cytochrome b gene and cytochrome oxidase subunit I gene. Korean J. Herpetol. 2012, 4, 31–41. [Google Scholar]
  23. Borzée, A.; Seliger, B. Dryophytes suweonensis (Suweon Treefrog). Herpetol. Rev. 2018, 49, 707. [Google Scholar]
  24. Poyarkov, J.N.; Che, J.; Min, M.-S.; Kuro-o, M.; Yan, F.; Li, C.; Lizuka, K.; Vieites, D. Review of the systematics, morphology and distribution of Asian Clawed Salamanders, genus Onychodactylus (Amphibia, Caudata: Hynobiidae), with the description of four new species. Zootoxa 2012, 3465, 1–106. [Google Scholar] [CrossRef]
  25. Maslova, I.V.; Portnyagina, E.Y.; Sokolova, D.A.; Vorobieva, P.A.; Akulenko, M.V.; Portnyagin, A.S.; Somov, A.A. Distribution of rare and endangered amphibians and reptiles in Primorsky Krai (far east, Russia). Nat. Conserv. Res. 2018, 3, 61–72. [Google Scholar] [CrossRef]
  26. Borzée, A.; Andersen, D.; Groffen, J.; Kim, H.-T.; Bae, Y.; Jang, Y. Climate change-based models predict range shifts in the distribution of the only Asian plethodontid salamander: Karsenia koreana. Sci. Rep. 2019, 9, 11838. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Bishop, P.J.; Angulo, A.; Lewis, J.P.; Moore, R.D.; Rabb, G.B.; Garcia Moreno, J. The amphibian extinction crisis-what will it take to put the action into the Amphibian Conservation Action Plan? Surv. Perspect. Integr. Environ. Soc. 2012, 5, 97–111. [Google Scholar]
  28. Lannoo, M.J. A perspective on amphibian conservation in the United States. Alytes 2012, 29, 133–144. [Google Scholar]
  29. Okada, Y. The tailless bactracians of the Japanese Empire. J. Imp. Agric. Exp. Stn. Nishigahara Tokyo 1931, 29, 1–215. [Google Scholar]
  30. Borzée, A.; Kyong, C.N.; Kil, H.K.; Jang, Y. Impact of water quality on the occurrence of two endangered Korean anurans: Dryophytes suweonensis and Pelophylax chosenicus. Herpetologica 2018, 74, 1–7. [Google Scholar] [CrossRef]
  31. Yeum, J.-H.; Lee, J.-H.; Lee, H.-J. Designation Properties of Ramsar Site for 20 Years in Korea. J. Environ. Sci. Int. 2018, 27, 319–332. [Google Scholar] [CrossRef]
  32. Son, E.; Borzée, A. Pelophylax chosenicus (Gold-spotted Pond Frog), Geographic distribution note. Herpetol. Rev. 2020, 51, 134. [Google Scholar]
  33. Borzée, A.; Bae, Y.H.; Seliger, B. Pelophylax chosenicus (Gold-spotted Pond Frog). Herpetol. Rev. 2018, 49, 709. [Google Scholar]
  34. Gunther, A. Catalogue of the Batrachia Salientia in the Collection of the British Museum; Collection of the British Museum; British Museum: London, UK, 1958. [Google Scholar]
  35. Kuzmin, S.; Maslova, I.; Matsui, M.; Liang, F.; Kaneko, Y. Dryophytes Japonicus (Amended Version of 2014 Assessment); A. S. Group, Ed.; The IUCN Red List of Threatened Species; IUCN: Gland, Switzerland, 2017. [Google Scholar]
  36. Roh, G.; Borzée, A.; Jang, Y. Spatiotemporal distributions and habitat characteristics of the endangered treefrog, Hyla suweonensis, in relation to sympatric H. japonica. Ecol. Inform. 2014, 24, 78–84. [Google Scholar] [CrossRef]
  37. Dufresnes, C.; Litvinchuk, S.N.; Borzée, A.; Jang, Y.; Li, J.-T.; Miura, I.; Perrin, N.; Stöck, M. Phylogeography reveals an ancient cryptic radiation in East-Asian tree frogs (Hyla japonica group) and complex relationships between continental and island lineages. BMC Evol. Biol. 2016, 16, 253. [Google Scholar]
  38. Fei, L.; Ye, C.-Y.; Huang, Y.-Z.; Liu, M. Atlas of amphibians of China; Henan Science and Technology Press: Zhengzhou, China, 1999. [Google Scholar]
  39. Liu, P.; Zhang, Z.; Zhao, S.; Zhao, W. Effect of female reproduction and mate choice on sexual size dimorphism in the northeast treefrog Hyla ussuriensis (Anura: Hylidae) in China. Asian Herpetol. Res. 2012, 3, 273–279. [Google Scholar]
  40. Matsui, M. Description of a New Brown Frog from Tsushima Island, Japan (Anura: Ranidae: Rana). Zool. Sci. 2014, 31, 613–620. [Google Scholar] [CrossRef] [Green Version]
  41. Günther, A.C.L.G. Description of a new frog from north-eastern Asia. Ann. Mag. Nat. Hist. 1876, 4, 387. [Google Scholar] [CrossRef]
  42. Yang, B.-T.; Zhou, Y.; Min, M.-S.; Matsui, M.; Dong, B.-J.; Li, P.-P.; Fong, J.J. Diversity and phylogeography of Northeast Asian brown frogs allied to Rana dybowskii (Anura, Ranidae). Mol. Phylogenet. Evol. 2017, 112, 148–157. [Google Scholar] [CrossRef]
  43. Boulenger, G.A. Notes sur les grenouilles rousses d’Asie. Bull. Société Zool. Fr. 1886, 11, 595–600. [Google Scholar]
  44. Song, J.-Y.; Matsui, M.; Chung, K.-H.; Oh, H.-S.; Zhao, W. Distinct specific status of the Korean brown frog, Rana amurensis coreana (Amphibia: Ranidae). Zool. Sci. 2006, 23, 219–224. [Google Scholar] [CrossRef] [Green Version]
  45. Kim, J.-B.; Min, M.-S.; Yang, S.-Y.; Matsui, M. Genetic relationships among Korean brown frog species (Anura, Ranidae), with special reference to evolutionary divergences between two allied species Rana dybowskii and R. huanrenensis. Zool. Sci. 2002, 19, 369–382. [Google Scholar] [CrossRef] [Green Version]
  46. Liu, M.-Y.; Zhang, S.-Q.; Liu, M. A new species of Ranidae fron Liaoning, China (Anura). Acta Zootaxonomica Sin. 1993, 18, 493–497. [Google Scholar]
  47. Yang, S.Y.; Kim, J.B.; Min, M.S.; Suh, J.H.; Kang, Y.J.; Matsui, M.; Fei, L. First record of a brown frog Rana huanrenensis (Family Ranidae) from Korea. Korean J. Biol. Sci. 2000, 4, 45–50. [Google Scholar] [CrossRef]
  48. Che, J.; Pang, J.; Zhao, E.-M.; Matsui, M.; Zhang, Y.-P. Phylogenetic relationships of the Chinese brown frogs (genus Rana) inferred from partial mitochondrial 12S and 16S rRNA gene sequences. Zool. Sci. 2007, 24, 71–80. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Zhou, Y.; Yang, B.-T.; Li, P.-P.; Min, M.-S.; Fong, J.J.; Dong, B.-J.; Zhou, Z.-Y.; Lu, Y.-Y. Molecular and morphological evidence for Rana kunyuensis as a junior synonym of Rana coreana (Anura: Ranidae). J. Herpetol. 2015, 49, 302–307. [Google Scholar] [CrossRef]
  50. Shannon, F.A. The Reptiles and Amphibians of Korea. Herpetologica 1956, 12, 22–49. [Google Scholar]
  51. Szyndlar, Z. A description of a small collection of amphibians and reptiles from the Democratic People’s Republic of Korea with notes on the distribution of the herpetofauna in that country. Acta Zool. Crac. 1984, 17, 1–11. [Google Scholar]
  52. Borzée, A.; Groffen, J.; Seliger, B. Bombina orientalis (Oriental Fire-bellied Toad). Herpetol. Rev. 2018, 49, 69. [Google Scholar]
  53. Stuart, S.N.; Chanson, J.S.; Cox, N.A.; Young, B.E.; Rodrigues, A.S.; Fischman, D.L.; Waller, R.W. Status and trends of amphibian declines and extinctions worldwide. Science 2004, 306, 1783–1786. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Wake, D.B. Facing extinction in real time. Science 2012, 335, 1052–1053. [Google Scholar] [CrossRef] [PubMed]
  55. Olson, D.M.; Dinerstein, E.; Wikramanayake, E.D.; Burgess, N.D.; Powell, G.V.N.; Underwood, E.C.; D’Amico, J.A.; Itoua, I.; Strand, H.E.; Morrison, J.C.; et al. Terrestrial Ecoregions of the World: A New Map of Life on Earth. BioScience 2013, 51, 933–938. [Google Scholar] [CrossRef]
  56. de Magalhães, R.F.; Rocha, P.C.; Santos, F.R.; Strüssmann, C.; Giaretta, A.A. Integrative taxonomy helps to assess the extinction risk of anuran species. J. Nat. Conserv. 2018, 45, 1–10. [Google Scholar] [CrossRef]
  57. Rohland, N.; Hofreiter, M. Ancient DNA extraction from bones and teeth. Nat. Protoc. 2007, 2, 1756. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Nguyen, H.Q.; Kim, Y.I.; Borzée, A.; Jang, Y. Efficient isolation method for high-quality genomic DNA from cicada exuviae. Ecol. Evol. 2017, 7, 8161–8169. [Google Scholar] [CrossRef] [Green Version]
  59. Campos, P.F.; Gilbert, T.M. Ancient DNA; Shapiro, B., Hofreiter, M., Eds.; Humana Press: New York, NY, USA, 2012; Volume 840, pp. 81–85. [Google Scholar]
  60. Gilbert, M.T.P.; Haselkorn, T.; Bunce, M.; Sanchez, J.J.; Lucas, S.B.; Jewell, L.D.; Van Marck, E.; Worobey, M. The isolation of nucleic acids from fixed, paraffin-embedded tissues—Which methods are useful when? PLoS ONE 2007, 2, e537. [Google Scholar] [CrossRef]
  61. Qiagen. Sample & Assay Technologies. In QIAamp DNA FFPE Tissue Handbook; Qiagen: Venlo, The Netherlands, 2012. [Google Scholar]
  62. Hua, X.; Fu, C.; Li, J.; Oca, A.N.M.D.; Wiens, J.J. A revised phylogeny of holarctic treefrogs (genus Hyla) based on nuclear and mitochondrial DNA sequences. Herpetologica 2009, 65, 246–259. [Google Scholar] [CrossRef]
  63. Jeong, T.J.; Jun, J.; Han, S.; Kim, H.T.; Oh, K.; Kwak, M. DNA barcode reference data for the Korean herpetofauna and their applications. Mol. Ecol. Resour. 2013, 13, 1019–1032. [Google Scholar] [CrossRef]
  64. Igawa, T.; Kurabayashi, A.; Usuki, C.; Fujii, T.; Sumida, M. Complete mitochondrial genomes of three neobatrachian anurans: A case study of divergence time estimation using different data and calibration settings. Gene 2008, 407, 116–129. [Google Scholar] [CrossRef]
  65. Li, J.-T.; Wang, J.-S.; Nian, H.-H.; Litvinchuk, S.N.; Wang, J.; Li, Y.; Rao, D.-Q.; Klaus, S. Amphibians crossing the bering land bridge: Evidence from holarctic treefrogs (Hyla, Hylidae, Anura). Mol. Phylogenet. Evol. 2015, 87, 80–90. [Google Scholar] [CrossRef]
  66. Borzée, A.; Didinger, C.; Jang, Y. The Complete Mitochondrial Genome of Dryophytes suweonensis (Anura Hylidae). Mitochondrial DNA Part B Resour. 2017, 1, 5–6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Lee, M.-Y.; Jeon, H.S.; Min, M.-S.; An, J. Sequencing and analysis of the complete mitochondrial genome of Hyla suweonensis (Anura: Hylidae). Mitochondrial DNA Part B 2017, 2, 126–127. [Google Scholar] [CrossRef] [PubMed]
  68. Sun, Q.; Xie, Y.; Zhao, W.; Liu, P. Sequencing and analysis of the complete mitochondrial genome of Hyla ussuriensis (Anura: Hylidae). Mitochondrial DNA Part A 2017, 28, 328–329. [Google Scholar] [CrossRef] [PubMed]
  69. Sumida, M.; Ueda, H.; Nishioka, M. Reproductive isolating mechanisms and molecular phylogenetic relationships among Palearctic and Oriental brown frogs. Zool. Sci. 2003, 20, 567–580. [Google Scholar] [CrossRef] [PubMed]
  70. Zhou, Y.; Wang, S.; Zhu, H.; Li, P.; Yang, B.; Ma, J. Phylogeny and biogeography of South Chinese brown frogs (Ranidae, Anura). PLoS ONE 2017, 12, e0175113. [Google Scholar]
  71. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [Green Version]
  72. Guindon, S.; Gascuel, O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst. Biol. 2003, 52, 696–704. [Google Scholar] [CrossRef] [Green Version]
  73. Okada, Y. Frogs in Korea. J. Chosen Biol. Soc. 1928, 6, 15–46. [Google Scholar]
  74. Milto, K.D.; Barabanov, A.V. An annotated catalogue of the amphibian types in the collection of the Zoological Institute, Russian Academy of Sciences, St. Petersburg. Russ. J. Herpetol. 2011, 18, 137–153. [Google Scholar]
  75. Berman, D.B.N.; Balan, I. The Siberian frog. Priroda 2017, 8. n.a.. [Google Scholar]
  76. Stejneger, L. Herpetology of Japan and Adjacent Territory; Smithsonian Institute: Washington, DC, USA, 1907. [Google Scholar]
  77. Borzée, A.; Park, S.; Kim, A.; Kim, H.-T.; Jang, Y. Morphometrics of two sympatric species of tree frogs in Korea: A morphological key for the critically endangered Hyla suweonensis in relation to H. japonica. Anim. Cells Syst. 2013, 17, 348–356. [Google Scholar] [CrossRef]
  78. Borzée, A.; Messenger, K.R.; Chae, S.; Andersen, D.; Groffen, J.; Kim, Y.I.; An, J.; Othman, S.; Ri, K.; Nam, T.Y.; et al. Yellow sea mediated segregation between North East Asian Dryophytes species. PLoS ONE 2020, 15, e0234299. [Google Scholar]
  79. Borzée, A.; Kong, S.; Didinger, C.; Nguyen, H.; Jang, Y. A ring-species or a ring of species? Phylogenetic relationship between two treefrog species, Dryophytes suweonensis and D. immaculatus, around the Yellow Sea. Herpetol. J. 2018, 28, 160–170. [Google Scholar]
  80. Borzée, A.; Kim, K.; Heo, K.; Jablonski, P.G.; Jang, Y. Impact of land reclamation and agricultural water regime on the distribution and conservation status of the endangered Dryophytes suweonensis. PeerJ 2017, 5, e3872. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Borzée, A.; Andersen, D.; Jang, Y. Population trend inferred from aural surveys for calling anurans in Korea. PeerJ 2018, 6, e5568. [Google Scholar] [CrossRef] [Green Version]
  82. Smith, G.W.; Nydegger, N.C. A spotlight, line-transect method for surveying jack rabbits. J. Wildl. Manag. 1985, 49, 699–702. [Google Scholar] [CrossRef]
  83. R Core Team. R Foundation for Statistical Computing; R Core Team: Vienna, Austria, 2017. [Google Scholar]
  84. Araya-Salas, M.; Smith-Vidaurre, G. WarbleR: An R package to streamline analysis of animal acoustic signals. Methods Ecol. Evol. 2017, 8, 184–191. [Google Scholar] [CrossRef]
  85. Sueur, J.; Aubin, T.; Simonis, C.; Lellouch, L.; Brown, E.C.; Depraetere, M.; Desjonqueres, C.; Fabianek, F.; Gasc, A.; Kasten, E.; et al. Package “Seewave”. 2018. Available online: https://cran.r-project.org/web/packages/seewave/index.html (accessed on 1 December 2020).
  86. Ligges, U.; Krey, S.; Mersmann, O.; Schnackenberg, S.; Guenard, G.; Preusser, A.; Thieler, A.; Mielke, J.; Weihs, C. Package “TtuneR”. 2008. Available online: https://cran.r-project.org/web/packages/tuneR/index.html (accessed on 1 December 2020).
  87. Park, S.; Jeong, G.; Jang, Y. No reproductive character displacement in male advertisement signals of Hyla japonica in relation to the sympatric H. suweonensis. Behav. Ecol. Sociobiol. 2013, 67, 1345–1355. [Google Scholar] [CrossRef]
  88. Won, H.K.; Choy, Y.G. Checklist of Korean Chordata: Amphibia, Reptilia; Publication of the Korean Academy of Science: Pyongyang, Democratic People’s Republic of Korea, 1967; pp. 167–245. [Google Scholar]
  89. Won, H.K. Amphibian and Reptilian Fauna of Korea; Publication of the Korean Academy of Science: Pyongyang, Democratic People’s Republic of Korea, 1971; p. 170. [Google Scholar]
  90. Nikolsky, A.M. Faune de la Russie et des Pays limitrophes; Izd. Ros. Akad. Nauk: Petrograd, Russia, 1918; p. 309. [Google Scholar]
  91. Slevin, J. Contributions to oriental herpetology II. Korea or Chosen. Calif. Acad. Sci. 1925, 115, 89–100. [Google Scholar]
  92. Okada, Y. A study on the distribution of tailess batrachians of Japan. Annot. Zool. Jpn. 1926, 11, 137–144. [Google Scholar]
  93. Slevin, J.R. Contributions to Oriental herpetology. V. Honshu Hondo Neighboring Isl. Sado Awaji Seven Isl. Idzu. Ibid 1937, 23, 175–190. [Google Scholar]
  94. Kang, Y.S.; Yoon, I.B. Amphibia and Reptilia; Illustrated Encyclopedia of Fauna and Flora of Korea; Ministry of Education: Seoul, Republic of Korea, 1975; Volume 17, p. 191. [Google Scholar]
  95. Lee, S.; Lee, K. The Yearbook of Natural Monument in D. P. R. Korea; Agriculture Press: Pyeongyang, Democratic People’s Republic of Korea, 1994. [Google Scholar]
  96. Pawlowski, J.; Tomek, T. Zoological expeditions to the North Korea organized in the years 1971–1992 by the Cracow Institute of Systematics and Evolutions of Animals of the Polish Academy of Science. Fragm. Faun. 1997, 40, 231–246. [Google Scholar] [CrossRef] [Green Version]
  97. Cantor, T. General features of Chusan, with remarks on the flora and fauna of that island. Ann. Mag. Nat. Hist. 1842, 1, 481–493. [Google Scholar] [CrossRef] [Green Version]
  98. Mészáros, F.; Zombori, L. Zoological collectings by the Hungarian Natural History Museum in Korea. 107. A report of the collecting of the fourteenth expedition. Folia Entomol. Hung. 1992, 52, 73–80. [Google Scholar]
  99. Lee, S.; Gil, J. A study of mammal and amphibians-reptile species around Shin-Po area, Hamkyungbuk-do, North Korea. Korean Wetl. Soc. 2007, 9, 13–20. [Google Scholar]
  100. Schmidt, P.Y. News from the Korean-Sakhalin expedition [Vesti is Koreysko-Sakhalinskoy ekspeditsii]. Izv. Imp. Rus. Geogr. Obs. St. Petersburg. 1900, 36, 463–501. [Google Scholar]
  101. Borkin, L. On rare Far Eastern toad species Bufo stejnegeri Schmidt, 1931. Proc. Zool. Inst. Leningr. 1986, 157, 192–196. [Google Scholar]
  102. Strauch, A. Mongolia and the Tangut Country: A Three Years’ Journey in the Eastern Highlands of Asia; Izd. Imp. Russkago Ob-va: Sankpeterburg, Russia, 1876; Volume 2. [Google Scholar]
  103. Boulenger, G. Annals and Magazine of Natural History (5), XIX; not of Audouin; Taill. Batr. Europe, Pt.: Gensan, Korea, 1887; p. 67. [Google Scholar]
  104. Borzée, A. First amphibian behavioural observation from the Democratic People’s Republic of Korea: Predation of a Dryophytes japonicus tadpole by Hydaticus sp. larvae. Nat. Conserv. Res. 2019, 4, 106–108. [Google Scholar] [CrossRef]
  105. Boulenger, G.A. A list of the reptiles and batrachians of Amoorland. Ann. Mag. Nat. Hist. 1890, 6, 137–144. [Google Scholar] [CrossRef] [Green Version]
  106. Zvegintsov, A. Travel in North Korea [Poezdka v severnuyu Koreyu]. Izv. Imp. Rus. Geogr. Obs. St. Petersburg 1900, 36, 502–518. [Google Scholar]
  107. Suk, J. Collection record of animal in Mt. Baekdusan area. J. Chosen Biol. Soc. 1934, 18, 129–136. [Google Scholar]
  108. California Academy of Sciences. Herpetology 32497; California Academy of Sciences, Ed.; California Academy of Sciences: San Francisco, CA, USA, 2017; Volume 2019. [Google Scholar]
  109. California Academy of Sciences. Herpetology 32483; California Academy of Sciences, Ed.; California Academy of Sciences: San Francisco, CA, USA, 2017; Volume 2019. [Google Scholar]
  110. California Academy of Sciences. Herpetology 32519; California Academy of Sciences, Ed.; California Academy of Sciences: San Francisco, CA, USA, 2017; Volume 2019. [Google Scholar]
  111. Barbour, T. Some new reptiles and amphibians. Bull. Mus. Comp. Zool. 1908, 51, 315–325. [Google Scholar]
  112. Kim, R.; Ju, J.; Sin, G.; Choi, G.; Park, H.; Go, Y.; Seok, S. The Illustrated Book of the Korean Natural Monuments 3: Animals Part; Seong, J., Kim, S., Eds.; Korean Culture Preservation Company: Pyongyang, Democratic People’s Republic of Korea, 2007; Volume 96. [Google Scholar]
  113. Choi, D.; Park, C.; Jang, U.; Jeong, S.; Go, I.; Park, Y.; Gang, C. Gwangmyeong Encyclopedia: Agriculture, Forestry, Fishery; Gwangmyeong Encyclopedia; Pyongyang Printing Press: Pyongyang, Democratic People’s Republic of Korea, 2009; Volume 18, p. 888. [Google Scholar]
  114. California Academy of Sciences. Herpetology 36028; California Academy of Sciences, Ed.; California Academy of Sciences: San Francisco, CA, USA, 2017; Volume 2019. [Google Scholar]
  115. California Academy of Sciences. Herpetology 32735; California Academy of Sciences, Ed.; California Academy of Sciences: San Francisco, CA, USA, 2017; Volume 2019. [Google Scholar]
  116. Bannikov, A.; Darevsky, I.; Ishchenko, V.; Rustamov, A.; Szczerbak, N. Guide to Amphibians and Reptiles of the USSR Fauna; Prosveshchenie: Moscow, Russia, 1977; p. 414. [Google Scholar]
  117. California Academy of Sciences. Herpetology 32757; California Academy of Sciences, Ed.; California Academy of Sciences: San Francisco, CA, USA, 2019; Volume 2019. [Google Scholar]
  118. Nikolskii, A.M. Rana emeljanovi sp. n. Annu. Musée Zool. L’academie Impériale Sci. St. Pétersbourg 1913, 18, 148–150. [Google Scholar]
  119. Hallowell, E. Report upon the Reptilia of the North Pacific Exploring Expedition, under command of Capt. John Rogers. Proc. Acad. Nat. Sci. Phila. 1861, 12, 480–510. [Google Scholar]
  120. California Academy of Sciences. Herpetology 32875; California Academy of Sciences, Ed.; California Academy of Sciences: San Francisco, CA, USA, 2017; Volume 2019. [Google Scholar]
  121. Boulenger, G. Description of a new tailed batrachian from Corea. J. Nat. Hist. 1887, 14, 67. [Google Scholar] [CrossRef]
  122. Kim, J. Natural monuments of our country. Chosun Nat. Conserv. Prot. Group Cent. Branch 2017, 4, 50–51. [Google Scholar]
  123. Kim, K.; Kim, M.; Kim, R.; Kim, K.; Kim, T.; Rim, C.; Pak, U.; Han, K. Red Data Book of Democratic People’s Republic of Korea; MAB national committee of Democratic People’s Republic of Korea: Pyongyang, Democratic People’s Republic of Korea, 2002. [Google Scholar]
  124. Mori, T. On a new Hynobius from Quelpaert Isl. J. Chosen Nat. Hist. Soc. 1928, 6, 47–52. [Google Scholar]
  125. Nadachowski, A.; Szyndlar, Z.; Tomek, T. Fauna kregowcow ladowych polwyspu Koreanskiego. Prz. Zool. 1982, 26, 79–92. [Google Scholar]
  126. Shin, Y.; Jang, Y.; Allain, S.J.R.; Borzée, A. Catalogue of herpetological specimens of the Ewha Womans University Natural History Museum (EWNHM), Republic of Korea. ZooKeys 2020, 965, 103–139. [Google Scholar] [CrossRef]
  127. Lu, Y.Y.; Wang, J.Q.; Dong, B.J.; Li, P.P. A New Record of Rana huanrenensis in Jilin Province, China. Sichuan J. Zool. 2005, 3, 3. [Google Scholar]
  128. Jiang, X. Thoughts on the development of ecotourism in Laoguiding Nature Reserve of Huanren. For. Prod. Spec. China 2015, 6, 89–91. [Google Scholar]
  129. Malyarchuk, B.; Derenko, M.; Berman, D.; Perkova, M.; Grzybowski, T.; Lejrikh, A.; Bulakhova, N. Phylogeography and molecular adaptation of Siberian salamander Salamandrella keyserlingii based on mitochondrial DNA variation. Mol. Phylogenet. Evol. 2010, 56, 562–571. [Google Scholar] [CrossRef]
  130. Min, M.S.; Yang, S.-Y.; Bonett, R.; Vieites, D.; Brandon, R.; Wake, D. Discovery of the first Asian plethodontid salamander. Nature 2005, 435, 87–90. [Google Scholar] [CrossRef] [PubMed]
  131. Elith, J.; Phillips, S.J.; Hastie, T.; Dudík, M.; Chee, Y.E.; Yates, C.J. A statistical explanation of MaxEnt for ecologists. Divers. Distrib. 2011, 17, 43–57. [Google Scholar] [CrossRef]
  132. Tytar, V.; Nekrasova, O.; Pupina, A.; Pupins, M.; Oskyrko, O. Long-term bioclimatic modelling the distribution of the fire-bellied toad, Bombina bombina (Anura, Bombinatoridae), under the influence of global climate change. Vestn. Zool. 2018, 52, 341–348. [Google Scholar] [CrossRef] [Green Version]
  133. Allouche, O.; Tsoar, A.; Kadmon, R. Assessing the accuracy of species distribution models: Prevalence, kappa and the true skill statistic (TSS). J. Appl. Ecol. 2006, 43, 1223–1232. [Google Scholar] [CrossRef]
  134. Maes, D.; Isaac, N.J.; Harrower, C.A.; Collen, B.; Van Strien, A.J.; Roy, D.B. The use of opportunistic data for IUCN Red List assessments. Biol. J. Linn. Soc. 2015, 115, 690–706. [Google Scholar] [CrossRef] [Green Version]
  135. Lamoreux, J.; Akçakaya, H.R.; Bennun, L.; Collar, N.J.; Boitani, L.; Brackett, D.; Bräutigam, A.; Brooks, T.M.; da Fonseca, G.A.; Mittermeier, R.A. Value of the IUCN red list. Trends Ecol. Evol. 2003, 18, 214–215. [Google Scholar]
  136. Rodrigues, A.S.; Pilgrim, J.D.; Lamoreux, J.F.; Hoffmann, M.; Brooks, T.M. The value of the IUCN Red List for conservation. Trends Ecol. Evol. 2006, 21, 71–76. [Google Scholar] [CrossRef] [PubMed]
  137. Keller, V.; Bollmann, K. From red lists to species of conservation concern. Conserv. Biol. 2004, 18, 1636–1644. [Google Scholar] [CrossRef]
  138. Fitzpatrick, U.; Murray, T.E.; Paxton, R.J.; Brown, M.J. Building on IUCN regional red lists to produce lists of species of conservation priority: A model with Irish bees. Conserv. Biol. 2007, 21, 1324–1332. [Google Scholar] [CrossRef] [Green Version]
  139. Simaika, J.P.; Samways, M.J. Reserve selection using Red Listed taxa in three global biodiversity hotspots: Dragonflies in South Africa. Biol. Conserv. 2009, 142, 638–651. [Google Scholar] [CrossRef]
  140. Butchart, S.H.; Akcakaya, H.R.; Kennedy, E.; Hilton-Taylor, C. Biodiversity indicators based on trends in conservation status: Strengths of the IUCN Red List Index. Conserv. Biol. 2006, 20, 579–581. [Google Scholar] [CrossRef]
  141. Mace, G.M.; Collar, N.J.; Gaston, K.J.; Hilton-Taylor, C.; Akçakaya, H.R.; Leader-Williams, N.; Milner-Gulland, E.J.; Stuart, S.N. Quantification of extinction risk: IUCN’s system for classifying threatened species. Conserv. Biol. 2008, 22, 1424–1442. [Google Scholar] [CrossRef]
  142. Gärdenfors, U.; Hilton-Taylor, C.; Mace, G.M.; Rodríguez, J.P. The application of IUCN Red List criteria at regional levels. Conserv. Biol. 2001, 15, 1206–1212. [Google Scholar] [CrossRef]
  143. Miller, R.M.; Rodríguez, J.P.; Fowler, T.A.; Bambaradeniya, C.; Boles, R.; Eaton, M.A.; Gärdenfors, U.; Keller, V.; Molur, S.; Walker, S. National threatened species listing based on IUCN criteria and regional guidelines: Current status and future perspectives. Conserv. Biol. 2007, 21, 684–696. [Google Scholar] [CrossRef]
  144. IUCN. Guidelines for Application of IUCN Red List Criteria at Regional and National Levels: Version 4.0; IUCN Species Survival Commission (SSC): Gland, Switzerland, 2012. [Google Scholar]
  145. Mace, G.M. Classifying threatened species: Means and ends. Phil. Trans. R. Soc. Lond. B 1994, 344, 91–97. [Google Scholar]
  146. Hermoso, V.; Kennard, M.J.; Linke, S. Evaluating the costs and benefits of systematic data acquisition for conservation assessments. Ecography 2015, 38, 283–292. [Google Scholar] [CrossRef]
  147. IUCN. International Union for Conservation of Nature Red List; IUCN, Ed.; IUCN: Gland, Switzerland, 2018. [Google Scholar]
  148. Kim, Y.I. Intraspecific Latitudinal Variation in Critical Thermal Maximum of Tadpoles in Two Anuran Species; Ewha Womans University: Seoul, Korea, 2018. [Google Scholar]
  149. Heo, K.; Kim, Y.I.; Bae, Y.; Jang, Y.; Amaël, B. First report of Dryophytes japonicus tadpoles in saline environment. Russ. J. Herpetol. 2019, 26, 87–90. [Google Scholar] [CrossRef]
  150. Borzée, A.; Purevdorj, Z.; Kim, Y.I.; Kong, S.; Choe, M.; Yi, Y.; Kim, K.; Kim, A.; Jang, Y. Breeding preferences in the Treefrogs Dryophytes japonicus (Hylidae) in Mongolia. J. Nat. Hist. 2019, 53, 43–44. [Google Scholar]
  151. Gosner, K.L. A simplified table for staging anuran embryos and larvae with notes on identification. Herpetologica 1960, 16, 183–190. [Google Scholar]
  152. Matsui, M. Pelophylax chosenicus. In The IUCN Red List Threatened Species; IUCN: Gland, Switzerland, 2004; Volume e.T58577A11806007. [Google Scholar] [CrossRef]
  153. Lee, S.-D.; Miller-Rushing, A.J. Degradation, urbanization, and restoration: A review of the challenges and future of conservation on the Korean Peninsula. Biol Cons 2014, 176, 262–276. [Google Scholar] [CrossRef]
  154. Borzée, A.; Struecker, M.-Y.; Yi, Y.; Kim, D.; Kim, H. Time for Korean wildlife conservation. Science 2019, 363, 1161–1162. [Google Scholar]
  155. The World Bank. World Development Indicators; The World Bank: Washington, DC, USA, 2013. [Google Scholar]
  156. Stone, R. Seeking cures for North Korea’s environmental ills. Science 2012, 335, 1425–1426. [Google Scholar] [CrossRef] [Green Version]
  157. UNEP. United Nations Environment Programme; Piddington, K., Ed.; UNEP: Klong Luang, Thailand, 2003. [Google Scholar]
  158. Democratic People’s Republic of Korea. National Biodiversity Strategy and Action Plan of DPR Korea; Democratic People’s Republic of Korea: Pyongyang, Democratic People’s Republic of Korea, 2007. [Google Scholar]
  159. Pimm, S.L.; Jenkins, C.N.; Abell, R.; Brooks, T.M.; Gittleman, J.L.; Joppa, L.N.; Raven, P.H.; Roberts, C.M.; Sexton, J.O. The biodiversity of species and their rates of extinction, distribution, and protection. Science 2014, 344, 987–997. [Google Scholar] [CrossRef]
  160. IUCN SSC Amphibian Specialist Group. Dryophytes Suweonensis. In The IUCN Red List of Threatened Species; IUCN SSC Amphibian Specialist Group; IUCN: Gland, Switzerland, 2017. [Google Scholar] [CrossRef]
  161. Kuzmin, S.; Pipeng, L.; Matsui, M.; Ishchenko, V.; Maslova, I. Bombina orientalis. In The IUCN Red List of Threatened Species; IUCN: Gland, Switzerland, 2004; Volume e.T54449A11146991. [Google Scholar] [CrossRef]
  162. Kuzmin, S.; Yuezhao, W.; Matsui, M.; Kaneko, Y.; Maslova, I. Bufo gargarizans. In The IUCN Red List of Threatened Species; IUCN: Gland, Switzerland, 2004; Volume e.T54647A11180910. [Google Scholar]
  163. Matsui, M.; Wenge, Z. Bufo stejnegeri. In The IUCN Red List of Threatened Species; IUCN: Gland, Switzerland, 2004; Volume e.T54766A11200818. [Google Scholar] [CrossRef]
  164. Koo, K.-S. Distribution Characteristics and Body Size Variation of Japanese Tree Frog (Hyla japonica) in South Korea; Gangwon National University: Chuncheon, Korea, 2018. [Google Scholar]
  165. Borzée, A. Why Are Anurans Threatened? The Case of Dryophytes suweonensis. Ph.D Thesis, Seoul National University, Seoul, Korea, 2018. [Google Scholar]
  166. Borzée, A.; Fong, J.J.; Nguyen, H.; Jang, Y. Large-scale hybridisation as an extinction threat to the Suweon treefrog (Hylidae: Dryophytes suweonensis). Animals 2020, 10, 764. [Google Scholar] [CrossRef] [PubMed]
  167. Groffen, J.; Kong, S.; Jang, Y.; Borzée, A. The invasive American bullfrog (Lithobates catesbeianus) in the Republic of Korea: History and recommendation for population control. Manag. Biol. Invasions 2019, 10, 517–535. [Google Scholar] [CrossRef] [Green Version]
  168. Borzée, A.; Kosch, T.A.; Kim, M.; Jang, Y. Introduced bullfrogs are associated with increased Batrachochytrium dendrobatidis prevalence and reduced occurrence of Korean treefrogs. PLoS ONE 2017, 12, e0177860. [Google Scholar]
  169. IUCN SSC Amphibian Specialist Group. Dryophytes suweonensis (amended version of 2014 assessment). In The IUCN Red List of Threatened Species; IUCN: Gland, Switzerland, 2017; Volume e.T55670A112715252. [Google Scholar] [CrossRef]
  170. Kuzmin, S.; Matsui, M.; Wenge, Z.; Kaneko, Y. Glandirana emeljanovi. In The IUCN Red List of Threatened Species; IUCN: Gland, Switzerland, 2004; Volume e.T58591A11793116. [Google Scholar] [CrossRef]
  171. Matsui, M.; Wenge, Z. Hynobius leechii. In The IUCN Red List of Threatened Species; IUCN: Gland, Switzerland, 2008; Volume e.T59095A11870831. [Google Scholar]
  172. Matsui, M.; Wenge, Z. The IUCN Red List of Threatened Species; IUCN: Gland, Switzerland, 2004; Volume e.T57849A11693398. [Google Scholar]
  173. Kuzmin, S.; Matsui, M.; Wenge, Z.; Maslova, I.; Tuniyev, B. Onychodactylus fischeri. In The IUCN Red List of Threatened Species; IUCN: Gland, Switzerland, 2004; Volume e.T39419A10234504. [Google Scholar] [CrossRef]
  174. Kuzmin, S.; Maslova, I.; Tuniyev, B.; Matsui, M.; Pipeng, L.; Kaneko, Y. Pelophylax nigromaculatus. In The IUCN Red List of Threatened Species; IUCN: Gland, Switzerland, 2004; Volume e.T58679A11809026. [Google Scholar]
  175. Litvinchuk, S.N.; Schepina, N.A.; Borzée, A. Reconstruction of past distribution for the Mongolian toad, Strauchbufo raddei (Anura: Bufonidae) using environmental modeling. PeerJ 2020, 8, e9216. [Google Scholar] [CrossRef] [PubMed]
  176. Kuzmin, S.; Matsui, M.; Gang, L.; Maslova, I. Strauchbufo raddei. In The IUCN Red List of Threatened Species; IUCN: Gland, Switzerland, 2004; Volume e.T54744A11198066. [Google Scholar]
  177. Kuzmin, S.; Maslova, I.; Matsui, M.; Wenge, Z. Rana amurensis. In The IUCN Red List of Threatened Species; IUCN: Gland, Switzerland, 2008; Volume e.T58542A11799162. [Google Scholar] [CrossRef]
  178. Song, J.-Y. Rana coreana. In The IUCN Red List of Threatened Species; IUCN: Gland, Switzerland, 2008; Volume e.T135860A4212618. [Google Scholar] [CrossRef]
  179. Kuzmin, S.; Ishchenko, V.; Maslova, I.; Ananjeva, N.; Orlov, N.; Matsui, M.; Feng, X.; Kaneko, Y. Rana dybowskii. In The IUCN Red List of Threatened Species; IUCN: Gland, Switzerland, 2004; Volume e.T58589A11792510. [Google Scholar] [CrossRef]
  180. Matsui, M.; Feng, X. Rana huanrensis. In The IUCN Red List of Threatened Species; IUCN: Gland, Switzerland, 2004; Volume e.T58619A11812585. [Google Scholar] [CrossRef]
  181. Fei, L.; Changyuan, Y.; Jianping, J. Colored Atlas of Chinese Amphibians and Their Distributions; Sichuan science and technology Press: Chendu, China, 2012. [Google Scholar]
  182. Korotkov, Y.M. To the Biology of Bombina Orientalis, Agkistrodon Blomhoffii and A. halis in Primorye Region; Zoologicheskie Problemy Sibiri: Novosibirsk, Russia, 1972; Volume 302. [Google Scholar]
  183. Korotkov, Y.M.; Korotkova, E.B. Some patterns of reproduction of Rana chensinensis in the Far East marine territory. Zool. Zhurnal 1977, 56, 1057–1061. [Google Scholar]
  184. Piao, R.Z. List of Wildlife in Heilongjiang China; The institute of wildlife in Heilongjiang province China: Harbin, China, 1990. [Google Scholar]
  185. Fok, K.W.; Wade, C.M.; Parkin, D.T. Inferring the phylogeny of disjunct populations of the azure–winged magpie Cyanopica cyanus from mitochondrial control region sequences. Proc. R. Soc. Lond. B Biol. Sci. 2002, 269, 1671–1679. [Google Scholar] [CrossRef] [Green Version]
  186. Zanini, F.; Klingemann, A.; Schlaepfer, R.; Schmidt, B.R. Landscape effects on anuran pond occupancy in an agricultural countryside: Barrier-based buffers predict distributions better than circular buffers. Can. J. Zool. 2008, 86, 692–699. [Google Scholar] [CrossRef] [Green Version]
  187. Ruiz, G.M.; Hines, A.H.; Jivoff, P. Biotic resistance to invasion: Native predator limits abundance and distribution of an introduced crab. Ecology 2005, 86, 3364–3376. [Google Scholar]
  188. Marco, A.; Quilchano, C.; Blaustein, A.R. Sensitivity to nitrate and nitrite in pond-breeding amphibians from the Pacific Northwest, USA. Environ. Toxicol. Chem. 1999, 18, 2836–2839. [Google Scholar] [CrossRef]
  189. Macias, D.A.; Groffen, J.; Jang, Y.; Borzée, A. Rana coreana (Korean Brown Frog) and R. uenoi (Ueno’s Brown Frog). Hibernaculum. Herpetol. Rev. 2018, 49, 121–122. [Google Scholar]
  190. Borzée, A.; Choi, Y.; Kim, Y.E.; Jablonski, P.G.; Jang, Y. Interspecific variation in seasonal migration and brumation behavior in two closely related species of treefrogs. Front. Ecol. Evol. 2019, 7, 55. [Google Scholar] [CrossRef] [Green Version]
  191. Wiens, J.J. Speciation and ecology revisited: Phylogenetic niche conservatism and the origin of species. Evol. Int. J. Org. Evol. 2004, 58, 193–197. [Google Scholar] [CrossRef]
  192. Vences, M.; Wake, D.B. Amphibian Biology; Heatwole, H., Tyler-Michael, J., Eds.; Surrey Beatty and Sons: Norton, Australia, 2007; Volume 7, pp. 2613–2671. [Google Scholar]
  193. Hendry, A.P. Ecological speciation! Or the lack thereof? Can. J. Fish. Aquat. Sci. 2009, 66, 1383–1398. [Google Scholar] [CrossRef] [Green Version]
  194. Hua, X.; Wiens, J.J. How does climate influence speciation? Am. Nat. 2013, 182, 1–12. [Google Scholar] [CrossRef] [Green Version]
  195. Jo, Y.-S.; Baccus, J.T.; Koprowski, J. Mammals of Korea; National Institute of Biological Resources: Incheon, Korea, 2018. [Google Scholar]
  196. Brazil, M. Birds of East Asia: China, Taiwan, Korea, Japan, and Russia; Princeton University Press: Trenton, NJ, USA, 2009; p. 528. [Google Scholar]
  197. Hewitt, G. The genetic legacy of the Quaternary ice ages. Nature 2000, 405, 907–913. [Google Scholar] [CrossRef]
  198. Veith, M.; Kosuch, J.; Vences, M. Climatic oscillations triggered post-Messinian speciation of Western Palearctic brown frogs (Amphibia, Ranidae). Mol. Phylogenet. Evol. 2003, 26, 310–327. [Google Scholar] [CrossRef]
  199. Avise, J.C. Phylogeography of Southern European Refugia; Weiss, S.F.N.E., Ed.; Phylogeography of southern European refugia; Springer: New York, NY, USA, 2007; p. 377. [Google Scholar]
  200. Knowles, L.L. Did the Pleistocene glaciations promote divergence? Tests of explicit refugial models in montane grasshopprers. Mol. Ecol. 2001, 10, 691–701. [Google Scholar] [CrossRef]
  201. Kent, M. Advanced Biology; Oxford University Press: New York, NY, USA, 2000. [Google Scholar]
  202. Nosil, P. Ecological Speciation; Oxford University Press: New York, NY, USA, 2012. [Google Scholar]
  203. Aizawa, M.; Kim, Z.-S.; Yoshimaru, H. Phylogeography of the Korean pine (Pinus koraiensis) in northeast Asia: Inferences from organelle gene sequences. J. Plant Res. 2012, 125, 713–723. [Google Scholar] [CrossRef] [PubMed]
  204. Kim, S.I.; Park, S.K.; Lee, H.; Oshida, T.; Kimura, J.; Kim, Y.J.; Nguyen, S.; Sashika, M.; Min, M.S. Phylogeography of Korean raccoon dogs: Implications of peripheral isolation of a forest mammal in East Asia. J. Zool. 2013, 290, 225–235. [Google Scholar] [CrossRef]
  205. Davis, M.B.; Shaw, R.G. Range shifts and adaptive responses to Quaternary climate change. Science 2001, 292, 673–679. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  206. Provan, J.; Bennett, K. Phylogeographic insights into cryptic glacial refugia. Trends Ecol. Evol. 2008, 23, 564–571. [Google Scholar] [CrossRef] [PubMed]
  207. Stöck, M.; Dufresnes, C.; Litvinchuk, S.N.; Lymberakis, P.; Biollay, S.; Berroneau, M.; Borzée, A.; Ghali, K.; Ogielska, M.; Perrin, N. Cryptic diversity among Western Palearctic tree frogs: Postglacial range expansion, range limits, and secondary contacts of three European tree frog lineages (Hyla arborea group). Mol. Phylogenet. Evol. 2012, 65, 1–9. [Google Scholar]
  208. Zhang, H.; Yan, J.; Zhang, G.; Zhou, K. Phylogeography and demographic history of Chinese black-spotted frog populations (Pelophylax nigromaculata): Evidence for independent refugia expansion and secondary contact. BMC Evol. Biol. 2008, 8, 21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  209. Borzée, A.; Santos, J.L.; Sanchez-Ramirez, S.; Bae, Y.; Heo, K.; Jang, Y.; Jowers, M.J. Phylogeographic and population insights of the Asian common toad (Bufo gargarizans) in Korea and China: Population isolation and expansions as response to the ice ages. PeerJ 2017, 5, e4044. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  210. Lee, J.E.; Yang, S.Y.; Lee, H.Y. Genetic relationships among six Korean Rana species (Amphibia; Ranidae) based on the mitochondrial cytochrome b gene. Korean J. Biol. Sci. 2000, 4, 117–121. [Google Scholar] [CrossRef]
  211. Qian, H.; Ricklefs, R.E. Large-scale processes and the Asian bias in species diversity of temperate plants. Nature 2000, 407, 180. [Google Scholar] [CrossRef] [PubMed]
  212. Serizawa, K.; Suzuki, H.; Iwasa, M.A.; Tsuchiya, K.; Pavlenko, M.V.; Kartavtseva, I.V.; Chelomina, G.N.; Dokuchaev, N.E.; Han, S.-H. A spatial aspect on mitochondrial DNA genealogy in Apodemus peninsulae from East Asia. Biochem. Genet. 2002, 40, 149–161. [Google Scholar] [CrossRef]
  213. Aizawa, M.; Yoshimaru, H.; Saito, H.; Katsuki, T.; Kawahara, T.; Kitamura, K.; Shi, F.; Kaji, M. Phylogeography of a northeast Asian spruce, Picea jezoensis, inferred from genetic variation observed in organelle DNA markers. Mol. Ecol. 2007, 16, 3393–3405. [Google Scholar] [CrossRef]
  214. Bai, W.N.; Liao, W.J.; Zhang, D.Y. Nuclear and chloroplast DNA phylogeography reveal two refuge areas with asymmetrical gene flow in a temperate walnut tree from East Asia. New Phytol. 2010, 188, 892–901. [Google Scholar] [CrossRef]
  215. Ding, L.; Gan, X.N.; He, S.P.; Zhao, E.M. A phylogeographic, demographic and historical analysis of the short-tailed pit viper (Gloydius brevicaudus): Evidence for early divergence and late expansion during the Pleistocene. Mol. Ecol. 2011, 20, 1905–1922. [Google Scholar] [CrossRef]
  216. Othman, S.N.; Putri, E.T.; Messenger, K.R.; Bae, Y.; Yang, Y.; Bova, T.; Reed, T.; Amin, H.; Chuang, M.-F.; Jang, Y.; et al. Impact of south-central Asian orogenesis on Kaloula spp. (Anura: Microhylidae) radiation and implication of Pleistocene glaciation on population extension to northern latitudes. Integr. Zool. 2021, (in press). [Google Scholar]
  217. Liu, K.; Wang, F.; Chen, W.; Tu, L.; Min, M.-S.; Bi, K.; Fu, J. Rampant historical mitochondrial genome introgression between two species of green pond frogs, Pelophylax nigromaculatus and P. plancyi. BMC Evol. Biol. 2010, 10, 201. [Google Scholar]
  218. Komaki, S.; Igawa, T.; Lin, S.M.; Tojo, K.; Min, M.S.; Sumida, M. Robust molecular phylogeny and palaeodistribution modelling resolve a complex evolutionary history: Glacial cycling drove recurrent mtDNA introgression among Pelophylax frogs in East Asia. J. Biogeogr. 2015, 42, 2159–2171. [Google Scholar] [CrossRef]
  219. Fairbanks, R.G. A 17,000-year glacio-eustatic sea level record: Influence of glacial melting rates on the Younger Dryas event and deep-ocean circulation. Nature 1989, 342, 637–642. [Google Scholar] [CrossRef]
  220. Park, Y.; Khim, B.; Zhao, S. Sea level fluctuation in the Yellow Sea Basin. J. Korean Soc. Oceanogr. 1994, 29, 42–49. [Google Scholar]
  221. Kim, J.-M.; Kennett, J.P. Paleoenvironmental changes associated with the Holocene marine transgression, Yellow Sea (Hwanghae). Mar. Micropaleontol. 1998, 34, 71–89. [Google Scholar] [CrossRef]
  222. Ijiri, A.; Wang, L.; Oba, T.; Kawahata, H.; Huang, C.-Y.; Huang, C.-Y. Paleoenvironmental changes in the northern area of the East China Sea during the past 42,000 years. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2005, 219, 239–261. [Google Scholar] [CrossRef]
  223. Nguyen, T.Q.; Bernardes, M.; Nguyen, T.T.; Ziegler, T. First records of Bufo gargarizans Cantor, 1842 and Odorrana lipuensis Mo, Chen, Wu, Zhang et Zhou, 2015 (Anura: Bufonidae, Ranidae) from Vietnam. Russ. J. Herpetol. 2016, 23, 103–107. [Google Scholar]
  224. Liu, W.; Lathrop, A.; Fu, J.; Yang, D.; Murphy, R.W. Phylogeny of East Asian bufonids inferred from mitochondrial DNA sequences (Anura: Amphibia). Mol. Phylogenet. Evol. 2000, 14, 423–435. [Google Scholar] [CrossRef] [Green Version]
  225. Fu, J.; Weadick, C.J.; Zeng, X.; Wang, Y.; Liu, Z.; Zheng, Y.; Li, C.; Hu, Y. Phylogeographic analysis of the Bufo gargarizans species complex: A revisit. Mol. Phylogenet. Evol. 2005, 37, 202–213. [Google Scholar] [CrossRef]
  226. Zhan, A.; Fu, J. Past and present: Phylogeography of the Bufo gargarizans species complex inferred from multi-loci allele sequence and frequency data. Mol Phylogenet. Evol 2011, 61, 136–148. [Google Scholar] [CrossRef] [PubMed]
  227. Garcia-Porta, J.; Litvinchuk, S.; Crochet, P.-A.; Romano, A.; Geniez, P.; Lo-Valvo, M.; Lymberakis, P.; Carranza, S. Molecular phylogenetics and historical biogeography of the west-palearctic common toads (Bufo bufo species complex). Mol. Phylogenet. Evol. 2012, 63, 113–130. [Google Scholar] [CrossRef]
  228. Tong, H.; Wo, Y.; Lioa, P.; Jin, Y. Phylogenetic, demographic and dating analyses of Bufo gargarizans populations from the Zhoushan archipelago and mainland China. Asian Herpetol. Res. 2017, 8, 165–173. [Google Scholar]
  229. Kong, W.S. Vegetational history of the Korean Peninsula. Glob. Ecol. Biogeogr. 2000, 9, 391–402. [Google Scholar] [CrossRef]
  230. Yan, F.; Zhou, W.; Zhao, H.; Yuan, Z.; Wang, Y.; Jiang, K.; Jin, J.; Murphy, R.W.; Che, J.; Zhang, Y. Geological events play a larger role than Pleistocene climatic fluctuations in driving the genetic structure of Quasipaa boulengeri (Anura: Dicroglossidae). Mol. Ecol. 2013, 22, 1120–1133. [Google Scholar] [CrossRef] [PubMed]
  231. Ryu, E.; Lee, S.-J.; Yang, D.-Y.; Kim, J.-Y. Paleoenvironmental studies of the Korean peninsula inferred from diatom assemblages. Quat. Int. 2008, 176, 36–45. [Google Scholar] [CrossRef]
  232. He, L.; Mukai, T.; Chu, K.; Ma, Q.; Zhang, J. Biogeographical role of the Kuroshio Current in the amphibious mudskipper Periophthalmus modestus indicated by mitochondrial DNA data. Sci. Rep. 2015, 5, 15645. [Google Scholar] [CrossRef] [Green Version]
  233. Wang, B.; Lin, H. Rainy season of the asian-pacific summer monsoon. J. Clim. 2002, 15, 386–398. [Google Scholar] [CrossRef] [Green Version]
  234. Lee, E.J.; Jhun, J.G.; Park, C.K. Remote connection of the northeast Asian summer rainfall variation revealed by a newly defined monsoon index. J. Clim. 2005, 18, 4381–4393. [Google Scholar] [CrossRef]
  235. Beals, E.W. Vegetational change along altitudinal gradients. Science 1969, 165, 981–985. [Google Scholar] [CrossRef]
  236. Stevens, G.C. The elevational gradient in altitudinal range: An extension of Rapoport’s latitudinal rule to altitude. Am. Nat. 1992, 140, 893–911. [Google Scholar] [CrossRef] [PubMed]
  237. Rosenzweig, M.L. Species Diversity in Space and Time; Cambridge University Press: Cambridge, UK, 1995. [Google Scholar]
  238. Gaston, K.J. Global patterns in biodiversity. Nature 2000, 405, 220. [Google Scholar] [CrossRef] [PubMed]
  239. Fischer, A.; Blaschke, M.; Bässler, C. Altitudinal gradients in biodiversity research: The state of the art and future perspectives under climate change aspects. For. Ecol. Landsc. Res. Nat. Conserv. 2011, 11, 35–47. [Google Scholar]
  240. Willig, M.R.; Presley, S.J. Biodiversity and metacommunity structure of animals along altitudinal gradients in tropical montane forests. J. Trop. Ecol. 2016, 32, 421–436. [Google Scholar] [CrossRef]
  241. Loveridge, A. Scientific Results of an Expedition to Rain Forest Regions in Eastern Africa: Zoogeography and Itinerary. Bull. Mus. Comp. Zool. 1937, 79, 481–541. [Google Scholar]
  242. Daniels, R.R. Geographical distribution patterns of amphibians in the Western Ghats, India. J. Biogeogr. 1992, 19, 521–529. [Google Scholar] [CrossRef]
  243. Poynton, J.; Loader, S.; Sherratt, E.; Clarke, B. Amphibian diversity in East African biodiversity hotspots: Altitudinal and latitudinal patterns. Biodivers. Conserv. 2006, 16, 277–292. [Google Scholar]
  244. Cogălniceanu, D.; Székely, P.; Samoilă, C.; Ruben, I.; Tudor, M.; Plăiaşu, R.; Stănescu, F.; Rozylowicz, L. Diversity and distribution of amphibians in Romania. ZooKeys 2013, 296, 35. [Google Scholar]
  245. Gaston, K.J.; Fuller, R.A. The sizes of species’ geographic ranges. J. Appl. Ecol. 2009, 46, 1–9. [Google Scholar] [CrossRef]
  246. Dodd, C.K., Jr. Biological diversity of a temporary pond herpetofauna in north Florida sandhills. Biodivers. Conserv. 1992, 1, 125–142. [Google Scholar] [CrossRef]
  247. Semlitsch, R.D. Biological Delineation of Terrestrial Buffer Zones for Pond-Breeding Salamanders. Conserv. Biol. 1998, 12, 1113–1119. [Google Scholar] [CrossRef]
  248. Natuhara, Y. Ecosystem services by paddy fields as substitutes of natural wetlands in Japan. Ecol. Eng. 2013, 56, 97–106. [Google Scholar] [CrossRef]
  249. Holzer, K.A. Amphibian use of constructed and remnant wetlands in an urban landscape. Urban Ecosyst. 2014, 17, 955–968. [Google Scholar] [CrossRef]
  250. Quesnelle, P.E.; Lindsay, K.E.; Fahrig, L. Relative effects of landscape-scale wetland amount and landscape matrix quality on wetland vertebrates: A meta-analysis. Ecol. Appl. 2015, 25, 812–825. [Google Scholar] [CrossRef]
  251. Bennett, A.F. Linkages in the Landscape: The Role of Corridors and Connectivity in Wildlife Conservation; IUCN: Gland, Switzerland, 1998. [Google Scholar]
  252. Forman, R.T.; Alexander, L.E. Roads and their major ecological effects. Annu. Rev. Ecol. Syst. 1998, 29, 207–231. [Google Scholar] [CrossRef] [Green Version]
  253. Geneletti, D. Using spatial indicators and value functions to assess ecosystem fragmentation caused by linear infrastructures. Int. J. Appl. Earth Obs. Geoinf. 2004, 5, 1–15. [Google Scholar] [CrossRef]
  254. Moriyama, H. What Does it Mean to Protect Paddy Fields; Rural Culture Association Japan: Tokyo, Japan, 1997. [Google Scholar]
  255. Machado, I.F.; Maltchik, L. Can management practices in rice fields contribute to amphibian conservation in southern Brazilian wetlands? Aquat. Conserv. Mar. Freshw. Ecosyst. 2010, 20, 39–46. [Google Scholar] [CrossRef]
  256. Naito, R. Perspectives of conservation of pond-breeding frogs (focusing on the Nagoya Daruma pond frog) in rice paddy areas in Japan. Landsc. Ecol. 2012, 17, 57–73. [Google Scholar]
  257. Naito, R.; Sakai, M.; Morimoto, Y. Negative effects of deep roadside ditches on Pelophylax porosa brevipoda dispersal and migration in comparison with Hyla japonica in a rice paddy area in Japan. Zool. Sci. 2012, 29, 599–603. [Google Scholar] [CrossRef]
  258. Naito, R.; Sakai, M.; Natuhara, Y.; Morimoto, Y.; Shibata, S. Microhabitat use by Hyla japonica and Pelophylax porosa brevipoda at levees in rice paddy areas of Japan. Zool. Sci. 2013, 30, 386–391. [Google Scholar] [CrossRef]
  259. Borzée, A.; Jang, Y. Impact of rice and bean harvests on the Suweon Treefrog (Dryophytes suweonensis). Int. J. Curr. Res. 2017, 9, 59620–59623. [Google Scholar]
  260. Borzée, A.; Heo, K.; Jang, Y. Relationship between agroenvironmental variables and breeding Hylids in rice paddies. Sci. Rep. 2018, 8, 8049. [Google Scholar] [CrossRef] [PubMed]
  261. Groffen, J.; Borzée, A.; Jang, Y. Positioning of two treefrog species within rice paddies in relation to different habitat borders. Anim. Cells Syst. 2018, 22, 205–211. [Google Scholar] [CrossRef] [Green Version]
  262. Borzée, A.; Jang, Y. Policy recommendation for the conservation of the Suweon Treefrog (Dryophytes suweonensis) in the Republic of Korea. Front. Environ. Sci. 2019, 7, 39. [Google Scholar] [CrossRef] [Green Version]
  263. Orchard, D.; Tessa, G.; Jehle, R. Age and growth in a European flagship amphibian: Equal performance at agricultural ponds and favourably managed aquatic sites. Aquat. Ecol. 2019, 53, 37–48. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Map of survey locations and geocoded points for amphibian species presence in the Democratic People’s Republic of Korea. See the Supplementary Materials Tables for exact GPS coordinates for species in DPR Korea and Supplementary Materials Table S1 for all georeferenced datapoints for the species present in DPR Korea. The layer with ecoregions is attributed to [55]. Map created in ArcGIS 10.5 (ESRI; Redlands, CA, USA).
Figure 1. Map of survey locations and geocoded points for amphibian species presence in the Democratic People’s Republic of Korea. See the Supplementary Materials Tables for exact GPS coordinates for species in DPR Korea and Supplementary Materials Table S1 for all georeferenced datapoints for the species present in DPR Korea. The layer with ecoregions is attributed to [55]. Map created in ArcGIS 10.5 (ESRI; Redlands, CA, USA).
Animals 11 02057 g001
Figure 2. Pictures of the holotype of Bufo sambangensis presented by [9]. Note the morphological similarities with B. stejnegeri: (A) detailed pictures of the head; (B) detailed picture of the hind limbs; (C) dorsal picture. Pictures taken by Mr. Tu Yong Nam from Institute of Zoology of the State Academy of Science in the Democratic People’s Republic of Korea.
Figure 2. Pictures of the holotype of Bufo sambangensis presented by [9]. Note the morphological similarities with B. stejnegeri: (A) detailed pictures of the head; (B) detailed picture of the hind limbs; (C) dorsal picture. Pictures taken by Mr. Tu Yong Nam from Institute of Zoology of the State Academy of Science in the Democratic People’s Republic of Korea.
Animals 11 02057 g002
Figure 3. Simplified maximum likelihood phylogenetic tree built with the PHYML plugin in Geneious for Dryophytes sp. samples. The results clarify the presence of D. japonicus in the vicinity of Pyongyang, DPR Korea. Hyla orientalis is used as outgroup in this analysis (GenBank accession numbers shown in Table 2). Branch distances represent nucleotide substitution rates and the scale bar represents the number of changes per nucleotide position. Samples originating from this study are shown in bold.
Figure 3. Simplified maximum likelihood phylogenetic tree built with the PHYML plugin in Geneious for Dryophytes sp. samples. The results clarify the presence of D. japonicus in the vicinity of Pyongyang, DPR Korea. Hyla orientalis is used as outgroup in this analysis (GenBank accession numbers shown in Table 2). Branch distances represent nucleotide substitution rates and the scale bar represents the number of changes per nucleotide position. Samples originating from this study are shown in bold.
Animals 11 02057 g003
Figure 4. Simplified maximum likelihood phylogenetic tree built with the PHYML plugin in Geneious for Rana sp. samples. The results clarify the presence of R. uenoi in the vicinity of Pyongyang, DPR Korea (GenBank accession numbers shown in Table 2). Branch distances represent nucleotide substitution rates and the scale bar represents the number of changes per nucleotide position. The species identification for individual Rana dybowskii from R Korea was arbitrarily reassigned to R. uenoi based on the study by Matsui (2014). These individuals are indicated by *. Samples originating from this study are shown in bold.
Figure 4. Simplified maximum likelihood phylogenetic tree built with the PHYML plugin in Geneious for Rana sp. samples. The results clarify the presence of R. uenoi in the vicinity of Pyongyang, DPR Korea (GenBank accession numbers shown in Table 2). Branch distances represent nucleotide substitution rates and the scale bar represents the number of changes per nucleotide position. The species identification for individual Rana dybowskii from R Korea was arbitrarily reassigned to R. uenoi based on the study by Matsui (2014). These individuals are indicated by *. Samples originating from this study are shown in bold.
Animals 11 02057 g004
Figure 5. Pictures of Rana coreana or Rana amurensis (A) and Dryophytes suweonensis (B) used for species identification in the vicinity of Pyongyang. The samples are stored at the ZISP (Zoological Institute of the Russian Academy of Sciences, St. Petersburg, Russia).
Figure 5. Pictures of Rana coreana or Rana amurensis (A) and Dryophytes suweonensis (B) used for species identification in the vicinity of Pyongyang. The samples are stored at the ZISP (Zoological Institute of the Russian Academy of Sciences, St. Petersburg, Russia).
Animals 11 02057 g005
Figure 6. Location of the sampling sites in Rason, North Hamgyong, DPR Korea. Dark blue pins indicate agricultural wetlands, green pins indicate natural wetlands, white pins indicate grassy hills and light blue pins indicate mountainous areas. Satellite views extracted from Google Earth Pro (Google Earth Pro imagery, v7.1.2.2041; Mountain View, CA, USA).
Figure 6. Location of the sampling sites in Rason, North Hamgyong, DPR Korea. Dark blue pins indicate agricultural wetlands, green pins indicate natural wetlands, white pins indicate grassy hills and light blue pins indicate mountainous areas. Satellite views extracted from Google Earth Pro (Google Earth Pro imagery, v7.1.2.2041; Mountain View, CA, USA).
Animals 11 02057 g006
Figure 7. Spectral call properties for Dryophytes suweonensis, recorded in Mundeok, DPR Korea, in May 2018. (A) Spectrogram with oscillogram underneath a representative call of D. suweonensis. The spectrogram highlights the call between the orange vertical dashed lines. (B) Frequency tracking of the same D. suweonensis call. Calls detected are highlighted between red dashed lines. Fundamental frequency measurements are shown as purple triangles and dominant frequency measurements as yellow crosses.
Figure 7. Spectral call properties for Dryophytes suweonensis, recorded in Mundeok, DPR Korea, in May 2018. (A) Spectrogram with oscillogram underneath a representative call of D. suweonensis. The spectrogram highlights the call between the orange vertical dashed lines. (B) Frequency tracking of the same D. suweonensis call. Calls detected are highlighted between red dashed lines. Fundamental frequency measurements are shown as purple triangles and dominant frequency measurements as yellow crosses.
Animals 11 02057 g007
Figure 8. Amphibian species richness for all amphibian species recorded in DPR Korea. This map combines MaxEnt models (Figure 9) using a sum of presence probabilities. Map was created in ArcGIS 10.5 (ESRI; Redlands, CA, USA).
Figure 8. Amphibian species richness for all amphibian species recorded in DPR Korea. This map combines MaxEnt models (Figure 9) using a sum of presence probabilities. Map was created in ArcGIS 10.5 (ESRI; Redlands, CA, USA).
Animals 11 02057 g008
Figure 9. MaxEnt models for all amphibian species recorded in DPR Korea. Maps were created in ArcGIS 10.5 (ESRI; Redlands, CA, USA).
Figure 9. MaxEnt models for all amphibian species recorded in DPR Korea. Maps were created in ArcGIS 10.5 (ESRI; Redlands, CA, USA).
Animals 11 02057 g009
Table 1. Amphibian collection at the Zoological Institute of the Russian Academy of Sciences, St. Petersburg, Russia, and the Zoological Museum of Moscow State University, Moscow, Russia, originating from DPR Korea and collected by the Chosun Academy of Science (DPR Korea Academy of Sciences). These datapoints are not reported later in the text in case of overlap with known combinations of species and localities.
Table 1. Amphibian collection at the Zoological Institute of the Russian Academy of Sciences, St. Petersburg, Russia, and the Zoological Museum of Moscow State University, Moscow, Russia, originating from DPR Korea and collected by the Chosun Academy of Science (DPR Korea Academy of Sciences). These datapoints are not reported later in the text in case of overlap with known combinations of species and localities.
Voucher IDSpeciesDescription
ZISP.2150Bombinator orientalisKorea, 1898, 2 sp., Imperial Russian Geographical Society (perhaps, A. Zvegintsov; description of the trip in Zvegintsov (1900); his route—Tumangan River, Purenga, Musan, Peyktusan volcano, Amnoka River, Kange, Togurion-Miti pass, Chin’-Chan’-Gan’ River, Tsynampo).
ZISP.2151Rana sp.”Korea, 1898, 1sp. (lost), Imperial Russian Geographical Society (the same).
ZISP.2835Pelophylax nigromaculatusWonsan (=Genzan), Korea, 1900, 4 sp., P.Yu. Schmidt.
ZISP.2875Pelophylax nigromaculatusWonsan (=Genzan), Korea, 1900, 5 sp., P.Yu. Schmidt.
ZISP.2876Pelophylax nigromaculatusWonsan (=Genzan), Korea, 5 July 1900, 8 sp., P.Yu. Schmidt.
ZISP.2912Bombina orientalisbetween Tonchen and Chogu-Chyen-Dogu Gulf, Korea, 1900, 4 sp., P.Yu. Schmidt
ZISP.2914Bufo stejnegeriKorea, 1900, 1 sp., P.Yu. Schmidt (description of the trip in Schmidt (1900); his route—03–24 July, Wonsan—25 July, Anbion—Suchumde Lake—27 July, Tonchen—28 July, Chogu-Chyen-Dogu Gulf—30 July, Singesa monastery, Diamond Mountains—01 August, Uonjon termal springs—02 August, Uonjon’on Pass—03 August, Chaansa monastery—04 August, Peansa monastery—05 August, Ammudo-Koge pass—05–06 August, Yuchomsa monastery—07–09 August, Kosyon Town—10–12 August, Kansen—13 August, Oridin village—14 August, Naksans monastery—15 August, Yak-Yan Mount—16 August, Kannyn Town—17 August, Kusan village—18 August, Tegulien Pass—19–20 August, Ol’chansa monastery—21 August, Kvanmul’ River—22 August, Pkhen-Khan Town—23 August, Pkhen-Khan River—24 August, Chkhe-Chkhen Mount—25 August, Iyonmusa village—26 August, Changnim village—27–28 August, Shuakori village—29 August, Chu-Yon-Che pass—30 August, Iyochen’ Town—01 August, Andon Town—02 August, Chingpo Town—04 August, Iondo Town—Khynkhe Town—Kyonju Town—12 August, Fuzan City).
ZISP.2915Hynobius leechiiKorea (according to date—Ammudo-Koge pass), 05 August 1900, 2 sp., P.Yu. Schmidt.
ZISP.2925Bombina orientalisKorea, 1900, 5 sp., P.Yu. Schmidt.
ZISP.3044Pelophylax nigromaculatusrice fields in Wonsan (Genzan), Korea, 05 July 1900, 4 sp., P.Yu. Schmidt.
ZISP.13968–13981Rana uenoiKorea, 14 sp., Chosun Academy of Science.
ZISP.13982–13990Pelophylax nigromaculatusKorea, 9 sp., Chosun Academy of Science.
ZISP.13991–13995Pelophylax chosenicusKorea, 5 sp., Chosun Academy of Science.
ZISP.13996–14000Glandirana emeljanoviKorea, 5 sp., Chosun Academy of Science.
ZISP.14001–14009Dryophytes japonicusKorea, 9 sp., Chosun Academy of Science.
ZISP.14010–14011Hynobius leechiiKorea, 2 sp., Chosun Academy of Science.
ZISP.14012Onychodactylus koreanusKorea, 1 sp., Chosun Academy of Science.
ZISP.14013–14016Bufo gargarizansKorea, 4 sp., Chosun Academy of Science.
ZISP.14017–14018Kaloula borealisKorea, 2 sp., Chosun Academy of Science.
ZISP.14019–14021DryophytessuweonensisPyongyang, Korea, 25–31 July 1947, 3 sp., V. Gnezdintsev.
ZISP.14022–14033Pelophylax nigromaculatusPyongyang, Korea, 25–31 July 1947, 12 sp., V. Gnezdintsev.
ZISP.14034–14035Rana uenoiPyongyang, Korea, 25–31 July 1947, 2 sp., V. Gnezdintsev.
ZISP.14036–14037Rana coreanaPyongyang, Korea, 25–31 July 1947, 2 sp., V. Gnezdintsev.
ZISP.14038–14044Kaloula borealisPyongyang, Korea, 25–31 July 1947, 7 sp., V. Gnezdintsev.
ZISP.14045–14047Bombina orientalisEiko, Korea, 29 August 1947, 3 sp., V. Gnezdintsev.
ZMMGU.871Pelophylax nigromaculatusKheyduzuo [=Haeju], 6 October 1947, D. I. Bibikov, 1 sp.
ZMMGU.877Rana dybowskii8 October 1947, Kheyduzuo, D. I. Bibikov, 3 sp.
ZMMGU.878Strauchbufo raddei6 October 1947, Kheyduzuo, D. I. Bibikov, 6 sp.
ZMMGU.1141Pelophylax nigromaculatusHamhung, 25 September 1970, R. Bielawski, 1 sp.
Table 2. List of samples with GenBank accession numbers used in this analysis. The species identification for individual Rana dybowskii samples from R Korea was arbitrarily reassigned to R. uenoi based on [40]. These individuals are indicated by *. References are indicated when available or left blank.
Table 2. List of samples with GenBank accession numbers used in this analysis. The species identification for individual Rana dybowskii samples from R Korea was arbitrarily reassigned to R. uenoi based on [40]. These individuals are indicated by *. References are indicated when available or left blank.
Accession NumberSpeciesLocalityReferences
Dryophytes
KY419887Dryophytes suweonensisPyeongtaek, R Korea[67]
KX854020Dryophytes suweonensisPyeongtaek, R Korea[66]
KP742630Hyla orientalisF. Russia[65]
KP742620Dryophytes japonicusTsushima, Japan
KP742616Dryophytes japonicusTaksimo, F. Russia
KP742598Dryophytes japonicusShenyang, PR China
NC010232Dryophytes japonicusHiroshima, Japan[64]
KT964710Dryophytes japonicusHeilongjiang, PR China[68]
Rana
KX024858Rana dybowskiiHuanran, PR China[42]
KX024857Rana dybowskiiFushun, PR China[42]
KX024869Rana huanrenensisHuanran, PR China[42]
KU310893Rana kukunorisPR China
DQ289092Rana chensinensisHuixian, PR China[48]
KX024861Rana uenoiTsushima, Japan[42]
KX024859Rana uenoiTsushima, Japan[42]
KX024860Rana uenoiTsushima, Japan[42]
KM876980Rana uenoi *Boeun, R Korea[42]
KM876979Rana uenoi *Boeun, R Korea[42]
AB058855Rana uenoiTsushima, Japan[42]
KX139719Rana maoershanensisGuangxi, PR China[70]
KX024866Rana piricaHokkaido, Japan[42]
KX024865Rana piricaHokkaido, Japan[42]
AB058857Rana ornativentrisAomori, Japan[69]
AB058856Rana ornativentrisHiroshima, Japan[69]
Table 3. AUC and TSS values of the 10-percentile training omission threshold generated by the MaxEnt models. The values were used to select the best fitting models.
Table 3. AUC and TSS values of the 10-percentile training omission threshold generated by the MaxEnt models. The values were used to select the best fitting models.
SpeciesAUC ± SDTSS ± SD
Bombina orientalis0.8651 ± 0.01160.4626 ± 0.0805
Bufo gargarizans0.8617 ± 0.01990.4570 ± 0.0568
Bufo stejnegeri0.9387 ± 0.00780.6905 ± 0.1167
Dryophytes japonicus0.8459 ± 0.01550.4958 ± 0.0525
Dryophytes suweonensis0.9637 ± 0.00640.8148 ± 0.0767
Glandirana emeljanovi0.8877 ± 0.03670.5080 ± 0.1041
Hynobius leechii0.8459 ± 0.01340.4538 ± 0.0544
Kaloula borealis0.9158 ± 0.01420.6230 ± 0.1133
Onychodactylus sp.0.9185 ± 0.01470.5924 ± 0.1078
Onychodactylus koreanus0.9568 ± 0.00950.6962 ± 0.1314
Pelophylax chosenicus0.9451 ± 0.00540.7137 ± 0.0810
Pelophylax nigromaculatus0.8324 ± 0.01050.4085 ± 0.0518
Strauchbufo raddei0.9881 ± 0.00520.9211 ± 0.0382
Rana amurensis0.9104 ± 0.00910.6929 ± 0.0161
Rana coreana0.8697 ± 0.01190.4652 ± 0.1158
Rana dybowskii0.9268 ± 0.02800.5747 ± 0.1993
Rana huanrenensis0.9913 ± 0.00300.8411 ± 0.1978
Rana uenoi0.8869 ± 0.00790.5874 ± 0.0315
Salamandrella tridactyla0.9716 ± 0.01430.4488 ± 0.0909
Table 4. Summary of population dynamics of amphibian populations in DPR Korea or closely related regions, extracted from the literature. W * indicates “workshop in Pyongyang in June 2019”. Here, [89] is cited as reported in [51].
Table 4. Summary of population dynamics of amphibian populations in DPR Korea or closely related regions, extracted from the literature. W * indicates “workshop in Pyongyang in June 2019”. Here, [89] is cited as reported in [51].
SpeciesPopulation DynamicsReference
Bufo gargarizansCommon, reported from urban districts[51,89]
Strauchbufo raddeiUncommon or only locally abundant[51,89]
Dryophytes japonicusPresent all over the nation[51,89]
Bombina orientalisPresent all over the nation[51,89]
Glandirana emeljanoviOne of the most common species in Central Korea[50]
Hynobius leechiiCommon in mountain forests[51]
OnychodactylusCommon in mountain forests[51]
Salamandrella tridactylaOnly in the high regions of the mountains[51]
Pelophylax nigromaculatusAbundant in rice paddies[113]
Rana dybowskii-R. uenoiLocally abundant[113]
Bombina orientalisAbundant in the lowlands[113]
Bufo gargarizansAbundant in the lowlands[113]
Dryophytes japonicusAbundant in the lowlands[113]
Pelophylax chosenicusDecline in abundance. Decline in North Pyongan province and Seoncheon because of droughts likely linked to climate changeW *
Rana coreanaDecline in abundanceW *
Bufo gargarizans:High density on the west coastW *
Glandirana emeljanoviNot abundant but widespreadW *
Kaloula borealisContinuous presence but not common, all the way to North Hamgyong provinceW *
Pelophylax nigromaculatusMost abundant amphibian speciesW *
Strauchbufo raddeiIsolated populations in Shinuiju, Amrok River, Chilbeok mountain, North Pyongan provinceW *
Rana dybowskiiAbundant at high altitude, populations >2100 m at Baekdu lakeW *
Hynobius leechiiAbundant <1000 m asl throughout its rangeW *
Lithobates catesbeianusIsolated large populations originating from escaped individuals from farmsW *
Rana sp.Isolated large populations originating from escaped individuals from farmsW *
Table 5. Description of presence models used for each species in DPR Korea. The surface area in km2 provided matches with the 10-percentile training and the 95% presence used in models. IUCN threat levels as of February 2019 (www.iucnredlist.org (accessed on 1 February 2019)).
Table 5. Description of presence models used for each species in DPR Korea. The surface area in km2 provided matches with the 10-percentile training and the 95% presence used in models. IUCN threat levels as of February 2019 (www.iucnredlist.org (accessed on 1 February 2019)).
SpeciesMinimum Suitable Area (km2)10-Percentile Omission Area (km2)IUCN StatusRecommendation at National LevelLocal Name
Bombina orientalis119,53130,729LCLC비단개구리
Bufo gargarizans101,62729,954LCLC두꺼비
Bufo stejnegeri45,30111,805LCVU B1ab(i,ii,iii)금강두꺼비
Dryophytes japonicus114,92451,290LCLC청개구리
Dryophytes suweonensis25,2149785ENEN B1ab(i,ii,iii)수원청개구리
Glandirana emeljanovi107,08530,201LCLC옴개구리
Hynobius leechii60,54141,443LCLC도롱룡
Kaloula borealis40,57127,098LCNT맹꽁이
Onychodactylus sp.57,96020,119
Onychodactylus koreanus86376DDVU B1ab(i,ii,iii)발톱도롱룡
Pelophylax chosenicus35,90813,810VUVU B1ab(i,ii,iii)금개구리
Pelophylax nigromaculatus102,32851,779NTNT참개구리
Strauchbufo raddei22,0267,775LCVU B1ab(i,ii,iii)작은두꺼비
Rana amurensis110,68717,271LCLC양용개구리
Rana coreana54,74026,249LCNT애기개구리
Rana dybowskii101,43030,360LCLC북방산개구리/기름개구리
Rana huanrenensis78,53334,021LCLC계곡산개구리
Rana uenoi36,0669579NEVU B1ab(i,ii,iii)
Salamandrella tridactyla15,9866294NENT합수도롱룡
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Borzée, A.; Litvinchuk, S.N.; Ri, K.; Andersen, D.; Nam, T.Y.; Jon, G.H.; Man, H.S.; Choe, J.S.; Kwon, S.; Othman, S.N.; et al. Update on Distribution and Conservation Status of Amphibians in the Democratic People’s Republic of Korea: Conclusions Based on Field Surveys, Environmental Modelling, Molecular Analyses and Call Properties. Animals 2021, 11, 2057. https://doi.org/10.3390/ani11072057

AMA Style

Borzée A, Litvinchuk SN, Ri K, Andersen D, Nam TY, Jon GH, Man HS, Choe JS, Kwon S, Othman SN, et al. Update on Distribution and Conservation Status of Amphibians in the Democratic People’s Republic of Korea: Conclusions Based on Field Surveys, Environmental Modelling, Molecular Analyses and Call Properties. Animals. 2021; 11(7):2057. https://doi.org/10.3390/ani11072057

Chicago/Turabian Style

Borzée, Amaël, Spartak N. Litvinchuk, Kyongsim Ri, Desiree Andersen, Tu Yong Nam, Gwang Hyok Jon, Ho Song Man, Jong Sik Choe, Sera Kwon, Siti N. Othman, and et al. 2021. "Update on Distribution and Conservation Status of Amphibians in the Democratic People’s Republic of Korea: Conclusions Based on Field Surveys, Environmental Modelling, Molecular Analyses and Call Properties" Animals 11, no. 7: 2057. https://doi.org/10.3390/ani11072057

APA Style

Borzée, A., Litvinchuk, S. N., Ri, K., Andersen, D., Nam, T. Y., Jon, G. H., Man, H. S., Choe, J. S., Kwon, S., Othman, S. N., Messenger, K., Bae, Y., Shin, Y., Kim, A., Maslova, I., Luedtke, J., Hobin, L., Moores, N., Seliger, B., ... Jang, Y. (2021). Update on Distribution and Conservation Status of Amphibians in the Democratic People’s Republic of Korea: Conclusions Based on Field Surveys, Environmental Modelling, Molecular Analyses and Call Properties. Animals, 11(7), 2057. https://doi.org/10.3390/ani11072057

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop