Stomatal CO2 responsiveness and photosynthetic capacity of tropical woody species in relation to taxonomy and functional traits.

Stomatal CO2 responsiveness and photosynthetic capacity vary greatly among plant species, but the factors controlling these physiological leaf traits are often poorly understood. To explore if these traits are linked to taxonomic group identity and/or to other plant functional traits, we investigated the short-term stomatal CO2 responses and the maximum rates of photosynthetic carboxylation (V cmax) and electron transport (J max) in an evolutionary broad range of tropical woody plant species. The study included 21 species representing four major seed plant taxa: gymnosperms, monocots, rosids and asterids. We found that stomatal closure responses to increased CO2 were stronger in angiosperms than in gymnosperms, and in monocots compared to dicots. Stomatal CO2 responsiveness was not significantly related to any of the other functional traits investigated, while a parameter describing the relationship between photosynthesis and stomatal conductance in combined leaf gas exchange models (g 1) was related to leaf area-specific plant hydraulic conductance. For photosynthesis, we found that the interspecific variation in V cmax and J max was related to within leaf nitrogen (N) allocation rather than to area-based total leaf N content. Within-leaf N allocation and water use were strongly co-ordinated (r 2 = 0.67), such that species with high fractional N investments into compounds maximizing photosynthetic capacity also had high stomatal conductance. We conclude that while stomatal CO2 responsiveness of tropical woody species seems poorly related to other plant functional traits, photosynthetic capacity is linked to fractional within-leaf N allocation rather than total leaf N content and is closely co-ordinated with leaf water use.

Introduction

Anthropogenic fossil fuel burning and land use change have increased the atmospheric pan class="Chemical">carbon dioxide concentration [pan class="Chemical">CO2] by over 40% (Ciais et al. 2013), with today’s concentration of over 400 μmol mol−1 being the highest in approximately 40 millions years (Frank et al. 2015). This large and rapid ongoing increase in [CO2] has profound impacts on land plants, which typically respond to altered [CO2] by increasing leaf photosynthesis (A) and decreasing stomatal conductance (g s; Ainsworth and Rogers 2007). However, the direct stomatal responses to a short-term increase in [CO2] vary greatly among species and experiments, from no change to 75% reduction in g s at doubled compared to ambient [CO2] (Morison 1998). Understanding this variation is important, since the short-term stomatal CO2 response appears to be an important determinant for the long-term effect of growth in elevated [CO2] on g s under field conditions. In so called free-air CO2 enrichment (FACE) experiments with trees, the interspecific variation in long-term (years) effects of growth under elevated [CO2] on g s was significantly and positively related to the variation in short-term stomatal CO2 responsiveness, measured as the effect of short (hours) interruptions in CO2 enrichment on sap flow (Cech et al. 2003; Keel et al. 2007; Tor-ngern et al. 2015) or leaf g s (Maier et al. 2008; Domec et al. 2009; Onandia et al. 2011) in some of these experiments (Fig. 1; r 2 = 0.81; P = 0.014). Tree species that exhibit a pronounced direct stomatal closure response to a short-term increase in [CO2] are thus also likely to develop a long-term g s reduction in the field. Therefore, a better understanding of the factors that control the large natural variation in short-term stomatal CO2 responses would provide indication of which plant species and groups that are likely to experience large increases in water-use efficiency in future higher [CO2] scenarios.

Fig. 1

Relationship between the short-term response of g s to elevated [CO2] [measured during short (hours) interruptions in CO2 enrichment] and the long-term effect of growth under elevated [CO2] on g s in temperate forest free-air CO2 enrichment (FACE) experiments. Regression statistics are shown in the figure. Based on data from Cech et al. (2003), Keel et al. (2007), Maier et al. (2008), Domec et al. (2009), Onandia et al. (2011) and Tor-ngern et al. (2015)

Relationship between the short-term response of g s to elevated [pan class="Chemical">CO2] [measured during short (hours) interruptions in CO2 enrichment] and the long-term effect of growth under elevated [CO2] on g s in temperate forest free-air CO2 enrichment (FACE) experiments. Regression statistics are shown in the figure. Based on data from Cech et al. (2003), Keel et al. (2007), Maier et al. (2008), Domec et al. (2009), Onandia et al. (2011) and Tor-ngern et al. (2015)

Both the short-term stomatal pan class="Chemical">CO2 response and the long-term effect of growth under elevated [CO2] on g s have been found to be weaker in gymnosperms compared to angiosperms (Medlyn et al. 2001; Brodribb et al. 2009). Brodribb et al. (2009) even suggested that gymnosperms lack a primary stomatal CO2 response, while others have found that this is not the case (Haworth et al. 2013). To our knowledge, there has been no study designed to explore possible taxonomic patterns in the large interspecific variation in stomatal CO2 responsiveness among different angiosperm taxa. Furthermore, possible functional determinants of the interspecific variation in stomatal CO2 responsiveness have also been poorly explored. If stomatal CO2 responses are linked to physiological, chemical or structural plant traits for which more knowledge is available, this could facilitate future trait-based modelling of plant water-use under rising [CO2] (Van Bodegom et al. 2012). For example, it is possible that (but unknown if) stomatal CO2 responsiveness is linked to stomatal density since, in general, plants with high stomatal density have more efficient stomatal regulation (Franks et al. 2009). It is also plausible that stomatal CO2 responsiveness is linked to hydraulic traits since plants with high hydraulic conductance often exhibit tight stomatal control over transpiration to avoid severe cavitation (Bond and Kavanagh 1999). Since stomatal regulation is an energy consuming process, CO2 sensitivity may also be linked to indicators of leaf metabolic activity, such as photosynthetic capacity.

Most land-surface and ecosystem models apply combined stomatal–photosynthesis models (Ball et al. 1987; Leuning 1995; Medlyn et al. 2011), in which plant pan class="Chemical">carbon and pan class="Chemical">water vapor fluxes are linked as they both pass through plant stomata. The empirical slope parameter of the combined stomatal–photosynthesis model (g 1), which is inversely proportional to leaf water use efficiency (at constant environmental conditions; Medlyn et al. 2011), has been determined for a large number of plant functional types and in different biomes (Lin et al. 2015). In angiosperm trees, it was found that g 1 was positively related to wood density, reflecting a higher cost of wood construction per unit water transported in species with higher wood density (Lin et al. 2015). However, the relationship with wood density was quite weak (r 2 = 0.21; Lin et al. 2015) and when gymnosperms were included it was not significant at all. There is thus a need to further explore the relationship of g 1 with wood density and other plant traits, such as hydraulic conductance. Hydraulic conductance and g s are often co-ordinated in trees (Brodribb and Jordan 2008) and since high hydraulic conductance essentially reduces the carbon cost of water use (Prentice et al. 2014) it is plausible that (but unexplored if) g 1 is positively related to hydraulic conductance. With respect to taxonomic patterns, there are typically distinct differences in water-use strategy between angiosperm and gymnosperm trees, with angiosperms typically exhibiting higher g s, hydraulic conductance and g 1 at the expense of considerably smaller hydraulic safety margin compared to gymnosperms (Choat et al. 2012; Lin et al. 2015).

There is also considerable uncertainty with respect to the variation in the maximum rates of photosynthetic carboxylation (V cmax) and electron transport (J max) among plant species in general (Ali et al. 2015) and tropical tree species in particular (Dusenge et al. 2015). Current vegetation models typically base their values of V cmax and J max on area-based total leaf N content (e.g., Zaehle et al. 2005; Rogers 2014). However, a recent global meta-analysis found that V cmax and J max were more closely related to within-leaf n class="Chemical">nitrogen (N) allocation than to total leaf N content (Ali et al. 2015), as also observed in two studies with tropical tree species (Coste et al. 2005; Dusenge et al. 2015). n>n class="Chemical">Part of the reason for why interspecific variation in photosynthetic capacity is often poorly related to total leaf nutrient content may be linked to differences in within-leaf N allocation among species with different successional strategies. It is well established that early-successional species invest more leaf N into N-rich molecules involved in photosynthesis and respiration compared to late-successional tree species, which typically have larger fractional N investments in rather N-poor structural compounds (Raaimakers et al. 1995; Valladares and Niinemets 2008). Dusenge et al. (2015) further suggested that there is a trade-off in how the photosynthetic leaf N is invested, with pioneer species making higher relative investments into compounds maximizing photosynthetic carboxylation and bioenergetics (i.e., V cmax and J max) and climax species investing more into compounds involved in light-harvesting. However, that study included a rather small number of species and this hypothesis, therefore, needs to be further evaluated by field data.

Tropical forests have a fundamental role in the global pan class="Chemical">carbon and pan class="Chemical">water cycles and in controlling the rate of climate change (Lewis 2006). They account for about a third of the global terrestrial primary production (Beer et al. 2010) and their evapotranspiration regulates regional temperature and precipitation and runoff patterns (Bonan 2008). Given all the challenges associated with research in tropical forests (e.g. large stature, biological complexity and logistics), much less is known about stomatal regulation and photosynthesis of trees in these ecosystems compared to temperate forests (Kattge et al. 2009). To improve the poor understanding of the phylogenetic and functional controls of the large interspecific variation in stomatal CO2 responsiveness and photosynthetic capacity among tropical woody species, we examined physiological, chemical and structural plant traits in an evolutionary broad cross-section of tropical woody seed plant species in a common garden experiment in Rwanda. The specific objectives were to answer the following three research questions: (1) Does the stomatal response to a short-term increase in [CO2] vary among major taxonomic groups? (2) What functional characteristics can explain the interspecific variation in stomatal behavior (i.e. in the short-term CO2 response and g 1)? (3) Is the interspecific variation in photosynthetic capacity controlled by differences in within-leaf nutrient allocation rather than by differences in total area-based leaf nutrient content?

Materials and methods

Study sites and plant material

From February 27 until March 25 in 2014, data were collected on mature (having reached reproductive stage) woody plants of 21 species growing at four locations in the surroundings of Butare, Huye district, Rwanda, central Africa: Rwanda Agriculture Board (RAB) Ruhande Arboretum (“Arboretum” hereafter; 2°36′55.2′′S, 29°44′53.8′′E; 1700 m asl; 17 species), Ruhande Fisheries Research Station (2°36′22.1′′S, 29°45′29.7′′E; 1650 m asl; one species), Marist Missionary Sisters of the Society of Mary’s garden (2°36′32.1′′S, 29°44′37.6′′E; 1700 m a.s.l.; two species) and RAB Rubona station (2°38′57.1′′S, 29°45′54.2′′E; 1690 m a.s.l.; one species). The Arboretum was established in 1934 and since then has gathered 227 tree species (50 native to Rwanda) planted, in most part, as replicated (n = 3) monospecific 50 × 50 m plots within its 200 ha plantation area. Out of the 21 species investigated five were native to Africa and 16 were exotic (Table 1).

Table 1

Tree species description and their plot location at the Ruhande Arboretum

Group		Order	Family	Species	Origin	Plot	Mean DBH (cm)	Leaf longevity	Leaf shape	Stand planted (year)
Gymnosperms		Pinales	Araucariaceae	Araucaria angustifolia (Bertol.) Kuntze	Exotic	47	19.2	EG	NL	1938
		Pinales	Cupressacaea	Cupressus lusitanica Mill.	Exotic	144	36.6	EG	NL	1948
		Pinales	Pinaceae	Pinus patula Schlecht. & Cham.	Exotic	48	66.9	EG	NL	1959
		Pinales	Podocarpaceae	Podocarpus latifolius (Thunb.) R.Br. ex Mirb.	Native	43	19.7	EG	NL	1948
		Pinales	Podocarpaceae	Podocarpus falcatus (Thunb.) R.Br. ex Mirb	Native	2	39.2	EG	NL	1955
Angiosperms
Monocots (Commelinids)		Arecales	Arecaceae	Phoenix reclinata Jacq.	Native	a	22.0	EG	BL	~2004
		Poales	Poaceae	Dendrocalamus giganteus Munro	Exotic	293	6.4	EG	BL	2005
		Poales	Poaceae	Bambusa vulgaris Schrad.	Exotic	26	5.0	EG	BL	1934
		Zingiberales	Heliconiaceae	Heliconia rostrata Ruiz & Pavon	Exotic	b	2.2	EG	BL	~2009
		Zingiberales	Musacea	Musa sapientum L.	Exotic	b	5.6	EG	BL	~2009
Dicots
Rosids		Malphigiales	Chrysobalanaceae	Macaranga kilirnandscharica (Pax.)	Native	224	8.0	EG	BL	2008
		Rosales	Rosaceae	Prunus caretta	Exotic	65	15.0	EG	BL	1943
		Myrtales	Myrtaceae	Eucalyptus maculata (Hook.) K.D. Hill & L.A.S. Johnson	Exotic	6	15.7	EG	BL	1934
		Sapindales	Meliaceae	Carapa grandiflora Sprage	Native	11	19.9	EG	BL	1948
		Sapindales	Meliaceae	Cedrela serrata Royle	Exotic	111–112	22.2	EG	BL	1945
		Malvales	Malvaceae	Brachychiton acerifolius (A.Cunn. Ex G.Don) Macarthur & C. Morre	Exotic	75	14.9	DE	BL	1943
Asterids		Lamiales	Bignoniaceae	Jacaranda rnirnosifolia D.Don	Exotic	8	42.4	DE	BL	1934
		Lamiales	Boraginaceae	Cordia alliodora (Ruiz & Pay.) Oken	Exotic	30	9.1	EG	BL	1984
		Lamiales	Oleaceae	Ligustrum lucidum W. T.Aiton	Exotic	67	5.1	EG	BL	1943
		Solanales	Solanaceae	Cyphornandra betaceae Cay.	Exotic	c	2.7	EG	BL	2012
		Asterales	Asteraceae	Tithonia diversifolia (Hemsl.) A.Gray	Exotic	37	1.8	EG	BL	~2008

DBH diameter at breast high, EG evergreen, DE deciduous, NL needle leaf, BL broadleaf)

a Fisheries Research Station Adjacent to the Arboretum

b Marist Missionary Sisters of the Society of Mary’s Garden

c RAB Rubona station

pan class="Chemical">DBH diameter at breast high, EG evergreen, DE deciduous, NL needle leaf, BL broadleaf)

The nearest meteorological station is located ca. 2 km from the Arboretum, showing for the period 2006–2013 mean day and night air temperatures of 20.8 and 17.1 °C, respectively, average humidity of 74%, and mean annual rainfall of 1231 mm. The climate of the region is characterized as tropical humid, with 1.5 °C difference between the warmest and coolest months. Precipitation exhibits some seasonal variation, with the highest and lowest amounts occurring in March–May and June–August, respectively. Additional meteorological information can be found in Nsabimana et al. (2009) and Vårhammar et al. (2015).

The 21 species were selected to represent a broad evolutionary cpan class="Chemical">ross-section of seed plant species, representing four major lineages of seed plants: gymnosperms, monocots [commelinids], rosids and asterids (Table 1). Of the selected species, 16 were silviculturally or agriculturally important exotics, while five were native tropical species. All species exhibited a woody stem and plants investigated were older than two years. In this study, we define wood as tissue produced of secondary xylem, including highly lignified xylem tissue in secondary vascular bundles in monocots. The focus on tropical woody species minimized sources of variation related to growth form and climate of origin, thereby increasing the chance of finding possible taxonomic patterns. For each species, measurements were done on five individual plants, measuring leaf gas exchange on one fully developed leaf of each individual (except for Phoenix reclinata, where two and three leaves were measured in each of two different plants).

Gas exchange measurements

Leaf gas exchange measurements were conducted in situ between 09:00 and 15:00 h using two portable photosynthesis systems (LI-6400 and LI-6400XT with the 2 × 3 cm standard leaf chamber and LI-6400-02B LED light source; Li-Cor Inc., Lincon, NE, USA). For each species, five fully developed, visually healthy, sun-exposed leaves (one per plant) were selected and measured for the response of net photosynthetic rate (A n) to eight different [pan class="Chemical">CO2] [so called A-C i curves; C i = leaf intercellular (pan class="Chemical">CO2)]. In large trees, the measured leaves were attached to lower sun-exposed branches by the edge of the stand. All measurements were conducted at a photosynthetic photon flux density (PPFD) of 1800 µmol m−2 s−1 and a leaf temperature of 25 °C. Following the A-C i curve measurements, direct responses of g s to a short-term change in leaf chamber [CO2] (600 vs. 400 μmol mol−1) were measured on the same leaves, at the same leaf temperature and PPFD. Measurements were taken every minute and the steady-state g s at each [CO2] was recorded after g s values had stabilized for at least 15 min and changed ≤2% within a 5 min period. The order of the [CO2] applied was randomized for the first leaf of a species, and then alternated for the remaining leaves of that species. For each leaf, the leaf-to-air vapour pressure deficit (D) was kept constant (± 0.03 kPa) around a leaf-specific value in the range 1.15–1.80 kPa during the entire measurement.

Shoots of needle-leaf species were also measured using the standard leaf chamber. Needles were removed at the points where the shoot axis passed the chamber gasket. Needles were scanned (Epson V200 Perfection, Epson America Inc., USA) and their projected surface area were estimated using WinRHIZO analysis software (WinRHIZO, Regent Instruments Inc., Canada), and leaf and needle temperatures were calculated using energy balance equations.

Photosynthesis model parameterization

The photosynthesis model by Farquhar et al. (1980), with modifications of photosynthetic temperature dependencies by Bernacchi et al. (2001), was used to parameterize V cmax and J max from A-C i curve data by the least squares method. The rates of V cmax-limited net photosynthesis (A c) and J max-limited net photosynthesis (A j) were calculated using:andwhere K c and K o are Michaelis–Menten constants for n class="Chemical">CO2 and n>n class="Chemical">O2, respectively; Γ* is the [CO2] at which the carboxylation reaction of Rubisco equals the oxygenation reaction; R d is the non-photorespiratory CO2 release in the light; and J is the rate of electron transport. For K c, K o and Γ*, the values at 25 °C as well as the temperature sensitivities were taken from Bernacchi et al. (2001). The internal leaf conductance for CO2 was not estimated and, therefore, ‘apparent’ V cmax and J max values are reported, based on C i rather than on the [CO2] at the chloroplast. The uncertainty of the values of the curvature of the light-response (0.9) and quantum yield of electron transport (0.3 mol electrons mol−1 photons) used when calculating J max from J has only a minor effect on the estimated value of J max (Medlyn et al. 2002). Values of V cmax, J max and R d were determined simultaneously with the only a priori restriction made to the A−C i fitting that data points with C i below 100 µmol mol−1 were forced to be V cmax-limited. The discontinuity of the transition from V cmax- to J max-limitation in the A−C i curve was not smoothed out. J max values were reported only if the A j limited part of the A−C i curve had at least two data points, or from one single data point if A j was at least 10% lower than A c at the C i value of that data point. The latter was the case in nine leaves only and the J max value of these leaves were similar to those of the other leaves. Light-saturated net photosynthesis at a C i of 280 µmol mol−1 (A n280; i.e. assuming a C i that is 70% of an ambient CO2 concentration of 400 µmol mol−1) was calculated based on the fitted photosynthesis model for each leaf.

Leaf-specific values of g 1 were calculated from the steady-state leaf gas exchange measurement at 400 μmol mol−1 pan class="Chemical">CO2 concentration of the air, according the model for optimal stomatal conductance given in Medlyn et al. (2011):where g 0 is the g s at zero A n, 1.6 is the ratio of the diffusivities of water vapour and CO2, and C a is the [CO2] of the air surrounding the leaf (µmol mol−1). Values of g 0 are often small (Medlyn et al. 2011) and were assumed to be zero in the present study. The model parameter g 1, which is the key parameter of the optimal stomatal conductance model and inversely related to water use efficiency, could thus be analytically calculated for each individual leaf.

Additional leaf traits

After the gas exchange measurements, epidermal impressions of around 2 cm2 of the center of the left part of the abaxial leaf surface were made from all measured leaves except pan class="Species">Tithonia diversifolia (due to its hairy leaf surface) to determine stomatal density and length. A thin layer of nail varnish was carefully peeled off after 10 min drying, and was gently placed over a microscope slide and covered and sealed by a cover slip. Six photos of evenly distributed areas of each peel where taken using a microscope (Zeiss, Axio Scope.A1, Germany) equipped with a digital camera (Zeiss, AxioCam MRm, Germany) using 100 × magnification. The photos were treated for higher contrast and definition using ImageJ 1.48v software before measurements of stomatal density and guard cell length. Stomatal density was estimated by calculating the average number of stomata per mm2 from the six photos taken from each peeling, except in Macaranga kilimandcharica where photos could be used from only one of the five peels. Stomatal length was estimated by calculating the average of 30 randomly chosen guard cell pairs from three of the six photos taken (10 from each photo). The estimates of stomatal density and length were used to calculate the maximal g s that was anatomically possible (g smax) according to Franks and Beerling (2009):where d w is the diffusivity of water vapour (0.26 m2 s−1), v is the molar volume of air (24.47 dm3 mol−1), ρ is the number of stomata per m2, a max is the area of a fully open stomatal pore (m2) and l is the depth of stomatal pore (m). Assumptions used in the calculation of stomatal pore length, width and depth in relation to guard cell length were the same as in Franks and Beerling (2009).

Leaf length, width and thickness were recorded. Leaf discs of known area were punched from each leaf, avoiding major veins, to determine leaf mass per unit area (pan class="Chemical">LMA) and leaf N and pan class="Chemical">phosphorus (P) content. Leaves of species with needles (Cupressus lusitanica, Pinus patula, Araucaria angustifolia, Podocarpus latifolius, Podocarpus falcatus) or small leaves (Jacaranda mimosifolia) were scanned (Epson Expression 1600, Epson America Inc., CA, USA) and projected leaf surface area was determined using WinRhizo software (2007, Regents Instruments, Canada). The collected leaf material was dried to constant weight at 70 °C before dry mass was recorded. The leaf discs and needles were milled using a ball mill (Retsch MM301, Haan, Germany) and leaf N concentration was determined using an elemental analyser (Europe EA-GSL, Sercon Ltd., Crewe, UK) coupled to a stable isotope ratio mass spectrometer (20-22, Sercon Ltd. Crewe, UK). Leaf P content was determined for three leaves per species by first treating samples with the Kjeldahl method and later analyzing them with a continuous flow automated analyzer (SEAL Analytical, AutoAnalyzer 3HR, Norderstedt, Germany). Leaf chlorophyll content was determined by extraction of leaf samples in 80% acetone and subsequent filtering and spectrophotometric absorbance measurement, as described by Uddling et al. (2007).

Immediately after the gas exchange measurements, leaves were collected and leaf pan class="Chemical">water potentials (Ψ leaf, MPa) were measured using a Scholander type pressure chamber (SAPS II Model 3115, Soilmoisture Equipment Corp., Santa Barbara, CA, USA). In addition, pre-dawn Ψ leaf measurements were done in the last day of the measurement campaign (except for Phoenix reclinata, Musa sapientum, Heliconia rostrata and Cyphomandra betaceae). The measured values of pre-dawn Ψ leaf were high, and this variable likely did not vary much over the 1 month measurement campaign since precipitation was quite high and temperature similar in both March (109 mm/19.5 °C; measurement month) and February 2014 (89 mm/19.2 °C; preceding month). The evening before the pre-dawn Ψ leaf measurements, one leaf per tree was kept in darkness inside a water vapour saturated Ziploc plastic bag covered in aluminium foil. At pre-dawn on the following day (ca. 05:00 h) the plastic bags were collected and pre-dawn Ψ leaf was measured and used as an estimate of soil water potential (Ψ soil). From data of Ψ soil (i.e., pre-dawn Ψ leaf), leaf transpiration (E), leaf temperature and Ψ leaf, leaf area-specific plant hydraulic conductance (K p) was determined according to:where ṽ is the relative temperature dependent kinematic viscosity of water (set to 1 at 20 °C; dimensionless). Values of ṽ were calculated using leaf-specific values of leaf temperature, which may be a fair simplification since the leaf accounts for about half of the whole-plant hydraulic resistance (Sack and Holbrook 2006). Soil water potential values for the three species lacking pre-dawn data were taken as the mean value across all other species.

A 3 cm piece of wood on a branch close to the measured leaf was cut and its volume (without the bark) was measured by immersion in a volumetric cylinder on an analytical balance (except for pan class="Species">H. rostrata and pan class="Species">M. sapientum which have pseudostem). All pieces of wood were later dried to constant weight at 70 °C and wood mass and density were determined.

The fractional leaf N investments were determined for the following components of the photosynthetic apparatus: Rubisco (NR); bioenergetics, including coupling factors, electron carriers except for photosystems, and Calvin-Benson cycle enzymes except for Rubisco (NB); and light-harvesting complexes and photosystems (NLH). The NR was estimated using the equation and parameters provided by Niinemets and Tenhunen (1997):where V cmax is the maximum rate of carboxylation, 0.16 converts Rubisco to N [g N in Rubisco (g Rubisco)−1] and V cr the specific activity of Rubisco at 25 °C [20.78 μmol n class="Chemical">CO2 (g Rubisco)−1 s−1]. The NB was estimated as:where it is assumed that N in bioenergetics is proportional to J max, that 156 is the ratio of electron transport to cytochrome f content in mol mol−1 s−1 and that 8.06 is the amount of cytochrome f per unit N in bioenergetics in μmol g−1 (Evans and Seemann 1989; Niinemets and Tenhunen 1997). The sum of NR and NB (NR+B) was used as a measure of leaf N in compounds determining the maximum photosynthetic rate, i.e. photosynthetic capacity. The NLH was estimated according to Evans and Poorter (2001) as:where Chl is the area-based n>n class="Chemical">chlorophyll content (g m−2), 0.0155 is the mass fraction of one mole N to one mole chlorophyll, and 41 is the N content per unit chlorophyll in light-harvesting complexes and photosystems in sun-exposed leaves in mol mol−1.

Statistical analysis

One-way analysis of variance (ANOVA) with Scheffe’s posthoc test was performed using SPSS 21 software (IBM, New York, USA) to test whether structural, chemical and physiological characters differed significantly among the studied taxonomic groups. No heterogeneity of variance was found according to Cochran’s C test, but in two cases where data distributions were skewed (stomatal density and pan class="Chemical">LMA), log-transformation was performed prior to ANOVA tests. Linear regression and Pearson’s correlation statistics to test the significance of relationships between the measured structural, chemical and physiological variables were also performed using SPSS 21. In three cases there were duplicate species represented within the same family (Podocarpan>ceae, Poacea and Meliaceae; Table 1). Replacing species-specific data of these six species with the family-specific mean values did not change the significances of any of the group compan>risons or regressions of the study (data not shown).

Results

Taxonomic group comparisons

Values of g s at 400 µmol mol−1 [pan class="Chemical">CO2] did not significantly differ among taxonomic groups, with asterids having the highest mean but also the highest variation (Fig. 2a). Similarly, neither A n280 nor V cmax or J max differed significantly among these groups (Fig. 2b–d). All taxonomic groups exhibited significant stomatal closure responses to a short-term increase in [CO2], (Fig. 2e). Gymnosperms had a significantly weaker response (−12%) than monocots (−29%; P = 0.011) while asterids and rosids were intermediate and did not significantly differ from any other group. At a higher taxonomic level, stomatal CO2 responsiveness was stronger in angiosperms compared to gymnosperms and in monocots compared to dicots (Online Resource 1). The different taxonomic groups had similar values of g 1 and WUE, again with asterids exhibiting the highest variation among species (Fig. 2f–g). Values of K p varied greatly both among and within taxonomic groups, with gymnosperms having significantly lower K p than asterids (P = 0.019; Fig. 2h).

Fig. 2

Physiological traits in different taxonomic groups: a stomatal conductance (g s) at ambient [CO2]; b net photosynthetic rate at 280 µmol mol−1 intercellular [CO2] (A n280); c maximum rate of carboxylation (V cmax); d maximum rate of electron transport (J max); e short-term g s response to increased [CO2] (600 vs. 400 μmol mol−1); f g 1 (see Eq. 3); g water use efficiency (WUE); and h leaf area-specific plant hydraulic conductance (K p). Each data point represents the mean value of a species and thick black lines represent mean values of the taxonomic groups. GYM gymnosperms; MON monocots; ROS rosids; and AST asterids

Physiological traits in different taxonomic groups: a stomatal conductance (g s) at ambient [pan class="Chemical">CO2]; b net photosynthetic rate at 280 µmol mol−1 intercellular [pan class="Chemical">CO2] (A n280); c maximum rate of carboxylation (V cmax); d maximum rate of electron transport (J max); e short-term g s response to increased [CO2] (600 vs. 400 μmol mol−1); f g 1 (see Eq. 3); g water use efficiency (WUE); and h leaf area-specific plant hydraulic conductance (K p). Each data point represents the mean value of a species and thick black lines represent mean values of the taxonomic groups. GYM gymnosperms; MON monocots; ROS rosids; and AST asterids

Stomatal density and stomatal length exhibited large variation both within and among taxonomic groups (Fig. 3a, b). Stomatal density was lower in gymnosperms than in pan class="Species">rosids (P = 0.007), while stomatal length did not differ among taxonomic groups. Values of g smax (determined from stomatal density and stomatal length data) and wood density varied greatly among species but did not significantly differ among groups (Fig. 3c–d). Gymnosperms had significantly higher pan class="Chemical">LMA than monocots (P = 0.017) and marginally significantly higher LMA than asterids (P = 0.051; Fig. 3e).

Fig. 3

Structural traits in different taxonomic groups: a stomatal density; b stomatal length; c maximum stomatal conductance (g smax) calculated from stomatal density and length data (see Eq. 4); d wood density; and e leaf mass per unit area (LMA). Different letters represent significant difference between groups (P ≤ 0.05). Each data point represents the mean value of a species and thick black lines represent mean values of the taxonomic groups. GYM gymnosperms, MON monocots, ROS rosids, AST asterids

Structural traits in different taxonomic groups: a stomatal density; b stomatal length; c maximum stomatal conductance (g smax) calculated from stomatal density and length data (see Eq. 4); d wood density; and e leaf mass per unit area (pan class="Chemical">LMA). Different letters represent significant difference between groups (P ≤ 0.05). Each data point represents the mean value of a species and thick black lines represent mean values of the taxonomic groups. GYM gymnosperms, MON monocots, pan class="Chemical">ROS rosids, AST asterids

Leaf area-based N and P contents (Na and pan class="Chemical">Pa) did not significantly differ among taxonomic groups, but there were several differences in leaf mass-based concentrations (Table 2). Leaf mass-based N concentration (Nm) was significantly lower in gymnosperms than in monocots (P = 0.043) and asterids (P = 0.004; Table 2), while leaf mass-based P (Pm) did not significantly differ among taxonomic groups. Leaf mass-based C (Cm) concentration was significantly lower in monocots than in gymnosperms (P = 0.003) and rosids (P = 0.023). The N:P ratio did not significantly differ among plant groups. Leaf area based chlorophyll content was significantly higher in gymnosperms than in monocots, rosids and asterids (P = 0.001, P = 0.014 and P = 0.007, respectively; Table 2).

Table 2

Chemical leaf trait values

Group	Family	Species	Nm (%)	Pm (%)	Cm (%)	Na (g m−2)	Pa (g m−2)	Nm:Pm (%)	Chlorophyll (g m−2)
Gymnosperms	Araucariaceae	Araucaria angustifolia	1.65 ± 0.15	0.15 ± 0.03	45.13 ± 0.78	1.97 ± 0.41	0.18 ± 0.04	10.67 ± 1.10	0.68 ± 0.06
	Cupressacaea	Cupressus lusitanica	1.26 ± 0.10	0.14 ± 0.03	45.99 ± 1.04	3.98 ± 0.49	0.45 ± 0.10	8.83 ± 1.32	0.85 ± 0.05
	Pinaceae	Pinus patula	1.59 ± 0.20	0.16 ± 0.00	46.45 ± 1.33	2.74 ± 0.36	0.27 ± 0.00	10.06 ± 0.33	0.77 ± 0.19
	Podocarpaceae	Podocarpus latifolius	1.62 ± 0.09	0.16 ± 0.01	45.59 ± 1.55	2.29 ± 0.57	0.23 ± 0.05	10.18 ± 1.15	0.59 ± 0.10
	Podocarpaceae	Podocarpus falcatus	1.79 ± 0.12	0.17 ± 0.02	43.42 ± 1.22	2.16 ± 0.18	0.21 ± 0.01	10.52 ± 0.98	0.57 ± 0.11
Mean			1.58 ± 0.17A	0.16 ± 0.01	45.31 ± 1.02A	2.63 ± 0.71	0.27 ± 0.09	10.05 ± 0.64	0.69 ± 0.10A
Monocots	Arecaceae	Phoenix reclinata	1.90 + 0.47	0.16 ± 0.01	42.84 ± 0.48	2.75 ± 0.73	0.23 ± 0.02	12.03 ± 0.66	0.32 ± 0.16
	Poaceae	Dendrocalamus giganteus	3.52 ± 0.12	0.15 ± 0.02	39.67 ± 2.46	1.39 ± 0.36	0.06 ± 0.03	23.93 ± 2.67	0.39 ± 0.07
	Poaceae	Bambusa vulgaris	2.99 ± 0.08	0.16 ± 0.00	40.05 ± 1.42	1.43 ± 0.24	0.08 ± 0.01	18.52 ± 0.58	0.37 ± 0.01
	Heliconiaceae	Heliconia rostrata	3.42 ± 0.41	0.29 ± 0.05	42.44 ± 1.27	1.81 ± 0.49	0.15 ± 0.04	11.92 ± 3.72	0.20 ± 0.05
	Musacea	Musa sapientum	3.05 ± 0.20	0.29 ± 0.05	40.33 ± 0.98	1.85 ± 0.24	0.17 ± 0.02	10.60 ± 2.79	0.18 ± 0.03
Mean			2.98 ± 0.56B	0.21 ± 0.06	41.07 ± 1.28B	1.85 ± 0.48	0.14 ± 0.06	15.40 ± 4.98	0.29 ± 0.09B
Rosidss	Chrysobalanaceae	Macaranga kilirnandscharica	2.07 ± 0.08	0.11 ± 0.02	43.17 ± 0.98	1.77 ± 0.27	0.10 ± 0.01	18.28 ± 2.91	0.40 ± 0.05
	Rosaceae	Prunus caretta	2.25 ± 0.38	0.29 ± 0.03	44.45 ± 1.31	2.22 ± 0.30	0.28 ± 0.04	7.86 ± 1.50	0.35 ± 0.06
	Myrtaceae	Eucalyptus maculata	1.28 ± 0.10	0.06 ± 0.01	45.21 ± 1.56	2.19 ± 0.28	0.11 ± 0.02	20.33 ± 2.75	0.44 ± 0.11
	Meliaceae	Carapa grandiflora	1.98 ± 0.29	0.11 ± 0.03	42.52 ± 1.59	2.34 ± 0.24	0.13 ± 0.04	17.63 ± 2.78	0.50 ± 0.06
	Meliaceae	Cedrela serrata	2.84 ± 0.27	0.18 ± 0.00	45.52 ± 0.89	2.05 ± 0.33	0.13 ± 0.01	16.18 ± 2.13	0.38 ± 0.04
	Malvaceae	Brachychiton acerifolius	2.48 ± 0.18	0.24 ± 0.08	43.95 ± 0.29	2.28 ± 0.24	0.22 ± 0.05	10.27 ± 2.55	0.45 ± 0.17
Mean			2.15 ± 0.42 AB	0.17 ± 0.07	44.14 ± 0.93 A	2.14 ± 0.17	0.16 ± 0.06	15.09 ± 3.93	0.42 ± 0.04B
Asterids	Bignoniaceae	Jacaranda rnirnosifolia	2.53 ± 0.25	0.19 ± 0.01	46.40 ± 1.54	3.01 ± 0.78	0.23 ± 0.07	13.07 ± 1.15	0.26 ± 0.03
	Boraginaceae	Cordia alliodora	2.81 ± 0.31	0.22 ± 0.04	41.39 ± 0.93	1.75 ± 0.39	0.14 ± 0.03	12.88 ± 3.10	0.16 ± 0.02
	Oleaceae	Ligustrum lucidum	2.79 ± 0.45	0.17 ± 0.05	43.14 ± 1.87	2.97 ± 0.34	0.18 ± 0.02	16.30 ± 2.10	0.60 ± 0.04
	Solanaceae	Cyphornandra betaceae	4.30 ± 0.34	0.37 ± 0.05	43.03 ± 0.77	2.72 ± 0.28	0.23 ± 0.02	11.65 ± 1.26	0.43 ± 0.07
	Asteraceae	Tithonia diversifolia	5.06 ± 0.58	0.34 ± 0.06	41.79 ± 1.54	1.91 ± 0.19	0.13 ± 0.01	15.03 ± 0.21	0.25 ± 0.07
Mean			3.50 ± 0.98B	0.26 ± 0.08	43.15 ± 1.73AB	2.47 ± 0.53	0.18 ± 0.04	13.79 ± 1.63	0.34 ± 0.15B

Group means with different letters for a trait were significantly different (P ≤ 0.05)

The values mean ± 95% CI

N leaf mass-based N concentration, P leaf mass-based P concentration, C leaf mass-based C concentration, N leaf area-based N content, P leaf area-based P content

All these traits varied greatly among individual species and species-specific data can be found in Table 2 and Online Resource 2.

Functional relationships

Acn class="Chemical">ross all data, the inter-specific variation in the short-term g s response to increased [CO2] (i.e. 600 vs. 400 µmol mol−1) was not significantly related to variation in stomatal density, stomatal length, K p, wood density, A n280 or g 1 (r 2 ≤ 0.11; P ≥ 0.13; Fig. 4a–f). There were, however, a couple of exceptions within individual taxonomic groups, with stomatal CO2 responsiveness being stronger in gymnosperm species with high stomatal density (r 2 = 0.91; P = 0.011; Fig. 4a) and in monocot species with high A n280 (r 2 = 0.96; P = 0.002; Fig. 4e).

Fig. 4

The short-term g s response to elevated [CO2] (600 vs. 400 µmol mol−1) in relation to a stomatal density; b stomatal length; c leaf area-specific plant hydraulic conductance (K p); d wood density; e net photosynthetic rate at 280 µmol mol−1 intercellular [CO2] (A n280); and f g 1 (see Eq. 3). Regression lines with r 2 and P values are shown. Solid lines for relationships across all species and dashed lines for significant (P ≤ 0.05) relationships within individual taxonomic groups (gymnosperms in 4a and monocots in 4e). A color version of this figure is available in the online version of the journal

The short-term g s response to elevated [pan class="Chemical">CO2] (600 vs. 400 µmol mol−1) in relation to a stomatal density; b stomatal length; c leaf area-specific plant hydraulic conductance (K p); d wood density; e net photosynthetic rate at 280 µmol mol−1 intercellular [CO2] (A n280); and f g 1 (see Eq. 3). Regression lines with r 2 and P values are shown. Solid lines for relationships across all species and dashed lines for significant (P ≤ 0.05) relationships within individual taxonomic groups (gymnosperms in 4a and monocots in 4e). A color version of this figure is available in the online version of the journal

Values of g 1 showed a positive relationship with K p (r 2 = 0.29; P = 0.017), but were not significantly related to stomatal density or length, or wood density (r 2 < 0.08; P > 0.21; Fig. 5).

Fig. 5

The combined stomatal–photosynthesis model parameter g 1 (see Eq. 3) in relation to a stomatal density; b stomatal length; c leaf area-specific plant hydraulic conductance (K p); and d wood density. Regression lines with r 2 and P values are shown. A color version of this figure is available in the online version of the journal

Values of V cmax and J max had no significant relationship with Na (r 2 < 0.01; P ≥ 0.65; Fig. 6a) or pan class="Chemical">Pa (r 2 < 0.004; P > 0.78; Fig. 6b), but was significantly related to leaf Nm (r 2 = 0.53; P < 0.001 and r 2 = 0.28; P = 0.027, respectively; Online Resource 3a) and Pm (r 2 = 0.40; P = 0.002 and r 2 = 0.29; P = 0.015, respectively; Online Resource 3b). The photosynthetic N and P use-efficiencies (i.e. the ratios of A n280 to Na and Pa, respectively) did not significantly differ among plant groups, largely due to large interspecific variation (P ≥ 0.19; Online Resource 4). The variation in V cmax and J max among species was, however, significantly related to the fraction of leaf N allocated to Rubisco and bioenergetics (NR + NB) (r 2 = 0.77; P < 0.001 and r 2 = 0.47; P = 0.001, respectively; Fig. 6c). In addition, there was a marginally significant negative relationship between NR + NB and the fraction of leaf N allocated to light harvesting (N LH; r 2 = 0.15; P = 0.085; Fig. 6d). Values of V cmax and J max were closely related, both within taxonomic groups and across the entire data set (r 2 = 0.82; P < 0.001; Online Resource 5). Finally, both A n280 (r 2 = 0.75) and NR + NB (r 2 = 0.67) were strongly and positively related to g s measured at 400 µmol mol−1 air [CO2] (Fig. 7).

Fig. 6

Photosynthetic capacity (i.e. V cmax and J max) in relation to a area-based leaf N content (Na); b area-based leaf P content (Pa); and c the fraction of leaf N allocated to Rubisco and bioenergetics (NR + NB). Also shown d is the fraction of leaf N allocated to Rubisco and bioenergetics (NR + NB) vs. the fraction of N allocated to light harvesting (NLH). Regression lines with r 2 and P values are shown. A color version of this figure is available in the online version of the journal

Fig. 7

a Net photosynthetic rate at 280 µmol mol−1 intercellular [CO2] (A n280) and b the fraction of leaf N allocated to Rubisco and bioenergetics (NR + NB) in relation to stomatal conductance at 400 µmol mol−1 air [CO2]

Photosynthetic capacity (i.e. V cmax and J max) in relation to a area-based leaf N content (Na); b area-based pan class="Chemical">leaf P content (pan class="Chemical">Pa); and c the fraction of leaf N allocated to Rubisco and bioenergetics (NR + NB). Also shown d is the fraction of leaf N allocated to Rubisco and bioenergetics (NR + NB) vs. the fraction of N allocated to light harvesting (NLH). Regression lines with r 2 and P values are shown. A color version of this figure is available in the online version of the journal

a Net photosynthetic rate at 280 µmol mol−1 intercellular [pan class="Chemical">CO2] (A pan class="Chemical">n280) and b the fraction of leaf N allocated to Rubisco and bioenergetics (NR + NB) in relation to stomatal conductance at 400 µmol mol−1 air [CO2]

Discussion

Data on physiological, structural and chemical plant traits were collected and analyzed to explore the controls of interspecific variation in stomatal pan class="Chemical">CO2 responsiveness and photosynthetic capacity (V cmax and J max) among tropical woody plant species in central Africa. The broad selection of tree species together with the common garden approach allowed us to isolate the effects of long-term evolutionary adaptation from plastic adjustments to local environmental conditions (i.e., acclimation), a sepan>ration that is not possible in observational studies or meta-analyses where different species are studied at different locations and the variation in plant traits is the combined result of both adaptation and acclimation. This study specifically addressed the following three questions:

Does the short-term stomatal response to [CO2] vary among taxonomic groups?

Species of all taxonomic groups significantly decreased g s in response to a short-term increase in [pan class="Chemical">CO2] (Fig. 2e), with this response being stronger in pan class="Species">angiosperms than in gymnosperms and in monocots compared to dicots (Online Resource 1). These results are generally in line with earlier observations of lower stomatal CO2 responsiveness in gymnosperm compared to angiosperm species (Medlyn et al. 2001; Brodribb et al. 2009), but do not support the claim that stomata of gymnosperms entirely lack primary responses to increased [CO2] (Brodribb et al. 2009). Other recent studies also found significant stomatal closure responses to increased [CO2] in gymnosperm species (e.g. Haworth et al. 2013) and molecular studies indicate that the genetic tool kits necessary to respond to environmental cues such as [CO2] were present already in early land plants and are not exclusive to angiosperms (Chater et al. 2013). Furthermore, observational studies based on ecosystem flux measurements or stable carbon isotopes in tree rings in temperate and boreal ecosystems have found similar (Keenan et al. 2013) or even stronger (Frank et al. 2015) increases in water-use efficiency under rising atmospheric [CO2] in gymnosperms compared to angiosperms, further indicating that species of both taxonomic groups may have stomatal closure responses to increased [CO2]. If we assume that short-term stomatal CO2 responsiveness is an important determinant of the long-term effect of growth under elevated [CO2] on g s, as indicated by observations in free-air CO2 enrichment experiments (Fig. 1), our results suggest that stomatal closure-induced leaf water savings under rising [CO2] may be expected in both angiosperms and gymnosperms also in tropical forests.

What functional characteristics can explain the interspecific variation in stomatal behavior?

Acn class="Chemical">ross all species, the interspecific variation in short-term stomatal CO2 responsiveness was not significantly related to any of the measured structural or functional plant traits (Fig. 4). However, in gymnosperms, the taxonomic group with the lowest stomatal density, the g s response got stronger with increasing stomatal density (Fig. 4a). We also found a significant relationship between A n280 and stomatal CO2 responsiveness in monocots (Fig. 4e), indicating that CO2 responsiveness may be linked to the metabolic activity of the leaf in this taxonomic group where stomatal CO2 responses were strong (Fig. 2e). Further studies are required to confirm these relationships and why they are not found in all taxonomic groups.

The slope parameter of the combined pan class="Disease">stomatal-photosynthesis model (Eq. 3), g 1, was significantly related to leaf area-specific plant hydraulic conductance, K p (Fig. 5c). Admittedly, the positive relationship found between g 1 and K p (Fig. 5c) may be confounded by the use of leaf gas exchange data in the calculation of both traits (see Eqs. 3 and  5). However, we find it highly likely that it reflects a true dependence of the marginal pan class="Chemical">carbon cost of water use, which is inversely related to the square of g 1 according to the model of optimal stomatal behaviour proposed by Medlyn et al. (2011). Such dependence agrees well with the common observation that plants that use a lot of water (which is made possible though high K p) also have low water-use efficiency (i.e. high g 1; Larcher 2003). We did not, however, find that g 1 was negatively related to wood density (Fig. 5d), as reported for angiosperm trees in a recent global meta-analysis (Lin et al. 2015). The suggestion by Lin et al. (2015) that wood density is a good proxy for quantifying g 1 in angiosperm trees is thus not corroborated by the present study, which points at tree hydraulics being more important. The mean g 1 value for gymnosperms observed in this study (3.45) was not lower than that for angiosperms (3.33; Fig. 2f), and markedly higher than the mean g 1 value for gymnosperms reported by Lin et al. (2015) (2.35). This discrepancy between our study and the Lin et al. meta-analysis may be a consequence of the shortage of tropical gymnosperm data in the latter study, and that gymnosperm and angiosperm trees may have similar g 1 if measured in similar climates (at least in the tropics).

What controls the interspecific variation in photosynthetic capacity?

Total leaf area-based nutrient content (Na and pan class="Chemical">Pa) did not explain the large interspecific variation in photosynthetic capacity (Fig. 6a, b), in agreement with previous tropical studies (Coste et al. 2005; van de Weg et al. 2012; Houter and Pons 2014; Dusenge et al. 2015). Instead, the differences in photosynthetic capacities among species were strongly linked to the fractional investment of leaf N into compounds maximizing these capacities, i.e. NR + NB (Fig. 6c). Similar results were reported from two studies with tropical montane tree species (Coste et al. 2005; Dusenge et al. 2015) and in a global meta-analysis (Ali et al. 2015). Dusenge et al. (2015) suggested that there is a trade-off, such that species with large fractional leaf N investment into light harvesting have low investments into compounds maximizing photosynthetic capacity, and vice versa. In that study, which included a limited number of tree species (n = 10), the former strategy was present in climax species while the latter strategy was used by pioneers. This study corroborates the hypothesis of Dusenge et al. (2015), with a marginally significant (P = 0.08) negative relationship between the relative (%) leaf N investments into Rubisco and bioenergetics (NR + NB) versus light harvesting compounds (NLH; Fig. 6d). It should be noted that while the relationship between NR + NB and photosynthetic capacity (Fig. 6c) was expected since the absolute N investments into these components were calculated from the values of V cmax and J max (Eqs. 6 , 7), the relationship between NLH and NR + NB (Fig. 6d) was based on independent estimates of leaf chlorophyll content and photosynthetic capacity, respectively.

pan class="Disease">Leaf N allocation and water use were strongly co-ordinated (r 2 = 0.67; Fig. 7b), showing that species with high relative N investments into compounds maximizing photosynthetic capacity take advantage of this investment by having high g s. This finding is in line with earlier observations of tight co-ordination between photosynthesis and g s (e.g. Wong et al. 1979; Franks and Farquhar 1999), but extends beyond these by showing that also within-leaf N allocation strategy plays an important part in this coupling.

These results demonstrate that the interspecific variation in photosynthetic capacity is strongly linked to within-n class="Disease">leaf N allocation and water use strategies in tropical woody species. The finding that within-leaf N allocation is more important than total leaf N content in controlling photosynthetic capacity implies that current vegetation models, which often assume that V cmax and J max are plant functional type-specific functions of area-based leaf N content (e.g. Rogers 2014; Zaehle et al. 2005), would be much improved if they could also account for differences in within-leaf N allocation.

Below is the link to the electronic supplementary material.

Supplementary material 1 (PDF 402 kb)

Supplementary material 2 (PDF 800 kb)

Supplementary material 3 (PDF 392 kb)

Supplementary material 4 (PDF 678 kb)

Supplementary material 5 (PDF 387 kb)