- Open Access
Do Marmorkrebs, Procambarus fallax f. virginalis, threaten freshwater Japanese ecosystems?
Aquatic Biosystemsvolume 8, Article number: 13 (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 . 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) . (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 [7–10]. 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  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 . 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 . This was a reasonable hypothesis, given that not all introduced species establish populations . As it happened, Marmorkrebs have established populations in Germany , indicating that there was a “lag phase”  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 . This species is endangered , due in part to the introduction of the North American crayfish species, Pacifastacus leniusculus[20–24]. 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 , although the extent of damage is not clear . 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 .
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 . As found in other studies of other species , 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 ). 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 . 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.
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 . 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 , 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  and cold  water, whereas Marmorkrebs appears to prefer shallow, warm waters . Thus, C. japonicus probably faces greater competition from another cold water crayfish species, Pacifastacus leniusculus, than from Marmorkrebs [20–24]. Marmorkrebs may pose a threat to C. japonicus if Marmorkrebs carry crayfish plague. Cambaroides japonicus is susceptible to crayfish plague , 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. 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[32–34].
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 , 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 , similar to what was found for the native and introduced distributions of signal crayfish in Japan . 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” .
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 , 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 , and water temperature affects growth and mortality in Marmorkrebs . 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, 41–43]. Further research is needed to determine if mismatches in predicted and actual distributions of aquatic species (e.g., signal crayfish ) 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 , 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 , 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  (spreadsheet of locations used in that model available at http://marmorkrebs.org). The main difference in this dataset and that used in  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 [7–9] 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 .
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 .
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.
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.
Data on paddy field rice production for “agricultural regions” was taken from . Each agricultural region is composed of one or more prefectures. We used GIS coordinates for prefectures from , with each prefecture being assigned the average rice production for its agricultural region.
The modeling process was similar to that used in . 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 . Several studies have used terrestrial climatic variables to model the broad-scale distribution of freshwater aquatic species across terrestrial landscapes [28, 35, 41–43].
We created models using MaxEnt v. 3.3.3e [26, 60] because it discriminates suitable and unsuitable areas better than other methods [26, 61–63]. 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 . 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 [64–66]. 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  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 . 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 . 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 , 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.
Biology of crayfish. Edited by: Kawai T, Takahata M. 2010, Hokkaido University Press, Sapporo
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.
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. 10.1002/jmor.10250.
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.
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. 10.1007/s10530-008-9334-y.
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. 10.1651/08-3125.1.
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. 10.3391/ai.2010.5.4.10.
Löwe K: Gefahr aus dem Aquarium. 2010, Halle: Mitteldeutsches Druck- und Verlagshaus GmbH & Co. KG, Mitteldeutsche Zeitung, http://www.mz-web.de/servlet/ContentServer?pagename=ksta/page&atype=ksArtikel&aid=1286541137817&openMenu=1013016724285&calledPageId=1013016724285&listid=1018881578312,
Privenau K: Marmorkrebs bringt Pest. 2010, Halle: Mitteldeutsches Druck- und Verlagshaus GmbH & Co. KG, Mitteldeutsche Zeitung, http://www.mz-web.de/servlet/ContentServer?pagename=ksta/page&atype=ksArtikel&aid=1286541132341&calledPageId=987490165154,
Stloukal E: Recent distribution of non-indigenuous (sic) crayfish species in Slovakia. Folia Faunistica Slovaca. 2009, 14: 119-122.
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. 10.3391/ai.2009.4.2.19.
Faulkes Z: The spread of the parthenogenetic marbled crayfish, Marmorkrebs (Procambarus sp.), in the North American pet trade. Aquat Invasions. 2010, 5: 447-450. 10.3391/ai.2010.5.4.16.
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-
Invasion of self-cloning crayfish alarms Madagascar. http://www.earthtimes.org/articles/news/339974,alarms-madagascar-feature.html,
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. 10.3391/ai.2010.5.4.09.
Williamson MH, Fitter A: The characters of successful invaders. Biol Conserv. 1996, 78: 163-170. 10.1016/0006-3207(96)00025-0.
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. 10.1111/j.1469-8137.2007.02207.x.
Okada Y: Some observations of Japanese crayfishes. Sci Rep Tokyo Bunrika Daigaku B. 1933, 14: 155-188.
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. 10.1016/j.biocon.2006.09.002.
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. 10.1651/S-2652.1.
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. 10.1046/j.1444-2906.2002.00491.x.
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. 10.1651/C-2411.
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. 10.1023/A:1014573915464.
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. 10.1127/0003-9136/2006/0166-0001.
Rodda GH, Jarnevich CS, Reed RN: Challenges in identifying sites climatically matched to the native ranges of animal invaders. PLoS One. 2011, 6: e14670-10.1371/journal.pone.0014670.
Phillips SJ, Anderson RP, Schapire RE: Maximum entropy modeling of species geographic distributions. Ecol Model. 2006, 190: 231-259. 10.1016/j.ecolmodel.2005.03.026.
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. 10.1890/1051-0761(2007)017[0181:FNZMIR]2.0.CO;2.
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. 10.3391/ai.2011.6.1.07.
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. 10.1890/10-1171.1.
Ministry of Agriculture, Forestry and Fisheries: Appendix 3: Rice Production. Mon Stat Agric For Fish. 2011, 694: 92-
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.
Unestam T: On the host range and origin of the crayfish plague fungus. Rep Inst Freshw Res Drottningholm. 1972, 52: 192-198.
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. 10.1111/j.1523-1739.2009.01240.x.
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. 10.1007/s10530-010-9828-2.
Larson ER, Olden JD, Usio N: Decoupled conservatism of Grinnellian and Eltonian niches in an invasive arthropod. Ecosphere. 2010, 1: art16-10.1890/ES10-00053.1.
Hobbs HH: The crayfishes of Florida. Univ Fl Publ Biol Ser. 1942, 3: 1-179.
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. 10.1023/A:1007904817619.
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.
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. 10.1111/j.1365-2427.2007.01860.x.
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. 10.1899/08-151.1.
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.
Januchowski-Hartley S, VanDerWal J, Sydes D: Effective control of aquatic invasive species in tropical Australia. Environ Manag. 2011, 1-9.
Reshetnikov A, Ficetola G: Potential range of the invasive fish rotan (Perccottus glenii) in the Holarctic. Biol Invasions. 2011
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. 10.1016/j.ecoleng.2011.06.014.
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.
Blanke D, Schulz H: Situation des Edelkrebses (Astacus astacus L.) sowie weiterer Flusskrebsarten in Niedersachsen. Tagungsbericht der Deutschen Gesellschaft für Limnologie (DGL), Braunschweig, 2002. 2003, DGL, Werder, 385-389.
Holdich DM, Pöckl M: Invasive crustaceans in European inland waters. Biological Invaders in Inland Waters: Profiles, Distribution, and Threats. Edited by: Gherardi F. 2007, Springer, Dordrecht, The Netherlands, 29-75. [Drake JA (Series Editor): Invading Nature: Springer Series in Invasion Ecology]
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. 10.1672/0277-5212(2000)020[0194:DARAOT]2.0.CO;2.
Andrews EA: Sperm-transfer organs in Cambaroides. Biol Bull. 1909, 17: 257-270. 10.2307/1536028.
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.
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. 10.1651/0278-0372(2001)021[0885:OOTMBA]2.0.CO;2.
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. 10.1651/0278-0372(2002)022[0532:BOJOTJ]2.0.CO;2.
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.
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.
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.
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. 10.1046/j.1463-6395.2002.00113.x.
Yamanaka K, Kuwabara R, Shio T: Larval development of a Japanese crayfish, Cambaroides japonicus (De Haan). Bull Mar Sci. 1997, 61: 165-175.
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.
Japan GIS Datasets V2. http://www.fas.harvard.edu/~chgis/japan/datasets.html,
Phillips SJ, Dudík M, Schapire RE: A maximum entropy approach to species distribution modeling. Proceedings of the 21st International Conference on Machine Learning; Banff. Edited by: Greiner R, Schuurmans D. 2004, ACM Press, , 655-662.
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. 10.1111/j.0906-7590.2008.5378.x.
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. 10.1111/j.2006.0906-7590.04596.x.
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. 10.1111/j.0906-7590.2006.04700.x.
Köppen WP: Das geographische system der klimate. 1936, Gebrüder Borntraeger, Berlin
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.
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. 10.1111/j.1472-4642.2011.00811.x.
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. 10.1002/joc.1276.
Elith J, Kearney M, Phillips S: The art of modelling range shifting species. Methods Ecol Evol. 2010, 1: 330-642. 10.1111/j.2041-210X.2010.00036.x.
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. 10.1111/j.1472-4642.2010.00725.x.
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. 10.1017/S0376892997000088.
The authors declare no competing interests.
ZF compiled data; TPF and JM ran models; ZF, TPF, and JM wrote the paper. All authors read and approved the final manuscript.