Abstract
Predation is a key driver of phenotypic diversification with prey having evolved sets of correlated anti-predator traits. Changes in anti-predator traits can be studied on an evolutionary as well as on a developmental timescale. Using a common garden setup, we studied inter- and intraspecific correlations of behavioural and morphological defences in four damselfly species that either occur in habitats dominated by predatory fish (fish habitats) or fishless habitats by raising larvae either with predatory fish or in a control treatment. We found inter- as well as intraspecific trait compensation (negative correlations) between behavioural and morphological defences. Compared to fishless habitat species, fish habitat species invested more in behavioural defences and less in morphological defences. This was mirrored by fish habitat species investing more in behavioural defences and less in morphological defences when reared with predatory fish whereas fishless habitat species invested less in morphological defences only. Our results emphasise the role of context-specific combinations of defensive traits to avoid predation. We suggest, considering changes in multiple correlated traits on different timescales when studying the evolution of anti-predator traits.
Introduction
Predators are well-known to structure phenotypic abundance and occurrence as well as to generate species divergence (Wellborn et al., 1996; Vamosi, 2005; Schluter, 2009; Piovia-Scott et al., 2017). In response to the strong selective impact of predation, prey evolved a multitude of defensive traits to avoid and evade predation (Edmunds, 1974; Tollrian and Harvell, 1999). It is apparent that prey animals typically recruit sets of correlated behavioural and morphological anti-predator traits to avoid predation [reviewed in DeWitt and Langerhans (2003) and Bourdeau and Johansson (2012)]. Trait correlations among defensive traits have been studied in aquatic organisms like snails [e.g., Cotton et al. (2004), Bourdeau (2009), Brönmark et al. (2012), and Ahlgren et al. (2015)], odonates [e.g., Mikolajewski and Johansson (2004) and Mikolajewski et al. (2010)], amphibians [e.g., Hossie et al. (2017) and Pacheco et al. (2019)], and fish [e.g., Hulthén et al. (2014), Kern et al. (2016), and Marshall and Wund (2017)], as well as terrestrial organisms such as moths (Vogelweith et al., 2014), snakes (Mayer et al., 2016), and rodents (Agnani et al., 2020). Functionally and mechanically independent behavioural and morphological defences can either be negatively (trait compensation) or positively (trait cospecialization) correlated (DeWitt et al., 1999; DeWitt and Langerhans, 2003). In the case of trait cospecialization, prey applies both traits simultaneously to enhance the overall defence. In the case of trait compensation, prey applies only one trait, most likely to trade-off costs of defensive traits (DeWitt et al., 1999; DeWitt and Langerhans, 2003).
Habitats differ in occurrence of predators, resulting in well-distinguished predator gradients (Gunzburger and Travis, 2004; Sorace and Gustin, 2009; Renner et al., 2012; Mori and Saitoh, 2014). The expression of behavioural and morphological defences strongly depends on the current predation regime (Laurila et al., 2004; Relyea, 2004a; Bourdeau, 2009; Mikolajewski et al., 2010; Marshall and Wund, 2017). Such different phenotypes could for instance arise through developmental plasticity [i.e., intraspecific differences; e.g., Relyea (2001), Hoverman and Relyea (2009), and Herzog et al. (2016)] or can be the result of evolutionary changes [i.e., interspecific differences; e.g., Mikolajewski et al. (2006), Stoks and McPeek (2006), Langerhans (2009), and Jiang and Mikolajewski (2018)]. In cases where closely related species occupy different levels of the predator gradient, we are able to study inter- and intraspecific phenotypic changes using different predator treatments in a common garden experiment and by this to compare developmental with evolutionary timescales (Mikolajewski et al., 2015a).
Because different predator species exhibit distinct hunting modes and foraging styles, selection on prey phenotypes varies along predator gradients (Werner and McPeek, 1994; McPeek et al., 1996; Mikolajewski et al., 2006; Stoks and McPeek, 2006; Marchinko, 2009). The freshwater predator gradient is an excellent model system to study predator-related changes in correlated defensive traits. Permanent freshwater habitats exhibit a well-described predator gradient from habitats being dominated by predatory fish (hereafter referred to as “fish habitat”) and habitats that lack predatory fish [hereafter referred to as “fishless habitat”; Wellborn et al. (1996)]. Because predatory fish are highly efficient visual predators, selection on morphological and behavioural traits in fishless habitats [that are usually dominated by predatory invertebrates like large dragonfly larvae, notonectids, and beetle larvae; Wellborn et al. (1996)], is often considerably reduced in comparison to fish habitats (Stoks and De Block, 2000; Mikolajewski et al., 2015a).
Here, we investigate correlations of behavioural and morphological anti-predator traits in four species of the damselfly genus Coenagrion. Two species mainly occur in fish habitats whereas the other two species occur mainly in fishless habitats (Kuhn and Burbach, 1998; Dijkstra and Lewington, 2006). Furthermore, previous research showed that Coenagrion damselfly larvae can detect and respond to chemical cues of predatory fish (Chivers et al., 1996). Thus, by rearing larvae of all four species either with or without predatory fish, we were able to study inter- as well as intraspecific trait correlations. As a behavioural trait we quantified larval activity. Reducing activity is a common strategy in damselfly larvae to avoid detection and encounters with predatory fish (McPeek, 1990; Stoks and Johansson, 2000; Strobbe et al., 2011), but is costly because of reduced encounters with food and decelerated growth rate (Stoks et al., 2005a,b; Biro et al., 2006). As a morphological trait, we measured larval caudal lamellae size. Larval damselflies use their caudal lamellae to evade invertebrate predators by accelerating burst escape swimming (McPeek et al., 1996; McPeek, 1997; Stoks et al., 1999; Stoks and De Block, 2000; Strobbe et al., 2010). However, burst escape swimming does not operate against predatory fish as larval burst escape swimming is not fast enough to evade predatory fish (Convey, 1988; Stoks and De Block, 2000; Gyssels and Stoks, 2005). As a result, damselfly species occurring in fish habitats express smaller caudal lamellae (McPeek, 1995; McPeek et al., 1996).
Comparing fish habitat species with fishless habitat species, we expected that (a) fishless habitat species invest less into behavioural defences than fish habitat species and by this show higher activity (see above). Furthermore, we expected that (b) fish habitat species invest less into morphological defences than fishless habitat species, and by this, express smaller caudal lamellae (see above). This would result in interspecific trait compensation of behavioural and morphological defences comparing fish habitat species with fishless habitat species. Comparing larvae of one species raised under the fish treatment with larvae raised under the control treatment, we expected that (c) only fish habitat species would plastically reduce activity in the fish treatment. Further, we expected that (d) only fish habitat species would develop smaller lamellae in the fish treatment. Thus, we also predicted intraspecific trait compensation of behavioural and morphological defences within fish habitat species but not within fishless habitat species. Prediction (c) and (d) are based on previous findings that showed that freshwater prey might show reduced or even lack responses to allopatric predator species (Richardson, 2001; Stoks et al., 2003; Åbjörnsson et al., 2004; Jiang and Mikolajewski, 2018).
Materials and Methods
Study Organisms and Maintenance
In this study we focus on four European Coenagrion species [Coenagrion lunulatum (Charpentier, 1840), Coenagrion mercuriale (C., 1840), Coenagrion ornatum (Selys, 1850), and Coenagrion pulchellum (Vander Linden, 1825)] that occupy permanent waters (Dijkstra and Lewington, 2006). All four species emerge in early May (Sternberg and Buchwald, 1999) with females laying eggs during late spring and summer. Larvae hatch approximately 2–3 weeks after oviposition. Larvae overwinter and finish development in the following year. Species show distinct differences with regards to predator occurrence. Whereas C. lunulatum and C. pulchellum can dominantly be found in small permanent ponds lacking predatory fish, C. mercuriale and C. ornatum dominantly occupy low order streams that contain small predatory fish species such as three-spined sticklebacks, Gasterosteus aculeatus and European perch, Perca fluviatilis (Sternberg and Buchwald, 1999).
To test for differences between fish habitat species and fishless habitat species as well as the effect of presence/absence of predatory fish on trait correlations, we reared larvae of all four species in a full factorial experiment. Therefore, egg clutches of all four species were collected in June 2017 following the method described in Mikolajewski et al. (2015b) [n = 4; 4; 7 and 8 for C. mercuriale, C. ornatum, C. pulchellum, and C. lunulatum, respectively; sampling location for C. mercuriale and C. ornatum: Artern (51°43′ N, 11°23′ E); sampling location for C. pulchellum and C. lunulatum: Wilmersdorfer Soll (58°86′ N, 42°95′ E)]. Egg clutches were stored separately in opaque plastic containers (7 × 4 × 10 cm) filled to 8 cm with aged tap-water. All containers were kept in a climate room at 21°C and a day-night light cycle corresponding to beginning of June. Larvae of all four species hatched after ∼2 weeks. Larvae received Artemia nauplii twice daily ad libitum until they were separated. The rearing experiment took place in six black 90 L plastic containers filled with 85 L of aged-tap-water (Mikolajewski et al., 2008). Half of the containers contained one perch, P. fluviatilis (body length approx. 10 cm). Perch is a well-known predator of odonate larvae (Rask, 1986). The fish were fed dead blood worms (Chironmidae) twice a week ad libitum. In early July 2017, we randomly distributed larvae of one species among opaque plastic vials (6 × 6 × 10 cm) floating in the tubs. Each tub had one vial per species that held 20 larvae of mixed heritage at the start of the experiment. Larval density within vials corresponded to densities being observed in nature (Banks and Thompson, 1987). We placed a 5 × 5 cm piece of plastic gauze in each vial for larvae to perch on. To ensure vials were floating, they were stabilised by a ring of Styrofoam board. To ensure water exchange between containers and vials, each vial had two circular holes (Ø 6 cm) covered with plastic mesh (mesh width 500 μm). Thus, larvae in the fish treatment could sense the fish by chemical cues but were safe from predation (Ferrari et al., 2010). Containers were stored on campus of Freie Universität Berlin. To prevent complete freezing during winter, we stored containers from December 2017 to March 2018 in a climate room at 5°C. Damselfly larvae were fed ad libitum with Daphnia spec. twice a week, ensuring that vials always held some daphnids as food. Larvae stayed in their respective vials until the start of the behavioural assays (i.e., approx. 9 months). For sample sizes used in this experiment see Supplementary Table 1.
Behavioural Assays
Behavioural assays were carried out during April and May 2018 using only larvae that had reached one of the last three developmental instars before metamorphosis. Larvae of all species were of similar size based on head width (F = 1.833, DF = 3/119, P = 0.145). All behavioural trials were conducted by the same person (MM) in a climate room at a constant temperature of 20°C. For the observations, we used circular plastic bowls (Ø 15 cm). The sides of the bowls were covered with black paper and bottoms of the bowls were roughened with sandpaper to enable larval foraging. Bowls were placed 70 cm underneath a 20 W light bulb and filled to a height of 5 cm with water mixed from the tubs of the corresponding treatment. Before usage in the trials, the water was sieved through a 500 μm Nitex screen to exclude larger items. We used separate bowls and Nitex screens for either water with or without chemical cues from predatory fish. This way, we ensured that water contained or lacked chemical cues by perch. 48 h before each observation, larvae were collected from the vials and held separately in opaque plastic cups (ø 14.5 cm, filled to a height of 5.5 cm) with water mixed from the three containers of the larvae’s respective treatment. We fed the larvae with ten D. pulex and only larvae that had consumed all D. pulex after the 48 h period were used in trials.
All replicates were carried out between 10:00 and 17:00. Ten minutes before behavioural trials, we added 25 D. pulex as prey. After each replicate we exchanged the water and replaced all remaining daphnids with new individuals. Each larva was observed for 10 min. We gently introduced one larva to the bottom of the bowl and recorded behaviour after a brief adaptation time (30 s). The following six variables were recorded during observations (Jiang and Mikolajewski, 2018): (i) orientation (= total number of orientations toward prey without advancing toward the prey); (ii) advance (=total number of directed advances toward prey); (iii) strike (= total number of strikes at prey with the extendable labium); (iv) capture (=total number of successfully caught prey items); (v) walk (=total number of walks defined by movements with contact to the bottom but not in the direction of prey); and (vi) swim (=total number of swims defined by movements with no contact with the bottom but not in the direction of prey). Walk was easily distinguishable from orient and advance because larvae would adapt a lurking position whenever targetting a daphnid as prey. After the trial, each larva was stored individually in 90% alcohol for subsequent measurements (see below).
Lamellae Morphometrics
To assess differences in lamellae size, we carried out morphological measurements of larval caudal lamellae using a stereomicroscope (Olympus Microscope SZX16) connected to a CCD-colour-camera (ColorView II). Damselfly larvae possess three caudal lamellae of which the two lateral lamellae are identical to each other. Therefore, photos of the medial and left lateral lamellae were taken with millimetre paper as a size reference. Subsequently, photos were analyzed using the image processing software “ImageJ” (Abràmoff et al., 2004). We quantified lamellae morphometrics that are likely to be related to lamellae size (McPeek, 1995; Strobbe et al., 2009). In short, for the two lamellae, we recorded the (i) area, (ii) perimeter, (iii) length of major axis, (iv) the largest width perpendicular to the major axis, (v) the distance between the upper end of the lamella base to the lamella tip, and (vi) the distance between the lower end of the lamella base to the lamella tip, respectively (Figure 1). Using the same protocol, we also measured head width of each larva to standardise the response variables for body size in subsequent analyses (see below). All morphological variables were measured in mm or mm2.
FIGURE 1
Statistical Analyses
All statistical analyses were performed using R 4.0.1 (R Core Team, 2020). Prior to data analysis, we explored the data set and did not find any unexplainable outliers (Zuur et al., 2010). Because of multicollinearity among our variables, we applied a principal component analysis (PCA) on the correlation matrices for the behavioural and morphological variables, using package “psych” (Revelle, 2017). Extracted principal components (PCs) are orthogonal to each other and therefore uncorrelated (Jolliffe, 2002). PC axes were rotated using varimax method (McPeek, 1990). We extracted the first three PCs (Table 1) as they had eigenvalues >1 (Jolliffe, 2002; Supplementary Figure 1). In all subsequent analysis, we used PC scores from the extracted PCs.
TABLE 1
| PC1 (61.3%) | PC2 (17.7%) | PC3 (7.1%) | |
| (A) Lamellae Morphometrics | |||
| Area (Medial Lamella) | 0.96 | 0.19 | –0.01 |
| Perimeter (Medial Lamella) | 0.97 | 0.14 | 0.01 |
| Length of major axis (Medial Lamella) | 0.95 | 0.14 | 0.00 |
| Width of major axis (Medial Lamella) | 0.85 | 0.27 | –0.06 |
| Distance upper lamella base to lamella tip (Medial Lamella) | 0.98 | 0.14 | 0.04 |
| Distance lower lamella base to lamella tip (Medial Lamella) | 0.98 | 0.14 | 0.04 |
| Area (Lateral Lamella) | 0.97 | 0.10 | –0.01 |
| Perimeter (Lateral Lamella) | 0.98 | 0.07 | 0.01 |
| Length of major axis (Lateral Lamella) | 0.96 | 0.08 | 0.00 |
| Width of major axis (Lateral Lamella) | 0.91 | 0.21 | 0.02 |
| Distance upper lamella base to lamella tip (Lateral Lamella) | 0.97 | 0.07 | 0.03 |
| Distance lower lamella base to lamella tip (Lateral Lamella) | 0.97 | 0.07 | 0.03 |
| (B) Behavioural variables associated with foraging activity | |||
| Orientation | 0.15 | 0.89 | 0.03 |
| Advance | 0.20 | 0.84 | 0.08 |
| Strike | 0.14 | 0.78 | –0.02 |
| Capture | 0.05 | 0.89 | –0.02 |
| (C) Behavioural variables associated with movement activity | |||
| Walk | –0.06 | –0.11 | 0.80 |
| Swim | 0.08 | 0.16 | 0.79 |
The principal component (PC) loadings for the variables describing (A) lamellae morphometrics and (B,C) behavioural variables examined in this study.
Given in parentheses after the PCs is the amount of variation explained by that PC.
To analyse differences between habitats and treatments in behavioural and morphological anti-predator traits, we fitted separate linear mixed-effect models with the three PCs as response variables, using packages “lmer4” (Bates et al., 2015) and “car” (Fox and Weisberg, 2011). PCs were log10(x + 3) transformed to meet normal distribution of residuals errors. As predictor variables, we used species (C. lunulatum, C. mercuriale, C. ornatum, or C. pulchellum) nested within habitat (fish habitat or fishless habitat), treatment (fish or control), and the interaction term. We added log-transformed head width as an additive independent variable to control for body size differences (Mikolajewski et al., 2015a,2016) and used tub number as a random effect to control for differences between larvae from different tubs but of the same treatment.
Results
Principal Component Analysis
The extracted PCs captured 86.1% variation of the original data set (Table 1). PC1 (hereafter referred to as “Lamella Size”) explained 61.3% of variation and was highly positively correlated with the variables describing medial and lateral lamellae size metrics (Table 1A). PC2 (hereafter referred to as “Foraging Activity”) explained 17.7% of variation and was highly positively correlated with the variables orient, advance, strike, and capture (Table 1B). PC3 (hereafter referred to as “Movement Activity”) explained 7.1% of variation and was highly positively correlated with the variables walk and swim (Table 1C).
Interspecific Differences
Overall, lamellae of fish habitat species were smaller compared to lamellae of fishless habitat species (Wald-χ2 = 126.789, DF = 1, P < 0.001; Figure 2A). Foraging activity (Wald-χ2 = 5,293, DF = 1, P = 0.021; Figure 2B) but not movement activity (Wald-χ 2 = 2.360, DF = 1, P = 0.124; Figure 2C) was lower in fish habitat species than in fishless habitat species. This resulted in a pattern of trait compensation of foraging activity with lamella size comparing fish- and fishless habitat species. Here, fish habitat species invest in reduced foraging but express smaller lamellae whereas fishless habitat species show lower investment in behaviour by expressing higher foraging activity but also exhibit larger lamellae. Within fishless habitat species, C. pulchellum had larger lamellae than C. lunulatum and within fish habitat species C. mercuriale had larger lamellae than C. ornatum (Wald-χ2 = 6.916, DF = 2, P = 0.031; Figure 2A). No other significant two- or three-way interaction has been detected (see Supplementary Table 2).
FIGURE 2
Intraspecific Differences
Overall, in the fish treatment individuals expressed smaller lamellae than in the control (Wald-χ2 = 9.289, DF = 1, P = 0.002; Figure 2A) which did not differ between individuals from fishless and fish habitats (Wald-χ2 < 0.001, DF = 1, P = 0.991; Figure 2A). Foraging activity did not differ between treatments (Wald-χ2 = 0.402, DF = 1, P = 0.526; Figure 2B) which was the case for individuals from fishless and fish habitats (Wald-χ2 = 0.219, DF = 1, P = 0.640, Figure 2B). However, individuals reduced movement activity in the fish compared to the control treatment (Wald-χ2 = 6.164, DF = 1, P = 0.013; Figure 2C), but this was only true for individuals stemming from fish habitats (Wald-χ2 = 4.994, DF = 1, P = 0.025; Figure 2C). This resulted in a pattern of trait compensation of movement activity with lamella size within fish habitat species. Here, individuals of fish habitat species in the fish treatment, show higher investment in behaviour by reducing movement, but also express smaller lamellae in the presence of predatory fish. This pattern was not present in fishless habitat species as individuals of these species did not show reduced movement activity. No other significant two- or three-way interaction has been detected (see Supplementary Table 2).
Discussion
Selection by predators on prey defensive traits plays a key role in promoting the diversity of prey phenotypes [reviewed in Langerhans (2007) and Schluter (2009)]. To avoid and/or repel predators, prey relies on sets of correlated defensive traits (DeWitt and Langerhans, 2003; Pigliucci, 2003; Relyea, 2004b). Here, we show that larvae of damselfly species, occurring in fish habitats, express lower foraging activity and smaller caudal lamellae than species from fishless habitats, resulting in interspecific trait compensation. This pattern of trait correlation was mirrored by intraspecific differences in individuals of fish habitat species, expressing lower movement activity and smaller caudal lamellae in the presence of predatory fish than in the control whereas individuals of fishless habitat species expressed smaller lamellae only. The convergence of developmental and evolutionary trajectories is not surprising given the identical selective environmental conditions, with predation being a major driver of phenotypic appearance.
By investigating trait correlations between closely related species, we were able to study trait correlation on an evolutionary timescale (Mikolajewski et al., 2016; Reilly et al., 2016; Messier et al., 2017; Coley et al., 2018). Our results indicate interspecific trait compensation in morphological and behavioural defences in Coenagrion species. Both trait compensation (Abrahams, 1995; Relyea, 2000; Rundle and Brönmark, 2001; Mikolajewski and Johansson, 2004; Vogelweith et al., 2014) and trait cospecialization (Smith and Van Buskirk, 1995; Relyea, 2001; Mikolajewski and Johansson, 2004; Dayton et al., 2005; Mikolajewski et al., 2010) have frequently been reported on an interspecific level. Cospecialization is expected when strengthening both defences is more effective in reducing mortality caused by predation than strengthening a single defence would be (Garcia and Sih, 2003). Conversely, compensation is expected when enhancing one defensive trait produces a strong defence while trading off costs for other defences (Mikolajewski and Johansson, 2004).
It is unsurprising that in fish habitat species, lamella size is traded-off for foraging activity. Increasing lamella size represents an ineffective defence as larval burst escape swimming is not fast enough to outswim predatory fish (Convey, 1988; Stoks and De Block, 2000; Gyssels and Stoks, 2005; Strobbe et al., 2010). Reducing foraging activity on the other hand is an effective defence as it reduces the encounter rate with predatory fish and predation by fish therefore imposes strong selection on foraging activity (Werner and Anholt, 1993; Stoks and McPeek, 2006; Strobbe et al., 2011; Jiang and Mikolajewski, 2018). Conversely, fishless habitat species trade-off foraging activity for lamella size. In fishless habitats, predatory invertebrates such as large dragonfly larvae are the main predators of damselfly larvae (Wellborn et al., 1996; Stoks and McPeek, 2006). Most invertebrate predators can be evaded by larval burst escape swimming (Stoks, 1998; Bose and Robinson, 2013). Therefore, the underlying morphological trait of burst escape swimming, lamella size, represents an effective defence in this habitat type and selection by invertebrate predators acts strongly on lamella size in fishless habitat species [e.g., McPeek (1995, 1997) and Strobbe et al. (2009)]. Investment into behavioural defences (i.e., reducing foraging activity), however, is lower in fishless habitat species. Reducing activity decreases encounter rates with predators (Werner and Anholt, 1993; Stoks and McPeek, 2006; Strobbe et al., 2011; Jiang and Mikolajewski, 2018), but carries the cost of reduced encounter rates with food and reduced growth rate (Stoks et al., 2005a,b; Biro et al., 2006). Increased activity is common in freshwater prey facing invertebrate predators compared to prey facing predatory fish [e.g., Richardson (2001), Stoks and McPeek (2006), and Jiang and Mikolajewski (2018)]. This prey usually relies on other defensive traits, such as larger appendages like fins or lamellae to accelerate burst escape swimming (Wellborn et al., 1996). While these morphological defences carry costs as well [e.g., Tollrian (1995), Van Buskirk and Schmidt (2000), Benard and Fordyce (2003), Flenner et al. (2009), Van Donk et al. (2011), and Vinterstare et al. (2019)], reducing foraging activity is likely costlier [e.g., Steiner and Pfeiffer (2009)] and fishless habitat species instead rely on lamella size, underlying burst escape swimming, instead of reducing foraging. Thus, differences in predator-driven selection and costs associated with defensive traits led to this pattern of interspecific trait compensation in Coenagrion species with fishless habitat species trading off behavioural avoidance for larger lamellae and vice versa in fish habitat species.
By rearing larvae of all four species either with or without a predatory fish, we were able to study intraspecific trait correlation mediated by developmental plasticity [e.g., Relyea (2001); Hoverman and Relyea (2009), and Herzog et al. (2016)]. Interestingly, only individuals of fish habitat species showed plastic reduction in movement activity whereas all species plastically expressed smaller caudal lamellae in the presence of predatory fish. This resulted in a pattern of intraspecific trait compensation for fish habitat species which mirrored the correlational pattern we found comparing species from different habitat types. Conversely, no pattern of intraspecific correlation was observed in fishless habitat species (see below).
Our results correspond to previous work reporting that prey species respond behaviourally to co-existing predator species while being oblivious of allopatric predator species [e.g., Stoks et al. (2003), Åbjörnsson et al. (2004), Jiang and Mikolajewski (2018), and Kadye et al. (2020)], supporting the prey naïvety-hypothesis. However, all Coenagrion species regardless of habitat plastically reduced lamellae size in the fish treatment, resulting in intraspecific trait compensation in fish habitat species and no observable trait correlation in fishless habitat species. This hints at functional and mechanistic differences in morphological and behavioural responses to predators. For example, Dijk et al. (2016) found that tadpoles with reduced morphological defences (i.e., smaller tail fins) but with no differences in behaviour were less likely to be attacked by a predator, possibly due to the “lure effect” [see also Crowder and Cooper (1982), Morin (1984a,b), and Robinson et al. (1991)]. Offering an alternative explanation, Relyea et al. (2021) showed that plasticity might evolve more quickly in morphological defences than in behavioural defences. Likewise, Hossie et al. (2017) found that the strength of morphological and/or behavioural anti-predator responses in tadpoles strongly differed with exposure to different predators and whether predator and prey shared a co-evolutionary history or not. In the case of Coenagrion damselfly larvae, fish habitat species increase behavioural avoidance while saving on a morphological defence. Conversely, fishless habitat species only save on a presumably maladaptive morphological defence in response to a predator that is only rarely encountered in the specific natural habitats (Kuhn and Burbach, 1998; Dijkstra and Lewington, 2006). Thus, we show that intraspecific correlations between anti-predator traits may vary among species, again driven by the predator gradient (Wellborn et al., 1996).
Our results indicate intraspecific trait compensation in fish habitat species, contrasting previous results of a meta-analysis by Hossie et al. (2017) on intraspecific anti-predator trait correlation in tadpoles. Hossie et al. (2017) showed that tadpoles most frequently cospecialise behavioural (i.e., activity reduction) with morphological (i.e., inducing deeper or larger tail fins) defences as a response to perceived predation risk by predatory fish. Both prey types show similar responses on the behavioural axis by reducing activity, most likely because visually hunting predators such as fish are attracted by movements (McPeek, 1990; Skelly, 1994; Wellborn et al., 1996; Stoks and De Block, 2000; Strobbe et al., 2011). However, they show diverging responses on the morphological axis. Tadpoles can escape predatory fish by swimming and thus, the underlying morphological trait: tail size, operates against predatory fish (Teplitsky et al., 2005). Damselfly larvae, however, cannot reach swimming speeds high enough to escape predatory fish by relying on burst escape swimming, and therefore, caudal lamellae size does not operate against predatory fish in damselfly larvae (Convey, 1988; Stoks and De Block, 2000; Gyssels and Stoks, 2005). Thus, functional differences in morphological defences might have resulted in diverging correlational patterns in both prey taxa with tadpoles showing trait cospecialization and damselfly larvae showing trait compensation.
Intriguingly, McPeek (1997) and Stoks et al. (1999) [see also Strobbe et al. (2009) and Strobbe et al. (2010)] reported the lack of plasticity in lamellae size in damselfly larvae of the genera Enallagma and Lestes, contrasting to our results. However, Black et al. (2019) reported functional plasticity in breaking joint sizes in two Enallagma-species, by investigating intraspecific differences driven by predator cues. Strikingly, the studies of McPeek (1997), Stoks et al. (1999), and this study is the experimental duration [1 week, 1 month, 9 months, for McPeek (1997), Stoks et al. (1999), and this study, respectively; the experimental duration in Black et al. (2019) was roughly 2 months]. Thus, experimental duration in the studies by McPeek (1997) and Stoks et al. (1999) might have been too short, to induce an effect of predator presence on lamellae size in damselfly larvae [see also DeWitt et al. (1998) and Diggle (2002)].
Conclusion
Phenotypic changes are frequently studied either on evolutionary timescales (i.e., comparing populations or species) or on developmental timescales (i.e., comparing individuals raised under different conditions). Thereby, research might ignore the interactive effect of natural selection and phenotypic plasticity on shaping phenotypic diversity [e.g., Schlichting (2021)]. Moreover, research on phenotypic changes is often limited only to single traits while most traits with a similar function likely change in concert and a single-trait focus could limit our understanding of phenotypic change [e.g., DeWitt and Langerhans (2003)]. Our research is important in linking different timescales with trait correlations in a freshwater predator gradient. Although our sample size was rather limited, our findings indicate that when studying the evolution of traits, it is important to consider a multiple timescale and trait approach to better understand how selection drives phenotypic divergence.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: Open Science Framework (OSF): https://osf.io/qxzuk/.
Author contributions
BJ, DM, and MM conceived and designed the study with assistance from OB. BJ and DM collected the egg clutches for this study with assistance from OB. MM and BJ performed the study and collected the data. MM, BJ, and DM analyzed the data and MM and DM wrote the draft of the manuscript with assistance from BJ. All authors contributed to the article and approved the submitted version.
Funding
The China Scholarship Council (CSC) provided a Ph.D.-stipend for BJ.
Acknowledgments
We would like to thank M. Friedrichs for collecting the perch. We acknowledge support by the OpenAccess Publication Fund of Freie Universität Berlin for paying the APC. Two reviewers helped improve the quality of this manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fevo.2022.874276/full#supplementary-material
References
1
ÅbjörnssonK.HanssonL. A.BrönmarkC. (2004). Responses of prey from habitats with different predator regimes: local adaptation and heritability.Ecology851859–1866. 10.1890/03-0074
2
AbrahamsM. V. (1995). The interaction between antipredator behaviour and antipredator morphology: experiments with fathead minnows and brook sticklebacks.Can. J. Zool.732209–2215. 10.1139/z95-261
3
AbràmoffM. D.MagalhãesP. J.RamS. J. (2004). Image processing with ImageJ.J. Biophoton.1136–42.
4
AgnaniP.ThomsonJ.SchradinT.CareauV. (2020). The fast and the curious II: performance, personality, and metabolism in Karoo bush rats.Behav. Ecol. Sociobiol.74:123. 10.1007/s00265-020-02908-y
5
AhlgrenJ.ChapmanB. B.NilssonP. A.BrönmarkC. (2015). Individual boldness is linked to protective shell shape in aquatic snails.Biol. Lett.11:20150029. 10.1098/rsbl.2015.0029
6
BanksM. J.ThompsonD. J. (1987). Regulation of damselfly populations: the effects of larval density on larval survival, development rate and size in the field.Freshw. Biol.17357–365. 10.1111/J.1365-2427.1987.TB01055.X
7
BatesD.MächlerM.BolkerB.WalkerS. (2015). Fitting linear mixed-effects models using lme4.J. Stat. Softw.671–48. 10.18637/jss.v067.i01
8
BenardM. F.FordyceJ. A. (2003). Are induced defenses costly? Consequences of predator-induced defenses in western toads, Bufo boreas.Ecology8468–78. 10.2307/3107998
9
BiroP. A.AbrahamsM. V.PostJ. R.ParkinsonE. A. (2006). Behavioural trade-offs between growth and mortality explain evolution of submaximal growth rates.J. Anim. Ecol.751165–1171.
10
BlackK. L.FudgeD.JarvisW. M.RobinsonB. W. (2019). Functional plasticity in lamellar autotomy by larval damselflies in response to predatory larval dragonfly cues.Evol. Ecol.33257–272. 10.1007/s10682-019-09979-y
11
BoseA. P. H.RobinsonB. W. (2013). Invertebrate predation predicts variation in an autotomy-related trait in larval damselfly.Evol. Ecol.2727–38. 10.1007/s10682-012-9581-3
12
BourdeauP. E. (2009). Prioritized phenotypic responses to combined predators in a marine snail.Ecology901659–1669. 10.1890/08-1653.1
13
BourdeauP. E.JohanssonF. (2012). Predator-induced morphological defences as by-products of prey behaviour: a review and prospectus.Oikos1211175–1190. 10.2307/23261003
14
BrönmarkC.LakowitzT.NilssonP. A.AhlgrenJ.LennartsdotterC.HollanderJ. (2012). Costs of inducible defence along a resource gradient.PLoS One7:e30467. 10.1371/journal.pone.0030467
15
ChiversD. P.WisendenB. D.SmithR. J. F. (1996). Damselfly larvae learn to recognize predators from chemical cues in the predator’s diet.Anim. Behav.52315–320. 10.1006/anbe.1996.0177
16
ColeyP. D.EndaraM. J.KursarT. A. (2018). Consequences of interspecific variation in defenses and herbivore host choice for the ecology and evolution of Inga, a speciose rainforest tree.Oecologia187361–337. 10.1007/s00442-018-4080-z
17
ConveyP. (1988). Competition for perches between larval damselflies: the influence of perch use on feeding efficiency, growth rate and predator avoidance.Freshw. Biol.1915–28. 10.1111/j.1365-2427.1988.tb00323.x
18
CottonP. A.RundleS. D.SmythK. E. (2004). Trait compensation in marine gastropods: shell shape, avoidance behavior, and susceptibility to predation.Ecology851581–1584. 10.1890/03-3104
19
CrowderL. B.CooperW. E. (1982). Habitat structural complexity and the interaction between bluegills and their prey.Ecology631802–1813. 10.1007/s00442-004-1644-x
20
DaytonG. H.SaenzD.BaumK. A.LangerhansR. B.DeWittT. J. (2005). Body shape, burst speed and escape behavior of larval anurans.Oikos111582–591. 10.1111/j.1600-0706.2005.14340.x
21
DeWittT. J.LangerhansR. B. (2003). Multiple prey traits, multiple predators: keys to understanding complex community dynamics.J. Sea Res.49143–155. 10.1016/S1385-1101(02)00220-4
22
DeWittT. J.SihA.HuckoJ. A. (1999). Trait compensation and cospecialization in a freshwater snail: size, shape and antipredator behaviour.Anim. Behav.58397–407. 10.1006/anbe.1999.1158
23
DeWittT. J.SihA.WilsonD. S. (1998). Costs and limits of phenotypic plasticity.Trends Ecol. Evol.1377–81. 10.1016/S0169-5347(97)01274-3
24
DiggleP. K. (2002). A developmental morphologist’s perspective on plasticity.Evol. Ecol.16267–283. 10.1023/A:1019680527788
25
DijkB.LaurilaA.OrizaolaG.JohanssonF. (2016). Is one defence enough? Disentangling the relative importance of morphological and behavioural predator-induced defences.Behav. Ecol. Sociobiol.70237–246. 10.1007/s00265-015-2040-8
26
DijkstraK.LewingtonR. (2006). Field Guide to the Dragonflies of Britain and Europe.London: British Wildlife Publishing.
27
EdmundsM. (1974). Defence in Animals: A Survey of Anti-Predator Defences.Harlow: Longman Publishing Group.
28
FerrariM. C.WisendenB. D.ChiversD. P. (2010). Chemical ecology of predator–prey interactions in aquatic ecosystems: a review and prospectus.Can. J. Zool.88698–724. 10.1139/Z10-029
29
FlennerI.OlneK.SuhlingF.SahlénG. (2009). Predator-induced spine length and exocuticle thickness in Leucorrhinia dubia (Insecta: Odonata): a simple physiological trade-off?Ecol. Entomol.34735–740. 10.1111/j.1365-2311.2009.01129.x
30
FoxJ.WeisbergS. (2011). An R Companion to Applied Regression, 2nd Edn. London: Sage.
31
GarciaT. S.SihA. (2003). Color change and color-dependent behavior in response to predation risk in the salamander sister species Ambystoma barbouri and Ambystoma texanum.Oecologia137131–139. 10.1007/s00442-003-1314-4
32
GunzburgerM. S.TravisJ. (2004). Evaluating predation pressure on green treefrog larvae across a habitat gradient.Oecologia140422–429. 10.1007/s00442-004-1610-7
33
GysselsF. G. M.StoksR. (2005). Threat-sensitive responses to predator attacks in a damselfly.Ethology111411–423. 10.1111/j.1439-0310.2005.01076.x
34
HerzogQ.RabusM.RibeiroB. W.LaforschC. (2016). Inducible Defenses with a” Twist”: Daphnia barbata abandons bilateral symmetry in response to an ancient predator.PLoS One11:e0148556. 10.1371/journal.pone.0148556
35
HossieT.LandoltK.MurrayD. L. (2017). Determinants and co-expression of anti-predator responses in amphibian tadpoles: a meta-analysis.Oikos126173–184. 10.1111/oik.03305
36
HovermanJ. T.RelyeaR. A. (2009). Survival trade-offs associated with inducible defences in snails: the roles of multiple predators and developmental plasticity.Funct. Ecol.231179–1188. 10.1111/j.1365-2435.2009.01586.x
37
HulthénK.ChapmanB. B.NilssonP. A.HollanderJ.BrönmarkC. (2014). Express yourself: bold individuals induce enhanced morphological defences.Proc. R. Soc. B Biol. Sci.281:20132703. 10.1098/rspb.2013.2703
38
JiangB.MikolajewskiD. J. (2018). Shift in predation regime mediates diversification of foraging behaviour in a dragonfly genus.Ecol. Entomol.43525–533. 10.1111/een.12530
39
JolliffeI. T. (2002). Principal Component Analysis, Second Edn. New York, NY: Springer.
40
KadyeW. T.LeighS.BoothA. J. (2020). Predator naïve minnows respond to their conspecific alarm substance but not the odour from a non-native predator.Afr. J. Ecol.58757–765. 10.1111/aje.12768
41
KernE. M. A.RobinsonD.GassE.GodwinJ.LangerhansR. B. (2016). Correlated evolution of personality, morphology and performance.Anim. Behav.11779–86. 10.1016/j.anbehav.2016.04.007
42
KuhnK.BurbachK. (1998). Libellen in Bayern.Stuttgart: Bayerisches Landesamt für Umweltschutz & Bund Naturschutz in Bayern.
43
LangerhansR. B. (2007). Evolutionary Consequences of Predation: Avoidance, Escape, Reproduction, and Diversification. Predation in Organisms.Berlin: Springer.
44
LangerhansR. B. (2009). Trade-off between steady and unsteady swimming underlies predator-driven divergence in Gambusia affinis.J. Evol. Biol.221057–1075. 10.1111/j.1420-9101.2009.01716.x
45
LaurilaA.Järvi-LaturiM.PakkasmaaS.MeriläJ. (2004). Temporal variation in predation risk: stage-dependency, graded responses and fitness costs in tadpole antipredator defences.Oikos10790–99. 10.1111/j.0030-1299.2004.13126.x
46
MarchinkoK. B. (2009). Predation’s role in repeated phenotypic and genetic divergence of armor in threespine stickleback.Evolution63127–138. 10.1111/j.1558-5646.2008.00529.x
47
MarshallC. A.WundM. A. (2017). The evolution of correlations between behavioural and morphological defence in Alaskan threespine stickleback fish (Gasterosteus aculeatus): evidence for trait compensation and co-specialization.Evol. Ecol. Res.18305–322.
48
MayerM.ShineR.BrownG. P. (2016). Bigger babies are bolder: effects of body size on personality of hatchling snakes.Behaviour153313–323. 10.1163/1568539x-00003343
49
McPeekM. A. (1990). Behavioral differences between Enallagma Species (Odonata) Influencing differential vulnerability to predators.Ecology711714–1726. 10.2307/1937580
50
McPeekM. A. (1995). Morphological evolution mediated by behavior in the damselflies of two communities.Evolution49749–769. 10.1111/j.1558-5646.1995.tb02311.x
51
McPeekM. A. (1997). Measuring phenotypic selection on an adaptation: lamallae of damselflies experiencing dragonfly predation.Evolution51459–466. 10.1111/j.1558-5646.1997.tb02433.x
52
McPeekM. A.SchrotA. K.BrownJ. M. (1996). Adaptation to predators in a new community: swimming performance and predator avoidance in damselflies.Ecology77617–629. 10.2307/2265635
53
MessierJ.LechowiczM. J.McGillB. J.ViolleC.EnquistB. J. (2017). Interspecific integration of trait dimensions at local scales: the plant phenotype as an integrated network.J. Ecol.1051775–1790. 10.1111/1365-2745.12755
54
MikolajewskiD. J.BlockM. D.RolffJ.JohanssonF.BeckermanA. P.StoksR. (2010). Predator-driven trait diversification in a dragonfly genus: covariaton in behavioral and morphological antipredator defense.Evolution643327–3335. 10.1111/j.1558-5646.2010.01078.x
55
MikolajewskiD. J.JohanssonF. (2004). Morphological and behavioral defenses in dragonfly larvae: trait compensation and cospecialization.Behav. Ecol.15614–620. 10.1093/beheco/arh061
56
MikolajewskiD. J.JohanssonF.WohlfahrtB.StoksR. (2006). Invertebrate predation selects for the loss of a morphological antipredator trait.Evolution601306–1310. 10.1111/j.0014-3820.2006.tb01208.x
57
MikolajewskiD. J.RüsenL.MauersbergerR.JohanssonF.RolffJ. (2015a). Relaxed predation results in reduced phenotypic integration in a suite of dragonflies.J. Evol. Biol.281354–1363. 10.1111/jeb.12658
58
MikolajewskiD. J.De BlockM.StoksR. (2015b). The interplay of adult and larval time constraints shapes species differences in larval life history.Ecology961128–1138. 10.1890/14-0262.1
59
MikolajewskiD. J.ScharnweberK.JiangB.LeichtS.MauersbergerR.JohanssonF. (2016). Changing the habitat: the evolution of intercorrelated traits to escape from predators.J. Evol. Biol.291394–1405. 10.1111/jeb.12879
60
MikolajewskiD. J.StoksR.RolffJ.JoopG. (2008). Predators and cannibals modulate sex-specific plasticity in life-history and immune traits.Funct. Ecol.22114–120. 10.1111/j.1365-2435.2007.01357.x
61
MoriT.SaitohT. (2014). Flood disturbance and predator–prey effects on regional gradients in species diversity.Ecology95132–141. 10.1890/13-0914.1
62
MorinP. J. (1984a). The impact of fish exclusion on the abundance and species composition of larval odonates: results of short-term experiments in a North Carolina farm pond.Ecology6553–60. 10.2307/1939457
63
MorinP. J. (1984b). Odonate guild composition: experiments with colonization history and fish predation.Ecology651866–1873. 10.2307/1937784
64
PachecoE. O.Almeida-GomesM.SantanaD. J.GuarientoR. D. (2019). Space use and phenotypic plasticity in tadpoles under predation risk.Hydrobiology83777–86. 10.1007/s10750-019-3962-3
65
PigliucciM. (2003). Phenotypic integration: studying the ecology and evolution of complex phenotypes.Ecol. Lett.6265–272. 10.1046/j.1461-0248.2003.00428.x
66
Piovia-ScottJ.YangL. H.WrightA. N. (2017). Temporal variation in trophic cascades.Annu. Rev. Ecol. Eol. S.48281–300. 10.1146/annurev-ecolsys-121415-032246
67
R Core Team (2020). R: A Language and Environment for Statistical Computing.Vienna: R Foundation for Statistical Computing.
68
RaskM. (1986). The diet and diel feeding activity of perch, Perca fluviatilis L., in a small lake in southern Finland.Ann. Zool. Fenn.2349–56.
69
ReillyS. M.MontuelleS. J.SchmidtA.KrauseC.NaylorE.EssnerR. L.Jr. (2016). Functional evolution of jumping in frogs: interspecific differences in take-off and landing.J. Morph.277379–393. 10.1002/jmor.20504
70
RelyeaR. A. (2000). Trait-mediated indirect effects in larval anurans: reversing competition with the threat of predation.Ecology812278–2289. 10.2307/177114
71
RelyeaR. A. (2001). Morphological and behavioral plasticity of larval anurans in response to different predators.Ecology82523–540. 10.1890/0012-9658(2001)082[0523:mabpol]2.0.co;2
72
RelyeaR. A. (2004a). Fine-tuned phenotypes: tadpole plasticity under 16 combinations of predators and competitors.Ecology85172–179. 10.1890/03-0169
73
RelyeaR. A. (2004b). “Integrating phenotypic plasticity when death is on the line,” in Phenotypic Integration: Studying the Ecology and Evolution of Complex Phenotypes, edsPigliucciM.PrestonK. A. (Oxford: Oxford University Press).
74
RelyeaR. A.StephensP. R.HammondJ. I. (2021). Phylogenetic patterns of trait and trait plasticity evolution: insights from tadpoles.Evolution752568–2588. 10.1111/evo.14338
75
RennerM.ArimitsuM. L.PiattJ. F. (2012). Structure of marine predator and prey communities along environmental gradients in a glaciated fjord.Can. J. Fish. Aquat. Sci.692029–2045. 10.1139/f2012-117
76
RevelleW. R. (2017). psych: Procedures for Personality and Psychological Research. R package version 2.2.3. Available online at: https://cran.r-project.org/web/packages/psych/index.html(accessed November 5, 2018).
77
RichardsonJ. M. L. (2001). A comparative study of activity levels in larval anurans and response to the presence of different predators.Behav. Ecol.1251–58. 10.1093/oxfordjournals.beheco.a000378
78
RobinsonJ. V.ShafferL. R.HagemierD. D.SmatreskN. J. (1991). The ecological role of caudal lamellae loss in the larval damselfly, Ischnura posita (Hagen)(Odonata: Zygoptera).Oecologia871–7. 10.1007/BF00323773
79
RundleS. D.BrönmarkC. (2001). Inter- and intraspecific trait compensation of defence mechanisms in freshwater snails.Proc. Biol. Sci.2681463–1468. 10.1098/rspb.2001.1682
80
SchlichtingC. D. (2021). Plasticity and Evolutionary Theory: Where We Are and Where We Should Be Going. Phenotypic Plasticity & Evolution.Boca Raton, FL: CRC Press.
81
SchluterD. (2009). Evidence for ecological speciation and its alternative.Science323737–741. 10.1126/science.1160006
82
SkellyD. K. (1994). Activity level and the susceptibility of anuran larvae to predation.Anim. Behav.47465–468. 10.1006/anbe.1994.1063
83
SmithD. C.Van BuskirkJ. (1995). Phenotypic design, plasticity, and ecological performance in two tadpole species.Am. Nat.145211–233. 10.1086/285737
84
SoraceA.GustinM. (2009). Distribution of generalist and specialist predators along urban gradients.Landsc. Urban. Plan.90111–118. 10.1016/j.landurbplan.2008.10.019
85
SteinerU. K.PfeifferT. (2009). Optimizing time and resource allocation trade-offs for investment into morphological and behavioral defense.Am. Nat.169118–129. 10.1086/509939
86
SternbergK.BuchwaldR. (1999). Libellen Baden-Württembergs.Stuttgart: Ulmer.
87
StoksR. (1998). Effect of lamellae autotomy on survival and foraging success of the damselfly Lestes sponsa (Odonata: Lestidae).Oecologia117443–448. 10.1007/s004420050679
88
StoksR.De BlockM. (2000). The influence of predator species and prey age on the immediate survival value of antipredator behaviours in a damselfly.Arch. Hydrobiol.147417–430. 10.1127/ARCHIV-HYDROBIOL/147/2000/417
89
StoksR.De BlockM.McPeekM. A. (2005a). Alternative growth and energy storage responses to mortality threats in damselflies.Ecol. Lett.81307–1316. 10.1111/j.1461-0248.2005.00840.x
90
StoksR.De BlockM.Van De MeutterF.JohanssonF. (2005b). Predation cost of rapid growth: behavioural coupling and physiological decoupling.J. Anim. Ecol.74708–715. 10.1111/j.1365-2656.2005.00969.x
91
StoksR.De BlockM.Van GossumH.De BruynL. (1999). Phenotypic shifts caused by predation: selection or life-history shifts?Evol. Ecol.13115–129. 10.1023/a:1006656616539
92
StoksR.JohanssonF. (2000). Trading off mortality risk against foraging effort in damselflies that differ in life cycle length.Oikos91559–567. 10.1034/j.1600-0706.2000.910318.x
93
StoksR.McPeekM. A. (2006). A tale of two diversifications: reciprocal habitat shifts to fill ecological space along the pond permanence gradient.Am. Nat.168S50–S72. 10.1086/509045
94
StoksR.McPeekM. A.MitchellJ. L. (2003). Evolution of prey behavior in response to changes in predation regime: damselflies in fish and dragonfly lakes.Evolution57574–585. 10.1554/0014-3820(2003)057[0574:EOPBIR]2.0.CO;2
95
StrobbeF.McPeekM.De BlockM.De MeesterL.StoksR. (2009). Survival selection on escape performance and its underlying phenotypic traits: a case of many-to-one mapping.J. Evol. Biol.221172–1182. 10.1111/j.1420-9101.2009.01733.x
96
StrobbeF.McPeekM. A.De BlockM.StoksR. (2010). Survival selection imposed by predation on a physiological trait underlying escape speed.Funct. Ecol.241306–1312. 10.1111/j.1365-2435.2010.01752.x
97
StrobbeF.McPeekM. A.De BlockM.StoksR. (2011). Fish predation selects for reduced foraging activity.Behav. Ecol. Sociobiol.65241–247. 10.1007/s00265-010-1032-y
98
TeplitskyC.PlenetS.LénaJ.-P.MermetN.MaletE.JolyP. (2005). Escape behavior and ultimate causes of specific induced defences in an anuran tadpole. J. Evol. Biol. 18, 180–190. 10.1111/j.1420-9101.2004.00790.x
99
TollrianR. (1995). Predator-induced morphological defenses: costs, life history shifts, and maternal effects in Daphnia pulex.Ecology761691–1705. 10.2307/1940703
100
TollrianR.HarvellC. D. (1999). The Ecology and Evolution of Inducible Defenses.Princeton: Princeton University Press.
101
VamosiS. M. (2005). On the role of enemies in divergence and diversification of prey: a review and synthesis.Can. J. Zool.83894–910. 10.1139/z05-063
102
Van BuskirkJ.SchmidtB. R. (2000). Predator-induced phenotypic plasticity in larval newts: trade-offs, selection, and variation in nature.Ecology813009–3028. 10.2307/177397
103
Van DonkE.IanoraA.VosM. (2011). Induced defences in marine and freshwater phytoplankton: a review.Hydrobiologia6683–19. 10.1007/s10750-010-0395-4
104
VinterstareJ.HegemannA.NilssonP. A.HulthénK.BrönmarkC. (2019). Defence versus defence: are crucian carp trading off immune function against predator-induced morphology?J. Anim. Ecol.881510–1521.
105
VogelweithF.ThiéryD.MoretY.ColinE.MotreuilS.MoreauJ. (2014). Defense strategies used by two sympatric vineyard moth pests.J. Insect. Physiol.6454–61. 10.1016/j.jinsphys.2014.03.009
106
WellbornG. A.SkellyD. K.WernerE. E. (1996). Mechanisms creating community structure across a freshwater habitat gradient.Annu. Rev. Ecol. Syst.27337–363. 10.1146/annurev.ecolsys.27.1.337
107
WernerE. E.AnholtB. R. (1993). Ecological consequences of the trade-off between growth and mortality rates mediated by foraging activity.Am. Nat.142242–272. 10.1086/285537
108
WernerE. E.McPeekM. A. (1994). Direct and indirect effects of predators on two anuran species along an environmental gradient.Ecology751368–1382. 10.2307/1937461
109
ZuurA. F.LenoE. N.ElphickC. S. (2010). A protocol for data exploration to avoid common statistical problems.Methods Ecol. Evol.13–14. 10.1111/j.2041-210X.2009.00001.x
Summary
Keywords
Coenagrion, defensive traits, phenotypic plasticity, predation, trait correlations
Citation
Mühlenhaupt M, Jiang B, Brauner O and Mikolajewski DJ (2022) Inter- and Intraspecific Trait Compensation of Behavioural and Morphological Defences in a Damselfly Genus. Front. Ecol. Evol. 10:874276. doi: 10.3389/fevo.2022.874276
Received
11 February 2022
Accepted
04 April 2022
Published
26 April 2022
Volume
10 - 2022
Edited by
Mauricio J. Carter, Universidad Andres Bello, Chile
Reviewed by
Andreas Walzer, University of Natural Resources and Life Sciences Vienna, Austria; Mark Galatowitsch, Centre College, United States
Updates
Copyright
© 2022 Mühlenhaupt, Jiang, Brauner and Mikolajewski.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Max Mühlenhaupt, max.muehlenhaupt@uni-bielefeld.de
This article was submitted to Ecophysiology, a section of the journal Frontiers in Ecology and Evolution
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.