Next Article in Journal
Strategies and Mechanisms of Thermal Compensation in Newborn Water Buffaloes
Previous Article in Journal
The Effect of Using Bovine Colostrum and Probiotics on Performance, Egg Traits, Blood Biochemical and Antioxidant Status of Laying Japanese Quails
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 13
Issue 13
10.3390/ani13132165
animals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Determining the Provenance of Traded Wildlife in the Philippines

by
Kate J. Brandis
1,*,
Phoebe Meagher
2,
Sabine Schoppe
3,
Kyle Zawada
4,
Indira Widmann
3,
Peter Widmann
3,
Roger G. Dolorosa
5 and
Roxane Francis
1
1
Centre for Ecosystem Science, School of Biological, Earth and Environmental Sciences, UNSW Sydney, Sydney 2052, Australia
2
Taronga Institute of Science and Learning, Taronga Conservation Society, Bradley’s Head Road, Mosman 2088, Australia
3
Katala Foundation Incorporated, Casoy Road, Purok El Rancho, Santa Monica, Puerto Princesa City 5300, Palawan, Philippines
4
Centre for Compassionate Conservation, University of Technology Sydney, Broadway, Sydney 2007, Australia
5
Puerto Princesa Campus, Western Philippines University, Santa Monica, Puerto Princesa City 5300, Palawan, Philippines
*
Author to whom correspondence should be addressed.
Animals 2023, 13(13), 2165; https://doi.org/10.3390/ani13132165
Submission received: 7 June 2023 / Revised: 26 June 2023 / Accepted: 28 June 2023 / Published: 30 June 2023
(This article belongs to the Section Wildlife)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

The illegal wildlife trade is a significant threat to biodiversity, with animals being taken from the wild and sold illegally. Often, species with the greatest demand are those that are already threatened, with illegal trade placing even more pressure on these populations. Knowing where traded animals or animal parts have come from is critically important for enforcement, education, and conservation actions. Using three highly traded Philippine species, the Palawan forest turtle, the Philippine cockatoo, and the Philippine pangolin, we developed a novel method of determining the geographic origin of animals or animal parts using a portable X-ray fluorescence device. This method can make significant contributions to combating the illegal wildlife trade through providing a rapid, non-destructive, cost-effective method for identifying the source of traded animals. End users of this tool include enforcement agencies, border force officials, wildlife rescue centres, and conservation researchers.

Abstract

The illegal wildlife trade is a significant threat to global biodiversity, often targeting already threatened species. In combating the trade, it is critical to know the provenance of the traded animal or part to facilitate targeted conservation actions, such as education and enforcement. Here, we present and compare two methods, portable X-ray fluorescence (pXRF) and stable isotope analysis (SIA), to determine both the geographic and source provenance (captive or wild) of traded animals and their parts. Using three critically endangered, frequently illegally traded Philippine species, the Palawan forest turtle (Siebenrockiella leytensis), the Philippine cockatoo (Cacatua haematuropygia), and the Philippine pangolin (Manis culionensisis), we demonstrate that using these methods, we can more accurately assign provenance using pXRF data ( x ¯   = 83%) than SIA data ( x ¯   = 47%). Our results indicate that these methods provide a valuable forensic tool that can be used in combating the illegal wildlife trade.

1. Introduction

The illegal wildlife trade is a significant threat to global biodiversity and a key driver of extinction risk, particularly as many targeted species are already endangered [1,2,3]. The interception of smuggled animals and animal parts is a common occurrence in the enforcement of wildlife trade legislation in many countries [4,5]. The ability to determine where the animal (or part) has come from, its provenance, is essential for conservation organisations, e.g., TRAFFIC, WWF, and IUCN, and enforcement agencies to address the source of the trade, allowing for targeted policy development, enforcement actions [6], repatriation of seized animals [7,8], education, and conservation actions [9,10].
The threat to biodiversity from illegal wildlife trade was formally recognised in 1973 with the signing of the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES), an intergovernmental agreement providing the framework for trade in flora and fauna that does not threaten survival of these species. Despite this, the illegal wildlife trade has continued to grow, facilitated by the globalisation of the world’s economy, flexible border arrangements, and improved transportation [11]. Significant declines in elephant Loxondonta africana and L. cyclotis numbers [12,13] and pangolins (Order: Pholidota) [14] and the extinction of the western black rhino Diceros bicornis longipes [15] are all attributed to demand on the illegal wildlife market. Ongoing demand for illegally traded flora and fauna is driven by traditional medicines, exotic pets, and animal parts, e.g., skins, furs [3,16].
In an attempt to alleviate the pressure on wild populations, the captive breeding of some species is permitted as defined in Article 1, paragraph (b) of the CITES convention, [17,18]. Under these circumstances, animals can be legally traded; however, ensuring that animals are not illegally taken from the wild and fraudulently traded as legal, i.e., laundered, is a problem that persists in many regions of the world [19,20,21]. Southeast Asia is recognised as a hotspot for illegally traded wildlife with its high levels of biodiversity, endemic species, and access to unregulated trade routes [11,22]. Our study region, Southern Palawan, Philippines, is a key source of animals in the trade of reptiles, birds, and pangolins [23,24]. Many of these animals are CRITICALLY ENDANGERED with very small populations remaining in the wild [25]. Captive breeding of some of these animals for release aims to bolster wild populations [26,27].
This study focused on three CRITICALLY ENDANGERED endemic Philippine species commonly traded in illegal wildlife markets. The Palawan forest turtle Siebenrockiella leytensis, a freshwater turtle endemic to the Province of Palawan [28,29], is the sixth most-confiscated species in the Philippines [30]. Demand for this species as pets and food have resulted in significant numbers being collected from the wild [29,31,32]. Seizure records from 2004–2018 indicated that 23 seizures collected 4723 turtles, with an estimated 1200 also sold illegally in China in 2015 [24]. This species is listed in CITES Appendix II [33], indicating it may become threatened with extinction if trade is not controlled. The Philippine cockatoo Cacatua haematuropygia is highly prized as a pet due to its mimicry skills [34]. This species has experienced severe population declines since the 1980s due to habitat loss, persecution, introduced infectious disease, and poaching for the pet trade [35,36]. It is listed in CITES Appendix I (CITES, 2021), indicating it is threatened with extinction and international trade is largely prohibited. The Philippine pangolin Manis culionensisis is restricted to the Palawan region, and is one of only eight species of its kind [14]. It is highly sought after for its meat and scales, used in traditional medicines [37]. This species is listed in CITES Appendix I [33]; however, trade is still occurring on a massive scale [37,38].
Using these three globally threatened species at risk from illegal trade, we aimed to test the efficacy of using a portable XRF device (pXRF) compared with stable isotopes (13C and 15N) as a method for determining provenance, both geographic and source (captive/wild). Portable instruments have the benefit of being more cost-efficient and enable the field collection of data, but they have yet to be applied in the determination of provenance to fight the illegal wildlife trade.
The pXRF collects data on the elemental composition of a sample for ~40 different elements. Provenance determination using element concentration data has been used in a range of fields, including food traceability and safety [39,40], archaeology [41], and conservation [42]. Methods for element measurement can be destructive and expensive, and require laboratories to conduct techniques such as inductively coupled plasma mass spectrometry (ICP-MS). Applications of instruments such as portable X-ray fluorescence (pXRF) devices are providing opportunities to determine animal provenance using rapid, non-destructive, cost-efficient techniques by measuring elemental abundances or concentrations within a sample. X-ray florescence technology includes laboratory-based instruments, e.g., ITRAX [43], or portable instruments [44,45]. Portable instruments have the benefit of being more cost-efficient, require minimal training for data collection, are battery powered, and can be used in a range of locations, e.g., Brazil [46] and Botswana [47], enabling field collection of data.
Elemental signatures, the unique combination of different element concentrations/abundances, are incorporated into tissue through diet, and are a representation of diet during the growth of the tissue. Keratinous tissues, e.g., feathers, scales, and scutes, provide a chemically inert record of diet over time [48] allowing for the study of provenance [48,49,50]. This is critical in determining the source of animals traded in illegal wildlife markets, and can contribute to mapping illegal wildlife trade networks [6,10] and the repatriation of animals [8].
Similarly, stable isotope ratios are determined by diet and can also provide information on provenance, allowing for the identification of geographic source location [51,52] and status as captive or wild animals [53,54]. Stable isotope analysis (SIA) is the measurement of ratios of different isotopes in a sample compared to natural occurrences. Biological studies typically use isotopes of carbon, oxygen, nitrogen, and hydrogen. Previous studies have used isotopes to distinguish between captive- and wild-bred fauna at high accuracies, including mink, frogs, and fish [55,56,57]. Commonly selected stable isotopes for the study of diet and provenance are 13C and 15N, due to their effectiveness in differentiating populations across a geographical range [53,58] and their direct relation to diet, indicating differences in terrestrial and aquatic diets, and trophic level and niche [59,60,61].
This study aims to demonstrate that pXRF is a valid tool that can be used in determining wildlife provenance, both geographic and source (captive/wild), with the hope that its portable nature will facilitate its use in combating the illegal wildlife trade with uptake by enforcement agencies for use in customs, wildlife markets, and pet stores. It aims to combat the illegal trade by identifying those animals or animal parts that have been fraudulently presented as legally obtained from captive breeding. It will also provide information for targeted enforcement actions [6], for repatriation of seized animals [7], and for education and conservation actions [9,10].
We hypothesise that pXRF data-derived models will perform better than SIA alone due to the quantity of data that can be collected through pXRF. The performance of geographic models compared to source (captive/wild) models will be influenced by the diet, environment, and life history traits of both the captive and wild species of interest.

2. Materials and Methods

2.1. Study Site and Sample Collection

Samples from the Palawan forest turtle, the Philippine cockatoo, and the Philippine pangolin were collected by Katala Foundation Incorporated (KFI), Western Philippines University (WPU), and Taronga Conservation Society across Palawan Island, Philippines (Figure 1, Table 1), in 2019. Samples were collected from a total of 45 Palawan forest turtles (24 captive, 21 wild) (Table 1). Captive animals were caught by hand from individual enclosures at the Katala Institute for Ecology and Biodiversity Conservation (KFI) in Narra, Palawan. Wild samples were collected opportunistically during survey field work undertaken by KFI in January 2019. Baited opera traps were set along creeks, semi-submerged, and secured to a nearby tree with rope [62]. Traps were re-checked during the evening and early morning. Scute samples from both captive and wild turtles were collected using a disposable biopsy punch. The edge of the scute was sampled using an upward scraping motion with uniform pressure, producing a long sample of an approximately 3 mm × 20 mm piece of keratin [63]. A single sample was taken from a marginal scute and stored refrigerated in a sealed, labelled cryovial.
Pangolin samples were a mix of scales (n = 6) and nails (n = 12), and all samples came from wild specimens. Scale samples were collected from live animals caught for population surveys. A small ‘v’-shaped notch was cut with a bone cutter and kept individually in a sample vial in a refrigerator. Samples were collected from the lateral, posterior scales, unless there was a damaged scale elsewhere on the animal that could be sampled. One sample per animal was collected. Nail samples were clipped from dead specimens that had been confiscated by the Palawan Council for Sustainable Development Staff (PCSDS) (Philippine government) and kept by the PCSDS in cold storage at their facilities in Irawan, Puerto Princesa City. Scales had been removed from these specimens and were not available for sampling.
Moulted feather samples were collected from within single and multiple bird aviaries at KFI (20 captive samples) and from wild individuals at roost and nest sites of the Philippine cockatoo (45 wild samples) (Table 1). Feathers were also collected opportunistically during routine health and population checks of wild birds. Where possible, feathers were collected from known individuals both captive and wild, i.e., single bird aviaries to avoid re-sampling of the same individuals and banded wild birds. For multi-bird aviaries, potential resampling was limited by only sampling a total of feathers for the total number of birds per aviary, e.g., 3 birds, 3 feathers. While this does not rule out potential resampling of the same individual, we attempted to limit the possibility and the impact it may have on data analyses. Single feathers from each wild roost/nest site were sampled to limit resampling of wild birds.
All samples were classified as being either captive or wild (i.e., animal had been bred or living in captivity for a period of time before sampling, or animal was collected from the wild), and labelled with their geographic provenance (i.e., wild animals collected from distinct regions, and captive animals assigned to their captive facility). Pangolin samples were classified into groups based on geographic location and body part. Although the pangolin sample size was small, we decided to include them in our analyses. Their elusive nature makes the sourcing of larger sample sizes difficult, they are very difficult to breed in captivity [64], and they are at high risk of extinction, making any samples valuable.
This study was conducted under the following permits: CITES permit WT2019–000849, University of New South Wales Animal Care and Ethics Committee number 18/127B, Gratuitous permit 2018–05 from the Palawan Council for Sustainable Development Staff (PCSDS).

2.2. X-ray Fluorescence

All samples were scanned for their elemental composition using an Olympus Vanta M-series portable XRF (pXRF, 50 kV, 80 uA) scanner using GeoChem (2) mode, with a 60 s scan time per beam (2 beams, 10 kV and 40 kV). Forty-two elements were measured and all were included in the models.). All samples were scanned within the Olympus Vanta Work Station. Samples were placed over the beam window ensuring the sample completely covered the scan area. This is a non-destructive analysis technique, and the same samples were then used for stable isotope analysis.

2.3. 13C and 15N Stable Isotope Analysis

Feather samples were cleaned in a 2:1 methanol:chloroform wash to remove any surface oils [65], then washed in mild detergent and rinsed three times in reverse osmosis (RO) water before being left to air-dry. Scutes, nails, and scales were cleaned of surface dirt and ground (using a ball grinder or mortar and pestle) to provide a homogenous sample. Samples were weighed into tin caps, and stable carbon (13C) and nitrogen (15N) isotopes were analysed at the Bioanalytical Mass Spectrometry Facility (BMSF) at the Mark Wainwright Analytical Centre (MWAC) University of New South Wales. A Delta V Advantage Isotope Ratio Mass Spectrometer and Flash 2000 Organic Elemental Analyzer fitted with a Conflo 4 were used.

3. Statistical Analyses

XRF beam spectra were analysed using the R package xrftools [66]. To build predictive models, we used the tree-based learning algorithm in the R package XGBoost [67]. We built models with two goals: predicting provenance (captive vs. wild) and predicting geographic location. These models were built for each species separately, in addition to separate models for the elemental pXRF and δ13C and δ15N isotope ratio data. XGBoost is an implementation of a gradient-boosting decision tree algorithm, which uses sequentially added models to correct errors made by existing models until no further improvements can be made [68]. We used default values of 6 for the maximum depth, 0.3 learning rate, a maximum of 10,000 training rounds, and the ‘multi:softprob’ objective. We used an 80% data split to create the training and test datasets, where the test set data were used to approximate ‘real-world’ model performance. Furthermore, the training data were split again to generate a watchlist dataset (20%) that the XGBoost algorithm assessed during training to avoid overfitting. The models produce probabilities for each class (wild vs. captive, geographical location), which we then converted to a prediction by taking the highest assigned probability as the prediction. To determine model accuracy, we compared the actual classification against the predicted class for the test set data. As model performance could vary based on the random split between training, watchlist, and test data, we ran the train–test process 50 times and aggregated the model performance into mean and standard deviation of accuracy. Variable importance was assessed by looking at the variables’ relative influence on the model outputs. As all pangolin samples were wild, we classified pangolins into 2 groups (Site F/claw or Site G/nail). We ran both geographic and captive/wild classification models on the turtle and cockatoo samples. We removed cockatoo samples that were labelled as wild but were housed within KFI for the testing of geographic origin.
To explore the isotopic differences between wild and captive turtles and cockatoos, we used simple linear models (excluding pangolins as we had no captive samples), separately testing for a significant difference (p < 0.05) in δ13C and δ15N based on their captive or wild origin. Again using separate linear models, we also tested for differences in δ13C and δ15N in sample geographic provenance for all three species, using the package emmeans to identify sites’ leading differences [69]. As pangolin samples were claws from one location and scales from another, we cannot be sure as to whether we are exploring geographic provenance or body part differences.

4. Results

4.1. X-ray Fluorescence

XGBoost predictive models were successful at differentiating between captive and wild specimens and at predicting geographic provenance using pXRF data (see Figure 2). Models predicted turtle samples as coming from a wild or captive specimen at 88 ± 12% accuracy, largely led by differences in Br, Pb, and Ni (Figure 2) based on scute elemental composition. Geographic provenance (n = 3 locations) was predicted at 94 ± 5% accuracy, led by differences in Br, Zr, and Y.
Philippine cockatoo feather classification models predicted geographic provenance (n = 3) at a 62 ± 9% accuracy rate. Key differences were found in the elements In, Mn, and Sr. Classification models also accurately assigned feathers into captive or wild classes with a mean accuracy of 78 ± 2%. Classification was largely led by differences in In, Sr, and Mn, with captive animals having higher concentrations (Figure 3).
Classification models using the elemental composition of pangolin samples were able to accurately predict groups (93 ± 14%). However, due to the samples available, we were unable to determine whether these differences are location- or body part-driven. Classification was largely led by differences in Ca, Zr, and V (Figure 4).

4.2. 13C and 15N Stable Isotopes

There were no significant differences in turtle scute δ15N (F(1,42) =0.81, p = 0.373, R2 = 0.02) or δ13C (F(1,43) = 0.04, p = 0.847, R2 = 0.0008) between wild and captive samples. There were no significant differences in turtle scute δ15N between sites (F(2,42) = 0.53, p = 0.59, R2 = 0.02). However, δ13C differed significantly between sites (F(2,42) = 3.95, p = 0.03, R2 = 0.16), led by differences between sites B and C (Figure 5a).
There were significant differences between captive and wild cockatoo feathers in both δ15N (F(1,28) = 10.336, p = 0.003, R2 = 0.26) and δ13C (F(1,28) = 12.97, p = < 0.001, R2 = 0.27) (Figure 5b). Within wild cockatoos, there were also significant differences in δ15N (F(2,22) = 4.49, p = 0.02, R2 = 0.29) when tested for geographic location. There were no significant differences in δ13C (F(2,22) = 0.21, p = 0.81, R2 = 0.02) (Figure 5b).
There were significant differences in Philippine pangolin δ15N groups (F(1,14) = 26.7, p = 0.0001, R2 = 0.65). δ13C also differed significantly between groups (F(1,14) = 9.10, p = 0.009, R2 = 0.35) (Figure 5c).

4.3. Predictive Stable Isotope Models

When using XGBoost predictive models and stable isotope results, we achieved 66 ± 16% accuracy predicting the captive or wild source of the Palawan forest turtle, led by δ15N, and 73 ± 14% when predicting geographic origin, led by δ13C. The models predicted the captive or wild status of the Philippine cockatoo at 28 ± 1% accuracy, led by differences in δ13C, and we achieved 27 ± 2% accuracy in predicting geographic provenance, led by differences in δ15N. Predictive models using Philippine pangolin stable isotopes performed poorly, with predicted geographic provenance at 41 ± 25% accuracy, led by δ15N.
Wild cockatoos from site E had a mean δ15N of 9.8‰ compared to site D’s 5.9‰ and captive animals’ 5.22‰ at KFI. Captive cockatoos had a less depleted mean δ13C (–19.03‰) than wild cockatoos, with δ13C −22.7‰ at site E and −21.93‰ at site D, indicating different sources of carbon.

5. Discussion

This study has demonstrated that it is possible to use data derived from a field-portable X-ray fluorescence device to accurately determine both geographic provenance and source (captive/wild) provenance of three Philippine species. Models using beam spectra data derived from pXRF were more accurate than those using stable isotope data (13C and 15N) for both geographic and source provenance models (Table 2). This suggests that pXRF technology shows promise as a low-cost, rapid, and portable solution for tackling the illegal wildlife trade.
The differences in prediction accuracy observed between the pXRF and SIA techniques are driven by the quantity of data collected through each technique. The ability of the pXRF to measure 42 elements in a single scan provides more sources of unique information compared to two stable isotopes, allowing models to utilise a richer dataset when training. It also uses a field-portable device, which does not require the use of large laboratory equipment, as is required for stable isotope analysis. It has the potential to be used on live animals, provided X-ray safety protocols are followed [70,71], reducing the need to invasively collect samples from animals, e.g., turtle carapaces may be scanned in situ without the need to remove scute samples. It can also be used by enforcement agencies in a range of locations including customs, wildlife markets, and pet stores.
Both the pXRF and SIA techniques measure values in tissues driven by diet [59,60], reflecting the habitat in which they grew the tissue. The species in this study represent a varied group of taxa, birds, mammals, and reptiles, and a range of feeding ecologies, from specialist pangolins [72], to omnivorous turtles [29], to frugivorous cockatoos [73]. These results are supported by other studies in a range of locations, including monotremes in Australia [54], terrestrial and marine waterbird species with diverse diets including piscivores, omnivores, and herbivores in Australia [49], and the southern Indian Ocean [74]. In addition, Buddhachat et al. (2016) successfully used pXRF to distinguish between elephant species, with results driven by differences in geographic location [75]. The study results have demonstrated that our technique is robust to varied feeding ecologies and taxa. The Palawan forest turtle model results (Table 2) were more accurate for the geographic provenance models than the captive/wild models. This was consistent for both pXRF- and SIA-derived data. This may reflect the relatively restricted home range individual turtles inhabit [76], ensuring that diet is representative of the local area only. However, both captive and wild Palawan forest turtles had similar ranges for δ15N and δ 13C values with large overlap (Figure 5), indicating similar sources of nitrogen and carbon in their diet. This highlights the benefits of using pXRF-derived data, providing a greater range of elements and enabling the detection of differences that otherwise might not be detected using SIA alone.
In contrast to the Palawan forest turtle, the Philippine cockatoo geographic source models performed less accurately than the captive/wild ones (Table 2). These birds have the ability to be wider-ranging than turtles, with the potential to be feeding across a larger geographic range and incorporating diet from a variety of locations and sources. The captive/wild models were more accurate, likely due to larger differences in diets between captive and wild animals, than between different wild locations. KFI cockatoo food sources included chicken feed, string beans, capsicum, and corn (S. Schoppe pers. comm.).
The Philippine pangolin models were limited to ‘Group’, representing both location and body part (claw or scale). While the pXRF predictive models for group performed well (93 ± 14%; Table 2), it is not possible to say whether this is due to locational differences or body part differences. The SIA data on their own did not result in accurate models (41 ± 25%; Table 2). Despite this limitation, the inclusion of these data is important as they represent the first collection of elemental data and add to the limited stable isotope data available for pangolins worldwide [77].
Other limitations of this this study include relatively small sample sizes, due to the limited population sizes of the globally threatened species we studied. Larger samples sizes would likely result in improved model accuracy and lower standard deviations. Another limitation is the potential resampling of individuals when using moulted feathers, particularly from multi-bird aviaries. As discussed in the methods section, we attempted to limit this. However, future studies should be aware of this restriction when using moulted feathers.
While this study achieved relatively accurate models for a range of taxa, the ability to do so for other taxa and species may be influenced by aspects that are diet-driven, e.g., life history traits, dietary breadth, dietary shifts, and home range sizes. Similarly, variations in captive diets and how closely they mimic wild diets will influence the ability to determine source provenance.
The further development of pXRF as a forensic tool, in addition to other tools such as SIA [78], genetics [4], and other molecular approaches [79], has significant potential and could achieve global application through the development of species reference libraries that can be collaboratively contributed to on an ongoing basis, providing data from a range of species and geographic locations.
The ability to predict geographic provenance at high accuracies using pXRF allows for the identification of animals and body parts, e.g., scales, horns, feathers, that are sourced from different geographic regions, e.g., Asian or African species [80]. This is of particular importance when only body parts are traded, which would otherwise require genetic testing, which can be expensive and time-consuming [81]. Knowledge regarding geographic and source provenance also enables anti-poaching and conservation efforts in specific geographic regions [10]. Predicting the captive or wild status of traded animals allows for differentiation between those individuals that were raised in captivity and those that were poached from the wild, facilitating the prosecution of illegal traders.
When combatting such a lucrative global problem as the illegal wildlife trade, it is important to have a range of forensic techniques to accurately identify wildlife as either captive or wild-sourced [54,82]. Furthermore, the ability to identify the location from which wild-sourced animals have been collected allows for the focusing of education [83,84], conservation [10], and enforcement actions [85]. Our study has demonstrated that pXRF provides a method that can achieve both these outcomes, and provides valuable additions to the knowledge of three highly traded Philippine species, contributing to the fight against the illegal wildlife trade and species conservation.

Author Contributions

Formal analysis, K.Z. and R.F.; Investigation, K.J.B., S.S., I.W. and P.W.; Data curation, K.J.B., P.M. and R.F.; Writing—original draft, K.J.B., P.M., K.Z. and R.F.; Writing—review & editing, K.J.B., S.S., K.Z. and R.G.D.; Visualization, R.F.; Supervision, S.S.; Funding acquisition, P.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

This study was carried out under the following permits: CITES permit WT2019-000849, University of New South Wales Animal Care and Ethics Committee number 18/127B, Gratuitous permit 2018–05 from the Palawan Council for Sustainable Development Staff (PCSDS).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the study animals being critically endangered.

Acknowledgments

The authors would like to thank the Palawan Council for Sustainable Development for the issuance of the permit, the access to specimens, and their continuous support. Thank you also to KFI staff who collected samples for this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zimmerman, M.E. The Black Market for Wildlife: Combating Transnational Organized Crime in the Illegal Wildlife Trade. J. Transnatl. Law 2003, 36, 1657–1689. [Google Scholar]
  2. Symes, W.S.; McGrath, F.L.; Rao, M.; Carrasco, L.R. The gravity of wildlife trade. Biol. Conserv. 2018, 218, 268–276. [Google Scholar] [CrossRef]
  3. Scheffers, B.R.; Oliveira, B.F.; Lamb, I.; Edwards, D.P. Global wildlife trade across the tree of life. Science 2019, 366, 71–76. [Google Scholar] [CrossRef]
  4. Rehman, A.; Jafar, S.; Ashraf Raja, N.; Mahar, J. Use of DNA Barcoding to Control the Illegal Wildlife Trade: A CITES Case Report from Pakistan. J. Bioresour. Manag. 2015, 2, 3. [Google Scholar] [CrossRef]
  5. Omifolaji, J.K.; Hughes, A.C.; Ibrahim, A.S.; Zhou, J.; Zhang, S.; Ikyaagba, E.T.; Luan, X. Dissecting the illegal pangolin trade in China: An insight from seizures data reports. Nat. Conserv. 2022, 46, 17–38. [Google Scholar] [CrossRef]
  6. Patel, N.G.; Rorres, C.; Joly, D.O.; Brownstein, J.S.; Boston, R.; Levy, M.Z.; Smith, G. Quantitative methods of identifying the key nodes in the illegal wildlife trade network. Proc. Natl. Acad. Sci. USA 2015, 112, 7948–7953. [Google Scholar] [CrossRef] [Green Version]
  7. Gaur, A.; Reddy, A.; Annapoorni, S.; Satyarebala, B.; Shivaji, S. The origin of Indian Star tortoises (Geochelone elegans) based on nuclear and mitochondrial DNA analysis: A story of rescue and repatriation. Conserv. Genet. 2005, 7, 231–240. [Google Scholar] [CrossRef]
  8. Gentile, G.; Ciambotta, M.; Tapia, W. Illegal wildlife trade in Galápagos: Molecular tools help the taxonomic identification of confiscated iguanas and guide their rapid repatriation. Conserv. Genet. Resour. 2013, 5, 867–872. [Google Scholar] [CrossRef]
  9. Kahler, J.S.; Gore, M.L. Beyond the cooking pot and pocket book: Factors influencing noncompliance with wildlife poaching rules. Int. J. Comp. Appl. Crim. Justice 2012, 36, 103–120. [Google Scholar] [CrossRef]
  10. Farine, D.R. Mapping illegal wildlife trade networks provides new opportunities for conservation actions. Anim. Conserv. 2020, 23, 145–146. [Google Scholar] [CrossRef] [Green Version]
  11. van Uhm, D. Illegal wildilfe trade to the EU and harms to the World. In Environmental Crime in Transnational Context: Global Issues in Green Enforcement and Criminology; Spapens, T., White, R., Huismen, W., Eds.; Routledge: Abingdon, UK, 2016. [Google Scholar]
  12. Poulsen, J.R.; Koerner, S.E.; Moore, S.; Medjibe, V.P.; Blake, S.A.; Clark, C.J.; Akou, M.E.; Fay, M.; Meier, A.; Okouyi, J.; et al. Poaching empties critical Central African wilderness of forest elephants. Curr. Biol. 2017, 27, R123–R138. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Maisels, F.; Strindberg, S.; Blake, S.; Wittemyer, G.; Hart, J.; Williamson, E.A.; Aba’a, R.; Abitsi, G.; Ambahe, R.D.; Amsini, F.; et al. Devastating decline of forest elephants in central Africa. PLoS ONE 2013, 8, e59469. [Google Scholar] [CrossRef] [PubMed]
  14. Challender, D.W.S.; Heinrich, S.; Shepherd, C.R.; Katsis, L.K.D. International trade and trafficking in pangolins, 1900–2019. In Biodiversity of the World: Conservation from Genes to Landscapes, Pangolins; Academic Press: Cambridge, MA, USA, 2020; pp. 259–276. [Google Scholar]
  15. Emslie, R. Diceros bicornis ssp. Longipes; The IUCN Red List of Threatened Species; IUCN: Gland, Switzerland, 2020; p. e.T39319A45814470. [Google Scholar] [CrossRef]
  16. ‘t Sas-Rolfes, M.; Challender, D.W.; Hinsley, A.; Veríssimo, D.; Milner-Gulland, E. Illegal wildlife trade: Scale, processes, and governance. Annu. Rev. Environ. Resour. 2019, 44, 201–228. [Google Scholar] [CrossRef] [Green Version]
  17. Thomson, J. Captive Breeding of Selected Taxa in Cambodia and Viet Nam: A Reference Manual for Farm Operators and Cites Authorities; TRAFFIC Southeast Asia, Greater Mekong Programme; TRAFFIC Southeast Asia: Ha Noi, Vietnam, 2008. [Google Scholar]
  18. Tensen, L. Under what circumstances can wildlife farming benefit species conservation? Glob. Ecol. Conserv. 2016, 6, 286–298. [Google Scholar] [CrossRef] [Green Version]
  19. Lyons, J.A.; Natusch, D.J.D. Wildlife laundering through breeding farms: Illegal harvest, population declines and a means of regulating the trade of green pythons (Morelia viridis) from Indonesia. Biol. Conserv. 2011, 144, 3073–3081. [Google Scholar] [CrossRef]
  20. Nijman, V.; Shepherd, C.R. Trade in non-native, CITES-listed, wildlife in Asia, as exemplifi ed by the trade in freshwater turtles and tortoises (Chelonidae) in Thailand. Contrib. Zool. 2007, 76, 207–211. [Google Scholar] [CrossRef] [Green Version]
  21. Moyle, B. Wildlife markets in the presence of laundering: A comment. Biodivers. Conserv. 2017, 26, 2979–2985. [Google Scholar] [CrossRef]
  22. Nijman, V. An overview of international wildlife trade from Southeast Asia. Biodivers. Conserv. 2009, 19, 1101–1114. [Google Scholar] [CrossRef] [Green Version]
  23. Cruz, R.M.; van den Beukel, D.; Lacerna-Widnmann, I.; Schoppe, S.; WIdmann, P. Wildlife Trade in Southern Palawan, Philippines. Banwa 2007, 4, 12–26. [Google Scholar]
  24. Sy, E.Y.; Schoppe, S.; Diesmos, M.L.L.; Lim, T.M.S.; Diesmos, A.C. Endangered by trade: Seizure analysis of the critically endangered Philippine Forest Turtle Siebenrockiella leytensis from 2004–2018. Philipp. J. Syst. Biol. 2020, 14, 1–8. [Google Scholar] [CrossRef]
  25. Schoppe, S.; Matillano, J.; Cervancia, M.; Acosta, D. Conservation Needs of the Critically Endangered Philippine Forest Turtle, Siebenrockiella leytensis, in Palawan, Philippines. Chelonian Conserv. Biol. 2010, 9, 145–153. [Google Scholar] [CrossRef]
  26. Chavez, L. Philippine Forest Turtles Stand a ‘Good Chance’ after First Wild Release. Available online: https://news.mongabay.com/2021/06/philippine-forest-turtles-stand-a-good-chance-after-first-wild-release/ (accessed on 28 June 2021).
  27. Griffiths, R.A.; Pavajeau, L. Captive breeding, reintroduction, and the conservation of amphibians. Conserv. Biol. 2008, 22, 852–861. [Google Scholar] [CrossRef] [PubMed]
  28. Asian Turtle Trade Working Group. Siebenrockiella Leytensis; The IUCN Red List of Threatened Species; IUCN: Gland, Switzerland, 2018. [Google Scholar] [CrossRef]
  29. Diesmos, A.; Buskirk, J.; Schoppe, S.; Diesmos, M.L.; Sy, E.; Brown, R. Siebenrockiella leytensis (Taylor 1920)—Palawan Forest Turtle, Philippine Forest Turtle. In Conservation Biology of Freshwater Turtles and Tortoises; IUCN: Gland, Switzerland, 2012. [Google Scholar]
  30. Sy, E.Y. The Palawan forest turtle. Traffic Bull. 2013, 25, 9. [Google Scholar]
  31. Schoppe, S.; Shepherd, C.R. The Palawan Forest Turtle: Under threat from international trade. TRAFFIC Bull. 2013, 25, 9–11. [Google Scholar]
  32. Schoppe, S.; Shepherd, C.R.; Beastall, C. The Palawan Forest Turtle: The story of a rare turtle found, then feared lost, only to be rediscovered—Now faces extinction. Tortoise 2013, 1, 108–117. [Google Scholar]
  33. CITES. Appendices I, II and III; IUCN: Gland, Switzerland, 2021. [Google Scholar]
  34. Boussekey, M. An integrated approach to conservation of the Philippine or Red-vented cockatoo: Cacatua haematuropygia. Int. Zoo Yearb. 2000, 37, 137–146. [Google Scholar] [CrossRef]
  35. Birdlife International. Cacatua Haematuropygia; The IUCN Red List of Threatened Species; IUCN: Gland, Switzerland, 2017. [Google Scholar] [CrossRef]
  36. Dickinson, E.C.; Kennedy, R.S.; Parkes, K.C.; McGowan, E. The Birds of the Philippines: An Annotated Checklist; British Ornithologists’ Union: Tring, UK, 1991; Volume 12. [Google Scholar]
  37. Schoppe, S.; Katsis, L.K.D.; Alvarado, D.; Acosta-Lagrada, L. Chapter 7—Philippine pangolin Manis culionensis (de Elera, 1915). In Pangolins: Science, Society and Conservation; Challender, D.W.S., Nash, J.C., Waterman, C., Eds.; Elsevier: Amsterdam, The Netherlands, 2020. [Google Scholar]
  38. Heinrich, S.; Wittmann, T.A.; Prowse, T.A.; Ross, J.V.; Delean, S.; Shepherd, C.R.; Cassey, P. Where did all the pangolins go? International CITES trade in pangolin species. Glob. Ecol. Conserv. 2016, 8, 241–253. [Google Scholar] [CrossRef] [Green Version]
  39. Aceto, M. The use of ICP-MS in food traceability. In Advances in Food Traceability Techniques and Technologies; Elsevier: Amsterdam, The Netherlands, 2016; pp. 137–164. [Google Scholar]
  40. Salvo, A.; Cicero, N.; Vadalà, R.; Mottese, A.F.; Bua, D.; Mallamace, D.; Giannetto, C.; Dugo, G. Toxic and essential metals determination in commercial seafood: Paracentrotus lividus by ICP-MS. Nat. Prod. Res. 2016, 30, 657–664. [Google Scholar] [CrossRef]
  41. Bonizzoni, L.; Galli, A.; Gondola, M.; Martini, M. Comparison between XRF, TXRF, and PXRF analyses for provenance classification of archaeological bricks. X-Ray Spectrom. 2013, 42, 262–267. [Google Scholar] [CrossRef]
  42. Natusch, D.J.; Carter, J.F.; Aust, P.W.; Van Tri, N.; Tinggi, U.; Riyanto, A.; Lyons, J.A. Serpent’s source: Determining the source and geographic origin of traded python skins using isotopic and elemental markers. Biol. Conserv. 2017, 209, 406–414. [Google Scholar] [CrossRef]
  43. Croudace, I.W.; Rindby, A.; Rothwell, R.G. ITRAX: Description and evaluation of a new multi-function X-ray core scanner. In New Techniques in Sediment Core Analysis; Rothwell, R.G., Ed.; Geological Society, London: London, UK, 2006; pp. 51–63. [Google Scholar]
  44. McGladdery, C.; Weindorf, D.C.; Chakraborty, S.; Li, B.; Paulette, L.; Podar, D.; Pearson, D.; Kusi, N.Y.O.; Duda, B. Elemental assessment of vegetation via portable X-ray fluorescence (PXRF) spectrometry. J. Environ. Manag. 2018, 210, 210–225. [Google Scholar] [CrossRef] [PubMed]
  45. Weindorf, D.C.; Bakr, N.; Zhu, Y. Advances in portable X-ray fluorescence (PXRF) for environmental, pedological, and agronomic applications. In Advances in Agronomy; Elsevier: Amsterdam, The Netherlands, 2014; Volume 128, pp. 1–45. [Google Scholar]
  46. Ribeiro, B.T.; Silva, S.H.G.; Silva, E.A.; Guilherme, L.R.G. Portable X-ray fluorescence (pXRF) applications in tropical Soil Science. Ciência e Agrotecnologia 2017, 41, 245–254. [Google Scholar] [CrossRef] [Green Version]
  47. Biagetti, S.; Alcaina-Mateos, J.; Ruiz-Giralt, A.; Lancelotti, C.; Groenewald, P.; Ibanez-Insa, J.; Gur-Arie, S.; Morton, F.; Merlo, S. Identifying anthropogenic features at Seoke (Botswana) using pXRF: Expanding the record of southern African Stone Walled Sites. PLoS ONE 2021, 16, e0250776. [Google Scholar] [CrossRef]
  48. Hobson, K.A.; Wassenaar, L.I. Tracking Animal Migration with Stable Isotopes; Academic Press: Cambridge, MA, USA, 2008. [Google Scholar]
  49. Brandis, K.J.; Mazumder, D.; Gadd, P.; Ji, B.; Kingsford, R.T.; Ramp, D. Using feathers to map continental-scale movements of waterbirds and wetland importance. Conserv. Lett. 2021, 14, e12798. [Google Scholar] [CrossRef]
  50. Hobson, K.A. Using stable isotopes to trace long-distance dispersal in birds and other taxa. Divers. Distrib. 2005, 11, 157–164. [Google Scholar] [CrossRef]
  51. Hobson, K.A. Tracing origins and migration of wildlife using stable isotopes: A review. Oecologia 1999, 120, 314–326. [Google Scholar] [CrossRef] [PubMed]
  52. Sung, Y.H.; Liew, J.H.; Chan, H.K.; Lee, W.H.; Wong, B.H.F.; Dingle, C.; Karraker, N.E.; Spencer, R.J.; Fong, J.J. 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. [Google Scholar] [CrossRef]
  53. Andersson, A.A.; Gibson, L.; Baker, D.M.; Cybulski, J.D.; Wang, S.; Leung, B.; Chu, L.M.; Dingle, C. Stable isotope analysis as a tool to detect illegal trade in critically endangered cockatoos. Anim. Conserv. 2021, 24, 1021–1031. [Google Scholar] [CrossRef]
  54. Brandis, K.J.; Meagher, P.J.; Tong, L.; Shaw, M.; Mazumder, D.; Gadd, P.; Ramp, D. Novel detection of provenance in the illegal wildlife trade using elemental data. Sci. Rep. 2018, 8, 15380. [Google Scholar] [CrossRef] [Green Version]
  55. Dempson, J.; Power, M. Use of stable isotopes to distinguish farmed from wild Atlantic salmon, Salmo salar. Ecol. Freshw. Fish 2004, 13, 176–184. [Google Scholar] [CrossRef]
  56. Hammershøj, M.; Pertoldi, C.; Asferg, T.; Møller, T.B.; Kristensen, N.B. Danish free-ranging mink populations consist mainly of farm animals: Evidence from microsatellite and stable isotope analyses. J. Nat. Conserv. 2005, 13, 267–274. [Google Scholar] [CrossRef]
  57. Dittrich, C.; Struck, U.; Rödel, M.O. Stable isotope analyses—A method to distinguish intensively farmed from wild frogs. Ecol. Evol. 2017, 7, 2525–2534. [Google Scholar] [CrossRef]
  58. Ziegler, S.; Streit, B.; Jacob, D.E. Assigning elephant ivory with stable isotopes. In Isotopic Landscapes in Bioarchaeology; Springer: Berlin/Heidelberg, Germany, 2016; pp. 213–220. [Google Scholar]
  59. DeNiro, M.J.; Epstein, S. Influence of diet on the distribution of carbon isotopes in animals. Geochim. Cosmochim. Acta 1978, 42, 495–506. [Google Scholar] [CrossRef]
  60. DeNiro, M.J.; Epstein, S. Influence of diet on the distribution of nitrogen isotopes in animals. Geochim. Cosmochim. Acta 1981, 45, 341–351. [Google Scholar] [CrossRef]
  61. Mumby, J.A.; Larocque, S.M.; Johnson, T.B.; Stewart, T.J.; Fitzsimons, J.D.; Weidel, B.C.; Walsh, M.G.; Lantry, J.R.; Yuille, M.J.; Fisk, A.T. Diet and trophic niche space and overlap of Lake Ontario salmonid species using stable isotopes and stomach contents. J. Great Lakes Res. 2018, 44, 1383–1392. [Google Scholar] [CrossRef]
  62. Schoppe, S.; Acosta, D. Philippine Freshwater Turtle Conservation Program (Pftcp) 2015 Annual Report; Katala Foundation: Puerto Princesa City, Philippines, 2016. [Google Scholar]
  63. Keller, J.M.; Pugh, R.S.; Becher, P.R. Biological and Environmental Monitoring and Archival of Sea Turtle Tissues (BEMAST): Rationale, Protocols, and Initial Collections of Banked Sea Turtle Tissues; U.S. Department of Commerce: Washington, DC, USA, 2014.
  64. Hua, L.; Gong, S.; Wang, F.; Li, W.; Ge, Y.; Li, X.; Hou, F. Captive breeding of pangolins: Current status, problems and future prospects. ZooKeys 2015, 507, 99–114. [Google Scholar]
  65. Paritte, J.M.; Kelly, J.F. Effect of cleaning regime on stable-isotope ratios of feathers in Japanese Quail (Coturnix japonica). The Auk 2009, 126, 165–174. [Google Scholar] [CrossRef]
  66. Dunnington, D. xrftools: XRF Tools for R. R package version 0.0.1.9000. 2021. Available online: https://github.com/paleolimbot/xrftools (accessed on 12 November 2021).
  67. Chen, T.; He, T.; Benesty, M.; Khotilovich, V.; Tang, Y.; Cho, H. Xgboost: Extreme Gradient Boosting. R package Version 1.3.2.1. 2021. Available online: https://rdrr.io/cran/xgboost/ (accessed on 12 November 2021).
  68. Chen, T.; Guestrin, C. Xgboost: A scalable tree boosting system. In Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, San Francisco, CA, USA, 13–17 August 2016; pp. 785–794. [Google Scholar]
  69. Lenth, R. Emmeans: Estimated Marginal Means, Aka Least-Squares Means; 2020. Available online: https://cran.r-project.org/web/packages/emmeans/emmeans.pdf (accessed on 12 November 2021).
  70. Gherase, M.R.; Mader, J.E.; Fleming, D.E. The radiation dose from a proposed measurement of arsenic and selenium in human skin. Phys. Med. Biol. 2010, 55, 5499–5514. [Google Scholar] [CrossRef]
  71. Specht, A.J.; Parish, C.N.; Wallens, E.K.; Watson, R.T.; Nie, L.H.; Weisskopf, M.G. Feasibility of a portable X-ray fluorescence device for bone lead measurements of condor bones. Sci. Total Environ. 2018, 615, 398–403. [Google Scholar] [CrossRef]
  72. Irshad, N.; Mahmood, T.; Hussain, R.; Nadeem, M.S. Distribution, abundance and diet of the Indian pangolin (Manis crassicaudata). Anim. Biol. 2015, 65, 57–71. [Google Scholar] [CrossRef]
  73. Sopsop, L.B. Floral survey in the coastal forest of Rasa wildlife sanctuary, Narra, Palawan, Philippines. J. Environ. Sci. Manag. 2011, 14, 71–76. [Google Scholar]
  74. Cherel, Y.; Hobson, K.A.; Weimerskirch, H. Using stable-isotope analysis of feathers to distinguish moutling and breeding origins of seabirds. Oecologia 2000, 122, 155–162. [Google Scholar] [CrossRef] [PubMed]
  75. Buddhachat, K.; Thitaram, C.; Brown, J.L.; Klinhom, S.; Bansiddhi, P.; Penchart, K.; Ouitavon, K.; Sriaksorn, K.; Pa-in, C.; Kanchanasaka, B.; et al. Use of handheld X-ray fluorescence as a non-invasive method to distinguish between Asian and African elephant tusks. Sci. Rep. 2016, 6, 24845. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Diesmos, A.C.; Gee, G.V.; Diesmos, M.; Brown, R.M.; Widmann, P.J.; Dimalibot, J.C. Rediscovery of the Philippine forest turtle, Heosemys leytensis (Chelonia; Bataguridae), from Palawan Island, Philippines. Asiat. Herpetol. Res. 2004, 10, 22–27. [Google Scholar]
  77. Pietersen, D.W.; Symes, C.T.; Woodborne, S.; McKechnie, A.E.; Jansen, R. Diet and prey selectivity of the specialist myrmecophage, Temminck’s ground pangolin. J. Zool. 2016, 298, 198–208. [Google Scholar] [CrossRef] [Green Version]
  78. Chesson, L.A.; Barnette, J.E.; Bowen, G.J.; Brooks, J.R.; Casale, J.F.; Cerling, T.E.; Cook, C.S.; Douthitt, C.B.; Howa, J.D.; Hurley, J.M.; et al. Applying the principles of isotope analysis in plant and animal ecology to forensic science in the Americas. Oecologia 2018, 187, 1077–1094. [Google Scholar] [CrossRef] [PubMed]
  79. Alacs, E.A.; Georges, A.; FitzSimmons, N.N.; Robertson, J. DNA detective: A review of molecular approaches to wildlife forensics. Forensic Sci. Med. Pathol. 2010, 6, 180–194. [Google Scholar] [CrossRef]
  80. Zhang, H.; Miller, M.P.; Yang, F.; Chan, H.K.; Gaubert, P.; Ades, G.; Fischer, G.A. Molecular tracing of confiscated pangolin scales for conservation and illegal trade monitoring in Southeast Asia. Glob. Ecol. Conserv. 2015, 4, 414–422. [Google Scholar] [CrossRef] [Green Version]
  81. Mwale, M.; Dalton, D.L.; Jansen, R.; De Bruyn, M.; Pietersen, D.; Mokgokong, P.S.; Kotze, A. Forensic application of DNA barcoding for identification of illegally traded African pangolin scales. Genome 2017, 60, 272–284. [Google Scholar] [CrossRef]
  82. van Uhm, D. Wildlife and laundering: Interaction between the under and upper world. In Green Crimes and Dirty Money, 1st ed.; Spapens, T., White, R., Huisman, W., van Uhm, D., Eds.; Routledge: London, UK, 2018. [Google Scholar]
  83. Van Song, N. Wildlife Trading in Vietnam; Situation, Causes and Solutions. J. Environ. Dev. 2008, 17, 145–165. [Google Scholar] [CrossRef]
  84. Thomas-Walters, L.; Veríssimo, D.; Gadsby, E.; Roberts, D.; Smith, R.J. Taking a more nuanced look at behavior change for demand reduction in the illegal wildlife trade. Conserv. Sci. Pract. 2020, 2, e248. [Google Scholar] [CrossRef]
  85. Kurland, J.; Pires, S.F.; Assessing, U.S. Wildlife Trafficking Patterns: How Criminology and Conservation Science Can Guide Strategies to Reduce the Illegal Wildlife Trade. Deviant Behav. 2017, 38, 375–391. [Google Scholar] [CrossRef]
Animals 13 02165 g001 550
Figure 1. Palawan forest turtles, Philippine cockatoos, and Philippine pangolins were sampled in the Province of Palawan (yellow) in the Philippines (green), situated in Southeast Asia (inset). Due to the ongoing illegal trade of these three species, fine-scale sampling locations are not provided.
Figure 1. Palawan forest turtles, Philippine cockatoos, and Philippine pangolins were sampled in the Province of Palawan (yellow) in the Philippines (green), situated in Southeast Asia (inset). Due to the ongoing illegal trade of these three species, fine-scale sampling locations are not provided.
Animals 13 02165 g001
Animals 13 02165 g002 550
Figure 2. Differences in beam spectra counts per second per element in Palawan forest turtle scute samples collected across the island of Palawan at three locations: KFI (orange, A) and two wild sampling locations (green, B, C).
Figure 2. Differences in beam spectra counts per second per element in Palawan forest turtle scute samples collected across the island of Palawan at three locations: KFI (orange, A) and two wild sampling locations (green, B, C).
Animals 13 02165 g002
Animals 13 02165 g003 550
Figure 3. Differences in beam spectra counts per second per element in Philippine cockatoo feather samples collected across the island of Palawan at three locations: KFI (orange, KFI) and two wild sampling locations (green, D,E).
Figure 3. Differences in beam spectra counts per second per element in Philippine cockatoo feather samples collected across the island of Palawan at three locations: KFI (orange, KFI) and two wild sampling locations (green, D,E).
Animals 13 02165 g003
Animals 13 02165 g004 550
Figure 4. Differences in beam spectra counts per second per element in wild Philippine pangolin scute and nail samples collected across two locations on the island of Palawan (F, G).
Figure 4. Differences in beam spectra counts per second per element in wild Philippine pangolin scute and nail samples collected across two locations on the island of Palawan (F, G).
Animals 13 02165 g004
Animals 13 02165 g005 550
Figure 5. Stable isotope δ13C and δ15N values for Philippine cockatoo, Philippine pangolin, and Philippine forest turtle with site and wild/captive group status identified.
Figure 5. Stable isotope δ13C and δ15N values for Philippine cockatoo, Philippine pangolin, and Philippine forest turtle with site and wild/captive group status identified.
Animals 13 02165 g005
Table
Table 1. Sample details of the Palawan forest turtle, Philippine cockatoo, and Philippine pangolin collected across Palawan Island, including captive animals cared for by Katala Foundation Incorporated. True locations are not provided in the interest of protecting these highly traded species.
Table 1. Sample details of the Palawan forest turtle, Philippine cockatoo, and Philippine pangolin collected across Palawan Island, including captive animals cared for by Katala Foundation Incorporated. True locations are not provided in the interest of protecting these highly traded species.
Common NameSpecies NameSampleProvenanceSiten
Palawan forest turtleSiebenrockiella leytensisScuteCaptiveKFI24
WildB18
C3
Philippine cockatooCacatua haematuropygiaFeatherCaptiveKFI20
WildD24
E21
Philippine pangolinManis culionensisClawWildF10
ScaleWildG6
Table
Table 2. Predictive model results differed in accuracy amongst the three study species, Palawan forest turtle, Philippine cockatoo, and Philippine pangolin, based on the datasets used (portable X-ray fluorescence (pXRF) or stable isotope analysis (SIA)) and based on the response variable (captive or wild status, geographic origin).
Table 2. Predictive model results differed in accuracy amongst the three study species, Palawan forest turtle, Philippine cockatoo, and Philippine pangolin, based on the datasets used (portable X-ray fluorescence (pXRF) or stable isotope analysis (SIA)) and based on the response variable (captive or wild status, geographic origin).
SpeciesResponse
Variable		pXRF Model
% Accuracy ± SD 		SIA Model
% Accuracy ± SD
Palawan forest turtleCaptive or wild88 ± 12%66 ± 16%
Geographic origin94 ± 5%73 ± 14%
Philippine cockatooCaptive or wild78 ± 2%28 ± 1%
Geographic origin62 ± 9%27 ± 2%
Philippine pangolinGroup93 ± 14%41 ± 25%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Brandis, K.J.; Meagher, P.; Schoppe, S.; Zawada, K.; Widmann, I.; Widmann, P.; Dolorosa, R.G.; Francis, R. Determining the Provenance of Traded Wildlife in the Philippines. Animals 2023, 13, 2165. https://doi.org/10.3390/ani13132165

AMA Style

Brandis KJ, Meagher P, Schoppe S, Zawada K, Widmann I, Widmann P, Dolorosa RG, Francis R. Determining the Provenance of Traded Wildlife in the Philippines. Animals. 2023; 13(13):2165. https://doi.org/10.3390/ani13132165

Chicago/Turabian Style

Brandis, Kate J., Phoebe Meagher, Sabine Schoppe, Kyle Zawada, Indira Widmann, Peter Widmann, Roger G. Dolorosa, and Roxane Francis. 2023. "Determining the Provenance of Traded Wildlife in the Philippines" Animals 13, no. 13: 2165. https://doi.org/10.3390/ani13132165

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

Article Metrics

Back to TopTop