According to theoretical predictions, endothermic homeotherms can be classified as either thermal specialists or thermal generalists. In high cost environments, thermal specialists are supposed to be more prone to using facultative heterothermy than generalists. We tested this hypothesis at the intraspecific level using male laboratory mice (C57BL/cmdb) fasted under different thermal conditions (20 and 10°C) and for different time periods (12–48 h). We predicted that variability of body temperature (Tb) and time spent with Tb below normothermy would increase with the increase of environmental demands (duration of fasting and cold). To verify the above prediction, we measured Tb and energy expenditure of fasted mice. We did not record torpor bouts but we found that variations in Tb and time spent in hypothermia increased with environmental demands. In response to fasting, mice also decreased their energy expenditure. Moreover, animals that showed more precise thermoregulation when fed had more variable Tb when fasted. We postulate that the prediction of the thermoregulatory generalist–specialist trade-off can be applied at the intraspecific level, offering a valid tool for identifying mechanistic explanations of the differences in animal responses to variations in energy supply.
Endothermic thermoregulation requires large amounts of energy to maintain body temperature (Tb) at a high and relatively constant level, which is particularly evident beyond the thermoneutral zone (TNZ) (Scholander et al., 1950). In active endothermic animals, the only way to balance high energy expenditure resulting from the need to maintain Tb is to increase energy intake (Humphries et al., 2003; McNab, 1974; Scholander et al., 1950). However, mammalian Tb varies continuously on a daily and seasonal basis (Geiser and Ruf, 1995; Refinetti and Menaker, 1992), and the degree of its variation differs among taxa (Boyles et al., 2013). According to theoretical predictions (Angilletta et al., 2006, 2010), thermoregulatory strategies of endothermic homeotherms range from specialists, which regulate their Tb precisely, to generalists, which tolerate greater variations of Tb. It was also proposed that in high cost environments, specialists spend more energy to maintain performance than generalists, which perform well in a wide spectrum of environments. Based on the above, Boyles and Warne (2013) further predicted that in the face of energy limitation, specialists should be more prone to using facultative heterothermy than generalists. Food limitation and cold increase the costs of maintaining high Tb and increase the frequency and lengthen duration of torpor bouts (Bozinovic et al., 2007; for rodent examples, see Ruf et al., 1993; Tomlinson et al., 2007). Hence, thermal generalists, or specialists that use facultative heterothermy, should be favored in environments with low or unpredictable energy supply (Lovegrove, 2000). Endotherms deal with energy supply deficit either by increasing activity and searching for food, or by decreasing it, and reducing energy expenditure, i.e. by entering torpor (Gutman et al., 2007). However, increasing activity may lead to resource exhaustion (Russell et al., 1987), and it increases the risk of predation (Lima and Dill, 1990); nevertheless, it brings about an opportunity to find food (Overton and Williams, 2004; Sakurada et al., 2000). Torpor, in turn, is a state of regulated decrease of Tb and metabolic rate (MR) (Heldmaier and Ruf, 1992; Ruf and Geiser, 2015; Snyder and Nestler, 1990), which brings about benefits when food is unavailable or when costs of foraging are too high (Hudson and Scott, 1979; Ruf and Heldmaier, 2000; Schubert et al., 2010; but see Humphries et al., 2003 and Wojciechowski et al., 2011 for a discussion of increased predation risk associated with torpor). In recent decades, several studies have focused on torpor use as a response to energy deficit (Bae et al., 2003; Gutman et al., 2006; Nespolo et al., 2010; Schubert et al., 2010, 2008). Wild house mice (Mus musculus), as well as different laboratory strains, enter daily torpor in the face of low food availability or low ambient temperature (Ta), or both (Morton, 1978; Tomlinson et al., 2007; Webb et al., 1982; Williams et al., 2002), indicating that facultative heterothermy is a response to increased energy demands (Dikic et al., 2008; Gordon, 2012; Hudson and Scott, 1979; Schubert et al., 2010).
To test the hypothesis that mice adjust their thermoregulation in response to environmental demands, we fasted C57BL6/cmdb mice for between 12 and 48 h at different Ta, and measured their Tb continuously. We predicted that variability of Tb and time spent with Tb below normothermy would increase with the increase of environmental demands (i.e. Ta and duration of fasting). As Tb of food-deprived mice showed large variations in response to food deprivation that could not be classified as torpor (see Results), the second aim of this work was to quantify the energetic consequences of such a pattern of heterothermic thermoregulation. To do so, we measured the MR of mice during 48 h food deprivation at constant, moderate Ta. Using data of Tb variability in fed and fasted mice, we tested the hypothesis that facultative heterothermy is more common in thermal specialists than in generalists (Boyles and Warne, 2013); specifically, we examined the prediction that at the intraspecific level, precise regulation of Tb during ad libitum food availability would be negatively correlated with greater variability of Tb during food deprivation.
MATERIALS AND METHODS
All experimental procedures were approved by the Local Committee for Ethics in Animal Research in Bydgoszcz, Poland (decision no. 6/2014 from 16/01/2014, and no. 32/2014 from 27/11/2014). We used two groups of male mice of the C57BL6/cmdb strain originating from the colony maintained at the Centre for Experimental Medicine of the Medical University of Białystok, Poland. In the first series of experiments we used eight individuals, and in the second series we used 12 individuals. Two-month-old mice were transported to the Nicolaus Copernicus University in Toruń and were acclimated for 2 months before starting the experiments. They were kept individually in standard rodent cages (model 1264, Tecniplast, Buguggiate, Italy) in a walk-in climate chamber at Ta=20±2°C and under a 12 h:12 h light:dark cycle, with food (Labofeed B, Kcynia, Poland) and water available ad libitum.
A week before each experiment (experiments 1 and 2; see below), mice were implanted intraperitoneally with miniaturized temperature-sensitive data loggers (modified iButton, model DS1922L, Dallas Semiconductor, Dallas, TX, USA), under ketamine (40 mg kg−1; Ketamine 10%, Biowet, Puławy, Poland) and xylazine (8 mg kg−1; Sedazin 2%, Biowet, Puławy, Poland) anesthesia. Before implantation, loggers were covered with paraffin wax (final mass: 1.5–2.2 g) and calibrated in a water bath against a high-precision thermometer. After implantation, mice were allowed to recover for 3 days and during that time they were supplemented with antibiotic in drinking water (Enflocyne 5%, Biowet, Puławy, Poland). During experiment 1, Tb was logged in 12-min intervals with resolution of 0.5°C, whereas during metabolic measurements (experiment 2), Tb was logged every 5 min with a resolution of 0.065°C.
Experiment 1: response of male mice to increasing environmental demands
In the first series of experiments we aimed to answer the question of whether the thermoregulatory response to food deprivation depends on environmental demands, namely, fasting and cold. Each mouse was deprived of food six times: for 12, 24, 36 and 48 h at 20°C, and then for 12 and 24 h at 10°C. During experiments at 20°C, mice were randomly divided into two groups of four mice each, which differed in the order of food deprivation trails, except for the 48 h fast, which was set as the last fast in both groups. This was done to minimize confounding effects of repeated food deprivation events on the animal’s response to fasting (McCue et al., 2017a). Because the order of trials did not affect Tb [linear mixed-effects (LME) models with order of trials as a fixed factor, body mass (mb) as a covariate and animal ID as a random factor; diurnal Tb: LME: F1,6=0.812, P=0.53, nocturnal Tb: LME: F1,12=1.11, P=0.32], the results from both groups were analyzed together. Thereafter, mice were transferred to 10°C, and after 3 days of habituation were deprived of food for 24 and 12 h. Because there is no information about mice being fasted in the cold for more than 24 h, we did not extend food deprivation for safety reasons. Before each trial, food was withdrawn ∼0.5 h before lights off, and cages, especially bedding material, were carefully checked for uneaten food. Water was available ad libitum during the entire experiment. Mice were weighed before and after each bout of food deprivation to the nearest 0.1 g (SPU402, OHAUS, Parsippany, NJ, USA); mass of the logger was always subtracted from a mouse's mb. After each food deprivation bout, mice were offered food and water ad libitum and were allowed to recover for between 3 and 6 days depending on fast duration, i.e. until their mb returned to the pre-fast level.
where Tb,mod is modal Tb (here, the most frequent Tb recorded when food was available ad libitum), Tbi is Tb recorded at time i, and n is the number of samples. HI during food deprivation (HIfast) was calculated for each trial separately. The HIs of mice fed ad libitum (HIfed) at both Ta were calculated from data obtained during 3-day recordings before the start of food deprivations. To compare HIfed and HIfast, we used LME models with duration of food deprivation included as a fixed factor; this analysis was done separately for 10 and 20°C. To account for repeated observations from the same individuals, in all LME models animal ID was included as a random factor.
Because the thermal state of fasted mice could not be unambiguously described as torpor (see Results; mice were also responsive to the stimuli), we defined it as hypothermia. According to the ‘Glossary of Terms for Thermal Physiology’ (IUPS Thermal Commission, 2001), heterothermy is ‘a pattern of Tb regulation that exceeds in range that characteristic for homeothermy’, while hypothermia is ‘the condition of a temperature regulator when core temperature is below its range specified for the normal active state of the species’, being either regulated or forced (pathological).
The lower limit for normothermic Tb was determined following Wojciechowski and Pinshow (2009) and was calculated separately for each animal. We assumed that normothermic Tb is normally distributed with a center on diurnal Tb,mod. Then we fitted the normal distribution curve shaped by the data equal to or higher than the Tb,mod. Hypothermic Tb (Tb,hypo) was accepted as Tb lower than diurnal (rest-phase) Tb,mod minus two standard deviations. Then we calculated the length of time for which animals were hypothermic during each trial.
Minimum Tb (Tb,min) was calculated separately for day and night as the mean of three consecutive lowest recordings during food deprivation that were taken 12 min apart. Body mass loss (Δmb) resulting from food deprivation was calculated as the difference between initial mb before food deprivation and final mb after food deprivation.
To analyze the effect of increasing environmental demands on HIfast, time with Tb,hypo and Tb,min, and Δmb of fasted mice, we used the LME model in which initial mb and Tb,mod were included as covariates. Duration of food deprivation and Ta were included as fixed factors and animal ID as a random factor. In all analyses, the interaction between duration of food deprivation and Ta was not significant and thus it was excluded from final analyses.
Experiment 2: metabolic rate and body temperature during 48 h food deprivation
To determine the energetic consequences of Tb variations recorded in experiment 1, the MR of food-deprived mice was measured by indirect calorimetry in an open-flow respirometry system. At the same time, we measured Tb of these mice as described above. Measurements were commenced ∼1 h before the active phase and lasted for ∼48 h. During measurements, mice were exposed to Ta=20°C and had access to water ad libitum. Mice were sealed for the duration of the experiment in 0.85 liter chambers constructed of translucent polypropylene containers (HPL 808, Lock&Lock, Hana Cobi, South Korea). To ensure that mice had access to water, ∼150 ml drinking bottle (model ACBT0152, Tecniplast) was mounted in the chamber's lid. Moreover, we put 3 g of sawdust from animal's home cage into each respirometry chamber; this amount of bedding was insufficient to build a nest but may have reduced the stress of a new environment.
We measured respiratory gas exchange of three animals simultaneously, using three parallel respirometry systems. In two systems, rates of O2 consumption (V̇O2) and CO2 production (V̇CO2) were measured, and in the third one only V̇O2. Air was pulled from outside the building using an air pump (5HCE-10-M553, Gast Manufacturing, Benton Harbor, MI, USA) and compressed in a balloon, then dried and scrubbed of CO2 with a PureGas Generator (Puregas, Westminster, CO, USA). The main air stream was split into three chambers and a reference gas stream. The air flow was regulated at ∼500 ml min−1 and measured upstream of animal chambers using a mass flow meter (FlowBar-8, Sable Systems International, North Las Vegas, NV, USA). Gases leaving respirometry chambers were selected sequentially with a computer-controlled multiplexer (Intelligent Multiplexer V3, Sable Systems). Then, air from each gas stream was subsampled at ∼100 ml min−1 and water vapor pressure of the subsampled air was measured with a water vapor analyzer (RH-300, Sable Systems). Air was then dried in a column of magnesium perchlorate (product number 11636.36, VWR International, Gdańsk, Poland) and subsequently the fractional concentrations of CO2 (FCO2) and O2 (FO2) were measured using a FoxBox-C integrated CO2 and O2 analyzer or with separate CO2 (CA-10) and O2 analyzers (FC-10a; both Sable Systems). All electronic elements of the respirometry system were connected to a PC via an analog-to-digital interface (UI2, Sable Systems) and data were acquired using ExpeData software (Sable Systems) at 0.5 Hz. Animal mb was measured before and after MR measurements to the nearest 0.1 g (Scout Pro 200, OHAUS).
Mean normothermic Tb was calculated from recordings obtained during 3 days before fasting, and mean Tb was calculated separately for the first and second day of food deprivation in a metabolic chamber.
For each mouse, its total energy expenditure (EE) during a 46 h food deprivation was calculated by integrating the area below the EE curve (we uniformly subtracted 1 h from the beginning each day during which we paused the recording to exchange desiccating columns and span gas analyzers; during that time, mice remained untouched in respirometry chambers; see below). Energy expenditure for the first day of food deprivation (EE1) was calculated from the beginning of dark phase (17:00 h) until 16:00 h the next day. Energy expenditure for the second day (EE2) was calculated in the same way. Total EE was the sum of EE1 and EE2.
To analyze the relationship between duration of food deprivation and Tb, we used an LME model with initial mb as a covariate, day of food deprivation as a fixed factor and animal ID as a random factor; in the analysis of the effect of food deprivation on EE, we included mean Tb and HI during each particular day as covariates.
Finally, to determine whether individuals that regulated their normothermic Tb more precisely showed greater variations of Tb during fast, we used an LME model with HIfast as a dependent variable, HIfed (here calculated from 2-day recordings before the experiment and used for analysis of the first and second day of food deprivation) and mb as covariates, day of food deprivation as a fixed factor and animal ID as a random factor. Because mb was not a significant covariate (F1,12=0.00, P=1.00), it was removed from the final model. We performed two analyses: first, for data obtained only during the active phase of the day (i.e. night), as suggested by Boyles and Warne (2013); and second, for data obtained for the whole 24 h, because fasted mice showed large variations in Tb both during active and inactive phases of the day.
All analyses were performed in SPSS Statistics v. 23 (IBM, Armonk, NY, USA). In all LME models, the restricted maximum likelihood method was used to estimate variance components. Assumptions of the linear modeling were checked post hoc by inspecting the distribution of residuals obtained from LME (check of histograms and quantile-quantile plots; Grafen and Hails, 2002). Data are presented as estimated marginal means±s.e.m., and statistical significance was accepted at α<0.05.
Experiment 1: response of mice to increasing environmental demands
Overall, when food deprived for 24 h, mice maintained Tb lower by ∼0.5°C than when they were fed ad libitum (LME with feeding status and Ta as fixed factors and animal ID as a random factor: F1,20=24.42, P<0.001). However, Tb was not maintained at the low level continuously, but after each drop it returned almost immediately to normothermy (Fig. 1). Body temperature was more variable during food deprivation than when mice were fed ad libitum, at both 20°C (F4,28=13.91, P<0.001) and 10°C (F2,28=20.13, P<0.001). HI during food deprivation also increased with the decrease in Ta (F1,28=7.68, P=0.01), and during 12 h fast at 20°C it equaled 1.56±0.17°C, whereas at 10°C it was 2.27±0.20°C (Table 1). Lengthening of food deprivation also resulted in higher HIfast (F3,28=3.36, P=0.03), and during 24 h food deprivation at 20°C HIfast was 1.63±0.17°C, while during 48 h it increased to 2.34±0.20°C (Table 1). The normothermic Tb of mice was 36.00±0.45°C and the lower limit for normothermy was 34.20±0.64°C at Ta=20°C and 33.24±0.55°C at Ta=10°C. Time spent with Tb,hypo increased with the duration of food deprivation (F3,28=12.13, P<0.001; Table 1). During the 12 h fast at Ta=20°C, mice spent 34.45±114.39 min with Tb,hypo, and during the 48 h fast they spent 1072.29±135.43 min with Tb,hypo (Table 1). Ambient temperature did not significantly affect time with Tb,hypo (F1,28=1.31, P=0.26; Table 1). Diurnal Tb,min decreased with Ta (F1,28=14.31, P=0.01) and with duration of food deprivation (F3,28=22.07, P<0.001; Table 1). Nocturnal Tb,min did not change with Ta (F1,28=2.32, P=0.14) but decreased as food deprivation lengthened (F3,28=10.68, P<0.001; Table 1). Mice lost more body mass (Δmb) as food deprivation was lengthened (F3,28=13.25, P<0.001) and Ta was decreased (F1,28=23.46, P<0.001; Table 1). Δmb did not correlate with initial mb (F1,28=0.01, P=0.94) or with Tb,mod (F1,28=0.07, P=0.79).
Experiment 2: metabolic rate and body temperature during 48 h food deprivation
Food deprivation led to a decrease of MR and Tb (Fig. 2). Mean normothermic Tb before food deprivation was 36.36±0.12°C. There was a significant effect of fasting duration on mean Tb (F2,8=130.15, P<0.001). During the first day, mean Tb decreased to 33.95±0.42°C, and during the second day to 31.58±0.99°C. Total EE during 46 h of food deprivation was 66.14±2.88 kJ. Energy expenditure was not significantly different between the two days of food deprivation (day 1: 34.74±155 kJ; day 2: 31.40±1.55 kJ; F1,11=1.41, P=0.26). Energy expenditure did not depend on initial mb (F1,11=0.22, P=0.64) or on HI (F1,11=2.32, P=0.15). Differences in energy expenditure during food deprivation were explained only by mean Tb (F1,11=5.51, P=0.04). The HI of food-deprived mice during their active phase did not correlate with active-phase HI of mice fed ad libitum (F1,13=2.43, P=0.14; Fig. 3A). However, when calculated for the entire 24 h, HIfast correlated negatively with HIfed (F1,13=5.41, P=0.04; Fig. 3B). It was true for both days of food deprivation, and HIfast on the second day of fast was higher than that on the first day (F1,13=27.04, P<0.001).
Any Ta lower than thermoneutral is an environmental stress for animals (reviewed in Kingma et al., 2012; Ravussin et al., 2012). We found that laboratory mice (strain C57BL6/cmdb) adjusted their thermoregulation in response to increased environmental demands. Their Tb was more variable and they spent more time below normothermy when they were deprived of food, and this effect was augmented by cold. Our findings give strong support to Angilletta and co-authors’ (2010) model of thermal physiology of endotherms, which posits that variations in Tb increase when food is limited and also when operative temperature (here, Ta) decreases. Moreover, we found that mice that showed more precise thermoregulation when fed had more variable Tb when fasted, supporting the prediction of the thermoregulatory generalist–specialist trade-off (Angilletta et al., 2010; Boyles and Warne, 2013) at the intraspecific level.
Thermoregulation depends on several factors such as Ta (Geiser, 2004; Overton and Williams, 2004; Williams et al., 2002), availability of bedding material and social interactions (Gordon, 2004), but also on sex (Geiser and Mzilikazi, 2011; Lovegrove and Raman, 1998). Variations of Tb that were observed in fasted male mice (Figs 1 and 2) could not be defined as torpor. Nevertheless, the use of heterothermy correlated positively with the duration of food deprivation. Such a correlation agrees with thermoregulatory adjustments observed in response to prolonged fasting or food restriction in rats (Wang et al., 2006; Yoda et al., 2000; McCue et al., 2017a), laboratory mice (McCue et al., 2017b), spiny mice (Acomys russatus; Gutman et al., 2006) and Chilean mouse-opossums (Thylamys elegans; Bozinovic et al., 2007). The lack of prolonged torpor in fasted male mice might be related to sex. In many species, females have been shown to be more prone to using torpor, whereas males often showed only slight variations of Tb (Geiser and Mzilikazi, 2011; Lovegrove and Raman, 1998). Swoap and Gutilla (2009) suggested that it was easier to induce torpor in female than in male C57BL/6 mice. This was supported by the present results; we used male mice and we did not find torpor bouts described for females of the species elsewhere (Dikic et al., 2008; Hudson and Scott, 1979; Schubert et al., 2010; Swoap and Gutilla, 2009). However, contrary to animals that use torpor spontaneously and lower their Tb only during rest phase (Ruf et al., 1993), fasted animals enter torpor regardless of the time of day (Dikic et al., 2008; Swoap and Gutilla, 2009), which is in line with the present results. Tb of mice started to decrease on the first day of the experiment, and both diurnal and nocturnal Tb,min decreased as food deprivation prolonged (Figs 1 and 2). Note, however, that Kanizsai and co-authors (2009), who fasted mice at Ta=28°C, found that the mice lowered normothermic Tb only at night, and only during the second day of fasting.
The present results show that Tb decreased with the duration of fasting (Table 1), which led to a decrease in daily energy expenditure. Hence, it seems that although mice did not enter deep torpor, they clearly could benefit from shallow and interrupted hypothermia in the face of increased environmental demands. Thus, even low variability of Tb and a concomitant decrease in MR seems to be an effective way to save energy in the face of a decrease in energy supply. Cold resulted in increased HI and Δmb, and in lower Tb,min, yet we did not find a significant effect of Ta on time spent with Tb below normothermy (Table 1). Williams and co-authors (2002) reported that food-restricted C57BL/6J mice reduced their energy expenditure much more at 23°C than at 30°C (TNZ), which could support the idea that variability of Tb increases with decreasing Ta (Gordon, 2009; but see Ravussin et al., 2012). Also, Overton and Williams (2004) showed that response to caloric restriction leads to greater changes in physiology at Ta values below TNZ, and in other studies, the effects of limited energy supply and decreased Ta on heterothermy use were found to be additive (Gordon, 2012; Nespolo et al., 2010; but see Ravussin et al., 2012; Williams et al., 2002).
Fasting or caloric restriction is a challenge, especially for small endothermic animals, and it significantly affects their physiology and behavior (Overton and Williams, 2004; Sakurada et al., 2000; Tucci et al., 2006). During fasting or prolonged food restriction in mice, heart rate, blood pressure, oxygen consumption and pulse are reduced gradually with increasing duration of negative energy balance (Swoap and Gutilla, 2009; Williams et al., 2002). Fasting may also lessen exploratory behavior and result in memory disorders (Tucci et al., 2006). However, it is known that regulating Tb at low levels during fasting provides energy savings (Hudson and Scott, 1979; Schubert et al., 2010), which may enhance survival in the face of adverse conditions (Geiser and Turbill, 2009), and lead to improved fitness of heterothermic animals compared with homeothermic ones (Angilletta et al., 2010). In line with that, we found a significant correlation between EE and mean Tb of fasted mice.
According to models of Gilchrist (1995) and Angilletta et al. (2010), animals can be classified as specialists or generalists, and they should differ in thermoregulatory responses to increasing environmental demands, such as food deprivation or cold (Angilletta et al., 2010; Boyles and Warne, 2013). Specialists should thermoregulate more precisely when energy supply is not limited, but in the face of high energy demands, their Tb should be more variable, leading to lower costs of maintaining homeothermy (Boyles and Warne, 2013). In support of this prediction, male mice that were characterized by lower HIfed showed higher HI during both the first and second days of fasting (Fig. 3B). This negative correlation between variability of Tb under feeding and fasting conditions is the first experimental support for a specialist–generalist trade-off at the intraspecific level (Boyles and Warne, 2013). However, the observed intra-individual variability of thermoregulatory responses to fasting (Fig. 3B) may indicate the proposed specialist–generalist trade-off is a continuum that includes a wide spectrum of thermal sensitivities, similar to what was inferred from interspecific analyses (Boyles et al., 2013). Moreover, intraspecific variability in the use of heterothermy might be important for predicting animal responses to changing conditions on the biogeographic scale (Bozinovic et al., 2011). Heterothermy could mitigate the effects of abiotic conditions on species distribution, enhancing their tolerance to increasing environmental demands (Boyles et al., 2011a). The presence of the proposed specialist–generalist trade-off at the intraspecific level also suggests that heterothermy could be correlated with selection for high mass-independent metabolic rates (Bozinovic et al., 2011; Boyles and Warne, 2013). Our results indicate that hypotheses testing trade-offs in the evolution of thermoregulatory strategies in endotherms can also be tested at the intraspecific level, offering a valid tool to seek mechanistic explanations of the observed differences in animal responses to variations in energy supply and environmental demands.
We thank Dr Justin G. Boyles for critical reading and comments on the manuscript. We also thank two anonymous reviewers for their constructive comments.
Conceptualization: A.S.P., J.B., M.S.W., M.J.; Methodology: A.S.P., J.B., M.S.W., M.J.; Formal analysis: A.S.P., M.S.W.; Investigation: A.S.P., M.S.W., M.J.; Writing - original draft: A.S.P., J.B., M.S.W., M.J.; Visualization: A.S.P.; Funding acquisition: J.B., M.J.
The study was partially funded by the Faculty of Biology and Environmental Protection of Nicolaus Copernicus University funds to M.J. and partially by a grant from the Faculty of Biology and Environmental Protection of Nicolaus Copernicus University awarded to J.S.B. (grant no. 1915-B).
The authors declare no competing or financial interests.