High genetic diversity and demographic history of captive Siamese and Saltwater crocodiles suggest the first step toward the establishment of a breeding and reintroduction program in Thailand
High genetic diversity and demographic history of captive Siamese and Saltwater crocodiles suggest the first step toward the establishment of a breeding and reintroduction program in Thailand
Sorravis Lapbenjakul 0 1
Watcharaporn Thapana 0 1
Panupon Twilprawat 0 1
Narongrit Muangmai 1
Thiti Kanchanaketu 1
Yosapong Temsiripong 1
Sasimanas Unajak 1
Surin Peyachoknagul 0 1
Kornsorn Srikulnath 0 1
0 Laboratory of Animal Cytogenetics and Comparative Genomics (ACCG), Department of Genetics, Faculty of Science, Kasetsart University , Chatuchak, Bangkok , Thailand , 2 Animal Breeding and Genetics Consortium of Kasetsart University (ABG-KU), Kasetsart University , Bangkok , Thailand , 3 Center for Advanced Studies in Tropical Natural Resources, National Research University-Kasetsart University (CASTNAR, NRU-KU), Kasetsart University , Bangkok , Thailand , 4 Department of Fishery Biology, Faculty of Fisheries, Kasetsart University , Bangkok , Thailand , 5 Division of Genetics, Department of Science, Faculty of Liberal Arts and Science, Kasetsart University (Kamphaeng Saen Campus), Kamphaeng Saen, Nakhon Pathom, Thailand, 6 R&D Center, Sriracha Moda Co., Ltd. , Sriracha, Chonburi , Thailand , 7 Department of Biochemistry, Faculty of Science, Kasetsart University , Bangkok , Thailand , 8 Department of Biology, Faculty of Science, Naresuan University , Muang Phitsanulok, Phitsanulok , Thailand , 9 Center of Excellence on Agricultural Biotechnology: (AG-BIO/PERDO-CHE) , Bangkok , Thailand
1 Editor: Michael A. Russello, University of British Columbia Okanagan , CANADA
The Siamese crocodile (Crocodylus siamensis) and Saltwater crocodile (C. porosus) are two of the most endangered animals in Thailand. Their numbers have been reduced severely by hunting and habitat fragmentation. A reintroduction plan involving captive-bred populations that are used commercially is important and necessary as a conservation strategy to aid in the recovery of wild populations. Here, the genetic diversity and population structure of 69 individual crocodiles, mostly members of captive populations, were analyzed using both mitochondrial D-loop DNA and microsatellite markers. The overall haplotype diversity was 0.924±0.971 and the mean expected heterozygosity across 22 microsatellite loci was 0.578±0.701 for the two species. This agreed with the star-like shaped topology of the haplotype network, which suggests a high level of genetic diversity. The mean ratio of the number of alleles to the allelic range (M ratio) for the populations of both species was considerably lower than the threshold of 0.68, which was interpreted as indicative of a historical genetic bottleneck. Microsatellite markers provided evidence of introgression for three individual crocodiles, which suggest that hybridization might have occurred between C. siamensis and C. porosus. D-loop sequence analysis detected bi-directional hybridization between male and female individuals of the parent species. Therefore, identification of genetically non-hybrid and hybrid individuals is important for long-term conservation man-
Data Availability Statement: All relevant data are
within the paper and its Supporting Information
Funding: This study was financially supported by
grants from Kasetsart University Research and
Development Institute (KURDI; grant no. 80.60),
Kasetsart University 72 year anniversary graduate
scholarship (no. 10/2557) from the Graduate
School, Kasetsart University, Professor Motivation
agement. Relatedness values were low within the captive populations, which supported
(PM; no. PM4/2558), Special Track Staff (STS; no.
STS1/2558), and the Science Research Fund
(ScRF; no. ScRF-S19-2558) from the Faculty of
Science, Kasetsart University, the Center for
Advanced Studies in Tropical Natural Resources,
National Research University-Kasetsart University
(CASTNAR, NRU-KU, Thailand; no. 10/2559), the
graduate scholarship provided by the National
Research Council of Thailand (NRCT) as of fiscal
year 2017, the Graduate School, Kasetsart
University, and the Science Achievement
Scholarship of Thailand (SAST; no. 5717400381)
from the Office of the Higher Education
Commission, Thailand. This research is also
partially supported by the Center of Excellence on
Agricultural Biotechnology, Science and
Technology Postgraduate Education and Research
Development Office, Office of Higher Education
Commission, Ministry of Education (AG-BIO/
PERDO-CHE). The funders had no role in study
design, data collection and analysis, decision to
publish, or preparation of the manuscript.
Yosapong Temsiripong is employed by Sriracha
Moda Co., Ltd. Sriracha Moda Co., Ltd provided
support in the form of salary for author YT, but did
not have any additional role in the study design,
data collection and analysis, decision to publish, or
preparation of the manuscript. The specific roles of
this author is articulated in the `author
Competing interests: Yosapong Temsiripong is
employed by Sriracha Moda Co., Ltd. There are no
patents, products in development or marketed
products to declare. This does not alter our
adherence to all the PLOS ONE policies on sharing
data and materials.
their genetic integrity and the viability of a breeding and reintroduction management plan.
This work constitutes the first step in establishing an appropriate source population from a
scientifically managed perspective for an in situ/ex situ conservation program and
reintroduction of crocodile individuals to the wild in Thailand.
The Siamese crocodile (Crocodylus siamensis) is a freshwater species that is commonly found
in swamps and sheltered portions of rivers and streams in Thailand [
], whereas the Saltwater
crocodile (C. porosus) inhabits more marine environments and typically resides in saline and
brackish mangrove swamps, estuaries, deltas, lagoons, and lower stretches of rivers [
However, habitat destruction and illegal hunting have resulted in fragmented crocodile populations
in the wild [
]. The severe decline of populations to fewer than 200 individuals of C.
siamensis and a report from a preliminary survey of only three C. porosus observed in the wild Thai
populations [4,5] led to both species being listed as Critically Endangered by the Crocodile
Specialist Group of the Species Survival Commission of the International Union for
Conservation of Nature (IUCN) and in Appendix I of the Convention on International Trade in
Endangered Species of Wild Fauna and Flora (CITES) [
]. Reintroductions of captive-bred
individuals and in situ/ex situ management are necessary for the recovery of wild populations
of C. siamensis and C. porosus, with decisions made at the national policy level. A large number
of C. siamensis individuals are maintained in 12 commercial captive facilities under the
auspices of the Crocodile Co-operatives of Thailand and CITES regulations for the leather and
food industries, whereas C. porosus has been managed as a captive population without
commercial use. However, no new individuals have been introduced into these captive populations
subsequent to their establishment and no historical information on the captive populations is
available, which might lead to the loss of genetic variation. This might consequently hamper
the adaptability of a captive-bred population to a new environment, the long-term
sustainability of the population, and increase the probability of species extinction . A captive
population should be fundamentally managed to retain maximum genetic variation by avoiding
inbreeding and outbreeding depression . Evaluation of the genetic diversity in captive
populations thus provides important information for prospective breeding and reintroduction
programs, and in situ/ex situ management.
Conservation management of animal wildlife is impacted by hybridization, which can result
from the possible invasion of one species into the habitat of another or involuntary mixing of
species in the same farm. Such hybridization might lead to the genetic extinction of the species
]. Hybridization has been widely reported in wild Crocodylus populations between Morelet's
crocodile (C. moreletti) and the American crocodile (C. acutus) and between the Cuban
crocodile (C. rhombifer) and American crocodile [10±15]. Hybrids between Siamese and Cuban
crocodiles or between Siamese and Saltwater crocodiles as a result of anthropogenic impacts
have been observed, despite the different chromosome constitutions of these species [16±20].
Most Crocodylus anthropogenic hybrids are a serious problem in the context of conservation
management because the hybrids are highly similar in morphology to the parental species and
might lead to introgression if included in a reintroduction program. Consequently, a genetic
tool must be developed to identify and exclude hybrids before reintroduction.
Assessment of genetic diversity in natural and captive populations is an important step to
understand population structure, history, and hybrid status better when developing breeding
2 / 21
programs for conservation management of threatened species [21±23]. Molecular genetic
markers such as mitochondrial DNA (mtDNA) and biparentally inherited nuclear DNA
microsatellites can facilitate the ability to characterize population diversity, assign possible
origins of individuals, and identify hybrids and their parents [
]. Although the 12
captive crocodile populations in Thailand contain both C. siamensis and C. porosus, the C.
siamensis captive-bred population is currently considered the most important source of
individuals for a reintroduction program. The captive crocodile populations were derived from twelve
major crocodile farms. These distribute crocodiles to sub-farms throughout the country to
increase the Thai captive crocodile population. However, the present lack of genetic
information impedes and undermines establishment of an effective long-term conservation
management plan. In the study described herein, we examined genetic diversity and relatedness
between captive-bred populations of C. siamensis and C. porosus, as well as wild C. siamensis
individuals in the Boraphet Wetland Wildlife Reserve. We used mtDNA D-loop and
microsatellite genotyping to assess the diversity of the gene pool to guide reintroductions to the wild
and the formulation of in situ/ex situ management recommendations. The admixture between
C. siamensis and C. porosus in captive populations is also discussed.
Materials and methods
Animal material and DNA extraction
Samples were collected from 47 adult C. siamensis and 17 adult C. porosus from 12 captive
populations in nine provinces of Thailand (Nakhon Ratchasima: 14Ê58050@N, 102Ê06000@E;
Lopburi: 14Ê48000@N, 100Ê37037@E; Bangkok: 13Ê45014@N, 100Ê30005@E; Ayutthaya: 14Ê20052@N,
100Ê33038@E; Chonburi: 13Ê21040.11@N, 100Ê59004.82@E; Chainat: 15Ê11010@N, 100Ê07024@E;
Chachoengsao: 13Ê41025@N, 101Ê04013@E; Saraburi: 14Ê31059@N, 100Ê55000@E; Ratchaburi: 13Ê
32008@N, 99Ê48048@E). Permission was granted by the farm owners and the Crocodile
Co-operatives of Thailand. Five C. siamensis were captured from the Boraphet Wetland Wildlife
Reserve (Nakhon Sawan: 15Ê42048@N, 100Ê08007@E), and all crocodiles were released
immediately in the same area after sample collection. This research was conducted under the authority
of the Department of Fisheries, Ministry of Agriculture and Cooperatives, Thailand. Detailed
information on the sampled individuals is presented in S1 Table and S1 Fig. Individuals were
classified as C. siamensis or C. porosus on the basis of morphology [
] and molecular
mtDNA markers . A piece of scale clipped from the tail of each specimen was collected as a
source of DNA [
]. Whole genomic DNA was extracted in accordance with the standard
salting-out protocol as described previously [
]. DNA quality and quantity were determined by
electrophoresis on 1% agarose gels and spectrophotometric analysis. Animal care and all
experimental procedures were approved by the Animal Experiment Committee, Kasetsart
University, Thailand (approval no. ACKU04959), and conducted in accordance with the
Regulations on Animal Experiments at Kasetsart University.
Mitochondrial D-loop DNA fragments were amplified using the primers mtCytbf2
(50TGCCATGTTCGCATCCATCC-30) and mt12srRNAr2 (50-CCAGAGGCTAGGCGTCGTGG-30),
which were designed based on five crocodilian mtDNA sequences: C. siamensis (GenBank
accession number: EF581857), C. porosus (GenBank accession number: AJ810453), C. niloticus
(GenBank accession number: AJ810452), Alligator mississippiensis (GenBank accession
number: Y13113), and Gavialis gangeticus (GenBank accession number: AB079596). PCR
amplification was performed using 15 μl of 1× ThermoPol buffer that contained 1.5 mM MgCl2, 0.2
mM dNTPs, 5.0 μM primers, 0.5 U of Taq polymerase (Apsalagen Co. Ltd., Bangkok,
3 / 21
Thailand), and 25 ng of genomic DNA. The PCR conditions were as follows: initial
denaturation at 94ÊC for 3 min, followed by 40 cycles of 94ÊC for 30 s, 55ÊC for 30 s, and 72ÊC for 1
min 30 s, and then a final extension at 72ÊC for 5 min. The PCR products were detected by
electrophoresis on 1% agarose gels. The PCR products were cloned using the pTG19-T vector
(Vivantis Technologies Sdn Bhd, Selangor Darul Ehsan, Malaysia). The nucleotide sequences
of the DNA fragments were determined by the DNA sequencing service of First Base
Laboratories Sdn Bhd (Seri Kembangan, Selangor, Malaysia). The BLASTn and BLASTx programs
(http://blast.ncbi.nlm.nih.gov/Blast.cgi) were used to search nucleotide sequences in the
National Center for Biotechnology Information (NCBI) database to confirm the identity of the
DNA fragments amplified in the present study. The sequences generated in this study were
deposited in the DNA Data Bank of Japan (DDBJ) (S1 Table).
D-loop sequence analysis
Multiple sequence alignment was performed for 69 sequences generated in this study (52 for
C. siamensis and 17 for C. porosus individuals), 20 sequences from other crocodile species
available in the GenBank database, and one sequence from a turtle Pelodiscus sinensis
(GenBank accession number: AY687385) as the outgroup, using the default parameters of the
Molecular Evolutionary Genetics Analysis 6 (MEGA6) software (Center for Evolutionary
Functional Genomics, The Biodesign Institute, Tempe, AZ, USA [
]). All unalignable and
gap-containing sites were carefully removed and trimmed from the data sets. Estimates of
haplotype (h) and nucleotide (π) diversity [
] were calculated based on mtDNA D-loop
sequences as implemented in DnaSP version 5 [
]. Tests of neutral sequence evolution,
namely Tajima's D [
] and Fu and Li's D and F tests [
], were performed using DnaSP
version 5, and Fu's Fs [
] was calculated using Arlequin version 184.108.40.206 [
]. Significance values
for differences among these test results were determined using 10,000 coalescent simulations
in accordance with the recommended parameters for the software. Phylogenetic analysis was
performed using Bayesian inference (BI) with MrBayes version 3.2.6 [
]. The best-fit model
of DNA substitution was determined for each genetic region using Kakusan4 [
Markov chain Monte Carlo process was used to run four chains simultaneously for one million
generations. After the log-likelihood value stabilized, a sampling procedure was performed
every 100 generations to obtain 10,000 trees, from which a majority-rule consensus tree with
average branch lengths was generated. All sample points prior to attaining convergence were
discarded as burn-in, and the Bayesian posterior probability in the sampled tree population
was calculated as a percentage. A statistical parsimony network of the consensus sequences
was constructed using the Templeton, Crandall and Sing (TCS) algorithm implemented in
PopART version 1.7. [
All 22 microsatellite primer sets that were developed originally from the Saltwater crocodile
(S2 Table) [
] were used for genotype determination in C. siamensis and C. porosus. Most
of the microsatellite loci were located in different linkage groups of the Saltwater crocodile
genome. PCR amplification was performed using 15 μl of 1× ThermoPol buffer that contained
1.5 mM MgCl2, 0.2 mM dNTPs, 5.0 μM primers, 0.5 U of Taq polymerase (Apsalagen Co. Ltd.,
Bangkok, Thailand), and 25 ng of genomic DNA. The PCR conditions were as follows: initial
denaturation at 95ÊC for 3 min; followed by four cycles of 95ÊC for 20 s, 65ÊC for 20 s, and
72ÊC for 30 s; followed by four cycles of 95ÊC for 20 s, 62ÊC for 20 s, and 72ÊC for 30 s;
followed by eight cycles of 95ÊC for 20 s, 60ÊC for 20 s, and 72ÊC for 30 s; followed by 24 cycles of
4 / 21
95ÊC for 20 s, 55ÊC for 20 s, and 72ÊC for 30 s; and a final extension at 72ÊC for 7 min. The
PCR products were separated on a 6% denaturing polyacrylamide gel (w/v) and visualized by
silver staining as described previously . Allele sizes were measured using the GelAnalyzer
software (http://www.gelanalyzer.com/index.html). DNA fragments were extracted from the
silver-stained gels and cloned into the pTG19-T vector (Vivantis Technologies). Nucleotide
sequencing was carried out to confirm the identity of the DNA fragments amplified in the
present study. Nucleotide sequences of all microsatellite loci for C. siamensis and C. porosus
were deposited in DDBJ (S3 Table).
Microsatellite data analysis
Allelic frequency, number of specific alleles, observed heterozygosity (Ho), expected
heterozygosity (He), Hardy±Weinberg equilibrium, and linkage disequilibrium were calculated for
each locus and for each captive population using Arlequin version 220.127.116.11. Shannon's
information index (I) was calculated for each locus for each species using GENALEX version 6.5
]. Polymorphic information content was estimated using the Excel Microsatellite Toolkit
]. Reduction in heterozygosity due to non-random mating (FST) was estimated to
determine pairwise population differentiation with corrected P values using Arlequin version
18.104.22.168. The state of heterozygosity excess and shift in allelic frequency distributions in
genetically bottlenecked populations was tested using BOTTLENECK version 1.2.02 [
The Wilcoxon signed rank test with a two-phased model of mutation (TPM) and stepwise
mutation model (SMM) was used to obtain probabilities for excess levels of heterozygosity
due to the small sample sizes of loci and small sample size. The TPM was carried out with
95% single-step mutations and 5% multistep mutations, and the variance among multiple
steps was set at 12 [
]. This test detects relatively short-term bottleneck events. To test for
relatively long-term bottleneck events, the M ratio test [
] was performed using Arlequin
version 22.214.171.124. The M ratio is the mean number of alleles in a population divided by the
allelic size range, and can indicate reductions in both recent and historical population sizes.
Phylogenetic analysis of all microsatellite loci was performed using the unweighted pair
group with arithmetic mean (UPGMA) clustering method using NTSYSpc version 2.1
(Exeter Software, New York, USA). Jaccard's coefficient was used to estimate genetic similarity
]. Bootstrap analysis was performed using FreeTree software with 500 replicates [
Mantel's test was used to determine the goodness of fit for a cluster analysis. The degree of
fit was interpreted as follows: 0.9 r (very good fit); 0.8 r < 0.9 (good fit); 0.7 r < 0.8
(poor fit); r < 0.7 (very poor fit). Principal component analysis (PCA) was performed to
visualize the overall relationship across individuals in the populations using GENALEX
The model-based clustering method implemented in STRUCTURE version 2.3.3 [
used to determine population structure. Run length was set to 100,000 Markov chain Monte
Carlo replicates after a burn-in period of 100,000 generations using correlated allelic
frequencies under a straight admixture model. The number of clusters (K) was varied from 1 to 28,
with 25 replicates for each value of K. The most likely number of clusters was determined by
plotting the log probability of the data (ln Pr(X|K) [
] across the range of K values tested and
selecting K at which the value of ln Pr(X|K) stabilized. The ΔK method [
] was applied with
STRUCTURE HARVESTER [
]. In addition, relatedness values (r) [
] were estimated
among individual crocodiles using GENALEX version 6.5. The individual and overall
inbreeding coefficients with 95% confidence intervals were calculated using Ritland's estimator [
implemented in COANCESTRY .
5 / 21
D-loop haplotype variation
After editing, sequences for the 1,110 to 1,528 bp fragment of both C. siamensis and C.
porosus D-loop into 263 bp in length, which corresponded to the fragment at positions 16,525
to 16,807 bp of the C. siamensis mitochondrial genome (GenBank accession number:
EF581859), were compiled into data sets for C. siamensis and C. porosus. This sequence
region was located in D-loop domain III, which is commonly used for the analysis of genetic
diversity and population structure [
]. We identified 152 polymorphic sites (comprising 84
transitions and 68 transversions) in the C. siamensis data set and 41 polymorphic sites
(comprising 22 transitions and 19 transversions) in the C. porosus data set. The number of
haplotypes was 35 for C. siamensis and 14 for C. porosus. The overall haplotype and
nucleotide diversities were 0.924 ± 0.031 and 0.031 ± 0.009, respectively, for C. siamensis and
0.971 ± 0.032 and 0.021 ± 0.006, respectively, for C. porosus (S4 Table).
The haplotype networks for both C. siamensis and C. porosus showed a star-like shaped
topology (Fig 1), which indicates a high level of haplotype diversity. All C. siamensis
haplotypes clustered in a group that was distinct from the C. porosus haplotypes by at least 34
mutational steps with a missing haplotype (Fig 1). The most common haplotype of C.
siamensis in the population sample (CSI44) differed from shared nucleotide sequences of the
haplotype of the mitochondrial genome of C. siamensis (GenBank accession number: EF581859)
that had been described previously [
] by only one mutational step. An additional common
haplotype of C. siamensis (CSI38) was distinguished from CSI44 by one mutational step.
However, one haplotype (CSI17) was distinguished from CSI44 by at least 50 mutational
steps with a missing haplotype. The most common haplotype of C. porosus differed from
the haplotype of the mitochondrial genome of C. porosus (GenBank accession number:
DQ273698) that had been described previously [
] by two mutational steps. Phylogenetic
analysis of a combined data set for the D-loop sequences from both C. siamensis and C.
porosus, together with those for 20 crocodile species obtained from the GenBank database,
indicated that most sequences of C. siamensis and all sequences of C. porosus each formed a
monophyletic clade. However, CSI17 was placed as a sister clade to C. rhombifer (GenBank
accession number: NC_024513) (Fig 2).
Tajima's D (−2.711, P < 0.001), Fu's Fs (−18.013, P < 0.001), Fu and Li's F (−5.099,
P < 0.02), and Fu and Li's D (−5.209, P < 0.02) were all negative and significant for the
Dloop sequences of C. siamensis. Similarly, Tajima's D (−2.240, P < 0.01), Fu's Fs (−6.188,
P < 0.02), Fu and Li's F (−3.322, P < 0.02), and Fu and Li's D (−3.123, P < 0.02) were all
negative and significant for the mtDNA D-loop sequences of C. porosus.
Genotypic variation, hybridization, and relatedness
Twenty-two microsatellite primer pairs were applied to genotype the crocodile individuals.
A total of 299 alleles, which comprised 138 private alleles in C. siamensis and 82 private
alleles in C. porosus, were detected among all loci, with a mean number of alleles per locus of
13.59 (S5 Table). The two crocodile species shared 79 alleles. Allelic frequencies showed
significant departures from Hardy±Weinberg expectations at five loci for the C. siamensis
population and eight loci for the C. porosus population with multiple lines of evidence for
linkage disequilibrium (S5±S7 Tables). However, the ability to detect significant departures
from Hardy±Weinberg equilibrium was low because of the small sample sizes. No consistent
patterns of deviation from Hardy±Weinberg equilibrium or linkage equilibrium were
detected across sites. Consequently, genetic analyses were then performed, based on all the
6 / 21
Fig 1. Haplotype network based on sequence data for the mitochondrial DNA D-loop region of
Siamese and Saltwater crocodiles constructed using statistical parsimony with the Templeton,
Crandall, and Sing (TCS) algorithm. The numbers of individuals that possessed a haplotype is indicated by
the different colors inside the circles. Inferred but unsampled haplotypes are indicated by slashes. Missing
haplotypes are indicated by a black circle. Detailed information for all crocodile individuals is presented in S1
7 / 21
Fig 2. Phylogenetic relationships among mitochondrial DNA D-loop region sequences inferred using Bayesian
inference analysis. Support values at each node are the Bayesian posterior probability. Detailed information for all crocodile
individuals is presented in S1 Table.
8 / 21
microsatellite loci. The polymorphic information content of both crocodile species ranged
from 0.038 to 0.937, and Shannon's information index ranged from 0.108 to 3.148 (S5
Table). The Ho values of C. siamensis ranged from 0.000 (CpP3001) to 1.000 (CpP4501)
(mean 0.486 ± 0.306) and He values ranged from 0.038 (CpP3008) to 0.949 (CpP501) (mean
0.578 ± 0.323). The Ho values of C. porosus ranged from 0.059 (CpP3001) to 0.941
(CpP1409) (mean 0.591 ± 0.253) and He values ranged from 0.059 (CpP3001) to 0.939
(CpP501) (mean 0.701 ± 0.221) (S5 and S8 Tables). After 110 permutations, estimates of FST
showed significant differences between captive and wild populations of crocodiles, which
indicated substantial genetic subdivision (S9 and S10 Tables). In the test for population
bottlenecks, SMM and TPM were 0.997 and 0.982, respectively, in C. siamensis (normal
Lshaped mode shift), and 0.999 and 0.996, respectively, in C. porosus (normal L-shaped mode
shift) as determined by Wilcoxon sign-rank tests (S11 Table). Although some populations
showed evidence of reduction in population size, overall results did not detect a significant
reduction in C. siamensis or C. porosus. However, the M ratio across all populations
averaged 0.264 ± 0.191 for C. siamensis and 0.252 ± 0.155 for C. porosus (S5 Table). These M
ratio values were lower than the 0.68 threshold identified by Garza and Williamson [
which indicated a historical population reduction.
Phylogenetic analysis of the 69 samples divided the crocodile individuals into two major
groups (C. siamensis and C. porosus) (Fig 3). The cophenetic coefficient of Mantel's test
indicated a very good fit for the cluster analysis (r = 0.959). The C. porosus group comprised 16
C. porosus individuals and two C. siamensis individuals (CSI05 and CSI06), whereas the C.
siamensis group contained 50 C. siamensis individuals and one C. porosus individual (CPO09)
(Fig 3). PCA revealed that the first, second, and third principal components accounted for
26.12%, 4.85%, and 4.47% of the total variation, respectively, and provided support for the
distinction between the two crocodile groups and admixture of the aforementioned three
individuals (CSI05, CSI06, and CPO09) (Fig 4).
Structure analysis revealed the highest posterior probability with one peak (K = 2) on the
basis of Evanno's ΔK (Fig 5a), with all crocodiles grouped into two clusters, α and β, which
corresponded to C. porosus and C. siamensis, respectively (Fig 5b). Two C. siamensis individuals
(CSI05 and CSI06) and one C. porosus individual (CPO09) showed high estimated levels of
admixture, whereas low levels of admixture were estimated in CPO14, CSI08, CSI17, CSI18,
CSI26, and CSI41. By contrast, STRUCTURE analysis based on the mean ln P(K) revealed one
peak (K = 13) (Fig 6a), which provided evidence for 13 clusters (Fig 6b). Clusters A±E
represented C. porosus, whereas clusters F±M comprised C. siamensis. A pairwise relatedness test
was performed to determine the level of relatedness between individuals in the study
population (Fig 7). The relatedness values (r) of 2,346 pairs of crocodiles among the 69 sampled
crocodiles were examined, and the mean pairwise value was −0.009. A total of 1,449 pairs showed
r < 0; there were 892 pairs with 0.5 < r < 0; and five pairs with 1< r < 0.5 (S12 Table). The
mean inbreeding coefficient was 0.186, with individual inbreeding coefficients ranging from
−0.016 to 0.824 (S13 Table).
The meat, clothing, and medicinal products that are derived from C. siamensis constitute 1%
of the agricultural income in Thailand [
]. Crocodile numbers are increasing in all Thai
crocodile farms, which is contrary to the decreasing number of individuals observed in the wild.
Consequently, CITES regulations permit only limited exportation of captive-bred crocodile
products until the reintroduction program and in situ/ex situ management are sustainable.
9 / 21
Fig 3. Microsatellite-based phylogenetic relationships for 52 Siamese crocodiles and 17 Saltwater crocodiles generated by the unweighted
pair group with arithmetic mean (UPGMA) clustering method. Support values at each node are bootstrap values. ª*º indicates a bootstrap
value < 50%. The genetic similarity matrices are shown by Jaccard's coefficient. Detailed information for all crocodile individuals is presented in S1
The first step toward establishment of a captive breeding and reintroduction program is the
evaluation of genetic diversity based on analysis of mtDNA D-loop sequences and
microsatellite markers as an indicator of the genetic status of captive crocodile populations in the main
10 / 21
Fig 4. Principal component analysis of 69 crocodiles (52 Siamese crocodiles and 17 Saltwater crocodiles) using
microsatellite data across 22 loci. Blue diamonds indicate Siamese crocodile samples. Red diamonds indicate
Saltwater crocodile samples. The yellow diamond indicates CSI17, which might be a Cuban crocodile or a hybrid between
Siamese and Cuban crocodiles.
Diversity and population structure
Although nucleotide diversity was low, analysis of mtDNA D-loop sequences indicated that
haplotype diversity was relatively high in the sampled individuals of C. porosus and C.
siamensis. This finding was probably attributable to the high haplotype heterogeneity in the study
population. Concomitantly, microsatellite analysis showed that the mean Ho and He were
0.486 and 0.578, respectively, for C. siamensis and 0.591 and 0.701, respectively, for C. porosus.
These values are similar to those of the Siamese crocodile population in Vietnam (Ho = 0.416
and He = 0.420) [
] and the Saltwater crocodile population in Palau (Ho = 0.570 and He =
], which suggests that the populations of both species show a state of high genetic
diversity. Differentiation among populations was also detected in both captive and wild
populations of C. siamensis. FST values were 0.073±0.240 (P < 0.05) from 50 pairwise captive/wild
population comparisons in C. siamensis and 0.200±0.356 (P < 0.05) from five pairwise captive
population comparisons in C. porosus, which might indicate high heterozygosity. This suggests
that captive-bred Thai populations of both C. siamensis and C. porosus can be recommended
for use in breeding programs and reintroduction plans.
Population history of captive C. siamensis and C. porosus
Network analyses revealed the existence of four different haplotype groups, of which two were
in the C. siamensis haplotype group. The most common haplotype (CSI44) was closely related
to the haplotype of three of the five individuals (CSI48, CSI49, and CSI52) that were sampled
from the Boraphet Wetland Wildlife Reserve. Many unique haplotypes were also detected,
which suggests that populations of C. siamensis had been isolated from each other. The
starlike shaped topology of the haplotype network of C. siamensis suggests a high level of genetic
heterogeneity. A similar pattern was observed for the C. porosus haplotype network, although
wild individuals were not available for study. This star-like network might also imply the
presence of a population bottleneck followed by a population expansion. Significant negative
values for neutrality statistics also indicate population expansion [
]. This finding suggests that
the captive crocodile populations originated from wild-caught crocodiles and were
subsequently bred within each captive breeding program to produce a large number of individuals.
11 / 21
Fig 5. Population structure of Siamese and Saltwater crocodiles from 69 crocodile individuals. (a) Evanno's ΔK graph. (b) STRUCTURE bar
plots depicting the model-based clustering results for inferred K = 2. Inferred genetic clusters are displayed as different colors. Each vertical bar on the
x-axis represents an individual, and the y-axis presents the proportion of membership (posterior probability) in each genetic cluster. Recovered
crocodile species, including clusters α and β, are superimposed on the plot, with black vertical lines indicating the boundaries. Detailed information for
all crocodile individuals is presented in S1 Table.
Unfortunately, the original source of most captive-bred populations is unknown. Although a
demographic bottleneck was not supported by the bottleneck test, the M ratio showed a signal
for historical population reduction for both C. siamensis and C. porosus. Collectively, these
results suggest that these captive crocodile populations underwent a recent population
expansion. Crocodile hunting to collect individuals has occurred over the last 50 years [
subsequent successful breeding and an increase in the size of captive populations. Low
haplotype heterogeneity is usually observed with rapid population expansion . However, in the
present study, captive populations of both C. siamensis and C. porosus showed high haplotype
heterogeneity. This suggests that individuals from various historically distinct lineages in the
wild might have contributed to the established populations, leading to the presence of
individuals with high levels of admixture. Range expansion and high haplotype diversity were aided
by intentional human introduction in each captive population.
12 / 21
Fig 6. Population structure of Siamese and Saltwater crocodiles from 69 crocodile individuals. (a) Mean Ln P(K) graph. (b) STRUCTURE
bar plots depict the model-based clustering results for inferred K = 13. Inferred genetic clusters are displayed as different colors. Each vertical bar
on the x-axis represents an individual, and the y-axis presents the proportion of membership (posterior probability) in each genetic cluster.
Recovered crocodile species, including clusters A±M, are superimposed on the plot, with black vertical lines indicating the boundaries. Detailed
information for all crocodile individuals is presented in S1 Table.
A missing haplotype was detected between most C. siamensis and C. porosus haplotype
groups, which determined that the two crocodile species were genetically distinct [
However, an additional missing haplotype was detected between CSI17 and all the remaining
C. siamensis haplotypes. This finding corresponded with the results of the microsatellite
phylogenetic analysis in which CSI17 was not placed in the cluster that contained most C. siamensis
individuals. Cuban crocodiles have also been identified in captive breeding programs in
Cambodia and Vietnam [
], and the admixture of crocodile individuals in Cambodia, Vietnam,
and Thailand might result from partial introgression in captive breeding programs over a long
]. The present results suggest that CSI17 is at least the result of hybridization
with a C. rhombifer, if not a pure Cuban crocodile.
13 / 21
Fig 7. Observed distribution of pairwise relatedness for 52 Siamese crocodiles and 17 Saltwater crocodiles plotted against expected
Hybridization between C. siamensis and C. porosus
The present microsatellite phylogenetic analyses and PCA indicated that three individuals
(CPO09, CSI05, and CSI06) might be the result of interspecific hybridization between C.
siamensis and C. porosus. Cluster analysis using the STRUCTURE software can help to determine
the degree of hybridization by aggregating individuals into a single cluster relative to
additional highly differentiated populations/species [
]. Using a posterior probability of 0.95 as
a criterion for assignment to a pure species as estimated with STRUCTURE , the three
individuals (CPO09, CSI05, and CSI06) showed evidence of admixture. Both CPO09 and
CSI05 were placed in strongly differentiated clusters (A and F) that were distinct from most
C. siamensis and C. porosus individuals, respectively, and possessed a high proportion of
private alleles. By contrast, CSI06 clustered with C. porosus. Observations on captive crocodiles
suggest that bi-directional hybridization occurs between C. siamensis and C. porosus [
conclusion is in agreement with the results of the D-loop sequence analysis. CPO09 shared a
mitochondrial haplotype with C. porosus, which suggests that hybridization probably occurred
between a female C. porosus and a male C. siamensis. By contrast, CSI05 and CSI06 shared
haplotypes with C. siamensis, which suggest that they are descended from a female C. siamensis
and a male C. porosus.
Natural hybridization often occurs in the genus Crocodylus [
]. However, hybrids
between C. siamensis and C. porosus are always observed in the presence of anthropogenic
impacts. This might reflect behavioral and geographic differences between the two species
]. Moreover, the occurrence of fertile hybrids might be rare as a consequence of
the different chromosome constitutions of the two species [
]. Hybridization is a serious
conservation concern when anthropogenic factors cause misclassification of individuals on the
basis of morphology, such as between Siamese and Saltwater crocodiles [
]. This might
lead to backcrossing, resulting in localized hybrids, widespread introgression, or complete
14 / 21
admixture . The removal of hybrids is beneficial to sustainable use programs because local
commercial captive operations can utilize the genetic material of hybrid crocodiles for
industry. Alternatively, we strongly suggest that efforts to avoid hybridization are taken into account
in conservation management and reintroduction programs because hybrids often show
superior survival and adaptive mechanisms under competition with non-hybrid individuals [
This scenario also elicits the likelihood of genetic pollution of the species by more abundant
The present D-loop and microsatellite data indicate that the captive populations of C.
siamensis and C. porosus are genetically divergent, with partial introgression and hybrids between
the two species. An alternative approach is required to identify hybrid or non-hybrid
individuals in addition to karyotyping [
]. A larger population sample and more detailed analysis
are required to estimate the degree of hybridization, such as analysis of F1, F2, or backcross
generations. The costs of molecular analytic methods have decreased greatly and the use of
such procedures should be a prerequisite for breeding management and the establishment of a
reintroduction program, rather than a karyological analysis approach. It is strongly
recommended that researchers conduct routine molecular analyses to identify individuals of hybrid
Implications for conservation and management
Inbreeding has negative impacts on both reproduction and survival [
]. The mean
inbreeding coefficient was 0.137 for C. siamensis and 0.334 for C. porosus, and a high level of
mean relatedness was indicated (0.008 and 0.109 for both species, respectively), implying that
the samples were not closely related . These results suggest that the populations of C.
porosus and C. siamensis were founded initially by individuals of unknown ancestry. Mate pairings
within representative captive crocodile populations might be determined to aid breeding and
reintroduction programs. In general, breeding that is based on the minimization of relatedness
generates maximum genetic diversity [
]. Prior to the current study, captive breeding and
reintroduction programs for C. siamensis and C. porosus in Thailand have proceeded in the
absence of genetic diversity data, and have relied solely on morphological, ethological,
demographic, and logistic information for the implementation of short-term management strategies
]. The present results provide an important genetic baseline for in situ/ex situ
management decisions. The populations of both species exhibited high heterozygosity, which could be
an indication of the sound establishment of the captive species populations. Specifically, 10
individuals of C. siamensis (19.231%) comprising six males and four females, and seven
individuals of C. porosus (43.750%) consisting of four males and three females, were highlighted as
genetically important for conservation and management, because they showed values of
relatedness that were lower than that of the overall captive population. This provides flexibility in
the implementation of breeding and reintroduction programs. However, one possible concern
for reintroduction is the likelihood of disrupting locally adapted genetic materials through
outbreeding among different source populations . This remains an inherent problem because
captive crocodile populations typically contain individuals that originate from widely scattered
populations. Consideration of the minimization of relatedness and relationship of haplotype
lineage with the source population are necessary to propose management action. Finding
suitable habitats for the introduction of the two species is an additional problem in Thailand [
Although the present study is preliminary, the results are important for ongoing conservation
and genetic management programs both locally and throughout the distribution ranges of the
species. We believe the substantial captive crocodile populations that were sampled in this
study to be broadly representative of captive-bred populations in Thailand as a whole. This
15 / 21
work constitutes the first step in establishing an appropriate source population from a
scientifically managed perspective for in situ/ex situ conservation and reintroduction programs in
Thailand. In addition, the present results provide reference data for further characterization of
cryptic diversity, which directly impacts on the conservation prioritization of C. siamensis and
C. porosus outside Thailand, such as in Southeast Asia.
The results of the present study indicate the status of genetic diversity in the extant Thai
captive crocodile population. Accurate information on captive populations for breeding
programs, reintroduction, or in situ/ex situ management will aid the management of subsequent
generations and maintain sustainable genetic diversity for long-term survival of the
population. Therefore, it is highly recommended that the genetic status of C. siamensis and C. porosus
populations are closely monitored using molecular genetic methods to improve in situ/ex situ
management. We are also convinced that the future conservation management of C. siamensis
will require the genetic identification of non-hybrid wild populations. Additional studies
involving genome-wide scans will be required for an improved understanding of overall gene
functions in different populations.
S1 Fig. Map showing the collection sites for the Crocodylus siamensis and C. porosus
specimens. Numbers indicate sample locality. Detailed information for all crocodile individuals is
presented in S1 Table.
S1 Table. Summary of crocodile specimens.
S2 Table. List of microsatellite primers and sequences used in the study.
S3 Table. Sequence accession numbers of 22 microsatellite loci from one representative
individual of the Siamese crocodile (Crocodylus siamensis) and one representative
individual of the Saltwater crocodile (C. porosus).
S4 Table. Mitochondrial DNA D-loop diversity based on a 263-bp fragment for the
Siamese crocodile (Crocodylus siamensis) and Saltwater crocodile (C. porosus). P < 0.001,
P < 0.02, P < 0.01, and ns = not significant.
S5 Table. Genetic diversity of 52 individuals of the Siamese crocodile (Crocodylus
siamensis) and 17 individuals of the Saltwater crocodile (C. porosus) based on 22 microsatellite
S6 Table. Pairwise differentiation of linkage disequilibrium among Siamese crocodile
(Crocodylus siamensis) individuals based on 22 microsatellite loci. The number indicates
P values, with 110 permutations.
S7 Table. Pairwise differentiation of linkage disequilibrium among Saltwater crocodile
(Crocodylus porosus) individuals based on 22 microsatellite loci. The number indicates
P values, with 110 permutations.
16 / 21
S8 Table. Comparison of observed and expected heterozygosity of the Siamese crocodile
(Crocodylus siamensis) and Saltwater crocodile (C. porosus) based on 22 microsatellite loci
in each captive/wild population. Detailed information for all crocodile individuals is
presented in S1 Table.
S9 Table. Pairwise genetic differentiation (FST) between Siamese crocodile (Crocodylus
siamensis) captive/wild populations based on 22 microsatellite loci. The number indicates P
values, with 110 permutations. Detailed information for all crocodile individuals is presented
in S1 Table.
S10 Table. Pairwise genetic differentiation (FST) between Saltwater crocodile (Crocodylus
porosus) captive populations based on 22 microsatellite loci. The number indicates P values,
with 110 permutations. Detailed information for all crocodile individuals is presented in S1
S11 Table. Test for genetic bottlenecks in the Siamese crocodile (Crocodylus siamensis) and
Saltwater crocodile (C. porosus) using BOTTLENECK version 1.2.02 and calculation of the
M ratio using Arlequin 126.96.36.199 for all populations. Detailed information for all crocodile
individuals is presented in S1 Table.
S12 Table. Pairwise genetic relatedness (r) for all 69 crocodile individuals. Detailed
information for all crocodile individuals is presented in S1 Table.
S13 Table. Pairwise inbreeding coefficients for all 69 crocodile individuals. Detailed
information for all crocodile individuals is presented in S1 Table.
We thank the crocodile farm members of the Crocodile Co-operatives of Thailand (COCOT),
Department of Fisheries, Ministry of Agriculture and Cooperatives (Thailand), and Faculty of
Veterinary Science, Mahidol University (Thailand) for advice on sample preparation. We are
also grateful to Siwapech Sillapaprayoon (Kasetsart University, Thailand) for helpful
discussions. This research is partially supported from graduate scholarship provided by the National
Research Council of Thailand (NRCT) as of fiscal year 2017, the Graduate School, Kasetsart
University and the Center of Excellence on Agricultural Biotechnology, Science and
Technology Postgraduate Education and Research Development Office, Office of Higher Education
Commission, Ministry of Education (AG-BIO/PERDO-CHE).
Conceptualization: Sorravis Lapbenjakul, Kornsorn Srikulnath.
Data curation: Sorravis Lapbenjakul, Kornsorn Srikulnath.
Formal analysis: Sorravis Lapbenjakul, Narongrit Muangmai, Thiti Kanchanaketu, Kornsorn
17 / 21
Funding acquisition: Sorravis Lapbenjakul, Watcharaporn Thapana, Yosapong Temsiripong,
Investigation: Sorravis Lapbenjakul, Narongrit Muangmai, Thiti Kanchanaketu, Yosapong
Temsiripong, Kornsorn Srikulnath.
Methodology: Sorravis Lapbenjakul, Narongrit Muangmai, Thiti Kanchanaketu, Yosapong
Temsiripong, Kornsorn Srikulnath.
Project administration: Kornsorn Srikulnath.
Resources: Sorravis Lapbenjakul, Yosapong Temsiripong, Kornsorn Srikulnath.
Software: Sorravis Lapbenjakul, Narongrit Muangmai, Thiti Kanchanaketu, Kornsorn
Supervision: Kornsorn Srikulnath.
Validation: Sorravis Lapbenjakul, Narongrit Muangmai, Thiti Kanchanaketu, Kornsorn
Visualization: Sorravis Lapbenjakul, Watcharaporn Thapana, Panupon Twilprawat,
Narongrit Muangmai, Thiti Kanchanaketu, Yosapong Temsiripong, Sasimanas Unajak, Surin
Peyachoknagul, Kornsorn Srikulnath.
Writing ± original draft: Sorravis Lapbenjakul, Kornsorn Srikulnath.
Writing ± review & editing: Sorravis Lapbenjakul, Watcharaporn Thapana, Panupon
Twilprawat, Narongrit Muangmai, Thiti Kanchanaketu, Yosapong Temsiripong, Sasimanas
Unajak, Surin Peyachoknagul, Kornsorn Srikulnath.
18 / 21
19 / 21
20 / 21
1. Platt SG , Lynam AJ , Temsiripong Y , Kampanakngarn M. Occurrence of the Siamese crocodile (Crocodylus siamensis ) in Kaeng Krachan National Park, Thailand. Nat. Hist. Bull. Siam Soc . 2002 ; 50 : 7± 14 .
2. Crocodile Specialist Group. Crocodylus porosus . The IUCN Red List of Threatened Species 1996 [cited 28 April 2017 ]. http://dx.doi.org/10.2305/IUCN.UK. 1996 .RLTS.T5668A11503588.en.
3. Ross JP . Crocodiles: Status survey and Conservation Action Plan. 2nd Edition . Gland, Switzerland and Cambridge, UK: IUCN/SSC Crocodile Specialist Group. IUCN. 1998 .
4. Ratanakorn P , Amget B , Ottley B . Preliminary surveys of crocodiles in Thailand . In Crocodiles. Proceedings of the 12th Working Meeting of the Crocodile Specialist Group . 1994 ; 35 ± 49 .
Wongsongsarn C. Annual Survey Report of the Siamese crocodiles (Crocodylus siamensis) in Bueng Boraphet Non-Hunting Area . Department of Fisheries. Bangkok, Thailand. 2010 .
6. Simpson BK , Bezuijen MR . Siamese crocodile Crocodylus siamensis . In Crocodiles: Status Survey and Conservation Action Plan , Mandis SC , Stevenson C . (eds). IUCN Crocodile Specialist Group: Darwin, Australia. 2010 ; 120 ± 126 .
7. Allendorf FW , Lundquist LL . Introduction: population biology, evolution, and control of invasive species . Conserv. Biol . 2003 ; 17 : 24 ± 30 .
Conserv . 2011 ; 20 : 1843 ± 1861 .
9. Allendorf FW , Leary RF , Spruell P , Wenburg JK . The problems with hybrids: setting conservation guidelines . Trends Ecol. Evol . 2001 ; 16 : 613 ± 622 .
10. Ray DA , Dever JA , Platt SG , Rainwater TR , Finger AG , McMurry ST , et al. Low levels of nucleotide diversity in Crocodylus moreletii and evidence of hybridization with C. acutus . Conserv. Genet . 2004 ; 5 : 449 ± 462 .
11. Cedeño-VaÂzquez JR , Rodriguez D , CalmeÂ S , Ross JP , Densmore LD , Thorbjarnarson JB . Hybridization between Crocodylus acutus and Crocodylus moreletii in the Yucatan Peninsula: I. Evidence from mitochondrial DNA and morphology . J. Exp. Zool. A. Ecol. Genet. Physiol . 2008 ; 309 : 661 ± 673 . https:// doi.org/10.1002/jez.473 PMID: 18626922
12. Rodriguez D , Cedeño-Vazquez JR , Forstner MRJ , Densmore LD . Hybridization between Crocodylus acutus and Crocodylus moreletii in the Yucatan Peninsula: II. Evidence from microsatellites . J. Exp. Zool. A. Ecol. Genet. Physiol . 2008 ; 309 : 674 ± 686 . https://doi.org/10.1002/jez.499 PMID: 18800373 Weaver JP , Rodriguez D , Venegas-Anaya M , Cedeño-VaÂzquez JR , Forstner MR , Densmore LD .
Genetic characterization of captive Cuban crocodiles (Crocodylus rhombifer) and evidence of hybridization with the American crocodile (Crocodylus acutus) . J. Exp. Zool. A. Ecol. Genet. Physiol . 2008 ; 309 : 649 ± 660 . https://doi.org/10.1002/jez.471 PMID: 18646197
14. MiliaÂn-GarcÂõa Y , Venegas-Anaya M , Frias-Soler R , Crawford AJ , Ramos-Targarona R , RodrÂõguezSoberoÂn R , et al. Evolutionary history of Cuban crocodiles Crocodylus rhombifer and Crocodylus acutus inferred from multilocus markers . J. Exp. Zool. A. Ecol. Genet. Physiol . 2011 ; 315 : 358 ± 375 . https://doi. org/10.1002/jez.683 PMID: 21462351
15. MiliaÂn-GarcÂõa Y , Ramos-Targarona R , PeÂrez-Fleitas E , Sosa-RodrÂõguez G , Guerra-Manchena L , Alonso-Tabet M , et al. Genetic evidence of hybridization between the critically endangered Cuban crocodile and the American crocodile: implications for population history and in situ/ex situ conservation . Heredity . 2015 ; 114 : 272 ± 280 . https://doi.org/10.1038/hdy. 2014 .96 PMID: 25335559
16. Cohen MM , Gans C. The chromosomes of the order Crocodilia . Cytogenet. Genome Res . 1970 ; 9 : 81 ± 105 . PMID: 5461083
17. Thang NQ. The status of Crocodylus rhombifer in the Socialist Republic of Vietnam . In Crocodiles, Proceedings of the 12th Working Meeting of the Crocodile Specialist Group . 1994 ; 141 ± 142 .
18. Chavananikul V , Suwattana D , Koykul W , Wattanodorn S , Sukkai J. A research report on karyotypes of freshwater crocodiles (Crocodylus siamensis), Saltwater crocodiles (Crocodylus porosus) and the interspecific hybrids by conventional and banding techniques . Chulalongkorn University, Bangkok, Thailand. 1998 .
19. Fitzsimmons NN , Buchan JC , Lam PV , Polet G , Hung TT , Thang NQ , Gratten J . Identification of purebred Crocodylus siamensis for reintroduction in Vietnam . J. Exp. Zool. A. Ecol. Genet. Physiol . 2002 ; 294 : 373 ± 381 . https://doi.org/10.1002/jez.10201 PMID: 12461816
20. Srikulnath K , Thapana W , Muangmai N. Role of chromosome changes in Crocodylus evolution and diversity . Genomics Inform . 2015 ; 13 : 102 ± 111 . https://doi.org/10.5808/GI. 2015 . 13 .4.102 PMID: 26865840
21. Moritz C . Conservation units and translocations: strategies for conserving evolutionary processes . Hereditas . 1999 ; 130 : 217 ± 228 .
22. Benavides E , Russello M , Boyer D , Wiese RJ , Kajdacsi B , Marquez L , et al. Lineage identification and genealogical relationships among captive GalaÂpagos tortoises . Zoo Biol . 2012 ; 31 : 107 ± 120 . https:// doi.org/10.1002/zoo.20397 PMID: 21674601
23. Meraner A , Cornetti L , Gandolfi A . Defining conservation units in a stocking-induced genetic melting pot: unraveling native and multiple exotic genetic imprints of recent and historical secondary contact in Adriatic grayling . Ecol. Evol . 2014 ; 4 : 1313 ± 1327 . https://doi.org/10.1002/ece3.931 PMID: 24834328
24. Funk SM , Fa JE . Phylogeography of the endemic St. Lucia whiptail lizard Cnemidophorus vanzoi: conservation genetics at the species boundary . Conserv. Genet . 2006 ; 7 : 651 ± 663 .
25. Ramirez O , Altet L , Enseñat C , Vilà C , Sanchez A , Ruiz A . Genetic assessment of the Iberian wolf Canis lupus signatus captive breeding program . Conserv. Genet . 2006 ; 7 : 861 ± 878 .
26. Ross FD , Mayer GC . On the dorsal armor of the Crocodilia . Advances in herpetology and evolutionary biology . 1983 ; 305 ± 331 .
27. Dinets V. Long-distance signaling in Crocodylia . Copeia. 2013 ; 2013 : 517 ± 526 .
28. Srikulnath K , Thongpan A , Suputtitada S , Apisitwanich S . New haplotype of the complete mitochondrial genome of Crocodylus siamensis and its species specific DNA markers: distinguishing C. siamensis from C. porosus in Thailand . Mol. Biol. Rep . 2012 ; 39 : 4709 ± 4717 . https://doi.org/10.1007/s11033-011 - 1263-7 PMID: 21947949
29. Barrow D , Halford A. SOP No 8.3: Biopsy pole collection of tissue samples from free-ranging crocodiles . Department of Parks and Wildlife Animal Ethics Committee. Australia , August 2016 .
30. Supikamolseni A , Ngaoburanawit N , Sumontha M , Chanhome L , Suntrarachun S , Peyachoknagul S , et al. Molecular barcoding of venomous snakes and species-specific multiplex PCR assay to identify snake groups for which antivenom is available in Thailand . Genet. Mol. Res . 2015 ; 14 : 13981 ± 13997 . https://doi.org/10.4238/ 2015 .October. 29 .18 PMID: 26535713
31. Tamura K , Stecher G , Peterson D , Filipski A , Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0 . Mol. Biol. Evol . 2013 ; 30 : 2725 ± 2729 . https://doi.org/10.1093/molbev/mst197 PMID: 24132122
32. Nei M. Molecular evolutionary genetics . New York: Columbia University Press; 1987 .
33. Librado P , Rozas J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data . Bioinformatics . 2009 ; 25 : 1451 ± 1452 . https://doi.org/10.1093/bioinformatics/btp187 PMID: 19346325
34. Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism . Genetics . 1989 ; 123 : 585 ± 595 . PMID: 2513255
35. Fu YX , Li WH . Statistical tests of neutrality of mutations . Genetics . 1993 ; 133 : 693 ± 709 . PMID: 8454210
36. Fu YX . Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection . Genetics . 1997 ; 147 : 915 ± 925 . PMID: 9335623
37. Excoffier L , Lischer HE . Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows . Mol. Ecol. Resour . 2010 ; 10 : 564 ± 567 . https://doi.org/10.1111/j. 1755- 0998 . 2010 . 02847 . x PMID : 21565059
38. Huelsenbeck JP , Ronquist F. MRBAYES : Bayesian inference of phylogeny . Bioinformatics . 2001 ; 17 : 754 ± 755 .
39. Tanabe AS . Kakusan4 and Aminosan: two programs for comparing nonpartitioned, proportional and separate models for combined molecular phylogenetic analyses of multilocus sequence data . Mol. Ecol. Resour . 2011 ; 11 : 914 ± 921 . https://doi.org/10.1111/j.1755- 0998 . 2011 . 03021 . x PMID : 21592310
40. Clement M , Snell Q , Walker P , Posada D , Crandall K. TCS: Estimating gene genealogies . Parallel and Distributed Processing Symposium , International Proceedings . 2002 ; 2 : 184 .
41. Miles LG , Isberg SR , Glenn TC , Lance SL , Dalzell P , Thomson PC , et al. A genetic linkage map for the saltwater crocodile (Crocodylus porosus) . BMC Genomics . 2009 ; 10 : 339 . https://doi.org/10.1186/ 1471 -2164-10-339 PMID: 19640266
42. Miles LG , Isberg SR , Moran C , Hagen C , Glenn TC. 253 Novel polymorphic microsatellites for the saltwater crocodile (Crocodylus porosus) . Conserv. Genet . 2009 ; 10 : 963 ± 980 .
43. Benbouza H , Jacquemin JM , Baudoin JP , Mergeai G . Optimization of a reliable, fast, cheap and sensitive silver staining method to detect SSR markers in polyacrylamide gels . Biotechnol. Agron. Soc . 2006 ; 10 : 77 ± 81 .
44. Peakall ROD , Smouse PE . GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and researchÐan update . Bioinformatics . 2012 ; 28 : 2537 ± 2539 . https://doi.org/10.1093/ bioinformatics/bts460 PMID: 22820204
45. Park S. Trypanotolerance in West African cattle and the population genetic effects of selection . Ph D thesis , University of Dublin; 2001 .
46. Cornuet JM , Luikart G. Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data . Genetics . 1996 ; 144 : 2001 ± 2014 . PMID: 8978083
47. Piry S , Luikart G , Cornuet JM . BOTTLENECK: a computer program for detecting recent reductions in the effective population size using allele frequency data . J. Hered . 1999 ; 90 : 502 ± 503 .
48. Garza JC , Williamson EG . Detection of reduction in population size using data from microsatellite loci . Mol. Ecol . 2001 ; 10 : 305 ± 318 . PMID: 11298947
49. Jaccard P. Nouvelles recherches sur la distribution florale . Bull. Soc. Vaud. Sci Nat . 1908 ; 44 : 223 ± 270 .
50. Hampl V , PavlÂõcek A , Flegr J . Construction and bootstrap analysis of DNA fingerprinting-based phylogenetic trees with a freeware program FreeTree: application to trichomonad parasites . Int. J. Syst. Evol. Microbiol . 2001 ; 51 : 731 ± 735 . https://doi.org/10.1099/ 00207713 -51-3-731 PMID: 11411692
51. Pritchard JK , Stephens M , Donnelly P . Inference of population structure using multilocus genotype data . Genetics . 2000 ; 155 : 945 ± 959 . PMID: 10835412
52. 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 . https://doi.org/10.1111/j. 1365 - 294X . 2005 . 02553 . x PMID : 15969739
53. Earl DA , vonHoldt BM . STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method . Conserv. Genet. Resour . 2012 ; 4 : 359 ± 361 .
54. Ritland K. Estimators for pairwise relatedness and individual inbreeding coefficients . Genet. Res . 1996 ; 67 : 175 ± 185 .
Wang J. COANCESTRY : a program for simulating, estimating and analyzing relatedness and inbreeding coefficients . Mol. Ecol. Resourc . 2011 ; 11 : 141 ± 145 . https://doi.org/10.1111/j.1755- 0998 . 2010 .
02885. x PMID : 21429111
56. Ray DA , Densmore L . The crocodilian mitochondrial control region: general structure, conserved sequences, and evolutionary implications . J. Exp. Zool. A Ecol. Genet. Physiol . 2002 ; 294 : 334 ± 345 . https://doi.org/10.1002/jez.10198 PMID: 12461813
57. Li Y , Wu X , Ji X , Yan P , Amato G. The complete mitochondrial genome of salt-water crocodile (Crocodylus porosus) and phylogeny of crocodilians . J. Genet. Genomics . 2007 ; 34 : 119 ± 128 . https://doi.org/10. 1016/S1673- 8527 ( 07 ) 60013 - 7 PMID: 17469784
58. Department of International Trade Promotion, Ministry of Commerce, Thailand. Foreign Buyers. Statistics. 2017 . http://www.ditp.go.th/main.php?filename=statistics
59. Russello MA , Brazaitis P , Gratten J , Watkins-Colwell GJ , Caccone A . Molecular assessment of the genetic integrity, distinctiveness and phylogeographic context of the Saltwater crocodile (Crocodylus porosus) on Palau . Conserv. Genet . 2007 ; 8 : 777 ± 787 .
60. Nwafili SA , Gao TX . Genetic diversity in the mtDNA control region and population structure of Chrysichthys nigrodigitatus from selected Nigerian rivers: implications for conservation and aquaculture . Arch. Pol. Fish . 2016 ; 24 : 85 ± 97 .
61. Thorbjarnarson J. Crocodiles: an action plan for their conservation . Messel H, King FW , Ross JP . (eds). IUCN/SSC Crocodile Specialist Group, Gland, Switzerland. 1992 .
62. Martin G , Pear JJ . Behavior modification: What it is and how to do it . Psychology Press; 2015 .
63. Meganathan PR , Dubey B , Batzer MA , Ray DA , Haque I. Molecular phylogenetic analyses of genus Crocodylus (Eusuchia, Crocodylia, Crocodylidae) and the taxonomic position of Crocodylus porosus . Mol. Phylogenet. Evol . 2010 ; 57 : 393 ± 402 . https://doi.org/10.1016/j.ympev. 2010 . 06 .011 PMID: 20601012
64. Jelden DC , Manolis C , Giam H , Thomson J , Lopez A . Crocodile Conservation and Management in Cambodia: a Review with Recommendations . IUCN Crocodile Specialist Group, Darwin, Australia; 2005 .
65. Starr A , Daltry JC , Nhek R . DNA study reveals C. siamensis at the Phnom Tamao Wildlife Rescue Centre, Cambodia. Newsletter of the IUCN/ SSC Crocodile Specialist Group. 2010 ; 28 : 5 ± 7 .
66. Youngprapakorn P. Hybrid CpsY . The Samutprakan Crocodile Farm and Zoo , Thailand, April 1990 , 3 ± 13 .
67. Lacy RC . Clarification of genetic terms and their use in the management of captive populations . Zoo Biol . 1995 ; 14 : 565 ± 578 .
68. Frankham R . Genetics and extinction . Biol. Conserv . 2005 ; 126 : 131 ± 140 .
69. Temsiripong Y. Reintroduction of the Siamese crocodile . Crocodile Specialist Group Newsletter. 2001 ; 20 : 10 ± 12 .
70. Temsiripong Y. Re-introduction of captive-raised Siamese crocodile in Thailand . Re-introduction Specialist Group Newsletter. 2007 ; 26 : 55 ± 57 .
71. Storfer A . Gene flow and endangered species translocations: a topic revisited . Biol. Conserv . 1999 ; 87 : 173 ± 180 .