E ﬀ ect of Temperature on the Size of Sedimentary Remains of Littoral Chydorids

: The body size of aquatic invertebrates is, to a great extent, dependent on ambient temperature, but size distributions are also determined by other factors like food supply and predation. The e ﬀ ect of temperature on organisms is formulated in the temperature–size hypothesis, which predicts a smaller body size with increasing temperature. In this study, the e ﬀ ect of temperature on the subfossil remains of three littoral Cladocera ( Alona a ﬀ nis , A. quadrangularis , and Chydorus cf. sphaericus ) was investigated. Exoskeletal remains of these species can be found in large numbers in lacustrine sediments and over a wide north–south range in Europe. The total length of both headshield and postabdomen for A. a ﬃ nis and A. quadrangularis and carapace length for C. cf. sphaericus were measured to observe their response to changes in latitude and temperature. A di ﬀ erent response to ambient temperature in the growth of body parts was observed. The size of the headshields of both Alona species and of the carapace of Chydorus was signiﬁcantly larger in colder regions as opposed to warm ones. It turned out that the postabdomen was not a good predictor of ambient temperature. While the sizes of all remains increased with latitude, the sizes of the Alona remains was smaller in the mountain lakes of the Southern Carpathians than in other cold lakes, in this case in Finland, a fact indicative of the importance of other factors on size distribution. This study demonstrates that a morphological response to climate is present in littoral cladocerans, and, therefore, changes in the length of headshield and carapace may be used as a proxy for climate changes in paleolimnological records.


Introduction
Crustacean cladocerans (Branchiopoda) are functionally adapted to diverse microhabitats in lakes. Some species live in open water (pelagic species, according to the strict definition of zooplankton), whereas others are attached to macrophytes (epiphytic) or abiotic surfaces, such as lake beds covered with silt, pebbles or sand. Cladocerans occupy a central position in the food web; their abundance and assemblage depend on bottom-up (food supply) and top-down (predation pressure) forcing to the same degree as on habitat availability. Therefore, any changes in cladoceran community composition and functional traits may well indicate changes in the prevailing environmental conditions [1]. The body size of cladocerans is a functional trait of great importance, because it has an influence on resource use, and may indicate predation patterns. This is because large-bodied zooplankton have greater grazing rates, they filter a wide range of food particles, and hence exert a stronger top-down effect on resources than small-bodied zooplankton do. It should also be noted that a larger body size tends to make them vulnerable to vertebrate predators, while a smaller size is preferred by invertebrate predators [2].
Both seasonal and latitudinal patterns can be observed in the body size distribution of zooplankton [3][4][5]. In warm water, the body size of zooplankton is generally smaller than that of those in cold water, and consequently zooplankters are larger in various waters in the cold season or in lakes at higher latitudes and/or altitudes, e.g., [5][6][7][8][9]. The direct effect of temperature on body size can be masked by other environmental factors that are related directly or indirectly to temperature, not to mention the fact that predation also has a significant influence on body size spectra [10]. In earlier studies where the relationship between planktonic cladoceran body size and latitude and temperature has been examined [5,11], the results have shown that the body size of planktonic cladocerans declines with temperature and increases with latitude. This decrease in body size in relation to temperature has also been found in littoral cladoceran species [12] and in the majority of ectotherm organisms [13,14].
Paleolimnological studies focus on the remains of organisms (e.g., algal pigments, chrysophyte cysts, diatoms, aquatic pollen, invertebrates and fish) and their use as proxies for paleoenvironmental reconstructions [15]. Due to taphonomic processes and the selective preservation of organisms and their remains, a certain degree of bias can be observed reconstructed communities are compared to modern ones [16][17][18][19]. Cladocerans are one of the most important paleolimnological proxies for discovering the trajectory of environmental changes [20]. Besides the species composition of cladoceran communities, their size structure also responds to changes in environmental conditions [21][22][23][24][25][26][27]. To date, comprehensive limnological studies of cladoceran body size distribution have concentrated on planktonic species [5,11], while studies of littoral habitats have revealed the horizontal migration of zooplankton between pelagial and littoral environments [28]. Although macrophytes provide shelter against pelagic predators, predation risk may nonetheless remain high due to the high abundance of invertebrate predators, and fish, which use the littoral bed as a way of avoiding predators [29]. Moreover, the density of plants also affects the predation capability of fish and invertebrate predators [30,31]. Despite the fact that the denser the macrophyte coverage, the better the refuge they offer, juvenile fish and invertebrates tend to suppress zooplankton communities in littoral zones/environments [32][33][34]. Furthermore, the structure of littoral zones affects the species distribution of littoral cladocerans, since cladocerans effectively share resources, allowing the co-existence of a high number of species [35][36][37]. Analyzing fish gut content, Adamchuk [36] demonstrated the high degree of selectivity in various fish species, and that fish generally prefer large-bodied littoral cladocerans. Basińska et al., [38] found the body size of littoral cladocerans living among macrophytes to be larger than in that of those in open water, and ascribed this to lower predation pressure. Predation on cladocerans in the littoral zone follows complex patterns. Predation by fish is more important in environments with low turbidity, since they may be regarded as generally visual predators, while invertebrates, which detect their prey by visual and tactile means, become an important factor in structuring littoral cladoceran communities in turbid environments [10,31]. The thickness of the carapace can also be seen as a defensive morphological trait against invertebrate predation [39][40][41].
Rizo et al. [42] investigated patterns of distribution and variability of body size among freshwater Cladocera taxa across different taxonomic levels, geographic distributions and habitat associations. They found that cladoceran body size is a phylogenetically conserved trait, but the observed variation and the varying means between biogeographical regions also suggests the effect of climate on size variability in cladocerans. This factor explains the high degree of the variance observed in size spectra in particular regions as a product of habitat variability, resource availability and the effect of predators. The same research also found that substrate-associated littoral species were smaller than to planktonic ones. They concluded that the smaller body is an adaptation for survival in the presence of high numbers of predators and low oxygen concentrations in preferred habitats.
In this study, the relationships between the size of sedimentary chydorid remains preserved in surface sediments and temperature at a large geographical scale from south to north of Europe are examined. Based on both correlation of the size of body parts with environmental conditions [23] and Atkinson's rule [13,14], the hypothesis that the mean size of chydorid remains will be larger in cold environments is adopted. In this hypothesis, there is an a priori assumption that benthic chydorids are less sensitive to top-down effects than planktonic species, and that this, in turn, highlights the temperature effect on the body size distribution. The focus is on the dominant chydorid species (Alona affinis (Leydig), Alona quadrangularis (O.F. Müller) and Chydorus cf. sphaericus), because the occurrence of their remains in sediments was frequent enough throughout the dataset to make statistically adequate measurements possible.

Study Sites
For the analysis, lakes in Finland, Hungary, Romania, Croatia and Spain were selected. The distribution of these lakes covers wide latitudinal and climatic gradients across Europe ( Figure 1; Table 1), from 40 • N in Spain to 70 • N in Finland. This latitudinal gradient also reflects a temperature gradient: the yearly average temperature in the most southern lake is 10.8 • C [19], while in the northernmost lake it is −1.9 • C [43]. The surface sediments of the lakes were sampled by means of sediment cores taken from the deepest part of the lakes, using gravity, piston and Russian peat corers. The uppermost portion of the sediment (0-1 cm, or in the case of very soft sediments, 0-2 cm) was used for the analysis. The exact sedimentation rate for the lakes was not known, but it was assumed that the uppermost centimeter represented a decadal time span, with the exception of Lake El Tobar (Spain). A short core was assessed from El Tobar, since it has very high accumulation rates, and the core represented several recent decades [19]. Additionally, sediment traps were also used at Lake Balaton [44]. In Finland, the sampling was conducted between 2003 and 2013 in winter, while at the other locations, it was performed between 2010 and 2014 in summer. Lake Balaton [44]. In Finland, the sampling was conducted between 2003 and 2013 in winter, while at the other locations, it was performed between 2010 and 2014 in summer.
Mean July air temperature data were used as a proxy for the water temperature of the lakes, since mean air temperature and epilimnetic water temperature correspond most closely in July [45,46]. Mean July temperatures were extracted from the WorldClim database at a resolution of 30 arc seconds (1 km 2 ) using dismo R packages [47][48][49].    Mean July air temperature data were used as a proxy for the water temperature of the lakes, since mean air temperature and epilimnetic water temperature correspond most closely in July [45,46]. Mean July temperatures were extracted from the WorldClim database at a resolution of 30 arc seconds (1 km 2 ) using dismo R packages [47][48][49].

Subfossil Analyses
For the analysis of subfossil cladoceran remains, between 2 and 8 g of sediment were subsampled and treated with 250 ml of 10% KOH solution, heated to 80 • C and kept at that temperature for at least 1 h, following the standard protocol for fossil Cladocera analysis [50]. Subsamples were gently stirred Water 2020, 12, 1309 5 of 21 with a plastic rod to deflocculate all organic matter without damaging the remains, and sieved through a 65-µm mesh with running tap water. Quantitative slides were prepared by pipetting 100-200 µL of each subsample on a microscope slide, and these were mounted using glycerol jelly stained with safranin. At least 20-30 headshields and postabdomens of Alona affinis and A. quadrangularis, and c. 30 carapaces of Chydorus cf. sphaericus were measured under a microscope at 100× and 200× magnifications to an accuracy of 1 µm using microphotography. Four morphological characteristics of the headshield (Figure 2a (headshield: 1. total length, 2. width, 3. length of distal part, 4. distance from the first pore to the tip of headshield)) and three of postabdomen (Figure 2b (postabdomen: 1. length, 2. length of ventral margin, 3. width of postabdomen)) were measured for the Alona species and the length of the carapace was measured for Chydorus cf. sphaericus (Figure 2c).
Water 2020, 12, x FOR PEER REVIEW 5 of 22

Subfossil Analyses
For the analysis of subfossil cladoceran remains, between 2 and 8 g of sediment were subsampled and treated with 250 ml of 10% KOH solution, heated to 80 °C and kept at that temperature for at least 1 h, following the standard protocol for fossil Cladocera analysis [50]. Subsamples were gently stirred with a plastic rod to deflocculate all organic matter without damaging the remains, and sieved through a 65-µm mesh with running tap water. Quantitative slides were prepared by pipetting 100-200 µL of each subsample on a microscope slide, and these were mounted using glycerol jelly stained with safranin. At least 20-30 headshields and postabdomens of Alona affinis and A. quadrangularis, and c. 30 carapaces of Chydorus cf. sphaericus were measured under a microscope at 100× and 200× magnifications to an accuracy of 1 µm using microphotography. Four morphological characteristics of the headshield (Figure 2a

Data Analyses
In the case of the fragmented remains of Alona species, linear regression was used to determine the relationship between the length of the headshield and other morphological characters (width, distal length and pore distance), and to establish the relationship between the length of the postabdomen and other morphological characters (width and ventral length). Due to the high proportion of fragmented remains, only those fragments were incorporated into the analysis upon which it was possible to perform at least one measurement.
Prior to linear regression analysis, outliers were determined and omitted from the regressions. Outliers included those with a 2 SD range of residuals of regression (in all, outliers constituted about 10% of total data). The linear model with the highest R 2 values was selected for the calculation of the inferred length, and this inferred length was then used as a substitute for missing data. The distributions of the headshield and postabdomen data were estimated by fitting normal, lognormal and gamma distributions, and from these, the best were selected according to the AIC score (Akaike's An Information Criterion), computed using the maximum-likelihood estimation method and the Kolmogorov-Smirnov goodness of fit test. The analyses were conducted using the 'fitdistrplus' R package.
Given that the sampling design was unbalanced (cold lakes outnumbered warm lakes), the following strategy was pursued in the grouping of the lakes. First, the question of whether the length of remains from cold environments differed from those from warm environments. At this point, lakes were grouped into "warm" and "cold" sets according to their yearly July mean air temperature. The Finnish and mountain lakes were considered "cold" lakes (range of Tmean July: 7-17 °C), and the others categorized as "warm" lakes (range Tmean July: 18-24 °C). Next, the cold lakes were examined to see if

Data Analyses
In the case of the fragmented remains of Alona species, linear regression was used to determine the relationship between the length of the headshield and other morphological characters (width, distal length and pore distance), and to establish the relationship between the length of the postabdomen and other morphological characters (width and ventral length). Due to the high proportion of fragmented remains, only those fragments were incorporated into the analysis upon which it was possible to perform at least one measurement.
Prior to linear regression analysis, outliers were determined and omitted from the regressions. Outliers included those with a 2 SD range of residuals of regression (in all, outliers constituted about 10% of total data). The linear model with the highest R 2 values was selected for the calculation of the inferred length, and this inferred length was then used as a substitute for missing data. The distributions of the headshield and postabdomen data were estimated by fitting normal, lognormal and gamma distributions, and from these, the best were selected according to the AIC score (Akaike's An Information Criterion), computed using the maximum-likelihood estimation method and the Kolmogorov-Smirnov goodness of fit test. The analyses were conducted using the 'fitdistrplus' R package.
Given that the sampling design was unbalanced (cold lakes outnumbered warm lakes), the following strategy was pursued in the grouping of the lakes. First, the question of whether the length of remains from cold environments differed from those from warm environments. At this point, lakes were grouped into "warm" and "cold" sets according to their yearly July mean air temperature. The Finnish and mountain lakes were considered "cold" lakes (range of T mean July : 7-17 • C), and the others categorized as "warm" lakes (range T mean July : 18-24 • C). Next, the cold lakes were examined to see if there comprised a homogeneous group, or whether there are any differences between latitudinal regions, and thus, "cold" lakes were further divided into regions. The Southern Carpathian formed a "Southern cold lakes" group (range of T mean July : 8-11 • C), and the Finnish lakes fell into a group named "Nordic cold lakes". Since the Finnish lakes distributed over a wide range of latitude (ca. 10 • N), they were further divided into subregions, namely, Southern Finnish (17 lakes, between 60-63 • N latitude), Central Finnish (37 lakes, between 63-66 • N) and Northern Finnish (24 lakes, above 66 • N).
To investigate the effect of depth on the size distribution of remains, the lengths of remains in deep and shallow lakes in the Southern Carpathians were compared. Two lake groups could be identified on the basis of the depth of the lakes. Since thermocline does not usually penetrate below 10 m, on this basis, three lakes could be assigned to the "deep lakes" group, while the others (n = 7) became the "shallow lakes" group.
Non-parametric Kruskal-Wallis and Wilcoxon tests were used for between-group comparisons, and the posthoc Nemenyi-test (package PCMCRplus [51]) was employed in the pairwise comparison of lake groups. Finally, a generalized linear model (GLM) was used to identify the effects of latitude and temperature on the length of the remains. All statistical analyses were performed in R [52].

Results
The various morphological characteristics examined turned out to provide an accurate basis for the estimation of the length of remains (Table A1). For the Alona species, the distal length of headshields and width of the postabdomen best corresponded to total length, and therefore these were used to infer total length for fragmentary remains ( Figures A1 and A2).
The lengths of the headshield and postabdomen of A. affinis were significantly greater in cold regions than in warm ones (Kruskall-Wallis test: χ 2 = 43.254, df = 1, p-value < 0.05), but the remains of A. quadrangularis displayed a slightly different pattern ( Figure 3). The average length of the headshield of A. quadrangularis was also significantly greater in cold regions (Kruskall-Wallis test: χ 2 = 12.134, df = 1, p-value = 0.05), but its postabdomen length did not differ between the climate regions ( Figure 3). The mean length of the carapace of C. cf. sphaericus was also greater in cold regions (Kruskall-Wallis test: χ 2 = 7.5481, df = 1, p-value < 0.05; Figure 3). there comprised a homogeneous group, or whether there are any differences between latitudinal regions, and thus, "cold" lakes were further divided into regions. The Southern Carpathian formed a "Southern cold lakes" group (range of Tmean July: 8-11 °C), and the Finnish lakes fell into a group named "Nordic cold lakes". Since the Finnish lakes distributed over a wide range of latitude (ca. 10°N), they were further divided into subregions, namely, Southern Finnish (17 lakes, between 60-63°N latitude), Central Finnish (37 lakes, between 63-66°N) and Northern Finnish (24 lakes, above 66°N).
To investigate the effect of depth on the size distribution of remains, the lengths of remains in deep and shallow lakes in the Southern Carpathians were compared. Two lake groups could be identified on the basis of the depth of the lakes. Since thermocline does not usually penetrate below 10 m, on this basis, three lakes could be assigned to the "deep lakes" group, while the others (n = 7) became the "shallow lakes" group.
Non-parametric Kruskal-Wallis and Wilcoxon tests were used for between-group comparisons, and the posthoc Nemenyi-test (package PCMCRplus [51]) was employed in the pairwise comparison of lake groups. Finally, a generalized linear model (GLM) was used to identify the effects of latitude and temperature on the length of the remains. All statistical analyses were performed in R [52].

Results
The various morphological characteristics examined turned out to provide an accurate basis for the estimation of the length of remains (Table A1). For the Alona species, the distal length of headshields and width of the postabdomen best corresponded to total length, and therefore these were used to infer total length for fragmentary remains ( Figure A1 and A2).
The lengths of the headshield and postabdomen of A. affinis were significantly greater in cold regions than in warm ones (Kruskall-Wallis test: 2 = 43.254, df = 1, p-value < 0.05), but the remains of A. quadrangularis displayed a slightly different pattern ( Figure 3). The average length of the headshield of A. quadrangularis was also significantly greater in cold regions (Kruskall-Wallis test: 2 = 12.134, df = 1, p-value = 0.05), but its postabdomen length did not differ between the climate regions ( Figure 3). The mean length of the carapace of C. cf. sphaericus was also greater in cold regions (Kruskall-Wallis test: 2 = 7.5481, df = 1, p-value < 0.05; Figure 3).  A. quadrangularis in the Nordic populations were not distinguishable from those in the warmer regions, and the smallest postambdomens were found in Carpathian lakes (Kruskal-Wallis test: χ 2 = 73.398, df = 2, p-value < 0.05). The largest Chydorus sphaericus carapaces were found in the Carpathian mountain lakes (Kruskall-Wallis test: χ 2 = 26.826, df = 2, p-value < 0.05).
The next task was to check how the climatic gradient of Finland affected size distribution of the remains, and therefore the Finnish lakes were grouped into three subregions. The size of remains from the Finnish subregions was very similar ( Figure A4), with the exception of Chydorus carapaces, which were smallest in the Southern Finnish lakes ( Figure A4). The size distribution of sedimental remains between subregions reflects a latitudinal effect, and an increasing trend can be seen in the mean size of remains towards higher latitudes ( Figure A4).
The size-frequency distribution of the remains of all species followed lognormal and normal distributions (Appendix A: Figure A2), and the lognormal distribution also applied to the headshield of A. affinis, and postabdomen of A. quadrangularis, the headshield of A. quadrangularis and carapace of C. sphaericus displayed a normal distribution These distributions were selected on the basis of the smallest AIC scores (Appendix A: Table A2). Latitude-size (L-S) GLM models revealed an increasing trend in the size of all remains, while temperature-size models (T-S) showed changes in the opposite direction in the length of Alona remains (Figures 4-6, Appendix A: Table A3). The body length of Chydorus sphaericus decreased with increasing temperature (Figure 6, Appendix A: Table A3). Since remains were generally smaller in lakes of the Southern Carpathians than in Finland in spite of the latter's cooler environment, the issue of how the performance of GLM models might change if Southern mountain lakes were excluded from the models was put to the test (Table A3). The L-S models of Alona remains displayed a significant positive relationship with latitude, and this did not change significantly when data from the Carpathian lakes were left out of the models. In contrast, the T-S model of length of all Alona remains showed a positive relationship with temperature when all Carpathian lakes were omitted from the model. The L-S model of C. spharicus did not correspond to temperature, but did display increasingly better performance when the Carpathian lakes were omitted.
The T-S model of carapace data also worked well, displaying an inverse correlation with temperature.   Table A3).        Table A3).   Table A3).   Table A3).

Southern mountain lakes,
Water 2020, 12, x FOR PEER REVIEW 8 of 22

Discussion
This is the first study in which the morphological response of littoral cladocerans has been examined in relation to the ambient temperature of the location of their sedimentary remains. Typically, paleolimnological inferences are based on assemblages of cladoceran species, but cladoceran size structure is also sensitive to changes in limnological conditions [53]. Mean cladoceran body size, a community metric, has been used mainly to infer fish and invertebrate predation regimes [54][55][56], but relationships between body size and climate have also been documented [5,7,57,58].
The comparison of inferred length to measured length reveals that the allometric growth of body parts was not to be observed, and inferred length could, therefore, be used to supplement missing data ( Figures A1 and A2).
The analyses revealed that the average length of headshield and carapace is smaller in cooler environments, as was expected. This result confirms the temperature-size hypothesis (TSR: [13,59]), and hence, the remains are generally smaller in warm lakes. However, the unexpected result of the analysis was that the length postabdomen of A. quadrangularis in warm lakes did not differ significantly from that of those found in cool environments. This can be explained by the distinctly different shapes of the postabdomen of A. affinis and A. angularis. The causes of these differences are difficult to unravel, but they may be related to habitat preferences and the ecological niche of the Alona species. The large Alona sp. (A. quadrangularis and A. affinis) are adapted to benthic habitats, and A. quadrangularis is specifically considered a mud dweller [60][61][62], while A. affinis prefers dense vegetation [61], but it can be found widely over the littoral zone [36]. Littoral cladocerans have developed many systems to exploit periphyton, including mechanisms of food selection, along with many other particular adaptations [63,64]. The postabdomenis is a key morphological feature of these cladocerans, as in the case of A. quadrangularis, which uses these parts for locomotion and to clean the food canal and toracopodes [63]. W might propose that in this, A. quadrangularis displays an evolutionary adaptation to a muddy bottom environment, with its widened carapace posteriorly. It is also possible that allometric body growth, in which the larger postabdomen has an advantage, is an evolutionary strategy for the maximization of niche exploitation. The axe shape (rounded ventral margin) of postabdomen of A. quadrangularis is very like that of other mud dwelling species (Ilyocryptus spp., Leydigia spp.), and thus can be regarded as an adaptation to a muddy environment. Low-oxygen conditions are also a factor to consider when dealing with benthic organisms. According to Van-Damme et al. [65] benthic Aloninae require a large "exopodite pump", which is probably a very specialized feature, along with the postabdomen. These factors together may well explain the differences found in the response of postabdomen length to ambient temperature in comparison with that observed in the case of headshields.
It was found that, in general, larger remains predominated in colder environments, confirming the hypothesized effect of climate on the size of the the distribution of littoral species, while the effect of climate on the size of planktonic species has previously been documented. Manca and Comoli [58] found larger Daphnia headshields coinciding with cold climate events in the late Holocene in the Nepalese Himalayas. In the regions in this study, Korponai et al. [25] found the occurrence of larger ephippia in Daphnia during the Late Glacial than those in the Holocene in the South Carpathians, and Nevalainen et al. [66] also documented that the size of the ephippia of Bosmina and Daphnia were larger during the cold periods of the Late Holocene in a Finnish lake. The physiological advantages of larger cladoceran body size in cold climatic conditions are manyfold. In a cold environment, individuals with a larger maturation size have a higher degree of fitness [67]. Fitness, in turn, depends on food supply, and therefore increased filtering efficiency increases individual fitness, and this has advantages in cases of low food concentration. The filtering setae number (FSN) is constant for the whole life-span of individuals [68], and FSN increases with body size [69]. Therefore larger individuals are efficient filterers of phytoplankton, due to their larger filtering screen, thus the large-bodied cladocerans are more effective grazers, since they have a lower threshold of food [70]. It seems reasonable to assume that a larger filtering screen is probably an adaptive response [71,72] in local benthic cladocerans as well. In addition, larger females produce larger ephippia and these ephippia have a larger hatching probability. Larger eggs then have the benefit of resulting in larger neonates, which in turn yield larger adults [73].
The clear increasing trend observed in the size of all remains with latitude is in agreement with the size-efficiency hypothesis, since latitudinal changes may describe a general pattern of body size changes articulated in both Bergmann's [74] and Atkinson's rules [13]. Generally, the T-S relationship would correspond to the L-S, and in aquatic species, this can be explained in terms of the oxygen supply, since warm water contains less oxygen [59,75,76]. Large body size has no advantage in a warm environment, while such species exhibit a low surface-area to body-mass ratio. Due to the increased oxygen demand of these species, they will suffer from oxygen depletion, and consequently, must sustain their growth in order to achieve a higher surface-area to body-mass ratio to increase oxygen uptake [77,78]. Therefore, smaller species are found in warm environments [76,79,80].
The strength of the associations between the T-S and L-S gradients can vary, depending on the particular species. Body size can be influenced by not only ambient temperature, but also season length, productivity, mortality and somatic plasticity, or by some combination of these [76]. Havens et al. [5] found that the body size of planktonic cyclopoid copepods had a strong relationship with temperature, while in the case of planktonic cladocerans this was weak and with calanoids it was insignificant in their comprehensive study. The discrepancy between the latitude-size (L-S) and temperature-size (T-S) models can be explained by the distribution of littoral species due to depth, hence in this research the question of whether the species would display any response to water depth was investigated. In deep and stratified lakes, the water temperature is lower, and its seasonal changes are smooth in the hypolimnion [81,82]. Therefore, populations of benthic species, which inhabit the region below the thermocline, live in a stable thermal environment. The thermal differences in the sediment surface between deep and shallow lakes may be reflected in the size of remains. Although the data in this study were unbalanced, since there were more shallow lakes in the data set than deep ones, the results show that thermal differences in littoral and profundal areas may not have a significant effect on the size distributions of the remains due to the accumulation of remains in the deepest areas. In his intra-lake studies, Luoto [83,84] (2010, 2012) found that the distribution of chironomid remains is closely related to water depth and suggested that the fact that fossils accumulate near each species' habitat may cause some bias in single-core paleolimnological studies.
A species-dependent response to temperature gradients was also observed. Direct T-S relationships were insignificant in the case of Alona species, due to the small size of the remains from the Southern Carpathians, but were strong and significant in the case of C. cf. sphaericus..
Remains of Alona species in the lakes of South Carpathians were generally smaller than those from Finland ( Figure A3). Since these Alona species are larger than C. cf. sphaericus, it was, therefore, reasonable to assume that in this case, the decrease in body size was the result of fish predation. Havens et al. (2015) demonstrated the effect of fish predation on cladoceran body size in the lakes of Greenland, with a fivefold increase in body size in the fish-free lakes than those with fish. The lengths of Alona remains from a subset of Finnish lakes inhabited by fish [85] were therefore compared with those from Southern Carpathian lakes. For the purposes of analysis, all southern mountain lakes were considered to be populated by fish, since a few lakes were naturally inhabited by trout (Salmo trutta), but almost all mountain lakes were stocked with trout between 1961 and 1977 [86]. Moreover, schools of common minnows (Phoxinus phoxinus) were seen from the shore of Carpathian lakes during the field sampling. Minnows are found in a variety of microhabitats in the littoral, but seem to prefer stony, shallow sites [87] as do most chydorids. Although minnow is the main prey of trout, the diets of young trout and minnow overlap, that is, they utilize the same food resources, and therefore both species have an effect on cladocerans [36,88,89]. Predation patterns in the littoral form a complex picture. Since both predatory fish and invertebrates utilize this zone, and fish also feed on invertebrate predators [31], the density of fish populations influences invertebrate predation pressure [29,31]. Both types of predation affect the structure of the littoral cladoceran community. Adamczuk [36] showed that fish prefer larger chydorids, and thus larger species tend to occupy dense vegetation. Basińska [38] found that Chydorus reached a large size in helophyte beds, a fact explicable in the light of the influence of macroinvertebrate predators. The scarce or absent macrophytes in the littoral zone of Carpathian lakes did not supply shelter for cladocerans, and thus predation pressure from fish explicitly affects the larger Alona species than the smaller C. sphaericus, causing a decrease in size in the Alona spp. population.
In arctic and subarctic lakes in Finland, zooplankton density is higher in lakes containing fish, and this is mainly due to the higher density of small-sized cladocerans and cyclopoids that escape grazing from high-latitude fish species. Moreover, large-sized cladocerans are found in lower densities in fish-free lakes [90,91]. Since fish prefer larger cladocerans, fish predation shifts the size spectrum of cladocerans towards the smaller range, and leads to a conclusion that fish predation enhances the effect of temperature on body size. Lakes in mid and southern Finland are exposed to nutrient enrichment [92,93], so in these lakes macrophyte density and species diversity increase. The macrophyte belt can provide shelter against predators for macroinvertebrates, and macroinvertebrate predation increases the body size of littoral cladocerans [34]. These top-down effects may well, therefore, explain the large headshield of Alonas in Southern Finnish lakes.
Alona spp. remains were smaller in the lakes of the Southern Carpathians ( Figure A3) than in fish-populated lakes in Finland, and in this case, the former was also fish-populated. Consequently, the incidence of smaller sized Alona remains could not be explained by the impact of fish. Manca and Comoli [58] studied fish-free mountain lakes in the Himalayas, where the number of moults and Daphnia body size depended on summer temperatures and on the duration of the ice-free period. They found that in cold years, when the productive season was short, the number of moults were low, and the range of body sizes was reduced. The range of size remains highly dependent on the proportion of juveniles, and thus in cases where the cladocera population finds suitable environmental conditions, a high number of juveniles will develop, causing a high degree of variance in the data. The weak correlation observed between the length of remains and ambient temperature was due to high variances. Complete T-S models of Alonas appear to contradict Atkinson's rule; however, if the Carpathian lakes are excluded, the models changed, decreasing to the point of a significant degree of agreement with the temperature-size hypothesis. We may, therefore, assume that the smaller size of the Southern mountain populations of Alona species is explained by factors other than temperature, factors which also affect size distribution. What is more, the smaller size is the result of the adaptation of a local population with a distinct genetic background.
Belyaeva and Taylor [94] revealed that C. sphaericus represents a species complex with a Holarctic distribution. They describe two European clades (A1 and A2), which look very similar, and hence cannot be distinguished through the use of carapace morphology, as the identifiable substantial difference is in gamogenetic individuals. Both clades (A1 and A2) occur in the Nordic populations of Chydorus sphaericus, but the A2 clade probably dominates at higher latitudes [94]. This A2 clade is similar to C. biovatus as described in Salmon Lake, Montana, USA [95]. They made the reasonable demand that, henceforth, further studies should be required to clarify the taxonomic status and the position of this taxon within the Chydorus sphaericus complex. Van Damme and Eggermont [96] also found that the average body size of high mountain populations of C. cf. sphaericus in the Rwenzori Mountains in the Democratic Republic of the Congo were larger than the Palearctic populations.
They suggested that the difference in size might be indicative of some degree of local adaptation and speciation (Löffler island).
In conclusion, this study reveals that littoral cladocerans exhibit a body-size response to climate over a wide geographical range across Europe. The length of headshields and carapaces may be regarded as a realistic proxy for climatic changes. Therefore, long-term changes in size distribution of their remains may provide information on climate oscillations in paleolimnological records. The analyses presented here also point to the fact that a regional heterogeneity in body size distribution which can obscure these relations may be expected to occur among populations. It may consequently be assumed that the genetic heterogeneity of species is responsible for size differences in littoral cladoceran populations.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study, in the collection, analyses, or interpretation of data, in the writing of the manuscript, or in the decision to publish the results.