ACTA BIOLOGICA SLOVENICA LJUBLJANA 2012 Vol. 55, Št. 1: 15-27 ЛВ Preliminary multispecies test of a model for non-lethal estimation of metabolic activity in freshwater crayfish Preliminarni test modela za oceno metabolne aktivnosti pri več vrstah potočnih rakov Tatjana Simčič*, Franja Pajk, Anton Brancelj, Al Vrezec Department of Freshwater and Terrestrial Ecosystems Research, National Institute of Biology, Večna pot 111, SI-1000 Ljubljana, Slovenia ^correspondence: tatjana.simcic@nib.si Abstract: We tested the applicability of electron transport system (ETS) derived from a single leg as a tool for non-lethal assessment of metabolic activity in freshwater crayfish. ETS activity of the whole body and of a leg was measured in four crayfish (Arthropoda, Crustacea, Decapoda) species: two European (Astacus astacus, Austro-potamobius torrentium), and two North American (Orconectes limosus, Pacifastacus leniusculus). Mass scaling of whole body ETS activity (ETSwhde) and leg ETS activity (ETSieg) was not significantly different for the European A. astacus and the American O. limosus. Therefore common models were constructed and tested on the remaining two species. The ratio ETSwhoie/ETSleg was significantly positively related to body mass. In the first model (model 1) ETSwhde was calculated from ETSleg multiplied by the ratio estimated from the known body mass. ETSwhde of A. torrentium was underestimated by this model, because they mature at smaller body size than the larger species. A direct relation between ETSleg and ETSwhde was therefore proposed as a general model (model 2), since they are correlated similarly in the studied species. The results show that model 2 is suitable for estimating the whole body ETS activity from leg ETS activity for the four investigated decapods. Keywords: electron transport system (ETS) activity, crayfish, size scaling, method Izvleček: V raziskavi smo testirali uporabnost modela za oceno metabolne aktivnosti pri potočnih rakih z merjenjem aktivnosti elektronskega transportnega sistema (ETS) le na eni nogi. Aktivnost ETS smo merili na celih osebkih in na nogi pri štirih vrstah potočnih rakov: dveh evropskih (Astacus astacus, Austropotamobius torrentium), in dveh severnoameriških (Orconectes limosus, Pacifastacus leniusculus). Ker se spreminjanje aktivnosti ETS celega telesa (ETSwhde ) in aktivnosti ETS noge (ETSleg) ni razlikovalo med jelševcem A. astacus in trnavcem O. limosus, smo naredili skupni model, ki smo ga testirali na preostalih vrstah. Razmerje ETSwhde/ETSleg je bilo v pozitivnem razmerju z maso telesa. Pri prvem modelu (model 1) smo ETSwhde izračunali iz ETSleg tako, da smo jo pomnožili z razmerjem, ki smo ga ocenili iz znane mase rakov. Izkazalo se je, da smo s tem modelom podcenili ETSwhde pri koščaku A. torrentium, verjetno zaradi nastopa zrelosti pri manjši velikosti kot pri večjih vrstah. Tako smo predlagali kot splošni model neposredno povezanost med ETSleg in ETSwhde, ki sta pri vseh vrstah rakov v podobnem razmerju. Rezultati so pokazali, da na podlagi modela 2 lahko na osnovi izmerjene aktivnosti ETS pri nogi napovemo aktivnost ETS pri vseh testiranih vrstah rakov. Ključne besede: aktivnost elektronskega transportnega sistema (ETS), potočni raki, velikost, metoda Introduction Freshwater crayfish are the largest freshwater macroinvertebrates. They are becoming recognized increasingly for their importance in the natural elimination of the dead organisms (Covich et al. 1999, Nyström 2002). In Europe, the increasing loss of freshwater habitats, coupled with the spread of non-indigenous North American crayfish species, and by the infection of Aphanomyces astaci that selectively kills the European species, populations of European crayfish has dramatically reduced (Holdich et al. 2009). The ecology of the native and introduced species and interactions between them have been studied intensively (e.g. Tamkevičiene 1988, Firkins and Holdich 1993, Holdich et al. 1995, Gil-Sànchez and Alba-Ter-cedor 2002, Paglianti and Gherardi 2004, Hudina et al. 2011). However, our knowledge about the role of crayfish species in the ecosystem in regard to energy flux and nutrient cycling through their metabolic activity is very limited. Crayfish are key energy transformers among the different trophic levels since animals show omnivorous feeding character. They make major sources of energy (detrital material, decaying wood, dead organisms and periphyton) available to higher trophic levels at a more rapid rate than any other consumers (Momot et al. 1978). An estimate of crayfish metabolic activity could provide useful information for ecophysiological studies that deal with energy flow through ecosystems with crayfish populations. However, many crayfish species are threatened with population decline or extinction (Taylor 2002), so the number of experimental animals that can be taken from the wild is usually limited. Therefore, an elaboration of a new method that could provide samples of crayfish without killing them, or even without taking them from the wild, would be revolutionary. Most comparative physiological studies on crayfish have been conducted on adult specimens (e.g. Demers et al. 2006, Styrishave et al. 2007). However, the effect of environmental pollution in the early development stages could be very different. For the study of this, the survey of the intra- and inter-specific, and size-dependent physiological reactions is very necessary. Knowledge of the metabolic activity of different sized crayfish would serve as a basis for such other ecophysio-logical population studies, in which the metabolic activity of entire crayfish populations has to be estimated, taking into account age and size. Most studies on the relationship between metabolic rate and body size in crustaceans have used respiration rate as a measure of metabolic activity (Buikema 1972, Ivleva 1980, Wheatly 1989, Glazier 1991, Marshall et al. 2003), which is impractical when dealing with large and endangered crayfish species. However, if organisms are basically similar in body size, size-related changes in most biochemical and physiological processes should parallel the scaling of metabolism (Peters 1983). Therefore, the use of enzyme activity to assess metabolic rates appears to be a better approach, as demonstrated in studies on several crustacean species which also took into account the body size effects (Berges and Ballantyne 1991, Berges et al. 1990, 1993, Simčič and Brancelj 2003). The test of the enzymatic respiratory electron transport system (ETS) activity is a useful tool for estimating metabolic potential in aquatic organisms, since the result indicates the amount of oxygen would be consumed if all enzymes functioned maximally (Muskó et al. 1995). The method is simple, rapid and sensitive and has been used extensively on zooplankton (e.g. Bamstedt 1980, 1988, Borgmann 1978, James 1987, G.-Tóth and Drits 1991, G.-Tóth et al. 1995a; Simčič and Brancelj 1997, 2004), various amphipod and isopod species (e.g. Muskó et al. 1995, Simčič and Brancelj 2006, 2007, Simčič et al. 2005, 2010, Mezek et al. 2010) and fish (G.-Tóth et al. 1995b). Nevertheless, it has been rarely used in decapods. Borgmann ( 1977) studied ETS activities in various tissues of the crayfish Orconectes pro- pinquus (Girard, 1852), while Simčič et al. (2012) used whole animal homogenization to measure the ETS activity of the whole animal (ETSwhole ) in the noble crayfish Astacus astacus (Linnaeus, 1758) in order to relate it to oxygen consumption. To avoid killing animals and to circumvent the impractical whole body homogenization procedure in further crayfish studies, Simčič et al. (2012) proposed a new approach to estimate the metabolic activity of whole crayfish based on measuring ETS activity of a leg (ETSleg). The leg can easily be removed in the field and regenerates afterwards without harmful effects on the animal. The method was however tested only for a single species, while it is essential to know whether a general model could be established and used for different crayfish species. The four species included in the present study belong to two crayfish families with natural distribution ranges in Europe (EU) and North America (NA): Astacidae (EU: Astacus astacus (Linnaeus, 1758), Austropotamobius torrentium (Schrank, 1803) NA: Pacifastacus leniusculus (Dana, 1852)), and Cambaridae (NA: Orconectes limosus (Rafinesque, 1817)). All species are distributed in lentic and lotic freshwater ecosystems in temperate regions (Holdich 2002). The aims of the present study were to determine size scaling of the relationship between ETSwhole and ETSleg in two crayfish species, in order to establish a general model that can be used for estimating the metabolic activity of a whole crayfish on the basis of measured ETSleg. The proposed model has been tested on additional two crayfish species. Material and Methods Four crayfish species were included in the study. Two species are indigenous to Europe, Astacus astacus (AA; n = 35), Austropotamobius torrentium (AT; n = 6), and two to North America, Orconectes limosus (OL; n = 12) and Pacifastacus leniusculus (PL; n = 5). One European (AA) and one North American (OL) species were used for model construction and the other two for validation of the models. All specimens used in the laboratory tests were sampled in natural water bodies in Slovenia and Italy in 2009 and 2010. Protected and endangered species were collected under special licence No. 35601-135/2010-10 (issued by the Slovenian Environment Agency) in limited numbers. Live crayfish were transported to the laboratory in thermo-isolated bags in order to reduce stress effects. In the laboratory all specimens were maintained, prior to use, in aerated dechlorinated tap water for three weeks at 10 °C and photo-period light:dark = 16:8 hours. They were fed a commercial food (Sera crabs natural) ad libitum. Before measurements, the animals were weighed to the nearest 0.1 mg. Electron transport system (ETS) activity was measured using the method originally proposed by Packard (1971) and improved by G.-Tóth (1999). The third walking leg or the whole crayfish was homogenized in liquid nitrogen using a mortar. A weighed amount (50-90 mg wet mass) was sonicated in 4 ml of ice-cold homogenization buffer (0.1 M sodium phosphate buffer pH = 8.4; 75 цМ MgSO4; 0.15% (w/v) polyvinyl pyrrolidone; 0.2% (v/v) Triton-X-100) for 20 sec (4710; Cole-Parmer) and centrifuged at 8500 x g for 4 min at 0 °C (Centrifuge Sigma). Three 0.5 ml samples from each homogenate were incubated for 30 min at 10 °C in 1.5 ml substrate solution (0.1 M sodium phosphate buffer pH = 8.4; 1.7 mM NADH; 0.25 mM NADPH; 0.2% (v/v) Triton-X-100) with 0.5 ml 2.5 mM 2-(p-iodophenyl)-3-(p-nitrophenyl)-5-phenyl tetrazolium chloride (INT) solution. The reaction was ended by addition of 0.5 ml of stopping solution (formalin: H3PO4 conc = 1:1 v/v). Blanks (1.5 ml substrate solution and 0.5 ml INT solution) were incubated and treated as for the samples, followed by addition of 0.5 ml of homogenate. Formazan production was determined spectrophotometrically from the absorbance of the sample at 490 nm against the control blank within 10 min of stopping the reaction (WTW PhotoLabSpektral). ETS activity was calculated according to Kenner and Ahmed (1975). Since body protein content measured in homogenates of whole animals is generally a constant fraction of organism composition and therefore often used in size scaling studies (Berges and Ballantyne 1991), we measured the protein concentration in whole body (PROTwhole) and legs (PROTleg), using a commercial BCA Protein assay kit (Pierce). 0.1 ml of homogenate, from either the leg or the whole body, was pipetted into a test tube with 2.0 ml of reagent and mixed well. Samples were incubated at 37 °C for 30 min. Absorbance was measured at 562 nm against the blank (distilled water) within 10 min and the protein concentration determined from a calibration curve constructed with bovine albumin as standard. For model construction we used two species, one indigenous (AA) and one non-indigenous (OL), for each of which we obtained a sufficiently large amount of data to distinguish between general and species-specific patterns. The effect of sex on body mass, ETSwhole, ETS^, PROTwhole, PROT^, and on the ratio ETSwhole/ETSleg was tested with separate paired t-tests for each species. Since there was no significant effect (P > 0.05), both sets of data were pooled in further analyses. Least-squares regression analyses were performed to establish the mass scaling of the ETS activity and protein content variables. Regressions were performed with a power model of the form y = axb, where y is ETS activity or protein content, x is wet body mass of a crayfish and a and b are regression coefficients. We tested whether a single mass scaling equation can be used for both AA and OL. The improved fit obtained by taking into account species identity was evaluated by multiple regression in which species identity was added with a binary dummy variable A. The variables were In transformed as necessary to fit a linear regression model: Iny = /n(a)+b/nx+cA+d(A*/nx). The coefficient c tests the significance of the difference in intercepts while d tests for the difference in slopes. Separate models were deemed necessary if either slope or intercept, or both, differed significantly between species. If not, a joint equation was calculated from the pooled AA and OL data. The statistically significant equations are shown in the Figures 1-4. Pearson correlation between ETS activity and protein content was calculated separately for whole bodies and for legs. The relationship between ETS activity and protein content was established, by power regression, separately for whole bodies and for legs. The difference between the regression models for AA and OL was tested as for mass scaling. We tested two models for predicting ETSwhole from ETSleg. The first model involves the use of the ratio ETSwhole/ETSleg, estimated from wet mass (WW). ETSwhole is thus obtained by multiplying the ETSleg value with the estimated ratio. The final model is: ETSwhole = aWWb ETS, leg (model 1) The alternative model predicts ETSwhole directly from ETSleg: ETSwhole = aETSleg b (model 2) Body mass was excluded from model 2, since it did not contribute significantly to the fit (P > 0.05). We tested for the difference between the models for AA and OL by including the dummy variable A. Both models ultimately contained the same number of parameters to be estimated, so they could be compared on the basis of the adjusted r2 values. Finally we evaluated the applicability of the joint AA and OL models for two other crayfish species (PL, AT). Values of ETSwhole measured in these species were compared with those predicted from the two joint AA and OL models with paired t-tests. The same method was used to test whether mass scaling equations for protein content and ETS activity constructed with AA and OL data function also for the other two crayfish species. All statistical analyses were conducted in SPSS 13.0. Results Size-structured samples (based on their body mass) of 58 crayfish specimens were used in laboratory experiments. Models were constructed on the basis of measurements on 47 specimens (AA, OL), and later evaluated on the basis of 11 specimens (AT, PL). The mass of specimens did not differ significantly between sexes in AA and OL (P > 0.05), so differences in investigated parameters between sexes were tested using t-test. Since sex had no effect on values of ETSwhole, ETSleg, PROTwhole, PROTleg, or on the ratio ETSwhole/ETSleg in AA and OL (t-tests, P > 0.05), the data for both sexes were pooled for further analysis. Values of mass scaling of ETSwhole for AA (0. 877) and OL (0. 907) did not differ significantly (Table 1), so a single equation for pooled data was established (Fig. 1a). Insignificant differences were also observed between the two species in the regression of ETSleg and wet mass (Table 1; Fig. 1b). Both ETSwhole and ETSleg decreased with increasing wet mass. The regression models for mass scaling of PROTwhole for AA and OL differed significantly (Table 1). PROTwhole of AA was independent of Figure 1: The relationship between wet mass and (a) electron transport system (ETS) activity of whole body and (b) ETS activity of a leg in different crayfish species (OL - Orconectes limosus, AA - Astacus astacus, PL - Pacifastacus leniusculus, AT -Austropotamobius torrentium). The function was derived based on AA and OL data only. Slika 1: Razmerje med svežo maso in (a) aktivnostjo elektronskega transportnega sistema (ETS) celega telesa in (b) aktivnostjo ETS noge za različne vrste potočnih rakov (OL - Orconectes limosus, AA - Astacus astacus, PL - Pacifastacus leniusculus, AT - Austropotamobius torrentium). Funkcija je bila narejena le na osnovi podatkov za AA in OL. (а) /о 6Ü CD I 5C (J ::: 4G P r- b; Щ ™зо о CD Ö 20 10 0 №3 12n s 100 Ш iTo > 8Q * AA lO о OL 0 о * Power (OL) * • • \ о О * * ' *• * * О * * * '"-OS о* о* • * • * -f.* • • у = 60.832х"с-194 г2 = 0.373 P <0.05 20 40 Wet mass (g) SO гл о п, & £ 60 -40 - 20 -0 * AA * О OL о • Power (AA&OL) V # л* oV 0 0 * л • * VVr У _ w * * * * = 93.S2Sx"D-321 r*== 0.753 P <0.001 * * 0 10 20 30 40 Wet mass (al 50 60 70 Figure 2: The relationship between wet mass and (a) protein content of whole body and (b) protein content of a leg in Orconectes limosus (OL) and Astacus astacus (AA). Slika 2: Razmerje med svežo maso in (a) vsebnostjo proteinov v celem telesu in (b) vsebnostjo proteinov v celem telesu pri Orconectes limosus (OL) in Astacus astacus (AA). body mass (r2 = 0.005, df = 28, P > 0.05), whereas that of OL decreased with body mass (Fig. 2a). On the other hand, PROTleg in AA and in OL decreased similarly with wet mass according to a power law (Fig. 2b). ETS activity correlated significantly with PROTwhole (r = 0.788, n = 42, P < 0.001) and in PROTleg (r = 0.870, n = 42, P < 0.001). Similar scaling exponents b were observed for AA and 1.2 сл Lu □i o.o - LO liJ 0.6 0.4 - 0.2 - у = 0.240xD-2C8 f = 0.447 P <0.001 A Ж * *• • * * ______ о • * _^ ^--- * A # ф * ** Л * AA £ш * О OL * Л PL Ж AT -Power (АА&0 L) 20 40 60 80 100 120 140 Wet mass Щ Figure 3: The relationship between wet mass and the ratio of electron transport system (ETS) activity of whole body to that of a leg, in different crayfish species (OL - Orconectes limosus, AA -Astacus astacus, PL - Pacifastacus leniusculus, AT -Austropotamobius torrentium). The function was derived based on AA and OL data only. Slika 3: Odnos med svežo maso in razmerjem med aktivnostjo elektronskega transportnega sistema (ETS) celega raka in noge pri različnih vrstah potočnih rakov (OL - Orconectes limosus, AA - Astacus astacus, PL - Pacifastacus leniusculus, AT - Austropotamobius torrentium). Funkcija je bila narejena le na osnovi podatkov za AA in OL. OL (0.973 and 1.117) when the ETSwhole was expressed in relation to protein mass of crayfish with a common scaling exponent of 0.961. The ratio ETSwhole/ETSleg scaled similarly in AA and OL (Table 1) and showed a significant positive correlation with wet mass (Fig. 3). The common model (model 1), predicting ETSwhole from ETSleg with the help of this estimated mass ratio, explained 45% of the variation in whole body ETS activity of AA and OL: ETSwhole = 0.240 WW0208 ETSleg (model 1) ETSwhole was estimated directly from ETSleg using a power equation (Fig. 4). The regression coefficients for AA and OL did not differ significantly, so a common model (model 2) was constructed: ETSwhole = 8.498 ETSleg0511 (model 2) This model explained 51% of the variance in ETSwhole of OL and AA. Mass scaling of ETSwhole was similar in all the crayfish species in this study. The observed values of ETSwhole of PL (t-test, t = 2.388, df = 4, P > 0.05), and AT (t-test, t = 0.305, df = 5, P > 0.05) did not differ significantly from those predicted with the mass scaling equation established from the AA and OL data. Similar results were obtained for ETSleg, except for the ETSleg of AT, where the observed and predicted values differed significantly (t-test, t = 4.270, df = 5, P < 0.01), indicating different mass scaling of ETSleg in AT from those in OL and AA. The observed values of ETSwhole of PL did not differ significantly from the predicted values obtained from model 1 (t-test, t = 0.062, df = 4, P > 0.05). However, the predicted values of ETSwhole of AT (t-test, t = 2.768, df = 5, P < 0.05) differed significantly from the observed values when model 1 was used. Figure 4: The relationship between electron transport system (ETS) activity in a leg and that of whole body in different crayfish species (OL - Orconectes limosus, AA - Astacus astacus, PL - Pacifastacus leniusculus, AT - Austropotamobius torrentium). The function was derived based on AA and OL data only. Slika 4: Razmerje med aktivnostjo elektronskega transportnega sistema (ETS) noge in celega raka pri različnih vrstah potočnih rakov (OL - Orconectes limosus, AA - Astacus astacus, PL - Pacifastacus leniusculus, AT - Austropotamobius torrentium). Funkcija je bila narejena le na osnovi podatkov za AA in OL. Table 1: Results of statistical testing (for details see methods) of the differences in intercept and slope between Orconectes limosus and Astacus astacus for the specified relationships (ETSwhole - whole body ETS activity, ETSleg - leg ETS activity, PROTwhole - whole body protein content, PROTleg - protein content of a leg). Tabela 1: Rezultati statističnega testiranja (za podrobnosti glej metode) razlik v odseku in naklonu med vrstama Orconectes limosus in Astacus astacus za podana razmerja (ETSwhole - aktivnost ETS celega telesa, ETSleg - aktivnost ETS noge, PROTwhole - vsebnost proteinov v celem telesu, PROTleg - vsebnost proteinov v nogi). Intercept Slope Relationship t P t ETSwhole - wet mass 0.655 0.516 0.218 ETSleg - wet mass 0.339 0.736 0.574 PROTwhole - wet mass 2.758 0.009 2.081 PROTleg - wet mass 0.068 0.946 0.366 ETSwhole/ETSleg - wet mass 1.200 0.237 0.426 ETSwhole - ETSleg 0.720 0.476 0.835 Predicted values of ETSwhoie obtained with model 2 for PL (t-test, t = 1.248, df = 4, P > 0.05), and AT (t-test, t = 1.840, df = 5, P > 0.05) did not differ significantly from the measured ETSwhole. Discussion It seems that estimating the metabolic activity of crayfish from a single leg only is an adjustable method and that the model for the relationship between ETSwhole and ETSleg can be applied to different crayfish species. The relation between ETS activity and wet mass is in agreement with the findings in previous investigations that ETS activity varies with body mass according to a power law (Cammen et al. 1990, Muskó et al. 1995, Simčič and Brancelj 1997, 2003; Simčič et al. 2012). In two crustaceans, Chirocephalus croaticus (Steuer, 1899) and Gammarus fossa-rum Koch, 1835, it was shown that ETS activity is related to body size in a manner typical of a metabolic function with the scaling exponent b being 0.787 and 0.651, respectively (Simčič and Brancelj 2000, 2003). Since mass scaling of ETS activity did not differ significantly between AA and OL, a common equation for pooled data was proposed. The established b-value was in the range reported for the metabolic rate of crustaceans in general (Wolvekamp and Waterman 1960, Ivleva 1980). Comparison of intra- versus inter-specific b exponents for oxygen consumption in aquatic crustaceans showed that the intraspecific slope was approximately 0.1 less than the slope of the overall collected data (Wheatly 1989). In crayfish we have found even smaller differences between intra- and inter-specific b exponents for ETS activity. Berges and Ballantyne (1991) reported that intra- and inter-specific exponents for whole body maximal enzyme activities in aquatic crustaceans, i.e. Macrobrachium rosenbergii (De Man, 1879), Artemiafranciscana (Kellogg, 1906) and Daphnia magna (Straus, 1820), were similar for enzymes such as citrate synthase, but significant differences between species were found for enzymes associated with pathways other than aerobic metabolism. However, several authors have assigned the decrease in mass-specific metabolic activity to an increasing proportion of metabolically inert mass as the animals grow (Glazier 1991, Simčič and Brancelj 2003). Thus, the proportion of metabolically inactive tissue differs with species and varies with size and developmental stage of the same species. The protein content per gram of whole crayfish mass did not change with increasing wet mass in AA, but decreased significantly in OL (Fig. 2a). However, when the ETS activity was expressed in relation to protein mass, similar exponents were observed for AA and OL. The common scaling exponent for ETS activity of the two species was 0.961, similar to the exponents reported for decapod species M. rosenbergii for a variety of enzymes, where whole animals were homogenized and the enzyme activity was expressed in relation to protein mass (Berges and Ballantyne 1991). A closely similar b-value (0.955) was found for the relationship between the amount of protoplasm and ETS activity per individual in G. fossarum (Simčič and Brancelj 2003). Berges and Ballantyne (1991) reported that an exponent close to 1.0 is characteristic of enzymes capable of functioning catabolically or anabolically. Moreover, for larger crustaceans, such as M. rosenbergii, that rely increasingly on anaerobic metabolism as body size increases, scaling exponents closer to 1.0 are expected for the enzymes that function in both anaerobic and aerobic processes. The exponent close to 1.0 obtained for ETSwhole in the present study was in accord with the findings of Berges and Ballantyne (1991), since ETS activity measures both aerobic and anaerobic metabolism (Packard 1985). ETSwhole and ETSleg exhibited different mass scaling exponents (Fig. 1). This was expected, because different tissues contribute to the sample for ETS measurements in the two cases. Borgmann (1977) found that ETS activities differed in the various tissues of the crayfish Orconectes propinquus. Thus, ETSwhole reflects the activity of the mixture of the large number of different metabolically active tissues, body storage materials and exoskeleton, while in a leg sample muscle tissue and exoskeleton material predominate. Moreover, the proportion of protein material in a leg decreased up to 10 g of body mass, while in larger animals it was relatively constant (Fig. 2b). The decrease of PROTleg with similar exponents in the two species could mean a lower variability in PROTleg than in PROTwhole. The reason for higher variability in the relation between protein content and body mass probably lies in variable amounts of different storage materials and other metabolically inert tissues in the crayfish during their inter-annual life history, as well as in their age. However, the ETSwhole and ETSleg correlated well with their protein content, indicating that the latter plays a key role in ETS activity. Since ETSwhole and ETSleg scaled differently with body mass (Fig. 1), the ratio ETSwhole/ETSleg showed a significant, positive power relationship with body mass (Fig. 3). Similar scaling of the ETSwhole/ETSleg ratio in AA and in OL indicates a greater influence of body size than of species-specific properties on the relationship between whole body and leg metabolic potential. Therefore estimation of the whole crayfish metabolic potential can be estimated from ETSleg using a common ratio for both species, but this ratio depends on body mass. Thus, calculation of the ETSwhole/ETSleg ratio for a crayfish of a given body mass on the basis of the equation in Fig. 3 provides a factor that can be used for the estimation of ETSwhole from ETSleg (model 1 ). The utility of this model for estimating ETSwhole in different crayfish species was tested by comparing observed and predicted values in two crayfish species not used in the model construction. For PL, observed did not differ from predicted values, but the ETSwhole of AT was underestimated by this model (Fig. 3). The reason probably lies in the different developmental stage of 2-5 g AT from that of other crayfish of this size. Wheatly (1989) reported that the comparison of organisms at different developmental stages is problematic due to the different physiology of immature and adult individuals. Small-sized crayfish species such as AT attain their maturity at a mass of 2 to 5 g, while large-sized species become mature at a mass of more than 5 g (Souty-Grosset et al. 2006). It means that the specimens of AT were actually in a mature developmental stage, but that the same sized individuals of other species were still immature, and therefore possessed different physiological characteristics. Berges et al. (1990) found that the activity per unit mass of the primary anabolic enzyme nucleoside diphosphate kinase (NDPK) decreases with size and that enzyme scaling is affected by differences in growth rate. Thus, different development stages probably contribute to the relatively low ETSleg in AT. To minimize the potential effect of different development stages on estimated ETS activity, we explored a model relating ETSwhole directly to ETSleg. A common model (model 2) was constructed for both species, since the species-specific models did not differ significantly (Fig. 4). High variability in physiological and biochemical variables, especially due to different moulting stages, reproduction cycle and fitness, resulted in scattered data and, consequently, low coefficients of determination. Nevertheless, the results confirmed our expectation that the metabolic potentials of a whole crayfish and of a leg scale similarly in different crayfish species, since the whole body ETS activity estimated from a leg did not differ significantly from the observed value in all species investigated. The results of this study suggested that the metabolic potential of whole crayfish could be estimated on the basis of the ETS activity measurement in a single leg, using a general model. Due to the non-lethal approach, the new method could allow larger sample sizes to be incorporated into metabolic activity studies on crayfish, which is essential in conducting studies on multi-population and interspecific levels. Direct measurement of respiration is time-consuming and impractical and may subject the animal to stress before and during measurement. The slow response of ETS activity to short-term variations in environmental factors or stress makes the method superior to direct respiratory measurements on incubated animals (Barnstedt 1980). Moreover, the results of the present study also revealed that the two species-specific models did not provide significantly better estimates of whole crayfish metabolic activity than the joint ones. Model 2, by which ETSwhole was related directly to ETSleg, seems the most appropriate for an approximate estimation of whole body metabolic activity in all tested species. Further studies taking into account more species and larger samples could contribute to a clearer picture of the generality of model use. The proposed model, in which only population density, size structure (as body mass) and ETSleg are needed for estimating whole population metabolic activity, could be applied, in particular, for the estimation of metabolic activity in endangered and rare crayfish species. Acknowledgements We thank Andrej Kapla and Martina Jaklič for assistance, Dr. Bruno Maiolini for providing the experimental animals and Dr. Roger Pain for linguistic improvement of the manuscript. We acknowledge the financial support of the Slovenian Research Agency (Project L1-2169). References Bamstedt, U. 1980. ETS activity as an estimator of respiratory rate of zooplankton populations. The significance of variations in environmental factors. J. Exp. Mar. Biol. Ecol., 42, 267-283. Bamstedt, U. 1988. Ecological significance of individual variability in copepod bioenergetics. Hydro-biologia, 167/168, 43-59. Berges, J. A., Ballantyne, J. S., 1991. Size scaling of whole-body maximal enzyme activities in aquatic crustaceans. Can. J. Fish Aquat. Sci., 48, 2385-2394. Berges, J. A., Roff, J. C., Ballantyne, J. S., 1990. Relationship between body size, growth rate, and maximal enzyme activities in the brine shrimp, Artemia franciscana. Biol. Bull., 179, 287-296. Berges, J. A., Roff, J. C., Ballantyne, J. S., 1993. Enzymatic indices of respiration and ammonia excretion: relationships to body sizes and food levels. J. Plankton Res., 15, 239-254. Borgmann, U., 1977. Electron transport system activity in Daphnia and crayfish. Can. J. Zool., 55, 847-854. Borgmann, U., 1978. The effect of temperature and body size on electron transport system activity in freshwater zooplankton. J. Zool., 56, 634-642. Buikema, A. L. Jr., 1972. Oxygen consumption of the cladoceran, Daphniapulex, as a function of body size, light and light acclimatization. Comp. Biochem. Physiol., 42A, 877-888. Covich, A. P., Palmer, M. A., Crowl, T. A., 1999. The role of benthic invertebrate species in freshwater ecosystems. BioScience, 49, 119-127. Cammen, L. M., Corwin, S., Christensen, J. P., 1990. Electron transport system (ETS) activity as a measure of benthic macrofaunal metabolism. Mar. Ecol. Prog. Ser., 65, 171-182. Demers, A., Souty-Grosset, C., Trouilhe, M.-C., Füreder, L., Renai, B., Gherardi, F., 2006. Tolerance of three European native species of crayfish to hypoxia. Hydrobiologia, 560, 425-432. Firkins, I., Holdich, D. M., 1993. Thermal studies with three species of freshwater crayfish. Freshwater Crayfish, 9, 241-248. Gil-Sanchez, J. M., Alba-Tercedor, J., 2002. Ecology of the native and introduced crayfishes Austro-potamobius pallipes and Procambarus clarkii in southern Spain and implications for conservation of the native species. Biol. Conserv., 105, 75-80. Glazier, D. S., 1991. Separating the respiration rates of embryos and brooding females of Daphnia magna: Implications for the cost of brooding and the allometry of metabolic rate. Limnol. Oce-anogr., 36, 354-362. G.-Tóth, L., Drits, A. V., 1991. Respiratory energy loss of zooplankton in Lake Balaton (Hungary) estimated by ETS-activity measurements. Verh. Internat. Verein. Limnol., 24, 993-996. G.-Tóth, L., Carrillo, P., Cruz-Pizarro, L., 1995a. Respiratory electron transport system (ETS)-activity of the plankton and biofilm in the high-mountain lake La Caldera (Sierra Nevada, Spain). Arch. Hydrobiol, 135, 65-78. G.-Tóth, L., Szabó, M., Webb, D. J., 1995b. Adaptation of the tetrazolium reduction test for the measurement of the electron transport system (ETS)-activity during the embryonic development of medaka. J. Fish. Biol., 46, 835-844. G.-Tóth, L., 1999. Aktivität des Elektronentransportsystems, In: von Tümpling W. and G. Friedrich (eds): Biologische Gewässeruntersuchung. Methoden der Biologischen Wasseruntersuchung 2, Gustav Fischer Verl., Jena, Stuttgart, Lübeck, Ulm, pp. 465-473. Holdich, D. M. (ed.), 2002. Biology of freshwater crayfish. Blackwell Sciences Ltd, Oxford. Holdich, D. M., Reader, J. P., Rogers, W. D., Harlioglu, M., 1995. Interactions between three species of crayfish (Austropotamiuspallipes, Astacus leptodactylus and Pacifastacus leniusculus). Freshwater Crayfish, 10, 46-56. Holdich, D. M., Reynolds, J. D., Souty-Grosset, C., Sibley P. J., 2009. A review of the ever increasing threat to European crayfish from non-indigenous crayfish species. Knowl. and Manag. Aquat. Ec., 394-395, 11. Hudina, S., Galić, N., Roessink, I., Hock, K., 2011. Competitive interactions between co-occuring invaders: identifying asymmetries between two invasive crayfish species. Biol. Invasions, 13, 1791-1803. Ivleva, I. V., 1980. The dependence of crustacean respiration rate on body mass and habitat temperature. Int. Rev. Gesamten Hydrobiol., 65, 1-47. James, M. R., 1987. Respiratory rates in cladoceran Ceriodaphnia dubia in lake Rotiongaio, a mono-mictic lake. J. Plankton Res., 9, 573-578. Kenner, R. A., Ahmed, S. I., 1975. Measurements of electron transport activities in marine phyto-plankton. Marine Biology, 33, 119-127. Marshall, D. J., Perissionotto, R., Holley, J. F., 2003. Respiratory responses of the mysid Gastrosaccus brevifissura (Peracarida: Mysidacea), in relation to body size, temperature and salinity. Comp. Biochem. Physiol., 134A, 257-266. Mezek, T., Simčič, T., Michael, T. A., Brancelj, A., 2010. Effect of fasting on hypogean (Niphargus stygius) and epigean (Gammarus fossarum) amphipods: a laboratory study. Aquat. Ecol., 44, 397-408. Momot, W. T., Gowing, H., Jones, P. D., 1978. The dynamics of crayfish and their role in ecosystems. Am. Midl. Nat., 99, 10-35. Muskó, I. B., G.-Tóth, L., Szàbo, E., 1995. Respiration and respiratory electron transport system (ETS) activity of two amphipods: Corophium curvispinum G. O. Sars and Gammarus fossarum Koch. Pol. Arch. Hydrobiol.,42, 547-558. Nyström, P., 2002. Ecology, pp. 190-235. In: Holdich, D. M. (ed.): Biology of Freshwater Crayfish. Blackwell Science, Oxford, UK. Packard, T. T.,1971. The measurement of respiratory electron-transport activity in marine phyto-plankton. J. Mar. Res., 29, 235-244. Packard, T. T.,1985. Measurement of electron transport activity of microplankton, pp. 207-261. In, Jannas, H. and P. J. LeB. Williams (eds), Advances in aquatic microbiology 3, Academic Press, London. Paglianti, A., Gherardi, F., 2004. Combined effects of temperature and diet on growth and survival of young-of year crayfish: A comparison between indigenous and invasive species. J. Crustacean Biol., 24, 140-148. Peters, R. H., 1983. The ecological implications of body size. Cambridge University Press, Cambridge. Simčič, T., Brancelj A., 1997. Electron transport system (ETS) activity and respiration rate in five Daphnia species at different temperatures. Hydrobiologia, 360, 117-125. Simčič, T., Brancelj, A., 2000. Energy exploitation in Chirocephalus croaticus (Steuer, 1899) (Crustacea:Anostraca): survival strategy in an intermittent lake. Hydrobiologia, 437, 157-163. Simčič, T., Brancelj, A., 2003. Estimation of the proportion of metabolically active mass in the am-phipod Gammarus fossarum. Freshwater Biol., 48, 1093-1099. Simčič, T., Brancelj, A., 2004. Respiratory electron transport system (ETS) activity as an estimator of the thermal tolerance of two Daphnia hybrids. J. Plankton Res., 26, 525-534. Simčič, T., Lukančič, S., Brancelj, A., 2005. Comparative study of electron transport system activity and oxygen consumption of amphipods from caves and surface habitats. Freshwater Biol., 50, 494-501. Simčič, T., Brancelj, A., 2006. Effects of pH on electron transport system (ETS) activity and oxygen consumption in Gammarus fossarum, Asellus aquaticus and Niphargus sphagnicolus. Freshwater Biol., 51, 686-694. Simčič, T., Brancelj, A., 2007. The effect of light on oxygen consumption in two amphipod crustaceans - the hypogean Niphargus stygius and the epigean Gammarus fossarum. Mar. Fresh. Beh. Physiol., 40, 141-150. Simčič, T., Pajk, F., Brancelj, A., 2010. Electron transport system activity and oxygen consumption of two amphibious isopods, epigean Ligia italica Fabricius and hypogean Titanethes albus (Koch), in air and water. Mar. Fresh. Beh. Physiol., 43, 149-156. Simčič, T., Pajk, F., Vrezec, A., Brancelj, A., 2012. Size scaling of whole-body metabolic activity in the noble crayfish (Astacus astacus) estimated from measurements on a single leg. Freshwater Biol., 57, 39-48. Souty-Grosset, C., Holdich, D. M., Noel, P. Y., Reynolds, J. D., Haffner, P., 2006. Atlas of Crayfish in Europe. Muséum National d'Histoire Naturelle, Paris. Styrishave, B., Bojsen, B. H., Witthafft, H., Andersen, O., 2007. Diurnal variations in physiology and behaviour of the noble crayfish Astacus astacus and the signal crayfish Pacifastacus leniusculus. Mar. Fresh. Beh. Physiol., 40, 63-77. Tamkevičiene, E. A., 1988. Growth and development of juveniles of the native and introduced species of freshwater crayfish. Freshwater Crayfish, 7, 396-400. Taylor, C. A., 2002. Taxonomy and conservation of native crayfish stocks, pp. 236-257. In: Holdich, D. M. (ed.), Biology of freshwater crayfish. Blackwell Science, Oxford. Wheatly, M. G., 1989. Standard rate of O2 uptake and body size in the crayfish Pacifastacus leniusculus (Dana, 1852) (Decapoda: Astacidae): Intra- versus interspecific relations in crustaceans. J. Crustacean Biol., 9, 212-216. Wolvekamp, H. P. I., Waterman T. H., 1960. Respiration, pp. 35-100. In: Waterman, T. H. (ed.): Physiology of Crustacea, Vol 1. Academic Press, New York.