UEF//eRepository
DSpace https://erepo.uef.fi
Rinnakkaistallenteet Luonnontieteiden ja metsätieteiden tiedekunta
2018
A Seven-Year Study of Phenolic
Concentrations of the Dioecious Salix myrsinifolia
Nissinen, Katri
Springer Nature
Tieteelliset aikakauslehtiartikkelit
© Springer Science+Business Media, LLC, part of Springer Nature 2018 All rights reserved
http://dx.doi.org/10.1007/s10886-018-0942-4
https://erepo.uef.fi/handle/123456789/6753
Downloaded from University of Eastern Finland's eRepository
A SEVEN-YEAR STUDY OF PHENOLIC CONCENTRATIONS
1
OF THE DIOECIOUS Salix myrsinifolia
2 3
*KATRI NISSINEN1
4
VIRPI VIRJAMO1
5
LAURI MEHTÄTALO2
6
ANU LAVOLA1
7
ANU VALTONEN1
8
LINE NYBAKKEN3
9
RIITTA JULKUNEN-TIITTO1
10
1Department of Environmental and Biological Sciences, University of Eastern Finland (UEF), 80101, Joensuu, Finland
11
2School of Computing, University of Eastern Finland, 80101, Joensuu, Finland
12
3Faculty of Environmental Sciences and Natural Resource Management, Norwegian University of Life Sciences, 1432
13
Ås, Norway
14 15
*Correspondence: Katri Nissinen, e-mail: katri.nissinen@uef.fi, Telephone +358 50 381 9326
16 17
ORCIDS – Katri Nissinen 0000-0002-6744-2773, Virpi Virjamo 0000-0002-2967-4813, Anu Valtonen 0000-0003-
18
1532-1563
19
20
Abstract–In boreal woody plants, concentrations of defensive phenolic compounds are expected to be at a high level
21
during the juvenile phase and decrease in maturity, although there is variation between plant species. Females of
22
dioecious species, like most of the Salicaceae, are expected to invest their resources in defense and reproduction, while
23
males are expected to be more growth-oriented. We studied age- and sex-dependent changes in leaf and stem phenolics,
24
and in height and diameter growth in a dioecious Salix myrsinifolia plants over a seven-year time period. In addition,
25
we registered flowering as well as rust damage in the leaves. From the first year and throughout ontogenetic
26
development from juvenile to adult phases, there was no significant change in the concentrations of any of the studied
27
compounds in the leaves of S. myrsinifolia. In the stems, the concentrations of six out of 43 identified compounds
28
decreased slightly with age, which may be partly explained by dilution caused by the increment in stem diameter with
29
age. The fairly steady chemistry level over seven years, accompanied by moderate genotypic phenolic variation,
30
indicates important roles of chemical defenses against herbivory for this early-successional species.
31 32
Key Words–Salix myrsinifolia, Dark-leaved willow, Melampsora, Sexual differences, Phenolic compounds.
33
34
INTRODUCTION
35 36
In plants, phenolic compounds, such as salicylates, flavonoids, phenolic acids and condensed tannins, function as
37
defensive compounds against various herbivores, as well as pathogens (Barbehenn and Constabel 2011; Boeckler et al.
38
2011). Flavonoids are also acknowledged to have antioxidant properties e.g. in photoprotection (Agati et al. 2012, 2013;
39
Ferreyra et al. 2012). Defense levels vary across ontogenetic phases (Barton and Koricheva 2010; Boege and Marquis
40
2005; Donaldson et al. 2006). Although there is variation between plant species, the concentration of different phenolic
41
compounds generally increases during the juvenile seedling stage of boreal woody plants (Barton and Koricheva 2010;
42
Donaldson et al. 2006). After the seedling phase, during the period of intensive growth, the concentration of total
43
phenolics often remains at constant level (Barton and Koricheva 2010). Later, during the reproductive phase, the
44
concentration usually decreases as resources are directed to flowering (Boege and Marquis 2005). Especially in boreal
45
areas with harsh winters, phenolic concentrations in woody species have been found to be lower in mature plants than in
46
juvenile ones (Barton and Koricheva 2010). Consequently, mature plants have been preferred for feeding by mammals
47
and birds (Boege and Marquis 2005, Tahvanainen et al. 1985). Proportions of different compounds and compound
48
groups may change during ontogeny, and as an example, salicylates can be replaced by other phenolic compounds such
49
as flavonoids and condensed tannins as the plant develops (Julkunen-Tiitto 1989a).
50
Females and males of dioecious species differ in their morphology, physiology, and defensive strategies
51
(Obeso 2002). According to ecological theory, the females of dioecious species should invest more in reproduction and
52
defense, whereas males are more growth-oriented (Obeso 2002). Experimental data are both in favor of and against this
53
theory (Nybakken and Julkunen-Tiitto 2013; Randriamanana et al. 2015a, 2015b; Robinson et al. 2014). Differences
54
between sexes in their adaptation to natural habitats may cause female- or male-biased sex ratios (Munné-Bosch 2015).
55
Male plants of different plant species are often more susceptible to herbivory (Cornelissen and Stiling 2006, Ågren et al.
56
1999). This difference is related to sexual differences in chemical composition, leaf morphology, vegetative phenology
57
as well as in the resulting nutrient content of the plants (Ågren 1999).
58
The rapid-growing Salix myrsinifolia is a dioecious species native to Europe and western Siberia. S.
59
myrsinifolia has variable growth habits, varying from a small-sized shrub to a small tree, thriving in many different
60
humid habitats (Cronk et al. 2015; Skvortsov 1999). As Salicaceae species in general, S. myrsinifolia is a dioecious
61
species (Skvortsov 1999). It is, as most of the Salix species, sexually female-biased 2:1 (Hughes 2010). In a three-year
62
field study, males of S. myrsinifolia had higher diameter and showed stronger growth response to enhanced temperature
63
during the first two years than did females, while females had more chlorogenic acids in leaves during the third year
64
compared to males (Nybakken et al. 2012; Randriamanana et al. 2015a). In a greenhouse experiment, where plants grew
65
in optimal conditions compared to those in field studies, female plants had clearly higher concentrations of most of the
66
studied phenolic compounds in twigs than males (Nybakken and Julkunen-Tiitto 2013).
67
Like most species belonging to the Salicaceae family, S. myrsinifolia contains high concentrations of
68
salicylates in its leaves and stems (Boeckler et al. 2011; Julkunen-Tiitto 1989a; Julkunen-Tiitto and Virjamo 2017). The
69
salicylates differ both qualitatively and quantitatively between Salicaceae species and among genotypes within the
70
species (Förster et al. 2008; Heiska et al. 2008; Nybakken et al. 2012; Nyman and Julkunen-Tiitto 2005). Usually,
71
salicylates and their degradation products repel generalist herbivores (e.g. Osier and Lindroth 2001; Pentzold et al.
72
2014; Volf et al. 2015). On the other hand, Salix specialist herbivores are attracted by the high concentration of
73
salicylates (Tahvanainen et al. 1985). Some specialist species can even benefit from salicylates and use the glucose
74
moiety as an energy source and the salicylaldehyde moiety in their own defense systems (Boeckler et al. 2000, Rowell-
75
Rahier and Pasteels 1986). By local outbreaks, these specialists can cause high seasonally varying stress on Salix plants
76
(Björkman et al. 2000). In addition to herbivores, pathogens, for example the Melampsora-species, may be detrimental
77
to the growth of Salix species (Dawson and McCracken 1994). Increased concentrations of some phenolic compounds
78
have been reported to decrease the infections of Melampsora-rust in S. myrsinifolia and S. triandra (Hakulinen 1998,
79
Hakulinen et al. 1999, Hjältén et al. 2007).
80
Concentrations of salicylates in many Salicaceae species change during they ontogeny (Nissinen et al.
81
2017, Nybakken et al. 2012, Randriamanana et al. 2015a). A high level of phenolics during the seedling phase is an
82
ancestral state, which some species from the subgenus Salix share with Populus (e.g. Julkunen-Tiitto 1989b). Populus is
83
a sister genus of Salix (Skvortsov 1999, Wu et al. 2015). For example, in P. tremuloides salicylates decreased with age,
84
while condensed tannins increased (Donaldson et al. 2006). On the other hand, in S. myrsinifolia, the salicylate level
85
can be at a high level even in mature plants (Julkunen-Tiitto 1989b). This may be one sign of S. myrsinifolia being a
86
derived species. However, phylogeny of the genus Salix is still under study (e.g. Lauron-Moreau et al. 2015, Wu et al.
87
2015). Skvortsov (1999) divided the genus Salix into three subgenera; Salix, Chamaetia and Vetrix, whereas e.g.
88
Ghahremaninejad et al. (2012), supports a unification of two subgenera, Chaemaetia and Vetrix, into one, based on leaf
89
morphological characters. S. myrsinifolia belongs to the subgenus Vetrix, which includes some of the most derived Salix
90
species (Lauron-Moreau et al. 2015, Skvortsov 1999).
91
Here we present a seven-year field study on the ontogenetic, age-related and sex-dependent patterns of
92
phenolic compounds, and on rust infestation of S. myrsinifolia. To our knowledge, this is the first long-term study on
93
the development of detailed leaf phenolic profiles over the ontogenetic phases in any tree species. In addition, we
94
studied changes in stem phenolics, height and diameter growth, and in flowering. Based on phytochemistry levels
95
measured from young (Nybakken et al. 2012) and adult S. myrsinifolia plants (Julkunen-Tiitto 1989b), we did not
96
expect S. myrsinifolia to undergo great changes over age in leaf and stem phenolics, neither qualitatively nor
97
quantitatively. On the other hand, we expected to find sex-related differences in leaf and stem phenolics, and also in
98
flowering and growth. We hypothesized that (1) males would be more growth-oriented, while females would invest
99
more in defensive compounds and reproduction. (2) Age would not influence the level of secondary compounds in the
100
leaves or stems of S. myrsinifolia. (3) Age would not influence the severity of rust infestation in the leaves, which we
101
expect to depend mainly on varying weather conditions between years.
102 103
METHODS AND MATERIALS
104 105
Plant Material and Experimental Design. In May 2007, we collected approx. 25 cm long cuttings from eight female and
106
nine male genotypes from the second-year-shoots of S. myrsinifolia plants. The samples were obtained from wild
107
populations sampled at nine different locations more than 8 kilometers apart in Eastern Finland. At each location, we
108
sampled a male and a female genotype, keeping a minimal distance of 500 m between the sampled individuals. The
109
same spring, we planted 40 cuttings of each genotype randomly in eight rows in an abandoned hay field in Luikonlahti
110
(62°54'N, 28°40'E) in Eastern Finland.
111 112
Growth Measurements and Sampling. Height and diameter growth measurements and sampling of leaves were done
113
every year at the end of the growing season, starting from the year of planting (Table S1). Sampling of stems was done
114
after the growing season (Table S1). The stem samples from 2009–2011 could not be used because the freezer, where
115
the samples were kept, broke down. Diameter was measured at the root collar, 1 cm above the ground level, using a
116
digital caliper. To detect possible differences between sexes in the starting time of flowering, we counted the number of
117
flowering plants and the number of the catkins in individual plants in 2009 and 2010.
118
Two to six of the youngest fully expanded leaves of the longest shoot were collected for chemical
119
analyses as well as for assessing leaf area, leaf dry weight, rust infection severity and herbivory damage. Leaves were
120
collected from six to ten individual plants of eight female genotypes and nine male genotypes per year (total 150–153
121
individuals per year, full total 1058). The leaves were air-dried in paper bags containing drying material (about 2 g of
122
silicagel, Oy FF-Chemicals Ab, Haukipudas) in a drying room at a relative humidity of 10 % (at room temperature) for
123
12 h. The dried leaves were weighed and their area was measured with a portable leaf area meter LI-3000C (LI-COR,
124
Lincoln, NE, USA). Specific leaf area (SLA, cm2 g-1 dw) was calculated by dividing leaf area by leaf dry weight. The
125
number of Melampsora-rust infection uredinia spots per leaf area was calculated by microscopic examination of the
126
same leaves that were used for chemistry analysis, in all the years 2007–2013. 10 cm long stem internode parts were
127
collected from a side shoot of three to ten individual plants of each genotype per year (total 72–143 individuals per
128
year, full total 494). Stem parts were kept at -20 °C before chemical analyses. Stems of one female genotype and four
129
male genotypes were not included in the study in 2007, as the number of plant individuals with side branches was too
130
low in these genotypes.
131 132
Extraction and Analysis of Leaf and Stem Phenolics. Leaf phenolics were extracted according to earlier published
133
methods (Nybakken et al. 2012). In short, approximately 6 mg of leaf discs was cut from the leaf blades of the sampled
134
2–6 leaves, avoiding leaf veins, mixed and homogenized for 30 s in 600 μl ice-cold 100 % methanol (MeOH) at 5500
135
rpm using a Precellys-homogenizer (Bertin Technologies, France). The samples were incubated in an ice bath for 15
136
min and centrifuged at 16700 g for 3 min at 4 °C (Eppendorf Centrifuge 5415 R, Hamburg, Germany). Supernatants
137
were collected immediately after centrifugation. The sample residues were re-extracted three times. Supernatants from
138
extractions were pooled and evaporated to dryness in a vacuum centrifuge (Eppendorf 270 Concentrator, Hamburg,
139
Germany) at +45 °C and stored at -20 °C. The pellets were dried in a fume hood at room temperature, stored at -20 °C,
140
and later used for quantification of the MeOH-insoluble condensed tannins.
141
Before extraction of stem phenolics, buds and a 2 cm piece from the stem tips were removed, and a 2.5
142
cm piece of stem was cut for chemical analyses. In addition, a 3 cm piece of the stem was cut for dry weight
143
measurement (dried in an oven at 105 °C for 48 hours). The 2.5 cm piece of stem was cut up into small pieces in a cold
144
room (4 °C) and they were put on ice. 20–25 mg of the stem pieces were homogenized for 30 s at 5500 rpm using a
145
Precellys-homogenizer (Bertin Technologies, France). After homogenization, the stem phenolics were extracted
146
similarly to the method used for leaf phenolics.
147
For analyses of leaf and stem phenolics, the dried leaf samples were re-dissolved in 1.2 ml
148
MeOH:MilliQ water (1:1), while the stem samples were re-dissolved in 1.0 ml MeOH:MilliQ water (1:1). Both sample
149
types were run on an HPLC-DAD system that consisted of a binary pump, a vacuum degasser, an autosampler with a
150
thermostat, a column oven and a diode array detector. A column Zorbax SB-C18 (4.6 x 75 mm, particle size 3.5 μm) was
151
used. The eluents used were A (1.5 % tetrahydrofuran + 0.25 % orthophosphoric acid in Milli-Q ultrapure water) and B
152
(100 % MeOH). The gradient for eluent A was as follows: 0–5 min 100 %, 5–10 min 100–85 %, 10–20 min 85–70 %,
153
20–40 min 70–50 %, 40–50 min 50 %, 50–52 min 50–0 %. The column temperature was 30 °C, the flow rate 2 ml min-
154
1, and the injection volume 20 μl.
155
The salicylates and other glucosides were quantified at 270 nm, flavonoids and phenolic acids at 320
156
nm. The concentrations of compounds (mg g-1 dw) were calculated according to the following standards: salicin
157
(Sigma-Aldrich Finland Oy, Helsinki, Finland) for salicin, salicyl alcohol and diglucoside of salicyl alcohol; salicortin
158
(Sigma-Aldrich Finland Oy, Helsinki, Finland) for salicortin, disalicortin, HCH-salicortin derivatives 1 and 2;
159
tremulacin (Apin Chemicals, Abingdon, UK) for tremulacin; (+)-catechin (Fluka Chemie AG, Buchs, Switzerland) for
160
(+)-catechin; quercetin 3-galactoside (Apin Chemicals, Abingdon, UK) for quercetin 3-galactoside, quercetin
161
diglucoside, quercetin 3-glucoside, quercetin glycoside derivative, quercetin 3-arabinoside, quercetin 3-rhamnoside and
162
quercetin rhamnetin derivative; myricetin 3-rhamnoside (Apin Chemicals, Abingdon, UK) for myricetin galactoside and
163
myricetin glucoside+glucuronide; luteolin 7-glucoside (Carl Roth, Karlsruhe, Germany) for luteolin 7-
164
glucoside+quercetin, luteolin glycoside derivatives 1 and 2, monomethyl luteolin 7-glucoside and luteolin derivative;
165
eriodictyol 7-glucoside (Carl Roth, Karlsruhe, Germany) for eriodictyol aglycone derivatives 1 and 2; apigenin 7-
166
glucoside (Carl Roth, Karlsruhe, Germany) for apigenin 7-glucoside; kaempferol 3-O-glucoside (Extrasynthése, Genay,
167
France) for coumaroyl-astragalin; chlorogenic acid (Sigma-Aldrich Finland Oy, Helsinki, Finland) for neochlorogenic
168
acid, chlorogenic acid, and chlorogenic acid derivatives 1, 2, 3, 4 and 5; cinnamic acid (Sigma-Aldrich Finland Oy,
169
Helsinki, Finland) for p-OH-cinnamic acid derivatives 1, 2, 3, 4 and 5; picein (Extrasynthése, Genay, France) for
170
picein; triandrin (from Prof. Beat Meier, Zürich, Switzerland) for triandrin; salidroside (from Prof. Beat Meier, Zürich,
171
Switzerland) for rosavin derivative; salireposide (from Prof. Beat Meier, Zürich, Switzerland) for salireposide.
172
The acid butanol assay (Hagerman 2011) was used to quantify condensed tannin concentration.
173
Concentrations of MeOH-soluble condensed tannins were determined from the re-dissolved MeOH:water extracts used
174
for HPLC-runs and MeOH-insoluble condensed tannins were determined from the dried pellet samples from MeOH
175
extractions. Purified crude tannin extracts from S. myrsinifolia leaves and stems were used as standards.
176 177
Statistical analyses. The effects of age and sex on the concentrations of leaf and stem total salicylates, total flavonoids,
178
total phenolic acids, total condensed tannins, soluble and insoluble condensed tannins and individual phenolic
179
compounds, height and basal diameter, leaf weight, leaf area and specific leaf area (SLA), for genotype i, plant
180
individual j and year k were analyzed using the linear mixed-effect model
181 182
𝑦𝑦𝑖𝑖𝑖𝑖𝑖𝑖=𝜷𝜷′𝒙𝒙𝒊𝒊𝒊𝒊𝒊𝒊+𝑢𝑢𝑖𝑖(1)+𝑢𝑢𝑖𝑖(2)𝑡𝑡𝑖𝑖𝑖𝑖𝑖𝑖+𝑣𝑣𝑖𝑖𝑖𝑖(1)+𝑣𝑣𝑖𝑖𝑖𝑖(2)𝑡𝑡𝑖𝑖𝑖𝑖𝑖𝑖+𝑝𝑝𝑖𝑖+𝑒𝑒𝑖𝑖𝑖𝑖𝑖𝑖 (1)
183
184
where the fixed part β’xijk includes the fixed predictors (age, sex and their interaction), ui(1) and ui(2) are random normally
185
distributed intercept and slope coefficient of age (tijk) at the genotype level and vij(1) and vij(2) are the random intercept
186
and slope at the plant individual level within genotype, pk is the random effect for calendar year and eijk is the residual
187
error for plant individual j in genotype i in year k. The age was coded as a continuous predictor, with an additional
188
binary predictor for the first year, when the plants were still at the juvenile phase.
189
In addition, we analyzed flowering and leaf rust infection severity by the linear mixed effect model.
190
Flowering and leaf rust infection severity were clearly a mixture of two random processes: i) the effects on presence
191
and absence of catkins or severe rust infection and ii) the effects on magnitude of the variable of interest. These data
192
were therefore analyzed in two parts. First, when analyzing the presence of catkins and the occurrence of severe
193
infection (≥ 1000 rust spots per leaf), we used a binary logistic mixed-effect model. In this model, the fixed part and
194
random effects were specified as in the previously introduced linear mixed-effect model (1). Second, we used this linear
195
mixed-effect model to analyze the number of catkins for those plants where these occurred and the number of rust spots
196
for leaves with non-severe rust infection.
197
The model selection was based on the conditional F tests and t tests, and on Likelihood ratio (Chi2)
198
tests, if F tests or t tests were not possible (Pinheiro and Bates 2000). The model fit was evaluated graphically and when
199
needed, a logarithmic (ln(x), ln(x+0.1) or ln(x+1)) transformation (Table S2) was used in the response variable for
200
testing purposes. For each variable, we selected the transformation that visually showed the most homoscedastic
201
residuals of the variable in question. However, the final parameter estimates are reported using a model fitted into the
202
original scale, to enable meaningful interpretation (Table S2). The approach was selected because the coefficients and
203
predictions from a model fitted in the original scale are unbiased, even though not fully efficient, whereas back-
204
transformed predictions from a transformed model would be biased. The statistical analyses were conducted using R
205
version 3.3.1 (R Core Team 2016) with the packages lme4 (Bates et al. 2015) and lmerTest (Kuznetsova et al. 2016) in
206
RStudio 0.99.903 (© 2009–2016 RStudio, Inc., Boston, USA) and IBM® SPSS® Statistics 22 (Armonk, NY, USA).
207
We did a graphic vector analysis (GVA) (Figure 1) to study the effects of age on the different phenolic
208
group levels; salicylates, flavonoids, phenolic acids and condensed tannins, in leaves of S.myrsinifolia (Haase and Rose
209
1995, Koricheva et al. 1999). The relative concentration and content values in GVA were calculated by using the
210
concentration (mg g-1) and content (mg) values of the phenolic group in the first-year plants as a reference. We plotted
211
the relative values of concentration (y) and content (x) of different phenolic compound groups in relation to leaf dry
212
weight of the leaves (z) used for phenolic analysis. In a GVA plot, an excess synthesis of a certain phenolic group
213
compared to the first year is revealed when both its relative content and concentration are increasing. Decreased relative
214
content and concentration implies reduced synthesis in the phenolic compound group of interest. A concentration effect
215
in the tissue is revealed if relative content decreases while relative concentration increases. On the other hand, an
216
increase in relative content and a decrease in relative concentration indicates dilution. The relative values used in a
217
GVA plot were calculated as studied year mean/first year mean x 100. The GVA plots were drawn using SigmaPlot
218
13.0 (Systat Software, Inc., Chicago, IL, USA). Because we did not collect the stem biomass samples, we did not have
219
stem biomass data for calculation of the relative content of phenolics for the GVA. Therefore, we did not analyze stem
220
phenolics with the GVA method.
221
We used a principal component analysis (PCA) to get an overview of the total leaf and stem phenolics
222
dataset (Figure S3, Figure S4). In the PCA we plotted all the observations according to all 33 different leaf and 43
223
different stem phytochemical compound concentrations to visualize the relationships of different genotypes in seven-
224
year time period. Sex, genotype and year were set as secondary observations identifiers. In addition, we scored the
225
leaves based on the amount of rust pustules on them (0–299=1, 300–799=2, >800=3). This rust score data was set as a
226
secondary observation identifier in the PCA to visualize the effect of rust infection severity on leaf phenolics. The PCA
227
was conducted in the SIMCA-P+-progman (Umetrics AB, Sweden). Furthermore, to obtain the proportion of variation
228
in leaf phenolic profiles explained by rust infection, age, sex and genotype, we used a permutational ANOVA-model
229
(PERMANOVA+ in PRIMER-E-software, version 6) (Anderson et al. 2008) to fit a hierarchical design for our data
230
(type III sums of squares and using Bray-Curtis as a measure of similarity). In PERMANOVA, rust infection and age
231
were added as continuous variables, sex as a fixed factor with two levels, genotype nested in sex and individual nested
232
in genotype as random factors.
233 234
RESULTS
235 236
Leaf and Stem Phenolics. In the leaves of S. myrsinifolia, we detected and quantified 31 different low molecular weight
237
phenolic compounds, in addition to soluble and insoluble high molecular weight condensed tannins (proanthocyanidins)
238
(Table 1). Salicylates were the main group of low molecular weight phenolics, representing 72 % of their total
239
concentrations in the leaves (Figure 2). Other groups were phenolic acids and flavonoids, representing shares of 19 %
240
and 9 %, respectively. Among the salicylates, salicortin and HCH-salicortin were the most abundant. Chlorogenic acids
241
were the most abundant phenolic acids, while quercetin 3-galactoside and (+)-catechin were the most abundant
242
flavonoids.
243
In the stems of S. myrsinifolia the identified number of different compounds, 43, was higher than in the
244
leaves (Table 2). Salicylates formed the most abundant low molecular weight phenolic group also in the stems.
245
Salicortin and salicin were present in the highest quantities. The stems contained phenolic glycosides not detected in the
246
leaves. Main phenolic glycosides were picein and a rosavin derivative. In contrast to the leaves, the stems had higher
247
concentrations of flavonoids than of phenolic acids, but the most abundant compounds were the same in both of these
248
phenolic groups in leaves and stems.
249
We found almost no significant developmental trends in the concentrations of the major phenolic
250
groups or of the individual phenolic compounds in the leaves of S. myrsinifolia. In the early development of the S.
251
myrsinifolia plants, between the first and the second year, we detected a significant increase in the concentration of only
252
three minor compounds of the analyzed leaf phenolics; in HCH-salicortin 2 (t(4.012)=-9.028 P<0.001), in monomethyl
253
luteolin 7-glucoside (t(3.997)=-3.453 P<0.05) and in luteolin 7-glucoside+quercetin glycoside (t(3.995)=-2.968 P<0.05)
254
(Table S2). No other significant effect of age was detected on leaf phenolics. The results of GVA analysis showed
255
increase in the concentrations of total phenolic acids and total salicylates only between 2007 and 2008. There was no
256
clear trend in the following years (Figure 1; Table S2). When visualizing the whole leaf individual phenolic compound
257
dataset in a PCA plot, years did not differ from each other (Figure S3). Based on PERMANOVA model, gneotype
258
explained the largest part of the variation in compositions of leaf phenolics (40.7%), while the proportion of variation
259
explaned by age was only 10.2% (Table 3).
260
We detected some developmental trends in the stem phenolics of S. myrsinifolia. Concentrations of six
261
out of the 43 low molecular weight phenolic compounds decreased significantly with increasing age (Table 2; Table
262
S2). These were salicin (F(2,29.880)=60.616 P<0.001), HCH-salicortin1 (χ2(2)=14.207 P<0.001), monomethyl 7-glucoside
263
(χ2(2)=15.061 P<0.001), methyl luteolin 7-glucoside (χ2(2)=7.116 P<0.05), chlorogenic acid derivative2
264
(F(2,27.069)=19.057 P<0.001) and salireposide (χ2(2)=9.830 P<0.01) (Table 2; Table S2). Two compounds, salicin
265
(t(422.2)=7.450 P<0.001) and chlorogenic acid derivative2 (t(363.9)=-3.823 P<0.001), decreased between the first and the
266
second year. In addition, phenolic glycosides as a group (χ2(2)=7.788 P<0.05) and MeOH-insoluble condensed tannins
267
(χ2(2)=9.916 P<0.01) decreased with age in stems. Furthermore, though not significant, there was a decreasing trend in
268
the concentrations of all the low molecular weight phenolic groups and most of the individual compounds in the stems
269
(Table 2; Table S2; Figure S4).
270
Generally, we found only minor differences in both leaf and stem phenolics between male and female
271
genotypes. In females, we detected significantly higher concentrations of the chlorogenic acid derivative2 (t(15.074)=-
272
2.560 P<0.05) in the leaves and of neochlorogenic acid (t(15.025)=-2.24 P<0.05) in the stems. On the other hand, the
273
concentrations of total flavonoids, (+)-catechin and condensed tannins in the leaves tended to be higher in males (Table
274
1; Table S2). The similarity in the phenolic profiles of the sexes can be seen in the PCA plot (Figure S3), and is also
275
supported by the results of the PERMANOVA model (Table 3).
276
In contrast to low sexual differences, there were large quantitative differences between genotypes in
277
both the leaf and the stem phenolics. This was proven by the PERMANOVA analysis, where the largest part of the
278
variation in the leaf phenolic profiles was explained by genotype (Table 3). High genotypic variation is seen in the PCA
279
plots also and in the high standard error values in the concentrations (Table 1; Table 2; Table S2; Figure S3; Figure S4).
280
In the stems, all the individual phenolic acids, phenolic glycosides and salicylates were detected in all of studied the
281
genotypes, while in the leaves only all the individual phenolic acid compounds were detected in all of the genotypes.
282
One individual leaf salicylate, tremulacin, and almost half of the total number of individual leaf and stem flavonoid
283
compounds, such as quercitrin in the leaves and quercetin 3-arabinoside in the stems, were detected in only some of the
284
genotypes. Condensed tannin concentration in the leaves of some of the genotypes was ten times higher than in the
285
others. In the PCA of leaf phenolic compounds, four genotypes stood out from the others quite clearly (Figure S3).
286
These distinctions were mainly explained by the high concentrations of (+)-catechin and condensed tannins in
287
genotypes 15 and 19, and tremulacin in genotypes 10 and 11.
288 289
Growth and Flowering. The height and diameter growth of the plants was steady over the years (Figure 3), and there
290
were no significant trends in the studied leaf parameters (leaf dry weight, leaf area and specific leaf area, SLA) with
291
increasing age (Figure 4). On the other hand, yearly weather conditions seemed to have an effect on leaf size; leaves
292
were largest and thickest during the warm summers 2010–2011 (mean temperatures during the months from June to
293
August 16.5 °C and 16.8 °C, respectively, compared to the mean for thirty years, 15.5 °C), and during a rainy summer
294
2012 (mean precipitation 301 mm from June to August, compared to the mean for thirty years, 75.5 mm). Females and
295
males did not differ in height and diameter growth, nor in any of the leaf size parameters (Figure 3, 5). S. myrsinifolia
296
started to flower at the age of three to four years, in 2009–2010. Number of catkins per shoot was higher in male
297
genotypes in 2009, while in 2010 there was no difference in flowering between sexes (significant year*sex interaction,
298
F (1, 663.5)=5.782 P<0.05; Figure 5).
299 300
Rust Damage. There was no significant effect of age or difference between females and males in the number of rust-
301
infected plants, nor in the severity of rust-infection in the leaves of S. myrsinifolia (Figure 6). There was, however,
302
some variation between years. Rust infections were severe in 2012, when the mean summer temperature was low (14.7
303
°C) and it rained frequently (301 mm) (Figure 6). On the other hand, rust infections were rare in 2010 and 2011, when
304
the mean summer temperature was high and it rained less, as mentioned above (156 mm and 196 mm, respectively)
305
(Figure 6). Neither the PERMANOVA (only 7.9% of explained variation) nor the PCA-analysis showed any strong
306
effects of the degree of rust-infection on the phenolic profile of the leaves in comparison to the effect of genotype
307
(Table 3; Figure S3).
308 309
DISCUSSION
310 311
This study gives a unique overview of the defensive phenolic profiles of leaves and stems of a Salicaceae species across
312
plant ontogenetic stages, and over a long-term study period (up to seven years). To our knowledge, this is the first time
313
that detailed phenolic profiles were analyzed from the same plant individuals and from the same microclimatic
314
environment for such a long period. As we hypothesized, age did not influence the concentrations of phenolic
315
compounds in the leaves or stems of S. myrsinifolia, which maintained constant defense levels over ontogeny after the
316
sapling phase and after reaching the reproductive phase. The concentration levels of different low molecular weight
317
phenolic groups and their relative proportions in the leaves and stems of S. myrsinifolia corresponded well with those
318
found in previous short-term studies of the same species in the field (Nybakken et al. 2012; Randriamanana et al.
319
2015a; Tegelberg et al. 2003). Yearly variation in the climatic conditions in the field may have reflected in changes in
320
leaf size, and may explain the detected slight variation in the concentrations of most of the studied compounds in the
321
leaves.
322
We detected only a few changes in stem phenolics and basically no changes in leaf phenolics with age
323
in S. myrsinifolia. The few significant changes that were detected with age in the concentrations of low molecular
324
weight phenolic compounds in the stems were directed towards lower amounts in older plants, contrary to what we
325
hypothesized. This may partly be a result of dilution (Koricheva 1999).
326
Salicylates, which represent a key group of the studied phenolics, have been shown to be largely
327
ineffective against specialized insect herbivores, such as the leaf beetle Phratora vitellineae (Sipura and Tahvanainen
328
2000). It is thus probable that the salicylates in S. myrsinifolia were primarily targeted at generalist insect herbivores
329
(Ruuhola et al. 2001; Tahvanainen et al. 1985; Volf et al. 2015) or small mammalian herbivores, such as hares (Lepus
330
ssp.) and voles (Microtus ssp.), which have been shown to be important herbivores of S. myrsinifolia (Heiska et al.
331
2008; Shaw et al. 2010; Tahvanainen et al. 1985). The constantly high levels of leaf phenolics in S. myrsinifolia are in
332
contrast to decreasing levels of phenolics in basal Salicaceae species. For example, salicylates are replaced by other
333
phenols, such as flavonoids and condensed tannins in the later ontogenetic phases of S. pentandra and P. tremula
334
(Donaldson et al. 2006, Julkunen-Tiitto and Virjamo 2017, Nissinen et al. 2017). These large tree species are less
335
exposed to herbivory by smaller mammals when mature, which may explain the decrease in salicylates with age. On the
336
other hand, the relatively smaller S. myrsinifolia remains exposed to mammalian herbivores even when mature, possibly
337
explaining its steady levels of phenolics. Indeed, larger Salix species have been reported to have lower concentrations of
338
phenolic compounds compared to the smaller sized S. daphnoides and S. purpurea (Förster et al. 2008, 2010;
339
Tahvanainen et al. 1985). Unfortunately, the phylogeny of the Salix species remains largely unresolved, and the
340
information of developmental changes in willow defenses is rather scarce (Lauron-Moreau et al. 2015; Liu et al. 2016;
341
Wu et al. 2015). We are thus not able to say whether the contrasting defensive strategies between S. myrsinifolia and the
342
basal Salicaceae species truly result from their different growth forms or from an evolutionary shift in defensive
343
strategies (Lauron-Moreau et al. 2015; Liu et al. 2016; Wu et al. 2015).
344
We hypothesized that males would be more growth-oriented, while females would invest more in
345
defensive compounds and reproduction. However, females and males of S. myrsinifolia differed little in the studied
346
parameters. Our results are in accordance with earlier field studies, but disagree with results from greenhouse studies
347
(Nybakken and Julkunen-Tiitto 2013; Nybakken et al. 2012; Randriamanana et al. 2014, 2015a). In greenhouses, plants
348
have optimal growing conditions, and thus the responses are not fully comparable to those found in the field. Several
349
concurrent stress factors (e.g. herbivory, plant deseases, competition by other plants, heat, heavy rain and strong wind)
350
in the field seem to level out possible sexual differences in general. Only one chlorogenic acid derivative was
351
significantly more abundant in the leaves, as was neochlorogenic acid in the stems of females. Chlorogenic acid
352
concentrations were also higher in females in field studies of one- to three-year-old S. myrsinifolia (Nybakken et al.
353
2012; Randriamanana et al. 2015a). Chlorogenic acids act as a deterrent against various leaf beetles and thrips (Leiss et
354
al. 2009; Matsuda and Senbo 1985). Chlorogenic acids may also serve as a precursor for the synthesis of flavonoids in
355
the phenylpropanoid pathway. In consequence, the higher concentrations of chlorogenic acids in females may indicate
356
that they rely more on inductive defenses than males do (Nybakken et al. 2012). Low number of significant sexual
357
differences detected may be also a consequence of the considerable differences in the concentrations of all the analyzed
358
compounds among the genotypes.
359
We found substantial genotypic differences in the concentrations of all analyzed phenolic compounds
360
in the leaves and stems of S. myrsinifolia. For example, tremulacin was found only in two genotypes. In addition, we
361
found large quantitative differences in the concentration of condensed tannins. The notable high genotypic variation in
362
phenolics among the Salicaceae species has been shown in several studies (e.g. Hakulinen and Julkunen-Tiitto 2000;
363
Paajanen et al. 2011; Randriamanana et al. 2014). An earlier study of S. myrsinifolia concluded that the genotypic
364
differences in phenolic concentrations were even greater than the differences between Salix species (Tegelberg et al.
365
2003). Genotypic variation allows plant species to survive in a larger range of environments and in changing conditions
366
(Liao et al. 2016). Furthermore, previous studies have demonstrated that closely related congeneric plants tend to
367
diverge in their defenses, which probably allows them to escape herbivory (Becerra 2007; Kursar et al. 2009; Salazar et
368
al. 2016; Volf et al. 2018). Similarly, intraspecific variation can reduce a pool of shared herbivores between S.
369
myrsinifolia genotypes.
370
As we hypothesized, age did not affect the severity of rust infections in S. myrsinifolia during the
371
seven-year period. However, there were substantial differences among years in the severity of rust infections. This is
372
probably a result of varying microclimate factors, such as humidity and temperature conditions. In general, for different
373
rust species, sufficient moisture is a prerequisite to thrive (Helfer 2014). In earlier studies of S. myrsinifolia and P.
374
tremula, rust severity was lower in rapid-growing genotypes and under elevated temperature (Heiska et al. 2007;
375
Nybakken et al. 2012; Randriamanana et al. 2015a, 2015b). Accordingly, we observed a clear relationship between
376
weather conditions and severity of rust infections. The leaves were heavily infected by rust in cold and rainy summers,
377
whereas during warmer and drier summers the infection rates were low. Many of the leaf phenolics have been found to
378
be inducible by different stress factors (Boege 2004; Coleman and Jones 1991; Karban and Myers 1989; Raupp and
379
Sadof 1991). The accumulation of both condensed tannins and salicin may have been induced by rust-damage in the
380
six-year old S. myrsinifolia plants (e.g. Hjältén et al. 2007), and the increase in total phenolics in the seven-year old
381
plants may suggest a delayed induction response to the severe rust attacks of the previous year. Severe damage may
382
have caused a regeneration process in the plants and thereby increased the synthesis of defensive compounds. For
383
example, Miranda et al. (2007) found a connection between rust-infection and the activity of genes encoding condensed
384
tannin synthesis in Populus trichocarpa x P. deltoides. On the other hand, high levels of salicin may also be a result of
385
the degradation of more complex salicylates, such as salicortin, which concentration decreased at the same time.
386 387
Conclusions. S. myrsinifolia has a high constitutive defense profile, which may have evolved as a consequence of being
388
constantly exposed to a high herbivore pressure. Concentrations of different phenolic compounds remained more or less
389
at a constant level over a seven-year period in both leaves and stems. The moderate yearly variation detected in
390
phenolic concentrations seemed to be highly affected by changing weather conditions as well as rust infestation levels.
391
The high recovered variability in polyphenol profiles among S. myrsinifolia genotypes makes this species an ideal
392
model for studying how genotypic variability in natural willow communities affect their chemical diversity and how
393
this is mediated to other trophic levels, such as rich communities of willow herbivores.
394 395
Acknowledgements–Our sincere thanks are owed to Sinikka Sorsa, Hannele Hakulinen, Riitta Pietarinen and Outi
396
Nousiainen for their help with laboratory analyses and Mervi Kupari, Teija Ruuhola, Anneli Salonen, Norul Sobuj,
397
Meeri Koivuniemi, Riia Hörkkä, Anu Ovaskainen and Henna-Riikka Leppänen for their help with field measurements,
398
leaf and stem sampling and in rust and herbivory damage assessment. We are also indebted to the reviewers for their
399
valuable and thorough comments on the manuscript. The English language of the article was checked by Rosemary
400
Mackenzie. This work was financed by the Academy of Finland (no. 267360).
401 402
REFERENCES
403 404
Agati G, Azzarello E, Pollastri S, Tattini M (2012) Flavonoids as antioxidants in plants: Location and functional
405
significance. Plant Sci. 196:67–76. doi:10.1016/jplantsci.2012.07.014
406
Agati G, Brunetti C, Di Ferdinando M, Ferrini F, Pollastri S, Tattini M (2013) Functional roles of flavonoids in
407
photoprotection: New evidence, lessons from the past. Plant Physiol. Biochem. 72:35–45.
408
doi:10.1016/j.plaphy.2013.03.014
409
Ågren JK, Danell T, Elmqvist T, Ericson L, Hjältén J (1999) Sexual dimorphism and biotic interactions, In: Geber MA,
410
Dawson TE, Delph L (ed), Gender and Sexual Dimorphism in Flowering Plants, Springer Verlag, Berlin, pp 217–246
411
Anderson MJ, Gorley RN, Clarke KR (2008) PERMANOVA+ for PRIMER: Guide to software and statistical methods.
412
PRIMER-E, Plymouth, UK.
413
Barbehenn RV, Constabel CP (2011) Tannins in plant-herbivore interactions. Phytochemistry 72:1551–1565.
414
doi:10.1016/j.phytochem.2011.01.040
415
Barton KE, Koricheva J (2010) The ontogeny of plant defence and herbivory: Characterizing general patterns using
416
meta-analysis. Am. Nat. 175(4):481–493. doi:10.1086/650722
417
Bates D, Maechler M, Bolker B, Walker S (2015). Fitting Linear Mixed-Effects Models Using lme4. J. Stat. Softw.
418
67(1):1-48. doi:10.18637/jss.v067.i01
419
Becerra JX (2007) The impact of herbivore-plant coevolution of plant community structure. PNAS 104(18):7483–7488.
420
doi:10.1073/pnas.0608253104
421
Björkman C, Bengtsson B, Häggström H (2000) Localized outbreak of a willow leaf beetle: plant vigor or natural
422
enemies? Popul. Ecol. 42:91–96. doi:10.1007/s101440050013
423
Boeckler GA, Gershenzon J, Unsicker SB (2011) Phenolic glycosides of the Salicaceae and their role as anti-herbivore
424
defenses. Phytochemistry 72:1497–1509. doi:10.1016/j.phytochem.2011.01.038
425
Boege K (2004) Induced responses in three tropical dry forest plant species – direct and indirect effects on herbivory.
426
Oikos 107:541–548. doi:10.1111/j.0030-1299.2004.13272.x
427
Boege K, Marquis J (2005) Facing herbivory as you grow up: the ontogeny of resistance in plants. Trends Ecol. Evol.
428
20(8):441–448. doi:10.1016/j.tree.2005.05.001
429
Coleman JS, Jones CG (1991) A phytocentric perspective of phytochemical induction by herbivores. In: Tallamy DW,
430
Raupp MJ (ed) Phytochemical induction by herbivores, John Wiley & Sons, New York, pp 3–45
431
Cornelissen T, Stiling P (2005) Sex-biased herbivory; a meta-analysis of the effects of gender on plant-herbivore
432
interactions. Oikos 111:488–500. doi:10.1111/j.1600-0706.2005.14075.x
433
Cronk Q, Ruzzier E, Belyaeva I, Percy D (2015) Salix transect of Europe: latitudinal patterns in willow diversity from
434
Greece to arctic Norway. Biodivers. Data J. doi:10.3897/BDJ.3e6258.
435
Dawson WM, McCracken AR (1994) Effect of Melampsora rust on the growth and development of Salix burjatica
436
‘Korso’ in Northern Ireland. Eur. J. For. Pathol. 24:32–39. doi:10.1111/j.1439-0329.1994.tb01320.x
437
Donaldson JR, Stevens MT, Barnhill HR, Lindroth RL (2006) Age-related shifts in leaf chemistry of clonal aspen
438
(Populus tremuloides). J. Chem. Ecol. 32:1415–1429. doi:10.1007/s.10886-006-9059-2
439
Ferreyra MLF, Rius SP, Casati P (2012) Flavonoids: biosynthesis, biological functions, and biotechnological
440
applications. Front. Plant Sci. 3 (222):1–15. doi:10.3389/fpls.2012.00222
441
Förster N, Ulrichs C, Zander M, Kätzel R, Mewis I (2008) Influence of the season on the salicylate and phenolic
442
glycoside contents in the bark of Salix daphnoides, Salix pentandra, and Salix purpurea. J. Appl. Bot. Food Qual.
443
82:99–102
444
Förster N, Ulrichs C, Zander M, Kätzel R, Mewis I (2010) Factors influencing the variability of antioxidative phenolic
445
glycosides in Salix species. J. Agric. Food Chem. 58:8205–8210. doi:10.1021/jf100887v
446
Ghahremaninejad F, Khalili Z, Maassoumi AA, Mirzaie-Nodoushan H, Riahi M (2012) Leaf epidermal features of Salix
447
species (Salicaceae) and their systematic significance. Am. J. Bot. 99(4):769–777. doi:10.3732/ajb.1100019
448
Haase DL, Rose R (1995) Vector analysis and its use for interpreting plant nutrient shifts in response to silvicultural
449
treatments. For. Sci. 41(1):54–66
450
Hagerman A.E. (2011) The Tannin Handbook. Department of Chemistry & Biochemistry, Miami University, Oxford,
451 452
OH.Hakulinen J (1998) Nitrogen-induced reduction in leaf phenolic level is not accompanied by increased rust trquency in
453
a compatible willow (Salix myrsinifolia)-Melampsora rust interaction. Physiol. Plant. 102:101–110. doi:10.1034/j.1399-
454
3054.1998.1020114.x
455
Hakulinen J, Sorjonen J, Julkunen-Tiitto R (1999) Leaf phenolics of three willow clones differing in resistance to
456
Melampsora rust infection. Physiol. Plant. 105:662–669
457
Hakulinen J, Julkunen-Tiitto R (2000) Variation in leaf phenolics of field-cultivated willow (Salix myrsinifolia) clones
458
in relation to occurrence of Melampsora rust. Eur. J. For. Pathol. 30:29–41
459
Heiska S, Tikkanen O-P, Rousi M, Turtola S, Tirkkonen V, Meier B, Julkunen-Tiitto R (2007) The susceptibility of
460
herbal willow to Melampsora rust and herbivores. Eur. J. Plant. Pathol. 118:275–285. doi:10.1007/s10658-007-9145-5
461
Heiska S, Tikkanen O-P, Rousi M, Julkunen-Tiitto R (2008) Bark salicylate and condensed tannins reduce vole
462
browsing amongst cultivated dark-leaved willows (Salix myrsinifolia). Chemoecology 17:245–253.
463
doi:10.1007/s00049-007-0385-9
464
Helfer S (2014) Rust fungi and global change. New Phytol. 201:770–780. doi:10.1111/nph.12570
465
Hjältén J, Niemi, L, Wennström A, Ericson L, Roininen H, Julkunen-Tiitto R (2007) Variable responses of natural
466
enemies to Salix triandra phenotypes with different secondary chemistry. Oikos 116:751–758.
467
doi:10.1111/j.2007.0030-1299.15365.x
468
Hughes FMR, Johansson M, Xiong S et al (2010) The influence of hydrological regimes on sex ratios and spatial
469
segregation of the sexes in two dioecious riparian shrub species in northern Sweden. Plant Ecol. 208:77–92.
470
doi:10.1007/s.11258-009-9689-x
471
Julkunen-Tiitto R (1989a) Phenolic compounds of the genus Salix: A chemotaxonomical survey of further Finnish
472
species. Phytochemistry 28:2115–2125
473
Julkunen-Tiitto R (1989b) Distribution of certain phenolics in Salix species (Salicaceae) Dissertation, University of
474
Joensuu
475
Julkunen-Tiitto R, Virjamo V (2017) Biosynthesis and roles of Salicaceae salicylates. In: Arimura G, Maffei M (ed)
476
Plant specialized metobolism. Genomics, biochemistry, and biological functions. Taylor & Francis Group, Boca Raton,
477
pp 65–83
478
Karban R, Myers JH (1989) Induced plant responses to herbivory. Annu. Rev. Ecol. Syst. 20:331–348.
479
doi:10.1146/annurev.es.20.110189.001555
480
Koricheva J (1999) Interpreting phenotypic variation in plant allelochemistry; problems with the use of concentrations.
481
Oecologia 119:467–473. doi:10.1007/s004420050809
482
Kursar TA, Dexter KG, Lokvam J, Pennington TR, Richardson JE, Weber MG, Murakami ET, Drake C, McGregor R,
483
Coley PD (2009) The evolution of antiherbivore defenses and their contribution to species coexistence in the tropical
484
tree genus Inga. PNAS 106(43):18073–18078. doi:10.1073/pnas.0904786106
485
Kuznetsova A, Brockhoff PB, Christensen RHB (2016). lmerTest:Tests in Linear Mixed Effects Models. R package
486
version 2.0-32. https://CRAN.R-project.org/package=lmerTest
487
Lauron-Moreau A., Pitre FE, Argus GW, Labrecque M, Brouillet L. (2015) Phylogenetic relationships of American
488
willows (Salix L., Salicaceae) PLoS One 10(4):e0121965. doi:10.1371/journal.pone.0121965
489
Leiss KA, Maltese F, Choi YH, Verpoorte R, Klinkhamer PGL (2009) Identification of chlorogenic acid as a resistance
490
factor for thrips in chrysanthemum. Plant Physiol. 150:1567–1575. doi:10.1104/pp.109.138131
491
Liao H, D’Antonio CM, Chen B, Huang Q, Peng S (2016) How much do phenotypic plasticity and local genetic
492
variation contribute to phenotypic divergences along environmental gradients in widespread invasive plants? A meta-
493
analysis. Oikos 125:905–917. doi:10.1111/oik.02372
494
Liu X, Wang Z, Wang D, Zhang J (2016) Phylogeny of Populus-Salix (Salicaceae) and their relative genera using
495
molecular datasets. Biochem. Syst. Ecol. 68:210–215. doi: 10.1016/j.bse.2016.07.016
496
Matsuda K, Senbo K (1986) Chlorogenic acid as a feeding deterrent for the Salicaceae-feeding leaf beetle, Lochmaeae
497
capreae cribrata (Coleoptera: Chrysomelidae) and other species of leaf beetles. Appl. Ent. Zool. 21(3):411–416
498
Miranda M, Ralph SG, Mellway R, White R, Heath MC, Bohlmann J, Constabel PC (2007) The transcriptional
499
response of hybrid poplar (Populus trichocarpa x P. deltoides) to infection by Melampsora medusa leaf rust involves
500
induction of flavonoid pathway genes leading to the accumulation of proanthocyanidins. Mol. Plant-Microbe Interact.
501
20(7):816–831. doi:10.1094/MPMI-20-7-0816
502
Munné-Bosch S (2015) Sex ratios in dioecious plants in the framework of global change. Environ. Exp. Bot. 109:99–
503
102. doi:10.1016/j.envexpbot.2014.08.007
504
Nissinen K, VirjamoV, RandriamananaT, Sobuj N, Sivadasan U, Mehtätalo L, Beuker B, Julkunen-Tiitto R, Nybakken
505
L. (2017) Responses of growth and leaf phenolics in European aspen (Populus tremula L.) to climate change during
506
juvenile phase change. Can. J. For. Res. doi:10.1139/cjfr-2017-0188
507
Nybakken L, Julkunen-Tiitto R (2013) Gender differences in Salix myrsinifolia at the pre-reproductive stage are little
508
affected by simulated climatic change. Physiol. Plant 147:465–476. doi:10.1111/j.1399-3054.2012.01675.x
509
Nybakken L, Hörkkä R, Julkunen-Tiitto R (2012) Combined enhancements of temperature and UVB influence growth
510
and phenolics in clones of the sexually dimorphic Salix myrsinifolia. Physiol. Plant 145:551–564.
511
doi:10.1111/j.1399.3054.2011.01565
512
Nyman T, Julkunen-Tiitto R (2005) Chemical variation within and among six northern willow species. Phytochemistry
513
66:2836–2843. doi:10.1016/j.phytochem.2005.09.040
514
Obeso JR (2002) The costs of reproduction in plants. New Phytol. 155:321–348. doi:10.1046/j.1469-
515
8137.2002.00477.x
516
Osier TL, Lindroth R. (2001) Effects of genotype, nutrient availability, and defoliation of aspen phytochemistry and
517
insect performance. J. Chem. Ecol. 27(7):1289–1313. doi:10.1023/A:1010352307301
518
Paajanen R, Julkunen-Tiitto R, Nybakken L, Petrelius M, Tegelberg R, Pusenies J, Rousi M, Kellomäki S (2011) Dark-
519
leaved willow (Salix myrsinifolia) is resistant to three-factor (elevated CO2, temperature and UV-B-radiation) climate
520
change. New Phytol. 190:161–168. doi:10.1111/j.1469-8137.2010.03583.x
521
Pentzold S, Zagrobelny M, Rook F, Bak S (2014) How insects overcome two-component plant chemical defense; plant
522
ß-glucosidases as the main target for herbivore adaptation. Biol. Rev. 89:531–551. doi:10.1111/brv.12066
523
Pinheiro JC, Bates DM (2000) Mixed-Effects Models in S and S-PLUS. Springer, New York
524
R Core Team (2016). R: A language and environment for statistical computing. R Foundation for Statistical Computing,
525
Vienna, Austria. URL https://www.R-project.org/
526
Randriamanana TR, Nybakken L, Lavola A, Aphalo P, Nissinen K, Julkunen-Tiitto R (2014) Sex-related differences in
527
growth and carbon allocation to defense on Populus tremula as explained by current plant defense theories. Tree
528
Physiol. 34(5):471–487. doi:10.1093/treephys/tpu034
529
Randriamanana TR, Nissinen K, Moilanen J, Nybakken L, Julkunen-Tiitto R (2015a) Long-term UV-B and temperature
530
enhancements suggest that females of Salix myrsinifolia plants are more tolerant to UV-B than males. Environ. Exp.
531
Bot. 109:296–305. doi:10.1016/j.envexpbot.2014.06.007
532
Randriamanana TR, Lavola A, Julkunen-Tiitto R (2015b). Interactive effects of supplemental UV-B and temperature in
533
European aspen seedlings: Implications for growth, leaf traits, phenolic defense and associated organisms. Plant
534
Physiol. Biochem. 93:84–93. doi:10.1016/j.plaphy.2015.03.001
535
Raupp MJ, Sadof CS (1991) Responses of leaf beetles to injury-related changes in their Salicaceous hosts. In: Tallamy
536
DW, Raupp MJ (ed) Phytochemical induction by herbivores, John Wiley & Sons, New York, pp 71–83
537
Robinson KM, Delhomme N, Mähler N, Schiffthaler B, Önskog J, Albrectsen BR, Ingvarsson PK, Hvidsten TR,
538
Jansson S, Street NR (2014) Populus tremula (European aspen) shows no evidence of sexual dimorphism. BMC Plant
539
Biol. 14:276. doi:10.1186/s12870-014-0276-5
540
Rowell-Rahier M, Pasteels JM (1986) Economics of chemical defense in Chrysomelinae. Journal of Chemical Ecology
541
12: 1189–1203. doi:10.1007/BF01639004
542
Ruuhola T, Tikkanen O-P, Tahvanainen J (2001) Differences in host use efficiency of larvae of a generalist moth,
543
Operophtera brumata on three chemically divergent Salix species. J. Chem. Ecol. 27(8):1595–1615.
544
doi:1023/A:1010458208335
545
Salazar D, Jaramillo MA, Marquis RJ (2016) Chemical similarity and local community assembly in the species rich
546
tropical genus Piper. Ecology 97(11):3176–3183. doi:10.1002/ecy.1536
547
Shaw RF, Pakeman RJ, Young MR, Iason GR (2010) Microsite affects willow sapling recovery from bank vole
548
(Myodes glareolus) herbivory, but does not affect grazing risk. Ann. Bot. 112:731–739. doi:10.1093/aob/mct126
549
Sipura M, Tahvanainen J (2000) Shading enhances the quality of willow leaves to leaf beetles – but does it matter?
550
Oikos 91:550–558. doi:10.1034/j.1600-0706.2000.910317.x
551
Skvortsov AK (1999) Willows of Russia and adjacent countries. Taxonomical and Geographical Revision. Joensuun
552
Yliopistopaino, Joensuu, Finland
553
Tahvanainen J, Helle E, Julkunen-Tiitto R, Lavola A (1985) Phenolic compounds of willow bark as deterrents against
554
feeding by mountain hare. Oecologia 65:319–323
555
Tegelberg R, Veteli T, Aphalo PJ, Julkunen-Tiitto R (2003) Clonal differences in growth and phenolics of willows
556
exposed to elevated ultraviolet-B radiation. Basic Appl. Ecol. 4:219–228
557
Volf M, Hrcek J, Julkunen-Tiitto R, Novotny V (2015) To each its own; differential response of specialist and
558
generalist herbivores to plant defense in willows. J. Anim. Ecol. 84:1123–1132. doi:10.1111/1365-2656
559
Volf M, Segar ST, Miller SE, Isua B, Sisol M, Aubona G, Šimek P, Moos M, Laitila J, Kim J, Zima Jr J, Rota J,
560
Weiblen GD, Wossa S, Salminen J-P, Basset Y, Novotny V (2018) Community structure of insect herbivores is driven
561
by conservatism, escalation and divergence of defensive traits in Ficus. Ecol. Lett. 21:83–92. Doi:10.1111/ele.12875
562
Wu J, Nyman T, Wang D-C, Argus GW, Yang Y-P, Chen J-H (2015) Phylogeny of Salix subgenus Salix s.l.
563
(Salicaceae): demilitation, biogeography, and reticulate evolution. BMC Evol. Biol. 15:31. doi:10.1186/s12862-015-
564
0311-7