INTRODUCTION
Elevational and environmental gradients are powerful tools for investigating the importance of habitat filtering and community assembly processes. Indeed, the potential value of environmental gradients for demonstrating how different factors drive patterns of plant diversity and vegetation structure has long been recognized [
1,
2] and previous studies of environmental gradients have provided evidence that many plant species have different distributions based on their tolerances to climate, topography, and soil characteristics [
3,
4,
5,
6].
Studies of elevational gradients can also provide important insights into how tropical species, communities, and ecosystems are responding to ongoing anthropogenic climate change [
7,
8,
9,
10,
11,
12]. Elevational gradients provide particular insight into the influence of temperature because air and soil temperatures typically decrease linearly with elevation [
12,
13]. In the Amazon-Andes region, elevational gradients have been used to assess the effects of temperature on ecosystem productivity and community structure [
8,
14,
15,
16,
17,
18] and to predict how species and communities will be affected by climate change [
19,
20]. However, much less is known about how changes in rainfall and water availability across elevation and/or through time affect tropical plant communities in these hyper-diverse areas. Some previous studies have shown that changes in precipitation are driving directional changes in the community composition of tropical forests [
21,
22], but other studies have shown little or no relationship between changes in species composition and changes precipitation [
11], and it has even been suggested that observed changes in composition in relation to drought tolerances may be driven primarily by changes in temperature and the inherent relationships between species’ heat and drought tolerances (especially when measured based on geographic occurrence locations) [
23]. Given that climate change is altering tropical precipitation regimes [
24,
25], this incomplete knowledge about how important water availability is for determining species’ current and future distributions greatly limits our ability to predict the fate of these tropical communities and ecosystems.
Since plant water availability is a complex phenomenon involving the interplay between multiple climate variables, soil and environmental factors, and intrinsic species traits related to water uptake and use, the characterization of water availability is nearly impossible for diverse species assemblages spread across large environmental gradients. Water availability is generally believed to decrease at high elevations in the tropics (due to reduced rainfall, increased rainfall seasonality, and thinner soil organic material) [
26,
27]; however, changes in water availability across elevations are not necessarily monotonic or universal. Lowland rainforests can have heavy rains and high temperatures, but these climate variables are modulated by the soil water capacity [
28]. In the sub-montane forests (above lowland rainforests) there is generally higher precipitation but VPD remains high, with increased seasonality at mid-elevations [
13]. Cloud forests, which occur above the sub-montane forests, can have high water availability due to high relative humidity and regular cloud and fog immersion, (although there can be seasonal periods with reduced fog, leading to high irradiance and increased vapor pressure deficit, VPD)[
13]. As a consequence of regular cloud immersion and high water availability, cloud forests generally have thicker soil organic layers, abundant mosses, and high densities of epiphytes [
29,
30,
31]. Above the cloud forests, the habitat transitions to high-elevation grasslands (in the southern Peruvian Andes these grasslands are referred to as
puna) where climate is drier and colder, precipitation is more seasonal, and soils tend to have thinner soil organic layer [
32]. Embedded in these high-elevation grasslands, there can be stands of shrubs and small patches of forests where the local climate is more like that of the cloud forests [
33]. These changes in water availability across elevations should lead to strong habitat filtering of plant species based on their specific habitat requirements and differences in drought tolerances.
As water availability decreases, plant hydraulic conductance may decrease through different biophysical and physiological mechanisms to prevent embolism or desiccation [
34]. In this study, we used turgor loss point, cuticular conductance and solute leakage as tools to characterize relative leaf drought tolerances of individual understory plants and understory plant communities. Turgor loss point (TLP), the leaf water potential at which cell turgor is lost and leaf wilting occurs, provides physiological information pertaining to cell wall integrity, stomatal closure and, more generally, the extent to which plants can maintain metabolism as water availability decreases due to soil drying [
35]. As such, plant drought tolerances should be negatively correlated with TLP (i.e., plants that can withstand greater negative water potentials before wilting are considered to be more drought tolerant than plants that wilt at less negative water potentials).
Another mechanism that can promote drought tolerance is to avoid desiccation and maintain high relative water content, thereby resisting the cell damage that could be caused by dehydration. Under dry conditions, plants close their stomata, relying on their cuticle to prevent further water loss. Cuticular conductance (G
min) is the conductance to vapor diffusion across the leaf epidermis (i.e., through the cuticle and any leaky stomata) once stomata are closed [
36]. Plants with low G
min are considered to be more drought tolerant than plants with high G
min [
37].
Responses at the cellular level may also prevent damage to plants during drought [
38,
39,
40]. The measurement of solute leakage (SL) from plant tissues is a proven method to detect cytorrhysis (damage of cell walls) under laboratory conditions causing desiccation. As such, measurements of SL from leaves exposed to standardized desiccating conditions is a potential metric of relative drought tolerances. Indeed, SL has previously been shown to be an accurate indicator of drought tolerance of crop plants [
41,
42,
43,
44]. However, SL has rarely been used in ecological field studies. This is surprising given that SL has many advantages that make it particularly well-suited for studying patterns of plant drought tolerances in biodiverse communities and in remote settings [
45,
46]. Specifically, the protocol is easy to follow and standardize, large numbers of treatments can be performed relatively quickly, the equipment needed for the experiments are low-cost relative to other physiological measures of drought tolerance, and measurements can be made in areas with little or no existing infrastructure (e.g., measuring SL does not require electricity or compressed gases, and only uses polyethylene glycol, which is widely available as over-the-counter laxatives).
In this study, we measured leaf TLP, G
min and SL in hundreds of individual understory saplings and shrubs and assessed community-level metrics of leaf drought tolerance in approximately 60 plots located along a 3000+ meter elevation and 5000 mm yr
-1 rainfall gradient from lowland rainforests to highland gallery forests and high elevation shrubs. We tested the relationships of these community drought tolerance metrics with elevation (a proxy for temperature and VPD; VPD decreases at higher elevations due to lower air temperatures) and mean annual precipitation. Specifically, we addressed the two following questions and corresponding
a priori hypotheses: A) Are the three metrics of drought tolerance coordinated or independent? We hypothesized that TLP, G
min, and SL would all be correlated due to a coupling of leaf-level mechanisms to resist and tolerate drought [
47]. And B) How is drought tolerance related to rainfall and elevation across the gradient? We hypothesized that since adaptation and/or acclimation to drought can incur significant costs, drought-tolerant species would be relatively infrequent in sites with high rainfall and that plant communities in the lower hotter elevations would have greater average drought tolerances than plants at higher elevations. Increasing our understanding of the role that water availability and drought tolerance plays in structuring understory plant communities will improve our ability to predict the responses of tropical ecosystems to environmental changes.
METHODS
Study site. We conducted our study of understory plants and plant communities along an elevational gradient in Manu National Park in Southeastern Peru (Cusco and Madre de Dios regions;
Figure 1). The study area ranged from 400 to 3600 m asl and included forests along the Trocha Union ridge in the Kosnipata Valley (900 to 3600 m asl) and around the Cocha Cashu Biological Station (400 m asl;
Table 1). Most of the study region is characterized as being “super humid"; the relative orientation of prevailing winds and topography create a zone of frequent cloud immersion when cold Andean winds collide with warm, moist air from the Amazonian lowlands [
48]. Mean annual temperature (MAT) ranges from ~23 ºC at 400 m asl to <9 ºC at 3400 m asl (
Table 1), with a lapse rate for the understory of 0.53 ºC per 100 m elevation [
13]. Mean annual precipitation (MAP) has been measured at a series of meteorological stations along the transect and ranges from ~2000 mm yr
−1 at low elevations to 5000 mm yr
−1 at the mid-elevations and <1000 mm yr
−1 at the very highest elevations. The MAP at each plot location was estimated based on the observed relationship of precipitation vs. elevation in nearby weather stations (adjusted R
2 = 0.99) and interpolating to the plot locations based on their elevations as measured with handheld GPS (
Table 1). Mean vapor pressure deficit (VPD) in the understory increases monotonically with elevation along the study gradient from 0.01 in the lowlands to 0.12 kPa in the highlands [
13].
Vegetation Sampling. We established 58 understory plots of 5 x 5 m (25m
2) each. Approximately four plots were installed every 250 m in elevation between 900 and 3600 m asl. Four additional plots were installed at approximately 400 m asl in the lowland rainforest near the Cocha Cashu Biological Station (
Figure 1), but these lowland plots were only used to measure SL (see below).
We defined the understory to include all woody shrubs and saplings with basal diameters ≥1 cm at ground level and heights between 0.5 and 5 m (plant height was measured as the vertical distance from ground to the top of the plant canopy; in some cases, stem length was >5 m if the plant was leaning). For every individual, we measured basal diameter, total height, and diameter at breast height (dbh) when height was >1.3m. Most individuals (90%) were identified to family level in the field, and some to genus level. Unidentified individuals were assigned a temporary name in the field. Subsequently, the unidentified individuals were compared with specimens in regional botanical collections and determined to species level or vouchered morphospecies. Vouchers were compared across all plots and taxonomy was standardized. Since all analyses were conducted at the plot-level using the averages of individual-level measurements, the high frequency of morphospecies does not affect the results.
One terminal branch was collected from each individual understory plant in each plot and placed in black plastic bags along with wet paper towels inside and stem ends submerged in plastic bags with water to avoid desiccation during transport to temporary field laboratories established near the top and bottom of the main elevation gradient at 3500 m asl and at 1300m asl, respectively, and at Cocha Cashu Biological Station at 400 m asl. Once in the “laboratory”, branches were placed in buckets, recut under water, and allowed to rehydrate covered with plastic. All traits were measured on leaves from the same branch. Leaves were only included if they looked healthy and were not discolored; leaves that were noticeably affected by pathogens were excluded.
Turgor loss point (TLP). Before measurements of leaf water potential, branches were allowed to rehydrate for 10 hours. After rehydration, three mature leaves were collected from each branch and allowed to slowly desiccate on a table over a period of 6-8 hours between 02:00 and 10:00am local time. During this desiccation period, leaf mass (g) and water potential (
Ψ leaf, MPa) were measured periodically. Leaf mass was measured with an analytical balance (±0.002 g, Ohaus, Scout SPX223) and
Ψ leaf was measured using a Scholander-type pressure chamber (model 1000, PMS Instruments, Corvallis, Oregon, USA). Plotting the inverse of leaf water potential against relative water content allows determination of the turgor loss point (TLP) as the point of transition between linear and nonlinear portions of the pressure-volume (PV) curve. Lower TLP (i.e., more negative water potentials) indicates higher resistance to leaf wilting [
49]. We were not able to construct PV curves or determine TLP for all individuals or species because of short petioles, sap presence in the leaves, or logistical difficulties.
Cuticular conductance (Gmin). For measuring the leaf G
min, we harvested one to three fully expanded leaves from each individual’s sample branch after ≥4 hours rehydration. We cut the leaf petioles and then sealed the petiole ends with paraffin. Leaves were dried at night using LED lanterns (PPFD 1.73 – 4.6 µmol m
-2s
-1) on a table for 1 hour in order to induce stomatal closure. The leaves were then massed six times at intervals of 30 minutes. Cuticular conductance was calculated as the rate of water loss (g) over time [
36]. At each weighing, we measured the air temperature (°C) and relative humidity (%) with a handheld sensor (Kestrel 2500NV, Nielsen-Kellerman USA) to calculate the saturation vapor pressure (VP
sat, kPa). Following standard protocols, the value of G
min (mmol m
-2 s
-1) was calculated as the water loss rate divided by VPD (mol mol
-1) and twice the total leaf area as measured on scanned leaf images (cm
2) [
50]. Lower G
min indicates high drought tolerance since lower cuticular conductance enables maintenance of hydration based on stored water.
Solute leakage (SL). For measuring solute leakage (SL), we harvested three to five leaves from the sample branches of each individual after ≥ 4 hours rehydration. We lightly cleaned each leaf to remove surface debris and cut three pieces from the apical, middle, and bottom area of the leaf with a 1.27 cm diameter hole punch. We submerged the leaf pieces in 20 mL of 50% Polyethylene glycol 3350 (PEG), a hypertonic solution with an osmotic potential of approximately -8.4 MPa, for 12 hours [
51]. After the PEG treatment, the leaf pieces were rinsed three times with distilled water (2 µS) to remove any PEG residue. We then transferred the treated leaf pieces to 50 mL falcon tubes with distilled water and measured the initial electrolyte conductivity (C
0) using an Oakton CON 6+ Meter conductivity meter. After 12 hours in the distilled water, we took a third conductivity measurement (C
1). The leaf pieces and water were then boiled in the falcon tubes for 20 minutes to rupture any remaining intact cells. Twelve hours after boiling, we took a second conductivity measurement (C
2). We calculated the relative solute leakage (%) as: (C
1 – C
0)/(C
2 – C
0) *100 [
45,
46]. Higher SL values indicate greater cytorrhysis when exposed to desiccating conditions and thus lower relative leaf drought tolerance. This method assumes that solute leakage is the result of membrane breakage due to water stress within cells, that leaked ions are representative of initial endogenous concentrations in leaf samples before PEG treatment [
52], and that differences in leaf nutrition do not influence the concentration of leaked ions [
51]. In addition, it assumes that variation in solute leakage is not caused by variation in the capacity of water transport by aquaporins and ion transport by channel and carrier proteins embedded in the plasma membrane. Given that SL is related to several key hydraulic and biochemical parameters, including stomatal resistance and tissue osmotic potential [
53], water use efficiency [
54], and tissue nitrogen levels [
55], we find this a valuable and accessible field tool for assessing relative drought tolerance across many individuals and species.
Statistical Analysis. We used three leaves per individual to calculate TLP and G
min and five leaves per individual to obtain SL. The values for all individuals in each plot were then averaged to obtain community-weighted mean (CWM) trait values. We measured SL on each individual plant in each plot. In the case of TLP and G
min, we were not able to sample every individual, and we therefore calculated the CWM trait values as the mean of species-level average trait values weighted by the number of stems per species. Species-level averages were always based on measurements of individuals within each focal plot; in other words, drought tolerance measurements were never generalized between plots. Using plot-specific measurements accounts for any intraspecific variation due to local adaptation or acclimation and eliminates the need for taxonomic standardization across plots. We conducted correlation analysis between the traits to determine possible relationships between the different metrics of physiological drought tolerance. Multiple regression and correlation analyses were used to assess the relationships between the plots’ CWM of SL, TLP, and G
min with the elevation and MAP. All analyses were conducted in R [
56].
RESULTS
Floristic composition. The 58 sample plots contained a total of 834 understory plants with 363 (morpho)species, representing 89 genera and 50 families (
Table 1). The most abundant families in the plots were
Melastomataceae (170 individuals),
Rubiaceae (155 individuals) and
Chloranthaceae (55 individuals). The most abundant genera were
Miconia (165 individuals),
Hedyosmum (55 individuals) and
Palicourea (45 individuals). Of the 834 individuals, 501 were saplings and 333 were shrubs (
Supplemental Table S1). On average, there were 14 individuals and 10 species per plot, but species richness and stem density varied greatly across elevations and between habitats. Plots at mid-elevations in the sub-montane forest habitat had the greatest density of individuals and the highest number of species (
Table 1;
Figure 2).
Drought Tolerance metrics. The values of SL, G
min, and TLP for each individual are reported in
Supplemental Table S1 (
Table S1 also includes other metrics derived from the pressure volume curves but not included in the analyses reported here). The community weighted mean values of SL, G
min, and TLP for each plot are reported in
Table 2. There were no significant correlations between the CWMs of SL, G
min, and TLP (
Figure 3). Plot mean TLP varied from -0.36 MPa to -1.42 MPa (lower values indicate higher drought tolerance) (
Table 2). Due to logistical constraints and the difficulties of measuring TLP in extremely remote locations lacking infrastructure, TLP was only measured in 28 plots from mid and high elevations (937 - 3604 m asl; i.e., only in sub montane forest, cloud forest, and gallery forest habitats). Mean G
min varied between plots from 4.18 to 99.11 mmol m
-2 s
-1 (lower values indicate higher drought tolerance) (
Table 2). As with TLP, G
min was only measured in mid and high elevations (1189 - 3604 m asl). Mean SL varied markedly among plots from 13% to 81% (
Table 2). SL was measured in all plots 400 - 3615 m asl.
There were no significant relationships between TLP or G
min with either elevation and/or mean annual precipitation (
Figure 4). In contrast, Solute leakage (SL) was significantly related to elevation (SL ~ 65.94 - 0.0067* Elevation (in m); P<0.00005; Adjusted R
2 = 0.23) and especially MAP (SL ~ 27.23 + 0.0058 * MAP (in mm); P<0.00005; Adjusted R
2 = 0.50) such that SL was lowest (indicating greater drought tolerance) at high elevations and in areas with lower precipitation (
Figure 4). If both elevation and MAP are included as explanatory variables, only MAP is significant. The significant relationships between SL and precipitation remains even if we limit the analysis to the subset of plots where we also have information on TLP and G
min (SL ~ 25.03 + 0.0062 * MAP; P<0.00005; Adjusted R
2 = 0.55).
DISCUSSION
We assessed changes in the relative drought tolerance of nearly 60 understory woody plant communities along a 3000+ m elevational gradient in the “super humid” tropics by testing the relationships of three different metrics of physiological leaf drought tolerance with precipitation and elevation. In total, we measured at least one measure of drought tolerance on almost 800 individual plants and measured all three measures of drought tolerance on almost 100 individual plants (
Supplemental Table S1). Notably, we did not find covariation between our three metrics of drought tolerance at either the plot level or the individual level (although there was a marginally significant correlation at the individual level between SL and G
min) (
Figure 3). This suggests that plants may develop different independent mechanisms of tolerating drought stress. Additional research is required to understand the interactions and tradeoffs in the different drought tolerance mechanisms employed by understory species and the implications of these mechanisms for plant performance.
The turgor loss point (TLP) values that we measured were generally higher (less drought tolerant) than those reported from other functional groups and habitats. Specifically the TLP of cloud forest vascular epiphytes (-2.71 to -1.67 MPa) [
30], understory shrubs from a tropical rainforest (- 2.55 to -1.43 MPa [
57]), cloud forest trees (-2.04 to -1.34 MPa [
58]), and trees from tropical rainforest (-2.18 to -1.02 MPa [
59]) are all lower than what we measured in understory saplings and shrubs along our tropical montane gradient. Other studies have also found TLP to be related to precipitation and VPD such that plants generally have lower TLP (indicating greater drought tolerance) in drier places [
60,
61]. The relatively high TLP and the absence of a relationship with either elevation or MAP along our gradient (
Figure 4) may indicate that our study area was wetter than those of previous studies either due to higher rainfall, frequent cloud immersion, or higher soil water retention. In addition, understories are generally cooler and more humid than canopies, which may explain the lower TLP values for trees and epiphytes. Alternatively, the plant communities that we studied may be relying more heavily on other mechanisms to tolerate periods of low water availability.
Similarly, we found a no relationship of cuticular conductance (G
min) with either precipitation or elevation (
Figure 4). The lowest values of G
min were in the high elevation gallery forests and the highest values of G
min occurred in the cloud forest plots at around 2500 m asl. Cloud forests are wet habitats with high rainfall and near-constant fog immersion – conditions that may not select for traits conferring drought tolerance. High cuticular conductance may also facilitate foliar water uptake in the cloud forests [
62].
Cell membranes are one of the first victims of many plant stresses and it is generally accepted that the maintenance of their integrity and stability under water stress conditions is a major component of resistance to environmental stresses in plants [
51]. As such, measurements of solute leakage (SL) under desiccating conditions should provide valuable information about plant strategies for dealing with low water availability and stress. Indeed, SL has been previously linked to plant tolerances to drought [
63], high temperatures [
64], and freezing [
65]. In our study, we found high SL in plots in wetter areas and at lower elevations, while understory plant communities located in areas with low precipitation and low temperature showed better membrane stability (low values of SL) (
Figure 4). High-elevation shrubs had the lowest SL (= highest drought tolerance) in the gradient, possibly due to the combination of the dry and cold (sometimes freezing) conditions that can occur in these habitats. Lowland and high elevation gallery forests (400 and 3548 - 3615 m asl, respectively) had intermediate SL. Gallery forests are located between or below mountain ridges and may receive less fog because of topographical configuration, resulting in relatively dry conditions and potentially explaining their high drought tolerance relative to other habitats. Lowland forests experience less rainfall and thus exhibit values of SL that are comparable with gallery forests >3000 m higher on the gradient. The highest SL value was exhibited by sub montane forests around 1000 m asl. The plots in these sub montane forests receive the highest precipitation (5362-5992 mm yr
-1) and can be extremely humid. The high SL of plants in the understory plants of the sub montane forests suggests that these communities are adapted to super wet conditions and may therefore suffer due to any future climate change-driven decreases in rainfall or humidity.
It is particularly noteworthy that SL had the strongest relationships with the environmental variables. Solute leakage has previously been used in laboratory and agricultural assays, but only a few very recent studies have used it to assess drought tolerance of plants in the field. Fadrique et al. (2022) worked along the same elevational gradient in Peru and measured SL in bamboos growing at different elevations. Like our understory dicots, bamboo SL decreased at higher elevations and under wetter conditions [
66]. Zuleta et al. (2022) analyzed intra and interspecific differences of SL in trees growing in different topographic habitats in the lowland Amazon rainforest. Consistent with the other results, Zuleta et al. found that trees growing in wetter microhabitats had higher SL than trees growing on the drier ridge tops [
46]. Based on these previous studies and the new results presented here, it does appear that differences in water availability lead to predictable changes in SL, indicating that SL may be a powerful tool for assessing plant water relations. The potential power of SL is especially exciting given that it is much cheaper and easier to measure in the field than TLP, G
min or other standard measures of drought tolerance (e.g., vascular conductance, leaf embolisms) since it does not require any specialized equipment or facilities.