All analyses were conducted in R v4.0.2 (http://www.R-project.org/). The structure of the final linear models and sample sizes break down by the levels of categorical explanatory variables are summarized in (Appendix S3). Non-significant interactions were removed from the models. Significance of the final models was evaluated using ANOVA tables using Type II and III sums of squares for models without and with significant interaction, respectively. Differences among groups were analysed using Tukey’s HSD post-hoc test. We have accounted for the ontogenetic effects by using fork length as a proxy of individuals age. Fork length has been used as explanatory variable in all models. When models contained other explanatory variables correlated to the fork length, we used residuals from a linear model between the log transformed variable of interest and fork length to prevent collinearity (see Appendix S3).