- Research article
- Open Access
Occurrence of the invasive Spanish slug in gardens: can a citizen science approach help deciphering underlying factors?
BMC Ecologyvolume 18, Article number: 23 (2018)
The Spanish slug (Arion vulgaris, also known as A. lusitanicus) is considered one of the most invasive species in agriculture, horticulture and private gardens all over Europe. Although this slug has been problematic for decades, there is still not much known about its occurrence across private gardens and the underlying meteorological and ecological factors. One reason for this knowledge gap is the limited access of researchers to private gardens. Here we used a citizen science approach to overcome this obstacle and examined whether the occurrence of Arionidae in Austrian gardens was associated with meteorological (air temperature, precipitation, global solar radiation, relative humidity) or ecological factors (plant diversity, earthworm activity). Occurrence of the invasive A. vulgaris versus the similar-looking native A. rufus was compared using a DNA-barcoding approach.
Slugs were collected from 1061 gardens from the dry Pannonian lowland to the wet alpine climate (altitudinal range 742 m). Slug abundance in gardens was best explained and negatively associated with the parameters “sum of the mean air temperature in spring”, “number of frost days in the previous winter” and “mean daily global solar radiation on the day of data collection”. Precipitation, plant diversity and earthworm activity were also related to slug abundance, but positively. Out of our genetic sampling of collected slugs, 92% belonged to A. vulgaris.
Our study showed that citizen science (i) is a feasible approach to record species occurrence in restricted areas across a wide geographical range and (ii) could be more widely employed in order to identify underlying environmental factors of species occurrence.
Invasive alien species have been recognized as a threat to biodiversity, human health and economy by many countries . Through global trade routes organisms are transported both deliberately and accidentally around the globe, sometimes finding a new home in areas they never could have reached on their own. In many cases these species remain unnoticed until they establish viable populations, spread quickly and outcompete native species, or become pestiferous in agriculture, horticulture or private gardens. At this stage residents of areas where the dispersal occurs are often among the first to notice this change in their environment [2, 3]. As a consequence, many environmental protection agencies, research facilities and municipalities use citizen science as a method of choice to detect invasive species at early stages [4, 5]. Citizen science describes a method where citizens and scientists work together in scientific projects. The citizens usually collect data, analyze data or are even involved in the whole scientific process .
In the current study we tested the citizen science approach to assess the occurrence of the Spanish slug Arion vulgaris Moquin-Tandon 1855 (formerly known as A. lusitanicus auct. non-Mabille) in private gardens in Austria. Urban biota can host a variety of native and non-native species  and gardens are particularly interesting areas, because they are ‘designed’ by humans  and therefore species numbers often vary from one garden to the next. In Austria A. vulgaris was described for the first time in 1971 in a village near Vienna . The origin of this slug is still controversial [10, 11], but it is nevertheless listed among the 100 most invasive species in Europe . In the meanwhile it has become one of the most common pest slugs in private gardens, horticulture and agriculture in central Europe [13,14,15]. The Spanish slug feeds on many vegetables and ornamentals as well as dead plant material [16, 17] and can survive low temperatures and droughts [18, 19]. Under normal circumstances, most Spanish slugs die after oviposition at the end of the year , but can live up to 3 years if they are unable to breed . They can adapt to a variety of environmental factors and therefore compete with native slug species . The activity of A. vulgaris is strongly correlated to humidity [19, 23, 24] and temperature [18, 21]. Additionally, the time of day [25, 26] and the availability of suitable shelters affect slug activity [25, 27]. Some studies also observed that earthworms [28, 29] or plant diversity  affected the extent of herbivory by slugs.
In many areas, A. vulgaris co-occurs with two similar looking native Arion species (the red slug A. rufus and the black slug A. ater). Recent studies have shown that no clear mating barrier exists between these three species and that fertile hybrids occur [31,32,33,34]. This makes traditional species identification difficult and genetic analyses a useful addition [35, 36].
Hence, the goals of the current study was (I) to test a citizen science approach in order to determine the occurrence and abundance of slugs in particular across gardens in Austria, (II) to determine any association of slug abundance with meteorological factors, plant diversity and earthworm activity in these gardens, and (III) to investigate the occurrence of the invasive A. vulgaris versus the native A. rufus by using DNA sequencing.
Study area and data collection
This study focused on private gardens across Austria. These gardens were distributed across all Austrian federal states, but most of them were located in the highly urbanized regions of Vienna and the greater Vienna area. Although the study covers a broad range of geographical conditions, there is more weight to the landscapes of the Vienna region: hilly terrains and flatlands with meadows, forests and cultivated areas. We chose a citizen science approach to independently assess the abundance and occurrence of slugs in Austrian gardens in the spring seasons 2014 and 2015. In 2014 the assessment took place from the 2nd of April to the 12th of June, in 2015 from the 3rd of April to the 15th of May. Participants were undergraduate students at the University of Natural Resources and Life Sciences, Vienna. All participants got 1 h training explaining the goal of this project and the sampling protocol. This method allowed us to get data on slug and plant diversity from private properties. To enable people to collect scientifically sound data, we designed a simple protocol. The participants used a cardboard quadrat (20 × 20 cm) that they randomly placed in their garden. This quadrat served as an artificial shelter for slugs during the day. No baits where used to attract slugs. The quadrats were left undisturbed on the same spot for 2 days, after which the slugs under the quadrats where counted, photographed and collected. Additionally, data were collected on altitude of the collection site, on the number of plant species within a radius of 2 m around the quadrat, number of earthworm casts under the quadrat as a measure for earthworm activity  and previous and current control methods against molluscs employed in the garden.
In spring 2014 participants used the smartphone app Epicollect and Epicollect+, respectively . Epicollect and Epicollect+ are open source smartphone applications available for Android and iOS for data collection in the field. The collected data were stored in a central database and can be exported in various formats for further analyses. In our study participants used the apps to submit pictures of slugs and data on location, altitude, number of plant species, number of earthworm casts and slug management methods. In 2015 the participants entered the data for the same parameters in a spreadsheet (Microsoft Excel) and uploaded it together with the associated photographs at the university’s online learning platform based on Moodle (https://learn.boku.ac.at/). The number of slugs was deduced from the spreadsheet and the photographs uploaded via app or on the Moodle platform, respectively, in the year 2015. In 2014 slugs were collected by the participants and preserved in 80% ethanol for further inspection by us. Slug numbers were then obtained by counting of the preserved slugs. All data were validated using the photographs the participants sent or the slug samples, respectively. Inconsistent or incomplete data sets were discarded and not included in further analyses.
In addition, we procured data on weather conditions from the nearest weather stations from each location of the year before the collection and the collection year from Austrian Federal meteorological and geophysical data provider (Zentralanstalt für Meteorologie und Geodynamik ZAMG). From these data we extracted various measures for vapour pressure (hPa), air temperature (°C), relative air humidity (%), precipitation (mm) and global radiation (J/cm2) from the collection day and periods before collection. See Table 1 for more details on the parameters considered for statistical analysis.
DNA barcoding approach
To identify and discriminate the invasive slug A. vulgaris from other slug species we used a DNA barcoding approach on 127 slug garden samples representing different altitudes and regions within Austria. The gardens were selected to cover all altitudinal levels and provinces represented in the whole data set.
DNA was extracted using the protocol of DNeasy® Blood & Tissue Kit, following the manufacturer’s protocol (Qiagen, Hilden, Germany). Incubation time was usually 1 h at 56 °C. Genomic DNA extracts were stored at − 20 °C until processing. The DNA barcode region of the cytochrome C oxidase subunit 1 (COI) gene (cf. ) was amplified using the LCO-1490 forward and the HCO-2198 reverse primers . Amplification was performed by taking 1 µl of genomic DNA extract and adding 16.3 µl double distilled H2O, 2.5 µl PCR Buffer, 2.5 µl MgCl2 (25 mM), 0.5 µl dNTP (2 mM), 0.5 µl of either forward or reverse primer (10 pmol/µl), 1 µl BSA and 0.2 µl Taq (5 U/µl). We used the following PCR protocol: initial step of 4 min at 94 °C, 5 cycles of 94 °C, 60 s/44 °C, 90 s/72 °C, 90 s, followed by 30 cycles of 94 °C, 60 s/49 °C, 75 s/72 °C, 75 s, and a post-treatment of 7 min at 72 °C. In order to check whether the DNA extraction and PCR amplification was successful, we run agarose gel electrophoresis. When the results of the gel electrophoresis were not convincing, we repeated DNA extraction and the amplification step and ran agarose gel electrophoresis again. Final cycle sequencing was processed by LGC Genomics (Berlin, Germany).
After a pre-check running the NCBI BLAST search algorithm (https://blast.ncbi.nlm.nih.gov/Blast.cgi?PROGRAM=blastn&PAGE_TYPE=BlastSearch&LINK_LOC=blasthome), we used eleven published reference DNA sequences of five different Arion species [41,42,43,44] and of the grey garden slug Deroceras reticulatum  for slug identification. At first, all sequences (n = 138) were manually aligned using BioEdit 7.2.5 . Subsequently, genetic distances among individual slug samples were quantified by using the K2P model , as is commonly done in non-phylogenetic, but species-identification DNA barcoding analysis [cf. 48]. The genetic distance matrix was used to generate a neighbor-joining (NJ) tree  with bootstrap support (BS > 60%) values indicated above species branches based on 10,000 replicates. These analyses were all performed running PAUP* (version 4.0b10; ). Different DNA sequences of the mitochondrial COI marker obtained in our analysis were deposited in the ‘European Nucleotide Archive’ under the accession numbers LS974172-LS974197.
All statistical analyses were done using R (Version 3.0.1) and R Studio (Version 1.0.136) . In order to identify which factors are strongly correlated we first calculated a heatmap using Spearman rank correlations with the R package “stats”. The resulting heatmap displays the strength of the relationships of the correlation matrix in colour codes and additionally groups the factors that are intercorrelated. Therefore, it is possible to visually identify factors that are rather independent from each other and have a combined effect on the response variable . This makes it possible to select most important factors based on the results of former studies on slug ecology and on the current results of our heatmap.
In a next step we used the R package MASS to calculate a generalized linear model and a stepwise forward and backward Poisson regression based on the identified factors from the heatmap and compared the Akaike Information Criteria (AIC) of both models to identify the one with the better fit. To account for potential multicollinearity we also calculated the variance inflation factors (VIF).
After consolidation of the received samples from citizen scientists, we obtained 1061 garden data sets, comprised of 687 in 2014 and 374 in 2015, sometimes including more than one slug, in 2014 and 2015, respectively (Fig. 1). In 459 of the 1061 data sets (43%) slug control measures were applied either in the year of data collection or in the year before.
The altitudinal range of the observations was from 134 to 876 m a.s.l. Altogether 1269 slugs were analysed, with a mean of 1.2 ± 2.5 SD (standard deviation) slugs and a maximum of 20 slugs under one sampling quadrat. On average, 3.9 ± 3.3 SD plant species were found within a radius of 2 m around the quadrat; the maximum was 31 different plant species within this area in one garden. Earthworm activity measured as the number of earthworm casts directly under the quadrat was 0.7 ± 1.6 SD earthworm casts with a maximum of 19 casts under one quadrat in one garden.
Heatmaps representing correlation of slug abundance and possible influencing factors are presented in Fig. 2. Groups of intercorrelated parameters are visible along the main diagonal. Overall, six different parameter groups may be distinguished. Out of these six groups we chose the sum of the mean temperature from the year previous to data collection (smetpy), the number of days below 0 °C from the winter before data collection (September to April; frostwint), the applied slug control methods in the previous year and the year of data collection (comc), the global radiation on the day of data collection (rad), the mean vapour pressure on the day of data collection (mvp), the number of days without precipitation from March to September from the year previous to data collection (drypy), the sum of the mean daily temperature from the beginning of the year until the day of data collection (smetcy), the precipitation on the day of data collection (prec) and the meters above sea level of the sampling site (alt) and additionally the number of plants in a radius of 2 m around the data collection point (plants) and the number of earthworm casts under the cardboard (worms) for further analyses (Table 1). We chose these parameters because on the one hand they represent their respective groups and on the other hand they are often found to affect slugs in the literature. Individual correlations with slug abundance were generally low.
To test the combined effect of these parameters on slug number we employed a generalized linear model (Table 2). A negative correlation was found between the number of slugs and the sum of the mean daily temperature from the beginning of the year until the day of data collection (smetcy), the number of days below 0 °C from the winter before data collection (September to April; frostwint), the global radiation on the day of data collection (rad), the mean vapour pressure on the day of data collection (mvp), the number of days without precipitation from March to September from the year previous to data collection (drypy), the sum of the mean temperature from the year previous to data collection (smetpy) and the meters above sea level of the sampling site (alt). A positive correlation was found with the number of earthworm casts under the cardboard (worms), the applied slug control methods in the previous year and the year of data collection (comc), the precipitation on the day of data collection (prec) and the number of plants in a radius of 2 m around the data collection point (plants). All variance inflation factors stayed sufficiently low (using a threshold value of 5 as a rule of thumb), so multicollinearity was not an issue. We checked the model for redundant information by computing a stepwise parameter fit. Based on AIC the stepwise regression model did not perform better (AIC for both models was 4308.1), therefore we chose the generalized linear model fitted to the full set of eleven parameters selected from the heatmap.
The length of the alignment of 138 individual slug COI sequences was 638 base pairs. The tree topology is displayed as a NJ phylogram based on K2P genetic distances between these samples (Fig. 3).
Our 127 newly generated individual slug sequences are representative of five different species: Arion vulgaris, A. fuscus, A. distinctus, A. fasciatus and Deroceras reticulatum. However, most of our samples (92%) represent A. vulgaris; and further species, like for instance the native A. rufus, were not detected within our data set. All species-branches had 100% bootstrap support indicating unambiguous species identification of slugs using this DNA barcoding approach. Grouping of individual sequences reflecting altitudinal and/or geographical clusters was not observed (i.e. genotypes observed are randomly distributed across Austria).
Our results indicate that meteorological parameters and plant diversity have a significant impact on the abundance of slugs in Austrian gardens. This confirms general results found in previous studies also on the level of individual gardens and indicates that citizen science is a suitable method for assessing the abundance of slugs on private properties.
In general, the use of a citizen science approach made it possible to collect data on a large geographical scale, covering all Austrian provinces in only 3 months. Additionally, we could get access to data from private properties, which we could not get otherwise. The project potentially raised the participants’ awareness for invasive species and ecosystem processes, although we did not measure this. The trade-off we had to accept was that we could not collect as extensive data as in a conventional project, in order to not overstrain participants and to keep motivation high.
We found a significant positive correlation between number of slugs observed and precipitation on collection day (prec), suggesting a high slug activity during moister conditions [24, 25, 54]. The significant negative association between slug abundance and mean vapour pressure (mvp) is connected to the vapour pressure deficit, the difference between the actual vapour pressure and the vapour pressure the air would have, when it is saturated with moisture, which has been described as one main factor affecting slug activity [19, 24].
The finding that the sum of the mean daily temperature from the sampling year (smetcy), the sum of the daily mean temperature of the year previous to data collection (smetpy) and the number of frost days from the previous winter (frostwint) influence slug abundance in gardens also supports the results of previous studies [18, 21, 23, 55], although in some respects other than expected:
The negative correlation between the number of slugs and the temperature sum of the sampling year (smetcy) and the year previous to data collection (smetpy) contradicts findings of colleagues, who demonstrated higher slug activity with rising temperatures [24, 56]. There could be multiple possible explanations for this finding.
Temperatures under the artificial slug shelters could have risen above tolerable limits, which made the slugs leave their shelter  seeking a cooler, shaded refuge under vegetation. Hommay et al.  recommended therefore inspections of slug shelters in the early morning before temperatures rise. Although we used cardboard quadrats that do not heat up very fast this could be a (partial) explanation for these results, since we did not check for the time when participants investigated the shelters. We therefore recommend including collection day time in future studies using a citizen science approach.
Another possible explanation could be that temperature is connected to global radiation (rad), which is also negatively correlated with slug numbers in our study and was reported to also affect slug activity [31, 58]. Moreover, in the time of late spring and early summer the temperatures can rise above 30 °C especially in Eastern Austria. As most of the data were collected in this region during rather warm temperatures the lower slug numbers could be the result of suboptimal temperatures for slug development [20, 21, 59]. High temperatures often are connected to drought, a parameter (drypy) correlated significantly negative to slug abundance in our study.
Another parameter in our results that negatively affected slug abundance was winter temperature (frostwint) or more specifically the number of frost days. Organisms have generally two strategies to overwinter: freeze avoidance or freeze tolerance. Some slug species, like A. vulgaris, have been shown to have a tolerance for frost [18, 60], others seek frost-free shelters. Particularly in gardens, slugs find good shelter in compost heaps or near heated buildings. However, the number of frost days in winter nevertheless seems to be an important parameter even in gardens, which could also be connected to the parameter altitude (alt), since at higher altitudes winters are usually harsher and we found a negative correlation between altitude and slug abundance.
Moreover, our analysis revealed that general slug abundance was positively correlated with plant diversity (plants) suggesting support of the “more individuals”-hypothesis  that states that basically more plant species result in higher productivity and therefore more herbivores. It is also likely that a higher diversity in gardens is associated with a higher richness in structures and potential shelters for slugs, which have been identified as key factor for slug activity, reproduction and survival [25, 62]. Particularly, big slugs are reliant on shelter in the vegetation, because they hardly find shelter in the soil . Additionally, higher plant species richness provides more diverse food sources and more shade for slugs, which could also have an effect on slug numbers.
Something that looked rather contradicting on first sight was the rise in the number of slugs when slug control methods were applied (as showed before, in 43% of the datasets slug control measures have been applied). This does not necessarily mean that control methods did not work properly, but maybe that people who experience more slug damage in their gardens were also more willing to apply control methods. However, this interpretation would need more research in the future.
Our current finding that slug numbers are positively correlated with earthworm activity (worms) suggests that earthworms’ effects on soil structure and moisture  benefits slugs in gardens. This seems to contradict other studies reporting that slug herbivory was lower when earthworms were active [28, 29]. However, those studies showed that earthworms stimulated the production of anti-herbivore compounds in the plants which could still take place in the gardens we investigated. Other explanations could be that slugs and earthworms prefer similar habitats, or that this correlation is a methodological artifact, because slugs could have been drawn due to the cardboard shelter to areas where earthworms already resided.
All these above mentioned factors were associated with slug abundance in gardens, in general. Using the citizen science approach we did not distinguish between slug species. However, our DNA analyses of representative samples showed that almost all slugs were identified as A. vulgaris (syn. A. lusitanicus), and no native A. rufus was among them. These results suggest that A. vulgaris is most likely the predominant slug species in (Austrian) gardens [65, 66]. Since in recent years more and more hybrids between the invasive A. vulgaris and the native A. rufus from Northern Europe and Germany have been reported [31,32,33] and individuals are difficult to distinguish, the use of DNA barcoding approaches seem to be helpful as additional identification tool, the more so as citizens prefer removal of slugs from the gardens and even juvenile individuals could be determined, exactly.
Although gardens are ‘designed’ by humans , our results indicate that meteorological factors like temperature or precipitation are important when studying slug abundance. Because our finding about the role of human affected parameters such as plant diversity on slug abundance contradicts other findings , further studies on the role of plant diversity on slug abundance would help to shed some light on these interactions. Additionally, interactions between earthworms and slug herbivory need to be investigated further in the field. Although previous studies showed negative effects of earthworms on slug herbivory in greenhouse experiments, our results indicate a positive relation between earthworm activity and slug abundance, which could also be a result of slugs being drawn to our artificial shelters. Accurate species identification may be difficult in a citizen science project and therefore genetic testing can help to overcome this issue. In summary, our study showed that a citizen science approach can help in gathering data on species occurrences from areas with restricted access (private properties) and can also be used to investigate relationships between species occurrences and meteorological as well as ecological parameters. As citizen science is integrated in early detection programs for invasive alien species by many institutions [4, 5] our results encourage other citizen science projects on early detection and environmental associations.
Sundseth K. Invasive alien species. Eur Comm Environ. 2014. http://www.physorg.com/tags/invasive+alien+species/. https://doi.org/10.2779/69473. Accessed 30 Sept 2017.
Dickinson JL, Shirk J, Bonter D, Bonney R, Crain RL, Martin J, et al. The current state of citizen science as a tool for ecological research and public engagement. Front Ecol Environ. 2012;10:291–7.
Pocock MJ, Chapman DS, Sheppard LJ, Roy HE. Developing a strategic framework to support citizen science implementation in SEPA. Final Report on behalf of SEPA. NERC Centre for Ecology & Hydrology. 2014. https://www.sepa.org.uk/media/163159/a-strategic-framework-to-support-the-implementation-of-citizen-science-for-environmental-monitoring.pdf. Accessed 4 June 2017.
Gallo T, Waitt D. Creating a successful citizen science model to detect and report invasive species. Bioscience. 2011;61:459–65.
Roy HE, Rorke SL, Beckmann B, Booy O, Botham MS, Brown PMJ, et al. The contribution of volunteer recorders to our understanding of biological invasions. Biol J Linn Soc. 2015;115:678–89.
Sanz FS, Holocher-Ertl T, Kieslinger B, García FS, Silva CG. White paper on citizen science for Europe. 2014. http://socientize.eu/?q=eu/content/white-paper-citizen-science. Accessed 12 Aug 2017.
Lososová Z, Horsák M, Chytrý M, Čejka T, Danihelka J, Fajmon K, et al. Diversity of Central European urban biota: effects of human-made habitat types on plants and land snails. J Biogeogr. 2011;38:1152–63.
Bergey EA, Figueroa LL. Residential yards as designer ecosystems: effects of management on land snail species composition. Ecol Appl. 2016;26:2538–47.
Reischütz PL, Stojaspal FJ. Bemerkenswerte Mollusken aus Ostösterreich. Mitteilungen der Zool Gesellschaft Braunau. 1972;1:339–43.
Pfenninger M, Weigand A, Bálint M, Klussmann-Kolb A. Misperceived invasion: the Lusitanian slug (Arion lusitanicus auct. non-Mabille or Arion vulgaris Moquin-Tandon 1855) is native to Central Europe. Evol Appl. 2014;7:702–13.
Zemanova MA, Knop E, Heckel G. Phylogeographic past and invasive presence of Arion pest slugs in Europe. Mol Ecol. 2016;25:5747–64.
Rabitsch W. DAISIE species factsheet Arion vulgaris. 2006. http://www.europe-aliens.org/speciesFactsheet.do?speciesId=52937. Accessed 6 June 2017.
Kozlowski J. The distribution, biology, population dynamics, and harmfullness of Arion lusitanicus in Poland. J Plant Prot Res. 2007;47:219–30.
Frank T, Friedli J. Laboratory food choice trials to explore the potential of common weeds to reduce slug feeding an oilseed rape. Biol Agric Hortic. 1999;17:19–29.
Frank T. Slug damage and numbers of the slug pests Arion lusitanicus and Deroceras reticulatum in oilseed rape grown beside sown wildflower strips. Agric Ecosyst Environ. 1997;67:67–78.
Briner T, Frank T. The palatability of 78 wildflower strip plants to the slug Arion lusitanicus. Ann Appl Biol. 1998;133:123–33.
Kozlowski J. Host plants and harmfulness of the Arion lusitanicus Mabille, 1868 slug. J Plant Prot Res. 2005;45:221–33.
Slotsbo S, Hansen LM, Holmstrup M. Low temperature survival in different life stages of the Iberian slug, Arion lusitanicus. Cryobiology. 2011;62:68–73.
Slotsbo S, Fisker KV, Hansen LM, Holmstrup M. Drought tolerance in eggs and juveniles of the Iberian slug, Arion lusitanicus. J Comp Physiol B. 2011;181:1001–9.
Kozłowski J. Reproduction of Arion lusitanicus Mabille, 1868 (Gastropoda: Pulmonata: Arionidae) introduced in Poland. Folia Malacol. 2000;8:87–94.
Slotsbo S, Damgaard C, Hansen LM, Holmstrup M. The influence of temperature on life history traits in the Iberian slug, Arion lusitanicus. Ann Appl Biol. 2013;162:80–8.
Knop E, Reusser N. Jack-of-all-trades: phenotypic plasticity facilitates the invasion of an alien slug species. Proc Biol Sci. 2012;279:4668–76.
Young AG, Port GR, Green DB. Development of a forecast of slug activity—validation of models to predict slug activity from meteorological conditions. Crop Prot. 1993;12:232–6.
Crawford-Sidebotham TJ. The influence of weather upon the activitiy of slugs. Oecologia. 1972;9:141–54.
Hommay G, Lorvelec O, Jacky F. Daily activity rhythm and use of shelter in the slugs Deroceras reticulatum and Arion distinctus under laboratory conditions. Ann Appl Biol. 1998;132:167–85.
Grimm B, Schaumberger K. Daily activity of the pest slug Arion lusitanicus under laboratory conditions. Ann Appl Biol. 2002;141:35–44.
Hommay G, Kienlen JC, Jacky F, Gertz C. Daily variation in the number of slugs under refuge traps. Ann Appl Biol. 2003;142:333–9.
Trouvé R, Drapela T, Frank T, Hadacek F, Zaller JG. Herbivory of an invasive slug in a model grassland community can be affected by earthworms and mycorrhizal fungi. Biol Fertil Soils. 2013;50:13–23.
Zaller JG, Parth M, Szunyogh I, Semmelrock I, Sochurek S, Pinheiro M, et al. Herbivory of an invasive slug is affected by earthworms and the composition of plant communities. BMC Ecol. 2013;13:1–10.
Fabian Y, Sandau N, Bruggisser OT, Kehrli P, Aebi A, Rohr RP, et al. Diversity protects plant communities against generalist molluscan herbivores. Ecol Evol. 2012;2:2460–73.
Roth SSS, Hatteland BA, Solhoy T, Solhøy T. Some notes on reproductive biology and mating behaviour of Arion vulgaris Moquin-Tandon 1855 in Norway including a mating experiment with a hybrid of Arion rufus (Linnaeus 1758) x ater (Linnaeus 1758). J Conchol. 2012;41:249–57.
Hatteland BA, Solhøy T, Schander C, Von Proschwitz T, Noble LR. Introgression and differentiation of the invasive slug Arion vulgaris from native A. ater. Malacologia. 2015;58:303–21.
Allgaier C. How can two soft bodied animals be precisely connected? A miniature quick-connect system in the slugs, Arion lusitanicus and Arion rufus. J Morphol. 2015;276:631–48.
Dreijers E, Reise H, Hutchinson JMC. Mating of the slugs Arion lusitanicus auct. non Mabille and A. rufus (L.): different genitalia and mating behaviours are incomplete barriers to interspecific sperm exchange. J Molluscan Stud. 2013;79:51–63.
Kałuski T, Soroka M, Kałuski T, Soroka M. Genetic studies on the invasive slug Arion lusitanicus Mabille, 1868 (Gastropoda: Pulmonata) in Poland. Folia Malacol. 2011;19:259–65.
Rowson B, Anderson R, Turner JA, Symondson WOC, Rowson B, Anderson R, et al. The slugs of Britain and Ireland: undetected and undescribed species increase a well-studied, economically important fauna by more than 20%. PLoS ONE. 2014;9:e91907.
Zaller JG, Arnone JA. Activity of surface-casting earthworms in a calcareous grassland under elevated atmospheric CO2. Oecologia. 1997;111:249–54.
Aanensen DM, Huntley DM, Feil EJ, Spratt BG. EpiCollect: linking smartphones to web applications for epidemiology, ecology and community data collection. PloS ONE. 2009;4(9):e6968.
Hawlitschek O, Morinière J, Lehmann GUC, Lehmann AW, Kropf M, Dunz A, et al. DNA barcoding of crickets, katydids and grasshoppers (Orthoptera) from Central Europe with focus on Austria, Germany and Switzerland. Mol Ecol Resour. 2015;17(5):1037–53.
Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994;3:294–9.
Hall T. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999;41:95–8.
Haro RJ, Gillis R, Cooper ST. First report of a terrestrial slug (Arion Fasciatus) living in an aquatic habitat. Malacologia. 2004;45:451–2.
Soroka M, Skujienė G. Species identification of slugs of genus Arion Férussac, 1819 (Mollusca, Pulmonata) on the basis of genetics studies. Ekologija. 2011;57:70–80.
Telfer A, Young M, Quinn J, Perez K, Sobel C, Sones J, et al. Biodiversity inventories in high gear: DNA barcoding facilitates a rapid biotic survey of a temperate nature reserve. Biodivers Data J. 2015;3:e6313.
Ahn S-J, Martin R, Rao S, Choi M-Y. The complete mitochondrial genome of the gray garden slug Deroceras reticulatum (Gastropoda: Pulmonata: Stylommatophora). Mitochondrial DNA Part B. 2017;2:255–6.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980;16:111–20.
Hebert PDN, Penton EH, Burns JM, Janzen DH, Hallwachs W. Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator. Proc Natl Acad Sci. 2004;101:14812–7.
Saitou N, Nei M. The neighbour-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4:406–25.
Swofford DL. Phylogenetic Analysis Using Parsimony. Options. 2002;42:294–307.
R Core Team. R: a language and environment for statistical computing. Vienna, Austria; 2013. https://www.R-project.org/Accessed 4 Jun 2017.
Van Loon AF, Laaha G. Hydrological drought severity explained by climate and catchment characteristics. J Hydrol. 2015;526:3–14.
Soroka M, Kozłowski J, Wiktor A, Kałuski T. Distribution and genetic diversity of the terrestrial slugs Arion lusitanicus Mabille, 1868 and Arion rufus (Linnaeus, 1758) in Poland based on mitochondrial DNA. Folia Biol (Praha). 2008;57:71–81.
Dahirel M, Ansart A, Madec L. Stage- and weather-dependent dispersal in the brown garden snail Cornu aspersum. Popul Ecol. 2013;56:227–37.
El-Danasoury H, Iglesias-Piñeiro J, Córdoba M. The effect of climate manipulations on the herbivory of the pest slug Deroceras reticulatum (Müller, 1774) (Pulmonata: Agriolimacidae). Int J Biometeorol. 2016;60:1501–7.
Grimm B, Kaiser H. Daily activity of the pest slug Arion lusitanicus Mabille. J Molluscan Stud. 2000;66:125–30.
Cordoba M, Iglesias J, Ribadulla P, Castillejo J. Performance of permanent refuge traps for the assessment of slug populations in pastureland. Ann Appl Biol. 2011;159:130–40.
Rollo CD. The regulation of activity in populations of the terrestrial slug Limax maximus (Gastropods: Limacidae). Res Popul Ecol. 1982;24:1–32.
Mellanby K. Slugs at low temperatures. Nature. 1961;189:944.
Slotsbo S, Hansen LM, Jordaens K, Backeljau T, Malmendal A, Nielsen NC, et al. Cold tolerance and freeze-induced glucose accumulation in three terrestrial slugs. Comp Biochem Physiol A Mol Integr Physiol. 2012;161:443–9.
Srivastava DS, Lawton JH. Why more productive sites have more species: an experimental test of theory using tree-hole communities. Am Nat. 1998;152:510–29.
South A. Terrestrial slugs—biology, ecology and control. London: Chapman & Hall; 1992. p. 196–202.
Eggenschwiler L, Speiser B, Bosshard A, Jacot K. Improved field margins highly increase slug activity in Switzerland. Agron Sustain Dev. 2012;33:349–54.
Blouin M, Hodson ME, Delgado EA, Baker G, Brussaard L, Butt KR, et al. A review of earthworm impact on soil function and ecosystem services. Eur J Soil Sci. 2013;64:161–82.
Kappes H, Kobialka H. Die Nacktschneckengesellschaften in NW-Deutschland (Gastropoda: Milacidae, Boettgerillidae, Limacidae, Agriolimacidae, Arionidae): ein Ergebnis der NRW-Kartierung. Natur und Heimat. 2009;69:73–94.
Noble LR, Jones CS. A molecular and ecological investigation of the large arionid slugs of North-West Europe: the potential for new pests. In: Symondson WOC, Liddell JE, editors. Ecology of agricultural pests: biochemical approaches. Springer, 1996. p. 93–131.
DD and JGZ conceived and designed the study; DD, MK and JGZ performed the study; GL analysed the data; DD, MK, GL and JGZ wrote the paper. All authors read and approved the final manuscript.
We want to thank all students from the course “Biology of terrestrial animals” from the years 2014 and 2015 participating in this project.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Consent for publication
This study was funded by the Federal Ministry of Agriculture, Forestry, Environment and Water Management (project number 100994).
Research involving animals
We followed all ethical standards based on Austrian law and controlled by the ethics commission of the University of Natural Resources and Life Sciences, Vienna.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.