Introduction

Light is one of the most important factors in plant growth and survival and has a direct influence on the level of many metabolites. The availability of light influences the amount of assimilated carbon and affects its allocation (Dudt and Shure 1994; Poorter et al. 2012). Carbon assimilated by tree seedlings during photosynthesis can be allocated to three main processes: growth (including reproduction), defense, and storage. Shade tolerance can determine the relative allocation of carbon among these three processes. It can also be assumed that the patterns of carbon allocation are optimized to maximize survival and growth (Imaji and Seiwa 2010). Compared with light-demanding species, shade-tolerant species allocate more carbon to defense and storage and less to growth. Many researchers have suggested the widespread existence of a trade-off between defense and growth (Zangerl and Bazzaz 1992; Haukioja 2003; Villar et al. 2006), but others have questioned the universality of such observations (Riipi et al. 2002; Gaylord et al. 2007). Koricheva (2002), using meta-analysis, described several potential sources of variation in the cost of exhibiting defense mechanisms, such as the presence or absence of herbivores, low or high resource availability, growth rate, and types of defense. She asserted that the failure to detect allocation trade-offs may be a consequence of the particular experimental designs used, and suggested that future studies should examine various aspects of the plant environment: level of competition, herbivore densities, or abiotic conditions.

All plants, including tree and shrub species, can be categorized according to their light requirements. Shade tolerance is usually defined as the capacity to grow in the shade (Niinemets and Valladares 2006), and is measured by the average lowest light levels at which a taxon can survive in a forest understory (Valladares and Niinemets 2008). It is assumed that light-demanding species (which are often early-successional) invest the majority of assimilated carbon in growth (Coley et al. 1985). This growth occurs at the expense of the production of secondary metabolites, which largely consist of compounds active in defense against herbivores and pathogens. Some light-demanding species, losing as many as half of their leaves to herbivory may not impact their capacity for CO2 assimilation because the remaining leaves are better insulated and better supplied with water and nutrients (Heichel and Turner 1983). In the case of trees, this means that the production of wood may not be reduced (Hoogesteger and Karlsson 1992). In contrast to light-demanding species, shade-tolerant species (which are often late-successional) invest a substantial part of their carbon resources in chemical defenses (Coley 1987), such as those provided by phenolic acids, flavonoids, and tannins (Bennett and Wallsgrove 1994). There are strong relationships between the content of phenolic compounds and resistance to insects (Barton and Koricheva 2010). Phenolic compounds are often negatively correlated with nitrogen concentration, and a decrease of phenolics is often observed after nitrogen fertilization (Hakulinen et al. 1995). This is explained by the carbon/nutrient balance hypothesis (Bryant et al. 1983), which states that phenolic compounds, as metabolites containing more carbon, are produced by plants under nutrient-poor conditions, mainly limited by the low level of nitrogen in soil. Under nutrient-rich conditions, carbon is used for the synthesis of proteins and other structural or energy-related compounds. If this hypothesis holds true, we can assume that the sensitivity of plants to herbivores may be connected to the benefits associated with leaves accumulating higher levels of proteins and carbohydrates alongside a lower level of phenolics (Mooney et al. 2009; Wright et al. 2010).

In terms of the association between the carbon/nutrient balance hypothesis and the use of different defense strategies, it is assumed that light-demanding species in contrast to shade-tolerant ones, defend less their leaves (Coley et al. 1985). Indeed, significant differences in the levels of antiherbivore metabolites were observed between birch and alder grazed by the alder beetle, Agelastica alni (Oleksyn et al. 1998). In the leaves of grey alder, the concentration of phenolic compounds was found to increase with increasing perforation of leaves caused by alder beetle. However, in the leaves of birch, this correlation was not observed. In our previous study, we observed that the time required for half of the leaves of a birch subjected to artificial perforation to fall was about 4 weeks earlier than that for control leaves. In the case of alder, this difference was limited to 1 week (Giertych et al. 2006). These observations support the idea that light-demanding and fast-growing species may have less developed leaf defenses.

Against this background, the objective of the present study was to examine the general hypothesis that light-demanding tree species differ from shade-tolerant ones in terms of their carbon allocation. A more specific hypothesis about differences in defense strategies between fast- and slow-growing species was proposed by Coley (1987) after investigation of tropical plants. This proposes that fast-growing species are usually light-demanding and do not invest heavily in chemical defense, whereas slow-growing species (which are often shade-tolerant) defend their leaves by producing secondary compounds.

In this study, we specifically ask the following questions: (1) To what extent do light conditions during growth and species’ capacity for shade tolerance alter the allocation of mass to different organs? (2) Does the level of carbon allocation to growth, defense, and storage vary depending on species’ shade tolerance, and is it affected by the light conditions during growth?

Materials and methods

Light demands of the species studied

The species used in this study differ in terms of their light demands (Table 1). We used indicators of shade tolerance (Ellenberg et al. 1991; Zarzycki 1984) and the tolerance scales range constructed by Niinemets and Valladares (2006) to select tree species. We chose two shade-tolerant species (Fagus sylvatica and Acer platanoides), two intermediate ones (Carpinus betulus and Quercus robur), and two light-demanding ones (Sorbus aucuparia and Betula pendula). In order to define an objective measure of shade impact, the ratio of the mean mass of seedlings from shade to the mean mass of seedlings from full light—shade mass effect was calculated (Table 1).

Table 1 Comparison of ecological light indicators for six examined tree species

Cultivation of seedlings

One-year-old silver birch (Betula pendula Roth, N = 60), European mountain ash (Sorbus aucuparia L., N = 32), pedunculate oak (Quercus robur L., N = 28), Norway maple (Acer platanoides L., N = 28), common hornbeam (Carpinus betulus L., N = 60), and European beech (Fagus sylvatica L., N = 28) seedlings were grown individually in 2000-cm3 pots filled with a 1:1 (v/v) mixture of forest soil and peat. Seedlings were obtained from commercial nursery (birch and hornbeam) or cultured from seeds in full-sunlight conditions. Before the experiments, seedlings were replanted to the pots and fertilized (2 kg/1000 dm3) with granulate Osmocote slow-release (5–6 months) fertilizer containing: 15 % N (7 % NO3 and 8 % NH4 +), 9 % P (P2O5), 12 % K (K2O), and 2.5 % Mg (MgO). At the beginning of the growing season (in April), the height and diameter of each seedling were measured (the initial height was later used as a covariant during the statistical analyses of growth parameters), and the seedlings of each species were randomly divided into two groups. The first group was placed into two identical shade houses covered with black polypropylene shade cloth, which allowed the transmittance of 5 % light (shade treatment). The second group of seedlings was grown in two plots with full exposure to sunlight. The sites where all treatments and blocks were carried out were no more than 10 m apart. Temperature and relative humidity during the experiments were monitored using HOBO data loggers (HOBO H08-032-08, Onset Computers, USA). Based the monthly average values, it was found that shade decreased temperature of about 2 °C but increased the relative humidity of dozen percent. This work was carried out at the experimental field of the Institute of Dendrology in Kórnik, Poland (52°14′N; 17°05′E; 75 m altitude). During growing season, the seedlings were watered as necessary and all weeds were periodically removed.

Growth parameters, morphological, and chemical analyses

For each seedling, the height and diameter above the root collar were measured at the beginning of the experiment, and three times during the growing season. In October, toward the end of the growing season, eight seedlings of each species from the two light treatments were randomly chosen for morphological and chemical analyses. In a few cases, owing to a small amount of available material, some analyses were carried out with fewer replications. The roots were rinsed, and each seedling was divided into four parts: coarse roots (diameter >2 mm), fine roots (diameter <2 mm), stem (with branches), and leaves. Subsamples of leaves and fine roots from each seedling were used to determine projected leaf area (PLA) using WinFolia Pro, and total surface root area (TRA) and fine root length using both WinRhizo softwares (Regent Instruments, Inc., Quebec, Canada). Specific leaf area (SLA: defined as the PLA divided by the leaf dry mass, cm2 g−1) and the total root area and length for each seedling were calculated. All parts of the seedlings were oven-dried (65 °C for 48 h), ground to a fine powder using a Mikro-Feinmühle Culatti mill (IKA Labortechnik, Staufen, Germany), and stored in plastic boxes. Descriptions and abbreviations of all measured parameters are given in Table 2.

Table 2 Parameters of plant morphology and chemistry, with their abbreviations and units of measurement

Nitrogen (N) and carbon (C) concentrations in the leaves were measured using the Elemental Combustion System CHNS-O (Costech Analytical Technologies Inc., Valentia, USA). The concentration of phenolic compounds (TPh) was measured colorimetrically using Folin and Ciocalteu’s Phenol Reagent (Sigma F-9252), following Johnson and Schaal (1957), as modified by Singleton and Rossi (1965). The content of total phenols was expressed as μmol chlorogenic acid per gram of dry mass. Total soluble carbohydrate and starch concentrations were determined by a modified version of the methods described by Hansen and Møller (1975) and Haissig and Dickson (1979). Carbohydrates were extracted from tissue powder in methanol–chloroform–water (12:5:3,v:v:v), and the tissue residue was used to determine the starch content. Concentrations of soluble carbohydrates and starch (with glucose as a standard) are expressed as  % of glucose per gram of dry mass.

It is clear that seedling mass depends on the amount of accumulated carbon. Accordingly, in order to clarify the distribution of carbon allocation, we expressed some of our results as the ratios of growth parameters or chemical compositions to seedling masses.

Statistical analyses

After initial testing of residuals normality (Shapiro–Wilk test) and equal variance (Levene test), analysis of covariance (ANCOVA) was used to evaluate the influence of species, light treatment, and their interaction on growth and morphological parameters. The initial height was used as a covariant (except for the final diameter, for which the initial diameter was used). Analysis of variance (ANOVA) was used to assess the influence of species, light treatment, and their interaction on allocation ratios, leaf morphology parameters, contents of leaf TPh, carbohydrate, N, and C. The results expressed in percent were arcsin-transformed for normality for ANOVA analyses. Relations between the shade tolerance (shade mass effect) and growth, morphological, and chemical parameters were expressed as a linear regression. Shade mass effect was used as independent variable in the regression analysis. All analyses were performed using JMP software (version 7.0.2; SAS Institute, Cary, NC, USA).

Results

Influence of light conditions and shade tolerance on growth parameters

The variability of reactions in examined species was high, and almost all tested morphological and chemical parameters demonstrated a significant statistical difference (Tables 3, 4, 5). The shade-tolerant species were characterized by smaller differences in entire-plant mass between shade and light treatment (higher values of shade mass effect). Regardless of light condition, the shade mass effect correlated negatively with aboveground to belowground mass ratio (Fig. 1). Light-demanding species allocate more biomass to aboveground parts of seedlings, and shade conditions even enhance the relationship.

Table 3 Summary of ANCOVA for seedling height, seedling diameter, seedling dry mass, stem dry mass, leaf dry mass, root dry mass, coarse root dry mass, fine root dry mass, fine root length, and surface root area with light conditions of growth, species, and interactions
Table 4 Summary of ANOVA for allocation to stem (SMF), allocation to leaves (LMF), allocation to roots (RMF), specific leaf area (SLA), mean leaf mass, mean leaf area, stem carbon, stem nitrogen, and stem C/N ratio with light conditions of growth, species, and interactions
Table 5 Summary of ANOVA for leaf soluble carbohydrates, leaf starch, root soluble carbohydrates root starch, leaf carbon, leaf nitrogen, leaf C/N ratio, leaf TPh, root carbon, root nitrogen, and root C/N ratio with species, light conditions of growth and interactions
Fig. 1
figure 1

Relationship between shade mass effect and aboveground to underground mass ratio for six species: asterisks Fagus sylvatica; box Acer platanoides; triangle Quercus robur; diamond Carpinus betulus; circle Betula pendula; square Sorbus aucuparia growing under shade (filled symbol) and full-sunlight (open symbol) conditions

Exposure to shaded conditions clearly reduced the growth of all species (Tables 3, 6). The extent of the reduction was correlated with the shade tolerance of the species, and was greater in light-demanding species. This became the basis for calculating shade mass effect. The light conditions clearly influenced seedling height, diameter, and mass. Among all of the growth-related traits, only plant height of the shade-tolerant species in shade conditions was similar or greater than that of the same species growing under full sunlight. Other growth-related traits (diameters above the root collars, and masses of leaves, stems, total roots, coarse roots, and fine roots) in all species showed greater values under full-sunlight conditions than under shade conditions (the exception is the greater mass of F. sylvatica leaves under shade). There were significant interactions between light conditions and tree species for all these traits (Table 3). It means that the shade conditions resulted reduction of the growth-related traits, particularly in the light-demanding species (Table 6; see Tukey test results). The reduction of mass associated with shade conditions was organ-specific, and was the greatest for roots and the smallest for leaves. The light conditions influenced mass allocation to stem (SMF), roots (RMF), and leaves (LMF), with shade-tolerant species allocating mass predominantly roots. Interactions between light conditions and tree species in terms of biomass allocation to different organs in seedlings were not statistically significant (Table 4).

Table 6 Mean (SD) of biometric parameters depending on species and light conditions

The light conditions and tree species affected almost all leaf parameters in all species. The seedlings grown in the shade were characterized by lower leaf mass (except F. sylvatica) and greater SLA and mean PLA values than those of plants grown in full sunlight (Tables 6, 7). There were also significant interactions between tree species and light conditions for SLA (Table 4). The largest increase of SLA and the largest decrease of leaf mass under shade conditions were exhibited by the light-demanding species (Table 7; see Tukey test results). The light-demanding species scored the largest leaf area ratio (LAR) under light conditions (23,7 and 79,3 for S. aucuparia and B. pendula, respectively) and the largest increase of LAR under shade conditions (221.7 and 216.5 for S. aucuparia and B. pendula, respectively). The relationship between shade mass effect and LAR was significant only in the case of shade condition (Fig. 2). Shade causes a high significant decrease in all measured root parameters such as total root mass, coarse root mass, fine root mass, and surface area (Tables 3, 6).

Table 7 Mean (SD) of allocation, leaf parameters, and stem chemical compositions depending on species and light conditions
Fig. 2
figure 2

Relationship between shade mass effect and leaf area ratio (cm2 leaves g−1 plant) for six species: asterisks Fagus sylvatica; box Acer platanoides; triangle Quercus robur; diamond Carpinus betulus; circle Betula pendula; square Sorbus aucuparia growing under shade (filled symbol) and full-sunlight (open symbol) conditions

Influence of light conditions and shade tolerance on chemical composition

The light conditions and tree species significantly influenced all the investigated chemical elements and compounds in leaves, stems, and roots. The leaves of species subjected to full-sunlight conditions had less nitrogen and more carbon, phenolic compounds and a higher carbon/nitrogen (C/N) ratio. In most cases, it was also a significant interaction between light conditions and species. The reason is that the observed differences between the contents of the examined compounds and elements were greater in the light-demanding species (Tables 5, 8).

Table 8 Mean (SD) of leaf and root chemical compositions depending on species and light conditions

The efficiency of carbon assimilation, when expressed as a ratio of carbon in the entire plant to leaf mass, showed significant dependence from shade tolerance (Fig. 3). This relationship is very similar for both light conditions. The highest efficiency of carbon assimilation showed a shade-tolerant species (F. sylvatica and A. platanoides) and the lowest light-demanding (B. pendula and S. aucuparia).

Fig. 3
figure 3

Relationship between shade mass effect and the ratio carbon content in the entire plant to leaf mass for six species: asterisks Fagus sylvatica; box Acer platanoides; ttriangle Quercus robur; diamond Carpinus betulus; circle Betula pendula; square Sorbus aucuparia growing under shade (filled symbol) and full-sunlight (open symbol) conditions

For stems, the contents of carbon and nitrogen, as well as the C/N ratio, depended on light conditions and on species (Table 4). The light conditions also altered the root contents of the examined chemical elements and compounds (Tables 5, 7, 8). The results were similar to those in leaves; however, the content of soluble carbohydrates was lower and that of starch was higher in roots compared with those in leaves.

The relationships between stem height increment and diameter increment to seedling mass, all of which may reflect how carbon is invested by a plant, were found to depend on the light conditions and the species’ tolerance of shade (Fig. 4). Under full-sunlight conditions, these relationships were similar for all groups of examined species. However, under shade conditions, the increases in the ratio height increment to seedling mass were the highest in the two light-demanding species (Fig. 4). Although the correlation coefficient was not significant, a clear trend is well noticed. The ratio of stem diameter increment to seedling mass for the shade-tolerant species was lower under shade conditions compared to the full sunlight. The relationships of carbohydrate concentrations in roots or leaves and total phenolic concentration in leaves to seedling mass also depended on species’ shade tolerance (shade mass effect) and light conditions during growth (Fig. 5). For full-sunlight conditions, these relationships were significant, light-demanding species under full sunlight invested less carbon to phenolic compounds and carbohydrates. Shade conditions increased the ratios of carbohydrate and phenolic contents to seedling mass in light-demanding species. For shade-tolerant species, shade conditions increased the ratios of carbohydrate to seedling mass, but decreased the ratio phenolic content to seedling mass (Fig. 5).

Fig. 4
figure 4

Relationship between shade mass effect and the ratio stem height increment to seedlings mass and stem radial increment to seedlings mass for six species: asterisks Fagus sylvatica; box Acer platanoides; ttriangle Quercus robur; diamond Carpinus betulus; circle Betula pendula; square Sorbus aucuparia growing under shade (filled symbol) and full-sunlight (open symbol) conditions

Fig. 5
figure 5

Relationship between shade mass effect and leaf soluble phenolics (TPh), leaf soluble carbohydrates (SC), root soluble carbohydrates (SC) to seedlings mass for six species: asterisks Fagus sylvatica; box Acer platanoides; triangle Quercus robur; diamond Carpinus betulus; circle Betula pendula; square Sorbus aucuparia growing under shade (filled symbol) and full-sunlight (open symbol) conditions

Discussion

For the seedlings of forest trees, during their first few years of life, development very often occurs in deep shade and shade tolerance plays a key role in the seedlings survival (Walters and Reich 1999). Our results showed that shade tolerance influences carbon allocation. The higher ratios of stem height increment and stem radial increment to seedling mass in light-demanding species (Fig. 4) indicate that these species invest proportionally more of their carbon into growth compared with shade-tolerant species. Higher allocation of carbon to growth processes in the case of light-demanding species was particularly pronounced under shade conditions. Niinemets (1998) maintains that, under conditions of low irradiance, shade-tolerant Acer platanoides and the intermediate species Quercus robur allocate more resources for promoting an increase in height. In our experiment, the mean heights (Table 6) of both of these species were also greater in the shade than under the full-sunlight conditions. This is connected with the so-called “lammas shoots,” which were observed in some seedlings growing under shade conditions and increased mean height relative to that of seedlings grown under high light intensities. Moreover, the pattern of carbon allocated to increase height in shade-tolerant and intermediate species was less dependent on light conditions than that in light-demanding species.

Rapid shoot elongation is an adaptation of light-demanding species to achieve high growth quickly when intense light conditions appear (Imaji and Seiwa 2010). In natural forests, conditions of intense light are found under gaps in the canopy, which are usually short-lived in mature forests (Madsen and Hahn 2008). Therefore, the seedlings of light-demanding species only have a few years to reach the forest canopy; otherwise, at strata nearer the ground, they are outcompeted by slower growing shade-tolerant species. Light-demanding species were previously shown to allocate more mass to the stem and less to the roots (Paz 2003), which was similar to the results in our experiment. Shade conditions increased the allocation to aboveground part of plant and decreased that of root (RMF), regardless of the shade tolerance of the species. However, this increase is higher in light-demanding species. Despite the increase of the LMF in the majority of examined species, the seedlings growing in the shade had a greater mean leaf area (Table 7) and greater LAR (Fig. 2). McCarthy and Enquist (2007) suggest that intraspecific variability in biomass allocation is correlated with environmental factors. An increase in leaf area is connected with the need to increase light capture. These results support the optimal partitioning theory (Thornley 1972; Bloom et al. 1985), which suggests that plants allocate biomass to components that will most facilitate the acquisition of those resources that are in limited supply at the time.

The influence of light conditions during growth on the chemical composition of plants was as expected. Under shade conditions, the carbon content in leaves, stem, and roots was lower, regardless of the shade tolerance of species. Walters and Reich (1999) suggested that young seedlings of shade-tolerant species do not have advantages over light-demanding species in terms of the efficiency of capture of carbon under low-light conditions. However, in our experiment, shade-tolerant species assimilated more carbon per leaf mass unit than intermediate or light-demanding species (Fig. 3). This confirms the suggestions of Teskey and Shrestha (1985) that shade-tolerant species may have more efficient systems for assimilating carbon dioxide.

The nitrogen content in leaves, stem, and roots of seedlings grown in the shade was approximately double than under full-sunlight conditions. A higher content of nitrogen in leaves upon growth in the shade was also reported by Niinemets (1997) and Baraza et al. (2004). However, the influence of shading on the concentration of nitrogen in roots and stem in our experiment differs from the findings from some other studies. Machado and Reich (2006) showed a significant decrease of nitrogen in the stem with increasing shade, but only for one (Pinus strobus) of three examined cold-temperate tree species. Villar-Salvador et al. (2004), who investigated the effect of nitrogen fertilization and shading on the functional quality of Quercus ilex seedlings, also found a significantly lower nitrogen contents in root and stem in plants grown under shade conditions than in those grown under full sunlight. These differences can be explained by the greater size of a plant upon growth in full-sunlight conditions, such that the available nitrogen is diluted across the greater mass of the plant.

According to the hypothesis of Coley (1987), shade-tolerant species allocate more carbon to metabolites involved in defense against pathogens and herbivores. In our experiment, the level of phenolic compounds in shade-tolerant species was also higher than that in light-demanding species, regardless of growth conditions. Shade-tolerant plants also have a longer leaf lifespan (Hallik et al. 2009); therefore, they are forced to establish defense systems for their leaves, especially in the case of plants growing in shaded conditions, which are usually more attractive to herbivores (Henriksson et al. 2003;). Other studies showed that shade-tolerant tree species have higher levels of defense metabolites than light-demanding species (Dudt and Shure 1994; Imaji and Seiwa 2010). We also noted the important interaction in phenolic content between light conditions during growth and shade tolerance. In this study, the decrease of leaf phenolic compounds under low-light conditions was greater in light-demanding species than in shade-tolerant ones. However, the amount of carbon allocated to phenolic compounds (phenolic compound content/seedling mass) in light-demanding species was significantly higher under shade conditions. In shade-tolerant species, the amount of carbon invested in defense metabolites (phenolics) was greater under full-sunlight conditions. This can be explained by the role of phenolics in protection against photodamage (Close and McArthur 2002). Despite the fact that the efficiency of carbon assimilation, expressed as the ratio of carbon in the entire plant to the leaf mass in shade-tolerant species, depends from shade mass effect (Fig. 3), shade-tolerant species contained more phenolic compounds and they allocated more carbon to defense than light-demanding or intermediate species (Fig. 5).

Our results suggest that the level of carbon allocation to carbohydrates differs between the various plant organs and light conditions. Higher ratios of leaf soluble carbohydrates/seedling mass and root soluble carbohydrates/seedling mass under shade conditions may reflect the lack of a trade-off relationship between storage and growth or storage and defense. According to Kobe (1997), allocation of carbon to stored carbohydrate may be an effective strategy in shade tolerance; however, many other researchers dispute the hypothesis (Machado and Reich 2006; Piper et al. 2009). A higher level of allocation of carbon in the form of carbohydrates to roots by shade-tolerant species under full-sunlight conditions and to leaves by light-demanding species under shade conditions (Table 8) may be characteristic only for the end of the growing season, that is, when the studies were carried out. There is thus a need to investigate these effects during the whole growing season.

In conclusion, this study demonstrated that the shade tolerance of tree seedlings influences carbon allocation. When expressed as the ratio of carbon in the entire plant to leaf mass, the efficiency of carbon assimilation correlated positively with shade mass effect (Fig. 3). In the case of shade-tolerant species, the findings suggest the existence of a trade-off mechanism. Under shade conditions, the shade-tolerant species restrict the allocation of carbon to roots, defense mechanisms (such as the accumulation of phenolics), and radial growth. Carbon assimilated by shade-tolerant seedlings of deciduous trees under shade conditions is instead allocated much more to the accumulation of storage compounds (nonstructural carbohydrates) and the promotion of plant height. In contrast, in light-demanding seedlings of deciduous trees grown under shade conditions, carbon is allocated not only to structural tissues, enabling a high level of growth, but also to nonstructural carbohydrates in leaves or roots. However, in the case of these light-demanding species, we did not note a trade-off mechanism. Under stress conditions (shade), the ratios of examined growth parameters and metabolites to seedling mass were always higher than under full-sunlight conditions. Only the allocation to roots was lower.

Author contribution statement

MJG, PK, and JO conceived and designed the experiments. MJG cultured plant materials, collected samples, performed the experiments and the data analysis, and drafted the earlier versions of the manuscript. PK and JO revised the manuscript.