Do Marmorkrebs, Procambarus fallax f. virginalis, threaten freshwater Japanese ecosystems?

  • Zen Faulkes1Email author,

    Affiliated with

    • Teresa Patricia Feria1 and

      Affiliated with

      • Jesús Muñoz2, 3

        Affiliated with

        Aquatic Biosystems20128:13

        DOI: 10.1186/2046-9063-8-13

        Received: 25 April 2012

        Accepted: 11 June 2012

        Published: 27 June 2012



        One marbled crayfish, Marmorkrebs, Procambarus fallax f. virginalis (Hagen, 1870), was discovered in a natural ecosystem in Japan in 2006. Because Marmorkrebs are parthenogenetic, they could establish a population from only a single individual, and thus pose a risk for becoming established in Japan, as they have in other countries. There are two major reasons to be concerned about the possibility of Marmorkrebs establishing viable populations in Japan. First, Japan’s only endemic crayfish, Cambaroides japonicus (De Haan, 1841), lives throughout Hokkaido and is endangered. Introduced Marmorkrebs are potential competitors that could further threaten C. japonicus. Second, Marmorkrebs live in rice paddies in Madagascar and consume rice. Marmorkrebs populations could reduce rice yields in Japan.


        We created five models in MaxEnt of the potential distribution of Marmorkrebs in Japan. All models showed eastern Honshu, Shikoku and Kyushu contain suitable habitats for Marmorkrebs. Hokkaido, the main habitat for C. japonicus, contained much less suitable habitat in most models, but is where the only Marmorkrebs in Japan to date was found.


        Marmorkrebs appear to be capable of establishing populations in Japan if introduced. They appear to pose minimal threat to C. japonicus, but may negatively affect rice production.


        In 2006, an unusual crayfish was collected in a river near Sapporo, and brought to the Sapporo Salmon Museum [1]. It had distinctive marbled colours, and all the offspring from this individual were later found to be female. This indicated that this crayfish was Marmorkrebs, a parthenogenetic [2, 3] crayfish that has been provisionally identified as Procambarus fallax f. virginalis (Hagen, 1870) [4]. (Here, we refer to the parthenogenetic form as “Marmorkrebs” and the sexual form as “P. fallax”).

        Marmorkrebs are unusual crayfish in two regards: they are the only known obligate parthenogenetic crayfish [2, 3], and the only known populations in natural ecosystems are the result of human introductions. Marmorkrebs have been introduced and established populations in Madagascar [5, 6] and Europe [710]. The discovery of Marmorkrebs in Hokkaido was the first well-documented case of an individual living in a natural ecosystem in Asia. Other reports of Marmorkrebs in Asia [11] have never been documented in the scientific literature or mainstream press. The most likely source of the Sapporo Marmorkrebs is a release or escape of a pet crayfish, as Marmorkrebs are widely circulated among aquarium hobbyists [12]. Marmorkrebs are parthenogenetic, and therefore a single female can initiate a stable population, resulting in an unwanted non-indigenous species.

        In Europe, single individual Marmorkrebs were discovered years before established populations were discovered [11, 13, 14]. Indeed, the delay was so long that it was questioned whether Marmorkrebs could establish populations in Europe [15]. This was a reasonable hypothesis, given that not all introduced species establish populations [16]. As it happened, Marmorkrebs have established populations in Germany [7], indicating that there was a “lag phase” [17] of several years between discovery of single individuals and establishment of populations. Thus, the discovery of the Sapporo Marmorkrebs may be a precursor to finding established populations of Marmorkrebs in Japan.

        We made quantitative models to assess the potential distribution of Marmorkrebs in Japan. There are at least two reasons to make such a threat assessment. First, Japan has only one native crayfish species, Cambaroides japonicus (De Haan, 1841). The historic range of C. japonicus is Hokkaido and the northern regions of Honshu [18]. This species is endangered [19], due in part to the introduction of the North American crayfish species, Pacifastacus leniusculus[2024]. The introduction of another exotic crayfish species could harm the remaining C. japonicus populations. Second, rice farming in Japan is economically important. Marmorkrebs have damaged rice paddies in Madagascar [14], although the extent of damage is not clear [5]. Marmorkrebs could become an agricultural pest in Japan if they become established. Thus, these models may help guide monitoring efforts, policy, and public information campaigns that could prevent further introductions or limit the spread.

        There are three common approaches to modeling the potential distribution of an exotic species: 1) extrapolate from the distribution of native populations only to the region of interest; 2) extrapolate from introduced populations in other regions to the region of interest; 3) extrapolate from a combination of native and introduced populations to the region of interest. Each of these methods has pros and cons, and it is not advisable to put too much weight on any single model [25].

        The first, and probably the most common approach, for developing models of the potential distribution of an exotic species is to describe the native distribution of the species, and then extrapolate from the climatic variables associated with those regions to find other geographic areas that have similar climatic features [25, 26]. Strictly speaking, it is not possible to use this method with Marmorkrebs because there are no known populations in the wild. We used the native distribution of P. fallax as a proxy for the native distribution of Marmorkrebs. We do so recognizing that, because the biotic conditions of each past introduction differ in overlap, the final fundamental niche space will vary depending on what dataset of presences is used to train the models [25]. As found in other studies of other species [27], invasive Marmorkrebs have a greater climatic range than wild P. fallax. There are many cases where parthenogens have quite distinct abilities and distributions than their sexual progenitors. This problem could be amplified by variation among clonal genotypes, which could be very distinct subsets of the range of abilities in the native sexual range. This has not been studied on Marmorkrebs yet, and using data from sexual P. fallax is a pragmatic approach.

        Other approaches to modeling potential distribution are to use established populations outside the native range, and a combination of the native and introduced populations. We use all three approaches (similar to [28]). We also created two additional models that included the location of the one Marmorkrebs individual discovered in Sapporo. These models including the Sapporo Marmorkrebs are speculative, because single individuals may not indicate that there is an established or even viable population [25]. Nevertheless, if it becomes apparent that there are established populations Marmorkrebs in Japan rather than single individuals, these models may help adjust the assessment of risk accordingly.


        According to the complexity of the models, measured as the number of parameters and lower AICc values (Table 1), and the visual inspection of response variables (Figure 1), we based all further analyses on models with a regularization multiplier of 1.5 and the reduced dataset including seven uncorrelated bioclimatic variables. The multivariate environmental similarity surfaces (MESS) showed that only a fraction of central Kyushu has novel conditions compared to the training dataset (hatched areas in Figures 2, 3, 4, 5 and 6), and thus only projections in that reduced area should be considered with caution. Procambarus fallax wild populations and Marmorkrebs have different environmental requirements, as shown by a Mann–Whitney U test (Table 2) and the response curves in Figure 1.
        Table 1

        Overview of MaxEnt models

        Trained on

        Sample size




        AICc score

        Training AUC

        Test AUC (30%)

        Procambarus fallax



        MTSPS = 0.284





        Procambarus fallax



        MTP = 0.284




        0.9766 (0.0043)

        Marmorkrebs populations



        MTSPS = 0.033


        not applicable



        Marmorkrebs populations



        MTP = 0.014




        0.9823 (0.0146)

        Marmorkrebs populations and Sapporo individual



        MTSPS = 0.045


        not applicable



        Marmorkrebs populations and Sapporo individual



        MTP = 0.030


        not applicable


        0.9723 (0.0212)

        P. fallax and Marmorkrebs populations



        MTSPS = 0.328





        P. fallax and Marmorkrebs populations



        MTP = 0.027




        0.9595 (0.0119)

        P. fallax, Marmorkrebs populations, and Sapporo individual



        MTSPS = 0.321





        P. fallax, Marmorkrebs populations, and Sapporo individual



        MTP = 0.020




        0.9537 (0.0148)

        Training AUC is the AUC obtained in the 1-replicate run used for model selection and MESS generation. Test AUC is the average AUC (and standard deviation) obtained in 10 runs that used a 30% of the presence points to calculate this statistic; only models considered superior based on reduced complexity without training AUC reduction were projected and test AUC calculated. MTSPS – minimum training sensibility plus specificity; MTP – minimum training presence; not applicable – AICc cannot be calculated if the model have more parameters than occurrence points [29].
        Figure 1

        Response curves of MaxEnt models. Response curves of MaxEnt models created using only the corresponding variable for (A) P. fallax; (B) Marmorkrebs in Madagascar and Germany; (C) Marmorkrebs in Madagascar, Germany, and Sapporo; (D) P. fallax and Marmorkrebs in Madagascar and Germany (E) P. fallax and Marmorkrebs in Madagascar, Germany, and Sapporo. Note that the variables used in this study are not correlated, and thus the graphs reflect the dependence of predicted suitability on the selected variable. The inclusion of the individual Marmorkrebs found in Sapporo has virtually no effect on the response curves (compare column B to C, or column D to E). Wild P. fallax and parthenogenetic Marmorkrebs show different optimal responses to variables, except for mean temperature of driest quarter. Such different responses would explain the bimodal curves of mean temperature of wettest quarter, the hump-shaped curves of precipitation seasonality or precipitation of wettest quarter, or the asymmetrically humped curves of precipitation of driest quarter in rows D and E.
        Figure 2

        Marmorkrebs distribution in Japan modeled using P. fallax. Potential distribution of Marmorkrebs in Japan as predicted by a model trained using native P. fallax distribution in United States. Legend in this and subsequent figures: Gray scale indicates similarity between known and predicted habitats (i.e., habitat suitability). Circles show presence of Cambaroides japonicus (see Table 3); square shows location of single Marmorkrebs found in Sapporo in 2006 [1]. Regions with red hatching in this and subsequent figures show areas with novel conditions compared to the training dataset according to multivariate environmental similarity surfaces (MESS) analysis. Inset in this and subsequent figures: Choropleth map of rice production of agricultural regions [30].
        Figure 3

        Marmorkrebs distribution in Japan modeled using introduced Marmorkrebs. Potential distribution of Marmorkrebs in Japan as predicted by a model trained using distributions of introduced Marmorkrebs populations in Madagascar and Europe.
        Figure 4

        Marmorkrebs distribution in Japan modeled using introduced Marmorkrebs and Sapporo Marmorkrebs. Potential distribution of Marmorkrebs in Japan as predicted by a model trained using introduced Marmorkrebs populations in Madagascar and Europe, and the single individual Marmorkrebs found in Hokkaido.
        Figure 5

        Marmorkrebs distribution in Japan modeled using P. fallax and introduced Marmorkrebs. Potential distribution of Marmorkrebs in Japan as predicted by a model trained using native P. fallax distribution in United States and introduced Marmorkrebs populations in Madagascar and Europe.
        Figure 6

        Marmorkrebs distribution in Japan modeled using P. fallax , introduced Marmorkrebs, and Sapporo Marmorkrebs. Potential distribution of Marmorkrebs in Japan as predicted by a model trained using native P. fallax distribution in United States, introduced Marmorkrebs populations in Madagascar and Europe, and the single individual Marmorkrebs found in Sapporo.

        Table 2

        Mean and standard deviation of bioclimatic variables

        Climatic variable

        Mean (s.d.) ofP. fallax, Florida populations

        Mean (s.d.) of Marmorkrebs populations, excluding Sapporo individual

        Mean difference

        Utest sig. (2-tailed)

        Mean diurnal range

        121.62 (11.994)

        107.78 (11.969)



        Temperature seasonality

        4799.3 (729.823)

        3092.17 (1600.579)



        Mean temperature of wettest quarter

        270.05 (2.697)

        206.61 (13.836)



        Mean temperature of driest quarter

        173.27 (14.327)

        131.17 (53.102)



        Precipitation seasonality

        45.27 (8.549)

        86.7 (28.919)



        Precipitation of wettest quarter

        540.62 (44.529)

        720.83 (233.981)



        Precipitation of driest quarter

        181.16 (24.443)

        48.43 (43.833)



        All models predict southeastern Japan – specifically, Shikoku, Kyushu and eastern Honshu – is suitable habitat for Marmorkrebs (Figures 2, 3, 4, 5 and 6). Models trained with P. fallax predict that Hokkaido, the main habitat for C. japonicus, is unlikely to contain suitable habitat for Marmorkrebs (Figure 2), but models that included the Marmorkrebs populations (Figures 3, 4, 5 and 6) predict that Hokkaido also contains suitable habitat. All models had high average performance scores, with an averaged Area Under the Curve (AUC) in Receiver Operating Characteristics (ROC) plots over 10 runs greater than 0.95 (Table 1).

        Model 1: Trained on P. fallax populations

        Marginal regions of eastern Shikoku and Kyushu, and south-central Honshu are likely to contain suitable habitat for Marmorkrebs under this model (Figure 2). This model predicted Japan provides less suitable habitat for Marmorkrebs than all other models. Similarly, the P. fallax-trained model of potential distribution of Marmorkrebs in Madagascar also predicted much less habitat than other models, and did not predict some regions where Marmorkrebs occurred [28]. The most important variables in this model were temperature of the wettest quarter, with a 46.1% contribution, and precipitation seasonality, with a 23.3% contribution.

        Model 2: Trained on Marmorkrebs populations

        Under this model, most of Shikoku is suitable for Marmorkrebs, as well as a substantial part of mainland Honshu, especially to the northeastern half north of Shizuoka Prefecture (Figure 3). Kyushu and Hokkaido show less area as suitable, mainly on the coastal regions.

        The climatic variables with the greatest predictive power were the same as in the previous model, although switched. The variables in the present study cannot be compared to the equivalent model in [28], because here we have removed correlated variables.

        Model 3: Trained on Marmorkrebs populations and Sapporo Marmorkrebs

        This model had the greatest variation in its predictive performance. This is probably due to most Marmorkrebs presences being located in one small geographic region in Madagascar, with a small number of outliers in Germany and Sapporo. Thus, it is expected that replicate tests of model performance would be the most susceptible to which data points are included or excluded in the training and testing phase (see Methods). The variable importance is virtually identical to the previous model, which is expected given that there is only one more presence in the training dataset.

        After applying the threshold rule, the regions predicted as suitable for the establishment of Marmorkrebs by this model are basically the same as predicted by the latter (Figure 1c): most of Shikoku, northern Kyushu, and a crescent area covering the east part of Honshu from southwest to northeast, although slightly displaced to the south. The addition of the data point from Hokkaido completely switches the area predicted present in this island to the western coastal regions, while the suitable area in the east disappears.

        Model 4: Trained on P. fallax and Marmorkrebs populations

        This model is similar to models 2–3, although the areas predicted as suitable are almost continuous across Kyushu, Shikoku, and the south half of Honshu, and have virtually disappeared from Hokkaido, except for a limited region in the southern central coast (Figure 5). The inclusion of Marmorkrebs populations from Madagascar and Germany has a dramatic effect, increasing the potential area in the southern Japan islands.

        The variables now more important in model construction are precipitation seasonality, as in all previous models, and mean temperature of driest quarter. The change in the variable importance can be an artefact due to the low number of Marmorkrebs localities known, and that most of the Madagascar localities are closely clustered and thus present very similar climate. Although variables are not correlated with regard to background area, they still show a strong relationship if only the Marmorkrebs localities are considered.

        Model 5: Trained on P. fallax, Marmorkrebs populations, and Sapporo Marmorkrebs

        This model resembles Model 3 in that the addition of the Sapporo Marmorkrebs resulted in more of Hokkaido being predicted as suitable habitat (Figure 1e). However, for the rest of the study area, this model is closer to Model 4, predicting much of Kyushu, Shikoku and south Honshu as suitable for Marmorkrebs. This suggests that the combination of both the native P. fallax populations and the stable Marmorkrebs populations generates a more useful model by extrapolating both the fundamental climate space of the wild P. fallax, which is conservative, and the expanded one of Marmorkrebs, which successfully established viable populations in regions with a broader climate space.


        Species distribution models predict that large regions of Japan are suitable habitat for Marmorkrebs. Although the most conservative model (Model 1, trained on P. fallax data) predicts only a few spots in Japan as suitable, the addition of established populations of Marmorkrebs highly increased the extension across Japan. The northern island of Hokkaido contains the least suitable habitat under models using existing populations; however, the addition of the temperate localities of Europe largely increased the area of potential presence in this island. These results indicate that Marmorkrebs pose a threat to freshwater ecosystems in Japan.

        Marmorkrebs may pose some threat to Japan’s endangered native crayfish, C. japonicus, through direct competition. Model 2 shows suitable habitat overlap with Marmorkrebs' known localities in eastern Hokkaido, Model 3 overlaps around Sapporo, and in Model 5 only central Hokkaido is not predicted as suitable, but most C. japonicus localities coincide with high suitability areas for Marmorkrebs. Nevertheless, there is evidence that C. japonicus prefers deep [31] and cold [21] water, whereas Marmorkrebs appears to prefer shallow, warm waters [7]. Thus, C. japonicus probably faces greater competition from another cold water crayfish species, Pacifastacus leniusculus, than from Marmorkrebs [2024]. Marmorkrebs may pose a threat to C. japonicus if Marmorkrebs carry crayfish plague. Cambaroides japonicus is susceptible to crayfish plague [32], but it is not clear if crayfish plague is widespread in Japan or if it has been a factor in the decline of C. japonicus[22]. There is no direct evidence that Marmorkrebs carry crayfish plague, but it is plausible that Marmorkrebs could act as a vector for the disease because every North American species tested has the ability to act as a chronic carrier for the causative agent of crayfish plague, the water mold Aphanomyces astaci[3234].

        Rice farming appears to be a more likely impact of a Marmorkrebs introduction. The eastern coastal areas of the most productive rice region, the Tohoku agricultural region in northernmost Honshu, and most of the Kanto-Tosan region, the next most productive, are consistently predicted to be highly likely to support Marmorkrebs populations (Figures 2, 3, 4, 5 and 6).

        As with models of potential habitat in Madagascar, North American, and Europe of other geographic regions [28], the model trained on P. fallax data makes the most conservative predictions of potential distribution. This can be due to the different requirements, or wider habitat breadth observed in Marmorkrebs (Table 2). Models that include Marmorkrebs data predict larger areas as potential habitat. There are several points to keep in mind in interpreting these models. The P. fallax-trained model was the most conservative in predicting the distribution of introduced Marmorkrebs, and vice versa [28], similar to what was found for the native and introduced distributions of signal crayfish in Japan [35]. Both studies interpreted these imperfect predictions indicating a niche shift for these crayfish species. The models predict lower habitat suitability in northern Japan, suggesting that the models pass the “eco-plausibility test” [25].

        There are two potential limitations to the models presented here. First, the locations of P. fallax have uncertainty associated with them, because the locations were estimated from records at the level of counties. The abundance and apparently uniform presence of P. fallax throughout its range, occupying “practically all of the streams” in its range [36], however, mitigates against the uncertainty inherent in using county level records. Second, terrestrial climatic variables were used to train the models rather than variables describing the aquatic habitats. Temperature and precipitation are the two major categories of climatic variables often used in species distribution models. Terrestrial air temperature correlates with lake water temperature [37], and water temperature affects growth and mortality in Marmorkrebs [38]. Similarly, precipitation levels (i.e., drought) affect wild populations of P. fallax[39, 40]. Thus, there are good reasons to expect that terrestrial climate variables can predict a substantial amount about the distribution of aquatic species; indeed, several other studies have done so [28, 35, 4143]. Further research is needed to determine if mismatches in predicted and actual distributions of aquatic species (e.g., signal crayfish [35]) might be explained by variables about aquatic habitats being omitted from models. Currently, however, there is no database for features of freshwater bodies equivalent to that for terrestrial climate.


        Controlling crayfish populations in other regions, such as Europe, has proved to be extraordinarily difficult. Physical barriers, such as dams [44], may hinder the spread of introduced crayfish, but there are no control agents that can target and remove crayfish once established. Drastic non-targeted measures, such as applying biocides to entire water bodies [45], have had mixed success. As the pet trade is the most likely source of Marmorkrebs introductions, education of pet owners and promoting responsible ownership may be a cost effective measure to reduce the chance of Marmorkrebs becoming established in Japan.


        Species distribution data

        We used the distributions of established populations of Procambarus fallax in the United States and established populations of Marmorkrebs in Madagascar that were used previously [28] (spreadsheet of locations used in that model available at http://​marmorkrebs.​org). The main difference in this dataset and that used in [28] was that we removed locations in Europe where only single Marmorkrebs had been reported [11, 15, 46, 47] and added two established Marmorkrebs populations that were subsequently found in Germany [79] and one new population found in Madagascar [Lake Alaotra; J.P.G. Jones, personal communication]. In all cases, the locations used to train the models represent populations of crayfish, not single individuals. Single individuals that are occasionally found outside of the normal range, but do not survive there, can distort models [25].

        The distribution of P. fallax was taken from previously published records of its presence in American counties [36, 48]. County names were searched though Google Earth v. 5.2.1. (http://​www.​google.​com/​earth/​), which provided latitude and longitude measurements near the center of each county. Although the available location data may be coarse, Hobbs noted that P. fallax is abundant throughout most of its distribution and could be found in “practically all of the streams” in its range [36].

        The distribution of Marmorkrebs in Madagascar was estimated from maps and published GPS coordinates [5, 6]. Published maps were superimposed on locations in Google Earth with the “Add Image Overlay” tool. The transparency of the map was adjusted so that both the published map and the Google Earth image simultaneously. Placemarks on Google Earth were put onto points indicated on the published map, from which latitude and longitude were recorded.

        The distribution of Marmorkrebs in Europe were taken from published GPS coordinates [7] or estimated from city locations [8, 9] using Google Earth.

        To visualize the potential overlap of Marmorkrebs and C. japonicus distributions, we searched the scientific literature for records of specific locations of C. japonicus (Table 3). In cases where GPS coordinates were not given, locations were estimated by nearest reference to a city or other geographic feature using Google Earth.
        Table 3

        Locations of Cambaroides japonicus described in the scientific literature








        Hakodate, Hokkaido





        Kushiro District, Hokkaido





        Lake Komadome, Hokkaido





        Otaru, Hokkaido





        Hamamasu, Hokkaido





        Atsuta, Hokkaido





        Shikotsu Lake, Hokkaido





        Lake Shikaribetsu, Hokkaido a





        Kikanko River, Obihiro City, Hokkaido a





        Aomori, Honshu





        Kusshara Lake, Hokkaido





        Otaru City, Hokkaido





        Oshoppu stream in Abashiri City, Hokkaido





        Sapporo, Hokkaido




        Rice production

        Data on paddy field rice production for “agricultural regions” was taken from [30]. Each agricultural region is composed of one or more prefectures. We used GIS coordinates for prefectures from [59], with each prefecture being assigned the average rice production for its agricultural region.

        Modeling process

        The modeling process was similar to that used in [28]. Models are trained using the existing distribution of the species of interest. The distribution of the species is matched to bioclimatic variables, and potential distributions are then predicted which finds other regions with bioclimatic conditions. Terrestrial climatic variables and those of non-ocean aquatic habitats are correlated in many ways. For example, terrestrial air temperature correlates well with water temperature in lakes [37]. Several studies have used terrestrial climatic variables to model the broad-scale distribution of freshwater aquatic species across terrestrial landscapes [28, 35, 4143].

        We created models using MaxEnt v. 3.3.3e [26, 60] because it discriminates suitable and unsuitable areas better than other methods [26, 6163]. While some parameters were set to defaults (convergence threshold = 10–5, or maximum iterations = 500), we used two regularization multiplier values (1.0 and 1.5). A regularization multiplier of 1.5 allows that variables’ average values in the projections spread from the empirical average of the background points (the situation if it is set to 1.0), thus avoiding model overfitting [29]. The 10,000 background points were randomly selected from the areas of the immediate adjacency of presences and pertaining to the same bioclimatic region according to the Köppen-Geiger classification scheme [6466]. To test different training strategies, we created five presence datasets: 1) the native distribution of P. fallax in the United States; 2) the established populations of Marmorkrebs in Madagascar and Germany; 3) the established populations of Marmorkrebs in Madagascar and Germany, plus the Sapporo Marmorkrebs; 4) a combination of the American distribution of P. fallax and Marmorkrebs established populations in Madagascar and Germany; and 5) the American distribution of P. fallax, the Marmorkrebs populations in Madagascar and Germany, plus the Sapporo Marmorkrebs.

        The nineteen Worldclim bioclimatic variables (http://​www.​worldclim.​org) were intersected with the presences and background datasets. To avoid multicollinearity and model overfitting, we ran a correlation analysis on the background dataset and eliminated one of the variables in each pair with a Pearson correlation value >0.8. The final datasets (Table 2) included mean diurnal range (mean of monthly (max temp - min temp)); temperature seasonality (standard deviation × 100); mean temperature of wettest quarter; mean temperature of driest quarter; precipitation seasonality (coefficient of variation); precipitation of wettest quarter, and; precipitation of driest quarter. These bioclimatic variables result from global land area interpolation of climate point data (years 1950–2000) at a spatial resolution of 30 arc-sec [67] which corresponds to a spatial resolution of approximately 0.008333° at the Equator (1 × 1 km2 grid cells).

        To confirm that combinations of novel climates were not cause of concern in the projections, we follow [68, 69] to generate the multivariate environmental similarity surfaces (MESS); this grid was reclassified and values below zero were masked to show areas of novel climate space relative to the range under which the model was fitted.

        To decide the combination of parameters that would be projected onto the final model, we used the small sample corrected Akaike’s Information Criterion (AICc) implemented in ENMTools [29]. The less complex model, i.e. lower AICc, was projected to the globe at 5 km pixel size, and to Japan at 1 km pixel size. To generate the presence/absence grids, we used the “minimum training presence” threshold, as in similar studies on invasive species [25]. For further testing of model transferability, we also checked the shape of the response curves looking for open-ended, unrealistic responses (Figure 1).

        To evaluate the predictive ability of these models, ten replicates of each model were run. In each replicate, 70% of the locations with known occurrences were selected randomly to train the model. The final models were in all cases the average of their 10 replicates for each set of training data; and the threshold used to create the binary grid the average threshold over the runs. The accuracy of each model was estimated using the Area Under the Curve (AUC) in Receiver Operating Characteristics (ROC) plots [70], calculated in MaxEnt on the 30% unused presences. ROC is a threshold–independent measure that evaluates the probability that the model produces a positive result in a positive locality (sensitivity) versus the probability that the model produces a negative result in a negative locality (specificity) when presented with new data. A ROC plot is obtained by plotting all sensitivity values on the y–axis against their equivalent (1–specificity) values for all available decision thresholds on the x–axis. An AUC score of 1 indicates perfect performance; a score of 0.5 indicates random performance.


        Authors’ Affiliations

        Department of Biology, The University of Texas-Pan American
        Real Jardín Botánico (RJB-CSIC)
        Universidad Tecnológica Indoamérica


        1. Kawai T, Takahata M (Eds): Biology of crayfish. Hokkaido University Press, Sapporo; 2010.
        2. Scholtz G, Braband A, Tolley L, Reimann A, Mittmann B, Lukhaup C, Steuerwald F, Vogt G: Parthenogenesis in an outsider crayfish. Nature 2003, 421:806–806.View Article
        3. Vogt G, Tolley L, Scholtz G: Life stages and reproductive components of the Marmorkrebs (marbled crayfish), the first parthenogenetic decapod crustacean. J Morphol 2004, 261:286–311.View Article
        4. Martin P, Dorn NJ, Kawai T, van der Heiden C, Scholtz G: The enigmatic Marmorkrebs (marbled crayfish) is the parthenogenetic form of Procambarus fallax (Hagen, 1870). Contrib Zool 2010, 79:107–118.
        5. Jones JPG, Rasamy JR, Harvey A, Toon A, Oidtmann B, Randrianarison MH, Raminosoa N, Ravoahangimalala OR: The perfect invader: a parthenogenic crayfish poses a new threat to Madagascar’s freshwater biodiversity. Biol Invasions 2009, 11:1475–1482.View Article
        6. Kawai T, Scholtz G, Morioka S, Ramanamandimby F, Lukhaup C, Hanamura Y: Parthenogenetic alien crayfish (Decapoda: cambaridae) spreading in Madagascar. J Crust Biol 2009, 29:562–567.View Article
        7. Chucholl C, Pfeiffer M: First evidence for an established Marmorkrebs (Decapoda, Astacida, Cambaridae) population in Southwestern Germany, in syntopic occurrence with Orconectes limosus (Rafinesque, 1817). Aquat Invasions 2010, 5:405–412.View Article
        8. Löwe K: Gefahr aus dem Aquarium. Halle: Mitteldeutsches Druck- und Verlagshaus GmbH & Co. KG, Mitteldeutsche Zeitung; 2010. http://​www.​mz-web.​de/​servlet/​ContentServer?​pagename=​ksta/​page&​atype=​ksArtikel&​aid=​1286541137817&​openMenu=​1013016724285&​calledPageId=​1013016724285&​listid=​1018881578312
        9. Privenau K: Marmorkrebs bringt Pest. Halle: Mitteldeutsches Druck- und Verlagshaus GmbH & Co. KG, Mitteldeutsche Zeitung; 2010. http://​www.​mz-web.​de/​servlet/​ContentServer?​pagename=​ksta/​page&​atype=​ksArtikel&​aid=​1286541132341&​calledPageId=​987490165154
        10. Stloukal E: Recent distribution of non-indigenuous (sic) crayfish species in Slovakia. Folia Faunistica Slovaca 2009, 14:119–122.
        11. Marzano FN, Scalici M, Chiesa S, Gherardi F, Piccinini A, Gibertini G: The first record of the marbled crayfish adds further threats to fresh waters in Italy. Aquat Invasions 2009, 4:401–404.View Article
        12. Faulkes Z: The spread of the parthenogenetic marbled crayfish, Marmorkrebs (Procambarus sp.), in the North American pet trade. Aquat Invasions 2010, 5:447–450.View Article
        13. Holdich DM, Reynolds JD, Souty-Grosset C, Sibley PJ: A review of the ever increasing threat to European crayfish from non-indigenous crayfish species. Knowl Manag Aquat Ecosyst 2009, 394–395:11.View Article
        14. Invasion of self-cloning crayfish alarms Madagascar. http://​www.​earthtimes.​org/​articles/​news/​339974,alarms-madagascar-feature.​html
        15. Martin P, Shen H, Füllner G, Scholtz G: The first record of the parthenogenetic Marmorkrebs (Decapoda, Astacida, Cambaridae) in the wild in Saxony (Germany) raises the question of its actual threat to European freshwater ecosystems. Aquat Invasions 2010, 5:397–403.View Article
        16. Williamson MH, Fitter A: The characters of successful invaders. Biol Conserv 1996, 78:163–170.View Article
        17. Theoharides KA, Dukes JS: Plant invasion across space and time: factors affecting nonindigenous species success during four stages of invasion. New Phytol 2007, 176:256–273.View Article
        18. Okada Y: Some observations of Japanese crayfishes. Sci Rep Tokyo Bunrika Daigaku B 1933, 14:155–188.
        19. Usio N: Endangered crayfish in northern Japan: distribution, abundance and microhabitat specificity in relation to stream and riparian environment. Biol Conserv 2007, 134:517–526.View Article
        20. Nakata K, Goshima S: Asymmetry in mutual predation between the endangered Japanese native crayfish Cambaroides japonicus and the North American invasive crayfish Pacifastacus leniusculus: a possible reason for species replacement. J Crust Biol 2006, 26:134–140.View Article
        21. Nakata K, Hamano T, Hayashi K-I, Kawai T: Lethal limits of high temperature for two crayfishes, the native species Cambaroides japonicus and the alien species Pacifastacus leniusculus in Japan. Fish Sci 2002, 68:763–767.View Article
        22. Nakata K, Goshima S: Competition for shelter of preferred sizes between the native crayfish species Cambaroides japonicus and the alien crayfish species Pacifastacus leniusculus in Japan in relation to prior residence, sex difference, and body size. J Crust Biol 2003, 23:897–907.View Article
        23. Usio N, Konishi M, Nakano S: Species displacement between an introduced and a ‘vulnerable’ crayfish: the role of aggressive interactions and shelter competition. Biol Invasions 2001, 3:179–185.View Article
        24. Usio N, Suzuki K, Konishi M, Nakano S: Alien vs. endemic crayfish: roles of species identity in ecosystem functioning. Arch Hydrobiol 2006, 166:1–21.View Article
        25. Rodda GH, Jarnevich CS, Reed RN: Challenges in identifying sites climatically matched to the native ranges of animal invaders. PLoS One 2011, 6:e14670.View Article
        26. Phillips SJ, Anderson RP, Schapire RE: Maximum entropy modeling of species geographic distributions. Ecol Model 2006, 190:231–259.View Article
        27. Loo SE, Nally RM, Lake P: Forecasting New Zealand mudsnail invasion range: model comparisons using native and invaded ranges. Ecol Appl 2007, 17:181–189.View Article
        28. Feria TP, Faulkes Z: Forecasting the distribution of Marmorkrebs, a parthenogenetic crayfish with high invasive potential, in Madagascar, Europe, and North America. Aquat Invasions 2011, 6:55–67.View Article
        29. Warren DL, Seifert SN: Ecological niche modeling in Maxent: the importance of model complexity and the performance of model selection criteria. Ecol Appl 2011, 21:335–342.View Article
        30. Ministry of Agriculture, Forestry and Fisheries: Appendix 3: Rice Production. Mon Stat Agric For Fish 2011, 694:92.
        31. Koga T, Kawai T: Water depth as a factor influencing growth rate and size of the Japanese endemic crayfish Cambaroides japonicus (Crustacea: Decapoda). Jpn J Benthol 2003, 58:11–15.
        32. Unestam T: On the host range and origin of the crayfish plague fungus. Rep Inst Freshw Res Drottningholm 1972, 52:192–198.
        33. Kozubíková E, Filipová L, Kozák P, Ďuriš Z, Martín MP, Diéguez-Uribeondo J, Oidtmann B, Petrusek A: Prevalence of the crayfish plague pathogen Aphanomyces astaci in invasive American crayfishes in the Czech Republic. Conserv Biol 2009, 23:1204–1213.View Article
        34. Aquiloni L, Martín M, Gherardi F, Diéguez-Uribeondo J: The North American crayfish Procambarus clarkii is the carrier of the oomycete Aphanomyces astaci in Italy. Biol Invasions 2011, 13:359–367.View Article
        35. Larson ER, Olden JD, Usio N: Decoupled conservatism of Grinnellian and Eltonian niches in an invasive arthropod. Ecosphere 2010, 1:art16.View Article
        36. Hobbs HH Jr: The crayfishes of Florida. Univ Fl Publ Biol Ser 1942, 3:1–179.
        37. Livingstone DM, Lotter AF: The relationship between air and water temperatures in lakes of the Swiss Plateau: a case study with palæolimnological implications. J Paleolimnol 1998, 19:181–198.View Article
        38. Seitz R, Vilpoux K, Hopp U, Harzsch S, Maier G: Ontogeny of the Marmorkrebs (marbled crayfish): a parthenogenetic crayfish with unknown origin and phylogenetic position. J Exp Zool A 2005, 303:393–405.View Article
        39. Dorn N, Trexler JC: Crayfish assemblage shifts in a large drought-prone wetland: the roles of hydrology and competition. Freshwater Biol 2007, 52:2399–2411.View Article
        40. Dorn N, Volin JC: Resistance of crayfish (Procambarus spp.) populations to wetland drying depends on species and substrate. J N Am Benthol Soc 2009, 28:766–777.View Article
        41. Kumar S, Spaulding SA, Stohlgren TJ, Hermann KA, Schmidt TS, Bahls LL: Potential habitat distribution for the freshwater diatom Didymosphenia geminata in the continental US. Front Ecol Environ 2008, 7:415–420.View Article
        42. Januchowski-Hartley S, VanDerWal J, Sydes D: Effective control of aquatic invasive species in tropical Australia. Environ Manag 2011, 1–9.
        43. Reshetnikov A, Ficetola G: Potential range of the invasive fish rotan (Perccottus glenii) in the Holarctic. Biol Invasions 2011.
        44. Dana ED, García-de-Lomas J, González R, Ortega F: Effectiveness of dam construction to contain the invasive crayfish Procambarus clarkii in a Mediterranean mountain stream. Ecol Eng 2011, 37:1607–1613.View Article
        45. Peay S, Hiley PD, Collen P, Martin I: Biocide treatment of ponds in Scotland to eradicate signal crayfish. Knowl Manag Aquat Ecosyst 2006, 380–381:1363–1379.View Article
        46. Blanke D, Schulz H: Situation des Edelkrebses (Astacus astacus L.) sowie weiterer Flusskrebsarten in Niedersachsen. In Tagungsbericht der Deutschen Gesellschaft für Limnologie (DGL), Braunschweig, 2002. DGL, Werder; 2003:385–389.
        47. Holdich DM, Pöckl M: Invasive crustaceans in European inland waters. In Biological Invaders in Inland Waters: Profiles, Distribution, and Threats. Edited by: Gherardi F. Springer, Dordrecht, The Netherlands; 2007:29–75. [Drake JA (Series Editor): Invading Nature: Springer Series in Invasion Ecology]View Article
        48. Hendrix AN, Loftus WF: Distribution and relative abundance of the crayfishes Procambarus alleni (Faxon) and P. fallax (Hagen) in southern Florida. Wetlands 2000, 20:194–199.View Article
        49. Andrews EA: Sperm-transfer organs in Cambaroides. Biol Bull 1909, 17:257–270.View Article
        50. Ichijo N, Inoue M, Ando T, Sasaki M, Arihara M, Yoshida K: Distribution of the Japanese crayfish (Cambaroides japonicus) and signal crayfish (Pacifastacus leniusculus), in the Kushiro District of Eastern Hokkaido, Japan –New challenge by high school students. Jpn J Conserv Ecol 2009, 14:243–248.
        51. Kawai T, Saito K: Observations on the mating behavior and season, with no form alternation, of the Japanese crayfish, Cambaroides japonicus (Decapoda, Cambaridae), in Lake Komadome, Japan. J Crust Biol 2001, 21:885–890.View Article
        52. Kawai T, Scholtz G: Behavior of juveniles of the Japanese endemic species Cambaroides japonicus (Decapoda: Astacidea: Cambaridae), with observations on the position of the spermatophore attachment on adult females. J Crust Biol 2002, 22:532–537.View Article
        53. Kawai T, Hamano T, Matsuura S: Survival and growth of the Japanese crayfish Cambaroides japonicus in a small stream in Hokkaido. Bull Mar Sci 1997, 61:147–157.
        54. Kawai T, Nakata K, Suzuki Y: Decrease of the habitat number of Japanese crayfish Cambaroides japonicus around Sapporo City. Bull Sapporo Salmon Mus 2001, 21–26.
        55. Niiyama H: The chromosomes of the crayfish, Cambaroides japonicus (De Haan). J Fac Sci Hokkaido Imp Univ Ser V I Zool 1934, 3:41–53.
        56. Scholtz G, Kawai T: Aspects of embryonic and postembryonic development of the Japanese freshwater crayfish Cambaroides japonicus (Crustacea, Decapoda) including a hypothesis on the evolution of maternal care in the Astacida. Acta Zool 2002, 83:203–212.View Article
        57. Yamanaka K, Kuwabara R, Shio T: Larval development of a Japanese crayfish, Cambaroides japonicus (De Haan). Bull Mar Sci 1997, 61:165–175.
        58. Yoshino M, Hisada M: The effects of some putative transmitters and biogenic amines on uropod abductor muscle in the crayfish Procambarus clarkii and Cambaroides japonicus. Comp Biochem Phys C 1984, 78:289–293.View Article
        59. Japan GIS Datasets V2. http://​www.​fas.​harvard.​edu/​~chgis/​japan/​datasets.​html
        60. Phillips SJ, Dudík M, Schapire RE: A maximum entropy approach to species distribution modeling. In Proceedings of the 21st International Conference on Machine Learning; Banff. Edited by: Greiner R, Schuurmans D. ACM Press, ; 2004:655–662.
        61. Phillips SJ: Transferability, sample selection bias and background data in presence-only modelling: a response to Peterson et al. (2007). Ecography 2008, 31:272–278.View Article
        62. Elith J, Graham C, Anderson RP, Dudik M, Ferrier S, Guisan A, Hijmans RJ, Huettmann F, Leathwick JR, Lehmann A, et al.: Novel methods improve prediction of species’ distributions from occurrence data. Ecography 2006, 29:129–151.View Article
        63. Hernandez PA, Graham CH, Master LL, Albert DL: The effect of sample size and species characteristics on the performance of different species distribution modelling methods. Ecography 2006, 29:773–785.View Article
        64. Köppen WP: Das geographische system der klimate. Gebrüder Borntraeger, Berlin; 1936.
        65. Kriticos DJ, Webber BL, Leriche A, Ota N, Macadam I, Bathols J, Scott JK: CliMond: global high-resolution historical and future scenario climate surfaces for bioclimatic modelling. Methods Ecol Evol 2011, 3:53–64.View Article
        66. Webber BL, Yates CJ, Le Maitre DC, Scott JK, Kriticos DJ, Ota N, McNeill A, Le Roux JJ, Midgley GF: Modelling horses for novel climate courses: insights from projecting potential distributions of native and alien Australian acacias with correlative and mechanistic models. Divers Distrib 2011, 17:978–1000.View Article
        67. Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A: Very high resolution interpolated climate surfaces for global land areas. Int J Climatol 2005, 25:1965–1978.View Article
        68. Elith J, Kearney M, Phillips S: The art of modelling range shifting species. Methods Ecol Evol 2010, 1:330–642.View Article
        69. Elith J, Phillips SJ, Hastie T, Dudík M, Chee YE, Yates CJ: A statistical explanation of MaxEnt for ecologists. Divers Distrib 2011, 17:43–57.View Article
        70. Fielding AH, Bell JF: A review of methods for the assessment of prediction errors in conservation presence/absence models. Environ Conserv 1997, 24:38–49.View Article


        © Faulkes et al.; licensee BioMed Central Ltd. 2012

        This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.