Research Article |
Corresponding author: Sarah R. Goodnight ( sarah.r.goodnight@gmail.com ) Academic editor: Ting Hui Ng
© 2025 Nicholas T. Lewis, Sarah R. Goodnight, Daya Hall-Stratton, Amy E. Fowler.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Lewis NT, Goodnight SR, Hall-Stratton D, Fowler AE (2025) Size-dependent desiccation tolerance in adult and juvenile introduced freshwater Japanese mystery snails (Cipangopaludina japonica, previously Heterogen japonica). Aquatic Invasions 20(2): 231-249. https://doi.org/10.3391/ai.2025.20.2.150239
|
Freshwater Japanese mystery snails (Cipangopaludina japonica, previously Heterogen japonica) were introduced to North America from Asia in the early 1900s and have colonized many lake and river systems across the United States (US). Tolerance to environmental stressors, such as desiccation, plays a large role in species’ invasion potential and persistence in novel environments. To characterize the desiccation tolerance of C. japonica snails, adults and juveniles from three eastern US populations were exposed to air for 13.5 weeks (adults, n = 650) or 48 hours (juveniles, n = 849) and their mortality assessed over time. Over 50% of adult snails from each population exposed to desiccation survived over 10 weeks of constant air exposure, while survival ranged from 10 to 64% at the end of the exposure experiment (13.5 weeks), depending on population, indicating exceedingly high resistance to desiccation mortality in adults of this species. In contrast, juvenile snails were much more vulnerable to desiccation, with over 70% mortality at just 24 hours of drying and only a single individual surviving 48 hours of desiccation stress. We found that the interaction between snail shell length and time affected survival for both adults and juveniles, where larger body sizes were associated with increased probability of survival as time of exposure increased (p < 0.001 for both juveniles and adults). Based on these data, juveniles cannot survive long-term air exposure, but the high desiccation tolerance of adults may facilitate survival and population persistence in stressful environments and allow for increased dispersal between water bodies. Therefore, both commercial and recreational users of water bodies containing introduced C. japonica should be aware of the risk of unintentional dispersal between water bodies via contaminated gear and/or boats, even if those materials are exposed to air for a significant amount of time.
Desiccation stress, management, mortality, invasive, shell length, Viviparidae
Invasive species are a global concern, threatening the biodiversity and function of many ecosystems (
Many successful aquatic invasive species can tolerate a broad range of common stressors, such as salinity, desiccation, and temperature stress (
Attempts to control invasive snails include a variety of strategies such as physical removal, molluscicides, or even biological control, with mixed results (
However, many aquatic snails have the ability to enter estivation, or dormancy, to survive prolonged drying events (
Two invasive freshwater snail species, the Japanese mystery snail, Cipangopaludina japonica (previously Heterogen / Bellamya / Viviparus japonica) and the Chinese mystery snail, Cipangopaludina chinensis (previously Bellamya chinensis), are native to parts of Asia and were introduced to North America over 100 years ago as a food source (
In spite of the extent of the mystery snail invasion in North American water bodies, there are few comprehensive studies identifying the specific aspects of their biology that afford them such success, and the mechanisms of dispersal between affected areas need further study (but see
In this study we evaluated the desiccation tolerance of both adult and juvenile Japanese mystery snails in a laboratory environment by placing snails in a temperature-controlled incubator and measuring survival over time. Based on previous work on Chinese mystery snails, we expected snails to survive at least eight weeks out of water. We also tested if snail size predicted survival across the experiment, with the expectation that larger snails would have higher tolerance to desiccation due to the ability to retain more water within their shells. This study represents the first desiccation challenge for invasive populations of C. japonica.
Adult Japanese mystery snails were collected in August and September of 2023 from three previously surveyed sites in central and northern Virginia, USA: Leesylvania State Park, on the Potomac River (38°35.1033'N, 77°15.7867'W), Lake Royal (38°48.2417'N, 77°17.455'W), a residential lake, and a stormwater pond in Maryland (39°10.6'N, 76°43.0467'W), hereafter designated “Maryland Pond”. These three populations were chosen based on previous genetic identification using the barcoding gene cytochrome c oxidase subunit 1 (COI), which confirmed they were Japanese mystery snails (
Live, fully formed juveniles were removed from the female’s brood-pouch, and only juveniles ≥ 3.8 mm shell length were used in experiments (n = 849, mean shell length = 5.64 mm, range = 3.8–8.2 mm; Table
Total number of adult and juvenile snails used in each experiment, with mean shell lengths in mm and standard deviations. Snails, both adults and juveniles, were significantly different in size between all three populations.
Population | Total snails used in experiment (n) | Shell length range (mm) | Mean shell length (mm) ± SD | |
---|---|---|---|---|
LEESYLVANIA | Adults | 234 | 17.8–63.0 | 52.16 ± 10.72 |
Juveniles | 400 | 3.8–7.8 | 5.46 ± 0.79 | |
MARYLAND POND | Adults | 207 | 15.2–55.5 | 40.06 ± 8.76 |
Juveniles | 200 | 4.3–8.1 | 5.98 ± 0.91 | |
LAKE ROYAL | Adults | 209 | 34.9–62.3 | 42.70 ± 4.84 |
Juveniles | 249 | 3.8–8.2 | 5.65 ± 0.90 |
For the experiment, juveniles were separated by population and all placed in the same incubator (daytime temperature = 32.2 °C; nighttime temperature = 21.1 °C; 12:12 light/dark cycle) on dry plates with moisture wicks immediately after removal from the female. Humidity was not explicitly measured during the experiment. After placement in the incubator, two (Maryland Pond and Lake Royal) or four (Leesylvania) juveniles per brood were randomly subsampled at 4-, 8-, 10-, 24-, and 48-hour time periods (for a total of 10 or 20 juveniles/brood/site sampled across five time periods = 50 or 100 juveniles per replicate experiment). At each time point after removal from the incubator, juveniles were individually placed in 60 mL of room temperature dechlorinated water in 18-well plastic bead boxes for a 24-hour rehydration period. After rehydration, juveniles were assessed for mortality by pushing the operculum with sharp forceps and looking under a dissecting microscope for movement within the shell. Juveniles were recorded as dead if they floated on the top of the water and/or if they had no response to their operculum being pushed in. The shell length of each juvenile was taken upon subsampling from the incubator. Snails were not fed before or during the experiment or during the rehydration period. One snail from Maryland Pond was accidentally crushed before being evaluated for survival or shell length, thus there are only 49 snails for that replicate trial rather than 50.
A total of 650 adults were used for the experiment (Leesylvania State Park, n = 234; Lake Royal, n = 207; Maryland Pond, n = 209) and were exposed to air concurrently using the same conditions as the juvenile experiment. Leesylvania snails were placed in the incubator on August 28th, 2023, while snails from the Maryland Pond and Lake Royal were placed in the incubator one week later. All snails were placed in the same incubator for the entire experiment.
Snails were kept in dry bins with moisture wicks and high edges, to prevent movement and escape, for 13.5 weeks (95 days). Snails were separated into bins by site of origin and placed on the two middle shelves of the incubator (Suppl. material
Data from adults and juveniles were analyzed separately. Statistical analyses and data wrangling were accomplished in the R programming environment (v4.4.0,
We then conducted two separate model-based analyses. The first approach focused on determining the impacts of source population and experimental timepoint on overall snail survival probability, and the second focused on the impact of shell length on individual snail survival over time. We conducted both analyses for the adult and juvenile experiments separately. For the first approach, we constructed generalized linear models (GLMs) or generalized linear mixed models (GLMMs) assuming binomial error distributions and logit link functions, where we modeled the number of live snails at each timepoint as the response and the number of dead snails as the binomial denominator. For both adult and juvenile experiments, we included timepoint (adults, week sampled; juveniles, hour sampled) and population (Leesylvania, Maryland Pond, or Lake Royal) as fixed effects. For the juvenile experiment, due to the replicate structure of the experiment (i.e., we conducted four or five independent trials), we also included replicate trial as a random effect.
For our second approach, to test for the impact of individual snail shell length on survival across time for both adults and juveniles, we constructed GLMMs treating each snail as an independent replicate with shell length (log scale) and timepoint (hours for juveniles, weeks for adults) as fixed effects. We also included population as a random effect in these models because we were primarily interested in determining the impact of shell size on individual survival across sites, which was not possible in our first approach. Additionally, there were site-level differences in shell size distributions (Table
We employed a model selection approach for both model-based analyses and tested all possible combinations of fixed effects. Models were selected based on the lowest AICc score. Models were fit using glmmTMB (
The three adult snail populations had significantly different mean snail shell lengths (in mm); snails sampled from the Leesylvania State Park were significantly larger than snails from the other two sites (p < 0.001, Table
As expected, the proportion of live juveniles decreased with increasing time of desiccation exposure (p < 0.001, Fig.
Survival probability of juvenile snails from three populations over a 48-hour period of desiccation stress. Plotted lines represent predicted values from a binomial generalized linear mixed model (GLMM), and shaded areas represent 95% confidence intervals plotted over raw data. Colors correspond to the population from which snails were collected.
When we investigated individual juvenile survival probability, the most parsimonious model included the interaction between time exposed (in hours) and shell length (in mm). Concordant with our previous analysis, individual probability of survival decreased with increasing hour of exposure (estimate = -0.679, SE = 0.158, p < 0.001; Table
Interaction between juvenile snail body size (shell length in mm) and timepoint sampled (in hours) predicting individual survival probability over 24 hours of desiccation stress. Fitted lines are Bernoulli generalized linear mixed model (GLMM) predictions with 95% confidence envelopes plotted over raw data points. Only one individual snail survived at the final timepoint of 48 hours, so data from that timepoint is excluded from this figure.
Juvenile individual survival table across all populations. Individual survival estimates and associated errors and confidence intervals are presented at all time points of the experiment. Predictions generated from a Bernoulli generalized linear mixed model (GLMM) with experimental replicate and population as random effects.
Time (hours) | Predicted survival probability | Standard error | Lower 95% CI | Upper 95% CI |
---|---|---|---|---|
4 | 84.7% | 77.6% | 54.8% | 96.2% |
8 | 70.9% | 76.9% | 35.1% | 91.7% |
10 | 61.8% | 76.7% | 26.5% | 87.9% |
24 | 8.4% | 80.4% | 1.9% | 30.7% |
48 | 0.1% | 101.9% | <0.1% | 0.5% |
Similar to results from our juvenile experiment, snail population and time exposed to desiccation (in weeks) interacted to determine adult snail overall survival (Fig.
Survival probability of adult mystery snails from three populations over 13.5 weeks of desiccation stress. Trendlines represent predictions from a binomial generalized linear model (GLM) and shaded areas represent 95% confidence intervals, plotted over jittered raw data. Colors correspond to population.
As for individual-level survival, predicted survival probabilities across populations ranged from 90–98% through five weeks of exposure, with many adults surviving to the end of the experimental period at 13.5 weeks (Table
Interaction between body size (shell length in mm) and week sampled predicting survival probability of adult snails. Fitted lines are Bernoulli generalized linear mixed model (GLMM) predictions with 95% confidence envelopes plotted over raw data points. Whole weeks of the experiment (1 through 13) are represented.
Adult individual survival table across all populations. Individual survival estimates and associated errors and confidence intervals are presented at every whole week of the experiment as well as at experiment end (13.5 weeks). Predictions generated from a Bernoulli generalized linear mixed model (GLMM) with population as a random effect.
Time (weeks) | Predicted survival probability | Standard error | Lower 95% CI | Upper 95% CI |
---|---|---|---|---|
1 | 98.5% | 30.2% | 94.3% | 99.7% |
2 | 97.8% | 29.4% | 92.4% | 99.5% |
3 | 96.7% | 28.6% | 90.0% | 99.2% |
4 | 95.3% | 28.0% | 87.1% | 98.7% |
5 | 93.4% | 27.4% | 83.5% | 98.0% |
6 | 91.1% | 27.0% | 79.4% | 97.0% |
7 | 88.2% | 26.7% | 74.6% | 95.6% |
8 | 84.6% | 26.5% | 69.2% | 93.8% |
9 | 80.4% | 26.5% | 63.2% | 91.5% |
10 | 75.5% | 26.5% | 56.8% | 88.7% |
11 | 70.1% | 26.7% | 50.1% | 85.3% |
12 | 64.1% | 27.0% | 43.3% | 81.4% |
13 | 57.8% | 27.5% | 36.7% | 76.9% |
13.5 | 54.6% | 27.7% | 33.4% | 74.5% |
Here we present the first study to our knowledge investigating the desiccation tolerance of invasive Japanese mystery snails (Cipangopaludina japonica). Mystery snails are likely to continue to spread to new waterbodies via anthropogenic (i.e., attachment to boats, equipment, or macrophytes) or “natural” (i.e., floating on the surface of the water, dispersal via bird or mammal predators, or flooding events) vectors (
Desiccation-resistant snails like C. japonica (for example, Pomacea spp. apple snails and Chinese mystery snails) are found in their native range in ecosystems like rice paddies, which are characterized by fluctuations in water level and frequent drying (
We also observed a positive relationship between snail shell length and survival probability at later stages of the adult experiment (Fig.
Our results are consistent with other studies that show high tolerance to desiccation for large freshwater snail species. For example, adults of the closely-related Chinese mystery snail (C. chinensis) can survive nine weeks of drying, and larger snails have higher survival rates than smaller snails, similar to our study (
In contrast to our findings in the adults, juvenile snails had high mortality when desiccated, representing a much lower probability of survival in dry environments as compared to fully developed adult snails. Specifically, predicted juvenile mortality at 24 hours of desiccation was nearly 90% when estimated across populations, while mortality at 48 hours was essentially 100% (Fig.
In accordance with our adult experiment, we again observed an increasingly positive relationship between juvenile snail size and survival as time exposed increased (Fig.
The three populations from which we collected snails had significantly different size distributions and tolerance to desiccation, suggesting that geographically isolated invasive populations have varying abilities to survive environmental fluctuations and colonize new areas (e.g.,
Further, we also hypothesize that effects of environmental stressors on adult snails likely cascaded to offspring, explaining the concurrent patterns we saw of low survival of both adults and juveniles from Lake Royal (Figs
An important factor in addition to survival probability when desiccated is the potential overhang effects of stress. That is, extended periods of desiccation and concurrent starvation likely have lasting impacts on snail fecundity, growth, and even future resistance to parasite infection (
Lastly, the degree to which this desiccation advantage for C. japonica adults compares to other species of potential native gastropod competitors commonly found in the same habitats (such as Elimia virginica, Physella/Physa spp., Campeloma decisum, and Planorbella trivolvis; (
This research highlights an important adaptation that facilitates C. japonica colonization of dynamic water bodies, dispersal to novel environments, and potential resistance to contemporary management methods. We found that large adult snails sourced from low-stress environments may represent a particularly high risk of spreading to other water bodies and may survive transport through dry, stressful conditions. Our study suggests that due to their high desiccation tolerance and reproductive strategy of bearing live young, large gravid females specifically may be able to establish introduced populations in dynamic environments previously considered unfit for stable snail populations. This work suggests that management of C. japonica populations must accommodate their enhanced ability to survive out of water.
NTL conducted the experiment and collected the data, interpreted the data, wrote the first draft of the manuscript, and reviewed all subsequent drafts. SRG assisted in curating data, conducted statistical analyses, and assisted in co-writing and reviewing/editing all manuscript drafts. DHS assisted in conceptualizing the project, collecting specimens and data, and reviewed and edited all manuscript drafts. AEF conceptualized the project, supervised project design and resources, assisted in data curation and analysis, and reviewed and edited all drafts of the manuscript.
George Mason University (GMU)’s Biology Research Semester funded this work and had no role in the planning and design of the study or publication of this manuscript.
Data and statistical code have been deposited in a public Figshare repository (accessible at https://figshare.com/s/ffdff14a9ed8b6a8e75b).
The authors would like to thank George Mason University and the Potomac Science Center for support in this project. We acknowledge Geraldine Grant, PhD, Arndt Laemmerzahl, PhD, Pratyush Jaishanker, and Durwood Moore for assistance and guidance in developing the project and working with the data. We also express our gratitude to the anonymous peer reviewers whose comments improved this manuscript.
Statistical tables and additional figure
Data type: docx
Explanation note: fig. S1. Diagram of snail spatial distribution in the incubator for the adult experiment. table S1. Top models and corresponding weights for overall juvenile survival analysis. table S2. Juvenile overall survival table separated by population. table S3. Top models and corresponding weights for individual juvenile survival analysis. table S4. Top models and corresponding weights for overall adult survival analysis. table S5. Adult overall survival table separated by population. table S6. Top models and corresponding weights for individual adult survival analysis.