Correction
10 Apr 2014: The PLOS ONE Staff (2014) Correction: Ecophysiological and Anatomical Mechanisms behind the Nurse Effect: Which Are More Important? A Multivariate Approach for Cactus Seedlings. PLOS ONE 9(4): e95405. https://doi.org/10.1371/journal.pone.0095405 View correction
Figures
Abstract
Background
Cacti establish mostly occurs under the canopy of nurse plants which provide a less stressful micro-environment, although mechanisms underlying this process are unknown. The impact of the combination of light and watering treatments on Opuntia streptacantha (Cactaceae) seedlings was examined.
Methods/Principal Findings
Ecophysiological [titratable acidity, osmotic potential (‘solute potential’, Ψs), relative growth rate (RGR) and their components (NAR, SLA, and LWR)], anatomical (chloroplast density, chloroplast frequency, and cell area), and environmental [photosynthetic photon flux density (PPFD) and air temperature] sets of variables were analyzed, assessing relationships between them and measuring the intensity of the relationships. Three harvests were carried out at days 15, 30, and 45. Ψs and acidity content were the most important responses for seedling establishment. The main anatomical and environmental variables were chloroplast density and water availability, respectively. Opuntia streptacantha seedlings establish better in the shade-watering treatment, due to higher Ψs and acidity, unaffected chloroplasts, and lower PPFD. In addition, the chloroplasts of cells under high-light and non-watering treatment were clumped closer to the center of the cytosol than those under shade-drought, to avoid photoinhibition and/or to better distribute or utilize the penetrating light in the green plant tissue.
Conclusions
Opuntia seedlings grow better under the shade, although they can tolerate drought in open spaces by increasing and moving chloroplasts and avoiding drastic decreases in their Ψs. This tolerance could have important implications for predicting the impact of climate change on natural desert regeneration, as well as for planning reforestation-afforestation practices, and rural land uses.
Citation: Delgado-Sánchez P, Yáñez-Espinosa L, Jiménez-Bremont JF, Chapa-Vargas L, Flores J (2013) Ecophysiological and Anatomical Mechanisms behind the Nurse Effect: Which Are More Important? A Multivariate Approach for Cactus Seedlings. PLoS ONE 8(11): e81513. https://doi.org/10.1371/journal.pone.0081513
Editor: Paul Adam, University of New South Wales, Australia
Received: July 15, 2013; Accepted: October 22, 2013; Published: November 28, 2013
Copyright: © 2013 Delgado-Sánchez et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This study was supported by Consejo Nacional de Ciencia y Tecnología-Fondos Sectoriales-Secretaría de Medio Ambiente y Recursos Naturales (2006-1-23818), Consejo Nacional de Ciencia y Tecnología-Fondo Ciencia Básica (156205 and 179982), Programa de Mejoramiento del Profesorado/103.5/08/5562 and Consejo Nacional de Ciencia y Tecnología-Fondos Sectoriales-Secretaría de Agricultura, Ganadería, Desarrollo Rural, Pesca y Alimentación (2004-C01-216). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Plant recruitment in arid environments often occurs only in years with above average rainfall or in safe sites under the canopy of nurse plants that provide shelter from high temperatures and low moisture. The nurse-protégé interaction is thought to be more frequent in harsh than in mild environments and appears to occur frequently in arid and semi-arid biomes [1]. The Cactaceae family includes the highest number of species protected by nurse plants [1]; for example, in the southern Chihuahuan Desert of Mexico, Opuntia streptacantha seedlings establish mostly under the canopy of nurse plants [2], [3]. However, physiological and anatomical mechanisms underlying the nurse effect on Opuntia seedlings are poorly known in arid environments [4] and the lack of this knowledge is impeding our understanding of these processes, e.g. there are no studies about the morphological [specific leaf area (SLA) and leaf weight ratio (LWR)] and physiological [net assimilation rate (NAR)] components of the relative growth rate (RGR), nor the anatomical responses such as chloroplast frequency, chloroplast density, and cell area.
Cactus plants generally exhibit CAM-photosynthesis, in which CO2 exchange occurs at night when the water vapor pressure difference between the air and the transpiring surfaces is lowest, resulting in high water-use efficiency in adult plants [5] but not necessarily during seedling stages and when water from rainfall is plentiful [6], [7]. The CAM pathway is very common in arid and semiarid environments, where plants are frequently exposed to low soil moisture and high solar radiation [6], [8]. There are numerous reports showing that incident photosynthetic photon flux density (PPFD) during daytime is positively correlated with nocturnal acid accumulation, in adult cacti [9], [10]. For cactus seedlings, Franco and Nobel (1989) [11] found that photosynthesis and growth of Ferocactus acanthodes under the canopy (diminished light) of nurse plants were lower than under direct sunlight. However, Gallardo-Vasquez and De la Barrera (2007) [12] found that nocturnal acid accumulation decreases for light - exposed cactus seedlings but increases for those shaded, regardless of their water treatment. Thus, the presence of CAM in shaded cactus seedlings can increase water-use efficiency and reduce the risk of photoinhibition [12]. It is possible that Opuntia seedlings established in open, highly illuminated spaces, show lower acid accumulation than those established in shade.
Diurnal movement of chloroplasts in the CAM plant Zygocactus truncates (epiphytic cactus) has also been reported under the combination of light and water stress [13], but not in Opuntia species or in cactus seedlings that require the protection of nurse plants. Chloroplasts are usually positioned on cell walls parallel to the leaf surface, but they can reorient along cell walls perpendicular to the leaf surface, parallel to the direction of incident light to avoid absorption [14]; [15].
Aridity is expected to increase in arid and semi-arid ecosystems due to global warming. Nobel (2010) [10] suggested that the future of succulents such as cacti is extremely bright under new climatic conditions expected in the 21st century. Thus, understanding how water availability and light intensity influence the ecophysiological and anatomical performance of cactus seedlings is crucial in having a better understanding of the organism's dynamics, adaptation, and plasticity [16], as well as of the impact of climate change on natural desert regeneration.
In this paper, the effect of the combination of light and watering treatments (shade and watering, shade and non-watering, light and watering, and light and non-watering) has been examined on O. streptacantha seedlings. We hypothesized that shade of nurse plants and watering favor seedling development. Opuntia streptacantha Lem. (Prickly Pear) is a perennial arborescent cactus of economic interest which has edible fruit and young pads. It is distributed throughout the semiarid lands of central Mexico [17], [18]. In order to corroborate our hypothesis, we assessed three types of response variables: ecophysiological [titratable acidity, osmotic potential (‘solute potential’, Ψs), and relative growth rate (RGR) and their components (NAR, SLA, and LWR)], anatomical (chloroplast frequency, chloroplast density, and cell area), and environmental (photosynthetic photon flux density (PPFD), air temperature, and soil water condition). In addition, we also analyzed the importance of ecophysiological and anatomical variables in order to know which is more important as mechanisms behind the nurse effect.
Materials and Methods
Study design
Opuntia streptacantha seeds were collected at Mexquitic de Carmona municipality of San Luis Potosi, Mexico (22°16′ N, 101°07′ W at 2,020 m asl), in the Southern Chihuahuan Desert, in 2007. The seed collection was carried out on a private land; the owner of the land gave permission to conduct the study on this site. This study did not involve endangered or protected species.
At May 15th, 2008, seeds were sown in plastic trays with Sunshine Mix # 3 substrate (Sun Gro Horticulture Canada Ltd), which were placed in a greenhouse and watered every other day. Four months later, on September 15th, 2008, the seedlings were transplanted into pots with soil from the O. streptacantha habitat and maintained at field capacity for 20 days, also in a greenhouse. All the experimental units were randomised within the glasshouse. Then, two treatments and two levels within each treatment were randomly assigned to 100 seedlings corresponding to a total of 25 replicates after treatment. Treatments included: neutral shade (25% light) with watering (SW) at field capacity every other day, neutral shade without watering (SNW), light (100% light) with watering (LW), and light without watering (LNW). The watering regime resembled natural conditions. The waterless treatments mimic dry periods with no rain. Treatments with water were intended to model years in which rains were frequent throughout the entire wet season. In addition, watering provided a good control of added water quantity to experimental units. Field capacity was determined in pots containing overwatered mixture and allowed to drain overnight. To reach field capacity, 200 mL water per pot was required. Treatments with high solar radiation simulated solar radiation in open spaces without vegetation, whereas shaded treatments represented the low solar radiation (≈75% of light reduction) found under shrubs in the area. All of these treatments allowed us to measure the ecophysiological and anatomical responses of seedlings to contrasting environmental conditions.
In order to simulate low solar radiation, plants were covered with a plastic film mesh which created shade; for high radiation treatments, plants were not covered. Seedlings in the low moisture treatments were not watered, while those in the high moisture treatments were watered every other day to field capacity during the entire course of the experiment.
The response variables were measured 15, 30, and 45 days after transplanting, and at these times three seedlings were harvested per treatment (light and watering) and the two levels within each treatment. At the same time, PPFD and temperature were measured using a portable fluorometer (pulse-amplitude modulated photosynthesis yield analyzer, Mini-PAM; Walz GmbH, Effeltrich, Germany). This instrument is equipped with a micro-quantum sensor to monitor PPFD and a thermocouple to measure temperature.
The following parameters were evaluated: osmotic potential (Ψs), titratable acidity content, chloroplast density, chloroplast frequency, and cell area. In addition, chloroplast arrangement was described.
Osmotic potential (Ψs)
Osmotic potential was evaluated as an estimation of plant water status, because the accumulation of malate in the vacuole of CAM species causes a reduction in sap osmotic potential, which may favour water absorption. Water status measurements were completed at predawn. The experimental unit consisted of one seedling per pot to obtain sufficient sap extract. Ex situ determinations were made of the Ψs of the sap obtained by mechanical compression [19]. Ψs was recorded in C-52 sampling chambers connected to a Wescor HR-33T dew-point microvoltmeter (Wescor Inc., Logan, Utah, U.S.A.). Sampling chambers were calibrated prior to use with standard NaCl solutions [19].
Titratable acidity content
Organic acid content, mainly malic acid, but also oxaloacetic acid and citric acid [6] in seedlings was estimated from titratable acidity, which was determined at 6:00 a.m., when its concentration is highest (6). Transverse samples of plant tissue were obtained from seedling shoots using a steel borer of cross section 0. 2 cm2. The plant material was sectioned and preserved in ethanol (80%) in 1.5 mL Eppendorf tubes. Titration was carried to neutrality with a 0.01 N NaOH solution to determine the organic acid concentration (mmol equiv. H+ m−2) per tissue according to Pearcy et al. (1991) [20].
Photosynthetic area
To determine the photosynthetic area of Opuntia streptacantha stems, seedlings were scanned with an HP 2400 ScanJet scanner. The images were analyzed with the software ImageJ 1.40 g (Wayne Rasband National Institutes of Health, USA).
Seedling growth
We measured relative growth rate (RGR) and its components at three time periods on harvested samples. The starting point for measuring changes in biomass was 15 days after transplanting. The calculations were realized using the formula Relative growth rate (RGR) = Net assimilation rate (NAR, g cm−2 day−1) * Specific leaf area (SLA, cm2 g−1) * leaf weight ratio (LWR, g g−l), following Shipley (2000) [21]. NAR represents an increase in plant total biomass (TB) per total leaf (or photosynthetic) area (TLA) unit and time (T) unit (NAR = (TB2−TB1)/(T2−T1). 2/(TLA1 + TLA2), mg/day/cm2); NAR is a physiological component because it is a measure of whole-plant daily net photosynthetic rate weighted by the rate of change in plant carbon content [22]. SLA is a morphological component which is determined by leaf dry matter concentration and leaf thickness [23], [24], [25]. In a succulent cactus without leaves, such as Opuntia, the SLA is the specific photosynthetic structure area (SPSA = photosynthetic structure area/photosynthetic structure dry mass at time t); LWR measures the allocation of biomass to leaves or photosynthetic organs vs. other plant parts (LWR = leaf (photosynthetic structure) dry mass at time t/plant dry mass at time t) [24], [25].
Light microscopy
Small samples of three O. streptacantha seedlings (five cells by seedling), including epidermis, hypodermis, and chlorenchyma were collected at 9:00 hours and immediately fixed with cold 4% glutaraldehyde in a phosphate buffer at pH 7.4, and dehydrated with graded ethanol series. Samples were hand-sectioned with a low profile microtome blade and the transverse sections were examined by light microscopy (Leica DM2000; Leica Microsystems Inc., Buffalo Grove, Il, USA) examination. Cell micrographs were obtained (LAS Imaging Software; Leica Microsystems Inc., Buffalo Grove, Il, USA) and digitized to determine the number of chloroplasts per cell and the cell cross-sectional area for every cell for five cells per micrograph (five micrographs per seedling per treatment). Chloroplasts were counted in one plane for each cell to estimate the number of chloroplasts on the top surface of each cell. Chloroplast density was calculated as the number of chloroplasts per 100 µm2 of cell area, and the chloroplast frequency as the number of occasions that chloroplasts are observed per grid unit within the cell, divided by the total grid number. Cell size was measured with the software ImageJ 1.40 g (Wayne Rasband National Institutes of Health, USA) using the microscope scale as reference.
Statistical analysis
A canonical discriminant analysis was used to understand the complex relationship between three groups of original variables (ecophysiological, anatomical, and environmental) and the relative contribution of these variables (within each group) to explain the effect of our experimental design, describing the linear combination of the original variable coefficients (canonical variables) that maximally discriminates among groups.
Data were analyzed d using a canonical discriminant analysis that allowed the identification of differences among the combined treatments (SW, SNW, LW and LNW) mentioned previously, using data obtained from variables from individuals belonging to each group and by utilizing linear functions from the quantified variables. Then the groups of individuals were separated to maximize the variance among treatments and to minimize the variance within them. The standardized canonical coefficients show the contribution of each of the joint variables which were analyzed to each of the canonical functions, and individually indicate the relative importance of each variable [26], [27]. In addition, a classificatory discriminant analysis was used to estimate the number and percent of entities classified correctly or incorrectly into each group of combined treatments. Both analyses were carried out using PROC CANDISC and PROC DISCRIM procedures.
The relationship among the three sets of variables (ecophysiological, anatomical, and environmental) was determined by performing a generalized canonical correlation analysis [28]. Only selected variables (P<0.005) by the canonical discriminant analysis were included for canonical correlation analyses. From the generalized canonical correlation results, a simple canonical correlation analysis was applied to the original variables related to their canonical variates (ecophysiological, anatomical, and environmental), arranged as two sets of plant and environmental variables. The goal of canonical correlation is to analyze the relationships between two sets of variables, and to elucidate the relationship between the two sets of variables. One set of variables consists of response variables and the other set consists of explanatory variables. A redundancy analysis was also performed to calculate the variance in a set of original variables explained by a canonical variate of another set [26].
In addition, three-way ANOVA were performed for osmotic potential, titratable acidity, specific photosynthetic structure area, leaf weight ratio, and chloroplast density, with light, watering and time as main factors, in order to determine the differences among treatments in the variables that contributed significantly to the canonical functions. Tukey tests were carried out to discern specific differences among treatments. All statistical analyses were performed with SAS software (SAS Institute Inc.).
Results
Photosynthetic photon flux density (PPFD)
The PPFD is a measurement unit that expresses the light quantum in photons of solar energy (specific to wavelength) related specifically to photosynthesis. The daily total PPFD received by O. streptacantha seedlings was 1059±5.5 µmol m−2 s−1 for the light treatments and 267±4 µmol m−2 s−1 for the shade treatments. Similar photo flux densities were observed in the three harvests.
Temperature
Temperatures for all treatments were of 39±0.55°C in the light treatments for the three harvests, and 33±0.38°C for the shade treatments.
Canonical discriminant analysis
We determined that two discriminant functions accounted for 99% of the data set total variation, contributing significantly to the separation among treatments (Wilks' λ: P<0.0001, n = 36). The first function (eigenvalue of 17.69; P<0.0001) explained 93% of the total variation, and the second (eigenvalue of 1.17; P<0.0001) explained 6% of the total variation. The remaining function accounted for 1% of total variation and was interpreted as trivial.
Variables that contributed significantly to the two canonical functions were present in the standardized canonical coefficients, and those that most contributed to centroid separation among treatments were Ψs, chloroplast density, specific photosynthetic structure area, titratable acidity, and leaf weight ratio in the first function, as well as cell area and chloroplast frequency in the second function (Figure 1).
The classificatory discriminant analysis showed that centroids of each group were significantly different (P<0.0001). All of the observations were correctly classified for all treatments.
Canonical correlation analysis
The generalized canonical correlation revealed that there was only one canonical variate. Vector defining canonical variates to produce the vector of the first canonical variates are shown in Table 1. The first eigenvalue (λ1 = 2.00) accounted for most of the variance of the first canonical variate, indicating that it is generated by only one common factor, and indicated that the three sets of variables are highly correlated and all had similar weight in determination of the first canonical variate. The second (λ2 = 0.68) and third eigenvalues (λ3 = 0.29) accounted for less of the variance. The first eigenvector elements (e1 = 0.63ecophysiological; 0.52anatomical; 0.56environmental) showed that the anatomical set had the lowest value. The results from the correlation matrix among canonical variates indicated that highest correlation between ecophysiological and environmental variates (0.65), followed by ecophysiological and anatomical (0.56), and the lowest between anatomical and environmental (0.19).
Table 2 shows correlations between the original variables and the canonical variates within sets. These correlation coefficients corroborated that Ψs and titratable acidity were highly correlated with ecophysiological (their first canonical) variate, which means that these variables were important in the resolution of the canonical variate. The highly associated variable with anatomical variables was chloroplast density. With respect to the environmental variables, air temperature in the harvests and PPFD were highly correlated. These were the variables selected for the canonical correlation analysis.
The first canonical correlation is 0.76 (84% of variance, P<0.004), however, the remaining correlation was not significant. Data for the first pair of canonical variates (cross-loadings) appear in Table 3. Percentage of variance (81%), redundancy (0.24), and canonical correlation (0.76) indicate that the first pair of canonical variates was highly correlated. The cross-loadings showed that Ψs was negatively correlated with PPFD, water treatment, titratable acidity, LWR, SPSA, and chloroplasts density. In addition, titratable acidity was correlated with chloroplast density, and both variables were correlated with PPFD and water treatment (Figure 1, Table 3).
ANOVA results
We determined the differences among treatments in the variables that contributed significantly to the canonical functions: osmotic potential, titratable acidity, specific photosynthetic structure area, leaf weight ratio, and chloroplast density.
The osmotic potential was more affected under drought and high-light intensity than under drought and low-light intensity at the three times measured (Figure 2). Tritatable acidity was lowest under drought and high-light intensity at days 15 and 30, but at day 45 the acidity content was lowest under both drought and high- and low-light intensity (Figure 2).
a) osmotic potential (MPa) and, b) titratable acidity (mmol equiv. H+ m−2). Bars represent the means ±SE (P<0.001; n = 3). Different letters indicated differences between treatments and days.
RGR and NAR were not important to the canonical functions, but they also were affected by drought, high solar radiation, and time (P<0.001; Figure 3). The lowest RGR was found under combined drought and both high- and low- solar radiation treatments at days 15, 30 and 45. This also provided the highest RGR under combined watering and low solar radiation treatments at days 15, 30 and 45, as well as under watering and high solar radiation at day 15. The highest NAR was found under watering and low-light intensity at days 30 and 45.
a) Relative growth rate (RGR; g g−l day−1); b) specific photosynthetic structure area (SPSA; cm2 g−1); c) leaf weight ratio (LWR; g g−l), and d) net assimilation rate (NAR; g cm−2 day−1). Bars represent the means ±SE (P<0.001; n = 3).
The highest specific photosynthetic structure area was found under combined watering and both high- and low- solar radiation treatments at days 15 and 30, as well as under low-light intensity and drought at day 45 (P<0.001; Figure 3). The highest leaf weight ratio was found under combined drought and both high- and low- solar radiation at days 30 and 45 (P<0.001; Figure 3).
Chloroplast density was higher under non-watering (both shade and light) than under watering treatments across measured times, but it was highest under combined light and non-watering treatment at day 30 (P<0.001; Figure 4).
a) chloroplast density (N° 100 µm−2); b) chloroplast frequency (%) and, c) cell area (µm2). Bars represent the means ±SE (P<0.001; n = 3).
Chloroplast frequency and cell area were not important to the canonical functions, the former was lightly affected by the treatments because it was lower under watering (both high- and low- light intensity) at day 15, as well as under all combined treatments at day 45 (P<0.001; Figure 4). The highest cell area was found at shade-watering treatment after 45 days of the experiment (P<0.001; Figure 4).
Chloroplasts arrangement
We found different chloroplast arrangements among the treatments in chlorenchyma cells. For cells from the watering treatment groups in both light and shade conditions, chloroplasts were always dispersed in the cytosol of the cells in the three harvests (Fig. 5a–f). The chloroplasts of cells from the LNW treatment group show a gradual aggregation toward the center of the cytosol for the three harvests, but this effect is most pronounced for the last harvest (Fig. 6b, d, f). A similar response of chloroplast in cells from the SNW group is observed, but chloroplasts are less aggregated than in cells of the LNW treatment group (Fig. 6a, c, e). In this sense, a similar behavior of chloroplast arrangement is observed when the seedlings were exposed to non-watering treatments, particularly at 45 days of the SNW treatment in comparison to LNW treatment at 30 days.
SW (a, c, e) and LW (b, d, f) at 15, 30 and 45 days, respectively. An arrow indicates the chloroplasts arrangement.
a) SNW at 15 days, chloroplasts clumped in one area of the cytosol (arrow); b) LNW at 15 days, chloroplasts clumped in the peripheral cytosol (arrow); c) SNW at 30 days, chloroplasts clumped in one area of the cytosol (arrow); d) LNW at 30 days, chloroplast clumped located close together at adjoining cells (arrow); e) SNW at 45 days, chloroplasts clumped in one area of the cytosol (arrow) and f) LNW at 45 days, chloroplasts densely clumped in the center of the cytosol (arrow).
Discussion
Ecophysiological and anatomical responses: Which is more important as mechanism behind the nurse effect?
In this study, ecophysiological responses are the most important for O. streptacantha seedling growth, followed by anatomical responses. In terms of seedling biology, ecophysiological variables such as Ψs and titratable acidity are very important for CAM plants, because acidity is regulated mainly by the water status of the plant [5], [8], [29], [30], [31], [32], [33], [34]. In this case, we found that Ψs of O. streptacantha seedlings is negatively correlated with titratable acidity. Interestingly, the Ψs was not particularly low (−1.5 MPa) for seedlings in the most stressful treatment (light and non-watering) at 45 days, which is similar to findings for adult Opuntia plants, i.e. during periods of high water stress, stem tissue water potentials are characteristically −1.2 to −1.6 MPa [32]. Our results show that Opuntia seedlings have the ability to osmotically adjust to tolerate drought, similar to findings for adult Opuntia plants, in which acidity content and Ψs decrease in cladodes under drought conditions [32], [35]. The importance of water status for the growth of O. streptacantha seedlings is also in agreement with studies on other cactus seedlings [12], [36].
Main ecophysiological mechanisms behind the nurse effect
Shade and higher soil humidity under nurse plants are two very important environmental factors facilitating seedling establishment [1]. Our results clearly show that for seedlings under shade the water status (in this case, the osmotic potential) and the acidity content did not decrease as rapidly as in high light conditions; this protection under nurse plants slows the effects of stress. Thus, Opuntia seedlings appear to have adapted to the protection of the shade of nurse plants, similar to findings for other cacti [12].
Growth frequently decreases under low water availability in a large number of species [37], [38], [39], including globose cactus seedlings [4]. However, RGR of O. streptacantha seedlings was not important inside the ecophysiological variables, probably because our study species is an obligatory CAM-photosynthesis [40] and other cactus seedlings could be CAM facultative [41], [42]. Similar findings were found under drought for Clusia hilariana seedlings, another obligatory CAM-photosynthesis plant [43]. Nevertheless, we found differences among treatments for RGR, having higher RGR in watering treatment. Interestingly, we also found higher cell size in the same treatment.
Main anatomical mechanisms behind the nurse effect
The chloroplasts in cells of many plants can be displaced and rearranged in accordance with the incidence of PPFD, with chloroplasts exhibiting maximum surface area at low PPFD and aggregated under high PPFDs [44]. Light-induced chloroplast movement is an adaptive response [45], because the aggregation of chloroplasts tends to avoid photoinhibition and/or to better distribute or utilize the penetrating light in the plant leaf tissue [46], [47], [48]. We suggest that Opuntia seedlings established in open, highly illuminated spaces show chloroplasts movement in order to tolerate light and/or water stress. In this study, we found that chloroplasts are highly clumped in the center of the cell under high solar radiation and drought, which is in agreement with findings that chloroplast aggregation occurs under combined light and water stress [13], [49]. This behavior is also observed under shade and drought treatment, although to a lesser degree, chloroplasts also aggregate. Indeed, chloroplast arrangement is similar for LNW treatment at day 30 and SNW treatment at day 45.
Under field conditions, O. streptacantha seedlings are primarily found in shade under nurse plants [2], [3]. Hence, we suggest that the chloroplasts of seedlings under nurse plants move to maximize photosynthesis, because light coming through the canopy provided by nurse plants is weak, diffuse, and transient owing to the movement of the sun, as proposed by Park et al. (1996) [50] for shade plants. Interestingly, seedlings under non-watering treatments, both in light and shade, showed higher chloroplast density, contrary to that described for other plants in which the number of chloroplasts decreases under water stress or when exposed to intense light [51], [52], [53]. This could be a mechanism used by O. streptacantha seedlings to maintain photosynthetic activity under stress conditions.
Under drought condition with lower Ψs, chloroplast clumping was also correlated with high chloroplast density and titratable acidity. Chloroplasts move towards a position of higher CO2 concentration for efficient photosynthesis [54], because higher concentration of malate in CAM plants is found in the vacuole. Thus, chloroplasts are closer to the vacuole for the transport of malate.
Conclusions
There are clearly complex interactions between drought and light that promote different ecophysiological and anatomical adaptations, which have direct consequences in the seedling growth under nurse plants or in open spaces. Ecophysiological variables are the most important for O. streptacantha seedlings, followed by anatomical and environmental variables. Ecophysiological responses are more sensitive to environmental conditions, and these would be mainly affecting anatomical responses, but at the same time, anatomical variables are affected in less degree by environmental conditions than ecophysiological variables. Water stress decreased the osmotic potential under drought and high-light intensity at the three measured times. Tritatable acidity was lowest under drought and high-light intensity at days 15 and 30, but at day 45 the acidity content was lowest under both drought and high- and low-light intensity. The lowest RGR was found under combined drought and both high- and low- solar radiation treatments at days 15, 30 and 45. The highest NAR was found under watering and low-light intensity at days 30 and 45. The highest specific photosynthetic structure area was found under combined watering and both high- and low- solar radiation treatments at days 15 and 30, as well as under low-light intensity and drought at day 45. The highest leaf weight ratio was found under combined drought and both hig- and low- solar radiation at days 30 and 45. Chloroplast density was higher under non-watering (both shade and light) than under watering treatments at days 15, 30 and 45, but it was highest under combined light and non-watering treatment at day 30. Chloroplast frequency was lower under watering (both high- and low- light intensity) at day 15, as well as under all combined treatments at day 45. The highest cell area was found at shade-watering treatment after 45 days.
We found that O. streptacantha seedlings grow well in exposed sites with watering, but these conditions are not so common in arid and semi-arid environments. Thus, these plants can tolerate high-light and drought in open spaces by increasing and moving chloroplasts and also avoiding drastic decrease in their osmotic potential, although they could suffer damage in longer water stress periods. Higher soil humidity under nurse plants could be longer during the year than in open sites, for this reason, O. streptacantha seedlings establish better in conditions simulating microhabitat under the shade of nurse plants with high soil humidity in the rainfall period, having higher Ψs and titratable acidity, as well as dispersed (unaffected) chloroplasts. These results could have important implications for predicting the impact of climate change on natural desert regeneration, as well as for planning reforestation-afforestation practices, and rural land uses.
Acknowledgments
The seed collection was carried out on a private land; the owner of the land gave permission to conduct the study on this site. This study did not involve endangered or protected species. We thank Teresa Terrazas and Enrique Jurado for the suggestions to improve the manuscript, as well as Paul Riley for grammatical review.
Author Contributions
Conceived and designed the experiments: PDS LYE JF. Performed the experiments: PDS. Analyzed the data: PDS LYE LCV. Contributed reagents/materials/analysis tools: PDS LYE JFJB LCV JF. Wrote the paper: PDS JFJB JF.
References
- 1. Flores J, Jurado E (2003) Are nurse-protégé interactions more common among plants from arid environments? J Veg Sci 14: 911–916.
- 2. Yeaton RI, Romero-Manzanares A (1986) Organization of vegetation mosaics in the Acacia schaffneri - Opuntia streptacantha association, southern Chihuahuan Desert, Mexico. J Ecol 74: 211–217.
- 3. Vargas-Mendoza MC, González-Espinoza M (1992) Habitat heterogeneity and seed dispersal of Opuntia streptacantha (Cactaceae) in nopaleras of Central Mexico. Southw Nat 37: 379–385.
- 4. Martínez-Berdeja A, Valverde T (2008) Growth response of three globose cacti to radiation and soil moisture: An experimental test of the mechanism behind the nurse effect. J Arid Environ 72: 1766–1774.
- 5. Winter K, Aranda J, Holtum JAM (2005) Carbon isotope composition and water-use efficiency in plants with crassulacean acid metabolism. Funct Plant Biol 32: 381–388.
- 6.
Nobel PS (2010) Desert wisdom/agaves and cacti: CO2, water, Climate change. New York: iUniverse, Inc.
- 7. Winter K, García M, Holtum JAM (2011) Drought-stress-induced up-regulation of CAM in seedlings of a tropical cactus, Opuntia elatior, operating predominantly in the C3 mode. J Exp Bot 62: 4037–4042.
- 8. Cushman JC (2001) Crassulacean acid metabolism: a plastic photosynthetic adaptation to arid environments. Plant Physiol 127: 1439–1448.
- 9. Nobel P (1983) Spine influences on PAR interception, stem temperature, and nocturnal acid accumulation by cacti. Plant Cell Environ 6: 153–159.
- 10. Sortibrán L, Tinoco-Ojanguren C, Terrazas T, Valiente-Banuet A (2005) Does cladode inclination restrict microhabitat distribution for Opuntia puberula (Cactaceae). Am J Bot 92: 700–708.
- 11. Franco AC, Nobel PS (1989) Effect of nurse plants on the microhabitat and growth of cacti. J Ecol 77: 870–886.
- 12. Gallardo-Vásquez JC, De la Barrera E (2007) Environmental and ontogenetic influences on growth, photosynthesis, and survival for young Pitayo (Stenocereus queretaroensis) seedlings. J Prof Assoc Cactus 9: 118–135.
- 13. Kondo A, Kaikawa J, Funaguma T, Ueno O (2004) Clumping and dispersal of chloroplasts in succulent plants. Planta 219: 500–506.
- 14.
Evert RF, Eichhorn SE (2006) Esau's plant anatomy: Meristems, cells, and tissues of the plant body. Their structure, function, and development. New Jersey: Wiley.
- 15. Li Z, Ahn TK, Avenson TJ, Ballottari M, Cruz JA, et al. (2009) Lutein accumulation in the absence of zeaxanthin restores nonphotochemical quenching in the Arabidopsis thaliana npq1 mutant. Plant Cell 21: 1798–1812.
- 16. Rosas U, Zhou RW, Castillo G, Collazo-Ortega M (2012) Developmental reaction norms for water stressed seedlings of succulent cacti. Plos One 7(3): e33936.
- 17. Littlejohn Jr RO, Williams GJ III (1983) Diurnal and seasonal variations in activity of crassulacean acid metabolism and plant water status in a northern latitude population of Opuntia erinacea. Oecologia 59: 83–87.
- 18. Winter K, Lesch M (1992) Diurnal changes in chlorophyll a fluorescence and carotenoid composition in Opuntia ficus-indica, a CAM plant, and in three C3 species in Portugal during summer. Oecologia 91: 505–510.
- 19.
Rundel PW, Jarrel WM (1991) Water in the environment. In: Pearcy RW, Ehleringer J, Mooney HA, Rundel PW, Plant physiological ecology. Field methods and instrumentation. London: Chapman and Hall, pp. 29–56.
- 20.
Pearcy RW, Ehleringer J, Mooney HA, Rundel PW (1991) Plant physiological ecology. Field methods and instrumentation. Chapman and Hall, London.
- 21. Shipley B (2000) Plasticity in relative growth rate and its components following a change in irradiance. Plant Cell Environ 23: 1207–1216.
- 22. McKenna MF, Shipley B (1999) Interacting determinants of interspecific relative growth: empirical patterns and a theoretical explanation. Ecoscience 6: 286–296.
- 23. Witkowski ETF, Lamont BB (1991) Leaf specific mass confounds leaf density and thickness. Oecologia 88: 486–493.
- 24. Shipley B (1995) Structured interspecific determinants of specific leaf area in 34 species of herbaceous angiosperms. Funct Ecol 9: 312–319.
- 25. Shipley B (2002) Trade-offs between net assimilation rate and specific leaf area in determining relative growth rate: relationship with daily irradiance. Funct Ecol 16: 682–689.
- 26.
McGarigal K, Cushman S, Stafford S (2000) Multivariate statistics for wildlife and ecology research. New York: Springer-Verlag.
- 27. Yáñez-Espinosa L, Terrazas T, López-Mata L (2006) Integrated analysis of tropical trees growth: A multivariate approach. Ann Bot 98: 637–645.
- 28. Ortega-Alcalá J, González-Cossío F, Castillo-Morales A, Martínez-Garza A (1992) Analysis of canonical correlation for M sets of variables. Agrociencia (México), Serie Matemáticas Aplicadas 3: 201–227.
- 29. North GB, Nobel PS (1992) Drought-induced changes in hydraulic conductivity and structure in roots of Ferocactus acanthodes and Opuntia ficus-indica. New Phytol 120: 9–19.
- 30. Pimienta-Barrios E, González del Castillo-Aranda ME, Nobel PS (2002) Ecophysiology of a wild Platyopuntia exposed to prolonged drought. Environ Exp Bot 47: 77–86.
- 31. Zañudo-Hernández J, González del Castillo-Aranda E, Ramírez-Hernández BC, Pimienta-Barrios E, Castillo-Cruz I, et al. (2010) Ecophysiological responses of Opuntia to water stress under various semi-arid environments. J Prof Assoc Cactus 12: 20–36.
- 32. Szarek SR, Johnson HB, Ting IP (1973) Drought adaptation in Opuntia basilaris. Significance of recycling carbon through crassulacean acid metabolism. Plant Physiol 52: 539–541.
- 33. Hanscom III Z, Ting IP (1978) Responses of succulents to plant water stress. Plant Physiol 61: 327–330.
- 34. Cushman JC, Borland AM (2002) Induction of crassulacean acid metabolism by water limitation. Plant Cell Environ 25: 295–310.
- 35. Goldstein G, Ortega JKE, Nerd A, Nobel PS (1991) Diel patterns of water potential components for the crassulacean acid metabolism plant Opuntia ficus-indica when well-watered or droughted. Plant Physiol 95: 274–280.
- 36. Franco AC, Nobel PS (1988) Interactions between seedlings of Agave deserti and the nurse plant Hilaria rigida. Ecology 69: 1731–1740.
- 37. Khurana E, Singh JS (2004) Germination and seedling growth of five-three species from tropical dry forest in relation to water stress: impact of seed size. J Trop Ecol 20: 385–396.
- 38. Galmés J, Cifre J, Medrano H, Flexas J (2005) Modulation of relative growth rate and its components by water stress in Mediterranean species with different growth forms. Oecologia 145: 21–31.
- 39. Bunker DE, Carson WP (2005) Drought stress and tropical forest woody seedlings: effect on community structure and composition. J Ecol 93: 794–806.
- 40. Nobel PS, Hartsock TL (1987) Drought-induced shifts in daily CO2 uptake patterns for leafy cacti. Physiol Plant 70: 114–118.
- 41. Altesor A, Ezcurra E, Silva C (1992) Changes in the photosynthetic metabolism during the early ontogeny of four cactus species. Acta Oecol 13: 777–785.
- 42. Loza-Cornejo S, Terrazas T, López-Mata L, Trejo C (2003) Características morfo-anatómicas y metabolismo fotosintético en plántulas de Stenocereus queretaroensis (Cactaceae): Su significado adaptativo. Interciencia 28: 83–89.
- 43. Cavalcante A, Braz MIG, de Mattos EA (2010) Germination biology and seedling growth of Clusia hilariana Schltdl., a dominant CAM-tree of drought-prone sandy coastal plains. Ecol Res 25: 781–787.
- 44. Nauš J, Rolencová M, Hlaváčková V (2008) Is chloroplast movement in tobacco plants influenced systemically after local illumination or burning stress? J Integr Plant Biol 50: 1292–1299.
- 45. Kagawa T, Sakai T, Suetsugu N, Oikawa K, Ishiguro S, et al. (2001) Arabidopsis NPL1: a phototropin homolog controlling the chloroplast high-light avoidance response. Science 291: 2138–2141.
- 46. Königer M, Delamaide JA, Marlow ED, Harris GC (2008) Arabidopsis thaliana leaves with altered chloroplast numbers and chloroplast movement exhibit impaired adjustments to both low and high light. J Exp Bot 59: 2285–2297.
- 47. Sharon Y, Beer S (2008) Diurnal movements of chloroplasts in Halophila stipulacea and their effect on PAM fluorometric measurements of photosynthetic rates. Aq Bot 88: 273–276.
- 48. Yamada M, Kawasaki M, Sugiyama T, Miyake H, Taniguchi M (2009) Differential positioning of C4 mesophyll and bundle sheath chloroplasts: Aggregative movement of C4 mesophyll chloroplasts in response to environmental stresses. Plant Cell Physiol 50: 1736–1749.
- 49. Königer M, Jessen B, Yang R, Sittler D, Harris GC (2010) Light, genotype, and abscisic acid affect chloroplast positioning in guard cells of Arabidopsis thaliana leaves in distinct ways. Photosynth Res 105: 213–227.
- 50. Park YII, Chow ES, Anderson JM (1996) Chloroplast movement in the shade plant Tradescanfia albiflora helps protect Photosystem II against light stress. Plant Physiol 111: 867–875.
- 51. Maxwell K, Marrison JL, Leech RM, Griffiths H, Horton P (1999) Chloroplast acclimation in leaves of Guzmania monostachia in response to high light. Plant Physiol 121: 89–95.
- 52. Rizhsky L, Liang H, Mittler R (2003) The water-water cycle is essential for chloroplast protection in the absence of stress. J Biol Chem 278: 38921–38925.
- 53. Nagy-Déri H, Péli ER, Georgieva K, Tuba Z (2011) Changes in chloroplast morphology of different parenchyma cells in leaves of Haberlea rhodopensis Friv. during desiccation and following rehydration. Photosynthetica 49: 119–126.
- 54.
Sato Y, Kadota A (2007) Chloroplast movements in response to environmental signals. In: Wise RR, Hoober JK, The structure and function of plastids. Dordrecht, The Netherlands: Springer. 527–537.