Introduction

The phenomenon of cannibalism is observed in many animal groups, inter alia, in protozoa, rotifers, snails, copepods, mites, insects, spiders, fish, amphibians, birds, and mammals (Fox, 1975). Increasing evidence indicates that cannibalism is not only common but that it is important to the ecology of many species. Fox (1975) concluded that cannibalism is a natural behaviour in 147 species, while Polis (1981) cites reports on approximately 1300 species. Smith & Reay (1991) define cannibalism more precisely as an act of killing and consuming wholly or in large part an individual belonging to the same species regardless of developmental stage. Similarly, Garcia & Zaniboni-Filho (2006) defines cannibalism as a special type of predation that depends on killing individuals of the same species for partial or total consumption.

Among fish, the phenomenon of cannibalism has been identified unequivocally in the literature in 36 of 410 teleost fish families (Nelson, 1984). It is likely, however, that cannibalism in fish is more widespread than a literature review can suggest (Dominey & Blumer, 1984) since the phenomenon has been described in the literature mainly in species that are of interest to humans as resources, are used as models in behavioural studies, or are popular species in fishkeeping (Smith & Reay, 1991). Pereira et al. (2017), who reviewed more than 1000 articles on fish cannibalism, conclude that this phenomenon among marine fishes occurs most frequently among Gadidae (106 records, Atlantic cod, Gadus morhua dominant), Gobiidae (54 records, common goby, Pomatoschistus microps, dominant), and Merlucciidae (29 records, genus Merluccius dominant), while among freshwater fishes it occurs among Percidae (140 records, European perch, Perca fluviatilis, and yellow perch, P. flavescens, dominant), Salmonidae (75 records, Arctic char, Salvelinus alpinus, dominant), and Esocidae (54 records, Northern pike, Esox lucius, dominant).

Although cannibalism occurs in both early developmental stages and adult fish, most instances of cannibalistic behaviour refer to the consumption of larvae or juvenile individuals (Pereira et al., 2017). This is true for the pike analysed in the current study in which cannibalism is particularly widespread in the first year of life (Craig, 1996). Feeding habits in young pike change as they grow and three feeding stages occur during development (Hunt & Carbine, 1951). This feeding sequence progresses with pike growth from lower crustaceans to insects, and then to fish. It is supplemented, even at an early developmental stage (15–20 mm), by attacks on larval fish, including cannibalistic attacks (Hunt & Carbine, 1951; Morrow et al., 1997; Mamcarz et al., 1998; Žiliukienė & Žiliukas, 2006).

Cannibalism is divided into seven classes according to the developmental stage of the prey, the genetic relationship between the cannibal and its prey, and the age relationship between them (Smith & Reay, 1991). According to the classification by Smith & Reay (1991), pike belongs to the class of “sibling intracohort cannibalism of post-hatching stages.” This class is characterized by sibling cannibalism or individuals belonging to the same age class, which is typical of piscivorous species in larval and juvenile early life stages (Naumowicz et al., 2017). Bearing in mind the methodology used in the present research and size range of analysed fish, type I cannibalism is considered. This occurs when the size of the prey is large relative to that of the predator, it can be ingested tail first (or less frequently head first) and consumed (Cuff, 1980); this is in contrast to type II cannibalism, which is characterized by the cannibal immediately and completely ingesting the prey (Baras & Jobling, 2002).

The ability to eat siblings in early life stages before offspring disperse, which occurs in pike, is facilitated by the pike’s large gape and well-developed dentition shortly after hatching (Giles et al., 1986). Cannibalistic larvae usually have augmented dentition, large gape size, and increased jaw musculature (Folkvord, 1997). In both adult and early life stages of pike, the long build of the body and its articulated jaw and backward-facing hinged teeth prevent prey from escaping, and the combination of the propellant force of the anal, caudal, and dorsal fins helps to obtain high acceleration when attacking prey (Craig, 1996; Skov & Nilsson, 2018). Additionally, the underslung pairs of pectoral and pelvic fins enable precise manoeuvring to sight prey stereoscopically down paired snout grooves (Craig, 1996).

In the natural environment, cannibalism is a key mechanism for regulating pike abundance that has significant implications for pike year class recruitment (Craig & Kipling, 1983) and also for the structure of given populations (Skov & Nilsson, 2018). It is especially important for over-exploited populations, which can be found not only in freshwater systems but also in many coastal areas of the Baltic Sea (Nilsson et al., 2014; Larsson et al., 2015; Skov & Nilsson, 2018). Although the main reasons for population declines are over-exploration and the disappearance of suitable spawning grounds and low recruitment (Skov & Nilsson, 2018), other factors related to human activities in coastal areas (e.g. wind farms and underwater cables) should also be considered (Fey et al., 2019). When there is a lack of recruitment from natural spawning, local populations are often supported by stocking with pike larvae or fry obtained from recirculating aquaculture systems (RAS), which have been found to be very efficient in producing pike stocking material (Szczepkowski et al., 2012). Unfortunately, cannibalism among pike early life stages plays a significant role not only in natural but also in aquaculture conditions because it can cause total losses that exceed 90% of total mortality (Szczepkowski, 2009). The significance and potential consequences of the occurrence of cannibalism during larval and juvenile stages mean it is extraordinarily important to acquire an understanding of this phenomenon.

There are several methods for determining the condition and fitness of fish specimens, including methods that evaluate Fulton’s condition factor (Jin et al., 2015), morphology and histology (Theilacker, 1986), lipid and protein composition (Håkanson, 1989; McCormick & Molony, 1993), RNA to DNA ratios (Clemmesen, 1988; Folkvord et al., 1996; Peck et al., 2015), and otolith microstructure (Campana & Jones, 1992). Alternatively, fitness is reflected in developmental instability (Zakharov, 1992). Developmental instability can then be reflected in fluctuating asymmetry (Palmer, 1994), which is defined as a random deviation in a perfect bilateral system (Palmer, 1994). Assuming environmental stress causes developmental instabilities, high fluctuating asymmetry can indicate the lower condition of larvae that experienced unfavourable environments. The link between fish condition and fluctuating asymmetry has been studied in both adult and larval fish, and a number of measurements have been analysed, including otolith size and shape (Escós et al.,1995; Somarakis et al., 1997a, b; Fey et al., 2020). However, the results obtained in different studies are inconsistent. Presumably, whether fish larvae belong to cannibal or prey groups within a cannibalistic population depends on differences in both growth rates and fish condition. The question posed in the present study is whether there are also differences in developmental instability that are reflected in fluctuating asymmetry in otolith size between groups of cannibals and prey. Although otolith growth and daily otolith increment deposition have recently been analysed for larval and early juvenile northern pike (Fey et al., 2018), no information on otolith asymmetry has been published to date.

The aim of the present study was to determine for larval pike (9–32.5 mm SL) held under controlled conditions if the water temperature (10, 15, 22 °C) had an effect on the age and size of predators during the first cannibalistic attacks, the intensity of the cannibalism phenomenon, and the proportion of prey to predator body size. The methodology used chosen to analyse type I cannibalism (the cannibal ingests its prey by the tail or head first). Additionally, whether belonging to predator or prey groups was related to the presence of higher or lower organism bilateral instability was also verified by analysing otolith size fluctuating asymmetry.

Materials and methods

Experimental conditions

The experiments were performed in nine freshwater tanks with volumes of 55 L each combined into three RAS providing three temperatures (mean ± SD): 10 ± 0.6°C, 15 ± 0.4 °C, 22 ± 0.5 °C, with three replicates of each temperature (3 temperatures × 3 replicates). Each of the systems was equipped with UV sterilization and a filtering-deposit container with a volume of 70 L providing both mechanical and biological filtration. The volume of biological filtration was 20 L. In each of the tanks a green plastic net (mesh size 0.5 cm) formed into a roller (20 cm in length × 8 cm in diameter) was used as a substitute for bottom substrate. The colour of the tank bottom was grey, and the glass walls were not covered. The tanks were located in a completely dark room in which light conditions (~ 100 Lux on the water surface, 12 h light/12 h dark, with light on from 8.00 to 20.00) were regulated automatically. Water temperature was measured three times daily with electronic temperature readers. The oxygen level during the experiment was measured daily and was 7.8 ± 0.9 mg l−1 in all tanks (mean ± SD) throughout the experiment.

The rationale for choosing the particular temperatures was to cover the wide range of temperatures that larvae encounter in natural conditions depending on the type of water body and depth. Temperatures typically range between 10 and 15 °C, but they can exceed 22 °C in shallow flooded fields where the water depth is only several centimetres. Two temperatures (15 °C and 22 °C) correspond to conditions that can occur in aquaculture systems.

The total time of experiment was 45 days, but this differed depending on the temperature: 17 days at 22 °C (no fish left beyond this point), 37 days at 15 °C (no fish left beyond this point), and 34 days at 10 °C (no cannibalism occurred). On day 34, the fish from the 10 °C treatment were transferred to the 15 °C treatment for 11 days (group 10–15 °C). The rationale behind the transfer from 10 to 15 °C was to test if cannibalism was triggered at the higher temperature and after how many days (see results for the information about no cannibalism at 10 °C). The scheme of the experimental setup is presented in Fig. 1.

Fig. 1
figure 1

Scheme of the experimental setup. The transfer of larvae from 10 to 15 °C on day 34 was a consequence of no cannibalistic attacks at 10 °C

Fish larvae rearing

Eyed-egg pike embryo were transported to the laboratory from a commercial hatchery, and 3000 eggs obtained from one female were placed in separate bags for transportation. The number of eggs was estimated by the volumetric method using a glass cylinder with a calibration scale (Szczerbowski, 2008). Mass hatching occurred during the first hour of transportation, and the hatched larvae were placed in tanks (one bag in one tank). The unhatched eggs left in the bags were counted to estimate the number of larvae in each tank. Hatching success was within a range of 92–95% in all of the transportation bags. The temperature was set to 10 °C in each tank at the beginning of the experiment. Thereafter, the temperature was increased to 15 °C within 24 h in one of the systems and to 22 °C within 48 h in the other one, and it was constant thereafter. Following yolk-sac absorption, feeding was conducted ad libitum every 1.5 h during periods with light with specialist fodder, as follows: ArtEX 2 (0.15 mm) (Aller Aqua) for the first 4 days after yolk-sac absorption; Perla Larva Proactive 4.0 (0.3-0.5 mm) (Skretting AS) from day 3 after yolk-sac absorption. Dead fish, unconsumed fodder, and excrement were cleaned from the tank bottoms every morning. The number of dead fish in each tank was counted to calculate the current density of larvae in each tank as well as the cannibalism intensity index (see data analysis section).

Cannibalism analysis

During each feeding (i.e. six times daily, every 1.5 h during the light period), the occurrence of cannibalistic pairs, i.e. prey caught by a cannibal inside its mouth (Fig. 2), was monitored in all tanks simultaneously for 15 min periods by two observers. Because only cases when prey was inside mouth of the cannibal were counted as cannibalistic attacks, this approach can be considered to be the analysis of successful cannibalism. Each pair was collected, and after placing it in 800 mg l−1 MS-222 anaesthetic fluid (ACROS Organics™, Belgium), it was transferred into an individual vial (one pair per vial) of 96% alcohol to preserve the material for otolith asymmetry analysis. Although the approach described did not provide information on the absolute number of cannibalistic attacks, it permitted calculating the relative type I cannibalism intensity index that is comparable among the three temperatures. A total of 71 cannibalistic pairs was collected.

Fig. 2
figure 2

Picture of a cannibalistic pair of pike (Esox lucius) representing type I cannibalism (without immediate ingestion). The cannibal holds the prey, which is not much smaller, by the tail (a). Once the difference in size between cannibals and prey increased, the attacks from the side and head were also observed (b). Only instances in which the prey’s body was inside the cannibal’s mouth were considered to be cannibalistic attacks in this study. All “connected” pairs were collected and analysed

The alcohol in each vial was changed after two weeks since it has been demonstrated that alcohol storage can damage otoliths and affect the results of otolith size measurements (Fey, 2018), which would influence the results of otolith size asymmetry analysis performed for the current study. Cannibal and prey body sizes (standard length, SL) were measured with an electronic calliper to the nearest 0.1 mm a few weeks after collection. For the purpose of background description (i.e. larvae size, both prey and cannibals, in a given tank in a given temperature), 10 larvae were collected every three days and placed first in the anaesthetic fluid and then in alcohol. The SL of these larvae were measured to the nearest 0.1 mm within a few weeks. The size of the preserved specimens was corrected for alcohol shrinkage, and unpreserved size was calculated with a correction factor (+0.7 mm) (Greszkiewicz & Fey, 2018). The size range of the larvae collected from cannibalistic pairs was between 11.8 and 32.5 mm SL.

Otolith fluctuating asymmetry analysis

Otoliths (left and right sagitta) were extracted from each larva and placed in DePeX mounting medium (Electron Microscopy Sciences, Fort Washington, PA) on microscope slides. The size of the otoliths was measured with an image analysis system (Image-Pro Premier, Media Cybernetics, Inc., Rockville, MD) under an Eclipse 80i transmitted light microscope (Nikon Corp., Tokyo, Japan).

The fluctuating asymmetry of otoliths (left otolith—L, and right otolith—R) was analysed according to methods described by Fey & Hare (2008). First, it was verified that the (L–R) was not different from zero (one-sample t test to 0) and that the distribution of differences was normally distributed (Shapiro–Wilk W-test) and free from antisymmetric (bimodality or platykurtosis) and directional symmetry (skew). Linear regression analysis was then employed to test if the absolute difference between the left and right otolith size (|L − R|) depended on trait size (i.e. otolith size). The largest measurement (left or right otolith) was used as the trait size since Sullivan et al. (1993) demonstrated that using the mean value as trait size can introduce a spurious correlation and should be avoided. The analysis was performed for cannibals and prey separately. Since the |L − R| was not dependent on trait size, the mean |L − R| was used as the fluctuating asymmetry index (FA). The comparison was based on the n = 71 cannibals and n = 71 prey sample size.

Data analysis

The time (dph) of the first cannibalistic attack and the size of the cannibals in the attacks were compared among the temperatures. How cannibalism intensity changed over time (i.e. with changing fish density in tanks and fish size) and what the effect of temperature was on the phenomenon were analysed with the cannibalism intensity index CI (%) that was calculated for each sampling day:

$${\text{CI}}\; = \;N_{\text{canib}} /N_{\text{all}} \cdot 100$$

where Ncanib is number of cannibalistic pairs collected on a given day in a given tank and Nall is number of all larvae in the tank.

The number of larvae in a given tank on a given day was estimated by subtracting the number of dead fish found during morning cleaning on that day and during previous days from the number of fish in that tank at the beginning of the experiment (i.e. number of eggs in the transportation bag minus unhatched eggs).

The CI was plotted against SL (mm) and larvae density (n·L−1) separately for the three temperature groups. The growth rate of the cannibals and the prey was compared with slope comparison (ANCOVA) of linear functions fitted to SL at age data. The differences in the FA between cannibals and prey were evaluated using the t test. All those analyses were conducted after data were pooled among the three subsample tanks for each temperature to obtain a sufficient number of observations for this calculation. Still, it was beneficial to have subsamples from three tanks as the CI index calculations were performed for the specific conditions (density and larval size) the larvae encountered in each of the tanks. Statistica, version 12.0 (TIBCO Software, Inc., Palo Alto, CA) software was used for data analysis. Differences were considered statistically significant at P < 0.05 (α = 0.05).

Results

First cannibalistic attack—cannibal size and age

In total, 71 cannibal-prey pairs representing larva size ranges from 11.8 to 32.5 mm SL were collected during the experiment. Thus, successful attacks were evaluated without considering unsuccessful attacks or bites only. The first successful cannibalistic attack was observed 5 dph at 22 °C, 10 dph at 15 °C, and no attacks were observed at 10 °C (Fig. 3). Cannibal size at the first attack was 13.0 mm at 22 °C (Fig. 3a) and 13.4 mm at 15 °C (Fig. 3b). When the fish were transferred from 10 to 15 °C (34 dph, size range 14.4–19.4 mm), the cannibalistic attacks started after 5 days at a cannibal size of 19.8 mm (Fig. 3c). Fish from the full-size range present in the tanks at a given time participated in the cannibalistic phenomenon with cannibals representing approximately the upper half of the size range of all the larvae in the tanks and the prey representing the lower half of the size range (Fig. 3).

Fig. 3
figure 3

Standard length (SL, mm) of larval pike, Esox lucius, cannibals (circle and solid line), prey (diamond and dashed line), and other specimens present in the tanks (filled square) on individual days after hatching at three different temperatures: 22 °C (A), 15 °C (B), 10 °C (C). Each point represents an individual fish. The data were pooled for three tanks at each temperature

The first attacks at 15 and 22 °C were from the tail only (Fig. 2a). At a cannibal size of approximately 17 mm, attacks from the side and head (Fig. 2b) were also noted and they represented 10% (side) and 14% (head) of all the attacks among larvae from the 17–32 mm size range.

Prey to cannibal size ratio

The prey size to cannibal size proportion ranged, on average, from 95% for 13 mm fish to 80% for 27 mm fish, with a mean for all specimens of 87% (Fig. 4). The smallest differences in size between the prey and cannibals (approximately 0.2 mm) occurred in the cannibal size range from 13.4 (difference = 1.5%) to 18.4 mm (difference = 1.1%). There was no temperature effect on the analysed proportion (Fig. 4).

Fig. 4
figure 4

Percentage values of the proportion of prey size to cannibal size of larval pike, Esox lucius, at three different temperatures: 22 °C (squares), 15 °C (circles), 10–15 °C (diamonds; fish transferred from 10 to 15 °C after 34 days of the experiment). The range of the highest similarities (range of 95–100%) in sizes between cannibals and prey is marked with a dashed line

Growth rate of cannibals and prey

The growth rates of cannibals and prey was estimated from the slopes of regression lines fitted to SL at age data. The growth rate of larvae at 22 °C was 0.57 mm·day−1 for cannibals and 0.32 mm·day−1 for prey (Fig. 3a). At 15 °C, the growth rate of cannibals was 0.55 and of prey it was 0.43 mm·day−1 (Fig. 3b). The differences between the growth rates of cannibals and prey were statistically significant for both 15 and 22 °C (ANCOVA, slope: P < 0.05). Since the slopes were different, no comparison was made for the elevations of the regression lines. In the 10–15 °C group, growth was not calculated because of the narrow size range of the larvae. However, the cannibals were clearly larger at age (i.e. at dph) than the prey (Fig. 3c), and the difference was statistically significant (t test, P < 0.01). This difference indicated that the growth rate of cannibals was faster than that of the prey.

Temperature effect on cannibalism intensity

Cannibalism intensity increased with time at each temperature (Fig. 5), and it corresponded to decreasing larval density in the tanks (Fig. 5) and to increasing fish size (Fig. 3). To account for these two factors—density and fish size—the effect of temperature on CI was analysed by plotting CI against SL (Fig. 6) and against fish density in the tanks (Fig. 7). In both comparisons, cannibalism intensity was higher at 22 °C compared to the 15 °C and 10–15 °C groups. Cannibalism intensity in the 15 °C groups was the same as in the 10–15 °C group. These differences were statistically significant (ANCOVA for log-transformed data, intercept P < 0.05).

Fig. 5
figure 5

Comparison of changes in the timing of cannibalism intensity (Cannibalism Intensity Index) and density of larvae in tanks for pike, Esox lucius, at three different temperatures: 22 °C (square and solid line), 15 °C (circle and dashed line), 10–15 °C (fish transferred from 10 to 15 °C after 34 days of the experiment) (diamonds, dotted line)

Fig. 6
figure 6

Cannibalism intensity (Cannibalism Intensity Index) depending on individual standard length for larval pike, Esox lucius, at three different temperatures: 22 °C (square and solid line), 15 °C (circle and solid line), 10–15 °C (fish transferred from 10 to 15 °C after 34 days of the experiment) (diamonds, no line was fitted)

Fig. 7
figure 7

Cannibalism intensity (Cannibalism Intensity Index) depending on fish density in tanks for larval pike, Esox lucius, at three different temperatures: 22 °C (square and solid line), 15 °C (circle and solid line), 10–15 °C (fish transferred from 10 to 15 °C after 34 days of the experiment) (diamonds, no line was fitted)

Cannibalism and otolith fluctuating asymmetry

The FA of otolith size (area, length, and width of sagitta) was compared between prey (n = 71) and cannibals (n = 71) (Table 1, Fig. 8). Although the FA for all the otolith size variables was higher for prey compared to cannibals, the differences were not statistically significant (paired t test, P > 0.05).

Table 1 Fluctuating asymmetry indexes (FA) (Mean ± SD) of sagittal otolith size calculated for cannibals (n = 71) and prey (n = 71) of larval pike, Esox lucius, kept under experimental conditions for 45 days
Fig. 8
figure 8

Comparison (Mean, SD, and non-outlier range) of fluctuating asymmetry (FA) for otolith area, otolith width, and otolith length for larval pike, Esox lucius, cannibals (n = 71) and prey (n = 71). The differences between cannibals and prey were not statistically significant (t test, P > 0.05)

Discussion

The results of the current study confirmed that cannibalism in pike is a common phenomenon, and it could possibly be a significant factor limiting survival in early developmental stages of the species. However, the cannibalism intensity values obtained in the present research, which were estimated during a few observations daily, were relative only. They were comparable among temperatures, but they did not provide an estimation of the absolute value of total mortality from cannibalism (for example during a 24 h period). This is especially so because total mortality from cannibalism resulted not only from the consumption of whole larvae, but also from death after bites even if the larvae attacked were not consumed (Szczepkowski, 2009). If we recalculate the mean daily mortality of 0.6% obtained in the present research during 90 min/day observations into a 24 h period, a value of approximately 10% daily is obtained. This does not include deaths resulting from unsuccessful attacks. However, we believe that such a calculation is too speculative to include it in the results section.

First cannibalistic attack—cannibal size and age

As the results of the present study show, pike cannibalism occurred as early as in the first days of life. The basic pike feeding pattern described by Hunt & Carbine (1951) can be divided into three stages. The first stage is from the moment pike begin exogenous feeding after yolk-sac resorption until they reach a length of about 25 mm when they feed mainly on lower crustaceans. In the second stage, which is during the period of growth from 26 to 50 mm in length, they feed mainly on immature aquatic insects. In the third stage, when pike are longer than 50 mm, the typical diet is composed almost exclusively of fish and other vertebrates. According to encyclopaedic information (Grabowska & Grabowski, 2013), cannibalism in this typical pattern only appears from a size of approximately 4 cm. Similarly, according to Wright & Giles (1987), pike larger than 30 mm feed exclusively on invertebrates. However, it is plausible that the exact timing (i.e. size and age) of the initiation of cannibalism depends on environmental conditions like temperature and feeding conditions.

During our study conducted under experimental conditions, the smallest cannibals observed were 13.0 mm (22 °C) or 13.4 (15 °C). However, no cannibalistic attacks occurred at 10 °C. These values are lower than those presented by other researchers to date regarding analyses of material obtained from natural conditions and during experiments. The occurrence of early cannibalistic attacks in the natural environment is confirmed, among others, by Žiliukienė & Žiliukas (2006), who report for pike that the smallest cannibalistic individuals were only 16 mm (mean water temperature of 15.7 °C). In the analyses conducted by Hunt & Carbine (1951), the size of the smallest cannibals was 21 mm, while Morrow et al. (1997) report a size of 25 mm. Giles et al. (1986) report that the size of the smallest cannibals under experimental conditions at a temperature of 12 °C was 30.3 mm. Wright & Giles (1987) report that cannibalistic attacks among pike cultured in ponds (no temperature information was provided) occurred at a size of 49 mm. Kucharczyk et al. (1997) observed attacks committed by smaller sized fish of approximately 18–24 mm during experiments conducted at 20 °C.

In reference to the timing, or larval age, of the occurrence of cannibalism, the first cannibalistic attack at the highest temperature analysed (22 °C) was observed as early as 5 dph, while at the temperature of 15 °C, the first attack occurred 10 dph. At the lowest temperature (10 °C) no cannibalistic attacks were noted for 31 days (following which the larvae were transferred to the temperature of 15 °C). Kucharczyk et al. (1997) analysed the phenomenon of cannibalism in young pike stages in experimental conditions at a temperature of 20 °C, and the first cannibalistic attacks, similarly to those in the current study at a similar temperature, occurred at 5–7 dph. However, Giles et al. (1986) noted the first cannibalistic fry (at a size of 30.3 mm) on 32 dph, but since the study was conducted at a lower temperature (12 °C), this could have been a factor limiting the phenomenon of cannibalism in smaller specimens that was similar to that in our study at 10 °C. Other studies report on experiments in which the phenomenon of cannibalism did not occur despite high temperatures (18–24 °C) (Wolnicki & Górny, 1997). It is possible that at some point, too high temperature could cause stress that inhibits cannibalistic behaviour, and the results by Wolnicki & Górny (1997) from 24 °C could confirm this hypothesis. Our results, however, indicate that if such a stressful level of high temperature exists, it occurs above 22 °C.

The timing (i.e. size and age of larvae or juveniles) of the first cannibalistic attacks depends of course not only on temperature but also on larval density and feeding conditions in tanks or the environment. These differed in the studies cited above. However, even if their results are not fully comparable, they provide some estimation of the timing of the first cannibalistic attack in different conditions, including different temperatures. The early onset of cannibalistic behaviour in our study could have resulted from feeding conditions. Even if the larvae were fed ad libitum, the commercial feed might not have been as efficient as live food (zooplankton, larvae of other fish species) in the natural environment. On the other hand, it could happen frequently that the availability of food in natural conditions is insufficient, and the feeding conditions could be similar to ours or worse. Thus, we assumed that the feeding conditions during the experiment were moderately good and were somewhere between very good and poor feeding conditions in nature.

Another factor affecting when cannibalistic behaviour begins is yolk-sac absorption, which can occur at various times depending on the temperature in which the larval pike are held. Fey et al. (2019) report that yolk-sac resorption in larval pike occurred 6 dph at 12 °C, while Demirci et al. (2017) report it occurring 10 dph at a temperature of 10 °C. In the current study, full yolk-sac resorption occurred 4 dph at a temperature of 22 °C, 6 dph at 15 °C, and 9 dph at 10 °C. The timing of yolk-sac resorption corresponded to the appearance of cannibalistic attacks, which were noted from 1 (at 22 °C) to 4 (at 15 °C) days after absorption. At 10 °C, no cannibalistic attacks were observed for the entire observation period, i.e. 25 days after yolk-sac resorption.

Thus, it can be concluded that factors facilitating the occurrence of cannibalism among larval pike was yolk-sac resorption and growth to a size of approximately 13 mm. This state can be obtained within 4 to 6 dph, depending on temperatures from 22 to 15 °C, respectively, which indicates a significant positive metabolism effect on the two factors, i.e. time to yolk-sac absorption and reaching the size of 13 mm. The impact of temperature is not restricted only to limiting the time until yolk-sac absorption and the larvae obtaining a given size. In the present study at a temperature of approximately 10 °C no attacks were noted even after yolk-sac resorption or when the larvae had achieved sizes considerably larger than 13 mm (up to 18 mm). This indicated that temperature affected the physiological state of the fish (Dabrowski, 1986). Clearly, metabolism and larval activity are reduced significantly and maintained at a minimal level below a certain temperature (Elliot, 1982).

Prey to cannibal size ratio

Nillson & Bronmark (2000) analysed the behaviour and morphology of pike and concluded that prey size is limited by its body depth and predator gape size. Body depth can be considered as both height and width since pike can attack from a lateral position. In the present work, the ratio of prey to predator ranged, on average, from 95% for 13 mm fish to 80% for 27 mm fish, with a mean for all specimens of 87%. The smallest size difference between predator and prey was just 0.2 mm and was observed at both 22 °C and 15 °C. Such a small difference was evident at the beginning of the occurrence of this phenomenon (e.g. 5 dph, predator 13 mm—prey 12.8 mm) and during later periods (e.g. 14 dph, predator 18.4 mm—prey 18.2 mm). Žiliukienė & Žiliukas (2006) report that prey length was 90.9% of pike predators measuring 17.6 mm. However, Hunt & Carbine (1951) report this proportion was 70% for pike measuring 23 mm. Bry et al. (1992) demonstrated that for much larger pike held in semi-natural conditions the relative size of the prey in relation to the predator ranged from 35 to 85.4% and fluctuated linearly from approximately 75% at a cannibal total length (TL) of 60 mm to approximately 60% for cannibals measuring 200 mm TL. Clearly, the values of the size proportion between prey and predator decreased (larger differences) with predator growth because of increasing differences in size between the smallest and largest larvae in a given tank over time.

Type of cannibalism—direction of the attack

The present results confirm that the first attacks observed for pike within first few weeks after hatching represent type I cannibalism, in which the prey is only partially ingested from the tail or head before it is completely digested. Cannibalism type II, in which prey are completely ingested, appears with increasing size differences between cannibals and prey (Baras & Jobling, 2002). This type of cannibalism was not observed during the present study.

Cannibalistic attacks among fish during the larval and juvenile periods can target the tail, head, or side of the prey (Baras & Jobling, 2002). In early developmental stages of fish, including pike, attacks are usually directed at the tail, but as size differentiation increases attacks directed at the head or from the side occur more frequently (Baras, 1999; Baras & Jobling, 2002; Colchen et al., 2019). Bry et al. (1992) confirm this with observations that 96% of the prey of pike measuring 60–200 mm SL were eaten from the head. Attacking larval pike from the head is a fairly common strategy under rearing conditions, and especially at high densities (Kucharczyk et al., 1997). Similar tendencies were observed in the present study when almost all of the attacks observed during the early larval stage occurred from the tail, but as the larvae grew, at a size of approximately 17 mm, attacks from side and the head began to appear and were common if the size difference between cannibal and predator was relatively large.

Temperature effect on cannibalism intensity

Undoubtedly, larval or juvenile density is a decisive factor in the intensity of cannibalism (Giles et al., 1986). Data from the literature on larval pike of a size range of 18–24 mm TL cultured for 12 days at densities of 9, 44, and 88 larvae·L−1 indicated there was increased cannibalism intensity at higher densities (Kucharczyk et al., 1997). In the present study, the inverse relationship was observed; cannibalism intensity increased temporally, which corresponded with decreasing larval density from approximately 34.5 to 0.4 n·L−1 (densities at which the cannibals were collected). However, the decrease in density was accompanied by a simultaneous increase in larval length (from approximately 14 mm to 25 mm). Separating these two factors of density and larval size that shifted temporally is impossible, which prevents fully comparing the current results to those obtained by Kucharczyk et al. (1997). Žiliukienė & Žiliukas (2006) report that the phenomenon of cannibalism in pike measuring from 16–22 mm occurred at both high (30 larvae·L−1) and low (2.3 larvae·L−1) densities, but they did not analyse the intensity of this phenomenon. Simultaneously, Žiliukienė & Žiliukas (2006) confirm that increased cannibalism intensity occurred with larval growth from 16.4 to 28.2 mm SL, which was also observed in the present study.

Presumably, the observed temperature effect on cannibalism intensity is a consequence of the well-known phenomenon of the positive impact water temperature has on fish metabolism (e.g. Jobling, 1981). Increased metabolism can, in turn, be responsible for increased feeding intensity (e.g. Johnston & Mathias, 1994) and also possibly on the intensity of cannibalism. However, since this subject was not studied in the present research, any attempt to explain the possible mechanisms would be too speculative.

Cannibalism and otolith fluctuating asymmetry

FA is considered to be an indicator of developmental instability in organisms (Zakharov, 1992), and as such it could impact differences in the stress fish experience and in growth rates. FA is reported to be a useful method to test in various fish species the environmental effects of El Nino conditions (Alados et al.,1993), temperature (Lu & Bernatchez, 1999), pollution (Hardersen, 2000), parasitism (Escós et al., 1995; Loot et al., 2001), salinity (Panfili et al., 2005), and fecundity (Hechter et al., 2000). FA is also useful as an indicator of larval fish condition (Grønkjær & Sand, 2003; Somarakis et al., 1997b) and stress (Gagliano et al., 2008; Lemberget & McCornick, 2009; Castillo-Hidalgo et al., 2018; Landaeta et al., 2018), even if the results obtained by Castillo-Hidalgo et al. (2018) were inconsistent. Most recently, Fey et al. (2020) showed for larval trout the significant effect of a magnetic field of 10 mT on otolith FA of larval trout. The opposite result of larval fish growth rate having no effect on FA is also available (Folkvord et al., 2000; Fey & Hare, 2008; Kristoffersen & Magoulas, 2009; Zenteno et al., 2014). Since negative results exist for growth rates and also other indicators (Diaz-Gil, 2015), the usefulness of FA as a fitness and stress indicator is still questioned. In the present study, prey, as opposed to cannibals, represent a group of larvae within a given cohort characterized first of all by slower growth. However, it cannot be excluded that the prey also experience higher stress from being surrounded by larger aggressive individuals. Indeed, for all the otolith size variables the FA index was higher in the present study for prey compared to cannibals, which indicated higher organism instability in the prey. Although these differences were not statistically significant, they can be treated as an indication of the possible effect of differences between cannibals and prey on otolith FA that could be explored in more detail in the future. However, considering the lack of support in the literature for the effect of growth rate differences on FA, the present differences in FA could be related more to the stress experienced than to growth rate differences between cannibals and prey. There is currently no other work that evaluates the FA issue for cannibal/prey groups for comparison.

Conclusions

Cannibalistic behaviour is a common phenomenon among early life stages of pike, starting as early as 5 dph at a larval size of 13 mm, i.e. as soon as the yolk-sac is absorbed. The phenomenon is affecting individuals with size difference as small as 0.2 mm, which, for larvae from the size range of 13–19 mm, corresponds to approximately 1.3%. The cannibalism could not occur in low temperature and intensify significantly with increased temperatures. No statistically significant differences in otolith FA between cannibals and prey provided no indication of higher organism instability in the prey.