Interannual Variation in Gas Exchange and Leaf Anatomy in Cenostigma pyramidale Is Exacerbated through El Niño and La Niña Climate Events

: This study sheds light on the remarkable physiological adaptations that the Cenostigma pyramidale makes, particularly during periods of extreme water scarcity, and their remarkable ability to rebound when a new rainy season arrives. C. pyramidale employs a survival strategy to mitigate the adverse effects of water shortage and then endures challenging environmental conditions and sustaining vital functions. Previously, this species was classified as deciduous since it retained at least 40% of its leaves to sustain basal gas exchange rates. The discrepancy in classification arises from differences in study methodologies, with this research being conducted in the field under natural conditions. This study demonstrates a negative correlation between rainfall and specific leaf area (SLA), highlighting that plants with smaller SLA are more drought-tolerant. Changes in leaf anatomy, including an increase in palisade parenchyma and reduction in spongy tissue, serve as adaptive strategies to enhance photosynthesis under water stress conditions. Hydraulic conductance plays a crucial role in plant adaptation to water scarcity. An intricate interplay between leaf anatomy and hydraulic conductance is observed, with adjustments in xylem characteristics influencing leaf gas exchange. The phenotypic plasticity is high in C. pyramidale , demonstrating the species’ ability to adapt to changing environmental conditions. In summary, this study illuminates the multifaceted strategies employed by plants to cope with water scarcity, from leaf shedding to anatomical and physiological adaptations, highlighting the resilience of native species in arid environments. These findings offer valuable insights into plant responses to environmental stress and their ability to thrive in challenging conditions.


Introduction
Industrial emissions of greenhouse gases contribute to global warming.As the Earth's temperature continues to rise, we witness shifts in climate patterns, including heightened temperature extremes and significant changes in both the timing and geographic distribution of rainfall, ultimately leading to prolonged periods of drought.Projections indicate a growing occurrence of temporary spikes in temperature surpassing the typical physiological limits and more frequent, intense drought periods in the future [1].The levels of climate-related risk frequently escalate in tandem with rising temperatures, and these risks are occasionally intertwined with various facets of climate change.Today, we are experiencing a pivotal example of the simultaneous occurrence of a substantial reduction in rainfall, along with isolated episodes of torrential rainfall in specific regions [2,3].This confluence of events can have profound implications for the broader climate landscape, exacerbating vulnerabilities and challenges associated with changing weather patterns and their associated impacts.With global climate change, it becomes increasingly imperative to gain deeper insights into the alterations in both plant physiological and morphological characteristics.This knowledge is not only pivotal for the conservation of dry tropical forests but also for addressing the broader environmental concerns we face today [4].In this context, El Niño is a climate phenomenon characterized by the periodic warming of ocean waters in the central and eastern tropical Pacific, leading to significant weather disruptions worldwide.As our planet's temperatures persist in their upward trajectory, the repercussions are becoming ever more conspicuous and all-encompassing [5,6].With the exacerbation of global climate changes, phenomena like El Niño and La Niña may become more frequent and affect larger geographic areas [7].Climate scenarios suggest that these global climate changes will disproportionately harm arid and semiarid regions, resulting in prolonged dry periods and posing significant risks to biodiversity [5].The significant increase in CO 2 emissions poses a grave threat to food and human security, amplifying the risk of severe, pervasive, and irreversible impacts on both ecosystems and people.This underscores the urgent need to address the challenges of global climate change, which can disrupt the delicate balance of our planet's well-being [4].
The interplay between drought conditions and intense solar radiation poses a substantial challenge to the efficiency of net photosynthesis (P N ) and the potential risk of leaf overheating, particularly when decreased transpiration restricts the leaf's ability to cool itself due to a scarcity of water [8,9].The drought stress is intricate, involving a spectrum of changes in their morphology, physiology, as well as primary and secondary metabolic pathways [10].In response to drought stress, plants activate intricate defense mechanisms to shield themselves against oxidative stress, which encompasses the loss of turgor, osmotic stress, and the closure of stomata to curtail water loss.Water stress also hampers photosynthesis, causing alterations in chlorophyll concentrations and causing damage to the photosynthetic apparatus [11].Additionally, a common occurrence is the reduction in enzyme activity within the Calvin Cycle [12][13][14].In tropical environments where water stress is prevalent, plants undergo significant adaptations.Their struggle for survival manifests through several pivotal adjustments in their photosynthetic processes.An initial indicator of stress is the disruption of stomatal regulation over leaf gas exchange, resulting in a diminished efficiency in carbon dioxide uptake and inhibition of photosynthesis [11,[15][16][17].The constraint on P N due to water stress can imperil the plant, especially under conditions of excessive irradiation.Failure to safely dissipate this excess energy can result in the overreduction of reaction centers and heightened production of reactive oxygen species (ROS) within the chloroplasts, ultimately causing harm to photosystem (PS)II [18].Consequently, the damage to photosynthetic membranes, crucial for the conversion of light energy into chemical energy, becomes an unfortunate outcome of the plant's struggle for survival [19].
The correlation between rising temperatures and the risks linked to climate change emphasizes the pressing need for a comprehensive approach to tackle this issue.There is a consensus in the scientific community that plant species thriving in arid and semiarid environments have evolved specific adaptations to endure periods of water scarcity [20,21].When faced with acute water shortages, regardless of whether they are exacerbated by increasing atmospheric temperatures or not, these plants often employ survival strategies.One such tactic involves shedding leaves as a means to reduce water loss to the atmosphere, enabling them to conserve precious moisture resources [22].In general, gas exchange is intricately governed by a spectrum of factors, encompassing both external environmental conditions and internal physiological aspects within the plant [23].During episodes of water stress, a typical response involves reduced water transport from the roots to the leaves, which consequently prompts a reduction in stomatal opening for CO 2 uptake [15,[24][25][26][27].Given that CO 2 serves as the primary substrate for photosynthesis, this diminished uptake often translates into a concurrent decrease in photosynthetic rates.Nevertheless, it's important to note that other contributing factors may exist, some of which could exert an equal or even greater influence on the regulation of photosynthesis rates.A simple rise in stomatal CO 2 absorption doesn't automatically lead to an increase in gas exchange.This is due to the fact that CO 2 must traverse several physical barriers before it can reach the RuBisCo carboxylation site [28].Nevertheless, various other factors can influence photosynthesis.In addition to physiological and biochemical responses, these plants have developed mechanisms to adapt to water stress through alterations in the anatomical structure of their leaves and roots, as well as modifications in specific characteristics [29,30].Leaf architecture plays a crucial role in facilitating CO 2 exchange within the mesophyll, making the study of plant leaf anatomy and its xylem vessels particularly intriguing.As the plant responds to decreased osmotic potential, it tends to curtail its water flux, impacting all metabolic processes, not limited to photosynthesis alone [27].
As sessile organisms, plants face a wide range of abiotic stresses, with water stress being a significant condition affecting their growth, development, and productivity, particularly in arid and semiarid regions [27].Caatinga has a very high plant diversity and rate of endemism but is also a region that is quite degraded because of the unsustainable use of its natural resources; human occupation is followed by abandonment and subsequent natural regeneration of secondary forests [31].Cenostigma pyramidale is considered one of the most frequent plant species in the Caatinga [32].In this scenario, physiological, biochemical, and morphological attributes play an important role in abiotic and biotic responses to environmental conditions [32].Moreover, plant sensitivity to adverse environmental conditions is majorly dependent on the severity of the stress as well as plant species, genotype, and developmental stage.Plant tolerance to periodic water deficit is expected to be of increasing importance as climate change intensifies inter-and intra-annual rainfall extremes [33][34][35].The responses of native plants to extreme climatic events, such as rising temperatures coupled with reduced rainfall, remain poorly understood [36], making their survival in one of the harshest environments truly remarkable [37].Substantial efforts have been made to better comprehend the drought tolerance strategies of the trees in Caatinga [12,22,31,[37][38][39][40].In this study, we evaluated the phenotypic plasticity in the physiological and anatomical characteristics of individuals of the Cenostigma pyramidale, a Fabaceae species.We hypothesize that C. pyramidale displays increased resilience when subjected to recurring water shortages.In this sense, this work can contribute to the understanding of the plant physiology of a symbolic species of the Brazilian Caatinga since the analyses were carried out over two years, encompassing a very dry year and another year with rainfall distributed according to the historical average and standard for those evaluated months.

Plant Species, Site Description, and Experimental Design
Cenostigma pyramidale (Tul.)Gagnon & G.P. Lewis (Fabaceae) is a Caatinga's symbol that predominates in all landscapes of the ecosystem.This tree is semideciduous and relevant as a pioneer species during the ecosystem succession.C. pyramidale, popularly known as "Catingueira," is a resilient tree thriving in the Caatinga, a challenging semiarid environment.Recognizable by its distinctive thorny branches, this deciduous species acts as vital in the ecological rejuvenation of the region.Its adaptability and cultural significance make it a key emblem of the Caatinga landscape.In this study, we evaluate the effects of climatic changes on the morphophysiological characteristics of this species over two distinct years.For this purpose, we selected ten adult trees (approximately 25 years of age) in the reproductive phase ranging from 2.5 to 4.0 m in height.All the plants were identified by a taxonomist and deposited in a plant collection at the Agronomic Institute of Pernambuco Herbarium (registered by IPA-88.494).
The Caatinga, a Brazilian semiarid forest covering an area of 862,818 km 2, constitutes approximately 10.1% of the country's landmass [41] (Figure 1).Despite being one of the most biodiverse dry forests globally [42], it remains one of the least studied and protected ecosystems within Brazil [31].Within the Caatinga, the annual rainy season typically extends for 4 to 6 months, marked by sporadic and meager rainfall patterns (Figure 2) with precipitation ranges around 250 and 700 mm per year, concentrated in 4 months, a few years with precipitation more than 1000 mm [43].These conditions result in high evaporation rates and limited water-holding capacity within the soil.Consequently, Caatinga trees must continually adapt to environmental challenges, evolving diverse survival tactics to vie for the scarce water supply [37].
ecosystems within Brazil [31].Within the Caatinga, the annual rainy season typically extends for 4 to 6 months, marked by sporadic and meager rainfall patterns (Figure 2) with precipitation ranges around 250 and 700 mm per year, concentrated in 4 months, a few years with precipitation more than 1000 mm [43].These conditions result in high evaporation rates and limited water-holding capacity within the soil.Consequently, Caatinga trees must continually adapt to environmental challenges, evolving diverse survival tactics to vie for the scarce water supply [37].ecosystems within Brazil [31].Within the Caatinga, the annual rainy season typical tends for 4 to 6 months, marked by sporadic and meager rainfall patterns (Figure 2 precipitation ranges around 250 and 700 mm per year, concentrated in 4 months, years with precipitation more than 1000 mm [43].These conditions result in high e ration rates and limited water-holding capacity within the soil.Consequently, Caa trees must continually adapt to environmental challenges, evolving diverse surviv tics to vie for the scarce water supply [37].As shown in Figure 2, the climate conditions in 2020 were very distinct from 2019 [43].Consequently, in 2019, we experienced a characteristic dry year, with a tendency to a rainy season from February to April, with a higher rainfall in March (181.2mm).In 2020, we verified a median rainfall as compared to expected rainfall from this area, with a rainy season extending from January to May, with the peak rainfall occurring in March 2020 (347.3 mm).Similarly, the dry season spanned from June to December 2020, with a tendency for the month to be less dry in November (89.6 mm) (Figure 2).These patterns were consistent not only in terms of precipitation but also with respect to global radiation and relative humidity.The maximum and minimum levels of global radiation were recorded in September 2020 (65.2 MJ m −2 day −1 ) and March 2020 (7.4 MJ m −2 day −1 ).Air humidity ranged from 17.3% (August 2020) to 51.7% (March 2020), reflecting the variations in moisture content in the atmosphere.Meanwhile, air temperatures fluctuated between 31.8 • C (March 2020) and 38.5 • C (August 2020), representing the seasonal temperature changes over the year.
The experiments were conducted in a completely randomized design in an open field environment, specifically within a Caatinga natural fragment in Serra Talhada (7 • 58 ′ 57" S; 38 • 18 ′ 16" W; 418 m a.s.l.), Pernambuco, Brazil (Figure 1).We selected a total of 10 Cenostigma pyramidale plants for our investigation.All of them were measured between 3 and 5 m in height, with breast diameter of 160 mm and plant age of at least 20 years.We collected data for the rainy (March) and the dry (August) seasons in both 2019 and 2020.In all experiments, to guarantee that the sampled leaves were developed in the rainy and dry seasons, monthly leaf primordia were marked with non-obstructive adhesive tape so that their leaf development could be observed and guaranteed between one collection and another.

Leaf Area and Specific Leaf Area
In January 2019, leaves primordia were marked with non-obtrusive adhesive tape in order to monitor the development of the leaves until its complete expansion.With each collection carried out, new leaf primordia were marked with adhesive tape to monitor their development.This procedure is necessary to ensure that the leaf that will be collected later was necessarily expanded during the evaluation period and to guarantee that the leaf that will be collected was not already on the mother plant before the treatments were imposed.Thus, in March and August 2019 and in March and August 2020, completely expanded leaves were collected from the plants after the gas exchange experiments described in topic 2.3 were completed.These leaves were inserted into zip-lock plastic packaging and transported under refrigeration to the laboratory.The leaves were digitalized using a scanner (HP Smart Suite, 1200 × 1200 dpi; Hewlett-Packard Company, Palo Alto, CA, USA), and the images were processed using Image-Pro ® Plus software (version 4.5, Media Cybernetics, Silver Spring, MD, USA).Then, five fragments 1 cm in diameter were collected, wrapped in an aluminum foil envelope, and placed in an oven at 70 • C for 72 h, when their dry mass was measured on an analytical balance (Sartorius Analytical Balance mod.ENTRIS224-1S, Bradford, MA, USA; accurate to 0.1 mg).The specific leaf area was obtained by the ratio of the dry mass to the leaf area of each leaf fragment, obtained by Equation (1): where DW is the dry weight of the leaf fragments, π is 3.14159, and r is the ray of the circumference, in this case 0.5 cm.

Predawn and Noon Leaf Water Potential (Ψw)
To assess leaf water status, we selected and collected one fully expanded leaf (i.e., from the third node from the apex) per tree.We measured the leaf water potential (Ψ w ) at predawn (5:00 a.m.) and noon (12:00 p.m.) with a Scholander-type pressure chamber [44].

Gas Exchange
In March and August, for both 2019 and 2020, leaf gas exchange was determined on the 3rd attached fully expanded leaf from the apex, using a portable open-flow infrared gas analyzer (LI-6400XT; LI-COR Inc., Lincoln, NE, USA).The net photosynthesis (P N , µmol CO 2 g −1 s −1 ), stomatal conductance (g s , mol H 2 O g −1 s −1 ), transpiration (E, Forests 2024, 15, 346 6 of 28 mol H 2 O g −1 s −1 ), internal-to-ambient CO 2 concentration ratio, leaf temperature (T leaf , • C), and vapor pressure deficit (VDP, kPa) were obtained on 20-year-old plants as described by Pinho-Pessoa et al. [36].All measurements were conducted from 6 a.m. to 6 p.m. (GMT-3) under a cloudless sky and with the leaf irradiance of saturation of 1200 µmol m −2 s −1 (as previously tested by light curves versus P N ; [36]), fixed CO 2 concentration in 390 µmol mol −1 and airflow of 400 µmol s −1 [25].Diurnal integrated CO 2 assimilation (integrated P N ) was obtained by integrating all measurements of P N measured throughout the day, as previously described in Mendes et al. [22].With gas exchange parameters, we made a correlation between P N and g s , P N and E, and P N and internal CO 2 concentration (Ci).WUEi was obtained by Equation (2).
Also, VPD was correlated with P N , g s , and E.

Leaf Anatomy
Leaf fragments 5 cm 2 of the 3rd attached fully expanded leaf from the apex (the same used in the previous topic) were sampled at the same time as described above.In each plant, two leaf fragments were collected, promptly immersed in FAA 50% for 48 h, and then stored in 70% (v/v) ethanol until analysis.Leaf samples were dehydrated in an ethylic series and embedded in plastic resin (HistoResin Leica Microsystems Nussloch, Heidelberg, Germany, part number 7592).Following this, all samples were processed as described in detail in Mendes et al. [40].For each sample, ten images were captured by a digital camera (Mikrosysteme Vertrieb GmbH, model ICC50 HD; Wetzlar, Germany) interfaced with a computer.To estimate the potential hydraulic conductivity (Kp), the principles of the Hagen-Poiseuille equation were implemented using Equation (1).
where Kp is the potential specific stem conductivity, h is the viscosity of water at 20 • C (1.002 × 10 −3 Pa), p w is the density of water at 20 • C (998.2 kg m −3 ), VD is the vessel density, and Dh is the hydraulically weighted vessel diameter (in m).Since the vessels were not exactly circular, the diameter of each vessel was calculated as the mean of the minimum and maximum diameters.The average Dh was calculated as well using the methodology described in Sterck et al. [45].As recommended by Scholz et al. [46], we measured at least 50 vessels per photomicrograph, resulting in 500 vessels (10 repetitions) per treatment.

Phenotypic Plasticity
The phenotypic plasticity ranges from 0 (any plasticity) to 1 (full plasticity) was calculated as the difference between the average maximum and minimum measurements for each sample divided by the maximum value [47] for all physiological and anatomical features.

Experimental Design and Statistical Analyses
The experiments were conducted in a completely randomized block design and ten replicates.All the data were analyzed by Two-way ANOVA, and means were compared using an SNK test (p < 0.05) by Statistic version 14.0 (StatSoft, Tulsa, OK, USA).The principal component analysis was estimated after a multivariate analysis for all analyzed features in Minitab 18.1 (Minitab, Inc., Chicago, IL, USA).Heat maps were used to compare the mean of each treatment, using the control (March, the rainiest month) as a reference.After the log 2 transformation, the false color method was used, including a color scale.

Leaf Area (LA) and Specific Leaf Area (SLA)
Tables 1 and 2 show that the specific leaf area (SLA) of C. pyramidale varies widely with the levels of rainfall received during its development, with a positive interaction for both year and seasonality (Table 1).Furthermore, we also show that the leaf area increased in the same proportion as rainfall increases, with leaves 40% and 67% larger in the rainy months (March, both in 2019 and 2020, respectively) compared to the driest months (August).The SLA ranges from 11.37 cm 2 g −1 to 15.07 cm 2 g −1 , while leaf area ranges from 1.59 cm 2 to 2.43 cm 2 (Table 2).  1.Each value denotes means (±SD).The values followed by different small case letters denote significant differences between seasons in the same year, and the means followed by uppercase letters denote significant differences between years in the same season.
(p ≤ 0.001).Knowing that SLA varies widely seasonally, we expect to find more photosynthetic tissues in thicker leaves, that is, in leaves with lower SLA (in the dry season).Therefore, presenting gas exchange data per unit of leaf area would not be the best option.So, in this study, we will present gas exchange data on a gram basis of leaf tissue and not by metric area.

Water Potential
The predawn Ψ w measurements from 2019 revealed significant plant water stress, with readings of −5.6 MPa and −4.2 MPa in March and August, respectively (Figure 3).In 2020, these values saw a notable improvement, with predawn Ψ w readings being 2.6 and 1.4 times higher for the same months in 2019.Interestingly, while predawn Ψ w values exhibited this substantial shift between the two years, the midday Ψ w values remained relatively consistent, showing less variability in both 2019 and 2020 (Figure 3).
1.4 times higher for the same months in 2019.Interestingly, while predawn Ψw values exhibited this substantial shift between the two years, the midday Ψw values remained relatively consistent, showing less variability in both 2019 and 2020 (Figure 3).2019) and a normally rain year (2020) in both seasonal conditions (rainy, March and dry, August).The measurements were made at predawn (5 a.m.) or at midday (12 p.m.).Means followed by small case letters denote statistical significance between seasons at the same time (5 a.m. or 12 p.m.), means followed by uppercase letters denote statistical significance between time and the same season (rainy and dry months) and the same year, and asterisks (*) denotes statistical difference between years in the same months and time.All values denote media ± SD. n = 10.

Gas Exchanges
In both 2019 and 2020, all gas exchange parameters fluctuated throughout the day.Net photosynthesis (PN), stomatal conductance (gs), and transpiration (E) start the day with higher values, and these were reduced throughout the day.In 2019, this pattern was very soft for PN, while for gs and E, the daily variation is not statistically different (Figure 4).From another perspective, in 2020, for both the rainy (March) and the dry (August) seasons the PN (at 6 a.m.) is measured at 213.8 µmol CO2 g −1 s −1 (Figure 4C) and 122.2 µmol CO2 g −1 s −1 (Figure 4D) and fall dawn to 44.2 µmol CO2 g −1 s −1 (−79.3%; Figure 4C) and 40.6 µmol CO2 g −1 s −1 (−66.8%; Figure 4D) at 6 p.m.Moreover, when the comparisons are made between seasons for both years, we show that the PN registered at 10 a.m. during the rainy season was 10.2-fold higher in 2020 than those values recorded in 2019.During the dry season, the PN at 10 a.m.2020 was 6-fold higher in 2020 than those values recorded in 2019.

Gas Exchanges
In both 2019 and 2020, all gas exchange parameters fluctuated throughout the day.Net photosynthesis (P N ), stomatal conductance (g s ), and transpiration (E) start the day with higher values, and these were reduced throughout the day.In 2019, this pattern was very soft for P N, while for g s and E, the daily variation is not statistically different (Figure 4).From another perspective, in 2020, for both the rainy (March) and the dry (August) seasons the P N (at 6 a.m.) is measured at 213.8 µmol CO 2 g −1 s −1 (Figure 4C) and 122.2 µmol CO 2 g −1 s −1 (Figure 4D) and fall dawn to 44.2 µmol CO 2 g −1 s −1 (−79.3%; Figure 4C) and 40.6 µmol CO 2 g −1 s −1 (−66.8%; Figure 4D) at 6 p.m.Moreover, when the comparisons are made between seasons for both years, we show that the P N registered at 10 a.m. during the rainy season was 10.2-fold higher in 2020 than those values recorded in 2019.During the dry season, the P N at 10 a.m.2020 was 6-fold higher in 2020 than those values recorded in 2019.
The daily fluctuations of g s and E were similar to P N .For g s , the values ranged from 151. 5 ), and from 20.9 mmol H 2 O g −1 s −1 to 7.4 mmol H 2 O g −1 s −1 (−64.8%).At 10 a.m., during March and August 2020, it measured 7-and 6.4-fold higher g s Forests 2024, 15, 346 9 of 28 than in 2019, while E registered 26.6-and 22.0-fold higher in 2020 when compared to the 2019 values (Figure 4).As described above, the comparisons of net photosynthesis are time-dependent.Thus, when comparing the daily integrated PN values between seasons and years (Figure 5), the differences are still more evident because those values were 6.6-fold higher in March 2020 as compared to March 2019 or 4.9-fold in August 2020 when compared to August 2019 (Figure 5).Between seasons, within the same year, we verified a significant effect only in 2020, when the daily integrated PN values were 37.1% larger in March than in August. in Cenostigma pyramidale grown in field conditions in the driest year (2019) and a normal rain year (2020) in both seasonal condition (rainy, March and dry, August).Means followed by small case letters denote statistical significance between months, regardless of year, and means followed by uppercase letters denote statistical significance between years in the same month.All values denote media ± SD. n = 10.
As described above, the comparisons of net photosynthesis are time-dependent.Thus, when comparing the daily integrated P N values between seasons and years (Figure 5), the differences are still more evident because those values were 6.6-fold higher in March 2020 as compared to March 2019 or 4.9-fold in August 2020 when compared to August 2019 (Figure 5).Between seasons, within the same year, we verified a significant effect only in 2020, when the daily integrated P N values were 37.1% larger in March than in August.
The correlations (Supplementary data file) of both g s and E could be modulated P N since correlations between P N and g s or P N and E were, respectively, 0.860 and 0.908 (both p < 0.00001).Therefore, to confirm that g s is not the preponderant factor mediating photosynthesis, we show in Figure 6A,B a regression between P N and g s and P N and E. The way of interpretation was season-dependent to contribute to correlations.In this case, we verify that all regressions between P N and E (Figure 6B) return higher regression coefficients than regressions between P N and g s (Figure 6A).As regression coefficients are high, we perform a regression between P N and Ci to infer the instantaneous ribulose 1,5-diphosphate carboxylase/oxygenase carboxylation efficiency.The analysis of this variable shows us that in the rainy months, the P N /Ci regression curve tended to be more responsive than in the dry months, although none of the curves were significant.However, when Ci is analyzed in an integrated manner with the other variables (multivariate analysis, Figure 10), Ci was highly responsive to photosynthesis and responsible for the separation of variables that led to the distinction of the August 2020 cluster (Details, see Figure 10).All this data confirms that in C. pyramidale, the P N is limited by stomatal conductance.This idea is confirmed by analysis of the Ci:Ca ratio throughout the day for all seasons (Figure 7A-D).We verify that all records from 2019 return a Ci:Ca ratio lower than 0.5 (dotted line in Figure 7A-D).In 2020, all values of the Ci:Ca ratio were higher than 0.5, with the exception of the record at noon in March 2020, during an effect that in plant physiology is called midday depression.The correlations (Supplementary data file) of both gs and E could be modulated PN since correlations between PN and gs or PN and E were, respectively, 0.860 and 0.908 (both p < 0.00001).Therefore, to confirm that gs is not the preponderant factor mediating photosynthesis, we show in Figure 6A,B a regression between PN and gs and PN and E. The way of interpretation was season-dependent to contribute to correlations.In this case, we verify that all regressions between PN and E (Figure 6B) return higher regression coefficients than regressions between PN and gs (Figure 6A).As regression coefficients are high, we perform a regression between PN and Ci to infer the instantaneous ribulose 1,5-diphosphate carboxylase/oxygenase carboxylation efficiency.The analysis of this variable shows us that in the rainy months, the PN/Ci regression curve tended to be more responsive than in the dry months, although none of the curves were significant.However, when Ci is analyzed in an integrated manner with the other variables (multivariate analysis, Figure 10), Ci was highly responsive to photosynthesis and responsible for the separation of variables that led to the distinction of the August 2020 cluster (Details, see Figure 10).All this data confirms that in C. pyramidale, the PN is limited by stomatal conductance.This idea is confirmed by analysis of the Ci:Ca ratio throughout the day for all seasons (Figure 7A-D).We verify that all records from 2019 return a Ci:Ca ratio lower than 0.5 (dotted line in Figure 7A-D).In 2020, all values of the Ci:Ca ratio were higher than 0.5, with the exception of the record at noon in March 2020, during an effect that in plant physiology is called midday depression.2019) and a normally rain year (2020) in both seasonal conditions (rainy, March, and dry, August).Means followed by small case letters denote statistical significance between months, regardless of year, and means followed by uppercase letters denote statistical significance between years in the same month.All values denote media ± SD. n = 10.
The intrinsic water use efficiency (WUEi; Figure 7E-H) fluctuates throughout the day, as do other gas exchange parameters.However, in the rainy season, the WUEi increased throughout the day, while in the dry season, the value measured at 6 p.m. was higher than that measured at 6 a.m.Regardless of season or evaluation time, the WUEi was higher in all analyses since this feature is dependent on net photosynthesis and stomatal conductance.This statement is confirmed by a negative correlation between rainfall and WUEi (r = −0.224)and a positive correlation between P N and WUEi (r = 0.215).
Along the day, the leaf temperature (Figure 7I-L) returns values dependent on air temperature and transpiration rate.As a consequence, T leaf shows a negative correlation between Ψ w(md) (r = −0.505),P N (r = −0.692),daily P N (r = −0.748),g s (r = −0.669),and E (r = −0.744).These patterns confirm that leaf transpiration decreases leaf temperature and improves net photosynthesis.
The vapor pressure deficit (VPD; Figure 8) also reflects the effect of air temperature with lower values in the early morning and late afternoon and higher values registered between 10 a.m. and 12 p.m. Commonly, VPD has an influence on the P N , g s , and E, because it is atmospheric water potential that governs VPD, and by analogy, P N , g s , and E are shown with a correlation between VPD and P N (r = −0.624),g s (r = −0.632),and E (r = 0.678).However, when correlations were made with these parameters (Figure 9), only E returned a weakly significative correlation of coefficients (Figure 9C).The intrinsic water use efficiency (WUEi; Figure 7E-H) fluctuates throughout the day, as do other gas exchange parameters.However, in the rainy season, the WUEi increased throughout the day, while in the dry season, the value measured at 6 p.m. was higher than that measured at 6 a.m.Regardless of season or evaluation time, the WUEi was higher in all analyses since this feature is dependent on net photosynthesis and stomatal conductance.This statement is confirmed by a negative correlation between rainfall and WUEi (r = −0.224)and a positive correlation between PN and WUEi (r = 0.215).
The vapor pressure deficit (VPD; Figure 8) also reflects the effect of air temperature with lower values in the early morning and late afternoon and higher values registered between 10 a.m. and 12 p.m. Commonly, VPD has an influence on the PN, gs, and E, because it is atmospheric water potential that governs VPD, and by analogy, PN, gs, and E are shown with a correlation between VPD and PN (r = −0.624),gs (r = −0.632),and E (r = 0.678).However, when correlations were made with these parameters (Figure 9), only E returned a weakly significative correlation of coefficients (Figure 9C).

Plant Anatomy
C. pyramidale leaflets consist of uniseriate epidermises, with anticlinal outer cell walls covered by a thin cuticle (Figure 10), and unicellular tectorial trichomes, stomata occurring on the abaxial epidermal surface, with glandular trichomes on the abaxial surface.The mesophyll is dorsiventra,l consisting of 1-2 layers of palisade parenchyma (PP) and 3-4 layers of spongy parenchyma (SP).Vascular bundles (VB) are collateral and associated with pericyclic fibers (Figure 9B).The secretory structures are located in the spongy parenchyma with fibers of the pericycle.
The midrib is biconvex, and its epidermis is uniseriate with a thin cuticle, also containing a few tectorial trichomes.
In this paper, we show that the leaf thickness varies seasonally.In 2019, there were no significant differences in mesophyll thickness.However, in 2020, the leaves sampled during the dry season were, on average, 25% thicker than the leaves sampled in the rainy season plants (Table 3).An impressive part of this increase is due to the greater thickness of the PP, which was 71% greater than during the dry season.Adaxial epidermis thickness (ADAE) (19.9%) and abaxial epidermis thickness (ABAE) (27.9%) contributed little to this thickness, while SP contributed negatively (−18.1%) to the increase in the mesophyll thickness (Table 3).

Phenotypic Plasticity
Phenotypic plasticity denotes the degree of modulation of the physiological and anatomical features when exposed to a specific environment situation.Table 5 shows us that physiological plasticity ranged from 0.077 (w(md)) to 0.995 (gs).

Phenotypic Plasticity
Phenotypic plasticity denotes the degree of modulation of the physiological and anatomical features when exposed to a specific environment situation.Table 5 shows us that physiological plasticity ranged from 0.077 (Ψ w(md) ) to 0.995 (g s ).
Table 5.The plasticity index of physiological and anatomical features in the leaves of Cenostigma pyramidale plants grows in a natural environment of semiarid Brazilian Caatinga in 2019/March to 2020/August in accordance with historical rainfall.Means followed by smallcase letters denote statistical differences between samples in the same feature (Newman-Keuls test p ≤ 0.001), and uppercase letters denote statistical differences between features (Bonferroni's test at p ≤ 0.001).

Physiological Features
March However, the high amplitude of the standard deviation provokes a non-significant mean value for other factors, but the plasticity index of March 2020 was 37.1% higher than those media presented in March 2019.An analogous situation occurred in anatomical plasticity, where the plasticity index ranged from 0.025 (SLA) to 0.900 (Kp).As a result, the amplitude of standard deviation does not provide a statistical difference between anatomical plasticity.Moreover, the mean value of the physiological plasticity was 0.567 ± 0.223, while an anatomical plasticity of 0.473 ± 0.187 shows a higher physiological plasticity (~20%) than anatomical plasticity (p ≤ 0.01).
Some physiological features present a strong ability to modulate when comparing the dry season with the rainy season, regardless of the year of evaluation.In this sense, the most plastic features are Ψ w(pd) , Ψ w(md) , Ci, Ci:Ca ratio, and T leaf .Anatomical features also showed a degree of plasticity but to a lesser extent.As mentioned above, the anatomical variables also modulate, and the most plastic is leaf area, ADAE, SP, MT, MXT, VA, NV, and Kp.

Discussion
Our findings highlight the pronounced seasonality of rainfall and drought in the Brazilian Caatinga environment [1].C. pyramidale showed a semideciduous pattern since it maintains at least 40% of the leaves on the plant, even during the dry season.This concept disagrees with Santos et al. [2], who classify this species as deciduous.The different patterns described for the same species probably arise from the methodology applied in the study.While Santos et al. [2] studied potted plants under artificial conditions, this study was conducted in the field with the plants in natural conditions.The discrepancy in this concept also is exacerbated by Santos et al. [2], which suspended 100% of the total irrigation of C. pyramidale seedlings for 5 days, but the plants still make photosynthesis, strict evidence of the presence of leaves since in this plants extremely drought stressed Water deficit could reduce leaf growth and, in turn, leaf area, because leaf area expansion depends on leaf turgor, temperature, and assimilating the supply for growth, which can be affected by drought as observed in this study and several other species [22,[48][49][50].In general, SLA is the main component of leaf water content.Leaf tissue density decreases as leaf water content increases and is influenced by precipitation [32].This statement is aligned with our results, which show a higher SLA in the driest year in comparison to the rainiest.The negative correlation between rainfall and SLA (r = −0.231)confirms our Table 2 data.Similar results were previously reported by Falcão et al. [51] in a study with this same species.The SLA of C. pyramidale was higher as the rain occurrence in both evaluated years in distinct forms were presented before [36], where SLA just increased in the rainy year and presented a similar pattern in the driest year.This finding supports our initial hypothesis that C. pyramidale displays increased resilience when subjected to recurring water shortages, given its status as a native species and its robust adaptation to the fluctuating environmental conditions throughout the seasons.This enduring nature was previously reported by Pagotto et al. [52] when this species was studied in the same biome as this study.These results lead us to argue that the increased SLA during the dry year may be indicative of compromised leaf biomass due to extremely low water potentials, which directly affects the turgor pressure required for cell elongation to make it more difficult.It is worth noting that SLA is influenced not only by leaf area but also by leaf density and thickness [53,54].These adaptive strategies permit the plant to conserve water, preventing significant drops in P N and increases in transpiration.Following the restoration of normal water availability in 2020, all plant physiological parameters rebounded strongly.The reduced g s in 2019 registered in the dry season were primarily due to a typical feed-forward response, where stomata closed in response to VPD, confirmed by a strongly significative correlation between VPD and E (r = 0.678) and gs (r = 0.632).Comparing 2020 with 2019, it's apparent that in 2019, the stomata tended to stay less open during the dry season, resulting in reduced daily photosynthesis.These findings suggest that in C. pyramidale, net photosynthesis is more constrained by stomatal conductance than by other parameters, as was previously reported [36].
Wright et al. [55] describe that C. pyramidale has more stable water use efficiency, which represents a background contribution to the evapotranspiration pattern over shorter timescales.Then, transpiration shifts were not detected, and responses were more buffered to VPD.In accordance with these authors, given this hydraulic response and because the phenological timing of this species is cued by rainfall, C. pyramidale is likely more vulnerable to acute drought.In situations of water scarcity, the carbon equilibrium within plant systems is disrupted as a protective measure [56].Consequently, a delicate balancing act unfolds for plants when they confront stress-inducing circumstances, primarily hinging on their ability to strike a balance between sustaining photosynthetic processes and preserving growth.Notably, our findings align with prior research documenting seasonal fluctuations in photosynthesis among Caatinga species and others [36,51,57,58].For instance, a study with Croton blanchetianus reported P N ranging from 15-20 µmol CO 2 m −2 s −1 during the rainy season to less than 5 µmol CO 2 m −2 s −1 in the dry season [22,40].Similarly, Pinho-Pessoa et al. [36] identified variations in the P N , with values ranging from 13.3 µmol CO 2 m −2 s −1 Forests 2024, 15, 346 20 of 28 during the rainy season to 6.7 µmol CO 2 m −2 s −1 in C. pyramidale under water-stressed conditions.
Similar results were described in other Caatinga species, such as Spondias tuberosa [20], Cenostigma pyramidale [36], Croton blanchetianus [22], Jatropha curcas [15], and Cynophalla flexuosa [59].According to Pinho-Pessoa et al. [36], Cenostigma pyramidale displays lower values of P N mediated by stomatal closure.This statement becomes clearer with the analysis of Figure 6, which shows that E is more responsive for P N than g s .The decrease in the Ci:Ca ratio indicates that constraints of g s are dominant under drought conditions, as further evidenced in the evaluations from the driest year.During this period, the C. pyramidale leaves displayed reduced Ci:Ca ratios coupled with diminished P N , implying limitations in carbon assimilation in a direct way.The confluence of drought conditions with an excess of incident radiation can result in a reduction in radiation use efficiency.Additionally, under such circumstances, there is an increased risk of leaf overheating, particularly when transpirational leaf cooling is compromised as a consequence of water deficits [8,9,60].Thus, we can argue that transpiration must modulate gas exchange more intensely than stomatal conductance.This is particularly evident during the dry months when high temperatures and stomatal closure prevail in the region [22,31,36,39,40].Furthermore, it must be considered that CO 2 , when captured by the stomata, must follow a very complex path until it reaches the rubisco carboxylation site, a route called mesophilic conductance [61][62][63].Furthermore, Rubisco and other enzymes of the Calvin-Benson Cycle must also be activated to receive this carbon, process it intrioses phosphates, and send it to the sink organs [64,65].High Ci:Ca ratio in the rainy seasons at the expense of the dry seasons may serve as a complement to our hypothesis that transpiration and stomatal conductance are mainly responsible for governing gas exchange in the dry months of the year.High Ci:Ca ratios also lead to a slower flow of CO 2 through the mesophyll (mesophilic conductance) and/or enzymatic restrictions in the Calvin Cycle.This is a critical point to be considered in future studies.
According to Falcão et al. [57], C. pyramidale can keep its photosynthetic rates even with low leaf water potential.This finding contradicts our results because a distinct reduction of net photosynthesis was found in the dry year compared to a year with a normal rainy season.However, net photosynthesis does not show a constant or fixed net value; instead, it ranges across the measurement of time.So, to give an accurate comparison, these values must be acquired at the same time or integrated in a full day.In this study, we verify that in March 2020, the net photosynthesis was increased in 11.1-, 11.1-, 10.2-, 6.1-, 4.4-, 2.7-, and 2.7-fold higher than the same times evaluated in 2019.For this reason, we preferred to present the data as integrated daily values, which were 6.6-and 4.9-fold higher in 2020 in comparison to 2019, respectively, in both the rainy and dry seasons.The transpiration rate, as part of the gas exchange, increased in March 2020 between 4.7-to 26.6-fold in comparison to evaluations conducted at the same times as in 2019, a range that was lower or higher than those presented only a singular evaluation (17-fold) by Falcão et al. [51].Falcão et al. [51] also describe that WUEi was higher in the driest year (63%) in comparison to the rainy year, a fact that contradicts our results, which, with the exception of those values registered at 8 a.m., all other values were 1-to 1.7-fold higher in the rainy year in comparison to the driest year.A reduction in transpiration via stomatal closure prevents excessive water loss in cases of low water availability and protects the hydraulic architecture of the plant, even under conditions of low leaf water potential [66,67].Wright et al. [55] describe that C. pyramidale appears to be operating near maximum hydraulic conductivity.The highest values of P N and g s from the first hours of the day agree with Dombroski et al. [37].Also, according to Roach et al. [68], drought conditions promote heat stress at the plant level.This is primarily due to the decrease in E, which, in turn, leads to a reduction in transpiration cooling, as demonstrated by a negative correlation between rain and T leaf (r = −0.902)and T leaf and E (r = −0.744)(Supplementary data file).Higher VPD conditions were observed to coincide with a significant increase in transpiration rates (r = 0.678) and transpiration reduce T leaf (r = −0.744).Smaller leaves come with distinct advantages, including large convection coefficients and reduced resistance to heat transfer across leaf boundary layers [69].This structural adjustment serves as a response to fine-tune leaf temperature in the face of drought stress.This hypothesis is corroborated by a positive correlation between leaf area and transpiration (r = 0.443), T leaf (r = 0.637), and VPD (r =0.474) or negative with WUEi (r = −0.436)maximizing the water use efficiency with less water loss to achieve the same leaf cooling.
Many scholars have shown that trees can increase their WUE when exposed to a water deficit [20,70,71].In this study, we show a reduction in P N due to low Ψ w(pd) (r = 0.788) as well Ψ w(md) (r = 0.851) or g s (r = 0.860).This phenomenon is notably evident through a strong correlation between rainfall and water potential (Ψ w ), with intermediate coefficients observed both in predawn (r = 0.518) and midday (r = 0.528), aligning with the findings of several other studies during the dry season in the Caatinga ecosystem, like Caesalpinia pyramidale [37], Cenostigma pyramidale [36], Spondias tuberosa [20].These findings underscore the remarkable resilience of these plants, demonstrating their capacity to endure and recover swiftly under severe water scarcity conditions.This underscores the remarkable adaptability of this particular species within arid environments.The very negative values of osmotic potential, reaching −7.7 MPa in plants evaluated at midday in August 2019, are consistent with other studies carried out with C. pyramidale, one on plants in a system similar to this one [36] and another with an artificial system in potted plants [72].Only plant species that are fully adapted to semiarid conditions can withstand very negative pressures [36].To do so, the walls of the plant vessels must be highly lignified in a way that avoids collapse or cavitation [69,73].Under water stress, the volume of water in the soil is not enough, so the water potential of the absorbing root cells can extract it from the soil and raise this water to the conductive vessels, which reduces water uptake from the soil [74] as well as its transmission to the leaves, a process that is called hydraulic conductivity (Kp) [75][76][77].An extracted relationship between water deficiency and tree physiology is based on hydraulic failure when the transpiration rate is higher than the root uptake, creating high tension in xylem vessels and leading to embolism and loss of conductivity of the tree transport system [78,79].This perception was also verified in our study when we correlated the rainfall with Kp (r = 0.535), Kp with P N (r = 0.389), g s (r = 0.364), E (r = 0.444), T leaf (r = 0.507), and VPD (r = 0.386) (Supplementary data file).
Plants are extraordinary organisms that have developed an adaptation modifying their Kp, which essentially governs the flow of water through the plant's vascular system [80].In response to limited water availability, plants finely tune their Kp to uphold a favorable water balance [74,81,82].Through this regulation of water flow, plants can efficiently distribute it, ensuring that vital processes like photosynthesis and nutrient uptake continue to function even during arid spells.Also, the Kp shows a significative correlation between rainfall (r = 0.635), leaf area (r = 0.320), Ψ w(pd) (r = 0.333), P N (r =0.389), E (r = 0.444), g s (r = 0.364), VPD (r = 0.386), and T leaf (r = 0.507).These results show an intimate relationship between evaporative demand (VPD), with transpiration and stomatal conductance in support of the current transpiration, assisted by the osmotic potential of the leaves, facts that translate into adjustments in hydraulic conductivity, which resulted in values 1.25-and 1.57-fold higher in the rainy season in comparison to the driest season.If we compare the average Kp of the rainy year and the dry year, the rainy year showed a 2.4-fold higher Kp compared to the dry year.These data are corroborated by a recent study that showed how hydraulic force is formed in Jatropha curcas plants evaluated over a period of 2 years in an arid climate under irrigation with saline water [27].Studies with Cyprus species [83] have shown that plants acclimatized in a given environment have a water memory to the point of regulating their Kp according to the need and the time of exposure to stressful conditions.In that study, the authors showed that in the first year, little significant effect of the moderate drought treatment on hydraulic conductivity was demonstrated; however, there were significant effects in the second year of the study.The same results were shown in Pinus sylvestris and Quercus pubescens [41].However, the reduction in the lumen vessel area emerges as a crucial strategy for countering and evading stresses [45,66,82].Furthermore, other pivotal factors exert their influence on altering resistances within the xylem, consequently leading to variations in hydraulic conductance.These factors encompass the characteristics of secondary wall thickening, perforation plates, vessel dimensions, and density [84].In a study conducted by Xu et al. [84], the intricate relationship between anatomical traits of vessels and water transport within J. curcas xylem was meticulously analyzed.Their findings underscored that despite xylem vessels offering a low-resistance conduit for water transport, even subtle changes in the inner diameter of the vessels had a significant impact on overall resistance.Recently, Pompelli et al. [85] showed by different methodologies that hydraulic conductance measured in the sugarcane stalk showed a fine-tuning with g s as a determining coefficient higher than 0.86 to non-irrigated plants and 0.94 with those that are irrigated.Typically, the atmosphere presents a negative water potential, often as low as −100 MPa.Consequently, leaves are forced to release water into the atmosphere through the process of transpiration [40].However, in order to prevent dehydration and potential harm, when faced with water scarcity, the leaf takes a defensive stance by closing its stomata, effectively minimizing or nearly halting transpiration, including the latent heat fluxes.Consequently, leaves with closed stomata, which are not actively transpiring, experience a rise in T leaf .This increase in leaf temperature is inversely related to photosynthesis, often leading to reduced photosynthetic activity, as shown in this study, where these features showed a negative correlation (r = −0.692) between them.On the other hand, plants inhabiting arid and xerophytic habitats tend to exhibit characteristics that promote the formation of smaller, thicker leaves [22,36,40,86,87].Larger leaves tend to possess reduced convection coefficients and exhibit greater resistance to heat transfer through the boundary layers of the leaves in comparison to their smaller counterparts [69].Consequently, these plants may undergo adjustments in leaf size as an adaptive strategy to optimize their T leaf .As emphasized by Glover [88], trichomes are prominent in saving water.Specifically, the non-glandular trichomes, such as those described herein C. pyramidale.
Oliveira and Almeida [89] demonstrated that leaves of C. pyramidale, which developed in the dry season, were more tender and sclerophyllous than leaves collected in the rainy season.In this study, we describe that in the rainy year, the less rainfall caused an increase in TLT, ADAE, ABAE, and PP, respectively in 24.9%, 19.7%, 27.9%, and 70.8% aligned with Taratima et al. [90] and Taratima et al. [91] which described that anatomical adaptations also do occur under drought stress with an increase in the cuticle, epidermis, and leaf thickness prevent water loss.An increase of 70.8% in PP could be related to an increase in the photosynthetic cells and then a higher number of chloroplasts, as well as a decrease in the thickness of spongy tissue, which facilitates CO 2 reaching chloroplasts in the PP.These anatomical alterations could be an adaptation strategy to facilitate the photosynthesis process under stress conditions [92,93].In fact, in this study, we described a negative correlation (r = −0.391) between the palisade and the spongy parenchyma (Supplementary data file).In accordance with Vogelman et al. [94], the arrangement of a palisade mesophyll composed of a single cell layer arrangement may be directly associated with the channeling of the light path within the leaf.
In accordance with Evert [54], a greater palisade parenchyma should be able to compensate for the reduction in leaf area, thus contributing to the increase in the total leaf thickness as a response to drought stress.As previously reported in the literature [54], xeric leaves displayed a notable increase in xylem thickness and total xylem area when contrasted with their counterparts in the forest.This adaptation appears to be a compensatory response aimed at accommodating the augmented transpiration rates associated with an elevated water supply.In this sense, it is probable that during the dry season, leaves exhibited a decreased xylem thickness in comparison to the wet season, aligning with findings from other investigations [95].These observations strongly suggest that, in these plant species, the hydraulic transport system is predominantly modulated by the availability of water.
Notably, all the attributes associated with the midrib displayed a positive correlation with net photosynthesis in our study.This finding suggests that alterations in vein xylem characteristics may play a pivotal role in regulating hydraulic mechanisms for gas exchange, particularly when adapting to water scarcity.It implies that the species under investigation adapt their midrib features, including the reduction in vessel area and vessel density, as a strategy to enhance their tolerance to reduced water availability during the dry season.These shifts in vein xylem conductivities can significantly influence leaf gas exchange and its responsiveness to changing environmental conditions.This assertion is substantiated by previously published correlation analyses between leaf hydraulic conductance and leaf anatomy [96], which describe the connection between leaf hydraulic conductance and leaf anatomy across six ecologically diverse species under both high and low irradiance levels.Their findings established that anatomical modifications statistically accounted for 40% of the observed variation in leaf hydraulic conductance across species.In line with these previous studies, our results highlighted substantial variations in midrib traits, particularly in the total midrib area, total xylem area, and vessel area.This aligns with findings from other investigations [96,97], further emphasizing the significance of these adaptations in leaf hydraulic regulation.
The ability of plants to alter their characteristics in response to changing environmental conditions is termed phenotypic plasticity, and it appears as morphological and/or physiological changes [98,99].Clearly, the variations in functional attributes depend on the plant species, the choice of the attributes to be analyzed, and the environment to which the plants are subjected [51].Higher plasticity index was previously reported in C. pyramidale [51,57].However, these authors do not conduct any phenotypic plasticity analysis to confirm their hypothesis.For clarity, phenotypic plasticity is an index ranging from 0 to 1 estimated using the maximum and minimum value for each analyzed feature estimation [47,100].Contrary to this, in our study, the phenotypic plasticity was computed for many features in several combinations of conditions.In accordance with Ribeiro et al. [32] and Gratani [98], higher standard deviations in the leaves of trees growing in some xerophytic environments could be related to the phenotypic plasticity of this species.Our results corroborate these presented data and add that due to a higher variance for data, any statistical difference was able to be distinguished in anatomical features, although physiological plasticity was higher than anatomical plasticity.

Conclusions
The data presented in this study leads us to believe that the main effect that led to the extremely low gas exchange in 2019 was the occurrence of the La Niña meteorological phenomenon, which cools the waters of the North Pacific, causing an extreme lack of rain in Northeast Brazil and other areas of the world.The primary action employed by C. pyramidale to mitigate desiccation and sustain essential photosynthetic rates involved lowering its foliar water potential.This reduction, aligned with transpiration chances and stomatal closure, serves to curtail excessive water loss during periods of limited water supply, safeguarding the plant's hydraulic structure even when the leaf water potential is low.Also, our dataset shows that C. pyramidale responds quickly to changes in the environment, such as variations in annual rainfall, VPD, transpiration, SLA, and leaf area.The lower Ψ w reduced the plant metabolism and, per analogy, net photosynthesis as a strategy to stress fugue or stress avoidance because the activation energy is lowered more than Ψ w .Thus, C. pyramidale can be considered an important species in the management and reforestation of degraded areas due to its phenotypic plasticity.This adjustment shows that C. pyramidale is a good indicator for the regrowth of Caatinga dry forest.

Figure 1 .
Figure 1.Brazilian ecosystems highlight the Brazilian Northeast, where the Caatinga ecosystem is predominant.The asterisk (*) denotes the Pernambuco state, which is represented in the base of the figure to show Serra Talhada city, where the experiments were conducted.

Figure 2 .
Figure 2. Rainfall (columns), relative humidity (blue line), mean temperature (green line), and global radiation (yellow line) were measured along the experiment in a restructuring Brazilian Atlantic Forest fragment.Source: INMET [43], free access.

Figure 1 .
Figure 1.Brazilian ecosystems highlight the Brazilian Northeast, where the Caatinga ecosystem is predominant.The asterisk (*) denotes the Pernambuco state, which is represented in the base of the figure to show Serra Talhada city, where the experiments were conducted.

Figure 1 .
Figure 1.Brazilian ecosystems highlight the Brazilian Northeast, where the Caatinga ecosys predominant.The asterisk (*) denotes the Pernambuco state, which is represented in the base figure to show Serra Talhada city, where the experiments were conducted.

Figure 2 .
Figure 2. Rainfall (columns), relative humidity (blue line), mean temperature (green line), and radiation (yellow line) were measured along the experiment in a restructuring Brazilian A Forest fragment.Source: INMET [43], free access.

Figure 2 .
Figure 2. Rainfall (columns), relative humidity (blue line), mean temperature (green line), and global radiation (yellow line) were measured along the experiment in a restructuring Brazilian Atlantic Forest fragment.Source: INMET [43], free access.

Figure 3 .
Figure 3. Water potential measured in Cenostigma pyramidale grown in field conditions in the driest year (2019) and a normally rain year (2020) in both seasonal conditions (rainy, March and dry, August).The measurements were made at predawn (5 a.m.) or at midday (12 p.m.).Means followed by small case letters denote statistical significance between seasons at the same time (5 a.m. or 12 p.m.), means followed by uppercase letters denote statistical significance between time and the same season (rainy and dry months) and the same year, and asterisks (*) denotes statistical difference between years in the same months and time.All values denote media ± SD. n = 10.

Figure 3 .
Figure 3. Water potential measured in Cenostigma pyramidale grown in field conditions in the driest year (2019) and a normally rain year (2020) in both seasonal conditions (rainy, March and dry, August).The measurements were made at predawn (5 a.m.) or at midday (12 p.m.).Means followed by small case letters denote statistical significance between seasons at the same time (5 a.m. or 12 p.m.), means followed by uppercase letters denote statistical significance between time and the same season (rainy and dry months) and the same year, and asterisks (*) denotes statistical difference between years in the same months and time.All values denote media ± SD. n = 10.

Forests 2024 ,
15,  x FOR PEER REVIEW 9 of 29 measured 7-and 6.4-fold higher gs than in 2019, while E registered 26.6-and 22.0-fold higher in 2020 when compared to the 2019 values (Figure4).

Figure 4 .
Figure 4. Net photosynthesis (A-D), stomatal conductance (E-H), and transpiration (I-L) were measured in March 2019 (A,E,I), August 2019 (B,F,J), March 2020 (C,G,K), and August 2020 (D,H,L) in Cenostigma pyramidale grown in field conditions in the driest year (2019) and a normal rain year (2020) in both seasonal condition (rainy, March and dry, August).Means followed by small case letters denote statistical significance between months, regardless of year, and means followed by uppercase letters denote statistical significance between years in the same month.All values denote media ± SD. n = 10.

Figure 4 .
Figure 4. Net photosynthesis (A-D), stomatal conductance (E-H), and transpiration (I-L) were measured in March 2019 (A,E,I), August 2019 (B,F,J), March 2020 (C,G,K), and August 2020 (D,H,L) in Cenostigma pyramidale grown in field conditions in the driest year (2019) and a normal rain year (2020) in both seasonal condition (rainy, March and dry, August).Means followed by small case letters denote statistical significance between months, regardless of year, and means followed by uppercase letters denote statistical significance between years in the same month.All values denote media ± SD. n = 10.

Forests 2024 , 29 Figure 5 .
Figure 5. Daily integrated CO2 assimilation (integrated PN) measured in Cenostigma pyramidale grown in field conditions in the driest year (2019) and a normally rain year (2020) in both seasonal conditions (rainy, March, and dry, August).Means followed by small case letters denote statistical significance between months, regardless of year, and means followed by uppercase letters denote statistical significance between years in the same month.All values denote media ± SD. n = 10.

Figure 5 .
Figure 5. Daily integrated CO 2 assimilation (integrated P N ) measured in Cenostigma pyramidale grown in field conditions in the driest year (2019) and a normally rain year (2020) in both seasonal conditions (rainy, March, and dry, August).Means followed by small case letters denote statistical significance between months, regardless of year, and means followed by uppercase letters denote statistical significance between years in the same month.All values denote media ± SD. n = 10.

Figure 7 .
Figure 7. Internal-to-ambient CO2 concentration (Ci:Ca ratio; A-D), intrinsic water use efficiency (WUEi; E-H), and leaf temperature (Tleaf; I-L) measured in March 2019 (A,E,I), August 2019 (B,F,J), March 2020 (C,G,K), and August 2020 (D,H,L) in Cenostigma pyramidale grown in field conditions in the driest year (2019) and a normal rain year (2020) in both seasonal condition (rainy, March and dry, August).Means followed by small case letters denote statistical significance between months, regardless of year, and means followed by uppercase letters denote statistical significance between years in the same month.All values denote media ± SD.

Figure 7 .
Figure 7. Internal-to-ambient CO 2 concentration (Ci:Ca ratio; (A-D)), intrinsic water use efficiency (WUEi; (E-H)), and leaf temperature (T leaf ; (I-L)) measured in March 2019 (A,E,I), August 2019 (B,F,J), March 2020 (C,G,K), and August 2020 (D,H,L) in Cenostigma pyramidale grown in field conditions in the driest year (2019) and a normal rain year (2020) in both seasonal condition (rainy, March and dry, August).Means followed by small case letters denote statistical significance between months, regardless of year, and means followed by uppercase letters denote statistical significance between years in the same month.All values denote media ± SD.

Forests 2024 , 29 Figure 8 .
Figure 8. Vapor pressure deficit (VPD) measured in March 2019 (A), August 2019 (B), March 2020 (C), and August 2020 (D) in Cenostigma pyramidale grown in field conditions in the driest year (2019) and a normal rain year (2020) in both seasonal condition (rainy, March and dry, August).Means followed by small case letters denote statistical significance between months, regardless of year, and means followed by uppercase letters denote statistical significance between years in the same month.All values denote media ± SD.

Figure 8 .
Figure 8. Vapor pressure deficit (VPD) measured in March 2019 (A), August 2019 (B), March 2020 (C), and August 2020 (D) in Cenostigma pyramidale grown in field conditions in the driest year (2019) and a normal rain year (2020) in both seasonal condition (rainy, March and dry, August).Means followed by small case letters denote statistical significance between months, regardless of year, and means followed by uppercase letters denote statistical significance between years in the same month.All values denote media ± SD.

Figure 8 .
Figure 8. Vapor pressure deficit (VPD) measured in March 2019 (A), August 2019 (B), March 2020 (C), and August 2020 (D) in Cenostigma pyramidale grown in field conditions in the driest year (2019) and a normal rain year (2020) in both seasonal condition (rainy, March and dry, August).Means followed by small case letters denote statistical significance between months, regardless of year, and means followed by uppercase letters denote statistical significance between years in the same month.All values denote media ± SD.

Forests 2024 , 29 Figure 11 .
Figure11.Principal component analysis showing the net photosynthesis (PN) and its correlated features.On the right side are presented all features responsive to PN, while on the left side, all features responsive to water deficit.Each feature is demonstrated with its vector and its strength via different line thicknesses.Ci, substomatal CO2 concentration; w (md), water potential measured at midday; w (pd), water potential measured at predawn; gs, stomatal conductance; E, transpiration; MA, midrib area; NV, number of vessels; Kp, potential specific stem conductivity; VPD, vapor pressure deficit; Tleaf, leaf temperature; MXA, midrib xylem area; VA, vessel area; MXT, midrib xylem thickness; LA, leaf area; ML, midrib length; SP, spongy mesophyll thickness; SLA, specific leaf area; MT, midrib thickness; ADAE, adaxial epidermis surface thickness; TLT, total leaf thickness; ABAE, abaxial epidermis surface thickness.

Figure 11 .
Figure11.Principal component analysis showing the net photosynthesis (P N ) and its correlated features.On the right side are presented all features responsive to P N , while on the left side, all features responsive to water deficit.Each feature is demonstrated with its vector and its strength via different line thicknesses.Ci, substomatal CO 2 concentration; Ψ w(md) , water potential measured at midday; Ψ w(pd) , water potential measured at predawn; g s , stomatal conductance; E, transpiration; MA, midrib area; NV, number of vessels; Kp, potential specific stem conductivity; VPD, vapor pressure deficit; T leaf , leaf temperature; MXA, midrib xylem area; VA, vessel area; MXT, midrib xylem thickness; LA, leaf area; ML, midrib length; SP, spongy mesophyll thickness; SLA, specific leaf area; MT, midrib thickness; ADAE, adaxial epidermis surface thickness; TLT, total leaf thickness; ABAE, abaxial epidermis surface thickness.

Forests 2024 , 29 Figure 12 .
Figure 12.Heatmap showing the degree of modulation in leaf characteristics of Cenostigma pyramidale evaluated as rainfall gradient.The sequence of reading is March 2020, March 2019, August 2019, and August 2020.The color code is followed by scale, and rainfall is scaled too.Asterisks (*) denote significant at p < 0.05.

Figure 12 .
Figure 12.Heatmap showing the degree of modulation in leaf characteristics of Cenostigma pyramidale evaluated as rainfall gradient.The sequence of reading is March 2020, March 2019, August 2019, and August 2020.The color code is followed by scale, and rainfall is scaled too.Asterisks (*) denote significant at p < 0.05.

Table 1 .
Summary of ANOVA for specific leaf area (SLA).Data were determined in two consecutive years (2019 and 2020) and two different months (March and August) in accordance with rainfall gradient.Degrees of freedom residuals (DF), sum of squares (SS), mean squares (MS), F value (F), and p value (p).

Table 2 .
Leaf area and specific leaf area (SLA) were determined in two consecutive years (2019 and 2020) and two different months (March and August) in accordance with rainfall gradient, as the ANOVA table showed in Table mmol H 2 O g −1 s −1 (2019 March, 6 a.m.) to 90 mmol H 2 O g −1 s −1 (2019 March, 6 p.m.).Measurements taken in August 2019, March 2020 and August 2020 show that g s ranged respectively from 95.6 mmol H 2

Table 3 .
Total leaf thickness, adaxial epidermis thickness, abaxial epidermis thickness, palisade mesophyll thickness, and spongy mesophyll thickness measured in Cenostigma pyramidale grown in field conditions in the driest year (2019) and a normal rain year (2020) in both seasonal conditions (rainy, March and dry, August).Means followed by small case letters denote statistical significance between months, regardless of year, and means followed by uppercase letters denote statistical significance between years in the same month.All values denote media ± SD. n = 10.

Table 4 .
Midrib thickness, midrib length, midrib area, midrib xylem thickness, midrib xylem area, vessel area, number of vessels, and potential specific stem conductivity measured in Cenostigma pyramidale grown in field conditions in the driest year (