Comparative study of acetylcholinesterase and glutathione S-transferase activities of closely related cave and surface Asellus aquaticus (Isopoda: Crustacea)

Abstract

The freshwater isopod crustacean Asellus aquaticus has recently been developed as an emerging invertebrate cave model for studying evolutionary and developmental biology. Mostly morphological and genetic differences between cave and surface A. aquaticus populations have been described up to now, while scarce data are available on other aspects, including physiology. The purpose of this study was to advance our understanding of the physiological differences between cave A. aquaticus and its surface-dwelling counterparts. We sampled two surface populations from the surface section of the sinking Pivka River (central Slovenia, Europe), i.e. locality Pivka Polje, and locality Planina Polje, and one cave population from the subterranean section of the sinking Pivka River, i.e. locality Planina Cave. Animals were sampled in spring, summer and autumn. We measured the activities of acetylcholinesterase (AChE) and glutathione S-transferase (GST) in individuals snap-frozen in the field immediately after collection. Acetylcholinesterase is likely related to animals’ locomotor activity, while GST activity is related to the metabolic activity of an organism. Our study shows significantly lower AChE and GST activities in the cave population in comparison to both surface A. aquaticus populations. This confirms the assumption that cave A. aquaticus have lower locomotor and metabolic activity than surface A. aquaticus in their respective natural environments. In surface A. aquaticus populations, seasonal fluctuations in GST activity were observed, while these were less pronounced in individuals from the more stable cave environment. On the other hand, AChE activity was generally season-independent in all populations. To our knowledge, this is the first study of its kind conducted in A. aquaticus. Our results show that among closely related cave and surface A. aquaticus populations also physiological differences are present besides the morphological and genetic. These findings contribute to a better understanding of the biology of A. aquaticus and cave crustaceans in general.

About the Authors:

Anita Jemec

* E-mail: [email protected]

Affiliation: University of Ljubljana, Biotechnical Faculty, Department of Biology, Ljubljana, Slovenia

ORCID http://orcid.org/0000-0003-3143-4461

David Škufca

Affiliation: University of Ljubljana, Biotechnical Faculty, Department of Biology, Ljubljana, Slovenia

Simona Prevorčnik

Affiliation: University of Ljubljana, Biotechnical Faculty, Department of Biology, Ljubljana, Slovenia

Žiga Fišer

Affiliation: University of Ljubljana, Biotechnical Faculty, Department of Biology, Ljubljana, Slovenia

Primož Zidar

Affiliation: University of Ljubljana, Biotechnical Faculty, Department of Biology, Ljubljana, Slovenia

Introduction

The freshwater isopod water louse, Asellus aquaticus (L.) (Isopoda: Crustacea), is a benthic crustacean with several important roles in freshwater ecosystems. It inhabits various freshwater habitats throughout most of Europe, including caves [1], [2]. The species exhibits strong genetic structuring in the southern and eastern part of its range [3], [4], [5]. Genetically distinct cave populations that resulted from polytopic and polychronous immigration to the cave environment still have their surface counterparts [4], [6], [7]. Pairs of surface and cave populations have gained increasing recognition as a model system to address questions of evolutionary parallelism and convergence [8], [9], [10].

The transition of organisms from surface to cave habitats presumably underlies adaptations at all levels of biological organization: morphology [9], [11], [12], [13], [14], physiology [15], [16], [17], [18], behaviour [16], [19], [20], and life histories [21], [22], [23]. Among these, mostly morphological [2], [24] and genetic differences [4], [6], [7], [8] between cave and surface A. aquaticus populations have been described, while scarce data are available on other aspects, including physiology. To our knowledge, only one physiological study comparing the cave and surface A. aquaticus has been carried out, and that was on respiration rate [25]. On the other hand, quite a few studies on presumable physiological adaptations of other cave crustaceans, such as different energy saving mechanisms, are available. These showed that many cave crustaceans have an amazing ability to endure starvation over prolonged periods of time due to their enhanced capacity for food storage [17], diminished locomotor activity [17], [18], and lower metabolic activity [16], [18], [26], [27], [28].

The purpose of this study was to advance our understanding of the physiological differences between the cave and surface A. aquaticus. While previous studies assessed the physiological state of some cave and surface crustaceans after a certain period of acclimation in the laboratory [18], [26], we measured the activities of two physiologically important enzymes in individuals snap-frozen immediately after collection. Acetylcholinesterase (AChE) plays a major role in cholinergic signal transmission in the sensory and neuromuscular systems and therefore appears to be likely related to animals’ locomotor activity [29], [30], [31]. Glutathione S-transferases (GSTs) are a family of multifunctional enzymes that play a central role in the detoxification of both endogenous and xenobiotic compounds and are also involved in intracellular transport, biosynthesis of hormones and protection against oxidative stress [32], [33], [34]. Some studies suggest that GST activity could be directly related to the metabolic activity of an organism. It is commonly assumed that higher metabolism leads to higher reactive oxygen species (ROS) production and results in higher GST activities [26], [35], [36], [37]. However, ROS levels and GST activities were also demonstrated to be increased under hypoxia (anoxia) and metabolic depression [38], [39], showing that the GST activity rate is not always representative of the basal metabolic rate.

The aim of this paper was to investigate whether AChE and GST activities differ between the cave A. aquaticus and its surface-dwelling counterparts, both living in well-oxygenated water. As cave animals are expected to have a reduced level of locomotor and metabolic activities, we hypothesized that generally there would be lower AChE and GST activities in cave compared to surface individuals. Due to more stable environmental conditions in caves, cave individuals were also expected to exhibit less season-dependent enzyme activities, in contrast to surface individuals, where considerable season-dependent enzyme activities were anticipated in accordance with evidence from studies on other crustaceans [40], [41].

Materials and methods

Chemicals

The following chemicals were purchased from Sigma (Germany): dibasic and monobasic potassium phosphate; 1-chloro-2,4-dinitrobenzene; L-glutathione (reduced form); 5,5-dithiobis(2-nitrobenzoic acid); sodium hydrogen carbonate; acetylthiocholine iodide. BCA Protein Assay Reagents A and B, cadmium chloride, and potassium dichromate were purchased from Pierce (U.S.A.). All chemicals were of the highest commercially available grade, typically 99% or higher.

Study system and field work

Closely related cave and surface populations of A. aquaticus, which inhabit subterranean and surface stretches of the sinking river Pivka (central Slovenia, Europe), provided an ideal ecological setup for the needs of our study. We sampled two surface populations from the river stretches flowing across two large enclosed karstic planes: (i) just prior to the river’s sink at Pivka Polje and (ii) after its resurgence on Planina Polje, and one cave population from (iii) the river’s subterranean section in the Planina Cave (Table 1, Fig 1). To account for seasonal differences in enzymatic activity, we conducted three samplings: spring, summer and autumn (Table 1). Each time 15 adult individuals regardless of sex, yet excluding ovigerous females, were collected from each population on the same day (hereafter referred to as a sample). Animals were snap-frozen on dry ice on site and immediately transferred to the laboratory where they were stored at -20°C until enzyme analysis. A reliable sex determination of A. aquaticus requires a careful examination of animal’s appendages (gonopods and peraeopod IV) under the stereo microscope. During this procedure we would risk that the specimens would warm up, which could endanger reliable enzyme measurements. We therefore did not determine the sex of each individual. During each sampling, we also measured physical and chemical properties of water at the sample site, i.e. temperature, dissolved oxygen concentration, conductivity, and pH, using a portable multimeter CyberScan 600 (Eutech Instruments).

[Figure omitted. See PDF.]

Fig 1. Sampling localities of Asellus aquaticus.

(Digital elevation model is based on LIDAR (Light Detecting and Ranging) data freely available from ARSO—Slovenian Environment Agency at http://gis.arso.gov.si/evode/profile.aspx?id=atlas_voda_Lidar@Arso).

https://doi.org/10.1371/journal.pone.0176746.g001

[Figure omitted. See PDF.]

Table 1. Information about sampling localities, taxa and dates of collection.

https://doi.org/10.1371/journal.pone.0176746.t001

Enzyme analysis

Prior to enzyme analysis, the fresh mass of all individuals was measured (S1 Fig). Whole body homogenates were prepared in 650 μL of 50 mM phosphate buffer (pH 7.0) with 0.5% triton X100 using a T10 IKA Ultra-turrax homogenizer. The homogenate was centrifuged for 15 min at 16000 g and 4°C. Enzyme activities were measured on freshly prepared supernatants for each sample in triplicate.

AChE activity was determined according to the method of Ellman et al. [42] using microtiter plates as described by Jemec et al. [43]. Kinetic measurements were performed using acetylthiocholine iodide as a substrate (final concentration 1 mM). We mixed 1 M acetylthiocholine iodide and 2.3 mM 5,5-dithiobis(2-nitrobenzoic acid) in 1:500 (v/v) ratio. A total of 100 μL of this mixture was applied to the microtiter plate, where 50 μL of 50 mM potassium phosphate buffer (pH 7.0) and 50 μL of protein supernatant had already been added. The reaction was followed spectrophotometrically at 405 nm and 25°C for 5 min using a microplate reader (Anthos, UK).

GST activity was measured according to the method by Habig et al. [44] and optimized for microtiter plates [43] using 1,2-dichloro-4-nitrobenzene (CDNB) as a substrate. The final concentrations of substrates and reagents were: 4 mM of CDNB and 1 mM of reduced glutathione. CDNB was dissolved in ethanol and further diluted in 50 mM potassium phosphate buffer (pH 7.0) to final concentrations. The concentration of ethanol in the final reaction solution was less than 1% (v/v). We added 50 μl of the protein supernatant to start the reaction, which was followed spectrophotometrically at 340 nm and 25°C for 3 min using a microplate reader (Anthos, UK).

Protein concentration of the supernatants for enzyme analysis was measured using the BCA^™ Protein Assay Kit, a modification of the bicinchoninic acid protein assay (Pierce, Rockford, IL, USA).

The AChE activity was calculated as nmoles of hydrolysed acetylthiocholine iodide/min/mg protein (extinction coefficient Ɛ₄₁₂ = 13600 M^-1cm^-1) and GST activity as nmoles of hydrolysed CDNB/min/mg protein (extinction coefficient Ɛ₃₄₀ = 9600 M^-1cm^-1). Hereafter, both of these are referred as an enzyme unit (EU). Protein concentration was calculated using bovine serum albumin as a reference.

Data analysis

Data were statistically analysed in R 3.3.2 [45]. As variance was non-homogenous between groups, we employed a robust two-way ANOVA implemented in the R package WRS2 [46] with a pbad2way() function; a modified one-step M-estimator based on Huber’s Psi was used as a robust measure of central tendency (est parameter set to “mom”). Note that the M-estimators of the central tendency values were close to the mean values in all cases. A separate ANOVA model was run for each enzyme, i.e. AChE and GST. Enzyme activity served as a continuous response variable while population (Planina Cave, Planina Polje, Pivka Polje) and season (spring, summer, autumn) were included as categorical explanatory variables, whose interaction was tested as well. Biologically meaningful pairwise comparisons among groups were tested using the pb2gen() function and p-values were adjusted according to Benjamini & Hochberg [47]. Variability of samples was assessed using a robust measure of variability analogous to the standard deviation, i.e. median absolute deviation (MAD). Hereafter, we use the term significant difference to refer to statistically significant difference. All plots were drawn with OriginPro 8.0.

Results

AChE and GST activities

Two-way robust ANOVA showed the same general pattern in AChE and GST: enzyme activities were noticeably lower in the cave population than in either surface population (Table 2, Figs 2A and 3A). Central values (see above) of AChE activity ranged between 3–4 EU in the cave population. In surface populations these values were 2.5 to 6-times higher, ranging between 10–11 EU and 11–18 EU in the Pivka Polje and Planina Polje populations, respectively. Similarly, central values of GST activity ranged between 23–33 EU in the cave population and were up to 2.5-times higher in both surface populations; between 45–75 EU and 30–80 EU in the Pivka Polje and Planina Polje populations, respectively (Table 2). Furthermore, between-individual variability in AChE and GST activity was generally lower in the cave population compared to either surface population during all seasons, with the exception of a slightlyhly variable GST activity in the spring sample from Planina Polje (S1 Table).

[Figure omitted. See PDF.]

Fig 2. A) AChE activity in cave and surface Asellus aquaticus populations. B) Seasonal fluctuation in AChE activity in cave and surface Asellus aquaticus populations.

Asterisks indicate statistically significant differences among samples (* 0.05 < p < 0.01, ** 0.01 < p < 0.001, *** p < 0.001). M-estimators of central tendency are shown as empty squares.

https://doi.org/10.1371/journal.pone.0176746.g002

[Figure omitted. See PDF.]

Fig 3. A) GST activity in cave and surface Asellus aquaticus populations. B) Seasonal fluctuation in GST activity in cave and surface Asellus aquaticus populations.

Asterisks indicate statistically significant differences among samples (** 0.01 < p < 0.001, *** p < 0.001). M-estimators of central tendency are shown as empty squares.

https://doi.org/10.1371/journal.pone.0176746.g003

[Figure omitted. See PDF.]

Table 2. Pairwise comparisons of AChE and GST activity among cave and surface Asellus aquaticus samples.

https://doi.org/10.1371/journal.pone.0176746.t002

Pairwise comparisons (Table 2) confirmed that AChE and GST activities were significantly lower in the cave population compared to both surface populations in all seasons. Only the insignificantly different GST activities of the Planina Cave and Planina Polje populations collected in spring diverged from this pattern. The enzyme activities also differed between both surface populations, although not in all seasons and with a smaller effect size. AChE activity in the Planina Polje population was significantly higher in summer and autumn, while GST activity in the Pivka Polje population was higher in the spring.

According to the two-way robust ANOVA, seasonal fluctuation in enzyme activity in surface populations was observed, yet it did not follow any clear trend. The interaction between population and season was significant in both AChE (p = 0.012) and GST (p < 0.001), indicating that seasons had different effects on the enzyme activity of each population. For both enzymes this was mainly due to their more pronounced seasonal fluctuation in the Planina Polje population (see S2 Table for details). Pairwise comparisons revealed that seasonal fluctuation was more pronounced in the GST than in the AChE activity in all populations (Table 2, Figs 2B and 3B). The latter was almost constant in the Planina Cave and Pivka Polje samples, while it was significantly lower in spring in the Planina Polje population. GST activity in the Planina Cave population was significantly lower in autumn than in both other seasons while in the Pivka Polje population the low autumn value significantly differed only from the higher summer value. On the other hand, GST activity in the Planina Polje population was significantly lower in spring, as well as significantly higher in summer than in autumn.

Physical and chemical parameters of water at sampling localities

Water temperature and dissolved oxygen concentration are shown in Fig 4, while all values of measured parameters are presented in S3 Table. As expected, variation in water temperature in the Planina Cave was negligible during the year. In spring and autumn, temperatures were similar at all three localities, while the summer temperatures at Pivka Polje and Planina Polje were considerably higher. Dissolved oxygen concentration was slightly lower in Planina Polje compared to Pivka Polje and Planina Cave, with mutually similar concentrations. Summer dissolved oxygen concentrations in Planina Cave and Planina Polje were lower than those measured in spring and autumn, while no seasonal fluctuation was observed in Pivka Polje. Water conductivity and pH were similar at all three localities throughout the year (see S3 Table).

[Figure omitted. See PDF.]

Fig 4. Water temperature and dissolved oxygen concentration at sampling localities.

https://doi.org/10.1371/journal.pone.0176746.g004

Discussion

Our study shows significantly lower AChE and GST activities in the cave population compared to both surface populations. These results support the suggested correlation between AChE activity and animals’ locomotor activity [29], [30] and the correlation between GST activity and animals' metabolic rate [36]. Namely, lower metabolic and locomotor activity of cave versus surface crustaceans has been reported in several comparative studies on locomotor [17], [18] and metabolic activities [16], [28], [48] of other cave and surface crustaceans. Decreased metabolic rates have also been shown in predatory stygobiotic amblyopsid fish compared to non-stygobiotic ones, presumably as a result of their adaptation to an energy-poor environment [49], [50].

Contrary, Mösslacher and Creuzé Des Châtelliers [25] observed higher locomotor activity and higher respiration rate in A. aquaticus from the chemoautotrophic Movile Cave (Romania) compared to surface individuals from Austria. However, the experimental design of their study was inappropriate to estimate animals’ locomotor activity since this behaviour was measured directly (i.e. by recording the time of movement), while animals simultaneously performed another type of behaviour (i.e. feeding on decaying leaves). Even the authors of the study have therefore acknowledged the high probability for biased results and speculative conclusions. Increased activity was also observed in characid Astyanax cavefish where it is a result of enhanced food searching behaviour [51], [52] and reduced sleep duration [53]. One has to consider though, that such behaviours likely provide an adaptive advantage to the predatory stygobiotic vertebrate, while it is harder to imagine their similar benefits for a grazing invertebrate such as A. aquaticus.

In surface A. aquaticus populations, evident seasonal fluctuations in GST activity were observed, while these were less pronounced in the cave population. There are a number of potential reasons for the observed results. Firstly, multiple environmental parameters, such as temperature [36], light intensity [27], [54], dissolved oxygen concentrations [38] and availability of food [17] are known to affect the GST activity of crustaceans. These environmental parameters are less variable in the stable cave environment compared to the surface [55], which may explain the lower GST activity fluctuation in cave individuals. Secondly, while mass of surface individuals differed significantly throughout the year, differences in mass were not as pronounced in the cave population (S1 Fig). In surface populations, we collected and compared three distinct generations over the three seasons [56], [57], [58], [59]. The overwintering generation included extremely large individuals collected in spring, the spring-brood generation included small individuals collected in summer, while a mixture of larger spring-brood and smaller summer-brood individuals was collected in autumn. Mentioned generations differ from one another according to reproductive stage, which is known to influence GST activity [60] via differences in metabolic activity. For example, it has been reported that the spring-brood generation in A. aquaticus matures within around 3 months, compared to the summer-brood generation that matures within about 6 months [58]. The described differences in life-histories are also the most plausible explanation for the single non-significant difference recorded between cave and surface populations used in our study, i.e. GST activity between the Planina Cave and Planina Polje spring sample. The overwintering individuals from Planina Polje were extremely large (S1 Fig) which may result in low GST activity as previously shown for other enzymes [41].

In contrast to the GST activity, the AChE activity at each locality did not vary throughout the seasons. There are mutually contradictory literature data for AChE dependence on environmental factors. While some authors report evident seasonal AChE fluctuation in crustaceans [40], [60], [61], [62], others state that some environmental parameters, such as temperature and salinity, do not influence AChE activity when studied in laboratory experiments [40], [63]. The only exception in our study was the surface population from Planina Polje, which had significantly lower AChE activity in the spring. The reason for this phenomenon remains unknown, but as already mentioned in the case of GST, we observed that the overwintering individuals at this locality were much larger and heavier than those from the other two seasons (S1 Fig). Similar observations were made by Xuereb et al. [41], where larger individuals of the freshwater amphipod Gammarus fossarum had a lower AChE activity than smaller ones.

Interestingly, larger variation in both enzyme activities was observed between surface than between cave individuals. Our study design does not allow the distinction between genetic and environmental causes for the observed differences. Nevertheless, the lower between-individual variability in the cave A. aquaticus may be due to strong directional and stabilizing natural selection for diminished locomotor and metabolic activity that supposedly increase fitness in an energy-poor environment. Just as likely, the higher variability between surface individuals could be due to the more variable environmental conditions at the surface. The same effect of the stable vs. changing environment on the variability in GST activities was observed in terrestrial isopods [43].

The employed enzymes are commonly used as biochemical biomarkers in routine environmental quality biomonitoring programmes [64], [65]. It has been already suggested that a number of environmental parameters and life-history traits of organisms could influence the enzyme activities [66]. The current study is an additional proof that the AChE and GST activities vary considerably with regard to locality and season, probably as a result of differences in environmental parameters and specific life-histories. Therefore, we suggest that potential future biomonitoring studies employing A. aquaticus should be designed closely in line with the results presented in this study.

To our knowledge, this is the first study of its kind conducted in A. aquaticus. In this model organism, AChE has so far been measured only as a biomarker of pollution in laboratory exposures [67], [68] and we found no records of GST measurements. Our future research will focus on verifying the actual correlation between both enzymes and their physiological roles. Confirmation of a direct link between A. aquaticus AChE and GST activity with its locomotor and metabolic activities would advance the use of biochemical approaches to studying both of these latter activities. The main advantage of such measurements is the assessment of the animals’ physiological condition in-situ, i.e. at the time and site of collection, as opposed to laboratory measurements that usually include an acclimation period. Namely, acclimation in the laboratory has been shown to alter the physiology of some crustaceans [43], [69].

In conclusion, our results show that considerable physiological differences exist between closely related cave and surface A. aquaticus populations. The lower AChE and GST activities of the cave population probably reflect physiological adaptations to the specific conditions of the cave environment. Seasonal fluctuation of GST enzyme activities was considerably more pronounced in surface populations and is likely a consequence of the joint effect of fluctuating environmental conditions and life histories. These findings contribute to a better understanding of both the biology of A. aquaticus and cave crustaceans in general.

Supporting information

[Figure omitted. See PDF.]

S1 Fig. Seasonal fresh mass fluctuation in cave and surface Asellus aquaticus populations.

https://doi.org/10.1371/journal.pone.0176746.s001

(DOCX)

S1 Table. Specific AChE and GST activities of cave and surface Asellus aquaticus samples.

https://doi.org/10.1371/journal.pone.0176746.s002

(DOCX)

S2 Table. Specific AChE and GST activities.

Robust two-way ANOVA post hoc tests comparisons of samples.

https://doi.org/10.1371/journal.pone.0176746.s003

(DOCX)

S3 Table. Physical and chemical parameters of water at sampling localities.

https://doi.org/10.1371/journal.pone.0176746.s004

(DOCX)

Acknowledgments

We are grateful to Cene Fišer for his contribution to the early versions of the manuscript. We acknowledge reviewers for their valuable suggestions.

Author Contributions

1. Conceptualization: AJ PZ.

2. Data curation: AJ PZ SP.

3. Formal analysis: ZF DS.

4. Investigation: AJ DS.

5. Methodology: AJ.

6. Resources: AJ.

7. Software: DS ZF.

8. Supervision: AJ PZ.

9. Validation: SP.

10. Visualization: ZF AJ.

11. Writing – original draft: AJ.

12. Writing – review & editing: SP PZ ZF.

Citation: Jemec A, Škufca D, Prevorčnik S, Fišer Ž, Zidar P (2017) Comparative study of acetylcholinesterase and glutathione S-transferase activities of closely related cave and surface Asellus aquaticus (Isopoda: Crustacea). PLoS ONE 12(5): e0176746. https://doi.org/10.1371/journal.pone.0176746

References

1. Sket B. Distribution of Asellus aquaticus (Crustacea: Isopoda: Asellidae) and its hypogean populations at different geographic scales, with a note on Proasellus istrianus. Hydrobiologia. 1994;287(1):39–47.

2. Prevorčnik S, Blejec A, Sket B. Racial differentiation in Asellus aquaticus (L.)(Crustacea: Isopoda: Asellidae). Arch. Hydrobiol. 2004;160(2):193–214.

3. Verovnik R, Prevorčnik S, Jugovic J. Description of a neotype for Asellus aquaticus Linné, 1758 (Crustacea: Isopoda: Asellidae), with description of a new subterranean Asellus species from Europe. Zoologischer Anzeiger-A Journal of Comparative Zoology. 2009;248(2):101–18.

4. Konec M, Prevorčnik S, Sarbu SM, Verovnik R, Trontelj P. Parallels between two geographically and ecologically disparate cave invasions by the same species, Asellus aquaticus (Isopoda, Crustacea). J Evol Biol. 2015;28(4):864–75. pmid:25728816

5. Sworobowicz L, Grabowski M, Mamos T, Burzyński A, Kilikowska A, Sell J, et al. Revisiting the phylogeography of Asellus aquaticus in Europe: insights into cryptic diversity and spatiotemporal diversification. Fresh.Biol. 2015;60(9):1824–40.

6. Verovnik R, Sket B, Prevorčnik S, Trontelj P. Random amplified polymorphic DNA diversity among surface and subterranean populations of Asellus aquaticus (Crustacea: Isopoda). Genetica. 2003;119(2):155–65. pmid:14620955

7. Verovnik R, Sket B, Trontelj P. Phylogeography of subterranean and surface populations of water lice Asellus aquaticus (Crustacea: Isopoda). Mol. Ecol. 2004;13(6):1519–32. pmid:15140095

8. Protas ME, Trontelj P, Patel NH. Genetic basis of eye and pigment loss in the cave crustacean, Asellus aquaticus. Proc Natl Acad Sci U S A 2011;108(14):5702–7. pmid:21422298

9. Protas M, Jeffery WR. Evolution and development in cave animals: from fish to crustaceans. Wiley Interdisciplinary Reviews: Dev. Biol. 2012;1(6):823–45.

10. Stahl BA, Gross JB, Speiser DI, Oakley TH, Patel NH, Gould DB et al. A Transcriptomic Analysis of Cave, Surface, and Hybrid Isopod Crustaceans of the Species Asellus aquaticus. PloS one. 2015;10(10):e0140484. pmid:26462237

11. Culver DC, Kane TC, Fong DW. Adaptation and natural selection in caves: the evolution of Gammarus minus. Harvard University Press; 1995.

12. Yamamoto Y, Byerly MS, Jackman WR, Jeffery WR. Pleiotropic functions of embryonic sonic hedgehog expression link jaw and taste bud amplification with eye loss during cavefish evolution. Dev. Biol. 2009;330(1):200–11. pmid:19285488

13. Bilandžija H, Ma L, Parkhurst A, Jeffery WR. A potential benefit of albinism in Astyanax cavefish: downregulation of the oca2 gene increases tyrosine and catecholamine levels as an alternative to melanin synthesis. PLoS One. 2013;8(11):e80823. pmid:24282555

14. Yoshizawa M, Yamamoto Y, O'Quin KE, Jeffery WR. Evolution of an adaptive behavior and its sensory receptors promotes eye regression in blind cavefish. BMC biology. 2012;10(1):1.

15. Hervant F. The activities of enzymes associated with the intermediary and energy metabolism in hypogean and epigean crustaceans. Comptes rendus de l'Académie des sciences. Série III, Sciences de la vie. 1996a; 319(12):1071–7.

16. Hervant F, Mathieu J, Garin D, Freminet A. Behavioral, ventilatory, and metabolic responses of the hypogean amphipod Niphargus virei and the epigean isopod Asellus aquaticus to severe hypoxia and subsequent recovery. Physiol. Zool. 1996b:1277–300.

17. Hervant F, Mathieu J, Barré H, Simon K, Pinon C. Comparative study on the behavioral, ventilatory, and respiratory responses of hypogean and epigean crustaceans to long-term starvation and subsequent feeding. Comp Biochem Physiol A Mol Integr Physiol. 1997;118(4):1277–83.

18. Hervant F, Mathieu J, Culver DC. Comparative responses to severe hypoxia and subsequent recovery in closely related amphipod populations (Gammarus minus) from cave and surface habitats. Hydrobiologia. 1999;392(2):197–204.

19. Parzefall J. Behavioural adaptations. In: Encyclopedia of caves. Editor: Culver D.C., White W.B. Elsevier Academic Press, 2012.

20. Fišer Ž*, Novak L*, Luštrik R, Fišer C. Light triggers habitat choice of eyeless subterranean but not of eyed surface amphipods. The Science of Nature. 2016;103(1–2):1–2. *equal contribution

21. Fišer C, Zagmajster M, Zakšek V. Coevolution of life history traits and morphology in female subterranean amphipods. Oikos. 2013;122(5):770–8.

22. Riesch R, Martin RA, Langerhans RB. Predation’s role in life-history evolution of a livebearing fish and a test of the Trexler-DeAngelis model of maternal provisioning. Am. Nat.;181(1):78–93. pmid:23234847

23. Taylor SJ, Krejca JK, Niemiller ML, Dreslik MJ, Phillips CA. Life history and demographic differences between cave and surface populations of the western slimy salamander, Plethodon albagula (Caudata: Plethodontidae), in Central Texas. Herpetol. Conserv. Biol. 2015;10 (2):740–52.

24. Turk S, Sket B, Sarbu Ş. Comparison between some epigean and hypogean populations of Asellus aquaticus (Crustacea: Isopoda: Asellidae). Hydrobiologia. 1996;337(1–3):161–70.

25. Mösslacher F, Creuzé des Châtelliers M. Physiological and behavioural adaptations of an epigean and a hypogean dwelling population of Asellus aquaticus (L.)(Crustacea, Isopoda). Arch. Hydrobiol. 1996;138(2):187–98.

26. Simčič T, Lukančič S, Brancelj A. Comparative study of electron transport system activity and oxygen consumption of amphipods from caves and surface habitats. Fresh. Biol. 2005; 50(3):494–501.

27. Simčič T, Brancelj A. The effect of light on oxygen consumption in two amphipod crustaceans–the hypogean Niphargus stygius and the epigean Gammarus fossarum. Mar. Freshw. Behav. Physiol. 2007;40(2):141–50.

28. Wilhelm FM, Taylor SJ, Adams GL. Comparison of routine metabolic rates of the stygobite, Gammarus acherondytes (Amphipoda: Gammaridae) and the stygophile, Gammarus troglophilus. Fresh. Biol. 2006;51(6):1162–74.

29. Xuereb B, Lefèvre E, Garric J, Geffard O. Acetylcholinesterase activity in Gammarus fossarum (Crustacea Amphipoda): linking AChE inhibition and behavioural alteration. Aquat.Toxicol. 2009a;94(2):114–22.

30. Jensen CS, Garsdal L, Baatrup E. Acetylcholinesterase inhibition and altered locomotor behavior in the carabid beetle Pterostichus cupreus. A linkage between biomarkers at two levels of biological complexity.). Environ.Toxicol. Chem. 1997;16(8):1727–32.

31. Wilczek G, Babczynska A, Migula P, Wencelis B. Activity of esterases as biomarkers of metal exposure in spiders from the metal pollution gradient. Polish Journal of Environmental Studies, 2003; 12(6): 765–772.

32. Hayes JD, Flanagan JU, Jowsey IR. Glutathione transferases. Annu. Rev. Pharmacol. Toxicol. 2005; 45: 51–88. pmid:15822171

33. Halliwell B, Gutteridge JMC. Free radicals in biology and medicine. NewYork, USA: Oxford, 2007.

34. Raza H. Dual localization of glutathione S‐transferase in the cytosol and mitochondria: implications in oxidative stress, toxicity and disease. FEBS Journal. 2011;278(22):4243–51. pmid:21929724

35. Janssens BJ, Childress JJ, Baguet F, Rees JF. Reduced enzymatic antioxidative defense in deep-sea fish. Journal of Experimental Biology, 2000; 203:3717–3725. pmid:11076735

36. Lushchak VI. Environmentally induced oxidative stress in aquatic animals. Aquat.Toxicol. 2011;101(1):13–30. pmid:21074869

37. Gismondi E, Beisel JN, Cossu-Leguille C. Influence of gender and season on reduced glutathione concentration and energy reserves of Gammarus roeseli. Environ. Res. 2012;118:47–52. pmid:22769238

38. Lawniczak M, Romestaing C, Roussel D, Maazouzi C, Renault D, Hervant F. Preventive antioxidant responses to extreme oxygen level fluctuation in a subterranean crustacean. Comp Biochem Physiol A Mol Integr Physiol. 2013;165(2):299–303. pmid:23545443

39. Welker AF, Moreira DC, Campos É G, Hermes-Lima M. Role of redox metabolism for adaptation of aquatic animals to drastic changes in oxygen availability. Comparative Biochemistry and Physiology Part A: Molecular & Integrative Physiology, 2013; 165(4): 384–404.

40. Menezes S, Soares AM, Guilhermino L, Peck MR. Biomarker responses of the estuarine brown shrimp Crangon crangon L. to non-toxic stressors: temperature, salinity and handling stress effects. J Exp Mar Bio Ecol. 2006;335(1):114–22.

41. Xuereb B, Chaumot A, Mons R, Garric J, Geffard O. Acetylcholinesterase activity in Gammarus fossarum (Crustacea Amphipoda): intrinsic variability, reference levels, and a reliable tool for field surveys. Aquat. Toxicol. 2009b;93(4):225–33.

42. Ellman GL, Courtney KD, Andres V, Featherstone RM. A new and rapid colorimetric determination of acetylcholinesterase activity. Biochem. Pharmacol. 1961;7(2):88–95.

43. Jemec A, Drobne D, Remškar M, Sepčić K, Tišler T. Effects of ingested nano‐sized titanium dioxide on terrestrial isopods (Porcellio scaber). Environ.Toxicol. Chem. 2008;27(9):1904–14. pmid:19086208

44. Habig WH, Pabst MJ, Jakoby WB. Glutathione S-transferases the first enzymatic step in mercapturic acid formation. J. Biol. Chem. 1974; 249(22):7130–9. pmid:4436300

45. R Core Team R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/. (21.6. 2016), 2015

46. Mair P, Schoenbrodt., Wilcox R. WRS2: Wilcox robust estimation and testing. 2016 https://cran.r-project.org/web/packages/WRS2/index.html

47. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol. 1995; 1:289–300.

48. Issartel J, Renault D, Voituron Y, Bouchereau A, Vernon P, Hervant F. Metabolic responses to cold in subterranean crustaceans. J. Exp. Biol. 2005;208(15):2923–9.

49. Poulson TL. Cave adaptation in amblyopsid fishes. Amer. Midl. Nat. 1963; 70: 257–290.

50. Poulson TL. Animals in aquatic environments: animals in caves. 1964; 749–771. In: Dill D.B. (ed.) Handbook of Physiology: Environment, Amer. Physiol. Soc., Washington, D.C.

51. Hüppop K. Food-finding ability in cave fish (Astyanax fasciatus). Int J Speleol 1987;16: 59–66.

52. Yoshizawa M, Goricki Š, Soares D, Jeffery WR. Evolution of a behavioral shift mediated by superficial neuromasts helps cavefish find food in darkness. Curr Biol 2010; 20: 1631–1636. pmid:20705469

53. Duboué ER, Keene AC, Borowsky RL. Evolutionary convergence on sleep loss in cavefish populations. Curr Biol 2011; 21: 671–676. pmid:21474315

54. Fanjul-Moles ML, Bosques-Tistler T, Prieto-Sagredo J, Castañón-Cervantes O, Fernández-Rivera-Río L. Effect of variation in photoperiod and light intensity on oxygen consumption, lactate concentration and behavior in crayfish Procambarus clarkii and Procambarus digueti. Comp Biochem Physiol A Mol Integr Physiol. 1998;119(1):263–9. pmid:11253792

55. Culver D. C. and Pipan T. The biology of caves and other subterranean habitats. Oxford University Press, USA. 2009

56. Steel EA. Some observations on the life history of Asellus aquaticus (L.) and Asellus meridianus Racovitza (Crustacea: Isopoda). In Proceedings of the Zoological Society of London 1961 (Vol. 137, No. 1, pp. 71–87). Blackwell Publishing Ltd.

57. Migliore L, Paluzzi R, Vitagliano-Tadini G. Reproductive activity in Asellus aquaticus L.(Crustacea, Isopoda) from southern Italy. Int. J. Inver. Rep. Dev. 1982;4(6):359–67.

58. Štrus J., and Blejec A.. Reproductive activity in Asellus aquaticus (Crustacea, Isopoda) from Ljubljansko Barje. Bioloski Vestnik 1983; 31: 83–92.

59. Iversen TM, Thorup J. A three years' study of life cycle, population dynamics and production of Asellus aquaticus L. in a macrophyte rich stream. Int. Rev. Gesamten Hydrobiol. 1988;73(1):73–94.

60. Leiniö S, Lehtonen KK. Seasonal variability in biomarkers in the bivalves Mytilus edulis and Macoma balthica from the northern Baltic Sea. Comp Biochem Physiol C Pharmacol Toxicol Endocrinol. 2005;140(3):408–21.

61. Dellali M, Barelli MG, Romeo M, Aissa P. The use of acetylcholinesterase activity in Ruditapes decussatus and Mytilus galloprovincialis in the biomonitoring of Bizerta lagoon. Comp Biochem Physiol C Pharmacol Toxicol Endocrinol. 2001;130(2):227–35.

62. Robillard S, Beauchamp G, Laulier M. The role of abiotic factors and pesticide levels on enzymatic activity in the freshwater mussel Anodonta cygnea at three different exposure sites. Comp Biochem Physiol C Pharmacol Toxicol Endocrinol. 2003;135(1):49–59.

63. Scaps P, Borot O. Acetylcholinesterase activity of the polychaete Nereis diversicolor: effects of temperature and salinity. Comp Biochem Physiol C Pharmacol Toxicol Endocrinol. 2000;125(3):377–83.

64. Cajaraville MP, Bebianno MJ, Blasco J, Porte C, Sarasquete C, Viarengo A. The use of biomarkers to assess the impact of pollution in coastal environments of the Iberian Peninsula: a practical approach. Sci. Total. Environ. 2000;247(2):295–311.

65. Damásio JB, Barata C, Munné A, Ginebreda A, Guasch H, Sabater S et al.Comparing the response of biochemical indicators (biomarkers) and biological indices to diagnose the ecological impact of an oil spillage in a Mediterranean river (NE Catalunya, Spain). Chemosphere. 2007;66(7):1206–16. pmid:16959291

66. Jemec A, Drobne D, Tišler T, Sepčić K. Biochemical biomarkers in environmental studies—lessons learnt from enzymes catalase, glutathione S-transferase and cholinesterase in two crustacean species. Environ. Sci.Poll. Res. 2010;17(3):571–81.

67. O'Neill AJ, Galloway TS, Browne MA, Dissanayake A, Depledge MH. Evaluation of toxicity in tributaries of the Mersey estuary using the isopod Asellus aquaticus (L.). Mar. Environ. Res. 2004;58(2):327–31.

68. Bouskill NJ, Handy RD, Ford TE, Galloway TS. Differentiating copper and arsenic toxicity using biochemical biomarkers in Asellus aquaticus and Dreissena polymorpha. Ecotox. Environ. Saf. 2006;65(3):342–9.

69. Lešer V, Drobne D, Vilhar B, Kladnik A, Žnidaršič N, Štrus J. Epithelial thickness and lipid droplets in the hepatopancreas of Porcellio scaber (Crustacea: Isopoda) in different physiological conditions. Zoology. 2008;111(6):419–32. pmid:18656333