Aquatic Botany 128 (2016) 48–57

Contents lists available at ScienceDirect

Aquatic Botany

journal homepage: www.elsevier.com/locate/aquabot

Differences in herbivory intensity between the seagrass Cymodocea

nodosa and the green alga prolifera inhabiting the same habitat

Lucia Del Río, Javier Vidal, Séfora Betancor, Fernando Tuya

¨ IU-ECOAQUA¨, Grupo en Biodivesidad y Conservación, Marine Sciences Faculty, Universidad de Las Palmas de Gran Canaria, Las Palmas, Canary Islands,

Spain

a r t i c l e i n f o a b s t r a c t

Article history: Seagrasses are frequently found mixed with other macrophytes, e.g., green macroalgae. We aimed to

Received 8 March 2015

assess whether the magnitude of herbivory differed between two coexisting macrophytes, the seagrass

Received in revised form

Cymodocea nodosa (Ucria) Ascherson and the green seaweed Caulerpa prolifera (Forsskål) Lamouroux, at

25 September 2015

Gran Canaria Island (eastern Atlantic). Both in situ (field) and aquaria experimentation demonstrated

Accepted 2 October 2015

a larger intensity of herbivory (between ca. 4–8 times) on C. prolifera than C. nodosa. At the scale of

Available online 8 October 2015

meadows, herbivorous fish abundance predicted the intensity of herbivory, in particular by the parrotfish

Sparisoma cretense. A plant physical attribute (“force-to-fracture”) negatively correlated with a larger

Keywords:

Seagrass consumption on C. prolifera, while differences in total phenolic compounds between both macrophytes

Herbivory were insignificant. Importantly, herbivory marks (bites) were significantly larger (ca. two times) on C.

Grazing nodosa leaves than in C. prolifera fronds, so differences in the magnitude of herbivory between C. nodosa

Grazers and C. prolifera were dependent on herbivorous size.

Green alga © 2015 Elsevier B.V. All rights reserved.

Leaf traits

1. Introduction epiphytes on seagrass leaves and the seagrass itself (Wressnig and

Booth, 2007), as well as flowers and seeds (Balestri and Cinelli,

Many organisms are associated with seagrasses, below the sed- 2003). Epiphytes are a key element in the relationship between

iment linked to the rhizomes, upon the leaves and stems, and over herbivores and seagrasses. Some studies suggest that epiphytic pro-

the seagrass canopy (epi- and suprabenthic organisms, Herrera duction may be elevated to exceed even that of seagrasses (Morgan

et al., 2014; Tuya et al., 2014b), which move throughout the and Kitting, 1984; Chiu et al., 2013). Temperate and sub-tropical

meadow and constitute the main consumers of seagrass and asso- fishes consuming seagrass material select seagrass leaves and parts

ciated vegetated material. Seagrasses are evolutionarily adapted of leaves with abundant epiphytic loads, whereas mesograzers usu-

to herbivory; this is indicated by the range of evolutionary adap- ally feed on algae attached to seagrass leaves (Valentine and Heck,

tions to mitigate the consequences of herbivory (Karban and Myers, 1999; Goecker et al., 2005). It has been postulated that internal

1989), including mechanical and chemical elements (Lucas et al., contents in N are a relevant factor mediating feeding preferences by

2000). Traditionally, it has been postulated that a small fraction marine herbivores (Vergés et al., 2007; Prado et al., 2010; Prado and

of seagrass production is directly consumed by marine herbivores Heck, 2011), despite in other occasions its influence is negligible

(Cebrián and Duarte, 1998; Valentine and Heck, 1999). The low con- (Lee et al., 2015).

sumption of seagrass by grazers has been explained by their poor is a seagrass distributed across the entire

nutritional quality (Prado and Heck, 2011), including a high content Mediterranean and the adjacent Atlantic coasts, from the south-

in cellulose that act as a structural deterrent. Recent studies; how- ern Iberian Peninsula to Senegal, including Madeira and the Canary

ever, have pointed out that herbivory over seagrasses has a larger Islands (Cunha and Araujo, 2009). Meadows constituted by C.

influence than previously considered (Tomas et al., 2005; Heck and nodosa are found on shallow soft substrates of Gran Canaria Island

Valentine, 2006; Doropoulos et al., 2009; Prado et al., 2007; Vergés (Tuya et al., 2014a), where it may form mixed meadows with

et al., 2011). Seagrasses offer herbivores two main food sources: green rhizophytic seaweeds of the genera Caulerpa, particularly

Caulerpa prolifera (Fig. A1). As a result of environmental deteriora-

tion, frondose C. nodosa meadows can turn into bottoms dominated

∗ by C. prolifera; this has been reported from the Mediterranean

Corresponding author. Fax: +34 928452900.

and the southern Iberian Peninsula (Ceccherelli and Cinelli, 1997;

E-mail address: [email protected] (F. Tuya).

http://dx.doi.org/10.1016/j.aquabot.2015.10.001

0304-3770/© 2015 Elsevier B.V. All rights reserved.

L. Del Río et al. / Aquatic Botany 128 (2016) 48–57 49

Lloret et al., 2005), as well as from the Canary Islands (Tuya et al., and fronds of C. prolifera. Firstly, we conducted an indirect approach

2013b). C. nodosa may be an important food source for macro- by estimating herbivory pressure as the number of bite marks left

herbivores (Cebrián et al., 1996a). In addition, leaves of this seagrass by herbivores on both C. nodosa leaves and C. prolifera fronds; since

are colonized by epiphytic assemblages that may provide food for the majority of fishes inhabiting these seagrass systems are small-

associated invertebrates (Vizzini et al., 2002; Tuya et al., 2013a). sized (Espino et al., 2011), complete removal of seagrass leaves is an

Accompanying macrophytes, e.g., green seaweeds, can also rep- unlikely process. The study was carried out at two times: October

resent an additional food source for herbivores. However, certain 2013 and May 2014 to test for the effect of seasonality on responses.

macrophytes have developed several mechanisms to minimize her- At each of the 4 meadows, 12 leaves of C. nodosa and 12 fronds of

bivory (Duffy and Hay, 1990). The primary deterrent substances in C. prolifera were haphazardly collected by SCUBA divers; adjacent

seagrasses and seaweeds are phenolic compounds (Arnold et al., leaves/fronds were >2 m apart. Samples were quickly transported

2012), which have been linked to a variety of functions, prevent- to the laboratory and preserved in ice until analysis. At the same

ing bacterial infections (Harrison and Chan, 1980), protecting algae time of collection, fish assemblages were counted at daylight hours

from high PAR and UV damage (Pavia et al., 1997) and deterring (between 9:00 and 12:00 a.m.) through underwater visual cen-

2

grazers (Van Alstyne and Paul, 1990). Yet, there is some contro- suses, following (n = 4) 25 × 4 m transects per meadow (100 m of

versy in the real effectiveness of phenols as grazer deterrents (Close observation per census); the abundance and size of each fish

and McArthur, 2002; Vergés et al., 2007). For algae within the gen- was annotated according to standard procedures implemented in

era Caulerpa, it has been largely hypothesized that the presence the study region for seagrass meadows (Tuya et al., 2006). No

of repulsive (toxic) secondary metabolites, e.g., caulerpenyne, may major herbivorous invertebrates (e.g., sea-urchins) were found in

also deter herbivores (Box et al., 2010). Preference for vegetated the study area. Once in laboratory, we measured the length of C.

material among herbivores is; however, not exclusively related to nodosa leaves (from the ligule to the upper tip of each leaf) and

chemical attributes (Hay and Kappel, 1994), but also to the physical C. prolifera fronds (from the base of the stipe to the upper tip of

structure and configuration of macrophytes, e.g., their resistance to the frond). Fronds with proliferations were not considered to avoid

breakage (Duffy and Hay, 1990; Lucas et al., 2000; Prado and Heck, confusion. Bite marks were recorded for each leaf/frond through

2011). image analysis (imageJ freeware); all material was then preserved

In mixed meadows (i.e., those constituted by seagrasses and in silica gel. Some bite marks were clearly crescent-shaped (Fig. A2),

green seaweeds), macro-herbivores have several choices of food, a clear indication of consumption by herbivorous fishes (Hay, 1984;

what may generate different patterns of vegetation consumption. Kirsch et al., 2002; White et al., 2011; Lee et al., 2015). In these cases,

The aim of this work was to compare the magnitude of herbivory we recorded each bite size, as the maximum diameter of the mark

between the seagrass C. nodosa and the green alga C. prolifera; these (cm). The cover of epiphytic material was also annotated by using a

two macrophytes inhabit the same habitat (mixed meadows on qualitative, visual, scale: 0 (cover: <1%), 1 (cover: 1–10%), 2 (cover:

shallow subtidal waters) at Gran Canaria Island (eastern Atlantic). 10–20%), 3 (cover: 20–40%), 4 (cover: 40–60%) and 5 (cover: >60%).

Differences in the intensity of herbivory were compared by com- When the apical part of either seagrass leaves or C. prolifera fronds

bining in situ assays, that assessed indirect (bite marks) and direct was damaged, we omitted to record these as bite marks, due to the

(rates of consumption of fresh material) measures of herbivory, difficult of ascertaining if these marks resulted from herbivory or

and an aquaria experiment that quantified rates of consumption other type of damage (e.g., currents and/or swells).

on fresh material under controlled laboratory conditions. We set We took measurements of the “force-to-fracture” (FTF), as a way

out these procedures to specifically test whether the intensity of to assess the physical resistance to breakage of both macrophytes;

herbivory differed between C. nodosa and C. prolifera. We addition- these measurements were calculated with a dynamometer. The tip

ally hypothesized that spatial and temporal variation in herbivory of each of n = 30 leaves and fronds of both C. nodosa and C. prolifera

intensity on these two macrophytes is connected with differences was attached to the pin of the dynamometer; the force (Newtons)

in the abundances of herbivorous fish. Finally, we analyzed differ- necessary to tear each leaf/frond was then annotated. All leaves

ences in phenolic compounds concentration (a chemical attribute) and fronds were collected at Gando meadow (May 2014), encom-

and leaf/frond resistance to breakage (a physical attribute) between passing the entire range of available sizes; measurements were

both macrophytes to help to explain differences in herbivory pat- taken from fresh material immediately after collection. On the 17th

ters. December 2013, we randomly collected leaves of C. nodosa and

thalli of C. prolifera (ca. 0.25 g FW each thalli, n = 9) from Gando

meadow (8–10 m) to analyze differences in total phenolic com-

2. Materials and methods pounds. All material was stored at −80 C until analysis. Once in the

laboratory, all thalli were initially cleaned and epiphytes removed.

2.1. Field observational approach In all cases, we selected the central parts of the thalli, with no

evidence of grazing activity. All samples were grounded with a

Four study sites were selected in mixed meadows constituted mortar and a pestle in sand at 4 C, and extracted overnight in cen-

by the seagrass C. nodosa and the rhizophytic seaweed C. prolif- trifuge tubes with 2.5 ml of 80% (v/v) methanol (Betancor et al.,

era at the east coast of Gran Canaria Island (Fig. 1); depth ranged 2014). The mixture was centrifuged at 4000 rpm for 30 min and the

between 8 and 12 m, all bottoms were sandy and proximity from supernatants were collected (Sigma 2-16PK, Göttingen, Germany).

−1

the adjacent coast varied between 150 and 250 m. These mixed Total phenolic compounds, expressed as mg GAE g DW (Gallic

meadows are permanent all year round (Tuya et al., 2013b, 2014b); Acid Equivalent), were determined using gallic acid as a standard

−2 ◦

the biomass of C. nodosa varies between 120 and 170 g DW m (Folin and Ciocalteu, 1927) after 120 min in darkness at 4 C. The

−2

and the biomass of C. prolifera between 0 and 70 g DW m across absorbance was then measured at 760 nm in a spectrophotometer

sites 10 s of meters apart (Tuya et al., 2013b). Alternative vegetation (Thermo Scientific Evolution 201, UV-visible, China).

is sparse and mainly restricted to epiphytes growing on seagrass

leaves. The seagrass shows a clear seasonal pattern, including a 2.2. Field experimentation

maximum in shoot density and biomass in summer and a mini-

mum in winter (Tuya et al., 2006). No information is available on Through an in situ experiment, we offered herbivores fresh C.

seasonal patterns of C. prolifera. We developed different types of nodosa leaves and C. prolifera fronds alternatively attached with

assays to evaluate the magnitude of herbivory on C. nodosa leaves clothespins to plastic mesh frames (Fig. A3). This was a way to

50 L. Del Río et al. / Aquatic Botany 128 (2016) 48–57

Fig. 1. Location of study area; including the 4 sampled meadows: C (Caballo), G (Gando), RA (Roque de Arinaga) and RV (Risco Verde).

directly estimate differences in consumption over both C. nodosa photoperiod. Only seagrass leaves/algal fronds without grazing

and C. prolifera by local herbivores (a double choice experiment). scars were selected. Each aquaria had two compartments sepa-

Fresh material was initially collected by SCUBA divers, before the rated by a mesh net (Fig. A4); one for the feeding trial including

experimental set up, which was located at Gando meadow at 10 m the fish, while the other compartment lacked fish and so acted as

depth (Fig. 1). A total of 14 plots (mesh frames) were directly a control to correct for possible autogenic changes in biomass not

attached to the bottom with metal bars and then covered with sand directly caused by grazing (Wressnig and Booth, 2007; Prado et al.,

to minimize any visual effect. Plots were placed at two randomly- 2011). A total of 5 seagrass leaves (with and without epiphytes)

selected areas within this meadow. The distance between adjacent and 5 fronds of C. prolifera were placed into each experimental

plots was approximately 0.5 m. All leaves/fronds were of the same compartment per aquaria. All vegetated material was distributed

length (15 cm for C. nodosa and 7 cm for C. prolifera fronds); we only randomly and secured to the bottom with wire stakes (Fig. A4,

used leaves with no previous herbivore marks and with no signs for a video clip visit the YouTube link: http://www.youtube.com/

of necrosis. Furthermore, to offer herbivores material under simi- watch?v=IhWCk62SZkE&list=UU9zh1SloG4bMsoJbW8y-hmw).

lar conditions, epiphytes were scraped off the leaves using a thin Prior to experimentation, the wet weigh of all vegetated material

blade. Three leaves of C. nodosa and 3 fronds of C. prolifera were per aquaria was obtained. We then conducted the feeding trial

used per plot. After 7 days underwater, all material was retrieved. over a 24 h period. After that time, the remaining material was

Once in the laboratory, the number of bites marks per leaf/frond removed from each aquaria, blot-dried, and weighed, what pro-

−1 −1

was recorded, as in the previous assay. Despite the number of bite vided consumption rates (g wet weight fish d ). The numbers

marks is not as critical as the mass of removed vegetation, it pro- of bite marks per leaf/frond was then annotated. To work out

vides an unconfounded measure of herbivory intensity; e.g., it is whether parrotfish had a preference for a particular macrophyte,

plausible that removal of vegetation biomass results, for example, we used Chesson’s index (Chesson, 1983):

from the action of waves. This assay was carried out at three dif-

ri/pi

  ferent times: October 2013, May 2014 and June 2014 to test for the ˛ =

r /p

temporal consistency of results at random times. i i

where ri is the proportional biomass removal due to consumption

on macrophyte i, and p is the proportional biomass of macrophyte

2.3. Aquaria experimentation i

i in the source sample (the aquaria). The value of ˛ varies between

0 and 1. When ˛ = 1/n, selective feeding does not happen (n is the

We tested for differences in consumption by the main herbivo-

total number of macrophytes). There were 3 macrophytes (n = 3,

rous fish in seagrass meadows of the study region, the parrotfish

1/n = 0.33). If ˛ > 1/n, it is then likely that parrotfish has a preference

Sparisoma cretense (Tuya et al., 2013a; this study: Appendix C),

for macrophyte i, whereas if ˛ < 1/n, the macrophyte i is probably

between C. nodosa seagrass leaves with, or without, epiphytes

avoided. Values of ˛ were calculated for each of the 3 macrophytes

and fronds of C. prolifera, under controlled laboratory conditions.

separately for every aquaria.

Fresh material of these 3 choices was offered to 3 S. cretense

per aquaria. Fish were collected in Gando meadow through a

seine net (May 2014) and immediately transported to the aquaria 2.4. Statistical analyses

facility under constant aeration. Fish were acclimatized to indoor

conditions for 7 days, with a lack of food provision prior to exper- Differences in the number of bite marks per thallus between the

imentation (Prado et al., 2011; Tuya and Duarte, 2012). All fish two macrophytes was tested using a 3-way ANOVA, which incorpo-

were small-sized (juveniles, total length between 7 and 10 cm). rated the factors: “Time” (fixed factor with 2 levels: October 2013

This experiment was replicated in 4 aquaria (84 l) with constant and May 2014), “Meadow” (random factor with 4 levels: the 4 study

aeration. The bottom of all aquaria was covered by sand; sea-water sites) and “Species” (fixed factor with two levels: C. nodosa vs. C.

temperature varied between 21 and 23 C under a 12:12 h natural prolifera). “Time” was considered as a fixed factor, since it reflected

L. Del Río et al. / Aquatic Botany 128 (2016) 48–57 51

Fig. 2. Herbivory intensity over Cymodocea nodosa leaves and Caulerpa prolifera fronds (mean number of grazing marks per thallus, n = 12 leaves/fronds) in (a) October 2013

and (b) May 2014. Error bars are +SE of means.

Table 1

Results of 3-way ANOVA testing for differences in the total number of grazing marks between times, meadows and macrophytes (C.n. = Cymodocea nodosa and C.p. = Caulerpa

prolifera).

df MS F-ratio P

Time 1 18.49 12.04 0.0404

Meadow 3 1.35 22.36 0.0001

C.n. vs. C.p. 1 52.9 38.11 0.0086

Time × meadow 3 1.5367 25.28 0.0001

Time × C.n. vs. C.p. 1 12.28 6.77 0.0803

×

Meadow C.n. vs. C.p. 3 1.3883 22.84 0.0001

Time × meadow × Cn. vs. C.p. 3 1.8161 29.87 0.0001

Residual 176 0.0608

potential seasonal differences. Data were square-root transformed significantly larger number of bite marks was recorded on C. prolif-

prior to analyses to avoid heterogeneous variances, as indicated era relative to C. nodosa. The overall number of bites was larger

by the Cochran’s test. Differences in the number of bite marks per in October (2013) than on May (2014) (Fig. 2; Table 1, “Time”

thallus after 1 week between both macrophytes in the in situ exper- P = 0.04); the magnitude of differences in the mean number of

iment was tested by a 3-way ANOVA, which included a similar bites per thallus varied from October (2013) to May (2014) (Fig. 2;

design relative to the previous analysis, but the factor “Meadow” Table 1, “Time × Cn. vs. C.p.” P = 0.083). In general, the total abun-

was replaced by “Plot” (also a random factor). A t-test tested for dance of herbivorous fishes (Appendix C) and the abundance of

significant differences in FTF between C. nodosa and C. prolifera at the most conspicuous herbivorous fish, the parrotfish S. cretense,

the range of sizes they typically occur, as well as the significance significantly predicted differences in herbivory intensity between

of differences in the size of grazing marks (bites) between both meadows (Fig. 3). In the field experiment, the number of bite marks

macrophytes. We ignored the contribution of spatial and temporal on C. prolifera fronds after 2 weeks exceeded those on seagrass

variation in differences in bite size between both macrophytes as leaves (Fig. 4; Table 2, “C.n. vs. C.p.” P = 0.003), particularly on time

a result of the large imbalances in the number of bite marks per 3 (June, 2014) (Fig. 4; Table 2, “Time × C.n. vs. C.p.” P = 0.008).

time and meadow; this was because only clear crescent-shaped Overall, C. nodosa leaves had a larger resistance to breakage

bite marks a doubtless indication of consumption by herbivorous (FTF) than fronds of C. prolifera, when considering the range of

fishes—were considered in the analysis. Differences in consump- sizes that both macrophytes typically reach in situ (Fig. 5; t-

tion in the aquaria trial were tested by means of a 2-way ANOVA, student = 1.9047, P = 0.061). The size of grazing marks (bites) on

which included the factors: “Treatment” (fixed factor with 3 levels: C. nodosa leaves were larger than those on C. prolifera fronds

C. nodosa with epiphytes, C. nodosa without epiphytes and C. prolif- (C. nodosa: 0.4413 cm ± 0.0238; C. prolifera: 0.2850 cm ± 0.041,

era) and “Aquaria” (random factor with 4 levels: aquaria 1, aquaria mean ± SE, n = 20; t-test = 6.0134, P < 0.00001).

2, aquaria 3 and aquaria 4). When appropriate, pair-wise compar- In the aquaria assay, consumption of C. prolifera fronds was

isons (via SNK tests) were performed. Significant differences in the larger than consumption of C. nodosa leaves either with or with-

content of phenols between both macrophytes were tested through out epiphytes (Fig. 6; Table 3, “Treatment”, P = 0.03), despite the

a t-test. A one-way ANOVA tested whether Chesson’s ˛ values dif- magnitude consumption among treatments varied from aquaria

fered between the 3 macrophytes. The significance level was set at to aquaria (Fig. 6; Table 3, “Aquaria × Treatment”, P = 0.01); in

the type I error rate of 0.05. all cases, independent SNKs tests for each aquaria demonstrated

a larger consumption on C. prolifera relative to C. nodosa with

or without epiphytes. In 3 out of the 4 aquaria, consump-

3. Results

tion on C. nodosa with epiphytes was larger than on C. nodosa

without epiphytes (SNK tests). Parrotfish clearly preferred C. pro-

The mean total number of bites per thallus differed between C.

lifera (˛ = 0.72 ± 0.04) over C. nodosa with (˛ = 0.16 ± 0.02) and

nodosa and C. prolifera (Fig. 2; Table 1, “C.n. vs. C.p.”, P = 0.0086); a

52 L. Del Río et al. / Aquatic Botany 128 (2016) 48–57

Fig. 3. Relationship between the mean number of bites on Cymodocea nodosa leaves (a and b) and Caulerpa prolifera fronds (c and d) and the abundance of the parrotfish,

Sparisoma cretense, and the total herbivorous fish abundance per meadow. Data were pooled for the two sampling times. Unfilled symbols: October 2013, filled symbols:

May 2014.

C. prolifera in Gran Canaria Island. The seaweed C. prolifera con-

sistently suffered higher rates of herbivory than C. nodosa, both

in the field, through direct and indirect measures, and in aquaria

conditions. Average feeding rates by S. cretense were estimated at

−1 −1

0.209 g wet weight ind d for C. prolifera and 0.086 g wet weight

−1 −1

ind d for C. nodosa, under laboratory conditions. Despite com-

parisons with other studies are complicated as a result of varying

methods and timing, the average feeding rates on C. nodosa leaves

−1 −1

obtained by our study (ca. 7.5 g seagrass kg fish d ) is within the

same order of magnitude as those reported for C. nodosa at other

zones. For example, Goldenberg and Erzini (2014) estimated the

average feeding rate by adults of the Sparid Sarpa salpa from the Ria

−1 −1

Formosa lagoon as ca. 32 g kg fish d ; obviously, this larger con-

sumption is explained because adult fish consume larger amounts

of vegetated material than juveniles.

Herbivorous fish abundance predicted the intensity of herbivory

over C. nodosa leaves and C. prolifera fronds through a positive cor-

relation between the average number of bite marks per thallus and

Fig. 4. Mean number of bites over Cymodocea nodosa leaves and Caulerpa prolifera

total fish abundance. This result is similar to patterns of seagrass

fronds after 1 week at each replicated assay (n = 14, T1 = October 2013, T2 = May 2014

consumption in tropical environments, where the abundances of

and T3 = June 2014). Error bars are +SE of means.

parrotfish is the main driver of removal of seagrass material (Lee

et al., 2015), despite in other occasions the abundance of megafauna

without epiphytes (˛ = 0.10 ± 0.02) (1-way ANOVA, F2,9 = 428.18,

(e.g., turtles) are the main contributors (Christianize et al., 2014).

P = 0.00001).

Herbivory intensity seemed to be seasonal, reaching maximum val-

The total content of phenols did not differ between both macro-

−1 ues in late summer (October 2013) and minimum values in spring

phytes (C. nodosa: 1.780 mg GAE g DW ± 0.1872; C. prolifera:

−1 (May 2014). There is a clear correlation between the annual vitality

1.8157 mg GAE g DW ± 0.1524; t-test = 0.14, P = 0.5737).

cycle of C. nodosa and the richness and abundance of associated

fish assemblages in the study region (Tuya et al., 2006; Espino

4. Discussion

et al., 2011), with maximum values in summer–autumn and mini-

mum in winter–spring. In particular, the abundance of small-sized

Our study has demonstrated that herbivory may remove

S. cretense, the main herbivorous fish in the study region, reach

amounts of vegetated material in mixed meadows of C. nodosa and

L. Del Río et al. / Aquatic Botany 128 (2016) 48–57 53

Table 2

Results of 3-way ANOVA testing for differences in the total number of grazing marks per thallus between times, plots and macrophytes (C.n. = Cymodocea nodosa and

C.p. = Caulerpa prolifera).

df MS F-ratio P

Time 1 1.32 8.217 0.015

Plot 13 0.4177 0.965 0.5

C.n. vs. C.p. 1 9.21 22.77 0.003

Time × plot 13 0.16 0.37 0.977

Time × C.n. vs. C.p. 1 1.7 11.46 0.008

Plot × C.n. vs. C.p. 13 0.4 0.93 0.536

×

Time plot × C.n. vs. C.p. 13 0.15 0.34 0.987

Residual 112 0.43

14 a) b) 12

10

8

6 Force (N)

4

2

0

0 5 10 15 20 25 30 35 2 4 6 8 10 12 14

Leaf length (cm) Frond len gth (cm)

Fig. 5. Differences in the “force-to-fracture” between (a) Cymodocea nodosa leaves and (b) Caulerpa prolifera fronds.

Table 3

Results of 3-way ANOVA testing for differences in mass loss between aquariums and treatments (Cymodocea nodosa leaves with epiphytes, C. nodosa leaves without epiphytes

and Caulerpa prolifera fronds).

df MS F-ratio P

Aquaria 3 23.34 17.29 0.0010

Treatment 2 27.14 6.07 0.0380

Aquaria × treatment 6 4.47 3.31 0.0100

Residual 48 1.35

maximum abundances in September–October (Espino et al., 2014). bite marks on seagrass blades (Fig. A5). Typically, adult and sub-

Our results are, moreover, in accordance with seasonal trends for adult fishes are big enough to move across meadows without the

herbivory in seagrass meadows from the Mediterranean, where necessity to hide from predators and so the entire canopy is avail-

maximum consumption rates occurred in summer (Tomas et al., able to them, while juveniles hide from predators within the dense

2005; Prado et al., 2007, 2010). Chiu et al. (2013) and Lee et al. canopy provided by seagrass leaves and accompanying green sea-

(2015) also demonstrated that leaf grazing rates were significantly weeds (Espino et al., 2014). As a result, small-sized fishes have a

greater in summer and autumn than in winter and spring. Sea- direct access to the shorter fronds of C. prolifera, which are softener

sonal differences in herbivory pressure can be related to migratory to feed. This hypothesis is consistent with the assumption that adult

(ontogenetic) movements, or varying seasonal feeding behavior, of fishes inhabiting seagrass meadows have a different trophic niche

herbivores (Prado et al., 2007). in comparison with juveniles (Livingston, 1982; Vizzini et al., 2002).

The different intensity of herbivory between C. nodosa and C. Many herbivores find difficulties to consume plant material

prolifera was herbivorous size-dependent, as bite marks were sig- (Watson and Norton, 1985). Calcification and toughness usually

nificantly larger on C. nodosa leaves than in C. prolifera fronds. correlate with low feeding preference (Litter et al., 1983; Hay, 1984;

The mean bite size on C. nodosa leaves normally exceeded 0.4 mm, Paul and Hay, 1986). Despite we lack data for our case study, C.

which likely corresponds with sub-adult fishes (e.g., S. cretense). In prolifera typically shows a higher total internal N content than C.

contrast, mean bite sizes on C. prolifera fronds were <0.3 mm, which nodosa and hence a lower C:N ratio (García-Sánchez et al., 2012).

likely correspond with juvenile fish stages. In other words, small- Plant quality (often typified as C:N ratios) has been largely demon-

sized herbivorous fishes preferably consume C. prolifera; once a strated to affect patterns in feeding preferences by herbivores in

certain size is reached, herbivorous fish may consume both macro- coastal habitats (Cebrián and Duarte, 1998), despite C:N rations

phytes. In this sense, it is plausible that the jaws of S. cretense may not explain feeding preferences over a suite of coexisting sea-

juveniles do not have enough force to ripe off pieces of C. nodosa grass species (Lee et al., 2015). In this sense, C. nodosa is a seagrass

leaves; this idea is consistent with the presence of unsuccessful with high fiber content (large amounts of cellulose in their cell

54 L. Del Río et al. / Aquatic Botany 128 (2016) 48–57

0.35 (McConnell et al., 1982) and a few reef fish (Targett et al., 1986;

C. nodosa +

Epiphytes Paul et al., 1987) have been deterred by caulerpenyne, despite most

C. nodosa - Epiphytes

0.30 reef fish were not deterred by Caulerpa extracts containing cauler-

Cau lerpa prolifera

penyne (Paul et al., 1987, 1990; Wylie and Paul, 1988). Macroalgal

fronds have often higher concentrations of caulerpenyne than

0.25

stolons, and this concentration may even change during the year

te

(Box et al., 2010) within and among species. In any case, the clas- Ra

0.20

sic idea that Caulerpa species deter herbivores seems not to be ion

t

supported by our data, and fits previous observations reporting

0.15 consumption of Caulerpa by a range of herbivores (Gavagnin et al.,

1994). The convergence of results from indirect and direct in situ

Consump 0.10 measurements and the aquaria assay reinforce this idea. With these

results, we hypothesize that caulerpenyne may have a higher effect

on small organisms (like epibiota or epifauna) than in vertebrates.

0.05

It is also possible that macro-herbivores inhabiting mixed seagrass

meadows may be adapted to allelochemicals produced in the sys-

0.00 tem (Prado et al., 2011; Goldenberg and Erzini, 2014). Seagrass

Aq uaria 1 Aquaria 2 Aquaria 3 Aquaria 4

leaves, C. nodosa in particular, are extensively colonized by a com-

−1 −1 plex epiphytic community (cyanobacteria, diatoms, crustose and

Fig. 6. Consumption rates of vegetated material (g wet weight ind d ) by the

ephemeral algae, invertebrates etc.), which provide food and habi-

parrotfish Sparisoma cretense (n = 5). Error bars are +SE of means.

tat for invertebrates and so increases the spatial complexity of

the habitat (Mazzella et al., 1992; Vizzini et al., 2002). Our results

walls). A high fraction of its internal C is used for the synthesis showed that C. nodosa leaves with epiphytes are preferred by her-

of structural carbohydrates that forms cell walls and fiber bun- bivores, at least by the parrotfish S. cretense in indoor conditions,

dles, which means that C. nodosa has a high mechanical resistance, over leaves devoid of epiphytes. This is consistent with previous

thereby exhibiting a low leaf nutritional value (Lucas et al., 2000; De findings, which confirm that fishes grazing on seagrass prefer epi-

los Santos et al., 2012). Fast-growing species (such as C. prolifera); phytes growing on seagrass leaves (Cebrián et al., 1996b; Wressnig

however, do not invest resources in leaf/frond toughening as much and Booth, 2007).

as large, long-lived, seagrass species (De los Santos et al., 2012). As In summary, this study has demonstrated that (i) C. prolifera

indicated by our data (through the “force-to-fracture” calculation), is more consumed by herbivores, particularly fish, than C. nodosa

C. nodosa is tougher than C. prolifera. This trait; however, was only when found intermixed; (ii) variation in herbivory intensity is con-

evaluated at one meadow and time, so we cannot totally infer our nected with variability in the abundance of herbivorous fishes; (iii)

results to all conditions. differences in the magnitude of herbivory between C. nodosa and

The polyphenols of seagrasses and seaweeds assist algae to C. prolifera were dependent on herbivorous size.

overcome oxidative stress (i.e., while subjected to high PAR and

UV levels), playing a putative role in their defense against herbi- Acknowledgements

vores (Altena and Steinberg, 1992), despite significant controversy

We acknowledge T. Sánchez and F. Espino for their help dur-

in their efficiency (Vergés et al., 2007). The presence of phenols in C.

ing fieldwork, laboratory assays and always useful comments, and

nodosa (García-Sánchez et al., 2012) and other seagrasses has been

R. Triay for his guidance during the aquaria assays. F. Tuya was

reported (Vergés et al., 2011). In our study, the total phenolic con-

supported by the MINECO ‘Ramón y Cajal’ program. Part of this

tents was similar for both macrophytes. This result was somehow

study was performed through the project ECOSERVEG (BEST ini-

unexpected, because a higher phenol content has been indicated for

tiative, Voluntary Scheme for Biodiversity and Ecosystem Services

C. prolifera than C. nodosa from the Mediterranean (García-Sánchez

in Territories of the EU Outermost Regions and Oversees Coun-

et al., 2012). This disagreement could result from different sam-

tries and Territories, Grant no 07.032700/2012/635752/SUB/B2).

pling depths between studies (1–2 meters in García-Sánchez et al.

The research staff was partially supported by the Campus Atlántico

(2012) vs. 8–10 m in our study); depth is known to affect the con-

Tricontinental.

tent of phenolic compounds, with higher levels for shallow-water

macrophytes (Betancor et al., 2015). Moreover, we lack temporal

Appendix A.

replication and so it is plausible that phenols concentrations may

change through seasons. As a result, the putative role of phenols

mediating differences in consumption between C. prolifera than C.

nodosa remains unresolved.

Siphonous, unicellular, green algae, e.g., Caulerpa spp., contain

secondary metabolites, which play a role during quick plug forma-

tion after injuries to avoid the release of internal cellular fluids.

The production and utilization of secondary metabolites has been

indicated to justify that plant chemistry is the central factor deter-

mining herbivore feeding choices (Ehrlich et al., 1988; Schultz,

1988). However, the relative importance of these factors proba-

bly varies from system to system depending on the identity of

herbivores, their abundances and sizes (Pennings and Paul, 1992;

Lee et al., 2015). Numerous studies reported that caulerpenyne,

the main secondary metabolite within the genus Caulerpa, actively

deter herbivores (e.g., gastropods and fishes), which is toxic to

larval and adult stages of several marine invertebrates and ver-

tebrates (Paul et al., 1987; Pedrotti and Lemée, 1999). Sea urchins Fig. A1. Mixed meadow of Cymodocea nodosa and Caulerpa prolifera.

L. Del Río et al. / Aquatic Botany 128 (2016) 48–57 55

Fig. A5. “Unsatisfactory” bite on a Cymodocea nodosa leaf.

Fig. A2. Crescent-shaped bite mark on a Cymodocea nodosa leaf.

Appendix B.

Average percentage of daily changes in Cymodocea nodosa leaves

and Caulerpa prolifera fronds as a result of autogenic processes.

Cymodocea nodosa Cymodocea nodosa Caulerpa prolifera

with epiphytes without epiphytes

Aquaria 1 0.377 2.432 0.069

Aquaria 2 1.877 0.860 1.848

Aquaria 3 0.631 1.263 1.322

Aquaria 4 0.068 −0.886 2.387

Total 2.199 3.670 5.628

Appendix C.

−2

Fish abundances at each meadow (100 m ± SE); species trophic

status are also included: CAR: carnivorous, OMN: omnivorous,

Fig. A3. Mesh frame with attached Cymodocea nodosa leaves and Caulerpa prolifera

fronds. The mesh and pegs were subsequently buried with sand. HER: herbivorous (*: only juveniles) at time 1 (October 2013) and

time 2 (May 2014).

Fig. A4. Aquaria with two separated compartments by a mesh (M), including a CON-

TROL without fish, and a FEEDING TRIAL with fish (F). The three treatments included:

C.n. +, Cymodocea nodosa leaves with epiphytes, C.n. −, C. nodosa leaves without

epiphytes, and C.p. Caulerpa prolifera fronds.

56 L. Del Río et al. / Aquatic Botany 128 (2016) 48–57

References

Altena, V., Steinberg, P.D., 1992. Are differences in the responses between North 0.85 3.97 1.03 1.03 0.25 0.25

American and Australasian marine herbivores to phlorotannins due to ± ± ± ± ± ±

1 0.71

differences in phlorotannin structure? Biochem. Syst. Ecol. 20, 493–499. ± ±

Arnold, T., Mealey, C., Leahey, H., Miller, A.W., Hall-Spencer, J.M., 2012. Ocean

T2 0.25 3.75 0.25 8.25 1 2.25 1.75 1

acidification and the loss of phenolic substances in marine plants. PLoS One 7,

e35107.

Balestri, E., Cinelli, F., 2003. Sexual reproductive success in Posidonia oceanica.

Aquat. Bot. 75, 21–32.

Betancor, S., Tuya, F., Gil-Díaz, T., Figueroa, F.L., Haroun, R., 2014. Effects of a 0.47 0.25 0.25

Verde submarine eruption on the performance of two brown seaweeds. J. Sea Res. 87, 0.86

25 ± ± ±

0.7 ±

68–78.

±

±

Betancor, S., Domínguez, B., Tuya, F., Figueroa, F.L., Haroun, R., 2015.

Risco T1 25 0.75 0.25 2.5 1 0.25

Photosynthetic performance and photoprotection of Cystoseira humilis

(Phaeophyceae) and Digenea simplex (Rhodophyceae) in an intertidal rock

pool. Aquat. Bot. 121, 16–25.

Box, A., Sureda, A., Tauler, P., Terrados, J., Marbà, N., Pons, A., Deudero, S., 2010.

Seasonality of caulerpenyne content in native Caulerpa prolifera and invasive C.

0.25 0.25

taxifolia and C. racemosa var. cylindracea in the western Mediterranean Sea. 1.26 0.87

± ±

1 Bot. Mar. 53, 367–375. ± ±

±

Cebrián, J., Duarte, C.M., Marbà, N., 1996a. Herbivory on the seagrass Cymodocea

T2 3.5 0.25 0.25 1.5 1

nodosa (Ucria) Ascherson in contrasting Spanish Mediterranean habitats. J.

Exp. Mar. Biol. Ecol. 204, 103–111.

Cebrián, J., Duarte, C.M., Marbà, N., Enríquez, S., Gallegos, M., Olesen, B., 1996b.

Herbivory on Posidonia oceanica: magnitude and variability in the Spanish

Mediterranean. Mar. Ecol. Prog. Ser. 130, 147–155.

Arinaga

17.5 0.47 0.25 0.25 0.25 Cebrián, J., Duarte, C.M., 1998. Patterns in leaf herbivory on seagrasses. Aquat. Bot.

0.28 4.41 0.65 0.5

5.7 ± ± ± ± ±

60, 67–82. ± ± ± ±

±

Ceccherelli, G., Cinelli, F., 1997. Short-term effects of nutrient enrichment of the

Roque T1 17.5 3.5 0.5 10 0.25 4.25 7.5 1.5 0.25 0.25

sediment and interactions between the seagrass Cymodocea nodosa and the

introduced green alga Caulerpa taxifolia in a Mediterranean bay. J. Exp. Mar.

Biol. Ecol. 217, 165–177.

Chesson, J., 1983. The estimation and analysis of preference and its relationship to

foraging models. Ecology 64, 1297–1304.

Chiu, S., Huang, Y., Lin, H., 2013. Carbon budget of leaves of the tropical intertidal

seagrass Thalassia hemprichii. Estuar. Coast. Shelf Sci. 125, 27–35.

0.25 0.25

Christianen, M.J.A., Herman, P.M.J., Bouma, T.J., Lamers, L.P.M., van Katwijk, M.M., ± ±

2.35 0.41

van der Heide, T., Mumby, P.J., Silliman, B.R., Engelhard, S.L., van de Kerk, M.,

± ±

Kiswara, W., van de Koppel, J., 2014. Habitat collapse due to overgrazing

T2 1.25 1 0.25

threatens turtle conservation in marine protected areas. Proc. R. Soc. B 281,

1777.

Close, D.C., McArthur, C., 2002. Rethinking the role of many plant

phenolics—protection from photodamage not herbivores? Oikos 99, 166–172.

Cunha, A.H., Araujo, A., 2009. New distribution limits of seagrass beds in West 0.33 0.16 0.16 0.54 3 0.16

Africa. J. Biogeogr. 36, 1621–1622. 0.22 0.22

± ± ± ± ±

± ±

De los Santos, C.B., Brun, F., Onoda, Y., Cambridge, M.L., Bouma, T.J., Vergara, J.J.,

Pérez-Lloréns, J.L., 2012. Leaf-fracture properties correlated with nutritional

Gando T1 0.67 1.5 0.16 0.16 0.83 0.5 0.16

traits in nine Australian seagrass species: implications for susceptibility to

herbivory. Mar. Ecol. Prog. Ser. 458, 89–102.

Doropoulos, C., Hyndes, G.A., Lavery, P.S., Tuya, F., 2009. Dietary preferences of two

seagrass inhabiting gastropods: allochthonous vs autochthonous resources.

Estuar. Coast. Shelf Sci. 83, 13–18.

0.65 0.86

Duffy, J.E., Hay, M.E., 1990. Seaweed adaptations to herbivory. Bioscience 40, 0.41 ± ±

368–375. ±

T2 1.5 3 1.5

Espino, F., Tuya, F., Brito, A., Haroun, R.J., 2011. Ichthyofauna associated with

Cymodocea nodosa meadows in the Canarian Archipelago (central eastern

Atlantic): community structure and nursery role. Cienc. Mar. 37, 157–174.

Espino, F., González, J.A., Haroun, R., Tuya, F., 2014. Abundance and biomass of the

parrotfish Sparisoma cretense in seagrass meadows: temporal and spatial

differences between seagrass interiors and seagrass adjacent to reefs. Environ.

0.75 0.63

Biol. Fish. 98, 121–133. ± ±

Folin, O., Ciocalteu, V., 1927. On tyrosine and tryptophane determinations in

Caballo T1 1.25 1.75 proteins. J. Biol. Chem. 12, 239–243.

Gavagnin, M., Marin, A., Castellucio, F., Villani, G., Cimino, G., 1994. Defensive

relationships between Caulerpa prolifera and its shelled sacoglossan predators.

J. Exp. Mar. Biol. Ecol. 175, 197–210.

García-Sánchez, M., Korbee, N., Pérez-Ruzafa, I., Marcos, C., Domínguez, B.,

Figueroa, F.L., Pérez-Ruzafa, A., 2012. Physiological response and

photoacclimation capacity of Caulerpa prolifera (Forsskíl) J.V. Lamouroux and status

Cymodocea nodosa (Ucria) Ascherson meadows in the Mar Menor lagoon (SE

Spain). Mar. Environ. Res. 79, 37–47.

Goecker, M.E., Heck, K.L., Valentine, J.F., 2005. Effects of nitrogen concentrations in HER CAR

Trophic OMN CAR OMN CAR HER* CAR CAR OMN CAR OMN CAR CAR CAR CAR

turtlegrass Thalassia testudinum on consumption by the bucktooth parrotfish

Sparisoma radians. Mar. Ecol. Prog. Ser. 286, 239–248.

Goldenberg, S.U., Erzini, K., 2014. Seagrass feeding choices and digestive strategies

of the herbivorous fish Sarpa salpa. J. Fish. Biol. 84, 1474–1489.

Harrison, P.G., Chan, A.T., 1980. Inhibition of the growth of micro-algae and

bacteria by extracts of eelgrass (Zostera marina) leaves. Mar. Biol. 61, 21–26.

Hay, M.E., 1984. Predictable spatial escapes from herbivory: how do these affect

cantharus the evolution of herbivore resistance in tropical marine communities?

hispidus

maderensis

capistrata

marmoratus Oecologia 64, 396–407. pavo

viridensis cretense

draco

novacula Hay, M.E., Kappel, Q.E., 1994. Synergisms in plant defenses against herbivores:

annularis

pastinaca sp.

podas surmuletus auriga pagrus interactions of chemistry, calcification, and plant quality. Ecology 75,

boops 1714–1726. Species Boop Sphoeroides Synodus Mullus Diplodus Sparisoma Xyrichtys Bothus Dasyatis Spondyliosoma Stephanolepis Trachinus Canthigaster Pagrus Pagrus Sphyraena Thalassoma

L. Del Río et al. / Aquatic Botany 128 (2016) 48–57 57

Heck, K.L., Valentine, J.F., 2006. Plant-herbivore interactions in seagrass meadows. Targett, N.M., Targett, T.E., Vrolijk, N.H., Ogden, J.C., 1986. Effect of macrophytes

J. Exp. Mar. Biol. Ecol. 330, 420–436. secondary metabolites on feeding preferences of the herbivorous parrotfish

Herrera, A., Landeira, J.M., Tuya, F., Packard, T., Espino, F., Gómez, M., 2014. Sparisoma radians. Mar. Biol. 92, 141–148.

Seasonal variability of suprabenthic crustaceans associated with Cymodocea Tomas, F., Turón, X., Romero, J., 2005. Seasonal and small-scale spatial variability of

nodosa seagrass meadows off Gran Canaria (eastern Atlantic). Contint. Shelf herbivory pressure on the temperate seagrass Posidonia oceanica. Mar. Ecol.

Res. 88, 1–10. Prog. Ser. 301, 95–107.

Karban, R., Myers, J.H., 1989. Induced plant responses to herbivory. Annu. Rev. Ecol. Tuya, F., Martín, J.A., Luque, A., 2006. Seasonal cycle of a Cymodocea nodosa

Syst. 20, 331–348. seagrass meadow and the associated ichthyofauna at Playa Dorada (Lanzarote

Kirsch, K.D., Valentine, J.F., Heck, K.L., 2002. Parrotfish grazing on turtlegrass Canary Islands, eastern Atlantic). Cienc. Mar. 32, 695–704.

Thalassia testudinum: evidence for the importance of seagrass consumption in Tuya, F., Duarte, P., 2012. Role of food availability in the bathymetric distribution of

food web dynamics of the Florida Keys National Marine Sanctuary. Mar. Ecol. the starfish Marthasterias glacialis (Lamk.) on reefs of northern Portugal. Sci.

Prog. Ser. 227, 71–85. Mar. 76, 1–6.

Lee, C.L., Huang, Y.H., Chung, C.Y., Hsiao, S.C., Lin, H.J., 2015. Herbivory in Tuya, F., Viera-Rodríguez, M.A., Guedes, R., Espino, F., Haroun, R., Terrados, J., 2013a.

multi-species, tropical seagrass beds. Mar. Ecol. Prog. Ser. 525, 65–80. Seagrasses responses to nutrient enrichment depend on clonal integration, but

Livingston, R.J., 1982. Trophic organization of fishes in a coastal seagrass system. not flow-on effects on associated biota. Mar. Ecol. Prog. Ser. 490, 23–35.

Mar. Ecol. Prog. Ser. 7, 1–12. Tuya, F., Hernandez-Zerpa, H., Espino, F., Haroun, R., 2013b. Drastic decadal decline

Lloret, J., Marín, A., Marín-Guirao, L., Velasco, J., 2005. Changes in macrophytes of the seagrass: Cymodocea nodosa at Gran Canaria (eastern Atlantic):

distribution in a hypersaline coastal lagoon associated with the development interactions with the green algae Caulerpa prolifera. Aquat. Bot. 105, 1–6.

of intensively irrigated agriculture. Ocean Coast. Manage. 48, 828–842. Tuya, F., Ribeiro-Leite, L., Arto-Cuesta, N., Coca, J., Haroun, R., Espino, F., 2014a.

Lucas, P.W., Turner, I.M., Dominy, N.J., Yamashita, N., 2000. Mechanical defenses to Decadal changes in the structure of Cymodocea nodosa seagrass meadows:

herbivory. Ann. Bot. London 86, 913–920. natural vs. human influences. Estuar. Coast. Shelf Sci. 137, 41–49.

Morgan, M.D., Kitting, C.L., 1984. Productivity and utilization of the seagrass Tuya, F., Png-González, L., Riera, R., Haroun, R., Espino, F., 2014b. Ecological

Halodule wrightii and its attached epiphytes. Limnol. Oceanogr. 29, 1066–1076. structure and function differs between habitats dominated by seagrasses and

Paul, V.J., Hay, M.E., 1986. Seaweed susceptibility to herbivory: chemical and green seaweeds. Mar. Environ. Res. 98, 1–13.

morphological correlates. Mar. Ecol. Prog. Ser. 33, 255–264. Valentine, J.F., Heck, K.L., 1999. Seagrass herbivory: evidence for the continued

Paul, V.J., Littler, M.M., Littler, D.S., Fenical, W., 1987. Evidence for chemical defense grazing of marine grasses. Mar. Ecol. Prog. Ser. 176, 291–302.

in tropical green alga Caulerpa ashmeadii (: ): Van Alstyne, K.L., Paul, V.J., 1990. The biogeography of polyphenolic compounds in

isolation of new bioactive sesquiterpenoids. J. Chem. Ecol. 13, 1171–1185. marine macroalgae: temperate brown algal defenses deter feeding by tropical

Paul, V.J., Nelson, S.G., Sanger, H.R., 1990. Feeding preferences of adult and juvenile herbivorous fishes. Oecologia 84, 158–163.

rabbitfish Siganus argenteus in relation to chemical defenses of tropical Vergés, A., Becerro, M.A., Alcoverro, T., Romero, J., 2007. Variation in multiple traits

seaweeds. Mar. Ecol. Prog. Ser. 60, 23–34. of vegetative and reproductive seagrass tissues influences plant-herbivore

Pavia, H., Cervin, G., Lindgren, A., Åberg, P., 1997. Effects of UVB radiation and interactions. Oecologia 151, 675–686.

simulated herbivory on phlorotannins in the brown alga Ascophylum nodosum. Vergés, A., Alcoverro, T., Romero, J., 2011. Plant defenses and the role of epibiosis in

Mar. Ecol. Prog. Ser. 157, 139–146. mediating within-plant feeding choices of seagrass consumers. Oecologia 166,

Pedrotti, M.L., Lemée, R., 1999. Effect of microalgae treated with natural toxins on 381–390.

the nutrition and development of filter-feeding sea-urchin larvae. Mar. Vizzini, S., Sarà, G., Michener, R.H., Mazzola, A., 2002. The trophic role of the

Environ. Res. 48, 177–192. macrophyte Cymodocea nodosa (Ucria) Asch: in a Mediterranean saltworks:

Pennings, S., Paul, V.J., 1992. Effect of plant toughness, calcification, and chemistry evidence from carbon and nitrogen. Bull. Mar. Sci. 71, 1369–1378.

on herbivory by Dolabella auricularia. Ecology 73, 1606–1619. Watson, D.C., Norton, T.A., 1985. The physical characteristics of seaweed thalli as

Prado, P., Alcoverro, T., Martínez-Crego, B., Vergés, A., Pérez, M., Romero, J., 2007. deterrents to littorine grazers. Bot. Mar. 28, 383–387.

Macrograzers strongly influence patterns of epiphytic assemblages in seagrass White, K.S., Westera, M.B., Kendrick, G.A., 2011. Spatial patterns in fish herbivory in

meadows. J. Exp. Mar. Biol. Ecol. 350, 130–143. a temperate Australian seagrass meadow. Estuar. Coast. Shelf Sci. 93, 366–374.

Prado, P., Romero, J., Alcoverro, R., 2010. Nutrient status: plant availability and Wressnig, A., Booth, D.J., 2007. Feeding preferences of two seagrass grazing

seasonal forcing mediate fish herbivory in temperate seagrass beds. Mar. Ecol. monacanthid fishes. J. Fish. Biol. 71, 272–278.

Prog. Ser. 409, 229–239. Wylie, C.R., Paul, V.J., 1988. Feeding preferences of the surgeonfish Zebrasoma

Prado, P., Heck, K.L., 2011. Seagrass selection by omnivorous and herbivorous flavescens in relation to chemical defenses of tropical algae. Mar. Ecol. Prog. Ser.

consumers: determining factors. Mar. Ecol. Prog. Ser. 429, 45–55. 45, 23–32.

Schultz, J.C., 1988. Many factors influence the evolution of herbivore diets, but

plant chemistry is central. Ecology 69, 896–897.