Research Article |
Corresponding author: Minh Duc Le ( minh.le.cres@gmail.com ) Corresponding author: Thomas Ziegler ( ziegler@koelnerzoo.de ) Academic editor: Franco Andreone
© 2024 Minh Duc Le, Dennis Rödder, Tao Thien Nguyen, Cuong The Pham, Truong Quang Nguyen, An Vinh Ong, Timothy E. M. McCormack, Thang Tai Nguyen, Mai Huyen Le, Hanh Thi Ngo, Thomas Ziegler.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Le MD, Rödder D, Nguyen TT, The Pham C, Nguyen TQ, Ong AV, McCormack TEM, Nguyen TT, Le MH, Ngo HT, Ziegler T (2024) Climatic niche modelling and genetic analyses highlight conservation priorities for the Spotted Softshell Turtle (Pelodiscus variegatus). Nature Conservation 55: 67-82. https://doi.org/10.3897/natureconservation.55.114746
|
The Spotted Softshell Turtle (Pelodiscus variegatus) has been recognised since 2019 from Vietnam and Hainan Island, China, but little information about its population status and distribution range is currently available. The species has been provisionally listed as Critically Endangered by the Turtle and Tortoise Working Group, although the status has not been officially accepted by the IUCN, due to the threats the species is facing, including habitat loss and degradation, overexploitation for food, competition with other non-native softshell turtles and pollution. To identify conservation priority sites for P. variegatus in mainland Indochina, this study combines molecular analyses and species distribution modelling. Our results show that, in Vietnam, Phong Nha-Ke Bang National Park has the largest suitable area and high probability of species occurrence, followed by Vu Quang National Park and Song Thanh and Ke Go Nature Reserves. In addition, the central provinces, from Thanh Hoa to Thua Thien Hue in Vietnam, constitute a key part of the species distribution and should be prioritised for conservation actions. According to the study’s findings, although P. variegatus is possibly found in Laos, the probability decreases sharply at the border between both countries and there is also a gap in the occurrence of wetlands, arguing for strong natural barriers. Unfortunately, to date, only part of the species potential distribution is protected, while no records are known from protected areas, highlighting the need for extended or even new reserves. To recover natural populations of the species and following the IUCN’s One Plan Approach to Conservation, breeding programmes have been established in Vietnam with a potential to expand to other facilities in the country and abroad. Once suitable sites are identified, offspring can be released into the protected areas to improve the current conservation status of this highly-threatened softshell turtle.
Cytochrome b, Laos, ND4, species distribution modelling, Vietnam
The Chinese Softshell Turtle, Pelodiscus sinensis, was described nearly 200 years ago and was believed to represent a morphologically variable, geographically widespread taxon (from the Russian Far East through the Korean Peninsula, eastern and central China to Vietnam). The Northern Chinese Softshell Turtle, P. maackii (Brandt, 1857), was described 23 years later, but was then thought to be a synonym of P. sinensis. Only 35 years ago, the populations from the northernmost part of the P. sinensis distribution range were shown to represent a distinct species, based on osteological differences (
Two additional species of this complex from central China were described in the 90s, based on morphological differences: the Hunan Softshell Turtle (P. axenaria) and the Lesser Chinese Softshell Turtle (P. parviformis) (
Thus, the Pelodiscus sinensis complex at the moment comprises seven species, with six species being distributed in China and four endemic to the country (
The exact range of the recently-described P. variegatus is still largely unknown, but historical records suggest that the species occupies lowland areas of central and northern Vietnam and parts of southern China, viz. Hainan Province (
As softshell turtles are common and prized as food where they occur, natural populations are threatened by local hunting, with further threats of habitat loss and competition with introduced softshell turtles (
To identify new samples collected in the field, we used a molecular approach. In total, 61 newly-collected samples were included in the analyses (Suppl. material
Primer | Sequence | Reference |
---|---|---|
Gludg (f) | 5′- TGACTTGAARAACCAYCGTTG - 3′ |
|
CB3 (r) | 5′- GGCAAATAGGAAATATCATTC - 3′ |
|
CB534 (f) | 5′- GACAATGCAACCCTAACACG- 3′ |
|
Tcytbthr (r) | 5′- TTCTTTGGTTTACAAGACC - 3′ |
|
ND4 672 (f) | 5′- TGACTACCAAAAGCTCATGTAGAAGC - 3′ |
|
Hist (r) | 5′- CCTATTTTTAGAGCCACAGTCTAATG - 3′ |
|
Data were then analysed using Maximum Likelihood (ML) as implemented in IQ-TREE v.1.6.7.1 (
This study only employed mitochondrial genes to provide taxonomic identification of samples collected from the wild. Although maternally inherited mitochondrial loci cannot help detect hybridisation events, interbreeding between different species of softshell turtles has only been reported in turtle farms (
For species records, we carefully checked potential records listed in recent studies, including,
For environmental predictors, we used a combination of weather station-derived precipitation data and remote sensing data. The full list of variables including their interpretation is provided in Table
Variables used for climatic niche modeling computation. NDVI = Normalized Difference Vegetation Index, EVI = Enhanced Vegetation Index.
Abbreviation | Remote sensing variable | Bioclimatic variable | Derived variable | Source |
---|---|---|---|---|
bio_12 | N/A | Annual Precipitation | N/A | Worldclim 2.1 |
bio_13 | N/A | Precipitation of Wettest Month | N/A | Worldclim 2.1 |
bio_14 | N/A | Precipitation of Driest Month | N/A | Worldclim 2.1 |
bio_15 | N/A | Precipitation Seasonality (Coefficient of Variation) | N/A | Worldclim 2.1 |
bio_16 | N/A | Precipitation of Wettest Quarter | N/A | Worldclim 2.1 |
bio_17 | N/A | Precipitation of Driest Quarter | N/A | Worldclim 2.1 |
bio_18 | N/A | Precipitation of Warmest Quarter | N/A | Worldclim 2.1 |
bio_19 | N/A | Precipitation of Coldest Quarter | N/A | Worldclim 2.1 |
ED1514_bio1 | MODIS V4 Band 14 Synoptic Months: NDVI | BIO1 = Annual Mean Temperature | Annual Mean of NDVI | EDENext |
ED1514_bio2 | MODIS V4 Band 14 Synoptic Months: NDVI | BIO2 = Mean Diurnal Range (Mean of monthly (max temp - min temp)) | Mean Diurnal Range of NDVI | EDENext |
ED1514_bio3 | MODIS V4 Band 14 Synoptic Months: NDVI | BIO3 = Isothermality (BIO2/BIO7) (×100) | Isothermaility (BIO2/BIO7) (*100) of NDVI | EDENext |
ED1514_bio4 | MODIS V4 Band 14 Synoptic Months: NDVI | BIO4 = Temperature Seasonality (standard deviation ×100) | Seasonality of NDVI | EDENext |
ED1514_bio5 | MODIS V4 Band 14 Synoptic Months: NDVI | BIO5 = Max Temperature of Warmest Month | Max NDVI of Monthly Scores | EDENext |
ED1514_bio6 | MODIS V4 Band 14 Synoptic Months: NDVI | BIO6 = Min Temperature of Coldest Month | Min NDVI of Monthly Scores | EDENext |
ED1514_bio7 | MODIS V4 Band 14 Synoptic Months: NDVI | BIO7 = Temperature Annual Range (BIO5-BIO6) | Annual Range of NDVI | EDENext |
ED1514_bio10 | MODIS V4 Band 14 Synoptic Months: NDVI | BIO10 = Mean Temperature of Warmest Quarter | Mean NDVI of Warmest Quarter | EDENext |
ED1514_bio11 | MODIS V4 Band 14 Synoptic Months: NDVI | BIO11 = Mean Temperature of Coldest Quarter | Mean NDVI of Coldest Quarter | EDENext |
ED1515_bio1 | MODIS V4 Band 15 Synoptic Months: EVI | BIO1 = Annual Mean Temperature | Annual Mean of EVI | EDENext |
ED1515_bio2 | MODIS V4 Band 15 Synoptic Months: EVI | BIO2 = Mean Diurnal Range (Mean of monthly (max temp - min temp)) | Mean Diurnal Range of EVI | EDENext |
ED1515_bio3 | MODIS V4 Band 15 Synoptic Months: EVI | BIO3 = Isothermality (BIO2/BIO7) (×100) | Isothermaility (BIO2/BIO7) (*100) of EVI | EDENext |
ED1515_bio4 | MODIS V4 Band 15 Synoptic Months: EVI | BIO4 = Temperature Seasonality (standard deviation ×100) | Seasonality of EVI | EDENext |
ED1515_bio5 | MODIS V4 Band 15 Synoptic Months: EVI | BIO5 = Max Temperature of Warmest Month | Max EVI of Monthly Scores | EDENext |
ED1515_bio6 | MODIS V4 Band 15 Synoptic Months: EVI | BIO6 = Min Temperature of Coldest Month | Min EVI of Monthly Scores | EDENext |
ED1515_bio7 | MODIS V4 Band 15 Synoptic Months: EVI | BIO7 = Temperature Annual Range (BIO5-BIO6) | Annual Range of EVI | EDENext |
ED1515_bio10 | MODIS V4 Band 15 Synoptic Months: EVI | BIO10 = Mean Temperature of Warmest Quarter | Mean EVI of Warmest Quarter | EDENext |
ED1515_bio11 | MODIS V4 Band 15 Synoptic Months: EVI | BIO11 = Mean Temperature of Coldest Quarter | Mean EVI of Coldest Quarter | EDENext |
ED15078_bio1 | MODIS V4 Band 07+08 Synoptic Months: day- & nighttime land surface temperature | BIO1 = Annual Mean Temperature | Annual Mean Temperature | EDENext |
ED15078_bio2 | MODIS V4 Band 07+08 Synoptic Months: day- & nighttime land surface temperature | BIO2 = Mean Diurnal Range (Mean of monthly (max temp - min temp)) | Mean Diurnal Range of Temperature | EDENext |
ED15078_bio3 | MODIS V4 Band 07+08 Synoptic Months: day- & nighttime land surface temperature | BIO3 = Isothermality (BIO2/BIO7) (×100) | Isothermality (Bio2/Bio7) (*100) | EDENext |
ED15078_bio4 | MODIS V4 Band 07+08 Synoptic Months: day- & nighttime land surface temperature | BIO4 = Temperature Seasonality (standard deviation ×100) | Seasonality | EDENext |
ED15078_bio5 | MODIS V4 Band 07+08 Synoptic Months: day- & nighttime land surface temperature | BIO5 = Max Temperature of Warmest Month | Max Temperature of Warmest Month | EDENext |
ED15078_bio6 | MODIS V4 Band 07+08 Synoptic Months: day- & nighttime land surface temperature | BIO6 = Min Temperature of Coldest Month | Min Temperature of Coldest Month | EDENext |
ED15078_bio7 | MODIS V4 Band 07+08 Synoptic Months: day- & nighttime land surface temperature | BIO7 = Temperature Annual Range (BIO5-BIO6) | Temperature Annual Range | EDENext |
ED15078_bio10 | MODIS V4 Band 07+08 Synoptic Months: day- & nighttime land surface temperature | BIO10 = Mean Temperature of Warmest Quarter | Mean Temperature of Warmest Quarter | EDENext |
ED15078_bio11 | MODIS V4 Band 07+08 Synoptic Months: day- & nighttime land surface temperature | BIO11 = Mean Temperature of Coldest Quarter | Mean Temperature of Coldest Quarter | EDENext |
Multi-collinearity amongst predictors may hamper successful model training and subsequent projection (
As the algorithm for climatic niche modelling development, we used Maxent v.3.4.0, which is specifically designed to derive potential distributions from presence-pseudoabsence data (
For the second land-cover niche modelling, we reclassified the potential distribution suggested by the first model applying the minimum training presence threshold, which was used as environmental background. Model selection followed again
We merged the occurrence data with existing protected areas (Reserves, National Parks etc.) in the country. Information of protected areas was obtained from the world dictionary of protected areas/protected planet (https://www.protectedplanet.net). We selected the targeted area in north-central and central Vietnam because most of distribution records were reported from the region. In total, there were 42 potential conservation units within the general area and, for each Reserve, we computed the number of suitable grid cells, the sum of probabilities and the mean probability in QGIS 3.14. Map resolution ca. 1 km (30 arc sec).
The molecular matrix contained 1921 aligned characters. Both BI and ML analyses showed that new samples belong to Pelodiscus sinensis and P. variegatus. The former species was only moderately supported (PP = 89%, UBF = 90%), while the latter received strong support from both BI and ML (PP = 100%, UBF = 99%). In total, 61 newly-collected and four GenBank samples were identified as P. variegatus (see Suppl. material
The niche modelling trained with only climatic variables had a good overall performance (AUCtest = 0.79; AUCtraining = 0.84, Suppl. material
The niche modelling, trained with microhabitat variables, had a good overall performance (AUCtest = 0.73; AUCtraining = 0.80, Suppl. material
When integrating the probabilities of occurrence derived from climatic and land-cover variables, the most suitable habitats for P. variegatus are near the Vietnam coastline, where extensive freshwater wetlands exist (Fig.
Only part of the potential distribution for P. variegatus in Vietnam is protected and no known occurrence is directly located within reserves, although they are close, such as nearby Ke Go Nature Reserve and Vu Quang National Park (Fig.
Protected areas in Vietnam with projected proper climatic conditions and land-cover (Average Probability > 0.1) for Pelodiscus variegatus, sorted by size of suitable area. A full list of protected areas in China, Laos and Vietnam is presented in the Suppl. material
Name | Type | IUCN category | Area [km2] | Average Probability | Suitable Area [km2] |
---|---|---|---|---|---|
Phong Nha-Ke Bang | National Park | II | 1222.825 | 0.303 | 424.578 |
Vu Quang | National Park | II | 591.661 | 0.275 | 168.536 |
Song Thanh | Nature Reserve | Not Reported | 890.589 | 0.110 | 117.079 |
Ke Go | Nature Reserve | IV | 239.724 | 0.364 | 91.145 |
Ben En | National Park | II | 141.816 | 0.123 | 85.311 |
Phong Dien | Nature Reserve | IV | 407.762 | 0.417 | 85.060 |
Bach Ma | National Park | II | 375.506 | 0.389 | 79.732 |
Nui Coc | Cultural and Historical Site | V | 90.263 | 0.200 | 60.395 |
Hoa Lu | Cultural and Historical Site | Not Reported | 65.819 | 0.126 | 59.431 |
Bac Huong Hoa | Nature Reserve | Not Reported | 235.057 | 0.162 | 56.869 |
Deo Ca-Hon Nua | Cultural and Historical Site | V | 215.611 | 0.122 | 56.472 |
Than Sa-Phuong Hoang | Nature Reserve | Not Reported | 136.517 | 0.112 | 35.724 |
Nui Thanh | Cultural and Historical Site | V | 64.139 | 0.134 | 34.971 |
Dakrong | Nature Reserve | Not Reported | 387.232 | 0.351 | 32.474 |
Dao Ho Song Da | Cultural and Historical Site | V | 78.889 | 0.130 | 31.908 |
Bac Me | Nature Reserve | Not Reported | 87.109 | 0.122 | 24.266 |
Hue Sao La | Nature Reserve | Not Reported | 377.864 | 0.422 | 17.208 |
Hon Ba | Nature Reserve | Not Reported | 198.213 | 0.146 | 14.560 |
Ngoc Linh (Quang Nam) | Nature Reserve | Not Reported | 192.496 | 0.146 | 10.710 |
Quy Hoa-Ghenh Rang | Cultural and Historical Site | Not Reported | 52.102 | 0.115 | 7.707 |
Ba Na-Nui Chua | Nature Reserve | IV | 268.922 | 0.226 | 6.019 |
Nui Chung | Cultural and Historical Site | Not Reported | 2.031 | 0.295 | 2.020 |
Bai Tu Long | National Park | Not Reported | 64.906 | 0.214 | 1.401 |
Son Tra | Nature Reserve | Not Reported | 38.368 | 0.198 | 1.213 |
Trimmed phylogram, based on the Bayesian analysis. Number above and below branches of major nodes are Bayesian posterior probabilities and ML ultrafast bootstrap values, respectively. Sample highlighted in bold and orange is the paratype of Pelodiscus variegatus.
Our habitat suitability analysis predicts that two most important protected areas for P. variegatus include Phong Nha-Ke Bang and Vu Quang National Parks in Vietnam. Other protected sites with the largest suitable sizes consist of Song Thanh and Ke Go Nature Reserves and Ben En National Park (Table
Potential distribution of Pelodiscus variegatus in Vietnam based on A climate B land cover C climate and land cover and D coverage with protected areas as number of suitable grid cells.
Prior to its discovery, P. variegatus was considered part of P. parviformis. The latter species was already assessed as threatened and included on appendix II of the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES). As its southern population has been split into a different species, P. variegatus, the overall population of each species becomes even smaller than previously thought.
As a first measure, based on the genetically-identified individuals, a conservation breeding programme has been established. This is following IUCN’s One Plan Approach to Conservation, developed by the Conservation Planning Specialist Group (CPSG), which combines in situ with ex situ conservation measures for the optimum protection of a given species (
To maximise positive outcomes and for security reasons, viz., to extend the conservation breeding network, another group of the genetically-identified P. variegatus was provided to another softshell breeding facility in northern Vietnam. Successful breeding has already been observed in the colony and offspring are ready for release to the original habitat sites. To extend the conservation breeding programme and, thus, contribute to the build-up of a stable assurance colony and conservation breeding network, a plan has been developed to transfer a limited number of surplus offspring to other facilities in Vietnam and overseas. In addition to these ex situ conservation measures already being in place, focus should now be directed to improving in situ conservation of this beautiful, but threatened softshell turtle species.
In the late 2023, 50 young and healthy spotted softshell turtles from the in-country breeding program initiated by the Institute of Ecology and Biological Resources (IEBR), Vietnam, together with the Cologne Zoo, Germany, were successfully released to a site in northern Vietnam.
We thank Jonathan Fong, Luca Luiselli and two anonymous reviewers for their insightful comments on the earlier version of the paper and John B. Iverson and Anders G. J. Rhodin for providing valuable distribution data.
The authors have declared that no competing interests exist.
No ethical statement was reported.
This study was partially supported by the Ministry of Education and Training of Vietnam (Project Code B2020-TDV-07) and the US Fish and Wildlife Service. Previous genetic screening was funded by Cologne Zoo and the European Union of Aquarium Curators (EUAC). Cologne Zoo is partner of the World Association of Zoos and Aquariums (WAZA): Conservation Projects 07011, 07012 (Herpetodiversity Research, Amphibian and Reptilian Breeding and Rescue Stations).
TZ, MDL, TTN, and TQN conceptualized the study; CTP, AVO, TEMM, and TTN conducted the fieldwork; DR, HTN, MDL, MHL, and TTN led the data analysis, data curation; TZ, MDL, DR led the writing and all authors edited and approved the manuscript.
Minh Duc Le https://orcid.org/0000-0002-2953-2815
Dennis Rödder https://orcid.org/0000-0002-6108-1639
Tao Thien Nguyen https://orcid.org/0000-0002-5640-4536
Cuong The Pham https://orcid.org/0000-0001-5158-4526
Truong Quang Nguyen https://orcid.org/0000-0002-6601-0880
An Vinh Ong https://orcid.org/0000-0003-3683-3832
Timothy E. M. McCormack https://orcid.org/0009-0003-7565-2345
Thang Tai Nguyen https://orcid.org/0009-0002-9088-3515
Mai Huyen Le https://orcid.org/0009-0002-1536-1274
Hanh Thi Ngo https://orcid.org/0000-0002-5283-6243
Thomas Ziegler https://orcid.org/0000-0002-4797-609X
All of the data that support the findings of this study are available in the main text or Supplementary Information.
Supplementary data
Data type: docx
Explanation note: table S1. Samples used in this study. table S2. Relative importance of protected areas in China, Laos and Vietnam for Pelodiscus variegatus in terms of climatic conditions and land-cover. The Table is sorted according to the suitable area in each Reserve. figure S1. Full phylogram based on the Bayesian analysis. Numbers above and below branches of selected nodes are Bayesian posterior probabilities and ML ultrafast bootstrap values, respectively. figure S2. Summary of the receiver operating characteristic curve of 100 Maxent models for P. variegatus trained with climatic variables. figure S3. Summary of the receiver operating characteristic curve of 100 Maxent models for P. variegatus trained with land-cover variables.