Present and Potential Future Distributions of Asian Horseshoe Crabs Determine Areas for Conservation
Stine Vestbo
Matthias Obst
Francisco J. Quevedo Fernandez
Itsara Intanai
Peter Funch- 1Genetics, Ecology, and Evolution, Department of Bioscience, Aarhus University, Aarhus, Denmark
- 2Department of Biological and Environmental Sciences, University of Gothenburg, Gothenburg, Sweden
- 3School of Computer Science and Informatics, Cardiff University, Cardiff, United Kingdom
- 4Department of Biology, Faculty of Science and Technology, Prince of Songkla University, Pattani, Thailand
Conservation of horseshoe crabs has recently received increasing attention as several populations are in decline. However, scarce information on their distributions in Southeast Asia is impairing conservation efforts. In this study, we sought to improve our understanding of the geographical range and distinct populations of the three Asian horseshoe crabs species in order to identify optimal conservation areas. We mapped the geographic range of Carcinoscorpius rotundicauda, Tachypleus gigas, and T. tridentatus using recent data from field work, literature, Global Biodiversity Information Facility (GBIF), and unpublished data from our scientific network. The data were correlated with 23 different environmental variables of potential ecological importance for horseshoe crabs using the openModeller webservices, including new tidal variables. Ecological niche models were generated using two algorithms, Maximum Entropy and support vector machine, for the three species under present conditions, and projected into a climate change scenario of 2050. The niches of the Asian horseshoe crabs were mostly determined by tidal regime, chlorophyll A concentrations, depth, distance to land, and sea surface temperature. According to our predictions, horseshoe crabs in Southeast Asia are not expected to experience any severe change in extent and distribution of suitable habitat in the future. In order to conserve Asian horseshoe crabs, we suggest establishing Marine Protected Areas at locations where distinct populations and several species occur, such as northern Vietnam, China, Borneo, and southern Japan.
Introduction
Coastal and intertidal areas are currently under threat globally due to a range of anthropogenic activities, including infrastructure development and coastal protection, as well as effects of climate change, such as rising sea levels pushing coastal areas closer to anthropogenic structures, resulting in “coastal squeeze” (Defeo et al., 2009). Coastal areas are characterized by high productivity and species richness (Ketchum, 1972; Ray, 1991), and humans depend strongly on a variety of species that are living and reproducing in the coastal zone, such as fishes (Allison et al., 2009; Bell et al., 2009), sea cucumbers (Purcell et al., 2013), and horseshoe crabs (Chatterji, 1994), thus the conservation of such species is of great importance.
During the past decade, horseshoe crab conservation has received increasing attention (Botton et al., 2015). Horseshoe crabs are marine chelicerates of the order Xiphosura, containing only four extant species. Three species occur in the coastal waters of India, Southeast Asia, China, and Japan [Carcinoscorpius rotundicauda (Latreille, 1802), Tachypleus gigas (Müller, 1785), and T. tridentatus (Leach, 1819)]. One species [Limulus polyphemus (Linnaeus, 1758)] occurs in the coastal waters of the eastern continental shelf of North America and in the Gulf of Mexico (Sekiguchi and Shuster, 2009). Populations of horseshoe crabs world-wide are currently in decline (e.g., Cartwright-Taylor et al., 2011; Kwan et al., 2016; Smith et al., 2017), and the IUCN Red List of Threatened Species now describes L. polyphemus as vulnerable (Smith et al., 2016). Although the Asian horseshoe crab species are thought to be in greater decline than L. polyphemus, these species are currently listed as data deficient (World Conservation Monitoring Centre, 1996a,b,c). While the geographic range of L. polyphemus is known (Sekiguchi, 1988; Anderson and Shuster, 2003; Sekiguchi and Shuster, 2009; Faurby et al., 2011), the precise distribution of the three Asian species remains less clear (World Conservation Monitoring Centre, 1996a,b,c), thus impairing conservation efforts.
The global decline in horseshoe crab populations is mainly due to anthropogenic activities. For L. polyphemus the decline has primarily been due to overharvesting and loss of breeding habitat (Smith et al., 2017), as this species breeds on beaches with slight slopes; a habitat commonly used for real estate construction and development (Nordstrom, 2004). Coastal armoring as a response to erosion is an additional factor reducing available breeding habitat (Jackson et al., 2015). Wild caught horseshoe crabs are used commercially in medicine, where substances from their hemolymph (Carcinoscorpius, Limulus, and Tachypleus Amoebocyte Lysate, CAL, LAL, and TAL) are used to test if drugs, blood products, and pharmaceutical devices are free from bacterial contamination (Levin and Bang, 1968; Rao and Bhagirathi, 1989; Levin et al., 2003). L. polyphemus is released after collection of LAL, but the two Tachypleus species are commonly used for chitin production after hemolymph collection, thus harvesting of these species results in 100% mortality (Gauvry, 2015). However, a synthetic, commercially available version of the CAL test has been developed (Ding et al., 1995; Lonza, 2016). The three Asian horseshoe crab species are fished for human consumption (Botton, 2001), and were previously imported to North America to be used as bait (Smith et al., 2016), although this import is now prohibited (IUCN, 2013). As for L. polyphemus, the expansion of urban infrastructure into coastal areas, as well as coastal armoring present additional anthropogenic habitat disturbances for the Asian horseshoe crabs, thus in some areas, such as Japan and Peninsular Malaysia, horseshoe crab breeding habitats are now almost completely lost due to coastal infrastructure development (Botton, 2001; Nelson et al., 2016). In response to the declining horseshoe crabs populations, Marine Protected Areas (MPAs) have been established, especially in the USA (Carl N. Shuster Jr. Horseshoe Crab Reserve, Delaware Bay) and Japan (Saikai National Park in Nagasaki). However, in most of Southeast Asia MPAs protecting critical horseshoe crab habitats are still scarce.
Climate change might have additional detrimental effects on the distribution and population sizes of horseshoe crabs. Rising sea levels could reduce suitable horseshoe crab breeding habitats, since the retreating shorelines would bring the horseshoe crabs closer to the human infrastructure near the coasts. This is the case for Delaware Bay, New Jersey, USA, where sea levels have risen at a high rate during the twentieth century and hence decreased horseshoe crab breeding habitat (Loveland and Botton, 2015). Climate change would also affect ecological factors that influence the developmental success of horseshoe crab eggs and larvae, such as salinity and temperature (Jegla and Costlow, 1982; Laughlin, 1983; Ehlinger and Tankersley, 2004; Zaleha et al., 2011), oxygen levels (Palumbi and Johnson, 1982; Funch et al., 2016), beach geochemistry, wave energy, and erosion (Botton et al., 1988; Penn and Brockmann, 1994; Jackson et al., 2008). One of the species, C. rotundicauda, is found in habitats characterized by mangroves, which are also vulnerable to the effects of climate change, primarily rising sea levels (Gilman et al., 2008).
Ecological niche modeling (ENM, also known as species distribution modeling) allows for predictions of suitable habitats for a given species by identifying the environmental factors driving the species' distribution, such as temperature, salinity, or depth (Peterson et al., 2011). Furthermore, ENM can be used to project the distribution of suitable habitats of a species into future scenarios (e.g., future climate change scenarios), allowing forecasts of possible species range shifts. For example, the method has been used to predict future distributions and range shifts for Mediterranean fishes (Albouy et al., 2013) as well as changes in suitable habitat for commercially important fish and invertebrates along the Atlantic coast of North America (Kleisner et al., 2017). ENM can be regarded as an important tool for conserving threatened species, and has proven useful for the assessment and planning of protected areas (Kremen et al., 2008; Stirling et al., 2016). Thus, the scattered information on Asian horseshoe crab distributions and their unknown conservation status would benefit from gathering recent occurrences as well as a modeling of suitable habitats in the region. The objectives of this study were to (1) identify important environmental factors that influence the distributions of the three Asian species of horseshoe crabs, (2) update the information on current species distributions, and (3) identify possible high priority regions for their conservation, where climate change has minimal effect and where species ranges overlap.
Materials and Methods
Species Data and Study Area
The geographic study area included all coastal waters between 10.0° S to 35.0° N and 77.0° W to 135.0° E, thus including the continental shelves of the Bay of Bengal, India and Indonesia and between Indonesia and Japan. This area, the geographical mask, was chosen based on the existing knowledge of horseshoe crab distribution limits in Asia (Sekiguchi, 1988), and used for geographical filtering of occurrence records, collecting background points, and building and projecting the models. We assembled occurrence records (presence-only) for all three species in this region from our own observations, collaborators, scientific networks as well as through publishing a scratchpad site at http://horseshoecrabs.myspecies.info/ (Table 1, Supplementary Table 1). For many species, numerous distribution records exist in the literature, and we manually geo-referenced additional occurrence data from these sources (Supplementary Table 2). Subsequently, we added data from the Global Biodiversity Information Facility (http://gbif.org) in September, 2015 (Supplementary Table 3). A taxonomic data refinement workflow (Mathew et al., 2014) was used to check the synonyms, download, visualize, filter, and integrate occurrence records for all species. We inspected all records and excluded those of dubious identity (e.g., on land, in deep water, or without clear reference of origin), as well as those older than 1995 and those outside the geographical mask. All new data (i.e., excl. GBIF records) were submitted to the Ocean Biogeographic Information System (http://www.iobis.org/) and are available under the link https://doi.org/10.14284/293.
Table 1. Occurrence data used in the analysis. Environmentally unique points (EUPs) were calculated from the compiled observations (records) using the respective geographical mask with a resolution of 300 arc seconds and 23 environmental variables.
Environmental Data
Based on available information on general ecology of horseshoe crabs (Bonaventura et al., 1982; Sekiguchi, 1988; Shuster et al., 2003), we included 23 environmental variables in the analysis, relating to sea surface temperature, chlorophyll A concentration, salinity, oxygen and nutrient concentrations, distance to land, water depth, pH, and tidal regime (Table 2). Chlorophyll A variables were included since chlorophyll levels acts as a proxy for phytoplankton primary production, which in turn reflects food availability for filter-feeding bivalves, which is one of the most important food items for horseshoe crabs (Chatterji, 1994). We used 13 marine layers from Bio-Oracle (http://www.bio-oracle.ugent.be/) with a resolution of five arc-minutes (Table 2; Tyberghein et al., 2012). These data layers are generated from monthly satellite data (Aqua-MODIS and SeaWiFS; https://oceancolor.gsfc.nasa.gov) as well as in situ measured oceanographic data from the World Ocean Database, 2009 (Boyer et al., 2009), and overlap with the species occurrence data in their temporal origin. In addition, we used five marine layers for both present and 2050 climate scenarios from AquaMaps (http://www.aquamaps.org/download/main.php) with a resolution of 30 arc-minutes (Table 2; Kaschner et al., 2008). Present-day (i.e., year 2010) datasets from AquaMaps were built from long-term averages of temporally varying environmental variables (Ready et al., 2010), while future data sets for 2050 were derived from the IPSL-CM4 A2 climate change scenario (IPCC, 2007).
Table 2. Environmental variables investigated for ecological niche modeling of the three Asian horseshoe crabs.
We additionally generated five present-day (i.e., year 2010) global tidal variables with the Finite Element Solution oceanographic model (FES2012), provided by Noveltis, Legos and CLS Space Oceanography Division and distributed by AVISO+ (http://www.aviso.altimetry.fr/) (Table 2, more detailed information on how the global tidal variables were created can be found in Supplementary Data Sheet 1). FES2012 is a fully revised version of the original global hydrodynamic tide solution model (Lyard et al., 2006) with improvements regarding longer altimeter time series, more accurate ocean bathymetry, major non-linear tides, and more accurate tidal currents. It includes overall 32 tidal constituents distributed on 1/16° grids (amplitude and phase), corresponding to 3.75 arc-minutes (Carrère et al., 2012). The tidal variable layers can be accessed through a webservice (https://www.biodiversitycatalogue.org/services/37), and is hence available for future research on distribution modeling of marine intertidal organisms. In addition, we deposited the raster files together with the algorithm for calling the FES program and creating the ASCII grids at the Swedish Environmental and Climate Data Repository (www.ecds.se) under the identifier 87196c86-8e36-4908-8f16-424e61aa313c.
Analysis of Environmental Variables
Environmental values were obtained using version 2 of the BioClim workflow (http://purl.ox.ac.uk/workflow/myexp-3725.2) available at the BioVeL portal, https://portal.biovel.eu (Hardisty et al., 2016). The workflow was used to retrieve environmentally unique points (EUPs) from the species occurrence files for a set of 23 environmental data sets (Table 2). Since all points within a raster cell share the same values for the environmental variables, EUPs are obtained by filtrating the occurrence points, so that each raster cell contains no more than a single point (Nix, 1986). The workflow returned environmental values from the specified layers that matched with the species occurrence records, and it was executed in batch mode (called data sweep function) to repeat all calculations with the same parameter settings for all species.
We analyzed differences in environmental tolerances between all three horseshoe crab species by applying a Welch's test for unequal variances to the retrieved environmental values (Supplementary Table 4). This test was chosen because the data did not conform to parametric test assumptions. In addition, pairwise comparisons of environmental tolerances between species were performed using the Games-Howell Post-Hoc test (Supplementary Table 5). Both statistical tests were performed using SPSS version 22 (IBM Corporation, 2013). We subsequently performed a correlation analysis (data not shown) as well as a principal component analysis (PCA) on the environmental values using the R statistical environment 3.0.2 (R Core Team, 2013). The analyses were performed in order to identify ecologically important variables that explain the variation in the data set, and which can be used as predictor variables (Metzger et al., 2005; Porfirio et al., 2014). The PCA was estimated as a 23-dimensional hypervolume and used to compare the environmental space occupied by the three species. We used all non-correlated variables that were represented (with >5%) in the first two components of the PCA to build the niche models, except for Mean SST, which was included in the future projection (Supplementary Tables 6–8).
Ecological Niche Modeling
We used version 20 of the ecological niche modeling (ENM) workflow (http://purl.ox.ac.uk/workflow/myexp-3355.20) in batch mode (called data sweep) to describe and compare the geographical space occupied by the three species and to estimate the distribution of suitable habitat with favorable biotic, environmental, and geographical conditions. Demographic or dispersal properties that may also be used in species distribution modeling were not considered (Peterson et al., 2011; Reiss et al., 2014). For more information on the ENM workflows see Holl et al. (2013), Leidenberger et al. (2015a,b), and De Giovanni et al. (2015). We executed parallel analyses with two ENM algorithms by means of the openModeller webservice suite (de Souza Muñoz et al., 2011). These are i) Maximum Entropy v. 1.0 (Phillips et al., 2006; Phillips and Dudík, 2008) and ii) support vector machine v. 0.5 (Schölkopf et al., 2001). These two methods were chosen because Maximum Entropy (Maxent) is a presence-background algorithm that has been shown to perform with high predictive accuracy even if there are moderate errors in the geo-referencing (Elith et al., 2006; Graham et al., 2008), and support vector machine (SVM) is a machine-learning algorithm that is not analyzing characteristics of statistical distributions and hence does not require independent observation data, and thereby overcomes potential problems with spatial autocorrelation in the species occurrence data (Drake et al., 2006).
Models were created using each species' maximum distribution range within the mask and a set of predictor variables identified in the PCA. Variables selected for present day projections (high resolution variables) and future trend analysis (low resolution variables) are listed in Table 2. Models were created based on EUPs (Table 1) with the following specifications. Maxent models were set to run with 10.000 background points (including input points) drawn from the mask. Feature selection was automated, allowing the algorithm to combine feature types when fitting a model, and perform 500 iterations. Tolerance for detecting model convergence was set to 0.00001, while sample threshold was set to 80 (product), 10 (quadratic), and 15 (hinge). SVM models were set to execute the C-SVC algorithm with radial basis kernels, gamma values 1/k (where k is the number of layers), and a cost value of 1. All models were set to produce a probabilistic, instead of binary output. Predictive models may encounter some problems with extrapolations (Zurell et al., 2012), i.e., when the environmental hyperspace of the scenario is not represented by true observations. These were minimized by using only few essential variables in the future predictions, thereby avoiding complexity and interaction between too many factors, as well as through using a consensus approach (Table 2; Pearson et al., 2006).
For each species, we ran both algorithms across two sets of environmental layers (specified above and in Table 2); one with high resolution (using Bio-Oracle and FES2012 layers) and one with low resolution (using AquaMaps layers). The high-resolution models were used for present day projections of suitable habitat, while the low-resolution models were used to deduct climate driven changes in suitable habitat. Models were tested using 10-fold cross-validation based on the area under the curve (AUC) value and omission error rate (false negative rate), and subsequently projected using present and (in case of AquaMaps layers) 2050 climate scenarios (Table 3). The results of the ENMs were visualized as maps showing the distribution of suitable habitat of a species in a particular geographic region at a particular time (de Souza Muñoz et al., 2011). The sweep function allowed automated batch processing of all species for a given algorithm and parameter set. Overall, we executed 12 niche models, i.e., two algorithms, three species, and two environmental datasets (Table 3).
Table 3. Results of model tests for all models of current distribution (based on Bio-Oracle and tidal variables) and future trend analysis (based on AquaMaps variables).
GIS Analysis
Raster maps created by the niche modeling algorithms were processed using the qGIS software package v. 2.6 Brighton (Quantum GIS Development Team, 2014). Maps of individual models are available in the Supplementary Material (Supplementary Images 1–6). For the consensus analysis presented in the maps, we transformed each individual model into a binary output using the lowest presence threshold (LPT) in order to include any signal of suitable habitat. Thereafter the maps of both models were merged into a single consensus model highlighting only areas where both algorithms agreed on the prediction of suitable habitat. Likewise, heat maps showing the predicted changes between 2050 and present-day distributions were created as consensus maps identifying only areas where both algorithms agreed on predicted loss or gain of suitable habitat. Finally, we compared high-resolution models (used for present day projections) with low-resolution models (for future change projections). We found that both model settings predict suitable habitat for similar biogeographic regions with few deviations mentioned in the discussion below.
Results
Environmental Factors Defining Horseshoe Crab Distributions
The PCA showed that C. rotundicauda and T. gigas largely share the same ecological niche, and that the niche of T. tridentatus to some extent overlaps with the niches of the other two species (Figure 1). Further, the PCA showed that the distributions of the three species are mainly driven by chlorophyll concentration, tidal regime, temperature, depth and distance to land (Figure 1, Table 2). The combination of the environmental tolerances represents the specific fundamental ecological niche for each of the three species. We found no difference between the three Asian species in any of the tidal variables (Figures 2A–C). All species experienced varied annual cycle deviation as well as cycle durations of varied lengths. The mean minimum sea surface temperature found for T. tridentatus was 17.0°C and hence 10 degrees lower than those for T. gigas and C. rotundicauda (Figure 2D). T. tridentatus experienced a larger mean temperature range of 12.7°C compared to T. gigas and C. rotundicauda, which experienced a mean range of 3.9 and 4.0°C, respectively (Figure 2E). The water oxygen levels followed the same pattern as sea surface temperatures, with a higher mean for T. tridentatus (5.0 mL/L) compared to T. gigas and C. rotundicauda (4.4 and 4.5 mL/L, respectively, Figure 2F). Neither the maximum chlorophyll A concentration nor the range of chlorophyll A concentration differed between the species (Figures 2G,H). We found a mean salinity of 33.5 PSU for T. tridentatus, which was slightly higher compared to T. gigas (31.8 PSU) and C. rotundicauda. (32.3 PSU, Figure 2I), although the narrow salinity range for C. rotundicauda was likely an artifact caused by inaccessible data from inland waters (see Discussion). The mean nitrate level was higher for T. tridentatus (1.9 μmol/L) than the other two species (0.6 μmol/L for both) (Figure 2J). Mean phosphate differed between all three species, with T. gigas experiencing the highest (0.28 μmol/L) and T. tridentatus the lowest (0.20 μmol/L) mean level (Figure 2K). T. tridentatus differed from the other two species in occupying a niche with higher mean pH value (Figure 2L).
Figure 1. Niche overlap analysis of the three Asian horseshoe crab species showing the principal component analysis plots of the environmental “hypervolume” for 23 predictor variables [CR, C. rotundicauda (
); TG, T. gigas (
); TT, T. tridentatus (
)]. Ellipses represent 68% of the hypervolume for each species; points represent presence of each species at environmentally unique locations. The positions of variable names have been adjusted to prevent overlapping.
Figure 2. Comparison of the three horseshoe crab species (CR, C. rotundicauda; TG, T. gigas; TT, T. tridentatus) showing the variation of 12 key environmental variables. Variable and unit are shown on the y-axis. Gray boxes indicate 50% of the sample points and are limited by the 1st (Bottom) and 3rd quartiles (Top). Black, horizontal lines within the gray boxes display the medians. Comparisons between C. rotundicauda and the two Tachypleus species for salinity were excluded because of absence of environmental data for inland waters. Asterisks show significant (<0.05) differences between species.
Species Distributions
The sampled distribution of C. rotundicauda includes the coasts of eastern India and Bangladesh, the coasts of Malaysia, Thailand, Cambodia, and southern Vietnam (Figure 3), but excludes the eastern coast of Vietnam, however, the distribution continues along the coast of northern Vietnam and southern China. C. rotundicauda is also present along the coasts of the Indonesian islands Sumatra and Java, as well as Borneo. The distribution of T. gigas is similar to that of C. rotundicauda, except that T. gigas seems to be absent from the coasts of Bangladesh, Cambodia, northern Vietnam, and China (Figure 4). T. tridentatus is distributed along the coasts of southern Japan, China, Taiwan, and northern Vietnam, in addition to the Philippines, and the islands of Borneo and Java (Figure 5).
Figure 3. Distribution maps for Carcinoscorpius rotundicauda with environmentally unique points, indicating the sampled distribution range for the species. Upper map (A) shows present-day projections of suitable habitat in red based on consensus of Maxent and SVM algorithms. Lower map (B) shows changes between the present-day and 2050 projections of suitable habitat based of consensus of Maxent and SVM algorithms, with green cells indicating loss of suitable habitat and red cells indicating gain of suitable habitat. Support values for the individual models are given in Table 3.
Figure 4. Distribution maps for Tachypleus gigas with environmentally unique points, indicating the sampled distribution range for the species. Upper map (A) shows present-day projections of suitable habitat in red based on consensus of Maxent and SVM algorithms. Lower map (B) show changes between the present-day and 2050 projections of suitable habitat based of consensus of Maxent and SVM algorithms, with green cells indicating loss of suitable habitat and red cells indicating gain of suitable habitat. Support values for the individual models are given in Table 3.
Figure 5. Distribution maps for Tachypleus tridentatus with environmentally unique points, indicating the sampled distribution range for the species. Upper map (A) shows present-day projections of suitable habitat in red based on consensus of Maxent and SVM algorithms. Lower map (B) show changes between the present-day and 2050 projections of suitable habitat based of consensus of Maxent and SVM algorithms, with green cells indicating loss of suitable habitat and red cells indicating gain of suitable habitat. Support values for the individual models are given in Table 3.
Distribution of Suitable Habitat
Both models, SVM and Maxent, demonstrated high predictive power, as seen by the high AUC values (all values >0.9) (Table 3). The Maxent model performed marginally better than SVM for all species except T. tridentatus. Omission errors were very low for both models (<3.5%), except for Maxent's future prediction for C. rotundicauda (6.2%).
Suitable habitats for C. rotundicauda are distributed along the coastline of mainland Southeast Asia, from the Bay of Bengal, including Sri Lanka, to southern China, except from the eastern coasts of Peninsular Malaysia and Vietnam (Figure 3A). Suitable habitats are also present in the Malacca Strait between Sumatra and Peninsular Malaysia, Southeast Sumatra, southern coast of Java, the coast of Kalimantan (Indonesian Borneo), and the coast of Sarawak (Malaysian Borneo). Additionally, small areas of suitable habitat are found in Sabah (Malaysian Borneo), West Papua, and in Manila Bay, Philippines. The distribution of suitable habitat differs from the sampled distribution, as it includes the coasts of Myanmar, southern Borneo, eastern Sumatra, and southern Java.
The distribution of suitable habitat for T. gigas is similar to that of C. rotundicauda, however, for T. gigas it does not go further east along the Southeast Asian mainland than to southeast Vietnam (Figure 4A), nor does it include the south coast of Java. Compared to the current known distribution, the distribution of suitable habitat includes the coast of Myanmar, Bangladesh, east coast of Sumatra, and southern Borneo.
The distribution of suitable habitat of T. tridentatus spans from northern Vietnam, along the coast of China, and into the Sea of Japan (Figure 5A). It also includes the coast of northwest Taiwan and Japan, including the Seto Inland Sea. Small areas of suitable habitat are present along the coasts of Southwest Myanmar, East Sumatra, South Java, Borneo, South Vietnam, Bay of Bangkok, and Manila Bay.
Changes in the Distribution of Suitable Habitat Under Climate Change Scenarios
We included three environmental variables, mean sea surface temperature, mean chlorophyll A concentration, and distance to land, for the modeling of the future distributions of suitable habitat (Table 2). Our models predicted that the habitat suitability for C. rotundicauda will increase in several coastal areas throughout Southeast Asia. More suitable habitats will arise in more eastern regions; the Philippines and along the coasts of several Indonesian islands: Sulawesi, Lombok, Sumbawa, and Flores (Figure 3B). Similarly, for T. gigas several areas along the coasts of Southeast Asia will become more suitable in the future (Figure 4B). The Bay of Bengal and North Sumatra will be less suitable. T. tridentatus will lose suitable habitat in the future, however, none of the lost areas are overlapping with the species' current occurrences (Figure 5B). In summary, none of the Asian horseshoe crabs seem to undergo a major range shift driven by prospective changes of the environmental factors tested here.
Discussion
Carcinoscorpius rotundicauda
We found that C. rotundicauda occurs along the coasts of Bangladesh, Cambodia, Vietnam, East Thai-Malay Peninsula, and China, in addition to the countries and coasts already surveyed by Sekiguchi (1988) (India, West Thai-Malay Peninsula, Singapore, and Indonesia). However, contrary to Sekiguchi (1988) we have not found evidence for its presence at Palawan, Philippines (Figure 3). According to the occurrence record, C. rotundicauda is absent along the coast of Myanmar, even though our models indicate suitable habitat in this area. Indeed, there is evidence that horseshoe crabs are present in this region, however, the available information points to L. polyphemus, and is likely a misidentification (Thapanand-Chaidee et al., 2010). We thus conclude that C. rotundicauda probably is present along the coast of Myanmar, and that missing records of horseshoe crabs are due to lack of studies in the area. The sampled distribution of C. rotundicauda shows a gap along the east coast of Vietnam, extending approximately 1,500 km in the high-resolution projections (based on Bio-Oracle variables, Figure 3) and approximately 400 km in the low-resolution projections (based on AquaMaps variables). The absence of the species in this region might be explained by the narrow continental shelf present along the coastline, which provides less horseshoe crab habitat, as well as water characterized by a high degree of mixing with deep water (Rojana-anawat et al., 2001). The distribution gap could also be related to food availability, as the east coast of central Vietnam has a lower macrozoobenthos biomass compared to the southern and northern coasts (Nguyen and Dao, 1995). However, the lack of mangroves along the east coast of central Vietnam may be the most plausible explanation for the distribution gap (Hong and San, 1993). C. rotundicauda is commonly known as the mangrove horseshoe crab, as it is often found in habitats characterized by mangroves and mudflats in contrast to the two Tachypleus species, which can be found at open sandy beaches (e.g., Cartwright-Taylor et al., 2011; Robert et al., 2014; Jawahir et al., 2017), thus a scarcity of mangroves would have a negative effect on the presence of C. rotundicauda.
When interpreting our model results, it is important to mention that species distribution models do have caveats, which may lead to false or uncertain projections of actual species distribution ranges (Jarnevich et al., 2015). We have tried to mitigate such caveats and uncertainties through an unbiased sampling of the biological data from independent sources, a careful selection of predictor variables, and through modeling with several environmental data sets and analytical algorithms. Nevertheless, the indication of a gap in the distribution of C. rotundicauda along the Vietnamese east coast should still be treated as a hypothesis that needs further validation with additional sampling and modeling.
The future distribution of suitable habitat of C. rotundicauda includes several coastal areas of Indonesian islands east of the sampled distribution. Since horseshoe crabs reside on the continental shelves (Sekiguchi and Shuster, 2009), it has been argued that the distributions of the Asian horseshoe crabs are delimited by the eastern edge of the continental Sunda Shelf (roughly corresponding to Wallace Line). Thus, even though areas of future suitable habitat are present east of the Sunda Shelf, we would not expect C. rotundicauda to be able to reside here, unless humans have mediated their dispersal. The models used in this study find an approximation of the fundamental niche of the horseshoe crabs, not the realized niche, thus geographical barriers are not necessarily detected.
We found that C. rotundicauda, as well as the other two Asian horseshoe crab species, inhabit coastal areas with varied tidal regimes. This is in contrast to L. polyphemus, as this species can be found in areas with mixed semi-diurnal tides (i.e., tidal cycles of approximately 12 h) (Rudloe, 1985; Barlow et al., 1986). The narrow range of salinity found for C. rotundicauda is most likely an artifact, as we know that the three Asian species occur in river estuaries, and C. rotundicauda has been found at salinities as low as 10 PSU (Chatterji, 1999; Mishra, 2009; Chen et al., 2015). Unfortunately, we could not include occurrences from inland rivers in Southeast Asia in our models, due to lack of environmental data from these areas. The narrow temperature range of C. rotundicauda reflects its tropical distribution, and can be related to its preference for temperatures between 20 and 30°C (Srijaya et al., 2014), which seems to be necessary for regular ecdysis (Lee and Morton, 2005). We find that C. rotundicauda is present in waters with lower pH levels compared to seawaters of normal pH level (around 8.2). This corresponds well with previous findings suggesting that pH levels found in horseshoe crab habitat are low due to the presence of estuaries (mixing of freshwater with pH levels of 7–7.5), photosynthesis by algae, as well as anoxic decomposition of organic material (Jawahir et al., 2017).
Tachypleus gigas
The distribution of T. gigas found in this study is consistent with that found by Sekiguchi (1988), with the exception that we also found occurrences of T. gigas along the east coast of the Thai-Malay Peninsula (Figure 4). Our high-resolution projections show limited suitable habitat at the southern extension range of the species (i.e., around Java), but the low-resolution projections clearly indicate suitable habitat along the northern coastline of Java.
As with C. rotundicauda, we suggest that T. gigas is present along the coast of Myanmar even though occurrence data are missing, due to the lack of research in the area. The presumed presence of T. gigas in this area needs to be confirmed with real observations in the future. Alternatively, there is the possibility that other factors not included in the models, such as biological interactions or substrate conditions, may actually refrain T. gigas form inhabiting this region. Our models predict additional future suitable habitats in more eastern regions of Indonesia, but again we argue that T. gigas is unable to disperse into these areas because of the deep waters east of the continental Sunda Shelf.
T. gigas is known to occur in areas with lower salinities than our results show (15–26 PSU) (Cartwright-Taylor et al., 2011). However, T. gigas is not as tolerant to low salinities as C. rotundicauda (Chatterji, 1999), and it has been shown that salinities lower than 20 PSU prevent its eggs from hatching (Zaleha et al., 2011). Apart from salinity, the ecological tolerances of T. gigas found in this study are similar to those of C. rotundicauda, which reflect their overlapping distributions.
Although we found overlapping niches and distributions for C. rotundicauda and T. gigas, the local habitat preferences are different between the species on a smaller spatial scale, i.e., C. rotundicauda is associated with mangroves and mudflats and T. gigas with sandy sediments at more exposed beaches, and adding an environmental factor such as sediment characteristics to our models might separate the niches of the two species (Jawahir et al., 2017).
Tachypleus tridentatus
In contrast to the T. tridentatus distribution reported by Sekiguchi (1988), we found no recent occurrence data on this species at the coasts of Sumatra, Sulawesi, South Vietnam, or Luzon, Philippines (Figure 5). However, we found several occurrences along the coast of South China. Although the coast of Myanmar is an area of suitable habitat for T. tridentatus, we would not expect it to occur here, as the area lies far west from any sampled occurrence sites. But according to our high-resolution model, the species might occur at the coast of Sumatra, as suggested by Sekiguchi (1988), even though we did not find any recent records of its presence here. There are no occurrences of T. tridentatus within the Yellow Sea, and we were unable to find historical evidence suggesting previous presence here, although the species is distributed at the coasts of the Korean Peninsula and Japan as well as in the East and South China Sea. Our high-resolution models indicated suitable habitat in the Yellow Sea, but this was not confirmed by the low-resolution models, leaving some degree of uncertainty in our projections for this area. The reason why T. tridentatus may be absent in the Yellow Sea could be related to a variety of factors, e.g., unsuitable sediment type or breeding habitat. If this discontinuity reflects a true distribution gap, it would suggest the presence of allopatric populations of T. tridentatus along each side of the Yellow Sea and future studies should aim to find more evidence for a distribution gap in this area. The future prediction of T. tridentatus shows less suitable habitat in several areas, however, none of these areas lie within the current known distribution of the species.
All Asian horseshoe crab species are known to occur at lower salinities than reported by this study, including T. tridentatus, as juveniles have been found at salinities of 8–24 PSU at a nursery beach in Hong Kong (Chiu and Morton, 2004). The temperature range of T. tridentatus is wider than of the two other Asian species, which reflects the wider latitudinal range of T. tridentatus, however, ecdysis in juveniles is halted at temperatures below 22°C (Lee and Morton, 2005). Thus, temperature is likely an important factor limiting the distribution of T. tridentatus to the north. Our results show that T. tridentatus experiences a high level and narrow range of pH, which might be explained by the interaction between temperature and pH. Lower temperatures result in higher pH (Zumdahl and Zumdahl, 2000), and thus the colder waters around Japan might have a higher pH compared to warmer waters surrounding more southern Asian countries. The chlorophyll A levels found for T. tridentatus does not differ from the two other Asian horseshoe crabs. Our results indicate that chlorophyll A concentration, in addition to temperature, is an important ecological factor determining the distributions of the Asian horseshoe crabs. This is supported by Hsieh and Chen (2009), who found that high densities of juvenile T. tridentatus were correlated with high chlorophyll A levels and high densities of polychaetes. They argued that chlorophyll A concentration reflected the amount of microalgae and hence the abundance of food available for higher trophic levels including polychaetes and that it was likely that polychaetes were an important food source for the juvenile horseshoe crabs. Thus, food resources for juvenile T. tridentatus increase with chlorophyll A concentration.
Potential Areas for Protection of Horseshoe Crabs
We suggest four criteria for identifying areas apt for conservation of horseshoe crabs in general, inspired from previous conservation studies (Petit et al., 1998; Hannah et al., 2007; Botero-Delgadillo et al., 2012). Marine Protected Areas (MPAs) with focus on horseshoe crab conservation should be prioritized in (1) geographic regions with disconnected distribution ranges, indicating physically distinct (i.e., allopatric) populations; (2) regions separated by large population breaks known from previous population genetic or phylogeographic studies; (3) areas with overlapping distributions of two or more species of horseshoe crabs; and (4) areas where the species experience little loss of suitable habitat under future climate scenarios. As previously noted, the potential distribution of suitable habitat, as well as the future potential distribution, found in this study should be considered as hypotheses of where Asian horseshoe crabs can be found and not as definitive results (Jarnevich et al., 2015). Thus, further investigations and models of the relationship between the species and the areas suggested here should be conducted prior to making decisions regarding conservation.
Our results show a wide gap in suitable habitat between C. rotundicauda populations in the southern and northern parts of the Vietnamese coast (Figure 3A). This, as well as the lack of mangrove habitat in the region (Hong and San, 1993), suggests the existence of allopatric populations of C. rotundicauda on the southern and northern coasts of Vietnam. Additionally, there is a gap in the distribution of suitable habitat at the northeast coast of Peninsular Malaysia, and there is evidence that genetically distinct populations of C. rotundicauda reside on each side of the Thai-Malay Peninsula (Obst et al., 2012; Adibah et al., 2015). The suitable habitat distribution also shows two separate distribution areas of C. rotundicauda on Borneo and Java, which could be allopatric to the three distinct distribution sites at mainland Southeast Asia. Consequently, based on our models, we hypothesize that there are probably five separated distribution ranges of C. rotundicauda in Southeast Asia, i.e., India to South Peninsular Malaysia, East Thailand to South Vietnam, North Vietnam to South China, Java, and Borneo, and none of these areas will become less suitable in the future according to our study. Based on our criteria 1, 2, and 4, we suggest the establishment of MPAs in coastal areas of West and East Thailand, North Vietnam and South China, and Borneo, in order to conserve C. rotundicauda in Southeast Asia (Figure 6).
Figure 6. Map of Southeast Asia and Japan showing regions which are optimal for the establishment of MPAs in order to conserve the three Asian species of horseshoe crabs. Rectangles with diagonal lines and green color indicate areas where MPAs for horseshoe crab conservation are absent, whereas rectangles with vertical lines and purple color indicate areas where MPAs for horseshoe crab conservation are already established.
Contrary to C. rotundicauda, we do not expect allopatric populations of T. gigas to exist across the Thai-Malay Peninsula as there is a continuity of suitable habitat for T. gigas in this region (Figure 4A). This assumption is supported by the short genetic distance between T. gigas collected from Vietnam and the Andaman Sea (Obst et al., 2012). The contrast in habitat continuity between T. gigas and C. rotundicauda can be explained by different habitat preferences for the two species. C. rotundicauda is more confined to mangroves and river deltas, while T. gigas is more frequent in truly marine environments, which is likely to result in higher dispersal along the coast. However, there is evidence of limited dispersal between the T. gigas populations in Northwest and Southwest Peninsula Malaysia (Rozihan and Ismail, 2011). Thus, based on our results, we hypothesize three distinct distribution areas for T. gigas: India to South Vietnam, Java, and Borneo. Additionally, C. rotundicauda and T. gigas are co-occurring at several locations along the coasts of the Southeast Asian mainland, hence according to criteria 1, 3, and 4 we suggest establishment of MPAs at the locations along the Southeast Asian mainland as wells as Borneo (Figure 6).
For T. tridentatus our results show four distinct distribution ranges; Japan to South Korea, China to North Vietnam, Borneo to Philippines, and Java (Figure 5A). T. tridentatus has overlapping distribution ranges with the two other Asian species; T. tridentatus and C. rotundicauda co-occur in South Vietnam as well as in the Hainan and Guangxi region of the Chinese coast, and all three Asian species co-occur along the coasts of Borneo and Java. According to criteria 1, 3, and 4, establishment of MPAs should hence be recommended in Japan, North Vietnam and South China, and Borneo (Figure 6). In Japan, MPAs are already established (http://mpatlas.org/explore/), and some of these include important breeding sanctuaries for horseshoe crabs, such as the Saikai National Park in Nagasaki. For areas where C. rotundicauda occurs sympatrically with T. gigas or T. tridentatus, it is important to consider the different habitat preferences for the species, i.e., low-salinity mangroves and river estuaries for C. rotundicauda, and high salinity sandy beaches for T. gigas and T. tridentatus.
Based on our analysis we propose that at least five regions should have high priority for horseshoe crab conservation: four in Southeast Asia, based on geographically separated areas with suitable habitat as well as the co-occurrence of two or more species, and one in Japan for the protection of distinct populations of T. tridentatus (Figure 6). We recommend that future MPAs are examined on site to ensure their suitability and that specific habitats, i.e., mangroves, mudflats, or gentle slope beaches, are present in order to provide the appropriate habitat for the horseshoe crab species in question (e.g., Kwan et al., 2016), as some factors of potential importance for horseshoe crab habitat suitability, such as sediment characteristics and presence of anthropogenic structures, could not be included in our analysis. The work presented here provides an initial step in the conservation of the Asian horseshoe crabs species, a research area in need of increasing attention considering the importance of these coastal organisms to human subsistence and health.
Future Studies
Destruction of breeding habitats pose a significant threat to horseshoe crabs (e.g., Mishra, 2009; Zaldívar-Rae et al., 2009; Nelson et al., 2016) and adult horseshoe crabs are harvested on the beaches during spawning for commercial exploitation. Since our study focused on the distribution range of adult and juvenile horseshoe crabs, we can only draw few conclusions about the distribution and preferred environmental conditions of horseshoe crab breeding sites. However, the approach implemented in this study could also be used to inventory and predict potential breeding sites using scientific networks and fieldwork, but should then include additional environmental variables. Horseshoe crabs nest on beaches with slight slopes and in mangroves and estuaries. Hence data on habitat types, beach inclination, sediment type and granularity, wave action, and organic contents of the sediment are essential variables that need to be measured and included for such analysis. Additionally, our study includes tidal regime environmental variables in the ecological niche modeling that can be used in future research on the distributions of other intertidal organisms, and thus could be of importance for the conservation of intertidal communities at a global scale.
Conclusions
The distributions of the three Asian horseshoe crab species have been assembled with this work, and marine areas of suitable habitat have been identified. Additionally, our models predict that C. rotundicauda and T. gigas in the future will experience increased habitat suitability. The information on ecological variables important for the three species, as well as the current and future distributions of suitable habitat provided here, can be used as an initial step in determining where MPAs should be established for conservation purposes and to determine the IUCN conservation status of the three Asian horseshoe crab species.
Author Contributions
SV, PF, and MO: designed the study; SV, MO, II, and PF: collected the data; MO and FQ: performed the modeling and analysis; SV: led the writing of the first draft in close collaboration with PF and MO. All authors contributed substantially to the revisions.
Funding
This study is supported by the EU's Seventh Framework Program project BioVeL (www.biovel.eu, grant no. 283359), the Horizon 2020 program ENVRI+ (grant agreement no. 654182), the Swedish LifeWatch program funded by the Swedish Research Council (grant no. 829-2009-6278), and by the Danish National Research Foundation Niels Bohr professorship project Aarhus University Research on the Anthropocene (AURA).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank Dr. Shin Nishida, Dr. Hwey-Lian Hsieh, and Joseph Klein for valuable feedback and help to improve this paper. SV and PF consider this work a contribution to the Danish National Research Foundation Niels Bohr professorship project Aarhus University Research on the Anthropocene (AURA).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2018.00164/full#supplementary-material
Supplementary Data Sheet 1. Explanation of the construction of global tidal variables.
Supplementary Image 1. Heat map for C. rotundicauda using Maxent.
Supplementary Image 2. Heat map for C. rotundicauda using SVM.
Supplementary Image 3. Heat map for T. gigas using Maxent.
Supplementary Image 4. Heat map for T. gigas using SVM.
Supplementary Image 5. Heat map for T. tridentatus using Maxent.
Supplementary Image 6. Heat map for T. tridentatus using SVM.
Supplementary Table 1. Data provider names for occurrence records obtained by own observations, collaborators, and scientific network.
Supplementary Table 2. Data references for occurrence records obtained by literature.
Supplementary Table 3. Data provider names for obtained GBIF occurrence records. The name (or acronym) in use by the institution having custody of the object(s) or information referred to in the record.
Supplementary Table 4. Test statistics for analysis of differences in environmental conditions between four horseshoe crab species using Welch's test for unequal variances.
Supplementary Table 5. Multiple comparisons of difference between means of environmental variables for four horseshoe crab species (CR – C. rotundicauda, TG – T. gigas, TT – T. tridentatus) using Games-Howell Post-Hoc test. Significant mean difference at the 0.05 level in bold.
Supplementary Table 6. Loadings for the different environmental variables shown for all PCA components for C. rotundicauda.
Supplementary Table 7. Loadings for the different environmental variables shown for all PCA components for T. gigas.
Supplementary Table 8. Loadings for the different environmental variables shown for all PCA components for T. tridentatus.
References
Adibah, A. B., Ng, W. L., and Tan, S. G. (2015). The Malay Peninsula as a barrier to gene flow in an Asian horseshoe crab species, Carcinoscorpius rotundicauda Latreille. Biochem. Syst. Ecol. 60, 204–210. doi: 10.1016/j.bse.2015.04.026
Albouy, C., Guilhaumon, F., Leprieur, F., Lasram, F. B. R., Aznar, S. S. R., Velez, L., et al. (2013). Projected climate change and the changing biogeography of coastal Mediterranean fishes. J. Biogeogr. 40, 534–547. doi: 10.1111/jbi.12013
Allison, E. H., Perry, A. L., Badjeck, M.-C., Adger, W. N., Brown, K., and Conway, D. (2009). Vulnerability of national economies to the impacts of climate change on fisheries. Fish Fish. 10, 173–196. doi: 10.1111/j.1467-2979.2008.00310.x
Almendral, M. A., and Schoppe, S. (2005). Population structure of Tachypleus tridentatus (Chelicerata: Merostomata) at a nursery beach in Puerto Princesa City, Palawan, Philippines. J. Nat. Hist. 39, 2319–2329. doi: 10.1080/00222930500063219
Anderson, L. I., and Shuster, C. N. Jr. (2003). “Throughout geologic time: where have they lived,” in The American Horseshoe Crab, eds C. N. Shuster Jr, R. B. Barlow, and H. J. Brockmann (Cambridge: Harvard University Press), 189–223.
Barlow, R. B. Jr., Powers, M. K., Howard, H., and Kass, L. (1986). Migration of Limulus for mating: relation to lunar phase, tide height, and sunlight. Biol. Bull. 171, 310–329. doi: 10.2307/1541674
Behera, S., Tripathy, B., Sivakumar, K., Choudhury, B. C., and Bhadury, P. (2015). “Distribution and abundance of two sympatric species of horseshoe crabs along the Odisha coast, India,” in Changing Global Perspectives on Horseshoe Crab Biology, Conservation and Management, eds R. H. Carmichael, M. L. Botton, P. K. S. Shin, and S. G. Cheung (New York, NY: Springer), 181–191.
Bell, J. D., Kronen, M., Vunisea, A., Nash, W. J., Keeble, G., Demmke, A., et al. (2009). Planning the use of fish for food security in the Pacific. Mar. Policy 33, 64–76. doi: 10.1016/j.marpol.2008.04.002
Bonaventura, J., Bonaventura, C., and Tesh, S. (eds.). (1982). Physiology and Biology of Horseshoe Crabs; Studies on Normal and Environmentally Stressed Animals Progress in Clinical and Biological Research, Vol. 81. New York, NY: Alan R. Liss, Inc.
Botero-Delgadillo, E., Páez, C. A., and Bayly, N. (2012). Biogeography and conservation of Andean and Trans-Andean populations of Pyrrhura parakeets in Colombia: modelling geographic distributions to identify independent conservation units. Bird Conserv. Int. 22, 445–461. doi: 10.1017/S0959270912000196
Botton, M. L., Shuster, C. N., Sekiguchi, K., and Sugita, H. (1996). Amplexus and mating behavior in the Japanese horseshoe crab, Tachypleus tridentatus. Zool. Sci. 13, 151–159. doi: 10.2108/zsj.13.151
Botton, M. L. (2001). “The conservation of horseshoe crabs: What can we learn from the Japanese experience?” in Limulus in the Limelight, ed J. T. Tanacredi (New York, NY: Springer), 41–51.
Botton, M. L., Carmichael, R. H., Shin, P. K. S., and Cheung, S. G. (2015). “Emerging issues in horseshoe crab conservation: a perspective from the IUCN Species Specialist Group,” in Changing Global Perspectives on Horseshoe Crab Biology, Conservation And Management, eds R. H. Carmichael, M. L. Botton, P. K. S. Shin, and S. G. Cheung (New York, NY: Springer), 369–381.
Botton, M. L., Loveland, R. E., and Jacobsen, T. R. (1988). Beach erosion and geochemical factors: influence on spawning success of horseshoe crabs (Limulus polyphemus) in Delaware Bay. Mar. Biol. 99, 325–332. doi: 10.1007/BF02112124
Boyer, T. P., Antonov, J. I., Baranova, O. K., Garcia, H. E., Johnson, D. R., Locarnini, R. A., et al. (2009). “World Ocean Database 2009,” in NOAA Atlas NESDIS 66, ed S. Levitus (Washington, DC: U.S. Gov. Printing Office), 1–217.
Cai, L., Li, H., Zhang, J., and Rao, Y. (2015). Characterization of twelve new microsatellite markers for the endangered Chinese horseshoe crab (Tachypleus tridentatus) through next-generation sequencing. J. Crustacean Biol. 35, 476–479. doi: 10.1163/1937240X-00002331
Carrère, L., Carrère, L., Lyard, F., Cancet, M., Roblou, L., and Guillot, A. (2012). “FES2012: A new global tidal model taking advantage of nearly 20 years of altimetry,” in Proceedings of Meeting “20 Years of Altimerty,” (Venice).
Cartwright-Taylor, L., Bing, Y. V., Chi, H. C., and Tee, L. S. (2011). Distribution and abundance of horseshoe crabs Tachypleus gigas and Carcinoscorpius rotundicauda around the main island of Singapore. Aquat. Biol. 13, 127–136. doi: 10.3354/ab00346
Cartwright-Taylor, L., Lee, J., and Hsu, C. C. (2009). Population structure and breeding pattern of the mangrove horseshoe crab Carcinoscorpius rotundicauda in Singapore. Aquat. Biol. 8, 61–69. doi: 10.3354/ab00206
Chatterji, A. (1999). New record of the sympatric distribution of two Asian species of the horseshoe crab. Curr. Sci. 77, 746–747.
Chatterji, A., Kotnala, S., and Mathew, R. (2004). Effect of salinity on larval growth of horseshoe crab, Tachypleus gigas (Müller). Curr. Sci. 87, 248–250. Available online at: http://www.jstor.org/stable/24108874
Chatterji, A., Mishra, J. K., and Parulekar, A. H. (1992). Feeding behaviour and food selection in the horseshoe crab, Tachypleus gigas (Müller). Hydrobiologia 246, 41–48. doi: 10.1007/BF00005621
Chatterji, A., and Parulekar, A. H. (1992). Fecundity of the Indian horse-shoe crab, Carcinoscorpius rotundicauda (Latreille). Trop. Ecol. 33, 97–102.
Chatterji, A., and Pati, S. (2014). Allometric relationship in the adult population of the Malaysian horseshoe crab (Tachypleus tridentatus; Leach). Int. J. Res. 1, 1378–1385.
Chen, C.-P., Yang, M.-C., Fan, L.-F., Qiu, G., Liao, Y.-Y., and Hsieh, H.-L. (2015). Co-occurrence of juvenile horseshoe crabs Tachypleus tridentatus and Carcinoscorpius rotundicauda in an estuarine bay, southwestern China. Aquat. Biol. 24, 117–126. doi: 10.3354/ab00641
Chiu, H. M. C., and Morton, B. (2003). The morphological differentiation of two horseshoe crab species, Tachypleus tridentatus and Carcinoscorpius rotundicauda (Xiphosura), in Hong Kong with a regional Asian comparison. J. Nat. History 37, 2369–2382. doi: 10.1080/00222930210149753
Chiu, H. M. C., and Morton, B. (2004). The behaviour of juvenile horseshoe crabs, Tachypleus tridentatus (Xiphosura), on a nursery beach at Shui Hau Wan, Hong Kong. Hydrobiologia 523, 29–35. doi: 10.1023/B:HYDR.0000033085.71861.63
Dao, H. V., Takata, Y., Sato, S., Fukuyo, Y., and Kodama, M. (2009). Frequent occurrence of the tetrodotoxin-bearing horseshoe crab Carcinoscorpius rotundicauda in Vietnam. Fish. Sci. 75, 435–438. doi: 10.1007/s12562-008-0041-5
Defeo, O., McLachlan, A., Schoeman, D. S., Schlacher, T. S., Dugan, J., Jones, A., et al. (2009). Threats to sandy beach ecosystems: a review. Estuar. Coast. Shelf Sci. 81, 1–12. doi: 10.1016/j.ecss.2008.09.022
De Giovanni, R., Williams, A. R., Ernst, V. H., Kulawik, R., Fernandez, F. Q., and Hardisty, A. R. (2015). ENM components: a new set of web service-based workflow components for ecological niche modelling. Ecography 39, 376–383. doi: 10.1111/ecog.01552
de Souza Muñoz, M. E., De Giovanni, R., de Siqueira, M. F., Sutton, T., Brewer, P., Pereira, R. S., et al. (2011). openModeller: a generic approach to species potential distribution modelling. Geoinformatica 15, 111–135. doi: 10.1007/s10707-009-0090-7
Ding, J. L., Navas, M. A. A. III., and Ho, B. (1995). Molecular cloning and sequence analysis of Factor C cDNA from the Singapore horseshoe crab, Carcinoscorpius rotundicauda. Mol. Marine Biol. Biotechnol. 4, 90–103.
Drake, J. M., Randin, C., and Guisan, A. (2006). Modelling ecological niches with support vector machines. J. Appl. Ecol. 43, 424–432. doi: 10.1111/j.1365-2664.2006.01141.x
Ehlinger, G. S., and Tankersley, R. A. (2004). Survival and development of horseshoe crab (Limulus polyphemus) embryos and larvae in hypersaline conditions. Biol. Bull. 206, 87–94. doi: 10.2307/1543539
Elith, J., Graham, C. H., Anderson, R. P., Dudík, M., Ferrier, S., Guisan, A., et al. (2006). Novel methods improve prediction of species distributions from occurrence data. Ecography 29, 129–151. doi: 10.1111/j.2006.0906-7590.04596.x
Faizul, M. I. M., Eng, H. T., Christianus, A., Abdel-Hadi, Y. M., et al. (2015). “Bacteria and fungi identified on horseshoe crabs, Tachypleus gigas and Carcinoscorpius rotundicauda in the laboratory,” in Changing Global Perspectives on Horseshoe Crab Biology, Conservation and Management, ed R. H. Carmichael, M. L. Botton, P. K. S. Shin, and S. G. Cheung (New York, NY: Springer), 303–311.
Faurby, S., Sauer, K., Nielsen, K., Bussarawit, S., Intanai, I., Cong, N., et al. (2011). Intraspecific shape variation in horseshoe crabs: the importance of sexual and natural selection for local adaptation. J. Exp. Mar. Biol. Ecol. 407, 131–138. doi: 10.1016/j.jembe.2011.05.025
Funch, P., Wang, T., Pertoldi, C., and Middelbo, A. B. (2016). Low oxygen levels slow embryonic development of Limulus polyphemus. Biol. Bull. 231, 113–119. doi: 10.1086/690091
Fusetani, N., Endo, H., Hashimoto, K., and Takahashi, K. (1982). Occurrence of potent toxins in the horseshoe crab Carcinoscorpius rotundicauda. Toxicon 20, 662–664. doi: 10.1016/0041-0101(82)90061-7
Gauvry, G. (2015). “Current horseshoe crab harvesting practices cannot support global demand for TAL/LAL: The pharmaceutical and medical device industries role in the sustainability of horseshoe crabs,” in Changing Global Perspectives on Horseshoe Crab Biology, Conservation and Management, eds R. H. Carmichael, M. L. Botton, P. K. S. Shin, and S. G. Cheung (New York, NY: Springer), 475–482.
Gilman, E. L., Ellison, J., Duke, N. C., and Field, C. (2008). Threats to mangroves from climate change and adaptation options: a review. Aquat. Bot. 89, 237–250. doi: 10.1016/j.aquabot.2007.12.009
Graham, C. H., Elith, J., Hijmans, R. J., Guisan, A., Peterson, A. T., Loiselle, B. A., et al. (2008). The influence of spatial errors in species occurrence data used in distribution models. J. Appl. Ecol. 45, 239–247. doi: 10.1111/j.1365-2664.2007.01408.x
Hannah, L., Midgley, G., Andelman, S., Araújo, M., Hughes, G., Martinez-Meyer, E., et al. (2007). Protected area needs in a changing climate. Front. Ecol. Environ. 5, 131–138. doi: 10.1890/1540-9295(2007)5[131:PANIAC]2.0.CO;2
Hardisty, A. R., Bacall, F., Beard, N., Balcázar-Vargas, M. P., Balech, B., Barcza, Z., et al. (2016). BioVeL: a virtual laboratory for data analysis and modelling in biodiversity science and ecology. BMC Ecol. 16:49. doi: 10.1186/s12898-016-0103-y
Holl, S., Garijo, D., Belhajjame, K., Zimmermann, O., De Giovanni, R., Obst, M., et al. (2013). “On specifying and sharing scientific workflow optimization results using research objects,” in Proceedings of the 8thWorkshop on Workflows in Support of Large-Scale Science (Denver, CO: ACM), 28–37.
Hong, R. F. (2004). Population and distribution of horseshoe crab Carcinoscorpius rotundicauda at the Kranji Nature Trail estuaries, Western Johor Straits. The National University of Singapore, Department of Biological Sciences, Project report (Singapore).
Hsieh, H.-L., and Chen, C.-P. (2009). “Conservation program for the Asian horseshoe crab Tachypleus tridentatus in Taiwan: characterizing the microhabitat of nursery grounds and restoring spawning grounds,” in Biology and Conservation of Horseshoe Crabs, eds J. T. Tanacredi, M. L. Botton, and D. R. Smith (New York, NY: Springer), 417–438.
Hu, M., Wang, Y., Chen, Y., Cheung, S.-G., Shin, P. K. S., and Li, Q. (2009). Summer distribution and abundance of juvenile Chinese horseshoe crabs Tachypleus tridentatus along an intertidal zone in southern China. Aquat. Biol. 7, 107–112. doi: 10.3354/ab00194
IBM Corporation (2013). IBM SPSS Statistics for Macintosh Version 22.0. Armonk, NY: IBM Corporation.
IPCC (2007). Climate Change 2007: Synthesis Report. Contribution of Working Groups I, II and III to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, eds Core Writing Team, R. K. Pachauri, and A. Reisinger (Geneva: IPCC).
Ismail, N., and Sarijan, S. (2011). Phylogenetic inference from 18S rRNA gene sequences of horseshoe crabs, Tachypleus gigas between Tanjung Dawai, Kedah & Cherating, Pahang, Peninsular Malaysia. Int. J. Biol. Food Vet. Agric. Eng. 5, 66–69. doi: 10.1999/1307-6892/8914
Itow, T., Loveland, R. E., and Botton, M. L. (1998). Developmental abnormalities in horseshoe crab embryos caused by exposure to heavy metals. Arch. Environ. Contam. Toxicol. 35, 33–40. doi: 10.1007/s002449900345
IUCN (2013). IUCN-led Intervention Stops Import of Asian Horseshoe Crabs Into USA. Available online at: http://www.iucnredlist.org/news/intervention-led-by-iucn-experts-stops-import-of-asian-horseshoe-crabs-into-usa (Accessed Sept 01, 2016).
Jackson, N. L., Nordstrom, K. F., Saini, S., and Smith, D. R. (2015). Influence of configuration of bulkheads on use of estuarine beaches by horseshoe crabs and foraging shorebirds. Environ. Earth Sci. 74, 5749–5758. doi: 10.1007/s12665-015-4592-3
Jackson, N. L., Smith, D. R., and Nordstrom, K. F. (2008). Physical and chemical changes in the foreshore of an estuarine beach: implications for viability and development of horseshoe crab Limulus polyphemus eggs. Mar. Ecol. Prog. Ser. 355, 209–218. doi: 10.3354/meps07211
Jarnevich, C. S., Stohlgren, T. J., Kumar, S., Morisette, J. T., and Holcombe, T. R. (2015). Caveats for correlative species distribution modeling. Ecol. Inform. 29, 6–15. doi: 10.1016/j.ecoinf.2015.06.007
Jawahir, A. R. N., Samsur, M., Shabdin, M. L., and Khairul Adha, A. R. (2017). Distribution of two species of Asian horseshoe crabs at west coast of Sarawak's Waters, East Malaysia. Egyptian J. Aquat. Res. 43, 135–140. doi: 10.1016/j.ejar.2017.03.002
Jeffries, W. B., Voris, H. K., and Yang, C. M. (1989). Observations on the incidence of the pedunculate barnacle, Octolasmis warwickii (Gray, 1825) on horseshoe crabs (Xiphosura) in the seas adjacent to Singapore. Raffles Bull. Zool. 37, 58–62.
Jegla, T. C., and Costlow, J. D. (1982). Temperature and salinity effects on developmental and early posthatch stages of Limulus. Prog. Clin. Biol. Res. 81, 103–113.
John, B. A., Jalal, K. C. A., Zaleha, K., Armstrong, P., and Kmaruzzaman, B.Y. (2011). Effects of blood extraction on the mortality of Malaysian horseshoe crabs (Tachypleus gigas). Mar. Freshw. Behav. Physiol. 44, 321–327. doi: 10.1080/10236244.2011.642505
Kamaruzzaman, B. Y., Akbar John, B., Aqilah Megat, M. H., and Zaleha, K. (2011). Bioaccumulation of heavy metals in horseshoe crabs (Tachypleus gigas) from Pekan, Pahang, Malaysia. Res. J. Environ. Toxicol. 5, 222–228. doi: 10.3923/rjet.2011.222.228
Kannan, K., Yasunaga, Y., Iwata, H., Ichihashi, H., Tanabe, S., and Tatsukawa, R. (1995). Concentrations of heavy metals, organochlorines, and organotins in horseshoe crab, Tachypleus tridentatus, from Japanese coastal waters. Arch. Environ. Contam. Toxicol. 28, 40–47. doi: 10.1007/BF00213967
Kaschner, K., Ready, J. S., Agbayani, E., Rius, J., Kesner-Reyes, K., Eastwood, P. D., et al. (2008). AquaMaps Environmental Dataset: Half-Degree Cells Authority File (HCAF). Available online at: http://www.aquamaps.org/download/main.php
Ketchum, B. H. (1972). The Waters Edge: Critical Problems of the Coastal Zone. Boston, MA: MIT Press.
Key, M. M. Jr., Jeffries, W. B., Voris, H. K., and Yang, C. M. (1996). Epizoic bryozoans, horseshoe crabs, and other mobile benthic substrates. Bull. Mar. Sci. 58, 368–384.
Key, M. M. Jr., Jeffries, W. B., Voris, H. K., and Yang, C. M. (2000). “Bryozoan fouling pattern on the horseshoe crab Tachypleus gigas (Müller) from Singapore,” in Proceedings of the 11th International Bryozoology Association Conference, eds A. Herrera and J. B. C. Jackson (Balboa: Smithsonian Tropical Research Institute), 265–271.
Kleisner, K. M., Fogarty, M. J., McGee, S., Hare, J. A., Moret, S., Perretti, C. T., et al. (2017). Marine species distribution shifts on the US Northeast Continental Shelf under continued ocean warming. Prog. Oceanogr. 153, 24–36. doi: 10.1016/j.pocean.2017.04.001
Kremen, C., Cameron, A., Moilanen, A., Phillips, S. J., Thomas, C. D., Beentje, H., et al. (2008). Aligning conservation priorities across taxa in Madagascar with high-resolution planning tools. Science 320, 222–226. doi: 10.1126/science.1155193
Kungsuwan, A., Nagashima, Y., Noguchi, T., Shida, Y., Suvapeepan, S., Suwansakornkul, P., et al. (1987). Tetrodotoxin in the horseshoe crab Carcinoscorpius rotundicauda inhabiting Thailand. Nippon Suisan Gakkaishi 53, 261–266. doi: 10.2331/suisan.53.261
Kwan, B. K. Y., Hsieh, H.-L., Cheung, S. G., and Shin, P. K. S. (2016). Present population and habitat status of potentially threatened Asian horseshoe crabs Tachypleus tridentatus and Carcinoscorpius rotundicauda in Hong Kong: a proposal for marine protected areas. Biodivers. Conserv. 25, 673–692. doi: 10.1007/s10531-016-1084-z
Kwan, B. K. Y., Cheung, S. G., and Shin, P. K. S. (2015). A dual stable isotope study for diet composition of juvenile Chinese horseshoe crab Tachypleus tridentatus (Xiphosura) on a seagrass-covered intertidal mudflat. Mar. Biol. 162, 1137–1143. doi: 10.1007/s00227-015-2647-3
Laughlin, R. (1983). The effects of temperature and salinity on larval growth of the horseshoe crab Limulus polyphemus. Biol. Bull. 164, 93–103. doi: 10.2307/1541193
Lee, C. N., and Morton, B. (2005). Experimentally derived estimates of growth by juvenile Tachypleus tridentatus and Carcinoscorpius rotundicauda (Xiphosura) from nursery beaches in Hong Kong. J. Exp. Mar. Biol. Ecol. 318, 39–49. doi: 10.1016/j.jembe.2004.12.010
Leidenberger, S., De Giovanni, R., Kulawik, R., Williams, A. R., Bourlat, S. J. and Maggs, C. (2015a). Mapping present and future potential distribution patterns for a meso-grazer guild in the Baltic Sea. J. Biogeogr. 42, 241–254. doi: 10.1111/jbi.12395
Leidenberger, S., Obst, M., Kulawik, R., Stelzer, K., Heyer, K., Hardisty, A., et al. (2015b). Evaluating the potential of ecological niche modelling as a component in marine non-indigenous species risk assessments. Mar. Pollut. Bull. 97, 470–487. doi: 10.1016/j.marpolbul.2015.04.033
Levin, J., and Bang, F. B. (1968). Clottable protein in Limulus; its localization and kinetics of its coagulation by endotoxin. Thromb. Diath. Haemorrh. 19, 186–197.
Levin, J., Hochstein, H. D., and Novitsky, T. J. (2003). “Clotting cells and Limulus Amebocyte Lysate: An amazing analytical tool,” in The American Horseshoe Crab, eds C. N. Shuster, R. B. Barlow, and H. J. Brockmann (Cambridge: Harvard University Press). 310–340.
Liao, Y., Chen, C.-P., Hsieh, H.-L., Cao, Y., and Chen, J. (2012). Sallow-skin horseshoe crabs (late juvenile Tachypleus tridentatus) as osmoconformers. J. Mar. Biol. Assoc. U.K. 92, 463–468. doi: 10.1017/S0025315411001196
Liew, P. L., Ng, W. L., and Tan, S. G. (2015). Levels and patterns of genetic variation in an Asian horseshoe crab species, Tachypleus gigas Müller, from the Malay Peninsula. Mar. Biol. Res. 11, 879–886. doi: 10.1080/17451000.2015.1024135
Lonza (2016). PyroGene™ Recombinant Factor C Assay Endotoxin Testing: Its Time to Embrace the Alternative. Lonza Walkersville Inc. Available online at: https://www.lonza.com/products-services/pharma-biotech/endotoxin-detection/endotoxin-detection-assays/recombinant-factor-c-assay.aspx
Loveland, R. E., and Botton, M. L. (2015). “Sea level rise in Delaware Bay, USA: adaptations of spawning horseshoe crabs (Limulus polyphemus) to the glacial past, and the rapidly changing shoreline of the bay,” in Changing Global Perspectives on Horseshoe Crab Biology, Conservation and Management, eds R. H. Carmichael, M. L. Botton, P. K. S. Shin, and S. G. Cheung (New York, NY: Springer), 41–63.
Lyard, F., Lefevre, F., Letellier, T., and Francis, O. (2006). Modelling the global ocean tides: modern insights from FES2004. Ocean Dyn. 56, 394–415. doi: 10.1007/s10236-006-0086-x
Mathew, C., Güntsch, A., Obst, M., Vicario, S., Haines, R., Williams, A. R., et al. (2014). A semi-automated workflow for biodiversity data retrieval, cleaning, and quality control. Biodiversity Data J. 2:e4221. doi: 10.3897/BDJ.2.e4221
Metzger, M. J., Bunce, R. G. H., Jongman, R. H. G., Mücher, C. A., and Watkins, J. W. (2005). A climatic stratification of the environment of Europe. Glob. Ecol. Biogeogr. 14, 549–563. doi: 10.1111/j.1466-822X.2005.00190.x
Mishra, J. K. (2009). “Horseshoe crabs, their eco-biological status along the northeast coast of India and the necessity for ecological conservation,” in Biology and Conservation of Horseshoe Crabs, eds J. T. Tanacredi, M. L. Botton, and D. R. Smith (New York, NY: Springer), 89–96.
Nelson, B. R., Satyanarayana, B., Moh, J. H. Z., Ikhwanuddin, M., Chatterji, A., and Shaharom, F. (2016). The final spawning ground of Tachypleus gigas (Müller, 1785) on the east Peninsular Malaysia is at risk: a call for action. PeerJ 4:e2232. doi: 10.7717/peerj.2232
Nguyen, V. C., and Dao, T. H. (1995). Zoobenthos fauna in the sea waters of Vietnam. Collect. Mar. Res. Works 6, 91–100.
Ngy, L., Yu, C. F., Takatani, T., and Arakawa, O. (2007). Toxicity assessment for the horseshoe crab Carcinoscorpius rotundicauda collected from Cambodia. Toxicon 49, 843–847. doi: 10.1016/j.toxicon.2006.12.004
Nishida, S., Kuroyanagi, K., and Koike, H. (2015). “Genetic features of Tachypleus tridentatus in Japan and an alien (non-native) population founded at Ise-Mikawa Bay, Chubu region, central Japan,” in Changing Global Perspectives on Horseshoe Crab Biology, Conservation and Management, eds R. H. Carmichael, M. L. Botton, P. K. S. Shin, and S. G. Cheung (New York, NY: Springer), 129–146.
Nix, H. A. (1986). “A biogeographic analysis of Australian elapid snakes,” in Atlas of Australian elapid snakes (Australian Flora and Fauna Series Number 7, ed R. Longmore. Canberra, ACT: Australian Government Publishing Service, 4–15.
Nordstrom, K. F. (2004). Beaches and Dunes of Developed Coasts. Cambridge: Cambridge University Press.
Obst, M., Faurby, S., Bussarawit, S., and Funch, P. (2012). Molecular phylogeny of extant horseshoe crabs (Xiphosura, Limulidae) indicates Paleogene diversification of Asian species. Mol. Phylogenet. Evol. 62, 21–26. doi: 10.1016/j.ympev.2011.08.025
Palumbi, S. R., and Johnson, B. A. (1982). A note on the influence of life-history stage on metabolic adaptation: the responses of Limulus eggs and larvae to hypoxia. Prog. Clin. Biol. Res. 81, 115–124.
Patil, J. S., and Anil, A. C. (2000). Epibiotic community of the horseshoe crab Tachypleus gigas. Mar. Biol. 136, 699–713. doi: 10.1007/s002270050730
Pearson, R. G., Thuiller, W., Araújo, M. B., Martinez-Meyer, E., Brotons, L., McClean, C., et al. (2006). Model-based uncertainty in species range prediction. J. Biogeogr. 33, 1704–1711. doi: 10.1111/j.1365-2699.2006.01460.x
Penn, D., and Brockmann, H. J. (1994). Nest-site selection in the horseshoe crab, Limulus polyphemus. Biol. Bull. 187, 373–384. doi: 10.2307/1542294
Peterson, A. T., Soberón, J., Pearson, R. G., Anderson, R. P., Martínez-Meyer, E., Nakamura, M., et al. (2011). Ecological Niches and Geographic Distributions. Princeton, NJ: Princeton University Press.
Petit, R. J., El Mousadik, A., and Pons, O. (1998). Identifying populations for conservation on the basis of genetic markers. Conserv. Biol. 12, 844–855. doi: 10.1046/j.1523-1739.1998.96489.x
Phillips, S. J., Anderson, R. P., and Schapire, R. E. (2006). Maximum entropy modeling of species geographic distributions. Ecol. Modell. 190, 231–259. doi: 10.1016/j.ecolmodel.2005.03.026
Phillips, S. J., and Dudík, M. (2008). Modeling of species distributions with Maxent: new extensions and a comprehensive evaluation. Ecography 31, 161–175. doi: 10.1111/j.0906-7590.2008.5203.x
Porfirio, L. L., Harris, R. M., Lefroy, E. C., Hugh, S., Gould, S. F., Lee, G., et al. (2014). Improving the use of species distribution models in conservation planning and management under climate change. PLoS ONE 9:e113749. doi: 10.1371/journal.pone.0113749
Purcell, S. W., Mercier, A., Conand, C., Hamel, J.-F., Toral-Granda, M. V., Lovatelli, A., et al. (2013). Sea cucumber fisheries: global analysis of stocks, management measures and drivers of overfishing. Fish Fish. 14, 34–59. doi: 10.1111/j.1467-2979.2011.00443.x
Quantum GIS Development Team (2014). Quantum GIS Geographic Information System. Open Source Geospatial Foundation Project. Available online at: http://qgis.org
Raman, N. J. B. A., Mohamad, S., Rahim, K. A. A., and Long, S. M. (2014). “Preliminary studies of horseshoe crabs in Sarawak, Malaysia,” International Conference on Marine Science and Aquaculture 2014: Ecosystem Perspectives in Sustainable Development (Kota Kinabalu: Universiti Malaysia Sabah), 67–74.
Rao, D. V., and Bhagirathi, B. (1989). Carcinoscorpius rotundicauda amoebocyte lysate for detection of endotoxins–its preparation, stability, sensitivity and comparison with Limulus amoebocyte lysate. Indian J. Exp. Biol. 27, 26–31.
R Core Team (2013). R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna.
Ready, J. S., Kaschner, K., South, A., Eastwood, P., Rees, T., Rius, J., et al. (2010). Predicting the distributions of marine organisms at the global scale. Ecol. Modell. 221, 467–478. doi: 10.1016/j.ecolmodel.2009.10.025
Reiss, H., Birchenough, S., Borja, A., Buhl-Mortensen, L., Craeymeersch, J., Dannheim, J., et al. (2014). Benthos distribution modelling and its relevance for marine ecosystem management. ICES J. Mar. Sci. 72, 297–315. doi: 10.1093/icesjms/fsu107
Robert, R., Muhammad Ali, S. H., and Amelia-Ng, P. F. (2014). Demographics of horseshoe crab populations in Kota Kinabalu, Sabah, Malaysia with emphasis on Carcinoscorpius rotundicauda and some aspects of its mating behaviour. Pertanika J. Trop. Agric. Sci. 37, 375–388. Available online at: http://www.pertanika.upm.edu.my/Pertanika%20PAPERS/JTAS%20Vol.%2037%20(3)%20Aug.%202014/08%20JTAS%200583-2013.pdf
Rojana-anawat, P., Pradit, S., Sukramongkol, N., and Siriraksophon, S. (2001). “Temperature, salinity, dissolved oxygen and water masses of Vietnamese waters,” Proceedings of the SEAFDEC Seminar on Fishery Resources in the South China Sea, Area IV: Vietnamese Waters, Vol. 4, 346–355.
Rozihan, M., and Ismail, E. (2011). Genetic structure and haplotype diversity of Tachypleus gigas population along the West Coast of peninsular Malaysia-inferred through mtDNA AT rich region sequence analysis. Biotechnology 10, 298–302. doi: 10.3923/biotech.2011.298.302
Rozihan, M., and Ismail, E. (2012). Impact of Malaysian continental drift on the genetic diversity of horseshoe crab inferred through mtDNA sequence analysis. Int. J. Biol. 4, 104–110. doi: 10.5539/ijb.v4n1p104
Rudloe, A. (1985). Variation in the expression of lunar and tidal behavioral rhythms in the horseshoe crab, Limulus polyphemus. Bull. Mar. Sci. 36, 388–395.
Sahu, A. C., and Dey, L. (2013). Spawning density and morphometric characteristics of the horseshoe crab Tachypleus gigas (Müller) on the Balasore coast of Bay of Bengal, India. Sci. Vis. 13, 76–84.
Schölkopf, B., Platt, J. C., Shawe-Taylor, J., Smola, A. J., and Williamson, R. C. (2001). Estimating the support of a high-dimensional distribution. Neural Comput. 13, 1443–1471. doi: 10.1162/089976601750264965
Sekiguchi, K., and Shuster, C. N Jr. (2009). “Limits on the global distribution of horseshoe crabs (Limulacea): lessons learned from two lifetimes of observations: Asia and America,” in Biology and Conservation of Horseshoe Crabs, eds J. T. Tanacredi, M. L. Botton, and D. R. Smith (New York, NY: Springer), 5–24.
Shakibazadeh, S., Christianus, A., and Hajeb, P. (2013). Blood cell histology of horseshoe crab, Tachypleus gigas. Asian J. Anim. Vet. Adv. 8, 275–283. doi: 10.3923/ajava.2013.275.283
Shigenaga, T., Muta, T., Toh, Y., Tokunaga, F., and Iwanaga, S. (1990). Antimicrobial tachyplesin peptide precursor. cDNA cloning and cellular localization in the horseshoe crab (Tachypleus tridentatus). J. Biol. Chem. 265, 21350–21354.
Shuster, C. N., Barlow, R. B., and Brockmann, H. J. (eds.). (2003). The American Horseshoe Crab. Cambridge, MA: Harvard University Press.
Smith, D. R., Beekey, M. A., Brockmann, H. J., King, T. L., Millard, M. J., and Zaldívar-Rae, J. A. (2016). Limulus polyphemus. The IUCN Red List of Threatened Species 2016. e.T11987A80159830. Available online at: http://www.iucnredlist.org/details/11987/0
Smith, D. R., Brockmann, H. J., Beekey, M. A., King, T. L., Millard, M. J., and Zaldívar-Rae, J. A. (2017). Conservation status of the American horseshoe crab, (Limulus polyphemus): a regional assessment. Rev. Fish Biol. Fish. 27, 135–175. doi: 10.1007/s11160-016-9461-y
Srijaya, T. C., Pradeep, P. J., Mithun, S., Hassan, A., Shaharom, F., and Chatterji, A. (2010). A new record on the morphometric variations in the populations of horseshoe crab (Carcinoscorpius rotundicauda Latreille) obtained from two different ecological habitats of Peninsular Malaysia. Our Nat. 8, 204–211. doi: 10.3126/on.v8i1.4329
Srijaya, T. C., Pradeep, P. J., Hassan, A., Chatterji, A., Shaharom, F., and Jeffs, A. (2014). Oxygen consumption in trilobite larvae of the mangrove horseshoe crab (Carcinoscorpius rotundicauda; Latreille, 1802): effect of temperature, salinity, pH, and light–dark cycle. Int. Aquat. Res. 6:60. doi: 10.1007/s40071-014-0060-z
Stirling, D. A., Boulcott, P., Scott, B. E., and Wright, P. J. (2016). Using verified species distribution models to inform the conservation of a rare marine species. Divers. Distributions 22, 808–822. doi: 10.1111/ddi.12447
Tan, A. N., Christianus, A., and Satar, M. K. A. (2011). Epibiont infestation on horseshoe crab Tachypleus gigas (Müller) at Pantai Balok in Peninsular Malaysia. Our Nat. 9, 9–15. doi: 10.3126/on.v9i1.5725
Tanu, M. B., and Noguchi, T. (1999). Tetrodotoxin as a toxic principle in the horseshoe crab Carcinoscorpius rotundicauda collected from Bangladesh. J. Food Hygienic Soc. Japan 40, 426–430. doi: 10.3358/shokueishi.40.6_426
Thapanand-Chaidee, T., Pe, U. M., and Chanrachkij, I. (2010). Status of marine fisheries resources in Myanmar waters: estimates from bottom trawl survey. Kasetsart Univer. Fish. Res. Bull. 34, 35–48. Available online at: https://www.tci-thaijo.org/index.php/JFE/article/view/81321
Tyberghein, L., Verbruggen, H., Pauly, K., Troupin, C., Mineur, F., and De Clerck, O. (2012). Bio-ORACLE: a global environmental dataset for marine species distribution modelling. Global Ecol. Biogeogr. 21, 272–281. doi: 10.1111/j.1466-8238.2011.00656.x
Vijayakumar, R., Das, S., Chatterji, A., and Parulekar, A. H. (2000). Morphometric characteristics in the horseshoe crab Tachypleus gigas (Arthropoda: Merostomata). Indian J. Mar. Sci. 29, 333–335. Available online at: http://nopr.niscair.res.in/handle/123456789/25499
Weng, Z., Xiao, Z., Xie, Y., Wang, Z., and Gui, J. (2012). Genetic difference of Chinese horseshoe crab (Tachypleus tridentatus) in southeast coast of China based on mitochondrial COI gene analysis. Acta Oceanol. Sin. 31, 132–137. doi: 10.1007/s13131-012-0213-0
World Conservation Monitoring Centre (1996a). Tachypleus gigas, The IUCN Red List of Threatened Species 1996: e.T21308A9266907. Available online: https://doi.org/10.2305/IUCN.UK.1996.RLTS.T21308A9266907.en
World Conservation Monitoring Centre (1996b). Tachypleus Tridentatus, The IUCN Red List of Threatened Species 1996: e.T21309A9267047. Available online: https://doi.org/10.2305/IUCN.UK.1996.RLTS.T21309A9267047.en
World Conservation Monitoring Centre (1996c). Carcinoscorpius rotundicauda, The IUCN Red List of Threatened Species 1996: e.T3856A10123044. Available online at: https://doi.org/10.2305/IUCN.UK.1996.RLTS.T3856A10123044.en
Xu, Q., Chen, F., Shin, P. K. S., Cheung, S. G., Chen, Y., and Ke, C. (2011). AFLP analysis of genetic variation among three natural populations of horseshoe crab Tachypleus tridentatus along Chinese coast. Chin. J. Oceanol. Limnol. 29, 284–289. doi: 10.1007/s00343-011-0066-y
Yang, K. C., and Ko, H. S. (2015). First record of tri-spine horseshoe crab, Tachypleus tridentatus (Merostomata: Xiphosurida: Limulidae) from Korean waters. Anim. Syst. Evol. Diversity 31, 42–45. doi: 10.5635/ASED.2015.31.1.042
Yang, M. C., Chen, C. A., Hsieh, H. L., and Chen, C. P. (2007). Population subdivision of the tri-spine horseshoe crab, Tachypleus tridentatus, in Taiwan Strait. Zool. Sci. 24, 219–224. doi: 10.2108/zsj.24.219
Yang, M.-C., Chen, C.-P., Hsieh, H.-L., Huang, H., and Chen, C. A. (2009a). “Phylogeography, demographic history, and reserves network of horseshoe crab, Tachypleus tridentatus, in the South and East China seaboards,” in Biology and Conservation of Horseshoe Crabs, eds J. T. Tanacredi, M. L. Botton, and D. R. Smith (New York, NY: Springer), 163–181.
Yang, M.-C., Chen, C.-P., Hsieh, H.-L., Huang, T.-S., Hsieh, H. J., Tsai, W.-S.W., et al. (2009b). Is a living fossil about to go locally extinct? No mitochondrial genetic variation in horseshoe crab juveniles Tachypleus tridentatus at Haomeiliao Nature Reserve, Budai, Taiwan. Zool. Stud. 48:737. Available online at: http://zoolstud.sinica.edu.tw/Journals/48.6/737.html
Yap, C. K., Chong, C. M., and Tan, S. G. (2011). Allozyme polymorphisms in horseshoe crabs, Carcinoscorpius rotundicauda, collected from polluted and unpolluted intertidal areas in Peninsular Malaysia. Environ. Monit. Assess. 174, 389–400. doi: 10.1007/s10661-010-1464-x
Zaldívar-Rae, J., Sapién-Silva, R. E., Rosales-Raya, M., and Brockmann, H. J. (2009). “American horseshoe crabs, Limulus polyphemus, in Mexico: open possibilities,” in Biology and Conservation of Horseshoe Crabs, eds J. T. Tanacredi, M. L. Botton, and D. R. Smith (New York, NY: Springer), 97–113.
Zaleha, K., Hazwani, I., Siti Hamidah, H., Kamaruzzaman, B. Y., and Jalal, K. C. A. (2011). Effect of salinity on egg hatching and early larvae of horseshoe crab Tachypleus gigas (Müller 1785) in laboratory culture. J. Appl. Sci. 11, 2620–2626. doi: 10.3923/jas.2011.2620.2626
Zaleha, K., Akbar John, B., Erni Atika, H., Kamaruzzaman, B. Y., and Jalal, K. C. A. (2012). Spawning and nesting behaviour of Tachypleus gigas along the east coast of Peninsular Malaysia. Int. J. Biol. 4, 102–111. doi: 10.5539/ijb.v4n2p102
Zhou, H., and Morton, B. (2004). The diets of juvenile horseshoe crabs, Tachypleus tridentatus and Carcinoscorpius rotundicauda (Xiphosura). from nursery beaches proposed for conservation in Hong Kong. J. Nat. History 38, 1915–1925. doi: 10.1080/0022293031000155377
Zumdahl, S. S., and Zumdahl, S. A. (2000). Chemistry, 5th Edn., Boston, NY: Houghton Mifflin Company.
Keywords: Carcinoscorpius, climate change, ecological niche modeling, conservation, marine protected area, Tachypleus, Xiphosura
Citation: Vestbo S, Obst M, Quevedo Fernandez FJ, Intanai I and Funch P (2018) Present and Potential Future Distributions of Asian Horseshoe Crabs Determine Areas for Conservation. Front. Mar. Sci. 5:164. doi: 10.3389/fmars.2018.00164
Received: 18 December 2017; Accepted: 24 April 2018;
Published: 14 May 2018.
Edited by:
Rochelle Diane Seitz, Virginia Institute of Marine Science, United StatesReviewed by:
Youji Wang, Shanghai Ocean University, ChinaFrancine Kershaw, Natural Resources Defense Council, United States
Copyright © 2018 Vestbo, Obst, Quevedo Fernandez, Intanai and Funch. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Stine Vestbo, vestbo@bios.au.dk