Free Access
Issue
Ann. For. Sci.
Volume 66, Number 7, October-November 2009
Article Number 709
Number of page(s) 12
DOI https://doi.org/10.1051/forest/2009062
Published online 02 October 2009

© INRA, EDP Sciences, 2009

1. INTRODUCTION

During ontogeny, trees increase both in size and structural complexity. It is now well established that leaf and whole plant traits change with tree size (mainly height) and age (e.g. Hoflacher and Bauer, 1982; Mediavilla and Escudero, 2003; Rijkers et al., 2000; Thomas and Winner, 2002). Tree height is known to influence structural and physiological leaf traits possibly due to hydraulic constraints imposed by increased transport distance (Niinemets, 1997a; 2002; Rijkers et al., 2000). In several tropical tree species, leaf mass-per-area ratio (LMA), area-based total leaf N content ( Na) and area-based light-saturated net CO2 assimilation rate (Asat a) increase with tree height whereas mass-based total leaf N (Nm), mass-based light saturated net assimilation rate (Asat m) and photosynthetic N-use efficiency remain stable (Rijkers et al., 2000).

On the other hand, relative irradiance is well known to dramatically modify leaf traits, eliciting a large phenotypic plasticity (Bradshaw, 1965). In general, shade leaves are thinner, less dense and display lower leaf mass-to-area ratios (LMA) than sun-exposed leaves (Givnish, 1988; Rozendaal et al., 2006; Valladares et al., 2000). They also exhibit lower Asat a and dark respiration per leaf area (Rd a) (Givnish, 1988; Rijkers et al., 2000). However, little is known about the modulation of irradiance-elicited plasticity by ontogeny. Many authors emphasized that, because traits vary during the course of growth and development, it is necessary to compare individuals at a common ontogenetic development stage (i.e., a common size rather than a common age) to identify environmentally induced variation in traits (Coleman et al., 1994; Poorter, 1999; Poorter et al., 2006; Veneklaas and Poorter, 1988; Wright and McConnaughay, 2002). Furthermore, growth and development rates themselves are highly plastic. Individuals grown under contrasted irradiance usually display different growth rates, and therefore reach largely different sizes at similar ages. Plasticity may therefore be confounded with ontogenetic responses (Veneklaas and Poorter, 1988). This was particularly emphasised in studies of shade tolerance (Lusk and Warton, 2007; Lusk et al., 2008; McConnaughay and Coleman, 1999; Niinemets, 2006).

The majority of studies addressing ontogenetic shifts of leaf traits used plant size (height, diameter or biomass) as a covariant (Poorter, 1999; Sterck and Bongers, 1998) or focused on a comparison between very young seedlings or saplings and mature trees and did not take into account any intermediate stage (Cavender-Bares and Bazzaz, 2000; Cornelissen et al., 2003; Holscher, 2004; Mediavilla and Escudero, 2003; Thomas and Winner, 2002). A meta-analysis of field measurements in saplings vs. mature trees (Thomas and Winner, 2002) showed that Asat a was larger in adult trees than in saplings, especially for temperate deciduous trees. An ontogenetic increase in LMA led to a decrease of Asat m (Mediavilla and Escudero, 2003; Roggy et al., 2005; Thomas and Winner, 2002).

The observed ontogenetic effects may also be induced by differences in tree architecture: an increased structural complexity, due to branching, crown development or occurrence of flowering may have important physiological consequences on functional leaf traits. Successive key stages of architectural development (ASD) have been identified (Barthelemy and Caraglio, 2007), on the basis of the degree of structural complexity and of the expression of processes like branching, reiteration or sexuality. The successive stages may be reached at different dimensions for individual trees, depending on environmental conditions (Nicolini, 1998; Nicolini et al., 2000; Roggy et al., 2005).

The main objectives of this study were to document these questions based on observations with young individuals from two tree species growing along a natural irradiance gradient in the tropical rain forest. Dicorynia guianensis Amshoff and Tachigali melinonii (Harms) Barneby (Fabaceae, Caesalpinioideae) were selected due to their ability to establish both in the understorey and in clearings. The two species belong to the same functional group (heliophilic non pioneer species, Molino and Sabatier, 2001). More specifically, our aim was: (1) to test whether structural and functional leaf traits differed among ASDs in these species; (2) to quantify the plasticity of leaf traits along a natural irradiance gradient and (3) to evidence potential interactions between irradiance and ASD, i.e., to test whether irradiance-elicited plasticity differed between ASDs.

In the present study, we assessed the variability of eleven morphological and physiological leaf traits related to photosynthesis over three different ASDs in forests around Kourou, French Guiana. Leaf structure was assessed from leaf mass-to-area ratio (LMA), leaf thickness, leaf density and nitrogen, carbon and chlorophyll contents. The semi-mechanistic photosynthesis model of Farquhar et al. (1980) was used to quantify photosynthetic capacity by assessing in vivo apparent maximal Rubisco activity (Vcmax), maximal light-driven electron flux (Jmax) and dark respiration (Rd). Photosynthetic N use efficiency was indirectly assessed by the ratios Vcmax /N; Jmax /N and the relative allocation of N to light interception was assessed from the ratio Chlorophylls/N.

2. MATERIAL AND METHODS

2.1. Study site

The study was carried out in coastal French Guiana in South America, in the lowland tropical rain-forest of the experimental site of Paracou and in nearby sites (5° 18′ N, 52° 55′ W). The climate is characterised by a wet season from December to August, interrupted during February by a short dry season, and followed by a long dry season from September to November. Average annual precipitation is 2875± 540 mm (mean ± SD from 1986 to 2005) with nearly two-third received between mid-March and mid-June (Baraloto et al., 2007; Bonal et al., 2000).

thumbnail Figure 1

Stages of architectural development (ASD) of Dicorynia guianensis (DG) and Tachigali melinonii (TM): ASD1: saplings with un-branched main stem; ASD2: saplings with sparsely branched main stem (order 2 axes remain un-branched); ASD3: immature trees (diamonds: sequential axes; crosses: reiterated axes). For a complete description of the different ASDs, see Material and methods.

2.2. Sampling

Dicorynia guianensis Amshoff and Tachigali melinonii (Harms) Barneby (Fabaceae, Caesalpinioideae) are two rather heliophilic but non-pioneer species with compound leaves (Molino and Sabatier, 2001). D. guianensis is generally considered to be relatively tolerant to shade (Favrichon, 1994; Roggy et al., 2005). Individuals were selected at three different stages of architectural development (ASD): (1) saplings with un-branched main stem (ASD 1); (2) saplings with a sparsely branched main stem (order 2 axes remain unbranched, ASD 2) and (3) immature trees with an abundantly and regularly branched main stem (ASD 3, Fig. 1). In the latter case, secondary axes are present; this stage follows the “architectural metamorphosis” during which the tree evolves from a system organised around a single main axis to a system organised around several major branches (Barthelemy and Caraglio, 2007). For more details, see descriptions in Roggy et al. (2005). 23 individuals were sampled for D. guianensis and 47 for T. melinonii. D. guianensis reaching ASD3 were only found under high irradiance and not under the full range of irradiance levels. Measurements were made from July to August 2000 for D. guianensis and during March and April 2003 for T. melinonii. The dimensions of individuals at the different ASDs differed severely between the two species, and in general D. guianensis displayed a larger height at any given ASD than T. melinonii (Fig. 2). For the tallest individuals, small scaffolds were used to access top canopy leaves from the last expanded cohort.

Daily photon flux (DPF) was measured with amorphous silicon quantum sensors for photosynthetically active radiation (PAR CBE 80 Solems S.A., Palaiseau, France). DPF was measured close to each sapling (or group of saplings) at the level of the measured leaf, and compared to a reference sensor located in a clearing; this reference sensor produced an estimate of the incident irradiance above the canopy. Measurements were made over 30 consecutive days from July 20 to August 18, 2000 for D. guianensis and over at least 3 d per individual for T. melinonii during March, April 2003 and January 2004. Values of DPF ranged from 0.05 to 36.7 mol m−2 d−1 among individuals. For each individual, relative DPF ( DPFrel , % of full irradiance) was computed as the ratio of the average local and above canopy DPF.

thumbnail Figure 2

Effects of successive stages of architectural development (ASD 1 to 3) on height (H, cm) for T. melinonii (white) and D. guianensis (black). Values are least square means (± standard error) adjusted to a common DPFrel of 26.1 %. Differences were tested with the multiple comparisons post-hoc test of Tukey (P < 0.05) and significant differences are indicated by different letters.

Table I

Primary parameters of the photosynthesis model used in this work. For details, see Dreyer et al. (2001).

2.3. Gas exchange measurements

Leaf gas exchange was measured with a portable gas exchange system (CIRAS-1, PP-Systems, Hitchin, U.K.) equipped with a Parkinson leaf chamber (2.5 cm2). All measurements were performed on young, fully expanded leaves selected on top of the main axis, usually on the second-last growth unit of the main shoot. A halogen lamp (Philips 12 V, 20 W) yielded a saturating photosynthetic photon flux density (PPFD) of 700 μ mol m−2 s−1. Air temperature in the leaf chamber and water vapor pressure deficit at the leaf surface were kept as stable as possible, but nevertheless varied significantly because of changes in the external microclimate (means ± SD, 29± 1.5°C and 1.6 ± 0.2 kPa, respectively). To compensate for differences in temperature during measurements, all parameters were scaled to a common reference temperature (Tref) of 25 °C (298.15 K). Because of the lack of data for tropical tree species, we used the temperature responses computed by Dreyer et al. (2001) for seedlings of temperate forest species (Tab. I). For each leaf, a full A - Ci curve was recorded, after induction under saturating irradiance for at least 20 min, by stepwise changes of CO2 mole fraction in the air between 1800 and 50 μ mol mol −1 (5-min stabilization at each step and three successive measurements at 1-min intervals to check for stability). At least 12 measurement steps were used for each response curve. Dark respiration, Rd, was recorded after 5 min in the dark by switching-off illumination. No correction was introduced for potential leaks in the gaskets of the chamber.

2.4. Estimating photosynthetic capacity

We used the Harley et al. (1992) version of Farquhar’s biochemical model of leaf photosynthesis (Farquhar et al., 1980) in which net CO2 assimilation rate ( A; μ mol CO2 m−2 s−1) is expressed as:

(1)
where Vc is Ribulose-bisphosphate (RuBP) carboxylation rate, Vo is RuBP oxygenation rate and Rd corresponds to CO2 production by non-photorespiratory respiration in the mitochondria. A - Ci curves were used to estimate: (1) Vcmax, the maximal rate of RuBP carboxylation, i.e., the potential activity of Rubisco in the absence of limitation by CO2 availability and photorespiration and (2) Jmax, the maximal light-driven electron flux. For a detailed description of the equation set, see Le Roux et al. (1999). Primary model parameters and their temperature dependencies are listed in Table I. The rate limitation by triose-phosphate availability was not estimated. Internal conductance to CO2 transfer, gi (Ethier and Livingston, 2004; Warren, 2006) was assumed to be infinite and the reported values of Vcmax (and to a lesser extent Jmax) should therefore be considered as apparent- Vcmax including a gi component. Real values of Vcmax were probably significantly larger. For each individual, a non-linear least-squares regression (Newton method, ProcNLIN, SAS v.8.1, SAS Institute Inc., Cary, NC) was used to estimate Vcmax and Jmax by fitting the A - Ci curves to model equations as described by Le Roux et al. (1999).

thumbnail Figure 3

Relationships between height (H, cm) and relative daily photon flux (logDPFrel) for the three stages of architectural development (ASD 1, white dots, ASD 2 grey diamonds and ASD 3 black dots) and for two species (T. melinonii and D. guianensis). Regression equations, level of significance (P) and coefficient of determination (r2) are indicated for statistically significant linear regressions. Significance levels: ns, P > 0.05;*, P < 0.05; **, P < 0.01 and ***, P < 0.001.

2.5. Leaf structure and chemistry

After gas exchange measurements, leaves were harvested and their area measured using an electronic area meter (LiCor 3000A, LiCor Inc., Lincoln, NE). Thickness of each T. melinonii leaf was obtained as the mean of three measurements with a digital micrometer (Digimatic micrometer, Mitutoyo, Japan), care being taken to avoid the main veins. For D. guianensis samples, leaf thickness was estimated from anatomical samples (see Roggy et al., 2005). Leaf dry mass was recorded after drying for 48 h at 60 °C. Total N (Nm) and C per unit leaf dry mass ( Cm) were measured with an elemental analyzer (NA 1500 NCS, Carlo Erba, Italy and SCA, CNRS Solaize, France for T. melinonii and D. guianensis, respectively). Values of LMA, Na and leaf density (LMA/thickness) were derived from these data. A portable chlorophyll meter SPAD-502 (Minolta, Osaka, Japan) was used to estimate chlorophyll content on an area basis (Chl a,). A calibration curve relating SPAD measurements and chlorophyll content was established for T. melinonii (Appendix 12). Chlorophyll content was not available for D. guianensis.

Indirect estimates of photosynthetic nitrogen-use efficiency (PNUE) were computed as the ratio between photosynthetic capacity and total leaf N content (in molar units), as Vcmax /N, Jmax /N. Relative allocation of leaf N to light interception was assessed from the ratio Chl/N.

2.6. Statistical analysis

All statistical analyses were performed with STATISTICA (Kernel Version 6.0, StatSoft, Tulsa, OK, USA). Data were log-transformed to respect the required data normality and variance homogeneity, as indicated in Table and Figure legends. One-way ANCOVA, with ASD as main factor and relative irradiance (DPF rel, %) as covariant, was applied for each of the two species. Homogeneity of regression slopes with relative irradiance was tested. When no significant difference was detected, then the interaction term was removed from the model. Inter-specific differences were assessed using adjusted means obtained after elimination of the covariant effect. For D. guianensis, we tested only the effect of ASD 1 and ASD 2. Results of an ANCOVA with ASD 1 and ASD 2 only are displayed in Appendix 22 for T. melinonii. All statistical relationships were considered significant at P < 0.05.

3. RESULTS

3.1. Size and architectural stage of development (ASD)

Height of the saplings ranged from 71 to 1950 cm for Dicorynia guianensis and from 30 to 732 cm for Tachigali melinonii. At all irradiance levels and ASDs, D. guianensis was higher than T. melinonii (Fig. 2). As expected, height (H, cm) and basal diameter (D, cm) increased with successive ASDs. Within a given ASD, larger relative daily photon flux (DPF rel) induced in general a decreased height (Fig. 3) and increased D, resulting in severely declining H / D in the two species (not shown). In ASD 1, the saplings of T. melinonii were very small and their dimension displayed no clear response to DPFrel .

Table II

Results of a one-way Ancova with stage of architectural development (ASD = 3) as main factor and relative irradiance (DPFrel, %) as covariant for Tachigali melinonii. F -values, levels of significance (P) and coefficients of determination (r2) of the model are given. DPFrel, H, D, LMA, Na and Jmax a were log-transformed prior to analysis. Cm and thickness were analyzed with a separate slope model. Vcmax m, Jmax m and Rd m were not linearly correlated with DPFreland thus they were analyzed with a one way ANOVA. Significance levels: ns, P > 0.05;*, P < 0.05; **, P < 0.01 and ***, P < 0.001. The corresponding adjusted means at a 26% DPFrel are displayed in Figure and numerical values are provided in the on-line only Appendix 2.

Table III

Results of a one-way Ancova with stage of architectural development (ASD 1 and ASD 2) as main factor and relative DPF (DPFrel, %) as covariant applied to Dicorynia guianensis. F -value, level of significance (P) and coefficient of determination (r2) of the model are given. DPFrel, H, D and Jmax a were log-transformed prior to analysis. Vcmax a was analyzed with a separate slope model. D, Cm,Jmax a, Jmax/ Vcmax, Vcmax/ N and Jmax/ N were not linearly correlated with DPFreland thus they were analyzed with a one way ANOVA. Significance levels: ns, P > 0.05; *, P < 0.05; **, P < 0.01 and ***, P < 0.001. The corresponding adjusted means at 26% DPFrelare displayed in Figure and numerical values are provided in the on-line only Appendix 2.

3.2. Effect of successive stages of ASD

Mean values of leaf traits adjusted to a common irradiance and for different ASDs are synthesized in Figure 4. In the two species, ASD had a significant impact on several traits, but the impact was larger in T. melinonii than in D. guianensis. This was not due to the fact that we had only 2 ASDs in the latter and 3 in the former case. Differences observed among ASDs were maintained in T. melinonii when ASD 3 was removed from the analysis (Appendix 32).

thumbnail Figure 4

Effects of stages of architectural development (ASD 1 to 3) on leaf traits for T. melinonii (white) and D. guianensis (black). Values are least square means (± Standard error) adjusted at a common relative daily photon flux (DPFrel=26.1 %). Significant differences were tested with the multiple comparisons post-hoc test of Tukey (P < 0.05) and indicated by different letters (a, b, c for T. melinonii and a, bc for D. guianensis).

thumbnail Figure 5

Relationships between leaf traits and relative daily photon flux (logDPFrel) for T. melinonii (white dots) and D. guianensis (black dots). Simple linear regressions were represented by solid lines for T. melinonii and by dotted lines for D. guianensis. Regression equation and r2 are indicated for statistically significant linear regressions (P < 0.05). Significance levels: ns, P > 0.05; *, P 60.05; **, P 60.01; *** and P 60.001.

In the case of T. melinonii, LMA, density and thickness increased from ASD 1 to 3 as did mass-based nitrogen content ( Nm). As a consequence, area-based photosynthetic capacity (Vcmax a and Jmax a) and respiration (Rd a) increased from ASD 1 to 3, as did, albeit to a lesser extent, mass-based photosynthetic capacity (Vcmax m and Jmax m). Photosynthetic nitrogen-use efficiency (PNUE, i.e., Vcmax /N and Jmax /N) was not affected by ASD, which shows that the increase of photosynthetic capacity was driven by increased N availability on a leaf area basis and not by a shift of N allocation towards photosynthesis (Tab. II, Fig. 4).

In the case of Dicorynia guianensis, leaf thickness, density, total leaf N (Na and Nm) and leaf C content (Cm) remained stable in the two ASDs. The ratio Jmax / Vcmax increased from ASD 1 to ASD 2, reflecting a slight increase of Jmax with a constant Vcmax. Photosynthetic N use efficiency (Vcmax /N and Jmax /N) remained stable (Tab. III, Fig. 4). The ASD induced variability was smaller than in T. melinonii and affect less variables.

3.3. Interactions between relative daily photon flux (DPFrel) and ASD

The interactions between DPF rel and ASDs were analysed separately in the two species with an ANCOVA assessing the effect of ASD (factor) and DPF rel

(covariant; Tab. II and III). For most traits, the impact of ASD was much smaller than that of DPF rel. The small level of interactions between species and DPF rel shows that irradiance-elicited plasticity was of similar extent in the two species (Fig. 5).

In the case of T. melinonii, the responses of leaf traits to DPFrel were as expected a large increase of LMA and its components density and thickness, an increase of area-based photosynthetic capacity and leaf N content (Fig. 5). Mass-based photosynthetic capacity did not change with DPFrel confirming that modulation of LMA and nitrogen content per area ( Na) were the main drivers of the effects of DPFrel on photosynthesis. As a consequence, PNUE increased slightly with increases in DPFrel (Tab. II, Fig. 5). Only few interactions were detected between ASD and DPFrel : thickness, Cm and Chl/N displayed such interactions that resulted in only minor changes in the relationship with DPFrel at a given ASD (Tab. II, Fig. 6).

In the case of Dicorynia guianensis, DPFrel affected all leaf traits except area based photosynthetic capacity and photosynthetic N use efficiency. LMA and to a lesser extent, leaf thickness, density and total leaf N content significantly increased along the DPFrel gradient (Fig. 5). Leaf C content (Cm) was unaffected by DPFrel and by ASD (Tab. III). The relationship between Vcmax a and DPFrel differed according to ASD: Vcmax a ranged from 21.5 to 31.5 μ mol m−2 s−1 in ASD 2 whereas no plasticity was detected in ASD 1 (Fig. 7). Mass-based photosynthetic capacities ( Vcmax m, Jmax m) increased with DPFrel while Rd m decreased. Photosynthetic N use efficiency (Vcmax /N and Jmax /N) remained unaffected by DPFrel (Tab. III).

4. DISCUSSION

Here we tested whether increasing structural complexity, as described by successive stages of architectural development (ASD) has an impact on individual leaf traits. We chose to compare three successive ASDs of young trees to evidence ontogenetic shifts in leaf traits. Nevertheless, we are aware that dimension, size and structural complexity are inter-related properties in trees, and that it is therefore uneasy to separate the effects of tree dimensions from those of structural complexity (i.e., ASDs). This is illustrated by the fact that height and basal diameter of the sampled individuals were always much larger in ASD 3 than in ASD 1. The time-course of successive ASDs is modulated by environment, and faster transitions from an ASD to the next one occur under favorable growth conditions (Barthelemy and Caraglio, 2007). Indeed, branching and crown development occur later and at a larger height in the understorey than in open areas. Consequently, a given ASD can be reached at different total heights and at different chronological ages depending on environmental conditions (Nicolini et al., 2000; Roggy et al., 2005). In our case, it was unfortunately impossible to retrospectively estimate the real age of the sampled young trees and we are unable to disentangle the effects of age, dimension and ASDs.

thumbnail Figure 6

Relationships between leaf thickness (μ m), leaf C content (Cm, mg g−1), ratio of chlorophyll content on N content (Chl/N, μ mol mmol−1) and relative daily photon flux (logDPFrel, %) for the three stages of architectural development (ASD 1, white dots, ASD 2 grey diamonds and ASD 3 black dots) of Tachigali melinonii. Regression equations, level of significance (P) and coefficient of determination (r2) are indicated for statistically significant linear regressions. Significant levels: ns, P > 0.05; * P < 0.05; ** P < 0.01 and *** P < 0.001.

thumbnail Figure 7

Relationships between Vcmax a (μ mol m− 2 s−1) and relative daily photon flux (logDPFrel, %) for the three architectural stage of development (ASD 1, white dots and ASD 2 grey diamonds) of Dicorynia guianensis. Regression equations, levels of significance (P) and coefficient of determination (r2) are indicated for statistically significance linear regressions. Significance levels: ns, P > 0.05; * P < 0.05; ** P < 0.01 and *** P < 0.001.

4.1. Differences in leaf traits among ASDs

Structural and physiological leaf traits were affected by ASDs (i.e., ontogeny) in saplings of Tachigali melinonii and Dicorynia guianensis. We concentrated on juvenile stages of development because they are of central importance for the regeneration of forest stands. At this stage, selection of individuals and species is the most severe during the life-cycle of a forest ecosystem. Moreover, ontogeny and responses to irradiance are major factors that influence the ability of a species to outcompete neighbouring species and to survive in older stages of the forest (Valladares and Niinemets, 2008).

Leaf mass-to-area ratio (LMA) is known to display large ontogenetic shifts as well as irradiance-elicited plasticity: it increases with age and with relative irradiance (e.g. Lusk and Warton, 2007; Niinemets, 1999a, 2006; Thomas and Winner, 2002). In the present study we observed a significant increase of LMA with the sequence of ASDs in juvenile T. melinonii and to a lesser extent D. guianensis. Such an increase can be attributed to different factors:

(i) Increased thickness, in particular of the palisade parenchyma as observed by Kenzo et al. (2006) in some dipterocarp species. Thiswas probably the case in T. melinonii (and to a lesser extent in D. guianensis) where leaf thickness increased with ASDs;

(ii) Increased leaf density due to enhanced investment into structural support required to sustain larger lamina (Lusk and Warton, 2007 and references therein); increased density was indeed recorded in T. melinonii but again not in D. guianensis;

(iii) Increased structural complexity due to the occurrence of branching may lead to restrictions of water supply to leaves due to hydraulic constrictions at branch bases; they may constrain leaf growth and result in smaller leaves (Niinemets,2002; Thomas and Winner, 2002). To our knowledge, the impact of successive ASDs on hydraulic properties has never been explicitly addressed.

In general, the ontogenetic increase in LMA leads in parallel to an increase in area-based and a stability of mass-based photosynthesis (Mediavilla and Escudero, 2003; Niinemets, 2002; Roggy et al., 2005; Thomas and Winner, 2002). Photosynthetic capacity increased in successive ASDs: area-based apparent maximal carboxylation rate (Vcmax a) increased from 25 to 41 and from 20 to 26 μ mol m−2 s − 1 for T. melinonii and D. guianensis, respectively. The effect on mass-based photosynthetic capacity was much smaller but still significant; this confirms that changes in area-based photosynthesis were mainly driven by variations in LMA. Parallel to photosynthetic capacity, total leaf N content increased with ontogenetic stages for T. melinonii whereas it remained constant for D. guianensis. As a consequence, photosynthetic N use-efficiency, measured here as the ratio of Vcmax /N and Jmax /N, was stable among ASDs. A slight decline of the ratio Chl/N was mainly due to increasing N content, whereas Chl m did not differ among ASDs (for T. melinonii).

Vcmax and Jmax increased in parallel, leading to a stability of the ratio among the three stages for T. melinonii. For D. guianensis, the ratio Jmax / Vcmax slightly increased from ASD 1 to ASD 2. Such a variability may be related to the method used to estimate the two parameters. We assumed an infinite value of internal leaf conductance (gi) whereas recent studies have shown that gi is not negligible and may result in a significant underestimation of Vcmax (Ethier and Livingston, 2004).

These results indicate that ontogenetic events like the development of branching do not necessarily induce similar effects on leaf traits in different species. Nevertheless, in the case of D. guianensis larger changes in leaf structure or physiology may appear at later ontogenetic stages, with the development of a complex crown for instance.

4.2. Difference in irradiance-induced plasticity of leaf traits among ASDs

As expected, irradiance affected a large majority of structural and physiological leaf traits, on-line with a number of published results (Givnish, 1988; Rozendaal et al., 2006; Valladares et al., 2000). However, irradiance-elicited plasticity could differ between ASDs (Wright and McConnaughay, 2002). LMA, leaf thickness and density increased as expected with increasing relative irradiance (Popma et al., 1992; Rijkers et al., 2000; Rozendaal et al., 2006; Valladares et al., 2000). Total leaf N decreased while a relative stability was expected from the literature (Ellsworth and Reich, 1993; Niinemets, 1997b; Niinemets and Kull, 1998; Rijkers et al., 2000). However this result fits with observations on co-occurring species of French Guiana grown in an experimental shade-house (Coste, 2008). Consequently, the increase of area-based nitrogen was mainly due to the large plasticity of LMA. Shade leaves generally display also higher Chlm (e.g. Poorter et al., 1995; Rijkers et al., 2000) which was confirmed here. Plasticity of leaf C content, Cm, differed between species. Thus Cm tended to increase for T. melinonii while it remained constant with irradiance (as with ontogeny) for D. guianensis. This result is consistent with the significant interspecific-diversity of irradiance-elicited plasticity of leaf specific construction costs (strongly correlated to Cm) put in evidence in co-occurring species (Coste, 2008).

It is generally accepted that area-based photosynthetic capacity increases with irradiance (Givnish, 1988; Rijkers et al., 2000), mainly as a consequence of structural adjustment (Niinemets,1999a; Roggy et al., 2005) which is confirmed by our results. A similar plasticity, was observed in the two species for a majority of leaf traits. However, mass-based photosynthetic capacity (Vcmax m, Jmax m and Rd m) was not affected by irradiance in T. melinonii whereas it significantly decreased with increasing irradiance for D. guianensis. Conversely,PNUE (i.e., both Vcmax/N and Jmax/N) increased with irradiance for T. melinonii whereas no effect was observed for D. guianensis.

Globally, irradiance-elicited plasticity remained constant among ASDs in the two species, i.e., only small interaction effects were detected between relative irradiance and ontogeny. This may not hold true when scaling up to branch and to crown level; it was indeed shown that ontogeny has an important effect on the plastic responses to light for traits at higher hierarchical scale (Valladares and Niinemets, 2008).

In T. melinonii leaf thickness and carbon content were independent of irradiance in ASD 1 while density increased; this contradicts partly the hypothesis that high leaf density may be partly associated with the presence of tissues with high lignin contents (Niinemets, 1999b; Niinemets and Kull, 1998).

We found only little interaction between ASDs and irradiance, and inter-specific differences were maintained along ontogenetic stages and irradiance. The two species are known to be rather shade-intolerant and the observed differences cannot be attributed to different ecological requirements. Moreover, to characterize shade tolerance as a trait, there is a need to scale up from leaf to at least branch level (if not canopy level) and to take into account potential effects of leaf life span (Lusk and Warton, 2007; Lusk et al., 2008; Niinemets, 2006; Valladares and Niinemets, 2008). Taking into account such integrative traits would be a target for future developments in the assessment of potential interactions between ontogeny and irradiance on tree physiology and carbon balance.

5. CONCLUSION

The present study, conducted in tropical rain forests, showed for two species that three successive ontogenetic stages of development induced visible changes in leaf traits, close to those described in other species (increased leaf-mass-to area ratio, increased area-based photosynthesis). The expected irradiance-elicited plasticity in these traits was also observed, with basically large increases in leaf-mass-to area ratio, and in area-based photosynthetic capacity. There was almost no interaction between the plastic response and the ontogenetic stage: we observed a similar plasticity at all stages and in the two species with the exception of a few minor traits. Such results also underline that leaf life-span and more integrative traits, like branch properties, should equally be considered when addressing the potential interactions between irradiance-elicited plasticity and ontogenic development.

Acknowledgments

S. Coste was supported by a PhD grant from the French ministry for higher education and research, and by a grant for technical cooperation in overseas regions of France. The authors are indebted to Jacqueline Marchand for carbon and nitrogen analyses and to Pascal Imbert and Julien Bachelier for their help throughout the experiment Particular thanks are due to Meriem Fournier, head of the UMR Ecofog at the time of this research, for continuous support and encouragement.


2

Available online only at www.afs-journal.org

References

  1. Baraloto C., Morneau F., Bonal D., Blanc L., and Ferry B., 2007. Seasonal water stress tolerance and habitat associations within four neotropical tree genera. Ecology 88: 478–489 [PubMed] [CrossRef].
  2. Barthelemy D. and Caraglio Y., 2007. Plant architecture: A dynamic, multilevel and comprehensive approach to plant form, structure and ontogeny. Ann. Bot. 99: 375–407 [PubMed] [CrossRef].
  3. Bonal D., Sabatier D., Montpied P., Tremeaux D., and Guehl J.M., 2000. Interspecific variability of $\delta^{13}$C among trees in rainforests of French Guiana: functional groups and canopy integration. Oecologia 124: 454–468 [CrossRef].
  4. Bradshaw A.D., 1965. Evolutionary significance of phenotypic plasticity in plants. Adv. Genet. 13: 115–155.
  5. Cavender-Bares J. and Bazzaz F.A., 2000. Changes in drought response strategies with ontogeny in Quercus rubra: implications for scaling from seedlings to mature trees. Oecologia 124: 8–18 [CrossRef].
  6. Coleman J.C., McConnaughay K.D.M., and Ackerly D.D., 1994. Interpreting phenotypic variation in plants. Trends Ecol. Evol. 9: 187–191 [CrossRef].
  7. Cornelissen J.H.C., Cerabolini B., Castro-Díez P., Villar-Salvador P., Montserrat-Martí G., Puyravaud J.P., Maestro M., Werger M.J.A., and Aerts R., 2003. Functional traits of woody plants: correspondence of species rankings between field adults and laboratory-grown seedlings? J. Veg. Sci. 14: 311–322 [CrossRef].
  8. Coste S., 2008. Diversité et plasticité des traits foliaires en forêt tropicale humide. Une analyse coût/bénéfice de l'acquisition de carbone, Ph.D. thesis, Nancy, AgroParisTech-Engref.
  9. Dreyer E., Le Roux X., Montpied P., Daudet F.A., and Masson F., 2001. Temperature response of leaf photosynthetic capacity in seedlings from seven temperate tree species. Tree Physiol. 21: 223–232 [PubMed].
  10. Ellsworth D.S. and Reich P.B., 1993. Canopy structure and vertical patterns of photosynthesis and related leaf traits in a deciduous forest. Oecologia 96: 169–178 [CrossRef].
  11. Ethier G.J. and Livingston N., 2004. On the need to incorporate sensitivity to CO2 transfer conductance into the Farquhar-von Caemmerer-Berry leaf photosynthesis model. Plant Cell Environ. 27: 137–153 [CrossRef].
  12. Farquhar G.D., von Caemmerer S., and Berry J.A., 1980. A biochemical model of photosynthetic CO2 assimilation in leaves of C3 species. Planta 149: 78–90 [CrossRef].
  13. Favrichon V., 1994. Classification des espèces arborées en groupes fonctionnels en vue de la réalisation d'un modèle de dynamique de peuplement en forêt guyanaise. Rev. Ecol. (terre et vie) 49: 379–403.
  14. Givnish T.J., 1988. Adaptation to sun and shade: A whole plant perspective. Aust. J. Plant. Physiol. 15: 63–92.
  15. Harley P.C., Thomas R.B., Reynolds J.F., and Strain B.R., 1992. Modelling photosynthesis of cotton grown in elevated CO2. Plant Cell Environ. 15: 271–282 [CrossRef].
  16. Hoflacher H. and Bauer H., 1982. Light acclimation in leaves of the juvenile and adult life phases of ivy (Hedera helix). Physiol. Plant. 56: 177–182 [CrossRef].
  17. Holscher D., 2004. Leaf traits and photosynthetic parameters of saplings and adult trees of co-existing species in a temperate broad-leaved forest. Bas. App. Ecol. 5: 163–172 [CrossRef].
  18. Kenzo T., Ichie T., Watanabe Y., Yoneda R., Ninomiya I., and Koike T., 2006. Changes in photosynthesis and leaf characteristics with tree height in five dipterocarp species in a tropical rain forest. Tree Physiol. 26: 865–873 [PubMed].
  19. Le Roux X., Grand S., Dreyer E., and Daudet F.A., 1999. Parameterization and testing of a biochemically based photosynthesis model for walnut trees (Juglans regia) and seedlings. Tree Physiol. 19: 481–492 [PubMed].
  20. Lusk C.H. and Warton D.I., 2007. Global meta-analysis shows that relationships of leaf mass per area with species shade tolerance depend on leaf habit and ontogeny. New Phytol. 176: 764–774 [PubMed] [CrossRef].
  21. Lusk C.H., Falster D.S., Jara-Vergara C.K., Jimenez-Castillo M., and Saldana-Mendoza A., 2008. Ontogenetic variation in light requirements of juvenile rainforest evergreens. Funct. Ecol. 22: 454–459 [CrossRef].
  22. McConnaughay K.D.M. and Coleman J.S., 1999. Biomass allocation in plants: ontogeny or optimality? A test along three resource gradients. Ecology 80: 2581–2593 [CrossRef].
  23. Mediavilla S. and Escudero A., 2003. Mature trees versus seedlings: Differences in leaf traits and gas exchange patterns in three co-occurring Mediterranean oaks. Ann. For. Sci. 60: 455–460 [CrossRef] [EDP Sciences].
  24. Molino J.-F. and Sabatier D., 2001. Tree diversity in tropical rain forests: A validation of the intermediate disturbance hypothesis. Science 294: 1702–1704 [PubMed] [CrossRef].
  25. Nicolini E., 1998. Architecture et gradients morphogénétiques chez de jeunes hêtres (Fagus sylvatica L. Fagaceae) en milieu forestier. Can. J. Bot. 76: 1232–1244 [CrossRef].
  26. Nicolini E., Barthélémy D., and Heuret P., 2000. Influence de la densité du couvert forestier sur le développement architectural de jeunes chênes sessiles, Quercus petraea (Matt.) Liebl. (Fagaceae), en régénération forestière. Can. J. Bot. 78: 1531–1344 [CrossRef].
  27. Niinemets Ü., 1997a. Distribution patterns of foliar carbon and nitrogen as affected by tree dimensions and relative light conditions in the canopy of Picea abies. Trees 11: 144–154 [CrossRef].
  28. Niinemets Ü., 1997b. Role of foliar nitrogen in light harvesting and shade tolerance of four temperate deciduous woody species. Funct. Ecol. 11: 518–531 [CrossRef].
  29. Niinemets Ü., 1999a. Components of leaf dry mass per area - thickness and density - alter leaf photosynthetic capacity in reverse directions in woody plants. New Phytol. 144: 35–47 [CrossRef].
  30. Niinemets Ü., 1999b. Energy requirement for foliage formation is not constant along canopy light gradients in temperate deciduous trees. New Phytol. 141: 459–470 [CrossRef].
  31. Niinemets Ü., 2002. Stomatal conductance alone does not explain the decline in foliar photosynthetic rates with increasing tree age and size in Picea abies and Pinus sylvestris. Tree Physiol. 22: 515–535.
  32. Niinemets Ü., 2006. The controversy over traits conferring shade-tolerance in trees: ontogenetic changes revisited. J. Ecol. 94: 464–470 [CrossRef].
  33. Niinemets Ü. and Kull O., 1998. Stoichiometry of foliar carbon constituents varies along light gradients in temperate woody canopies: implications for foliage morphological plasticity. Tree Physiol. 18: 467–479 [PubMed].
  34. Poorter L., 1999. Growth responses of 15 rain-forest tree species to a light gradient: the relative importance of morphological and physiological traits. Funct. Ecol. 13: 396–410 [CrossRef].
  35. Poorter L., Oberbauer S.F., and Clark D.B., 1995. Leaf optical properties along a vertical gradient in a tropical rain forest canopy in Costa Rica. Am. J. Bot. 82: 1257–1263 [CrossRef].
  36. Poorter L., Bongers L., and Bongers F., 2006. Architecture of 54 moist-forest tree species: Traits, trade-offs, and functional groups. Ecology 87: 1289–1301 [PubMed] [CrossRef].
  37. Popma J., Bongers F., and Werger M.J.A., 1992. Gap-dependence and leaf characteristics of trees in a tropical lowland rain forest in Mexico. Oikos 63: 207–214 [CrossRef].
  38. Rijkers T., Pons T.L., and Bongers F., 2000. The effect of tree height and light availability on photosynthetic leaf traits of four neotropical species differing in shade tolerance. Funct. Ecol. 14: 77–86 [CrossRef].
  39. Roggy J.-C., Nicolini E., Imbert P., Caraglio Y., Bosc A., and Heuret P., 2005. Links between tree structure and functional leaf traits in the tropical forest tree Dicorynia guianensis Amshoff (Caesalpiniaceae). Ann. For. Sci. 62: 553–564 [CrossRef] [EDP Sciences].
  40. Rozendaal D.M.A., Hurtado V.H., and Poorter L., 2006. Plasticity in leaf traits of 38 tropical tree species in response to light; relationships with light demand and adult stature. Funct. Ecol. 20: 207–216.
  41. Sterck F.J. and Bongers F., 1998. Ontogenetic changes in size, allometry, and mechanical design of tropical rain forest trees. Am. J. Bot. 85: 266–272 [CrossRef].
  42. Thomas S.C. and Winner W.E., 2002. Photosynthetic differences between saplings and adult trees: an integration of field results by meta-analysis. Tree Physiol. 22: 117–127 [PubMed].
  43. Valladares F. and Niinemets U., 2008. Shade tolerance, a key plant feature of complex nature and consequences. Annu. Rev. Ecol. Evol. Syst. 39: 237–57 [CrossRef].
  44. Valladares F., Wright S.J., Lasso E., Kitajima K., and Pearcy R.W., 2000. Plastic phenotypic response to light of 16 congeneric shrubs from a Panamanian rainforest. Ecology 81: 1925–1936.
  45. Veneklaas E.J. and Poorter L., 1988. Growth and carbon partitioning of tropical tree seedlings in contrasting light environments. In: H. Lambers, H. Poorter and M. Van Vuuren (Eds.), Inherent variation in plant growth: physiological mechanisms and ecological consequences, Backhuys publishers, Leiden, The Netherlands, pp. 337–355.
  46. Warren C.R., 2006. Estimating the internal conductance to CO2 movement. Funct. Plant Biol. 33: 431–442 [CrossRef].
  47. Wright S.D. and McConnaughay K.D.M., 2002. Interpreting phenotypic plasticity: the importance of ontogeny. Plant Species Biol. 17: 119–131 [CrossRef].

Online material

Download PDF file.

Appendix

thumbnail Appendix 1

Relationships between chlorophyll content (Chla, μ molChl m− 2) and SPAD values (SPAD-502, Minolta, Osaka, Japan). Chlorophyll content was measured on 2.3 cm2 leaf discs after extracting in DMSO (n = 35).

Appendix 2

Mean values of leaf traits in the different stages of architectural development for Tachigali melinonii and Dicorynia guianensis. Values of adjusted means (± standard error) at a common DPFrel of 26% are given for the different ASDs.

Appendix 3

Result of a one-way Ancova with stage of architectural development (ASD = 2) as main factor and relative irradiance (DPFrel, %) as covariables applied to Tachigali melinonii (Harms) Barneby. F -value, level of significance (P) and coefficient of determination (r2) of the model are given. DPFrel, LMA and Jmax a were log transformed prior to analysis.Values of adjusted means (± standard error) were also given (covariate mean = 21%) for the three stages of Architectural development (ASD). Cm and thickness were analyzed with a separate slope model. Vcmax m, Jmax m and Rd m were not linearly correlated with DPFrel and thus they were analyzed with a one way ANOVA. Significant levels: ns, P > 0.05; * P < 0.05; ** P < 0.01 and *** P < 0.001.

All Tables

Table I

Primary parameters of the photosynthesis model used in this work. For details, see Dreyer et al. (2001).

Table II

Results of a one-way Ancova with stage of architectural development (ASD = 3) as main factor and relative irradiance (DPFrel, %) as covariant for Tachigali melinonii. F -values, levels of significance (P) and coefficients of determination (r2) of the model are given. DPFrel, H, D, LMA, Na and Jmax a were log-transformed prior to analysis. Cm and thickness were analyzed with a separate slope model. Vcmax m, Jmax m and Rd m were not linearly correlated with DPFreland thus they were analyzed with a one way ANOVA. Significance levels: ns, P > 0.05;*, P < 0.05; **, P < 0.01 and ***, P < 0.001. The corresponding adjusted means at a 26% DPFrel are displayed in Figure and numerical values are provided in the on-line only Appendix 2.

Table III

Results of a one-way Ancova with stage of architectural development (ASD 1 and ASD 2) as main factor and relative DPF (DPFrel, %) as covariant applied to Dicorynia guianensis. F -value, level of significance (P) and coefficient of determination (r2) of the model are given. DPFrel, H, D and Jmax a were log-transformed prior to analysis. Vcmax a was analyzed with a separate slope model. D, Cm,Jmax a, Jmax/ Vcmax, Vcmax/ N and Jmax/ N were not linearly correlated with DPFreland thus they were analyzed with a one way ANOVA. Significance levels: ns, P > 0.05; *, P < 0.05; **, P < 0.01 and ***, P < 0.001. The corresponding adjusted means at 26% DPFrelare displayed in Figure and numerical values are provided in the on-line only Appendix 2.

Appendix 2

Mean values of leaf traits in the different stages of architectural development for Tachigali melinonii and Dicorynia guianensis. Values of adjusted means (± standard error) at a common DPFrel of 26% are given for the different ASDs.

Appendix 3

Result of a one-way Ancova with stage of architectural development (ASD = 2) as main factor and relative irradiance (DPFrel, %) as covariables applied to Tachigali melinonii (Harms) Barneby. F -value, level of significance (P) and coefficient of determination (r2) of the model are given. DPFrel, LMA and Jmax a were log transformed prior to analysis.Values of adjusted means (± standard error) were also given (covariate mean = 21%) for the three stages of Architectural development (ASD). Cm and thickness were analyzed with a separate slope model. Vcmax m, Jmax m and Rd m were not linearly correlated with DPFrel and thus they were analyzed with a one way ANOVA. Significant levels: ns, P > 0.05; * P < 0.05; ** P < 0.01 and *** P < 0.001.

All Figures

thumbnail Figure 1

Stages of architectural development (ASD) of Dicorynia guianensis (DG) and Tachigali melinonii (TM): ASD1: saplings with un-branched main stem; ASD2: saplings with sparsely branched main stem (order 2 axes remain un-branched); ASD3: immature trees (diamonds: sequential axes; crosses: reiterated axes). For a complete description of the different ASDs, see Material and methods.

In the text
thumbnail Figure 2

Effects of successive stages of architectural development (ASD 1 to 3) on height (H, cm) for T. melinonii (white) and D. guianensis (black). Values are least square means (± standard error) adjusted to a common DPFrel of 26.1 %. Differences were tested with the multiple comparisons post-hoc test of Tukey (P < 0.05) and significant differences are indicated by different letters.

In the text
thumbnail Figure 3

Relationships between height (H, cm) and relative daily photon flux (logDPFrel) for the three stages of architectural development (ASD 1, white dots, ASD 2 grey diamonds and ASD 3 black dots) and for two species (T. melinonii and D. guianensis). Regression equations, level of significance (P) and coefficient of determination (r2) are indicated for statistically significant linear regressions. Significance levels: ns, P > 0.05;*, P < 0.05; **, P < 0.01 and ***, P < 0.001.

In the text
thumbnail Figure 4

Effects of stages of architectural development (ASD 1 to 3) on leaf traits for T. melinonii (white) and D. guianensis (black). Values are least square means (± Standard error) adjusted at a common relative daily photon flux (DPFrel=26.1 %). Significant differences were tested with the multiple comparisons post-hoc test of Tukey (P < 0.05) and indicated by different letters (a, b, c for T. melinonii and a, bc for D. guianensis).

In the text
thumbnail Figure 5

Relationships between leaf traits and relative daily photon flux (logDPFrel) for T. melinonii (white dots) and D. guianensis (black dots). Simple linear regressions were represented by solid lines for T. melinonii and by dotted lines for D. guianensis. Regression equation and r2 are indicated for statistically significant linear regressions (P < 0.05). Significance levels: ns, P > 0.05; *, P 60.05; **, P 60.01; *** and P 60.001.

In the text
thumbnail Figure 6

Relationships between leaf thickness (μ m), leaf C content (Cm, mg g−1), ratio of chlorophyll content on N content (Chl/N, μ mol mmol−1) and relative daily photon flux (logDPFrel, %) for the three stages of architectural development (ASD 1, white dots, ASD 2 grey diamonds and ASD 3 black dots) of Tachigali melinonii. Regression equations, level of significance (P) and coefficient of determination (r2) are indicated for statistically significant linear regressions. Significant levels: ns, P > 0.05; * P < 0.05; ** P < 0.01 and *** P < 0.001.

In the text
thumbnail Figure 7

Relationships between Vcmax a (μ mol m− 2 s−1) and relative daily photon flux (logDPFrel, %) for the three architectural stage of development (ASD 1, white dots and ASD 2 grey diamonds) of Dicorynia guianensis. Regression equations, levels of significance (P) and coefficient of determination (r2) are indicated for statistically significance linear regressions. Significance levels: ns, P > 0.05; * P < 0.05; ** P < 0.01 and *** P < 0.001.

In the text
thumbnail Appendix 1

Relationships between chlorophyll content (Chla, μ molChl m− 2) and SPAD values (SPAD-502, Minolta, Osaka, Japan). Chlorophyll content was measured on 2.3 cm2 leaf discs after extracting in DMSO (n = 35).

In the text