Plan
Comptes Rendus

Stratigrahy, sedimentology (Palaeoenvironment)
Estimating variation in stomatal frequency at intra-individual, intra-site, and inter-taxonomic levels in populations of the Leonardoxa africana (Fabaceae) complex over environmental gradients in Cameroon
[Estimation de la variation de la fréquence stomatique aux niveaux intra-individuel, intra-site, et inter-taxonomique chez des populations du complexe Leonardoxa africana (Fabaceae) selon des gradients environnementaux au Cameroun]
Comptes Rendus. Géoscience, Volume 345 (2013) no. 7-8, pp. 350-359.

Résumés

Variation of stomatal frequency (stomatal density and stomatal index) includes genetically-based, potentially-adaptive variation, and variation due to phenotypic plasticity, the degree of which may be fundamental to the ability to maintain high water-use efficiency and thus to deal with environmental change. We analysed stomatal frequency and morphology (pore length, pore width) in leaves from several individuals from nine populations of four sub-species of the Leonardoxa africana complex. The dataset represents a hierarchical sampling wherein factors are nested within each level (leaves in individuals, individuals in sites, etc.), allowing estimation of the contribution of different levels to overall variation, using variance-component analysis. SI showed significant variation among sites (“site” is largely confounded with “sub-species”), being highest in the sub-species localized in the highest-elevation site. However, most of the observed variance was accounted for at intra-site and intra-individual levels. This variance could reflect great phenotypic plasticity, presumably in response to highly local variation in micro-environmental conditions.

Supplementary Materials:
Supplementary material for this article is supplied as a separate file:

La variation de la fréquence stomatique (densité de stomates [SD] et index stomatique [SI]) inclut des variations potentiellement adaptatives et des variations liées à la plasticité phénotypique, dont le degré peut être fondamental dans la régulation et pour le maintien d’une haute « water-use efficiency » et donc pour faire face aux changements environnementaux. Nous avons analysé les fréquences et la morphologie (longueur et largeur des pores) de stomates de feuilles provenant de neuf populations de quatre sous-espèces du complexe Leonardoxa africana. Le jeu de données représente un échantillonnage hiérarchisé dans lequel les facteurs sont emboîtés à chaque niveau (feuilles dans individus, individus dans site/populations, etc.). En utilisant une « variance-component analysis », on a estimé la contribution de chaque niveau à la variation globale. SI est significativement différent entre les sites et atteint les valeurs les plus élevées dans le site/sous-espèce localisé à la plus haute altitude. Néanmoins, la plus grande variance est enregistrée aux niveaux intra-sites et intra-individus. Cela pourrait être le reflet d’une grande plasticité phénotypique qui répond à une forte variation locale des conditions micro-environnementales.

Compléments :
Des compléments sont fournis pour cet article dans le fichier séparé :

Métadonnées
Reçu le :
Accepté le :
Publié le :
DOI : 10.1016/j.crte.2013.06.003
Keywords: Stomatal frequency, Phenotypic plasticity, Variance-component analysis, Leonardoxa, Central Africa
Mot clés : Fréquence stomatique, Plasticité phénotypique, Variance-component analysis, Leonardoxa, Afrique centrale

Walter Finsinger 1 ; Thibaut Dos Santos 2 ; Doyle McKey 2, 3

1 UMR 5059 CNRS/UM2/EPHE, Centre de bio-archéologie et d’écologie (CBAE), Institut de Botanique, 163, rue A.-Broussonnet, 34090 Montpellier, France
2 UMR 5175 CNRS, Centre d’écologie fonctionnelle et évolutive (CEFE), 1919, route de Mende, 34293 Montpellier cedex 05, France
3 Institut universitaire de France, 103, boulevard Saint-Michel, 75005 Paris, France
@article{CRGEOS_2013__345_7-8_350_0,
     author = {Walter Finsinger and Thibaut Dos Santos and Doyle McKey},
     title = {Estimating variation in stomatal frequency at intra-individual, intra-site, and inter-taxonomic levels in populations of the {\protect\emph{Leonardoxa} africana} {(Fabaceae)} complex over environmental gradients in {Cameroon}},
     journal = {Comptes Rendus. G\'eoscience},
     pages = {350--359},
     publisher = {Elsevier},
     volume = {345},
     number = {7-8},
     year = {2013},
     doi = {10.1016/j.crte.2013.06.003},
     language = {en},
}
TY  - JOUR
AU  - Walter Finsinger
AU  - Thibaut Dos Santos
AU  - Doyle McKey
TI  - Estimating variation in stomatal frequency at intra-individual, intra-site, and inter-taxonomic levels in populations of the Leonardoxa africana (Fabaceae) complex over environmental gradients in Cameroon
JO  - Comptes Rendus. Géoscience
PY  - 2013
SP  - 350
EP  - 359
VL  - 345
IS  - 7-8
PB  - Elsevier
DO  - 10.1016/j.crte.2013.06.003
LA  - en
ID  - CRGEOS_2013__345_7-8_350_0
ER  - 
%0 Journal Article
%A Walter Finsinger
%A Thibaut Dos Santos
%A Doyle McKey
%T Estimating variation in stomatal frequency at intra-individual, intra-site, and inter-taxonomic levels in populations of the Leonardoxa africana (Fabaceae) complex over environmental gradients in Cameroon
%J Comptes Rendus. Géoscience
%D 2013
%P 350-359
%V 345
%N 7-8
%I Elsevier
%R 10.1016/j.crte.2013.06.003
%G en
%F CRGEOS_2013__345_7-8_350_0
Walter Finsinger; Thibaut Dos Santos; Doyle McKey. Estimating variation in stomatal frequency at intra-individual, intra-site, and inter-taxonomic levels in populations of the Leonardoxa africana (Fabaceae) complex over environmental gradients in Cameroon. Comptes Rendus. Géoscience, Volume 345 (2013) no. 7-8, pp. 350-359. doi : 10.1016/j.crte.2013.06.003. https://comptes-rendus.academie-sciences.fr/geoscience/articles/10.1016/j.crte.2013.06.003/

Version originale du texte intégral

1 Introduction

The need to understand how species respond to environmental variation has become critical as environmental changes have demonstrable ecological effects in many natural systems (Parmesan, 2006; Walther et al., 2002). Determining how individuals base key life-history decisions on environmental cues is important to predict how environmental changes will influence them. Phenotypic plasticity, defined as the ability of a single genotype to modify its phenotype under heterogeneous environmental conditions (Houston and McNamara, 1992), is fundamental to a plant's ability to cope with environmental change at various temporal and spatial scales.

Stomata regulate the permeability of leaves of terrestrial plants to gases (gas phase conductance), so that the leaves can absorb CO2 for photosynthesis without losing excessive water vapour (Raven, 2002). There is at present sufficient evidence to suggest that perennial terrestrial plants regulate their gas exchange (and thus CO2 uptake and water loss) by varying the opening of the stomatal pores (Young et al., 2006) and by changing the stomatal frequency on leaves (Gagen et al., 2011; Hetherington and Woodward, 2003). Thus, to optimize their resource use and evade costs of having excess stomata, plants could be expected to adjust stomatal frequency (Roth-Nebelsick, 2007). Amongst several environmental factors, atmospheric CO2 has often been observed to have a strong influence in the regulation of stomatal frequency (Chen et al., 2001; Hetherington and Woodward, 2003; Wagner et al., 1996; Woodward, 1987). However, most scholars agree with the view that responses are species-specific (Marchi et al., 2004; Tognetti et al., 2000) and include cases in which species do not respond at all or respond to increasing atmospheric CO2 by increasing, rather than decreasing, stomatal frequency (Hetherington and Woodward, 2003).

Although strong plastic (developmental) responses have been observed in stomatal frequency (stomatal density and stomatal index [SI]) (Royer, 2001), there is clear evidence of genetic control of the CO2 response (Jordan, 2011). Thus, in Arabidopsis thaliana, the HIC gene codes for CO2 responses in both stomatal density and index (Gray et al., 2000). The latter evidence strongly suggests that stomatal frequencies are determined (at least to some extent) by the genotype and implies that plants can regulate their stomatal frequency by phenotypic plasticity, allowing them to adapt to changing environmental conditions. For example, it is often observed that plants have higher stomatal density when growing at lower CO2 concentrations (Hetherington and Woodward, 2003; Woodward, 1987) or at higher elevation (where the partial pressure of CO2 is presumably lower) (McElwain, 2004). The negative relationship between CO2 concentration (or partial pressure) and stomatal frequency has often been used in palaeoecological investigations to reconstruct temporal changes of atmospheric CO2 concentrations using (sub)fossil leaf cuticles collected from sedimentary archives (Beerling, 1999; Beerling and Chaloner, 1993; McElwain and Chaloner, 1995; McElwain et al., 2002; Steinthorsdottir et al., 2013). However, such reconstructions are based on the assumption that environmental changes at smaller scales (i.e., among individuals or within individuals of the same population) are negligible compared to the stomatal-frequency changes caused by varying CO2 concentrations (Finsinger and Wagner-Cremer, 2009). Although estimates of the stomatal-frequency variation at intra-population and intra-individual level might be necessary to fully account for uncertainties in paleo-CO2 reconstructions, the analysis of modern leaves and of these sources of variation is rarely conducted (Chen et al., 2001).

Here we aimed to compare stomatal frequency within the Leonardoxa africana complex (Fabaceae–Caesalpinioideae) at different hierarchical levels: among sub-species, among populations belonging to the same sub-species, among individuals of the same population, and within individuals. We were particularly interested in investigating at which hierarchical level the greatest changes in the SI were observed and in illustrating the role of environmental variables as potential predictors for cuticle-morphology characteristics (stomatal frequency and stomatal pore size). Our investigation focuses principally on variations of SI because this variable allows a more appropriate comparison among leaves, since the SI is less biased by external factors than is stomatal density (Beerling, 1999).

2 Material and methods

Leonardoxa africana (Baill.) Aubréville (Fabaceae: Caesalpinioideae) is a species complex native to Atlantic central Africa (Fig. 1) that comprises four distinct, mostly allopatric subspecies (L. a. africana, L. a. letouzeyi, L. a. gracilicaulis, and L. a. rumpiensis). These have been described and identified based on morphological characters, particularly those related to the plants’ interactions with ants (McKey, 2000).

Fig. 1

Location of the studied populations and contemporary distributions of Leonardoxa africana subspecies.

Emplacement des populations étudiées et distribution actuelle des sous-espèces de Leonardoxa africana.

Adapted from McKey, 2000.

These four taxa of small- to medium-sized understory trees are located in lower Guinea forests from extreme south-eastern Nigeria to Gabon (Fig. 1). Leonardoxa africana ssp. gracilicaulis McKey is a non-myrmecophyte involved in loose associations with opportunistic ants attracted to foliar nectaries that is most common on hills (McKey, 2000). Leonardoxa africana ssp. rumpiensis McKey has hollow internodes, inhabited by a variety of opportunistic ants (Chenuil and McKey, 1996; McKey, 2000). It appears to be very localized in submontane forests in the Rumpi Hills. Leonardoxa africana ssp. letouzeyi McKey is found in lowland forests in extreme western Cameroon and south-eastern Nigeria. Although its juvenile individuals can shelter several species of arboricolous ants in their internodes, mature trees of L. a. letouzeyi are associated specifically with a single ant species, Aphomomyrmex afer Emery (Gaume and McKey, 1998; McKey, 1991, 2000). Subspecies rumpiensis and letouzeyi may be considered as morphologically and ecologically ‘transitional’ myrmecophytes (Brouat and McKey, 2000; Fiala et al., 1999). The fourth sub-species of the complex, L. africana ssp. africana, is morphologically the most specialized myrmecophyte. It occurs in lowland coastal forests in southern Cameroon. Leonardoxa a. africana is associated, as early as the seedling stage, with one specific mutualist ant species, Petalomyrmex phylax Snelling (McKey, 1984, 2000). All four subspecies form pinnately compound leaves composed of two to four leaflet pairs.

2.1 Sampling design, leaf collection and cuticle analysis

Mature leaves of the four Leonardoxa subspecies were collected and dried at room temperature. All leaves were labelled with the subspecies name, the population, and the individual to which they belonged. We analysed two sets of leaves of Leonardoxa to investigate variation of cuticle-morphology characteristics at different hierarchical levels. The larger set comprised at least two leaves (sometimes three) for each individual growing in the sampled populations (Table 1). In the smaller set (hereafter named the ‘leaflet-comparison dataset’), all leaflets from one compound leaf were analysed for four individuals per species (Table 2). For the leaflet-comparison dataset (which included leaves of L. a. africana and L. a. gracilicaulis with three leaflet pairs), leaflet position within the leaves (proximal, middle, distal leaflet) was also noted.

Table 1

Moyennes et écart-types des paramètres de morphologie cuticulaire mesurées sur des feuilles provenant de différentes populations des sous-espèces du complexe L. africana.

L. a. gracilicaulis L. a. letouzeyi L. a. africana L. a. rumpiensis
Population Boundé (84 m) Ebogo (669 m) Mt Kala (931 m) Iboti (1009 m) Nguti (229 m) Korup (351 m) Boundé (84 m) Ebodjé (86 m) Dikome (1148 m)
SD (# mm−2) 243.3 ± 63.4 298.3 ± 108.3 239.6 ± 75.3 236 ± 61.7 233.9 ± 82.4 288.5 ± 102.4 229.5 ± 63.2 259.2 ± 72.4 354.4 ± 97.2
ED (# mm−2) 3446.6 ± 437.6 3339.2 ± 404.0 3233.8 ± 483.2 3533.3 ± 471.4 3492.7 ± 696.0 3071.5 ± 445.4 3200.9 ± 392.3 3116.0 ± 423.8 3531.8 ± 455.6
SI (%) 6.7 ± 1.9 8.3 ± 3.0 6.9 ± 2.2 6.3 ± 1.6 6.3 ± 1.9 8.6 ± 2.9 6.7 ± 1.8 7.7 ± 2.1 9.3 ± 3.1
PL (μm) 13.5 ± 1.2 13.2 ± 1.2 14.3 ± 1.4 13.5 ± 1.3 13.6 ± 1.1 13.4 ± 1.3 13.2 ± 1.0 12.0 ± 1.2 13.8 ± 1.3
PW (μm) 7.4 ± 0.1 7.7 ± 1.27 8.9 ± 1.2 7.3 ± 1.4 7.4 ± 0.9 8.3 ± 1.1 7.5 ± 1.0 6.7 ± 1.0 7.6 ± 1.0
# individuals 10 10 10 8 10 10 10 10 8
# leaves sampled 25 25 22 16 23 25 21 21 18
Table 2

Moyennes et écart-types des paramètres de morphologie cuticulaire mesurées sur des feuilles provenant de différentes populations des sous-espèces du complexe L. africana.

L. a. africana L. a. gracilicaulis
Population Boundé (84 m) Ebodjé (86 m) Boundé (84 m)
Leaflet position Distal Middle Proximal Distal Middle Proximal Distal Middle Proximal
SD (# mm−2) 284 ± 47.7 312 ± 46.0 312 ± 22.8 293 ± 76.3 292 ± 55.4 285 ± 65.2 348.0 ± 71.8 272.0 ± 17.9 328.8 ± 86.2
ED (# mm−2) 3064 ± 188.3 2864.0 ± 28.5 3140 ± 104.9 3224 ± 325.8 3173.3 ± 241.6 3570 ± 705.8 3421.6 ± 411.2 3520.0 ± 232.8 3564 ± 354.8
SI (%) 8.49 ± 1.40 9.82 ± 1.38 9.04 ± 0.69 8.36 ± 2.08 8.49 ± 1.82 7.63 ± 2.14 9.25 ± 1.78 7.19 ± 0.63 8.42 ± 1.95
PL (μm) 12.4 ± 1.0 11.9 ± 0.5 12.4 ± 0.8 12.7 ± 0.9 12.6 ± 0.8 12.6 ± 1.0 13.6 ± 1.3 13.0 ± 0.9 13.1 ± 1.2
PW (μm) 6.6 ± 0.6 7.1 ± 0.7 6.8 ± 0.4 7.1 ± 0.8 6.9 ± 0.8 7.3 ± 0.7 7.7 ± 0.8 7.8 ± 0.8 7.5 ± 0.9

For cuticle analysis, a fragment of 0.5 cm2 was cut with a knife, bleached in a 5% solution of NaHClO2 at ∼60 °C for 5–6 h and then washed in distilled water. The bleaching caused the separation of the abaxial and adaxial leaf epidermis, allowing the stomata-bearing abaxial epidermis to be easily mounted on microscopic slides with glycerine jelly (Fig. 2). Cuticle analysis was performed using an Olympus BX 60 transmitted-light microscope (Hamburg, Germany) at ×500 enlargements connected to a Leica DFC 300 FX R2 video camera (Münster, Germany). The Leica Application Suite v3.7 image-analysis software allowed measurement, on count fields of 0.05 mm2, of stomatal density (SD) and epidermal-cell density (ED). Stomatal frequency (SI) was derived using the formula SI = (SD/[ED + SD] × 100) after Salisbury (1927). For each leaf, up to seven count fields were analysed and in addition 20 stomata were measured for pore length (PL) and pore width (PW), following Finsinger and Wagner-Cremer (2009). Two analysts performed the measurements of cuticle morphology, with the larger set being analysed by TDS (Dos Santos, 2011) and the smaller leaflet-comparison dataset analysed by WF.

Fig. 2

Images of leaf cuticles of the four subspecies of the Leonardoxa complex analysed in this study. A, B. Leonardoxa africana subsp. gracilicaulis from Ebogo (A) having higher SI and SD than leaf from Mt. Kala (B). C. Leonardoxa africana subsp. letouzeyi. D. Leonardoxa africana subsp. africana. E. Leonardoxa africana subsp. rumpiensis. F. Image of compound leaf of L. a. africana (color available online).

Image des cuticules de feuille des quatre sous-espèces de Leonardoxa analysées dans cette étude. A, B. Leonardoxa africana subsp. gracilicaulis provenant d’Ebogo (A) et de Mt. Kala (B). C. Leonardoxa africana subsp. letouzeyi. D. Leonardoxa africana subsp. africana. E. Leonardoxa africana subsp. rumpiensis. F. Image de feuille de L. a. africana.

2.2 Environmental variables

For each population, 20 environmental variables (Appendix S1) were extracted with the Diva-GIS software from the WorldClim database (Hijmans et al., 2005) that summarizes climate indices over a period of 50 years (AD 1950–2000). The dataset has a spatial resolution of 30 arc s (∼1 km) and allows obtaining accurate estimates of temperature and precipitation and derived climate indices (seasonality of temperature and precipitation, isothermality) in mountainous regions with large climatic gradients over small horizontal scales.

2.3 Numerical analyses

The dataset represents a hierarchical sampling wherein factors are nested within each category (populations in subspecies, individuals in populations, etc.). To investigate at which hierarchical level most of the variation in SI originates, we analysed the dataset using a variance-components analysis (Crawley, 2005, p. 638; Sims et al., 2006). We fitted to the L. africana SI dataset a linear mixed-effect model that included each of the categorical factors as a random effect plus a random error term (ɛ). The linear mixed model was fitted to these data using restricted maximum likelihood (REML) as implemented in the nlme v3.1–107 library (Pinheiro et al., 2013) in R v2.14 (R Core Development Team, 2012).

The relationships between climate variables are illustrated with a standardized principal component analysis (PCA; based on a correlation matrix with focusing on inter-subspecies correlations and standardizing subspecies scores by standard deviations, with centering and standardization by subspecies). This choice is particularly appropriate if the variables are measured in different units, e.g. when the data are actually environmental variables (Leps and Šmilauer, 2003). The relationship between stomata data (SI, PL and PW) and climate variables was illustrated with PCA (based on a covariance matrix with focusing on inter-subspecies correlations, subspecies scores not divided by standard deviation, centering by subspecies) and climate data added passively. All ordinations were calculated in Canoco 4.5 (Ter Braak and Šmilauer, 2003).

3 Results

Means and standard deviations of SD, ED, PL, and PW are listed in Tables 1 and 2 by subspecies, population, and leaflet position within compound leaves. In the nine populations of the four subspecies (Table 1), SD ranged between 229 and 354 stomata mm−2. All populations had means lower than 300 stomata mm−2 except that of L. a. rumpiensis at Dikome. Mean values of ED ranged between ∼3000 and ∼3500 cells mm−2 without notable extreme values. Similarly to mean SD, mean SI values were relatively homogeneous for most of the populations (with mean SI values ranging from 6 to 7%) but three populations had higher values: L. a. gracilicaulis at Ebogo (8.3%), L. a. africana at Ebodjé (7.7%), and L. a. rumpiensis at Dikome (9.3%). The variability of SI values, as estimated based on the standard deviations, accounts for 20–30% of the mean SI for all populations.

Mean values of S.D., ED, and SI measured on different leaflets of compound leaves from L. a. africana and L. a. gracilicaulis from Boundé (Table 2) differed from those listed in Table 1, likely indicating an observer effect. For SI, leaflet position was not a statistically significant factor for L. a. africana (ANOVA, p = 0.298), but was significant for L. a. gracilicaulis (ANOVA, p = 0.0489) (Fig. 3).

Fig. 3

Boxplot illustrating the similarity between SI values in different leaflets (D: distal, M: median, P: proximal) of L. a. africana (LA from populations Boundé and Ebodjé) and L. a. gracilicaulis (LG from Boundé).

Boxplot illustrant la similarité entre les valeurs de SI sur différentes folioles (D : distal, M : médian, P : proximal) de L. a. africana (LA provenant des populations Boundé et Ebodjé) et L. a. gracilicaulis (LG provenant de Boundé).

SI and SD were strongly positively correlated in both datasets (r = 0.86–0.88; p < 0.001), indicating that variation in SI and SD was mainly due to differences in stomatal initiation rather than in the lateral expansion of epidermal cells.

The estimated variance components for SI were sometimes higher than the residual variance (Table 3). Variance among leaflets was always smaller than the intra-site (among-individual) variance (∼10% vs ∼40–50% of the variance, respectively), reflecting the marked difference of SI among individuals in comparison to the differences among leaflets of the same individual. Intra-site variance was mostly in the order of magnitude of residual variance except for L. a. gracilicaulis (51% vs 36% of the variance). Also in the full dataset (where among-leaflet variance was not taken into account), the intra-site variance was in most cases larger than variance at other hierarchical levels.

Table 3

Pourcentage de variance expliquée par niveau hiérarchique pour l’indice stomatique (SI).

Inter-taxon Inter-site Intra-site Intra-individual Inter-leaflet Residual % variance
Leaflet dataset
 All 8.0 × 10−7 8.1 × 10−9 45.6 / 12.0 42.4
L. a. africana / 1.4 × 10−6 47.9 / 10.2 41.9
L. a. gracilicaulis / / 51.1 / 12.9 36.0
Full dataset
 All 6.2 × 10−6 12.1 40.6 18.9 / 28.3
L. a. letouzeyi 30.2 41.5 7.4 / 20.9
L. a. africana 8.3 17.0 39.4 / 35.3
L. a. gracilicaulis / 9.0 39.6 19.4 / 31.9
L. a. rumpiensis / / 55.6 22.4 / 22.0

The first two Principal Component Analysis (PCA) axes of the climate data (Fig. 4a) for each of the nine populations explained ∼84% of the variance in the climate data (46 and 38%, respectively). The climate variables with the highest loading on the first PCA axis were altitude and mean annual temperature, whereas seasonality of temperature and precipitation, and isothermality had higher loadings on the second PCA axis.

Fig. 4

A. PCA illustrating main patterns of variation of climate indices at the sites where the nine studied populations of the L. africana complex occur. B. PCA illustrating the relationships between cuticle-morphology variables (SI, PL, and PW) in individuals of L. africana and climate indices (isothermality, altitude, seasonality of T and P, and mean annual P). Climate variables were added passively in this PCA (color available online).

A. ACP illustrant les principales tendances de variation des indices climatiques selon les sites où se trouvent les neuf populations étudiées du complexe L. africana. B. ACP illustrant les relations entre les variables de morphologie cuticulaire (SI, PL, PW) chez les individus de L. africana et les indices climatiques (isothermalité, altitude, saisonalité de la température et des précipitations, et précipitation annuelle moyenne). Les variables climatiques ont été additionnées passivement dans cet ACP.

The PCA of cuticle morphology parameters (Fig. 4b) placed SI with highest loading on the first axis (which explained ∼75% of the variance in the stomata data), whereas PL and PW were strongly correlated with each other and with the second PCA axis. Although the overall absence of distinct clusters of populations in the PCA biplot illustrates readily the high intra-population variability of stomatal morphology and frequency, some patterns seem to emerge. Populations having leaves with high SI values (Dikome [L. a. rumpiensis], Korup [L. a. letouzeyi], and Ebogo [L. a. gracilicaulis]) are mostly present towards the left part of the biplot. Also, whereas leaves from Ebodjé (characterised by low PL and PW values) are mostly located in the lower part of the biplot, those from Mt Kala are located in the upper part of the biplot. When added passively to the PCA with the stomata data, climate variables more strongly correlated with the first PCA axis are those related to the seasonality of temperature and precipitation, and isothermality.

4 Discussion

Regional environmental conditions only account for a relatively small part of the intra-taxon variability, as shown by higher variance of SI at the intra-site level than variance at higher and lower hierarchical levels (intra-individual, inter-taxon, and inter-site). Inter-site and intra-site variances of SI were smaller than intra-individual variance only for L. a. africana, of which the two populations studied are located at similar altitude and experience similar climatic conditions (Fig. 4a). This suggests that in most subspecies leaves within individuals were less different than leaves among individuals, and that increasing the number of replicate leaves per individual would have little impact on the precision of the variance-component estimates (Solomon, 2005). Inter-site variance was highest (although it did not exceed intra-site variance) for L. a. letouzeyi, whose two studied populations also experience similar climatic conditions (Fig. 4a) but are located at different altitudes, with leaves from the higher-altitude population having higher SI values (6.3 ± 1.9% vs 8.6 ± 2.9%).

Several herbarium studies have demonstrated SI-CO2 relationships in ferns, Gingkoaceae, conifers and angiosperms (Beerling et al., 1998; Kouwenberg et al., 2003; Wagner et al., 2007). However, stomatal-frequency datasets based on studies of herbarium specimens may invariably introduce significant variation into any SI-altitude or SI-CO2 relationship, owing to variation in the environmental conditions experienced by the sampled individuals, e.g., in habitats and locations of the individual trees sampled, or in light environments associated with different positions of leaves within the canopy (Kürschner, 1997; Smith et al., 2010; Woodward and Bazzaz, 1988).

Although the CO2 mixing ratio in air remains constant over altitudinal gradients, the CO2 partial pressure (pCO2, as measured in Pascals [Pa]) is lower at higher elevation because of the decrease in atmospheric pressure. Several C3 plants show stomatal frequency responses to altitude-controlled changes in pCO2 when CO2 mixing ratio is unaltered (Johnson et al., 2005; Kouwenberg et al., 2003; Woodward and Bazzaz, 1988). In this respect, our data show contrasting results: although we analyzed only two populations of L. a. letouzeyi, which cover a small altitudinal range below 1000 m asl (230 vs 350 m asl), the observed stomatal-frequency increase could be related to a pCO2 decrease. Similarly, the fact that the highest SI values amongst the nine Leonardoxa populations occurred in L. a. rumpiensis growing at Dikome (1148 m asl), the highest-altitude site, could also be related to differences in pCO2. By contrast, the four populations of L. a. gracilicaulis that span the greatest altitudinal gradient (∼920 m) amongst the four subspecies do not show an inverse relationship between stomatal frequency and altitude. In L. a. gracilicaulis, the highest SI values (8.3%) were observed in the Ebogo population (669 m asl) and the lowest SI values were observed at Iboti (1009 m asl). Overall, these results confirm the notion that stomatal-frequency differences can be negligible when comparing populations growing at altitudes below 1000 m asl (McElwain, 2004) or when the CO2 gradient is too shallow, and suggest that other factors overriding the effect of CO2 may have an influence on the SI values and thus on the plants’ ecophysiology.

Part of the intra-site variability is certainly related to variation among leaflet positions. In fact, our analysis of leaflets from pinnately compound leaves of L. a. gracilicaulis indicated statistically significant differences in SI among leaflet positions (Fig. 2). Although we did not attempt to estimate the intra-leaf (or the intra-leaflet) variance with multiple readings on each leaf (leaflet) lamina, as Poole et al. (2000) did in their extensive study of Alnus glutinosa leaves, our results are consistent with the results of those studies in showing that intra-leaf variability can add significant noise to stomatal frequency-environment relationships.

On the other hand, stomatal frequencies can also be influenced at intra-site and intra-individual levels by other factors, including the timing of leaf maturation, developmental stage of leaves (Marchi et al., 2008), length of shoots, and position in the canopy (Chen et al., 2001). Although the sampling strategy adopted to collect the Leonardoxa leaves did not take into account most of these factors, it is likely that the high variance at intra-site level reflects phenotypic plasticity within each of the L. africana subspecies, presumably in response to highly local variation in micro-environmental conditions. In fact, L. africana trees often grow in the shaded understorey, although they sometimes also grow in less shaded riverside settings (McKey, 2000), and are therefore exposed to varying degrees of light intensity and microclimate conditions, which make them liable to experience different light regimes. It will be essential to estimate the influence of these micro-environmental factors if stomatal responses to environmental changes are further studied in this genus.

Inter-site variance seems to be related to regional environmental (climatic) conditions, as suggested by the PCA (Fig. 4). In L. a. gracilicaulis, individuals growing in a strongly seasonal climate (Iboti, Fig. 4a) were characterized by the lowest stomatal frequencies, whereas populations growing in the least seasonal and most isothermal climate (Ebogo) had the highest stomatal frequencies (Table 1). The high SI values in L. a. rumpiensis, the highest observed in this study, could also be related to environmental conditions, as this subspecies occurs at high altitude and under low seasonality, in a climate characteristic of tropical cloud forests. Long leaf lifespan has been proposed as a factor that may reduce the selective pressure to regulate stomatal initiation (and thus stomatal density and index) (Haworth et al., 2011) based on conditions experienced during the leaf-formation season to optimise gas exchange in response to major inter-seasonal and inter-annual environmental changes such as variation in water availability. In comparison, plants with shorter leaf lifespans (such as deciduous species) may be better adapted to optimise for the conditions of the growing season and for year-to-year fluctuations of environmental variables (Miller-Rushing et al., 2009; Wagner et al., 1996). Because Leonardoxa trees possess leaves with a relatively long lifespan (>3 years) (McKey, 1984), our results may suggest that Leonardoxa tend to regulate their stomatal frequencies with respect to the least favourable climate conditions at the growing site (hydric stress, temperature). In addition to highlighting the extent of plastic responses, our study also shows that the two different subspecies growing at the same site (L. a. gracilicaulis and L. a. africana at Boundé), exhibit similar stomatal frequency values, although they are genetically well separated (Léotard et al., 2008). This may reflect a high degree of adaptation of the two subspecies to the environmental conditions at Boundé.

5 Conclusions

Our analysis of stomatal parameters at different hierarchical levels in four Leonardoxa subspecies focussed on the SI to investigate its variation among subspecies, among populations of the same subspecies, among individuals of the same population and among leaves of the same individual. Only in few cases did we detect clear relationships between stomatal frequency and altitude that could reflect the well-known negative SI–CO2 relationship. Although variance of SI at the intra-site level always exceeded the variance at the inter-site level, perhaps owing to highly variable micro-environmental conditions within sites, we could highlight the potential influence of environmental factors (e.g. seasonality of precipitation) on the SI at the inter-site level. The long life span of Leonardoxa leaves may be a key factor that causes these plants to optimise for year-to-year fluctuations of the environment with respect to the least favourable climate conditions at the growing site (hydric stress, temperature). What emerges, therefore, is that in modern Leonardoxa trees the stomatal frequency could reflect great phenotypic plasticity presumably in response to local variation in micro-environmental conditions (e.g., position in the canopy, exposure).

Acknowledgements

We thank the Ministry of Scientific Research and Innovation of the Republic of Cameroon for permission to carry out this study. We thank Xavier Garde and the IRD in Yaoundé for providing logistic help in the field, and the residents of the villages near our study sites for their hospitality. We also thank Geo Coppens (CEFE) for help with using the WorldClim database. This study was funded by two grants from the French Agence Nationale de la Recherche to DM: one from the “Biodiversity” programme (IFORA project) and one from the “Sixth Extinction” programme (C3A project). We are grateful to Giovanna Battipaglia for comments on an earlier version of this manuscript and to two reviewers, Margret Steinthorsdottir and Roberto Tognetti, for thoughtful and constructive comments.


Bibliographie

[Beerling, 1999] D.J. Beerling Stomatal density and index: theory and application (T.P. Jones; N.P. Rowe, eds.), Fossil Plants and Spores: Modern Techniques, Geological Society, London, 1999, pp. 251-256

[Beerling and Chaloner, 1993] D.J. Beerling; W.G. Chaloner Stomatal density responses of Egyptian Olea europaea L. leaves to CO2 change since 1327 BC, Ann. Bot., Volume 71 (1993), pp. 431-435

[Beerling et al., 1998] D.J. Beerling; J.C. McElwain; C.P. Osborne Stomatal responses of the “living fossil” Ginkgo biloba L. to changes in atmospheric CO2 concentrations, J. Exp. Bot., Volume 49 (1998), pp. 1603-1607

[Brouat and McKey, 2000] C. Brouat; D.B. McKey Origin of caulinary ant domatia and timing of their onset in plant ontogeny: evolution of a key trait in horizontally transmitted ant-plant symbioses, Biol. J. Linn. Soc., Volume 71 (2000), pp. 801-819

[Chen et al., 2001] L.Q. Chen; C.S. Li; W.G. Chaloner; D.J. Beerling; Q.G. Sun; M.E. Collinson; P.L. Mitchell Assessing the potential for the stomatal characters of extant and fossil Ginkgo leaves to signal atmospheric CO2 change, Am. J. Bot., Volume 88 (2001), pp. 1309-1315

[Chenuil and McKey, 1996] A. Chenuil; D.B. McKey Molecular phylogenetic study of a myrmecophyte symbiosis: did Leonardoxa/ant associations diversify via cospeciation?, Mol. Phylogenet. Evol., Volume 6 (1996), pp. 270-286

[Crawley, 2005] M.J. Crawley Statistics: An Introduction using R, Wiley, Chichester, UK, 2005

[Dos Santos, 2011] T. Dos Santos Variation dans la taille et la densité des stomates dans le complexe Leonardoxa africana (Fabaceae: Caesalpinioideae) : un pas vers l’étude écophysiologique d’un groupe d’« arbres à fourmis », Université Montpellier 2, France, 2011 (unpublished MSc thesis)

[Fiala et al., 1999] B. Fiala; A. Jakob; U. Maschwitz; K.E. Linsenmair Diversity, evolutionary specialization and geographic distribution of a mutualistic ant-plant complex: Macaranga and Crematogaster in South East Asia, Biol. J. Linn Soc., Volume 66 (1999), pp. 305-331

[Finsinger and Wagner-Cremer, 2009] W. Finsinger; F. Wagner-Cremer Stomatal-based inference models for reconstruction of atmospheric CO2 concentration: a method assessment using a calibration and validation approach, Holocene, Volume 19 (2009), pp. 757-764

[Gagen et al., 2011] M. Gagen; W. Finsinger; F. Wagner-Cremer; D. McCarroll; N. Loader; I. Robertson; R. Jalkanen; G. Young; A. Kirchhefer Evidence of changing intrinsic water-use efficiency under rising atmospheric CO2 concentrations in Boreal Fennoscandia from subfossil leaves and tree ring delta 13C ratios, Global Change Biol., Volume 17 (2011), pp. 1064-1072

[Gaume and McKey, 1998] L. Gaume; D.B. McKey Protection against herbivores of the myrmeciphyte Leonardoxa africana (Baill.)Aubrév. T3 by its principal ant inhabitant Aphomomyrmex afer Emery, C.R. Acad. Sci. Paris, Ser. III, Volume 321 (1998), pp. 593-601

[Gray et al., 2000] J.E. Gray; G.H. Holroyd; F.M. Van der Lee; A.R. Bahrami; P.C. Sijmons; I. Woodward; W. Schuch; A.M. Hetherington The HIC signalling pathway links CO2 perception to stomatal development, Nature, Volume 408 (2000), pp. 713-716

[Haworth et al., 2011] M. Haworth; C. Elliott-Kingston; J. McElwain The stomatal CO2 proxy does not saturate at high atmospheric CO2 concentrations: evidence from stomatal index responses of Araucariaceae conifers, Oecologia, Volume 167 (2011), pp. 11-19

[Hetherington and Woodward, 2003] A.M. Hetherington; F.I. Woodward The role of stomata in sensing and driving environmental change, Nature, Volume 424 (2003), pp. 901-908

[Hijmans et al., 2005] R.J. Hijmans; S.E. Cameron; J.L. Parra; P.G. Jones; A. Jarvis Very high resolution interpolated climate surfaces for global land areas, Int. J. Climatol., Volume 25 (2005), pp. 1965-1978

[Houston and McNamara, 1992] A.I. Houston; J.M. McNamara Phenotypic plasticity as a state-dependent life-history decision, Evolut. Ecol., Volume 6 (1992), pp. 243-253

[Johnson et al., 2005] D.M. Johnson; W.K. Smith; M.R. Silman Climate-independent paleoaltimetry using stomatal density in fossil leaves as a proxy for CO2 partial pressure. Comment and Reply, Geology, Volume 33 (2005), p. e82-e83 | DOI

[Jordan, 2011] G.J. Jordan A critical framework for the assessment of biological palaeoproxies: predicting past climate and levels of atmospheric CO2 from fossil leaves, New Phytol., Volume 192 (2011), pp. 29-44

[Kouwenberg et al., 2003] L.L.R. Kouwenberg; J.C. McElwain; W.M. Kürschner; F. Wagner; D.J. Beerling; F.E. Mayle; H. Visscher Stomatal frequency adjustment of four conifer species to historical changes in atmospheric CO2, Am. J. Bot., Volume 90 (2003), pp. 610-619

[Kürschner, 1997] W.M. Kürschner The anatomical diversity of recent and fossil leaves of the durmast oak (Quercus petraea Lieblein/Q. pseudocastanea Goeppert) – Implications for their use as biosensors of palaeoatmospheric CO2 levels, Rev. Palaeobot. Palynol., Volume 96 (1997), pp. 1-30

[Léotard et al., 2008] G. Léotard; A. Saltmarsh; F. Kjellberg; D. McKey Mutualism, hybrid inviability and speciation in a tropical ant–plant, J. Evol. Biol., Volume 21 (2008), pp. 1133-1143

[Leps and Šmilauer, 2003] J. Leps; P. Šmilauer Multivariate analysis of ecological data using CANOCO, Cambridge University Press, Cambridge, UK, 2003

[Marchi et al., 2008] S. Marchi; R. Tognetti; A. Minnocci; M. Borghi; L. Sebastiani Variation in mesophyll anatomy and photosynthetic capacity during leaf development in a deciduous mesophyte fruit tree (Prunus persica) and an evergreen sclerophyllous Mediterranean shrub (Olea europaea), Trees, Volume 22 (2008), pp. 559-571

[Marchi et al., 2004] S. Marchi; R. Tognetti; F.P. Vaccari; M. Lanini; M. Kaligarič; F. Miglietta; A. Raschi Physiological and morphological responses of grassland species to elevated atmospheric CO2 concentrations in FACE-systems and natural CO2 springs, Funct. Plant Biol., Volume 31 (2004), pp. 181-194

[McElwain, 2004] J. McElwain Climate-independent paleoaltimetry using stomatal density in fossil leaves as a proxy for CO2 partial pressure, Geology, Volume 32 (2004), pp. 1017-1020 | DOI

[McElwain and Chaloner, 1995] J.C. McElwain; W.G. Chaloner Stomatal density and index of fossil plants track atmospheric carbon dioxide in the Palaeozoic, Ann. Bot., Volume 76 (1995), pp. 389-395

[McElwain et al., 2002] J.C. McElwain; F.E. Mayle; D.J. Beerling Stomatal evidence for a decline in atmospheric CO2 concentration during the Younger Dryas stadial: A comparison with Antarctic ice core records, J. Quater. Sci., Volume 17 (2002), pp. 21-29

[McKey, 1984] D. McKey Interaction of the ant-plant Leonardoxa africana (Caesalpiniaceae) with its obligate inhabitants in a rainforest in Cameroon, Biotropica, Volume 16 (1984), pp. 81-99

[McKey, 1991] D.B. McKey Phylogenetic analysis of the evolution of a mutualism: Leonardoxa (Caesalpiniaceae) and its associated ants (C.R. Huxley; D.F. Cutler, eds.), Ant-plant interactions, Oxford Univ. press, Oxford, 1991, pp. 310-334

[McKey, 2000] D. McKey Leonardoxa africana (Leguminosae: Caesalpinioideae): a complex of mostly allopatric subspecies, Adamsonia, Volume 22 (2000), pp. 71-109

[Miller-Rushing et al., 2009] A.J. Miller-Rushing; R.B. Primack; P.H. Templer; S. Rathbone; S. Mukunda Long-term relationships among atmospheric CO2, stomata, and intrinsic water use efficiency in individual trees, Am. J. Bot., Volume 96 (2009), pp. 1779-1786

[Parmesan, 2006] C. Parmesan Ecological and evolutionary responses to recent climate change, Ann. Rev. Ecol. Evol. Syst., Volume 37 (2006), pp. 637-669 | DOI

[Pinheiro et al., 2013] J. Pinheiro; D. Bates; S. DebRoy; D. Sarkar; R Development Core Team nlme: Linear and Nonlinear Mixed Effects Models R package version 3. 1-109, 2013

[Poole et al., 2000] I. Poole; T. Lawson; J.D.B. Weyers; J.A. Raven Effect of elevated CO2 on the stomatal distribution and leaf physiology of Alnus glutinosa, New Phytol., Volume 145 (2000), pp. 511-521

[R Core Development Team, 2012] R Core Development Team R: A language and environment for statistical computing, reference index version 2.14, R Foundation for Statistical Computing, Vienna, Austria, 2012

[Raven, 2002] J.A. Raven Selection pressures on stomatal evolution. Tansley Review no. 131, New Phytol., Volume 153 (2002), pp. 371-386

[Roth-Nebelsick, 2007] A. Roth-Nebelsick Computer-based studies of diffusion through stomata of different architecture, Ann. Bot., Volume 100 (2007), pp. 23-32

[Royer, 2001] D.L. Royer Stomatal density and stomatal index as indicators of paleoatmospheric CO2 concentration, Rev. Palaeobot. Palynol., Volume 114 (2001), pp. 1-28

[Salisbury, 1927] E.J. Salisbury On the causes and ecological significance of stomatal frequency, with special reference to the woodland flora, Phil. Trans. Roy. Soc. B, Volume 216 (1927), pp. 1-65

[Sims et al., 2006] M. Sims; S. Wanless; M.P. Harris; P.I. Mitchell; D.A. Elston Evaluating the power of monitoring plot designs for detecting long-term trends in the numbers of common guillemots, J. Appl. Ecol., Volume 43 (2006), pp. 537-546

[Smith et al., 2010] R.Y. Smith; D.R. Greenwood; J.F. Basinger Estimating paleoatmospheric pCO2 during the Early Eocene Climatic Optimum from stomatal frequency of Ginkgo, Okanagan Highlands, British Columbia, Canada, Palaeogeogr. Palaeoclimatol. Palaeoecol., Volume 293 (2010), pp. 120-131

[Solomon, 2005] P.J. Solomon Variance Components, Encyclopedia of Biostatistics, John Wiley & Sons Ltd, Chichester, 2005

[Steinthorsdottir et al., 2013] M. Steinthorsdottir; B. Wohlfarth; M.E. Kylander; M. Blaauw; P.J. Reimer Stomatal proxy record of CO2 concentrations from the last termination suggests an important role for CO2 at climate change transitions, Quater. Sci. Rev., Volume 68 (2013), pp. 43-58

[Ter Braak and Šmilauer, 2003] C.J.F. Ter Braak; P. Šmilauer CANOCO reference manual and user's guide to CANOCO for Windows Version 4, Centre for Biometry Wageningen, Wageningen, The Netherlands, 2003

[Tognetti et al., 2000] R. Tognetti; A. Minnocci; J. Penuelas; A. Raschi; M.B. Jones Field water relations of three Mediterranean shrub species co-occurring at a natural CO2 vent, J. Exp. Bot., Volume 51 (2000), pp. 1135-1146

[Wagner et al., 1996] F. Wagner; R. Below; P. De Klerk; D.L. Dilcher; H. Joosten; W.M. Kürschner; H. Visscher A natural experiment on plant acclimation: Lifetime stomatal frequency response of an individual tree to annual atmospheric CO2 increase, Proc. Natl Acad. Sci. U. S. A., Volume 93 (1996), pp. 11705-11708

[Wagner et al., 2007] F. Wagner; H. Visscher; W. Kürschner; D.L. Dilcher Influence of ontogeny and atmospheric CO2 on stomata parameters of Osmunda regalis, Courier Forschungsinstitut Senckenberg, Volume 258 (2007), pp. 183-189

[Walther et al., 2002] G.R. Walther; E. Post; P. Convey; A. Menzel; C. Parmesan; T.J.C. Beebee; J.M. Fromentin; O. Hoegh-Guldberg; F. Bairlein Ecological responses to recent climate change, Nature, Volume 416 (2002), pp. 389-395

[Woodward, 1987] F.I. Woodward Stomatal numbers are sensitive to increases in CO2 from pre-industrial levels, Nature, Volume 327 (1987), pp. 617-618

[Woodward and Bazzaz, 1988] F.I. Woodward; F.A. Bazzaz The responses of stomatal density to CO2 partial pressure, J. Exp. Bot., Volume 39 (1988), pp. 1771-1781

[Young et al., 2006] J.J. Young; S. Mehta; M. Israelsson; J. Godoski; E. Grill; J.I. Schroeder CO2 signaling in guard cells: Calcium sensitivity response modulation, a Ca2+-independent phase, and CO2 insensitivity of the gca2 mutant, Proc. Natl Acad. Sci. U. S. A., Volume 103 (2006), pp. 7506-7511


Commentaires - Politique