Sex demographics alter the effect of habitat structure on predation by a temporary pond specialist

Habitat structure can profoundly influence interaction strengths between predators and prey. Spatio-temporal habitat structure in temporary wetland ecosystems is particularly variable because of fluctuations in water levels and vegetation colonisation dynamics. Demographic characteristics within animal populations may also alter the influence of habitat structure on biotic interactions, but have remained untested. Here, we investigate the influence of vegetation habitat structure on the consumption of larval mosquito prey by the calanoid copepod Lovenula raynerae, a temporary pond specialist. Increased habitat complexity reduced predation, and gravid female copepods were generally more voracious than male copepods in simplified habitats. However, sexes were more similar as habitat complexity increased. Type II functional responses were exhibited by the copepods irrespective of habitat complexity and sex, owing to consistent high prey acquisition at low prey densities. Attack rates by copepods were relatively unaffected by the complexity gradient, whilst handling times lengthened under more complex environments in gravid female copepods. We demonstrate emergent effects of habitat complexity across species demographics, with predation by males more robust to differences in habitat complexity than females. For ecosystems such as temporary ponds where sex-skewed predator ratios develop, our laboratory findings suggest habitat complexity and sex demographics mediate prey risk.


Introduction
Predation is a key structuring force within ecosystems and often regulates the stability and diversity of communities (Dayton, 1971;Paine, 1980;Sih et al., 1985). The implications of predator-prey dynamics are multifaceted, with both direct lethal (densitymediated) and non-lethal (trait-mediated) effects of predators identified as important ecosystem-level processes (Dodson, 1974;Lima, 1998;Alexander et al., 2013b). Predation pressure can, however, vary markedly over space and time, yet the mechanistic drivers of such variations remain poorly understood (Delclos & Rudolf, 2011;Cuthbert et al., 2018dCuthbert et al., , 2019c. In particular, habitat structure may differ considerably across spatio-temporal gradients and is known to be an important determinant of interaction strengths through, for example, the provision of prey refugia and mediation of predatorpredator interference (Crowder & Cropper, 1982;Barrios-O'Neill et al., 2015;Cuthbert et al., 2018e, 2019a. Furthermore, heterogeneity within populations (e.g. ontogenic or sex ratio change) can alter the intensity of predator-prey interactions (Delclos & Rudolf, 2011;Alexander et al., 2013a;Cuthbert et al., 2019c). However, whether such intraspecific demographic heterogeneity alters the influence of environmental context, such as habitat structure, on interaction strengths remains unclear.
The structuring of habitats in temporary freshwater ecosystems is particularly dynamic, owing to natural variability in water depth and vegetation structure over the hydroperiod. Temporary pond ecosystems provide important aquatic habitats for a variety of amphibians, birds and invertebrates (McCulloch et al., 2003;Ferreira et al., 2012) and are characteristically inhabited by specialist, often endemic, species adapted to withstand extended dry periods (De Meester et al., 2005). In comparison with permanent freshwaters, temporary wetland ecosystems are poorly understood and remain at risk of degradation or destruction as a result of multiple anthropogenic stressors (Dalu et al., 2017a). Accordingly, there is an urgent need to examine biotic interactions in these systems in order to better understand their stability in the context of ongoing environmental change. Moreover, understanding how habitat heterogeneity affects different population demographics is important for understanding population success and ecosystem stability (Kiørboe, 2006).
Animal population demographics are extremely variable over the course of the hydroperiod in temporary ponds , with the potential to alter offtake rates towards lower trophic groups (Cuthbert et al., 2019c). For copepods, sexskewed ratios often emerge in natural populations due to processes such as predation and physico-chemical variation (Gusmão et al., 2009;Wasserman et al., 2018b). Zooplankters such as calanoid copepods are numerically dominant in temporary wetland ecosystems for much of the hydroperiod, given their internal recruitment following submersion from egg banks deposited within substrate (Wasserman et al., 2016a). As a result, large predatory copepods can occupy the top trophic levels in these systems in the early stages of the hydroperiod (Dalu et al., 2016), with other higherorder predation extensively dependent on external recruitment dynamics at the landscape scale (e.g. Wasserman et al., 2018b). Intraspecific behavioural variations between sexes in copepod populations are marked, with males often more motile than females due to mate-searching processes and time limitations to reproduction (Kiørboe, 2006;Gusmão et al., 2009). There has, however, been no consideration for whether habitat heterogeneity effects (i.e. habitat structure) manifest differently among demographic subsets of populations (i.e. according to sex variations). Consequently, our understanding of the context dependency of biotic interaction strengths in temporary ponds is lacking, particularly given the documented development of sex-skewed populations in these ecosystems .
Functional responses are fundamental to derivations of consumer-resource interactions (Solomon, 1949;Holling, 1959). Through quantifications of resource utilisation as a function of resource density, applications of the functional response offer insights into biotic interaction strengths in the context of population-level stability (Alexander et al., 2012;Dick et al., 2014). Furthermore, the functional response approach allows for explicit consideration of multiple context dependencies which may alter consumer-resource interactions, including temperature (Wasserman et al., 2018a), species demographics (Cuthbert et al., 2019c) and multiple predator interactions (Wasserman et al., 2016b;Sentis & Boukal, 2018). Functional responses are typically characterised into three forms (Hassell, 1978), with each, theoretically, pertaining to different consumer-resource outcomes. Type I functional responses are linear and density-independent and are considered mechanistically exclusive to filter feeding organisms (Jeschke et al., 2004). Type II functional responses are hyperbolic and inversely density-dependent, characterised by high consumption rates at low resource (e.g. prey) densities, whereby most, if not all, prey are consumed. Consumption rates then decelerate towards an upper asymptote. Type III functional responses are sigmoidal and density-dependent and are often driven by interactions which require significant search time when resources are rare. Here, the consumption rate initially increases before again reaching a plateau, similar to the Type II functional response. Importantly, the functional response form of consumers may be an important indicator for the stability of interactions, whereby Type II functional responses are typically considered destabilising and Type III, conversely, considered stabilising for resource populations. Moreover, it is known that environmental context (e.g. habitat structure) can impart stability to populations through mediation of more stable functional response forms (e.g. Alexander et al., 2012).
The aim of this study is to use predatory functional responses to quantify combined abiotic and biotic context dependencies of interaction strengths by an important temporary pond specialist. We focus on interactions between the copepod Lovenula raynerae, which is known to be a top predator for much of the hydroperiod before other predators (e.g. notonectids) colonise (Dalu et al., 2017b), and mosquito larval prey, which are known to externally colonise temporary aquatic systems. Previous research has demonstrated the potential role of mosquito prey in the diet of L. raynerae (Cuthbert et al., 2018b(Cuthbert et al., , 2019d. Specifically, we aim to decipher the effects of habitat structure and population sex demographics on biotic interactions, and whether the influence of habitat structure is dependent on the internal demographic characteristics within species. Through a factorial laboratory-based approach, we provide comparative insights into the trophic dynamics of poorly understood temporary pond ecosystems.

Animal collection and maintenance
The focal predators, adult L. raynerae were collected from a temporary pond close to Peddie (33°12 0 46.9 00 S 27°20 0 11.2 00 E), Eastern Cape, South Africa, using a 64-lm mesh zooplankton net and transported in source water to a constant environment (CE) room at Rhodes University, Grahamstown. Adult females and males (total length ± SE: females, 4.39 ± 0.08 mm; males, 4.13 ± 0.08 mm) were maintained in the CE room at 25 ± 1°C and under a 12:12 light:dark regime for 48 h and starved in continuously aerated 30-l tanks containing dechlorinated tap water. The focal prey, larval Culex pipiens complex mosquitoes were cultured using egg rafts collected from artificial containers within the Rhodes University campus and reared to the desired size class (see below) in the same laboratory using a diet of crushed rabbit pellets (Agricol, Port Elizabeth), supplied ad libitum.

Experimental design
Either gravid adult females or males were used in the experiment under three habitat complexity treatments. Habitat complexity was created by positioning rinsed stalks of the bulrush Schoenoplectus brachyceras (Hochst. ex A. Rich.) Lye (4-6 mm dia.) in a uniform array in experimental arenas. Low complexity comprised no stalks, medium complexity two stalks and high complexity four stalks (as per Cuthbert et al., 2019a). This complexity variation represents that displayed throughout the hydroperiod in temporary systems. Second and third instar C. pipiens complex larvae (total length ± SE: 2.87 ± 0.19 mm) were established at five prey densities (2, 4, 8, 16, 32; n = 4 per density) in 80 mL arenas of 5.6 cm diameter containing dechlorinated tap water from a continuously aerated source, 2 hours before the addition of predators. Once predators were added, they were allowed to feed undisturbed for 6 h, after which they were removed and remaining prey counted to derive numbers killed. Controls consisted of four replicates at each density and habitat treatment without predators.

Data analyses
All statistical analyses were undertaken in R v3.5.1 (R Core Development Team, 2018). Generalised linear models (GLMs) assuming a Poisson error distribution and log link were used to analyse overall prey consumption with respect to the 'habitat complexity', 'predator type' and 'prey supply' factors and their two-and three-way interactions. We did not find any evidence for residual overdispersion in the model. An information theoretic approach was followed to evaluate variables of most importance in influencing numbers of prey eaten, with model averaging via second-order Akaike's Information Criterion (AICc; adjusted for small sample size, Burnham & Anderson, 2002) applied to identify models which minimised information loss (Bartoń, 2015). The relative variable importance (RVI) was also calculated, based on the sum of AICc weights which included the focal predictor variable. Models with DAICc \ 2 were considered interchangeable. In the top model, effect sizes were derived through analysis of deviance, with Tukey's tests used post hoc for multiple pairwise comparisons using the 'multcomp' package in R (Hothorn et al., 2008). Statistical significance was considered at the 95% confidence interval.
Functional response analyses were undertaken phenomenologically (Pritchard et al., 2017). Accordingly, given that we did not empirically measure the functional response parameters (Jeschke et al., 2002), functional responses were not interpreted mechanistically and were instead used for comparative purposes among factorial experimental treatments (Alexander et al., 2012). Logistic regression considering the proportion of prey consumed as a function of prey density was used to infer FR types. A Type II functional response is characterised by a significantly negative first-order term, whilst a Type III functional response is characterised by a significantly positive first-order term followed by a significantly negative second-order term (Solomon, 1949;Holling, 1959;Juliano, 2001). As prey were not replaced as they were consumed, we applied Rogers' random predator equation for depleting prey densities (Rogers, 1972;Trexler et al., 1988;Juliano, 2001): where N e is the number of prey eaten, N 0 is the initial density of prey, a is the attack constant, h is the handling time and T is the total experimental period. The Lambert W function was used to aid model fitting, owing to the implicit nature of Eq. 1 (Bolker 2008). The attack rate parameter corresponds to the initial slope of functional response curves (Hassell & May, 1973;Jeschke et al., 2002), and thus, predators which consume more prey at low prey densities should have a higher attack rate. The handling time parameter, reciprocally (1/h), corresponds to the maximum feeding rate of a predator and the height of functional response curves (Jeschke et al., 2002). Functional response parameters (attack rate, handling time) were non-parametrically bootstrapped (n = 2000) and parameters were compared based on overlapping of 95% bias-corrected confidence intervals across treatments. Given that bootstrapping allows for data to be considered in population terms, rather than at the sample level, a lack of convergence in confidence intervals allows for differences in functional response parameters to be inferred without additional statistical testing.

Results
Survival in predator-free controls exceeded 99% and so all experimental mortality was attributed to copepod predation. The 'habitat' (RVI = 1.00), 'predator' (RVI = 0.78), 'supply' (RVI = 1.00) and 'habitat 9 predator' (RVI = 0.40) terms were included in the top model (Table 1). Overall, predation increased significantly under greater prey supplies (v 2 = 152.52, df = 4, P \ 0.001) and was significantly affected by habitat complexity (v 2 = 16.12, df = 2, P \ 0.001). Significantly greater predation was exhibited under the low complexity treatment compared to the high complexity (P \ 0.001) and medium complexity treatment (P \ 0.01); there was no significant difference between the medium and high complexity treatments (P = NS). Overall predation was not affected by predator sex (v 2 = 3.38, df = 1, P = NS); however, gravid female copepods tended to consume more than males overall (Fig. 1). The effects of habitat complexity were consistent between predator types as there was no significant 'habitat 9 predator' interaction effect (v 2 = 4.86, df = 2, P = NS); however, consumption by female copepods was generally most affected by increases in habitat complexity (Fig. 1). Indeed, post hoc contrasts indicated that consumption by males was not significantly affected by habitat complexity (all P = NS), whilst consumption by females was significantly greater under low as compared to medium (P \ 0.01) and high (P \ 0.001) complexities. Type II functional responses were observed across each habitat complexity treatment and predator sex (Table 2; Fig. 1). Functional response parameters of female and male copepods showed subtle differences across the habitat complexity gradient (Fig. 2). For female copepods, attack rates were similar across the different levels of habitat; however, these tended to be Models with high importance (DAICc \ 2) are shown here df degrees of freedom, logLik log likelihood, DAICc difference between model Akaike Information Criterion (AICc, adjusted for small sample size) and minimum AICc value, w i probability of model i being the best of this set of candidate models given data at hand greater under higher complexities (Table 2; Fig. 1). This was further evidenced by overlapping of confidence intervals between these treatment groups (Fig. 2a). Similarly, whilst habitat complexity had little effect on the attack rates of male copepods ( Fig. 1), attack rates trended towards being higher under more complex treatments (Fig. 2a). Generally, there were negligible differences between attack rates of female and male copepods across the habitat types, as indicated by extensive confidence interval overlap.
Handling times and, inversely, maximum feeding rates of female copepods varied across the habitat  Fig. 1). Handling times were lowest, and therefore maximum feeding rates were highest, under the low complexity treatment as compared to medium and high habitat complexities (Fig. 1). This significant difference was further evidenced by a lack of overlap in confidence intervals between the low and high complexity treatments for female copepods, whilst medium complexity confidence intervals bridged these treatment groups (Fig. 2b). For male copepods, handling times were similar across the habitat gradient, and thus, maximum feeding rates were relatively unaffected for this predator group (Table 2; Fig. 1). Whilst handling times tended to be longer under higher complexity treatment groups, confidence intervals overlapped across all treatments and thus differences were nonsignificant for this copepod sex (Fig. 2b). At low complexities, handling times of female copepods were substantially lower than males and generally became more similar as habitat complexity increased; however, confidence intervals of males overlapped with all female copepod complexity groups, owing to extensive interval limits.

Discussion
The present study demonstrates the combined influence of habitat structure and population sex demographics on predator-prey interaction strengths between temporary pond specialist biota. Temporary ponds are particularly dynamic aquatic habitats, with substantial habitat heterogeneity present over short timescales due to characteristic hydroperiod variation (Wasserman et al., 2016c). However, hitherto, there has been little consideration for the influence of environmental variables on interactions in these systems (but see Wasserman et al., 2016c;Cuthbert et al., 2019a), and so, our understandings of temporary wetland ecology and the potential effects of environmental change on species interactions are limited. Furthermore, sex demographics are known to shift radically over the course of the hydroperiod due to selective processes such as predation in temporary ponds , and the implications of sex ratio variations for trophic interactions in such wetlands remain unclear (but see Cuthbert et al., 2019c). Here, we demonstrate that habitat complexity alters the interaction strength between temporary pond specialists and their prey, but that the effects of habitat structure on these interactions manifest differently between sexes. Both males and females of the predaceous calanoid copepod L. raynerae were capable of capturing and handling larval mosquito prey across the habitat complexity gradient in the present study. However, consumption was significantly affected by habitat complexity. Specifically, increasing habitat complexity tended to reduce interaction strengths towards larval mosquito prey, yet only significantly for female copepods. Whilst results of laboratory experiments must be viewed with caution when making real-world inferences, comparative laboratory feeding studies can provide useful and robust insights into the effects of multiple treatments in controlled environments which corroborate field patterns (Dick et al., 2014). Complex habitats can provide refugia for prey items and act as a barrier to predator movements (Sih et al., 1992;Barrios-O'Neill et al., 2015). However, the presence of habitat can also intensify interaction strengths, depending on the specific prey capture strategies of predators (Wasserman et al., 2016c). This may be particularly true for ambush predators, which perform better in structured habitats owing to sit-and-wait strategies (James & Heck, 1994). Therefore, habitat implications for biotic interactions such as predation are system-specific. Copepods are often reliant on hydromechanical rather than visual cues for prey detection (Hwang & Strickler, 2001;Cuthbert et al., 2018a), and L. raynerae predation has been previously shown to be unaffected by variations in water clarity (Cuthbert et al., 2018c). Considering the active searching strategy utilised by L. raynerae through the water column, it is likely that negative consumptive effects driven by habitat structure resulted through interference by structures which impeded prey detection. Further, gravid female L. raynerae tended to be more voracious than males in the present study, an effect likely driven by differences in energy demands in relation to the development of progeny. This intraspecific difference corroborates with results from other zooplankton systems (e.g. Laybourn-Parry et al., 1988;Marten & Reid, 2007;Cuthbert et al., 2019c).
Both female and male L. raynerae exhibited a Type II functional response across the habitat complexity gradient, owing to high rates of consumption at low prey densities. Contrastingly, a multitude of other studies has demonstrated shifts in functional response form associated with environmental context (e.g. Koski & Johnson, 2002;Alexander et al., 2012;South et al., 2018). Type II functional responses are indicative of potential population-destabilising effects and high ecological impact, owing to a lack of consumptive refuge where prey are scarce (Dick et al., 2014). Accordingly, it is possible that predatory zooplankters such as L. raynerae can extirpate prey types such as larval mosquitoes, even where habitats are structurally complex. Despite a lack of effect on functional response types, emergent differences in attack rate and handling time parameters were evident across habitat complexities and between copepod sexes. Attack rates in both sexes were generally similar across the habitat gradient, yet tended to peak under higher complexity levels. Attack rates are classically described as the scaling coefficient of functional responses and therefore influence the initial slope of the functional response curve (Hassell & May, 1973;Jeschke et al., 2002). Accordingly, high attack rates are conducive to high predatory impact at low prey densities (Cuthbert et al., 2018e), and this further supports the lack of habitat effect on low-density prey refugia in the present study. Similarly, Wasserman et al. (2016c) investigated emergent effects between habitat structure and temperature in aquatic systems, wherein notonectid attack rates increased with structural complexity under certain conditions.
Conversely, handling times in the present study differed according to habitat complexity. The handling time of a predator corresponds to the time taken to handle and digest a prey item (Jeschke et al., 2002). Inversely, handling times correspond to the maximum feeding rate (i.e. asymptote) of functional response curves. In this study, for gravid female L. raynerae, handling times were significantly reduced under low complexity treatments as compared to high complexity treatments. Accordingly, higher maximum feeding rates (i.e. functional response asymptotes) were displayed under these conditions by the female copepods. It is likely that lower encounter rates with prey resulted in increased handling times in complex environments, and therefore, maximum feeding rates were reduced. This finding corroborates with Cuthbert et al. (2019a), wherein handling times of notonectids tended to increase, and maximum feeding rates thus fall, under greater habitat complexity levels.
Copepod sex demographics, however, mitigated the influence of habitat complexity on interaction strengths. In the absence of habitat complexity, female copepod's maximum feeding rates were substantially higher than males, which saturated at a lower rate. These effects may be driven by differences in satiation levels between the two predator sexes, with energy demands by females often higher in copepods (e.g. Laybourn-Parry et al., 1988;Cuthbert et al., 2019c). For male L. raynerae, the effects of habitat structure on handling times were less prevalent, and therefore maximum feeding rates remained relatively unaffected. Male copepods are often characterised by mate-searching processes which may impart greater vulnerability to higher-order predation (Kiørboe, 2006, Gusmão et al., 2009, and motility rates in male L. raynerae have been shown to be particularly high and unaffected by higher predator cues (Cuthbert et al., 2019b). Accordingly, it is possible that high motility rates in males reduce the influence of habitat structure on handling times via negations of physical interference by habitat structures. However, further behavioural examinations are required to elucidate this, particularly in natural settings. Furthermore, whilst gravid female L. raynerae are significantly more voracious in the absence of habitat (Cuthbert et al., 2019c), in the present study we demonstrate emergent effects of habitat structure which nullify these demographic differences. Therefore, the influence of habitat structure may be, in turn, dependent on the internal characteristics of populations. Future research to examine how emergent habitat complexities alter interactions between multiple predators would provide further insight into trophic dynamics in temporary ponds, owing to the importance of intraspecific predator-predator interactions for prey risk .
Robust quantifications of trophic interactions within ecosystems are imperative to understanding the influence of environmental change on ecosystem stability, especially in vulnerable temporary wetlands (Dalu et al., 2017a). In particular, studies such as the present which examine the interactive influence of multiple abiotic and biotic context dependencies are crucial to understanding how predation pressure is affected by environmental context. We therefore demonstrate the importance of considering demographic variations within populations when quantifying the influence of environmental changes, such as habitat alterations, for trophic interactions and population-level stabilities.
Acknowledgements This study was funded by the Department for the Economy, Northern Ireland. We extend gratitude to Rhodes University for the provision of laboratory facilities. We acknowledge the use of infrastructure and equipment provided by the South African Institute for Aquatic Biodiversity (SAIAB) Research Platform and the funding channelled through the National Research Foundation-SAIAB Institutional Support system. This study was partially funded by the National Research Foundation -South African Research Chairs Initiative of the Department of Science and Technology (Inland Fisheries and Freshwater Ecology, Grant No. 110507). We also acknowledge the Natural Environment Research Council.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.