Introduction

Freshwater crayfish, mostly endemic to the continents of North America, Australia-Oceania, and Europe1, account for 1.71 million tons of global aquaculture production with a worth of 14.46 billion € as of 20182. Presently they contribute a negligible fraction in the global aquaculture scenario (~ 3.5% of total freshwater aquaculture production) but having great potential ahead. During the last half-decade alone (2013–2018), freshwater crayfish production, and its commercial valuation have tripled2. In terms of crayfish nutrition research, efforts have been quite limited compared to other commercially important crustaceans (like penaeids and palaemonids)3,4. Therefore, screening of novel feedstuffs, that too for crayfish, is somewhat ambiguous or problematic. A brief prologue in this regard is provided in the supplementary text. On the other hand, aquaculture nutrition research has focused on developing feed substitution strategies with a minimal supply of fishmeal and fish oil in recent times. One potential ingredient could be a microbial biomass meal from biofloc technology systems (BFT)4. BFTs are a modern, intensive aquaculture system that evolved from the classic ‘activated-sludge based sewage bioremediation’ in wastewater treatment plants. The system essentially operates on the rationale of maintaining an optimum C: N ratio (6:1 to 15:1) by daily purging with carbohydrate (carbon) source5,6. It is done to support the blooming of microbial biomass (flocs). These microbial flocculants, known as ‘bioflocs’ bioremediate the nitrogenous wastes generated by fish and uneaten feed into consumable microbial protein for cultured animals7,8. Although they are consumed by the fish or shrimp stock, the biofloc biomass (as measured in Imhoff cones) or total suspended solids (TSS) may often exceed the recommended values for fish (25–50 ml L−1; TSS up to 1000 mg L−1) and shrimp (10–15 ml L−1; TSS 400–600 mg L−1)—posing problems for the cultured animals7,9,10,11. It is advisable to drain part of the biofloc biomass daily through sedimentation or fractionation of biofloc system water10,12,13,14. Such thinning (filtering) of culture water generates a large amount of biofloc biomass as waste, quite frequently. This drained biofloc is often of limited use. In general, they can be used as an alternative to synthetic polymers for wastewater treatment15, fertilizer, or inoculum to start a new system16.

Our research intervenes in recycling this waste for aquatic animal nutrition. Since conventional protein sources in aquafeed (e.g., fishmeal) are becoming expensive and scarce, there has been a growing impetus in testing biofloc as an unconventional protein source for aquatic animals8,17,18,19. Few commercial floc meals are generically marketed under ‘single-cell protein (SCP)’ or ‘microbial protein’ category—Profloc (Nutrinsic), FeedKind (Calysta), and Novacq/OBM (Ridley, Maritech) with pricing (as of 2018) between 1.1–3.3 USD kg−117,18. One of these is listed in IAFFD (international aquaculture feed formulation database), with complete nutrient spectrum data, including essential amino acids20. So far, crayfish are not included in these mentioned researches. The novelty here is its potential use as a feedstuff (protein source) in the crayfish diet. In general, the protein (12–49%), lipid (0.5–12.5%), and ash (13–46%) contents in biofloc can vary substantially depending on several factors (reviewed by22). To the best of our knowledge, nutritional evaluation of biofloc as a feedstuff ingredient for artificial crayfish diets has not been done so far. Although rearing of crayfish in BFT system, where the animals co-fed on commercial feed pellets (primarily) and bioflocs suspended in the system, are recently being explored23,24. Our objective was to understand—(a) nutritional optima of freshwater crayfish from the available literature in the absence of centralized recommendations (see supplementary material); (b) growth trajectory and nutritional dependencies in crayfish (supplementary material); (c) response of red swamp crayfish to biofloc meal in their diet, in terms of nutrition, growth, and survivability; (d) the risk of heavy metals bioaccumulation or mineral stress in crayfish from feeding on biofloc, and; (e) evaluate nutritional strengths and bottlenecks associated with using biofloc meal in crayfish diet. The first two objectives (a and b) were rather a methodological and necessary step (placed in supplementary material) to the second part of our research related to the use of biofloc meal for crayfish (objectives c to e).

Results and discussion

Nutritional optima, growth trajectory, and nutritional dependencies of crayfish

Based on our meta-analyses, crayfish' optimum dietary nutritional requirement is tabulated as crayfish standards in Table 1. It is also compared with established standards of penaeid shrimps, often assumed as a template for most crustacean diets. Detailed information in this regard can be found in the supplementary material. In terms of crayfish growth trajectory, their thermal growth coefficient (TGC) may vary from 0.07–1 unit (interquartile range, IR 0.32–0.64 units). Results suggest any TGC in the range of 0.5–0.64 units may be regarded as ‘reasonably good growth’ in crayfish. Further insights into crayfish growth trajectory and its nutritional dependencies are presented in detail in the supplementary material. The information synthesized and approach used may serve as a template for future researchers exploring three less-established or unknown dimensions simultaneously (as in the present study)—novel feedstuff, optimum nutrition, and data-deficient (nutritionally) animals.

Table 1 Optimum dietary nutritional requirement of freshwater crayfish and its comparison with NRC (2011) standards for penaeid shrimps (usually adopted as status quo).

Growth response of crayfish to biofloc protein

Following a 9-week growth trial with graded BM levels in the diet, differential growth response by crayfish was realized (Fig. 1). Except for control and BM33 groups, crayfish' final body weight showed a significant deviation from the normal distribution. Further examining the skewness of final body weight distribution in BM66 and BM100 groups, it was apparent that these groups were dominated by runts (smaller sized individuals) with large size deviations from the handful of bigger individuals. The size heterogeneity showed a significant and negative correlation with BM inclusion in the diet (Pearson's 2-tailed r =  − 0.63, p < 0.05). Size heterogeneity in crayfish may aggravate community aggression25. However, the diet-driven size heterogeneity was not significantly correlated with mortality. The dietary treatments did not cause significant differences (p > 0.05) in survivability, confirmed by post-hoc analyses. Overall, the survivability remained > 70% through the experimental period in all groups (Table 2). It implies—BM does not pose a significant mortality risk to crayfish stocks irrespective of inclusion levels, but it has implications on growth (presented below).

Figure 1
figure 1

Body weight distribution in red swamp crayfish Procambarus clarkii fed graded level of biofloc meal (BM) in diets over 9 weeks of experimental duration. Measured on 20th, 38th and 63rd days post stocking. ‘Baseline’ indicates stocked stage-3 juveniles (0.007–0.008 g individual−1). Size heterogeneity (measured by coefficient of variance, CV) seems maximum and comparable in control (mean CV = 67%), BM33 (mean CV = 67.5%) and BM66 (mean CV = 63.4%) groups but significantly suppressed (p < 0.05) in BM100 (mean CV = 51%). BM100 showed poor size throughout the experiment.

Table 2 Response of the red swamp crayfish Procambarus clarkii (initial body weight 7–8 mg) under 9-week growth trial (21.8 °C) fed experimental diets. Values presented in interquartile range with mean ± standard deviation in parentheses.

The growth in terms of TGC, live-weight gain (LWG), and body weight (BW) were significantly depressed (p < 0.05) in the BM100 fed group. In contrast, the growth in control, BM33, and BM66 groups were higher with insignificant differences among them (p > 0.05) (Table 2, Fig. 1, 2). A statistically insignificant dampening of growth rate over time (p > 0.05) was observed in groups BM66 and BM100 (Fig. 2). At the end of culture (63 days), the realized TGC in crayfish fed on BM100 was on an average two times lower (p < 0.05) than the growth exhibited on control, BM33, or BM66 diets (Table 2, Fig. 2). In terms of feed utilization, the feed conversion ratio (FCR) and protein efficiency ratio (PER) were linearly related to increasing BM inclusion in the diet. The FCR increased with increasing share of BM in diet: FCR = 1.156 + 0.006 × BM (Adj. R2 0.95, p < 0.05). The PER decreased with an increasing BM inclusion: PER = 1.922 − 0.006 × BM (Adj. R2 0.95, p < 0.05). It means, for every 10% inclusion of BM, FCR increased by + 0.06 units, and PER decreased by − 0.066 units (Fig. 3). The results from the growth trial are summarized in Table 2, and the relationship of feed utilization parameters in response to BM inclusion is depicted in Fig. 3. Interestingly, the calculated FCR(s) of our respective diets, when multiplied with the dietary arginine content, seem to ‘hit the target’ of arginine requirements by crayfish (e.g., FCR of BM100 × Arginine in BM100 = Fulfillment of arginine requirement).

Figure 2
figure 2

Growth pattern (TGC: thermal growth coefficient) of red swamp crayfish Procambarus clarkii fed different experimental diets over 9 weeks. A dampening of growth over time gradually setting-in at higher BM inclusion in the crayfish diet (from BM66 to BM100). At the end of culture, BM100 resulted in twice less growth (p < 0.05) than achievable on other diets (control or BM33 and BM66—statistically comparable TGC).

Figure 3
figure 3

Feed utilization pattern (FCR in red and PER in blue) of red swamp crayfish Procambarus clarkii in response to the level of biofloc meal (indicated by BM.Inclusion, in %) in the diet. More feed is required per unit weight gain of crayfish with an increasing share of BM in the diet because protein utilization is lowered at higher BM inclusion.

As per the crayfish growth trajectory (quantified in the previous section), the group fed on 100% BM failed to show reasonably good growth. They were dominated by smaller-sized runts, poorest of the FCR and PER, but no significant mortality. Among the limited studies testing flocculated microbial meals in crustacean diets [reviewed in 8], BM inclusions were mostly up to 10–30% (of the total diet) or 30% (of fishmeal replacement). Good results in terms of growth were usually obtained at the maximum inclusion levels8. Like the present study, two previous studies had tested BM (on Litopenaeus vannamei) at a broader inclusion level from 17 to 84% of the total diet17,26. Despite different target species, the results seem close to that of the present study. Above 41–53% BM inclusion, the growth advantages were gradually lost17,26. Looking deeper into the aspects of our BM100-protein compared to control, BM33, or BM66-protein, the arginine seems to be a bottleneck for reasonably good growth (Tables 2, 4). Other EAAs, which could also be critical (e.g., methionine and lysine), were comparable-to-higher in BM100 than in other diets (Table 4). Although methionine and lysine levels in diets fell short of our formulated crayfish nutritional standard (Table 1), at least it fulfilled penaeid EAA standards of NRC4. It hints that NRC’s penaeid EAA standards cover well for most of the EAA requirements in crayfish, except for arginine (and tryptophan could not be judged). Arginine levels in BM (Table 4) neither fulfilled crayfish nor penaeid standards (Table 1).

Biofloc has been previously criticized for being partly deficient in arginine27,28,29. The arginine coefficient (proportion of total protein, in %) of biofloc meals, be it commercial ones like Novacq (2.38%19), FeedKind (2.54%20), or in the present study (2.73%) seem to have close resemblance (CV 5.5%). If we consider the mean arginine coefficient of BM from these data (2.55%) and tally it to fulfill the optimum arginine requirement of crayfish (minimum 1.8%), the crude protein level of such BM should be at least 70%. It is beyond the expected range of ordinary bioflocs22. BM harvested from high TSS systems (due to infrequent sedimentation or water exchange) can have lower protein content10. For example, the crude protein content of a biofloc can drop by − 34.5% if the TSS of the system is let to increase from ≤ 200 mg L−1 to 800–1000 mg L−110. As such, BM harvested from a low TSS system would have higher arginine (0.72%) compared to a high TSS system (0.47% arginine) (recalculated from10; using mean arginine coefficient = 2.55% of total protein). Even with aging biofloc, the content of arginine (also other EAAs) may decline. For example, from the 10th day to the 30th day of a biofloc culture, the arginine levels can decrease by 25–41% (recalculated from27). However, some specially produced commercial flocculated meals can have a high arginine coefficient (e.g., 5.3% of the protein in ProFloc17). Among all the EAAs, arginine content in red swamp crayfish seems maximum3,21,30,31, indicating a supposedly higher arginine demand in crayfish. The same is true for marbled crayfish Procambarus virginalis32. Arginine is perhaps the most limiting EAA in most crustacean diets and is required between 1.6–2.7% of diet33. Due to the poor activity of the urea cycle in crustaceans, arginine is indispensable for growth33,34. Arginine functions as a phosphagen in crustaceans, being the only amino acid providing amidino group for the synthesis of creatine—a major reserve of high-energy phosphate for ATP regeneration33,35.

Risk of heavy metals bioaccumulation or mineral stress from biofloc meal

The contents of heavy metals in BM were below the critical pollution limits. No critical limits were breached in the crayfish body that could qualify BM as a feedstuff capable of inducing unsafe heavy metal biomagnification, rendering them unfit for consumption. Content of Cd and Mn were mostly below the detection limits (Table 3). Except for mercury, hepatopancreas contained a higher amount of heavy metals (and minerals) than muscle. Hepatopancreas of crayfish, like most crustaceans, have been reported to be major storage of minerals, including heavy metals3,37. With increasing BM fraction in the diet, the concentration of Hg significantly increased in hepatopancreas (control → BM33 and BM66 → BM100; p < 0.05), while other metals did not show any significant trend (Table 3). Except for Cd, all metals were significantly higher (p < 0.05) in the hepatopancreas of BM100 fed crayfish compared to the control group. Such accumulation of heavy metal in hepatopancreas is capable of impairing metabolism in crayfish37. The concentration of Fe exhibits a rather ‘bell curve’ pattern, peaking at BM33 and receding thereafter, only in the muscle (Table 3). Cd and Zn did not exhibit any pattern as such. The heavy metal contents in crayfish and BM are given in Tables 3 and 4, respectively.

Table 3 Heavy metals and mineral content (mean ± SE; dry matter basis) in the tail muscle and hepatopancreas of red swamp crayfish Procambarus clarkii fed graded levels of biofloc meal.
Table 4 Proximate composition of biofloc meal, basal and treatment diets (dry matter basis).

Globally, the total ash content in biofloc may range between 13–46% (reviewed by22), also applicable in our case. The problem of high ash content in most biofloc, limiting its inclusion in diets (despite good protein content), has been briefly discussed in Sabry Neto et al.38. One previous study, which studied BM at a high enough inclusion level, attributed high ash and probable toxic effects of trace minerals to retarded growth in Litopenaeus vannamei fed > 60% BM in a diet26. Owing to high ash content in BM, mineral stress seems plausible in the present study as well (see Tables 1, 4). By mineral stress, we imply even if 10% of the ash or minerals from BM are digested by crayfish, it is potentially much higher ‘bioavailable minerals’ in the body than their optimum physiological limits. Information on this aspect have been limited for shrimps [reviewed in 39, 40] and none for crayfish3,36. In shrimps (Penaeus monodon, P. japonicus), retarded growth was observed when excessive mineral premixes were supplemented in a practical diet39, or more specifically, when trace minerals like Fe and Mn exceeded levels of 0.01% each in the diet40. The BM100 had all these factors (ash, Fe, and Mn) in excess (Table 4). Heavy metal stress could also be plausible. Any significant absorption of Hg in the body (presented above) is capable of impairing crayfish metabolism37, provoking hyper-osmoregulation in crustaceans41, with repercussions on aggravated energy expenditure42. Our metadata derived models show TGC in crayfish deteriorates at dietary ash levels > 14% (also in BM100), during which the retention of ash is merely < 10% of total dietary intake (see supplementary material and Fig S2, S3). Thus ≥ 90% of the ingested ash (exceeding physiological limits) are excreted through digestive and osmoregulatory (metabolic) pathways. It has its own energy cost, which could have been utilized for protein-sparing or growth42.

Recycling biofloc waste as a novel feedstuff for crayfish: Strengths and bottlenecks

Comparing the nutritional standards for crayfish with observed performance in growth trials, few strengths and bottlenecks of BM were realized (Tables 1 and 4). In terms of advantages: (a) BM has a high crude protein content (43.9%); (b) crude fiber content in BM (4.9%) was in the optimum range for crayfish, and; (c) BM is a rich supplier of minerals. However, there are more bottlenecks than limited advantages. BM has excessive total ash detrimental to crayfish growth, with probable manifestations on hyper-osmoregulation and energy expenditure (discussed above). A mediocre crude lipid content (4.5%) is another bottleneck for supplying non-protein energy. These, in combination, render the non-protein energy: protein ratio (NPE: P = 3.7 cal non-protein energy per 1 mg protein) in BM insufficient for effective protein sparing (≈growth). At such low NPE:P, the proteins are catabolized for meeting energy demand (even after oxidizing carbohydrates and lipids), rather than building biomass42. It is further compounded by arginine deficiency in BM (~ 14–20% less than an optimum requirement)—probably the most critical essential amino acid for crayfish (discussed above).

A retrospective evaluation of BM100 or BM (as a feedstuff for crayfish) applying our metadata derived ‘growth-retention models’ (supplementary Fig S3, S4, S6) could explain few nutrient utilization scenarios behind low growth in BM100. The ash, protein, and lipid retentions from BM should be less than 5%, 10%, and 3% of dietary intakes, respectively (predicted). For control, BM33, and BM66 diets, these retentions were well above the identified thresholds qualifying for reasonably good growth in crayfish (refer to supplementary material). Comprehensively, the retarded growth problem with solely feeding on biofloc biomass could be a synergistic effect of— (a) arginine deficiency, (b) mineral and heavy metal stress, and, (c) low non-protein energy to protein ratio.

Methods

Calculation of crayfish nutritional standards, growth trajectory, and its nutritional dependencies

In the absence of centralized nutrition recommendations for freshwater crayfish species, unlike other commercially important crustaceans (e.g., penaeid shrimps, see NRC4), available literature was meta-analyzed. Peer-reviewed and published articles (in English or at least with English abstract) were searched online (search engines: Web of Science, Scopus, and Google Scholar) using keywords like ‘growth trials’, ‘crayfish’, ‘nutrition’, ‘proximate composition’, ‘body composition’, ‘amino acids’, ‘heavy metals’, ‘optimum requirement’ were used in different combinations (depending on target information). Altogether 27 articles were sourced and data extracted for meta-analyses. Detailed methodology on each meta-analysis (i.e., formulation of nutritional standards, calculation of growth trajectory and feed utilization parameters, quantification of nutritional dependencies on growth) are provided in the supplementary material.

Collection of biofloc biomass

Biofloc biomass was obtained from a well-established indoor, freshwater biofloc system, stocked with Nile tilapia Oreochromis niloticus at a stocking density of 35 kg m−3. Commercial pellets (TILAPICO 3 mm, Coppens, The Netherlands) were used as standard feed for fish. Fish feed was given twice daily based on a feed amount equivalent to 2.5% of the fish body weight. Wheat flour (35.56% C; 2.38% N) served as a carbon source which was applied daily with feed (22.05% C; 7.07% N) in a ratio of 1:0.6 (feed: flour). Assuming a 30% retention of nutrients from feed to fish, the projected C: N ratio was ≈6:1. Such a low C: N ratio favored frequent harvest of young and N-rich wet biofloc biomass6,10 to be converted to dry matter for the ensuing experiment. Biofloc biomass was drained daily through a pump and a vortex separation device so that the suspended solids level stayed between 25 and 50 ml L−1 in the system. After separation, biofloc was filtered through a nylon screen (mesh size 60 μm) to drain the excess water. The filtrate was then dried at 80 °C to obtain a material of solid consistency. After obtaining enough dried biofloc, the samples were grounded by a hammer mill to yield finer particles and hereinafter referred to as the biofloc meal (BM).

Preparation of experimental feed

Commercial pellets (TILAPICO 3 mm, Coppens, The Netherlands) were used as the basal diet due to its similar protein content with our test ingredient (BM). The commercial ‘fish feed’ was chosen due to a lack of established ‘crayfish feeds’ in the market. Even the available ones appeared to be random feed mixtures targeted for ornamental crayfish keeping. Inclusion of BM by replacing basal diet was done on a weight by weight basis. All feeds were isonitrogenous. The graded inclusion levels were 0% (basal diet = control diet), 33% (67% basal + 33% BM; diet BM33), 66% (34% basal + 66% BM; diet BM66) and 100% (only BM; diet BM100). Feed pellets (pellet size 2 mm) were cold extruded, dried (12 h; 45 °C), vacuum sealed, and stored at 4 °C till further use. The diet samples were analyzed in an accredited third-party laboratory (AGRO-LA, spol. s.r.o., https://www.agrola.cz/zemedelske-a-potravinarske-sluzby/) employing analytical methods (ISO verified and certified protocols in the Czech Republic) for proximate composition, essential amino acids (EAAs; except tryptophan due to analytical error), heavy metals, and essential mineral contents. Detailed composition of basal diet, treatment diets and the biofloc meal are summarized in Table 4.

Crayfish keeping

A total of 120 juvenile red swamp crayfish (Procambarus clarkii; conservation status: least concern) having a mean weight of 7.8 ± 0.7 mg at the onset of exogenous feeding (developmental stage 3), were used as experimental animals (10 individuals per tank; 4 group x triplicate). The experiment lasting for nine weeks was conducted in a series of indoor glass aquaria (54 × 36 × 30 cm, volume 46 L) with aeration and attached to a recirculating aquaculture system. Two baked clay bricks (28.5 × 13.5 × 6.5 cm), each with 39 cross holes (26 and 13 holes with a profile of 1 × 3 cm and 1 × 1 cm, respectively), were placed in each aquarium to provide shelters/refugia for the stocked crayfish43. After three weeks, a block of joined polypropylene tubes containing five tubes (length 10 cm, inner diameter 35 mm) was added to each aquarium as an additional shelter for on-growing animals. The bases were represented by three longitudinally joined tubes with a further two tubes positioned pyramidal in the second layer44. Altogether, 12 tanks were used and subjected to stable indoor climatic conditions with natural photoperiod (12L:12D).

Growth trial and feed utilization parameters

Crayfish were fed twice a day to apparent satiation (roughly corresponding 5–6% of the body weight) with the abovementioned diets for nine weeks. Uneaten feed, feces, and other wastes were siphoned out manually every morning. Dissolved oxygen (7.9 ± 0.3 mg L−1), pH (7.6 ± 0.2), and temperature (21.8 ± 0.3 °C) were measured daily using Oxi 3205 and pH 720 m (WTW GmbH, Weilheim, Germany), respectively. Every three weeks, the body weight was measured using an electronic balance (lowest sensitivity 1 mg) and the number of survivors counted. The feed rationing was revised accordingly. Body weight measurements were taken before feeding. After the trial, final body weight and total length were recorded, including the number of survivors. The animals were not fed before the day of the final measurement.

The food conversion ratio (FCR, units), protein efficiency ratio (PER, units), and survivability (%) were determined for each diet following the formulas in Cortes-Jacinto et al.45. Live weight gain (LWG) was calculated applying the formula, LWG = final—initial weight (in mg)/ days reared. Coefficient of variance (CV) of body weight (standard deviation × 100 ∕ mean) was calculated as a measure of size heterogeneity. To eliminate statistical biasedness in the data due to hierarchical size distribution in crayfish groups, other measures of central dispersion like interquartile range (IR) and median were included besides the mean. The abovementioned parameters were calculated from the IR, median, and mean estimates of each treatment. All graphical models were generated using the ggplot2 package in R. Statistically significant differences (α level set at 0.05) in body weight, growth, and survivability of crayfish fed on different dietary treatments were tested. The grouped data were first subjected to a Shapiro–Wilk’s normality test; then following the p value, either one-way ANOVA with post-hoc Tukey HSD (parametric test), or, Kruskal–Wallis post-hoc Dunn’s test with Bonferroni correction (non-parametric test) was selected. The tests were performed using default commands in RStudio v1.2.5042.

Assessment of heavy metals risk from biofloc biomass

At the end of the experiment, tail muscle and hepatopancreas samples from representative crayfish of each group were collected and frozen (−20 °C). Selected heavy metals (Hg, Cd, Zn; following high bioaccumulation affinity realized in Kouba et al.46) and some additional minerals (Fe, Mn) were analyzed from these samples in the same accredited third-party laboratory. Body (muscle + hepatopancreas) heavy metal levels were compared with maximum permissible limits (Cd or Hg 0.5 mg kg−1 wet weight basis) given in the European Commission47 for aquatic meat products (in the context of safety for consumption). In the context of agricultural use safety (as fertilizers), the heavy metal content of biofloc meal was determined and compared with Czech EPA limits (Cd 5 mg kg−1, Hg 4 mg kg−1 dry matter basis) (Decree of Ministry of Environmental of the Czech Republic No. 437/2016 on the Code, 2016).

Ethics approval

All procedures performed in studies involving animals (Oreochromis niloticus and Procambarus clarkii) were in accordance with the ethical standards approved by the institutional ethics committee (Jihočeská univerzita v Českých Budějovicích Fakulta rybářství a ochrany vod).