Population genomics identifies genetic structure and admixture in the endangered Beale’s Eyed Turtle (Sacalia bealei), and implications for aquatic ecology and ex situ breeding

Abstract

Freshwater turtles are one of the most threatened animal groups in the world, especially in Asia. The Beale’s Eyed Turtle (Sacalia bealei) is a highly endangered Chinese species at risk of extinction due to overexploitation for the food and pet trades. Nonetheless, Hong Kong still houses a handful of relatively healthy populations of this species, providing an increasingly rare opportunity to study the population genetics of wild populations and, in turn, guide conservation action. In this study, we used genome-wide double-digest restriction-site associated DNA sequencing (ddRADseq) markers to examine the population structure of S. bealei from all known wild individuals with samples suitable for DNA sequencing, representing a subset of its distribution in Hong Kong (two localities: HK1, HK2) and Fujian Province. Although genetic diversity is relatively low, we recovered four genetic clusters, with three corresponding to each of the three known localities (HK1, HK2, Fujian), and admixture between the clusters. Captive-breeding colony individuals showed mixed geographic origins, with approximately half clustering with Fujian populations and half with Hong Kong populations. These data were used to evaluate the genetic diversity and infer the geographic origin of an ex situ breeding colony with individuals of unknown provenance. These results provide important baseline information on the population structure of wild S. bealei and the potential geographic origins of captive turtles, which directly contribute to in situ and ex situ conservation efforts of the species.

Keywords

Conservation genomics / ddRADseq / Ex situ conservation / Population structure / Sacalia bealei

References

[1] Anderson EC, Thompson EA. A model-based method for identifying species hybrids using multilocus genetic data. Genet. 2002, 160: 1217-1229.

[2] Asian Turtle Trade Working Group. Sacalia bealei (errata version published in 2016). IUCN Red List Threat Species. 2000, 2000: eT19796A97377507

[3] Baym M, Kryazhimskiy S, Lieberman TD, Chung H, Desai MM, Kishony R. Inexpensive multiplexed library preparation for megabase-sized genomes. PLoS ONE. 2015, 10e0128036.

[4] Catchen J, Hohenlohe P, Bassham S, Amores A, Cresko W. Stacks: an analysis tool set for population genomics. Mol Ecol. 2013, 223124-3140.

[5] Cheung SM, Dudgeon D. Quantifying the Asian turtle crisis: market surveys in southern China, 2000–2003. Aquat Conserv Mar Freshw Ecosyst. 2006, 16: 751-770.

[6] Çilingir FG, Rheindt FE, Garg KM, Platt K, Platt SG, Bickford DP. Conservation genomics of the endangered Burmese roofed turtle. Conserv Biol. 2017, 31: 1469-1476.

[7] Çilingir FG, Seah A, Horne BD, Som S, Bickford DP, Rheindt FE. Last exit before the brink: conservation genomics of the Cambodian population of the critically endangered southern river terrapin. Ecol Evol. 2019, 9: 9500-9510.

[8] Converse PE, Kuchta SR, Hauswaldt JS, Roosenburg WM. Turtle soup, prohibition, and the population genetic structure of diamondback terrapins (Malaclemys terrapin). PLoS ONE. 2017, 12: e0181898.

[9] Do C, Waples RS, Peel D, Macbeth GM, Tillett BJ, Ovenden JR. NeEstimator v2: re-implementation of software for the estimation of contemporary effective population size (N_e) from genetic data. Mol Ecol Resour. 2014, 14209-214.

[10] Dresser CM, Pierson TW, Fitzpatrick BM. Isolation by distance, local adaptation, and fortuitous coincidence of geo-political boundaries with spatial-genetic clusters in southern Bog Turtles. Glob Ecol Conserv. 2018, 16: e00474

[11] Earl DA, von Holdt BM. Structure harvester: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012, 4359-361.

[12] Endler J, Mazer S, Williams M, Sandoval C, Ferren W. Problems associated with the introduction of non-native genotypes on NRS reserves. 2010, California, University of California Natural Reserve System17

[13] Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software structure: a simulation study. Mol Ecol. 2005, 14: 2611-2620.

[14] Fong JJ, Qiao G. New localities of endangered Chinese turtles from museum specimens and the practical and ethical challenges using and reporting natural history collection data. Zootaxa. 2010, 2395: 59-68

[15] Frankham R. Genetics and extinction. Biol Conserv. 2005, 126: 131-140.

[16] Frankham R, Ballou J, Eldridge MDB, Lacy RC, Ralls K, Dudash MR, Fenster CB. Predicting the probability of outbreeding depression. Conserv Biol. 2011, 25: 465-475.

[17] Gong SP, Chow AT, Fong JJ, Shi HT. The chelonian trade in the largest pet market in China: scale, scope and impact on turtle conservation. Oryx. 2009, 43: 213-216.

[18] González-García MDP, García-Sais JR, García-Moliner G, Schizas NV. RAD-Seq-derived SNPs reveal no local population structure in the commercially important deep-sea queen snapper (Etelis oculatus) in Puerto Rico. Mar Life Sci. 2025, 7: 594-605.

[19] He B, Liu Y, Shi H, Zhang J, Hu M, Ma Y, Fu L, Hong M, Wang J, Fong JJ, Parham JF. Captive breeding of the four-eyed turtle (Sacalia quadriocellata). Asian Herpetol Res. 2010, 1: 111-117

[20] Hong X, Zhang X, Liu X, Wang Y, Yu L, Li W, Zhu X. Status and analysis of artificial breeding and management of aquatic turtles in China. Biol. 2022, 11: 1368.

[21] IUCN/SSC (2013) Guidelines for reintroductions and other conservation translocations. Version 1.0. Gland, Switzerland: IUCN Species Survival Commission.

[22] Janes JK, Miller JM, Dupuis JR, Malenfant RM, Gorrell JC, Cullingham CI, Andrew RL. The K= 2 conundrum. Mol Ecol. 2017, 263594-3602.

[23] Kopelman NM, Mayzel J, Jakobsson M, Rosenberg NA, Mayrose I. Clumpak: a program for identifying clustering modes and packaging population structure inferences across K. Mol Ecol Res. 2015, 15: 1179-1191.

[24] Laikre L, Schwartz MK, Waples RS, Ryman N. Compromising genetic diversity in the wild: unmonitored large-scale release of plants and animals. Trends Ecol Evol. 2010, 25: 520-529.

[25] Lau M, Shi HT. Conservation and trade of terrestrial and freshwater turtles and tortoises in the People’s Republic of China. Chelonian Res Monogr. 2000, 2: 30-38

[26] Lawson DJ, Hellenthal G, Myers S, Falush D. Inference of population structure using dense haplotype data. PLoS Genet. 2012, 8: e1002453.

[27] Le MD, McCormack TEM, Hoang HV, Duong HT, Nguyen TQ, Ziegler T, Nguyen HD, Ngo HT. Threats from wildlife trade: the importance of genetic data in safeguarding the endangered four-eyed turtle (Sacalia quadriocellata). Nat Conserv. 2020, 4191-111.

[28] Lin L, Hu Q, Fong JJ, Yang J, Chen Z, Zhou F, Wang J, Xiao F, Shi H. Reproductive ecology of the endangered Beal’s-eyed turtle, Sacalia bealei. PeerJ. 2018, 6: e4997.

[29] Lin L, Chen HQ, Gaillard D, Shi HT, Luo SJ. Phylogeography of the endangered “eyed” turtles (genus Sacalia) and the discovery of a lineage derived from natural interspecific hybridization. Ecol Evol. 2022, 12e9545.

[30] Malinsky M, Trucchi E, Lawson DJ, Falush D. RADpainter and fineRADstructure: population inference from RADseq data. Mol Biol Evol. 2018, 351284-1290.

[31] Meirmans PG. Genodive version 3.0: easy-to-use software for the analysis of genetic data of diploids and polyploids. Mol Ecol Res. 2020, 20: 126-1131.

[32] Nash HC, Wirdateti LGW, Choo SW, Chong JL, Semiadi G, Hari R, Sulaiman MH, Turvey ST, Evans TA, Rheindt FE. Conservation genomics reveals possible illegal trade routes and admixture across pangolin lineages in Southeast Asia. Conserv Genet. 2018, 19: 1083-1095.

[33] Nelson HV, Farquharson KA, Georges A, McLennan EA, DeGabriel JL, Giese M, Ormond C, McFadden M, Skidmore A, Prangell J, Belov K, Hogg CJ. A genomic framework to assist conservation breeding and translocation success: a case study of a critically endangered turtle. Conserv Sci Pract. 2024, 6: e13204.

[34] Peterson BK, Weber JN, Kay EH, Fisher HS, Hoekstra HE. Double digest RADseq: an inexpensive method for de novo SNP discovery and genotyping in model and non-model species. PLoS ONE. 2012, 7e37135.

[35] Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genet. 2000, 155945-959.

[36] Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira M, Bender D, Maller J, Sklar P, de Bakker P, Daly M, Sham P. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81: 559-575.

[37] Quinzin MC, Sandoval-Castillo J, Miller JM, Beheregaray LB, Russello MA, Hunter EA, Gibbs JP, Tapia W, Villalva F, Caccone A. Genetically informed captive breeding of hybrids of an extinct species of Galapagos tortoise. Conserv Biol. 2019, 33: 1404-1414.

[38] Rivera-Colón AG, Catchen J. Population genomics analysis with RAD, reprised: Stacks 2. Marine genomics: methods and protocols. 2021, New York, Springer99149

[39] Schmidt DJ, Espinoza T, Connell M, Hughes JM. Conservation genetics of the Mary River turtle (Elusor macrurus) in natural and captive populations. Aquat Conserv Mar Freshw Ecosyst. 2018, 28: 115-123.

[40] Shelmidine N, Murphy B, Massarone K. Husbandry and propagation of the Chinese big-headed turtle (Platysternon megacephalum) at the Wildlife Conservation Society’s Prospect Park Zoo. Zoo Biol. 2016, 35: 174-179.

[41] Shi H, Fong JJ, Parham JF, Pang J, Wang J, Hong M, Zhang YP. Mitochondrial variation of the “eyed” turtles (Sacalia) based on known-locality and trade specimens. Mol Phylogenet Evol. 2008, 49: 1025-1029.

[42] Sung YH, Fong JJ. Assessing consumer trends and illegal activity by monitoring the online wildlife trade. Biol Conserv. 2018, 227: 219-225.

[43] Sung YH, Karraker NE, Hau BCH. Demographic evidence of illegal harvesting of an endangered Asian turtle. Conserv Biol. 2013, 27: 1421-1428.

[44] Sung YH, Liew JH, Chan HK, Lee WH, Wong BHF, Dingle C, Karraker NE, Spencer RJ, Fong JJ. Assessing the diet of the endangered Beale’s eyed turtle (Sacalia bealei) using faecal content and stable isotope analyses: implications for conservation. Aquat Conserv Mar Freshw Ecosyst. 2021, 31: 2804-2813.

[45] Trizio I, Crestanello B, Galbusera P, Wauters LA, Tosi G, Matthysen E, Hauffe HC. Geographical distance and physical barriers shape the genetic structure of Eurasian red squirrels (Sciurus vulgaris) in the Italian Alps. Mol Ecol. 2005, 14: 469-481.

[46] Turtle Conservation Coalition. Turtles in trouble: the world’s 25+ most endangered tortoises and freshwater turtles-2011. 2011, Lunenburg, Turtle Conservation Coalition

[47] Turtle Taxonomy Working Group. Turtles of the world: annotated checklist and atlas of taxonomy, synonymy, distribution, and conservation status (9th Ed.). Chelonian Res Monogr. 2021, 8: 1-472

[48] van Dijk PP. The status of turtles in Asia. Chelonian Res Monogr. 2000, 2: 15-23

[49] Waples RS, Do C. Linkage disequilibrium estimates of contemporary N_e using highly variable genetic markers: a largely untapped resource for applied conservation and evolution. Evol Appl. 2010, 3: 244-262.

[50] Weissenbacher A, Preininger D, Ghosh R, Morshed AGJ, Praschag P. Conservation breeding of the Northern river terrapin Batagur baska at the Vienna Zoo, Austria, and in Bangladesh. Int Zoo Yearb. 2015, 49: 31-41.

[51] Wringe BF, Stanley RRE, Jeffery NW, Anderson EC, Bradbury IR. HYBRIDDETECTIVE: a workflow and package to facilitate the detection of hybridization using genomic data in R. Mol Ecol Res. 2017, 17: e275-e284.

[52] Yuan X, Hu Q, Bu R, Yang J, Lin L, Shi HT. Microhabitat selection by the Beale’s eyed turtle (Sacalia bealei) and conservation implications. Glob Ecol Conserv. 2025, 60: e03607

[53] Zheng X, Levine D, Shen J, Gogarten S, Laurie C, Weir B. A high-performance computing toolset for relatedness and principal component analysis of SNP data. Bioinformatics. 2012, 283326-3328.

Funding

Lingnan University

RIGHTS & PERMISSIONS

The Author(s)