Email updates

Keep up to date with the latest news and content from Aquatic Biosystems and BioMed Central.

Open Access Highly Accessed Research

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

Zen Faulkes1*, Teresa Patricia Feria1 and Jesús Muñoz23

Author Affiliations

1 Department of Biology, The University of Texas-Pan American, Edinburg, TX, USA

2 Real Jardín Botánico (RJB-CSIC), Plaza de Murillo 2, 28014, Madrid, Spain

3 Universidad Tecnológica Indoamérica, Bolívar 2035 y Quito, Ambato, Ecuador

For all author emails, please log on.

Aquatic Biosystems 2012, 8:13  doi:10.1186/2046-9063-8-13

The electronic version of this article is the complete one and can be found online at: http://www.aquaticbiosystems.org/content/8/1/13


Received:25 April 2012
Accepted:11 June 2012
Published:27 June 2012

© 2012 Faulkes et al.; licensee 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.

Abstract

Background

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.

Results

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.

Conclusions

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.

Background

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 [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 [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[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 [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.

Results

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

thumbnailFigure 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.

thumbnailFigure 2. Marmorkrebs distribution in Japan modeled usingP. 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].

thumbnailFigure 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.

thumbnailFigure 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.

thumbnailFigure 5. Marmorkrebs distribution in Japan modeled usingP. fallaxand 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.

thumbnailFigure 6. Marmorkrebs distribution in Japan modeled usingP. 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

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.

Discussion

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 [20-24]. 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[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 [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,41-43]. 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.

Conclusion

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.

Methods

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 webcite). 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 [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 [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/ webcite), 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 ofCambaroides japonicusdescribed in the scientific literature

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,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 [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 [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 webcite) 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.

Competing interests

The authors declare no competing interests.

Authors’ contributions

ZF compiled data; TPF and JM ran models; ZF, TPF, and JM wrote the paper. All authors read and approved the final manuscript.

References

  1. Kawai T, Takahata M (Eds): Biology of crayfish. Hokkaido University Press, Sapporo; 2010. OpenURL

  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. PubMed Abstract | Publisher Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text OpenURL

  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. OpenURL

  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. Publisher Full Text OpenURL

  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. Publisher Full Text OpenURL

  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. Publisher Full Text OpenURL

  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 webcite

    OpenURL

  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 webcite

    OpenURL

  10. Stloukal E: Recent distribution of non-indigenuous (sic) crayfish species in Slovakia.

    Folia Faunistica Slovaca 2009, 14:119-122. OpenURL

  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. Publisher Full Text OpenURL

  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. Publisher Full Text OpenURL

  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. OpenURL

  14. Invasion of self-cloning crayfish alarms Madagascar.

    http://www.earthtimes.org/articles/news/339974,alarms-madagascar-feature.html webcite

    OpenURL

  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. Publisher Full Text OpenURL

  16. Williamson MH, Fitter A: The characters of successful invaders.

    Biol Conserv 1996, 78:163-170. Publisher Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text OpenURL

  18. Okada Y: Some observations of Japanese crayfishes.

    Sci Rep Tokyo Bunrika Daigaku B 1933, 14:155-188. OpenURL

  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. Publisher Full Text OpenURL

  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. Publisher Full Text OpenURL

  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. Publisher Full Text OpenURL

  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. Publisher Full Text OpenURL

  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. Publisher Full Text OpenURL

  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. Publisher Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Phillips SJ, Anderson RP, Schapire RE: Maximum entropy modeling of species geographic distributions.

    Ecol Model 2006, 190:231-259. Publisher Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text OpenURL

  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. Publisher Full Text OpenURL

  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. PubMed Abstract | Publisher Full Text OpenURL

  30. Ministry of Agriculture, Forestry and Fisheries: Appendix 3: Rice Production.

    Mon Stat Agric For Fish 2011, 694:92. OpenURL

  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. OpenURL

  32. Unestam T: On the host range and origin of the crayfish plague fungus.

    Rep Inst Freshw Res Drottningholm 1972, 52:192-198. OpenURL

  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. PubMed Abstract | Publisher Full Text OpenURL

  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. Publisher Full Text OpenURL

  35. Larson ER, Olden JD, Usio N: Decoupled conservatism of Grinnellian and Eltonian niches in an invasive arthropod.

    Ecosphere 2010, 1:art16. Publisher Full Text OpenURL

  36. Hobbs HH Jr: The crayfishes of Florida.

    Univ Fl Publ Biol Ser 1942, 3:1-179. OpenURL

  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. Publisher Full Text OpenURL

  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. OpenURL

  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. Publisher Full Text OpenURL

  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. Publisher Full Text OpenURL

  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. OpenURL

  42. Januchowski-Hartley S, VanDerWal J, Sydes D: Effective control of aquatic invasive species in tropical Australia.

    Environ Manag 2011, 1-9. OpenURL

  43. Reshetnikov A, Ficetola G: Potential range of the invasive fish rotan (Perccottus glenii) in the Holarctic.

    Biol Invasions 2011. OpenURL

  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. Publisher Full Text OpenURL

  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. OpenURL

  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. OpenURL

  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]

    OpenURL

  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. Publisher Full Text OpenURL

  49. Andrews EA: Sperm-transfer organs in Cambaroides.

    Biol Bull 1909, 17:257-270. Publisher Full Text OpenURL

  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. OpenURL

  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. Publisher Full Text OpenURL

  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. Publisher Full Text OpenURL

  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. OpenURL

  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. OpenURL

  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. OpenURL

  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. Publisher Full Text OpenURL

  57. Yamanaka K, Kuwabara R, Shio T: Larval development of a Japanese crayfish, Cambaroides japonicus (De Haan).

    Bull Mar Sci 1997, 61:165-175. OpenURL

  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. OpenURL

  59. Japan GIS Datasets V2.

    http://www.fas.harvard.edu/~chgis/japan/datasets.html webcite

    OpenURL

  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. OpenURL

  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. Publisher Full Text OpenURL

  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. Publisher Full Text OpenURL

  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. Publisher Full Text OpenURL

  64. Köppen WP: Das geographische system der klimate. Gebrüder Borntraeger, Berlin; 1936. OpenURL

  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. OpenURL

  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. Publisher Full Text OpenURL

  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. Publisher Full Text OpenURL

  68. Elith J, Kearney M, Phillips S: The art of modelling range shifting species.

    Methods Ecol Evol 2010, 1:330-642. Publisher Full Text OpenURL

  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. Publisher Full Text OpenURL

  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. Publisher Full Text OpenURL