A comparative study on the leaf anatomical structure of Camellia oleifera in a low-hot valley area in Guizhou Province, China

The leaf serves as an important assimilation organ of plants, and the anatomical structure of leaves can reflect the adaptability of the plant to the environment to a certain extent. The current study aimed to cultivate superior local cultivars, and 35 healthy individual plants were selected from the Camellia oleifera germplasm resource nursery for a comparative study of the leaf structure. In July 2019, the leaves were collected from 35 selected healthy C. oleifera plants, and the leaf structure was observed by using the paraffin section method. Healthy individual plants were screened using variance analysis, correlation analysis and cluster analysis. The representative indices were selected according to the cluster membership, correlation indices and coefficient of variation (C/V) for a comprehensive evaluation of drought resistance via the membership function. There were extremely significant differences in 11 indices of leaf structure for these 35 healthy plants. C18 had the greatest leaf thickness, C7 the largest spongy tissue, and C38 the largest ratio of palisade tissue thickness to spongy tissue thickness (P/S). The clustering results of the healthy individual plants differed significantly. The membership function showed that the drought resistance of 35 C. oleifera plants was divided into five categories. C18 had very strong drought resistance, and C3, C7 and C40 had strong drought resistance. There were significant differences in terms of the upper epidermis, P/S ratio and spongy tissue among the C. oleifera plants. C18, C3, C7 and C40 exhibited satisfactory drought resistance. Although C39 and C26 had moderate drought resistance, their P/S ratios were high, which might be used to cultivate high-yield and drought-resistant C. oleifera varieties. The leaf P/S ratio of C. oleifera from low-hot valley areas was high. Among various leaf structures, spongy tissue, upper epidermis, P/S ratio and cuticle constitute the drought resistance evaluation indices for C. oleifera grown in low-hot valley areas.


Introduction
The leaf is an assimilation organ of the plant whose photosynthesis involved in its structural response to drought stress is complex, involves interactions between different structural levels [1,2], and directly affects photosynthetic factors such as water retention in the leaf, CO 2 stomatal conductance, and mesophyll conductance [3,4]. Water and CO 2 are substrates for photosynthesis. Under drought conditions, stomatal closure is an early adaptation of plants to cope with water deficit but also limits CO 2 uptake by leaves [5]. At this time, how other tissue structures maintain photosynthesis under drought conditions is of crucial importance for plants [6].
Under drought stress, to maintain its biological function to ensure the normal operation of photosynthesis, the leaf is often required to change its anatomical structure to adapt to the environment. Therefore, the same species may evolve different structural characteristics when faced with varying degrees of drought stress. There are differences among individuals within the population. Additionally, facing the same degree of drought stress, the same species may also exhibit different structural characteristics [7,8], which is also referred to as leaf plasticity; the plasticity of leaf anatomical structure and physiology is an important guarantee for plants to adapt to adverse environment [9]. Among the various anatomical characteristics, the ratio between palisade tissue and spongy tissue (P/S) is positively correlated with the net photosynthetic rate in a significant (r = 0.9) or a very significant (r = 0.985) manner [10,53]. Leaf anatomical structure reflects important photosynthetic physiological characteristics and is closely associated with function, and therefore, leaf structure is widely used to evaluate the drought resistance of different varieties of the same species [11][12][13].
Guizhou Province is the center of karst topography in Southwest China. The soil in karst areas is rich in calcium and characterized by severe erosion. The direct leakage of water content in soil underground is the main form of soil erosion. The growth and development of plants in karst areas are often accompanied by environmental conditions featuring water deficits and high calcium [14][15][16][17]. In Guizhou Province, soil erosion is mainly distributed in Wangmo County and Ceheng County in Southwest Guizhou [18]. In addition, due to gradual reductions in annual precipitation in recent years [19], droughts occur frequently as a result. The Beipan River flowing through Guizhou Province forms a low-hot valley area, where a distribution of Camellia oleifera is found. However, related research has not been conducted in this area. However, the seeds of Quercus sichourensis distributed in this area show higher drought resistance than those in other distribution areas [20]. In Hunan Province, the distribution center of Camellia oleifera, the annual precipitation is approximately 1450 mm, whereas that in low-hot valley areas is approximately 1200 mm; in the meantime, for karst topography in low-hot valley areas, water loss is more likely to occur. Therefore, C. oleifera in the low-hot valley area may have the potential for undeveloped germplasm resources with desirable drought resistance.
Camellia oleifera is widely distributed in southern China, but not all distribution areas can be called the most suitable cultivation area. One of the possible reasons for this is that C. oleifera is a calciphobous plant [21]. Therefore, the National Forestry and Grassland Administration issued a development plan for the C. oleifera industry in China and classified areas suitable for cultivation of C. oleifera according to its distribution areas into the most suitable cultivation area, suitable cultivation area and relatively suitable area (arranged in descending order). The habitat of the last two areas is worse than that of the most suitable one. As a result, C. oleifera distributed in the last two areas is prone to more abiotic stress, which also makes the local C. oleifera cultivars better adapted to these distribution areas. Therefore, it is one of our tasks to select local C. oleifera cultivars [22]. For this purpose, we selected 35 healthy individual plants in this study (to be published) as materials for seed breeding according to the yield criterion.
At present, the seed breeding of C. oleifera in the low-hot valley area is mainly to select healthy individual plants. However, currently, the adaptative mechanism of C. oleifera in the low-hot valley area is still poorly understood because the leaf structure of the plant reflects the plant's adaptation to the environment to a certain extent. Therefore, we conducted a comparative study to investigate the leaf anatomical structure of individual Camellia oleifera plants grown in the low-hot valley area in Guizhou Province, China. The following issues would be clarified by studying the leaf structure of C. oleifera. (1) What are the characteristics of anatomical structure of C. oleifera leaves in the low-hot valley area, and are there any differences between the anatomical structure of its leaves and those in other cultivation areas reported in the literature? (2) How does the leaf structure of C. oleifera in the low-hot valley adapt to water deficit and high calcium? (3) Are there any differences in drought resistance in terms of leaf anatomical structure for the 35 selected healthy C. oleifera plants? Our research is helpful in explaining how the leaf structure of C. oleifera, a calciphobous plant, adapts to drought and high-calcium environments in karst areas and will provide a theoretical basis for the selection of C. oleifera breeding materials in the low-hot valley area.

Sample collection
Leaf samples of healthy individual plants were collected in 2019 at the Guizhou University Camellia Oleifera Research Station, located in Ceheng County in southwestern Guizhou Province, China (24.71˚-24.94˚N, 105.79˚-106.05˚E), which has an average annual temperature of 19.7˚C, an average temperature of 27.2˚C in the hottest month and 10.5˚C in the coldest month, an annual accumulation temperature >10˚C of 6,348˚C, an annual precipitation of 1,197 mm, heat and precipitation in the same seasons, precipitation from November to April accounting for only 16% of the annual precipitation, a tendency of rising temperature, little rainfall and frequent droughts in March and April, a rainy season starting from May, and an annual sun exposure of 1,257 h. Biyou Town and Yangba Town are the two primary distribution areas of C. Oleifera in Ceheng County. At these two towns, the planting area reached up to 5786.67 hm 2 , which accounted for 53.90% of that of the whole county. These two towns are also the main excellent C. Oleifera seed supply sources in the low-hot valley area in Guizhou Province. In this study, the study area was mainly centered on Biyou Town. In the study area, the soil was the typical yellow soil in karst regions of Southwest China. Severe soil erosion leads to a shallow soil layer, and the resulting drought problem is a great challenge faced by C. oleifera planting in this area ( Fig 1A).
The nursery was established as follows. In the study area, C. oleifera plants with satisfactory economic characteristics were recorded (S2 Table). Their seeds were collected and then sowed in 2007 and 2008. In the sowing year, the land was prepared before sowing. To guarantee no less than one plant from the same preliminarily selected individual survived to blossom and bear fruit (normally eight years), three repetitions were set. The planting density was 2000 plants per hectare, which was slightly higher than the recommended density in the literature (900-1,800 plants per hectare) [23]. Due to differences in individual development, the C. oleifera plants also differed somewhat in canopy height, which led to inconsistent sunlight exposure among them or, even worse, disadvantageous tree death during the long-term photosynthesis-biomass accumulation process [24]. This phenomenon indicates that differences in the competitiveness to sunlight exist among the individuals.
At the time of sample collection, the trees aged 12-13 years, and they all reached fruit productive age [23]. The C. oleifera resource nursery is on the south slope. The soil is yellow soil with the soil layer thickness of 51.6±8.2 cm (Fig 1B). Between 2013 and 2018, the economic characteristics of the individual plants at the nursery were consecutively determined, and the plants growing at the stand edge were excluded. Within the involved six years, yield was observed within the late three continuous years. All these individuals were planted within two years, and their flowering and fruiting time had a difference of 1-2 years. This long-term determination was mainly based on the considerations that the growth differences among the trees led to the differences in the flowering and fruit setting periods and that long-term observation could help obtain more stable information about the trees.
We selected 45 healthy Camellia oleifera plants in total. The individual plants were named (C-number) according to the sequence of the healthy individual plants identified. Their yields per unit crown width were recorded (S1 Table).
For the current purpose, we excluded individuals with unstable yields and/or severe insect and pest diseases (original number was used). Insect and pest diseases mainly occurred between 2010 and 2011. The primary reason was that drought persisted in Southwest China in 2019 and 2010 [25], which weakened the resistance of young C. oleifera trees in Ceheng county [26]. Moreover, the high temperature accompanied by drought might bring the first appearance and peak of pests to earlier dates, which had adverse effects on C. oleifera [27]. Additionally, C. oleifera is plant of flowering-fruiting synchronization, and it normally takes 18-19 months from bud differentiation to fruit maturation. The number of differentiated flower buds in one year (it is the basis of fruit yield in the second year) is negatively correlated with fruit development in the same year, which leads to instable yields in successive years [22]. Although instable yields can be remedied though additional cultivation measures, such as artificial flower and fruit thinning and water and fertilizer management, these measures increase the production cost. Finally, 35 healthy plants were selected from the resource nursery for a comparative study of leaf structure.
There is an overlapping period between flower bud differentiation and fruit growth of C. oleifera, particularly in July and August each year. This period is the key time of oil conversion in C. oleifera fruit as well as the time of the C. oleifera leaf with the highest photosynthetic efficiency; the photosynthetic products of the leaves close to the base of the current-year branch are mainly provided for fruit development, whereas those with the top leaves are mainly provided for flower buds [28]. Considering that the competitiveness of the excellent plants differed somewhat, the leaves directly exposed to sunlight were taken as the samples. Therefore, in July 2019, three branches in the upper and middle parts of the east-facing canopy were selected from 35 healthy plants, and three leaves, the 5th to 7th healthy and mature leaves in the C. oleifera phyllotaxis, were selected from each branch. A total of nine leaf replicates were selected from each plant. Then, 0.5 � 0.5 mm leaf tissue was removed from the middle of each leaf while avoiding the midrib, which was immediately placed into carnoys fixative (95% ethanol: glacial acetic acid (V/V) = 3:1) for a fixation of 6 h. The leaves of the same plant were fixed in the same bottle. The bottle was exhausted of air after the fixation was completed, transferred to 70% alcohol solution and stored in a 4˚C refrigerator.

Precipitation in the study area
The percentage of precipitation anomaly (Pa) can visually reflect the drought condition confronted by plants in specific areas, and Pa is calculated as follows [29]: where � Prepresents the average monthly precipitation and P is the precipitation in a specific month. Pa > 0.4 indicates drought stress. The data used for precipitation analysis in this study were from the China Meteorological Information Center (http://data.cma.cn/).

Observation of leaf anatomical structure
The fixed material was dehydrated in ethanol, made transparent in xylene, embedded in paraffin, and then sliced using a microtome (Leica RM 2235, Germany) with a section thickness of 8 μm. Safranine O-fast green staining was performed [30] after the specimen was dried. Then, it was sealed using neutral resin and photographed with an optical microscope (Leica DM 3000, Germany). Data were measured by ImageJ (avoiding leaf veins and heterolayer epidermis). Three measurements were performed for each leaf, and an average was obtained. The measured parameters included leaf thickness (LT), cuticle thickness, thickness of upper epidermis (TU), thickness of palisade tissue (TP), thickness of spongy tissue (TS), thickness of lower epidermis (TL), and thickness of first layer palisade tissue (TFP). The organizational tightness (CTR), tissue porosity (SR), upper and lower epidermis thickness ratio (U/L), and the P/S ratio were calculated. CTR = CTR = TP/LT � 100; SR = SR = TS/LT � 100; U/L ratio = upper epidermis thickness to lower epidermis thickness; P/S ratio = TP/ TS. Coefficient of variation (C/V) = standard deviation/mean ×100%, the plasticity index (PI) of leaf tissue structures = the minimum mean of an index of an individual plant/the maximum mean of the index of the individual plant [31].

Data processing
Excel 2010 was used to sort out measurement results, Photoshop cs6 for plotting, SPSS19.0 for variance analysis, and Duncan for multiple comparison and hierarchical clustering of the data. The leaf tissue structure of an individual plant is indicated by the mean±SD. The rest is described by the mean. The mean leaf tissue of an individual plant was selected for clustering analysis.
According to the clustering of organizational structure, the typical indices were calculated by the correlation index [32], whose formula is: where R 2 i is the correlation index for each index in each category, n is the number of indices in each category, and r is the correlation coefficient between one index and other indices in the same category. When there is only one index in the same category, the value of the correlation index is 1. A larger correlation index between a certain index and other indices in the same category better represents the characteristics of this category. With reference to the correlation index and C/V for index selection, the drought resistance of leaves of 35 C. oleifera plants was evaluated in terms of the leaf anatomical structure by using the membership function in fuzzy mathematics [33]: is the value of the drought resistance membership function; x i is the mean leaf tissue structure of individual plants; and x min and x max are the average minimum and maximum of the values measured for the leaf tissue structure of all individual plants, respectively. If a certain index is negatively correlated with drought resistance, then the inverse membership function is used for conversion. The calculation formula is: After calculating the average membership degree of the indices, the healthy individual plants with a higher average degree of membership have a better stress resistance.

Precipitation at the leaf development stage of C. oleifera in Ceheng in 2019
After C. oleifera shooting period, the precipitation in Ceheng County between March and April in 2019 were comparable to those in the same time period between 2000 and 2018. In May, the precipitation (1182 mm) was lower than the normal range (1636 ± 383 mm), with a Pa value of 0.28, although this precipitation did no fall within the drought range on the monthly scale (Fig 2). However, further observation showed that the precipitation in May mainly concentrated in late May (Fig 3), which indicated that precipitation on the monthly scale might underestimate the drought condition in Ceheng in early summer. To keep consistent with the Pa index, we analyzed the precipitation data between April 25 and May 24 (totally, 30 d), and found that the precipitation in Ceheng in early summer of 2019 was 838 mm (the average in the considered years was 1387 mm), with a Pa value of 0.396. The precipitation during that period in 2019 was even lower than that in 2009 and in 2010 (these two years were drought years); that is, Ceheng suffered climate drought during that period. Nevertheless, the drought time lasted for one month, and the precipitation in June returned to normal.

Anatomical structural characteristics of C. oleifera leaves in the low-hot valley
The leaf structure of C. oleifera consisted of cuticle, upper epidermis, palisade tissue, spongy tissue, lower epidermis and vein (Fig 4). Most upper epidermis contained monolayer cells with double-layer cells (multiple epidermis), and some of the second layer of upper epidermis cells contained crystals. The palisade tissue was mainly composed of three layers (Fig 4A), and some leaves had four layers (Fig 4B). The palisade cells in the upper layer were mainly columnar-shaped cells with a dense distribution, while the palisade cells in the lower layer contained short-type cells with a relatively sparse distribution. The palisade cells in each layer were arranged in an orderly manner, and the staining using the safranine O-fast green staining method was relatively uniform. The thickness of spongy tissue was slightly higher than that of palisade tissue, where crystal distribution (P/S ratio = 0.95) was observed. The adaxial surface  of the leaf midrib was covered with thin hair (Fig 4C), while the upper and lower epidermis were hairless, with the upper epidermis thicker than the lower epidermis (U/L = 1.27).

Differences in leaf structure of healthy individual plants
Palisade tissue and spongy tissue were the main components (90.34%) of the C. oleifera leaf structure. These tissues were associated with the growth and development of plants in terms of light energy fixation, and they also contained important indices of plant stress resistance. There were very significant differences in different leaf tissues of healthy C. oleifera plants from the low-hot valley (Tables 1 and 2).
The leaf thickness of C18 was the largest, reaching 642.77 μm. Its palisade tissue and upper epidermis were also the thickest, reaching 297.83 μm and 37.56 μm, respectively. At the same time, the U/L ratio was the greatest, measuring 1.68. The leaf thickness and palisade tissue of C8 were the smallest, measuring 388.76 and 388.76 μm, respectively. The U/L ratio of C2 was the smallest, measuring 1.01. The spongy tissue of C7 was thick, reaching 301.45 μm. The spongy tissue and upper epidermis of C38 were the thinnest, at 166.27 and 20.57 μm, respectively. The porosity of C5 was the largest, while the CTR and P/S ratio were the smallest, which were 38.08, 52.02 and 0.73 μm, respectively. The porosity of C26 was the smallest, while the CTR and P/S ratio were the largest, measuring 38.85, 53.01 and 1.36, respectively. In addition, the thickness of the lower epidermis was the thinnest, at 18 μm ( Fig 5).
LT, palisade tissue, spongy thickness, upper epidermis, TFP and cuticle exhibited extremely significant correlations, while a significant correlation was found between lower epidermis and other tissues except for the cuticle, showing a coevolution of leaf tissue structure to adapt to the environment. There was a significant correlation between porosity and leaf thickness but no correlation between tightness and leaf thickness. There was no correlation between the P/S ratio and cuticle/upper epidermis thickness but a significant correlation between the P/S ratio and other tissues. The U/L ratio was extremely significantly correlated with cuticle thickness and TFP in addition to its extremely significant correlation with the upper and lower epidermis. Additionally, it was significantly correlated with palisade tissue thickness but was not correlated with other issues (Table 3).
According  The P/I of leaf tissue might be one of the important characteristics for C. oleifera to adapt to the low-hot valley habitat.

Evaluation of the leaf anatomical structure of healthy C. oleifera plants
In contrast to animals, forest trees are unable to avoid bad site conditions. The term "environment" can be used to describe the sum of the ambient conditional factors that affect the growth and development of animals and plants, while "site conditions" are specifically used for the environment in which forest trees grow. Plant leaves are assimilating organs, and the structural characteristics of their different internal tissues enable them to adapt to different site conditions. In this study, a hierarchical cluster analysis was performed according to the average values of 11 leaf anatomical indices of 35 C. oleifera plants. The tissue structure indices were classified by Euclidean measurement, in which the cuticle thickness at 1.5 and the thickness of the lower epidermis were classified in the same category, inconsistent with the result that the cuticle and the lower epidermis were not correlated in the correlation analysis. Therefore, cosine was adopted for measurement in the cluster analysis of tissue structure, while Euclidean measurement was used for the cluster analysis of individual plants.
According to the clustering analytical results, 11 indices (a distance of 15) and 35 healthy individual plants (a distance of 15) were clustered into 4 categories, and the distance between these categories was relatively large (Figs 6 and 7). The first category of leaf structure indices included LT, spongy thickness, lower epidermis and SR; the second category included upper epidermis, U/L ratio and TFP; the third category included cuticle; and the fourth category included P/S ratio, CTR and palisade thickness.
The cluster analysis of healthy individual plants revealed that the first category contained the best individual plants, with 23 plants clustering in one category; the second category contained C2, C8, C21, C26, C33, C38, and C39; the third category contained C3, C13 and C22; and the fourth category contained only C7 and C18. The characteristics of each category are shown in Table 4.
According to the clustering results of tissue structure, the item with the greatest correlation index in the same category (Table 5) was selected. If the difference between the two items with the largest correlation index in the same category was not sufficiently large, the selection was made with reference to C/V (Tables 1 and 2). In the fourth category, the difference between the correlation indices of CTR (0.66) and the P/S ratio (0.62) was not sufficiently large, and the C/V of the P/S ratio was the largest among various indices, reaching 22.00%, and CTR was 10.98%. Therefore, the P/S ratio was selected as the evaluation index for the fourth category.
Spongy tissue thickness, upper epidermis thickness, cuticle and P/S ratio were selected to evaluate and order the membership functions (Table 6). According to the ranking of the average membership degree of drought resistance and using the scoring difference of 0.2 to define one subrange (Table 5), the healthy individual plants of C. oleifera in the low-hot valley were divided and ranked into five categories: super strong > strong > medium > weak > none in terms of drought resistance. Only C18 had very strong drought resistance, and C3, C7 and C40 fell in the category of strong drought resistance. Camellia oleifera plants with medium drought resistance accounted for the majority, and included C29, C22, C6, C23, C30, C39, C11, C5, C13, C31, C14, C19, C26, C25, C10, C41, C17, C21, C15, C24, C35, C37 and C12. Plants with weak drought resistance included C36, C9, C4, C33, C27, C38 and C8. Only C2 had no drought resistance. In terms of distribution, the classification results have a normal distribution and therefore can be evaluated for significance (Table 6).

Discussion
Plant leaves are the organs most susceptible to environmental regulation [34]. Their structural characteristics vary with environmental changes. This property is also known as the structural plasticity of leaves [9]. Plasticity is of great significance for immobile plants to adapt to abiotic stresses of different types and varying degrees [35,36]. Tree body mainly provides water for leaves, which receive more solar radiation compared with other plant organs [37]. In this study, we selected the leaves fully exposed to sunlight for investigation. The palisade tissue (0.48) and the spongy tissue (0.45) both displayed high plasticity, and these indices determined leaf thickness (r = 0.724 and 0.819, respectively). Under short-term drought stress, Canarium album and Platanus orientalis maintain moisture by increasing leaf thickness [3,38] while Lycopersicon esculentum decreases leaf thickness under long-term drought stress [39]. Under short-term drought stress, the decrease in the water content of leaves exhibits a certain pattern: The water content of non-functional leaves (lower leaves, morphologically) continuously decreases, whereas that of functional leaves (upper leaves) is maintained to keep biomass accumulation [40]. After rehydration, functional leaves maintain their strategy for drought. Platanus orientalis maintains moisture by increasing leaf thickness while Platanus orientalis exerts antioxidant protection by increasing non-enzymatics [38], which indicate the diversity of different species in dealing with even the same stress. In our study, the excellent individual plants underwent rehydration after short-term drought, and their drought resistance was assessed based on leaf anatomic structure. The plasticity of leaf thickness was 0.4. Presumably, the main reason is that the involved 35 plants were from the karst area where the habitat is complex, and therefore, their strategies to cope with drought stress differ. In the future, the physiological functions of the plants remain to be assessed.
The karst area has a high calcium environment. The Ca content of plants in this area is significantly higher than that in nonkarst areas. The Ca content in the aboveground part of plants is significantly higher than that in the belowground part [41]. Plants can mostly adapt to this site condition. Lonicera confusa leaves in karst areas store or remove excess calcium through glandular hairs and stomata under high Ca conditions [42]. When cultivated under high calcium conditions, a larger area of calcium oxalate crystals was formed beneath the cuticle of Dracaena sanderiana compared with low calcium culture [43]. At present, among the studies on the upper epidermis of C. oleifera leaves, the upper epidermis cells of different cultivars are uniformly composed of one layer of cells [22,44,45], while in the low-hot valley, a heterolayer phenomenon was observed in the upper epidermis cell layers. Calcium oxalate crystals were

PLOS ONE
occasionally visible in the lower epidermis cells. An excessive supply of calcium and water loss contributes to the formation of calcium oxalate [46,47]. Calcium oxalate, closely related to the maintenance of ion balance in plants, is widely present in plants [48,49]. In addition, calcium oxalate crystals have a certain significance in taxonomy [50,51]. In the low-hot valley, a heterolayer phenomenon was observed in the upper epidermis in C. oleifera leaves, and crystals were distributed in the lower epidermis cells, which have not been discovered in C. oleifera leaves in other studies.
When plants are under drought conditions, maintaining photosynthesis is of essential importance for them to survive drought stress. In photosynthesis, CO 2 diffuses into mesophyll cells (palisade tissue) through stomata, passes through multilayer cells and reaches chloroplasts in palisade tissue before being fixed there [52]. Palisade tissue is the main site where chlorophyll is distributed. Plants with short transport distances from stomata to chloroplasts (rather thin spongy tissue and thick palisade tissue-P/S ratio) may consume less energy and have a higher net photosynthetic rate. Therefore, there is a significant positive correlation between the net photosynthetic rate of C. oleifera and the P/S ratio, yet no significant correlation between the net photosynthetic rate and palisade tissue [53,54]. For C. oleifera in the most suitable cultivation area, its P/S ratio ranged from 0.6 to 0.8 [54], lower than that (P/S ratio = 0.94) in healthy C. oleifera plants in the low-hot valley. These results may be responsible for the higher net photosynthetic rate of C. oleifera leaves in the low-hot valley than in other regions, making it easier for C. oleifera in the low-hot valley to survive the drought period.
The cuticle, covering the upper epidermis, serves as the boundary between the plant leaves and the environment, protecting plants from biotic or abiotic stresses, and it is the main barrier to limit nonstomatal water loss, thereby promoting the survival of plants under drought conditions [55][56][57][58]. Arabidopsis thaliana, whose cutin synthesis was impeded, displayed intolerance to drought at its seedling stage, and its leaves showed signs of necrosis [59]. Research showed that the synthesis of the main structural components of the cuticle was induced by drought [60]. A thick cuticle and upper epidermis are also characteristics of plants in arid areas [61]. These descriptions might explain why the cuticle was one drought resistance evaluation factor according to this study.
The carbon starvation hypothesis for plant mortality due to drought [62] holds that under drought conditions, the stomata of plants are closed by hydraulic pressure to prevent water loss, which diminishes the absorption of carbon. However, the respiratory demand caused by high temperature accompanying drought aggravates the metabolic demand for carbon and causes plants to die of carbon starvation. To cope with drought conditions, plants often reduce aboveground biomass and increase belowground biomass to increase the probability of water availability. Therefore, the root system is well-developed. Only with this adaptation, however, CO 2 required for photosynthesis is still insufficient. Through carbon starvation culture, calcium oxalate crystals in Amaranthus hybridus and Colobantus quitensis leaves significantly decreased. The degradation of crystals accompanied by the enhancement in oxalate oxidase activity might be attributable to the release of CO 2 from the decomposition of calcium oxalate crystals to function as the substrate of photosynthesis to maintain photosynthesis [63][64][65]. In relevant research on the structure of C. oleifera leaves, the distribution of crystals was observed only in a few varieties [66], and crystals were only found in the leaf margin and parenchymatous tissue under the main vein [22]. In the C. oleifera leaves from the low-hot valley, crystals widely existed in the spongy tissue, which might be because the soil in the karst area contained a large amount of Ca. The content of Ca in the C. oleifera leaves was the highest among the contents of minerals in the region [67], which was an indication that plants in the karst region adapted to the environment, manifested in the storage of calcium oxalate crystals in C. oleifera from the low-hot valley. Greater leaf thickness is conducive to the storage of calcium oxalate crystals, while the determination coefficient of spongy tissue thickness and leaf was the largest (r2 = 0.671), greater than that of palisade tissue (r 2 = 0.524). On the other hand, the conduction of CO 2 in leaves is divided into stomatal conductance and mesophyll conductance. Drought reduces mesophyll conductance and causes the stomata to diminish in size [68][69][70]. Thicker spongy tissue may increase CO 2 binding sites [3,71], promote the diffusion of gaseous CO 2 , and alleviate the reduction in CO 2 caused by the closure of the stomata. However, in most studies on the effect of drought stress on leaf structure, spongy tissue is inversely proportional to the degree of stress. Contrary to this study, this may be because we used yield as the indicator to select healthy individual plants. The environmental stress on the selected healthy individual plants did not reach the extent to which it reduced the thickness of spongy tissue. The phenomenon of soil erosion is serious in the low-hot valley, which results in the formation of a complex habitat. A single evaluation index cannot reflect the adaptability of C. oleifera to different habitats. Therefore, it is one-sided. Multiple indices jointly adopted for evaluation can effectively evaluate the interaction between the indices. There have been quite a few studies on the evaluation of plant stress resistance using the membership function method. Guo et al. [72] evaluated the cold resistance of different varieties of Prunus persica through the membership function, and the evaluation results were basically consistent with the field survey results. Zou [73] evaluated the drought resistance of Gossypium hirsutum varieties using physiological indices and the membership function and carried out evaluation verification of replanting under drought conditions. The results were consistent with the evaluation by the membership function, which has been widely adopted for the evaluation of plant stress resistance. In this study, we selected the cuticle, upper epidermis, spongy tissue and P/S ratio for a comprehensive evaluation and discussed these indices in the previous section. However, further study is still needed to verify these physiological indices and results.
In this study, we found that the leaf anatomical structure of C. oleifera from the low-hot valley was significantly different from that in other distribution areas. The C. oleifera leaf from the low-hot valley had a higher P/S ratio, and a distinct distribution of calcium oxalate crystals was observed in the mesophyll, which reflects the adaptability of the C. oleifera leaf structure to this area. Through the evaluation of the membership function, we screened four healthy C. oleifera plants with desirable drought resistance, namely, C18, C3, C7 and C40. Although the drought resistance of C39 and C26 was medium, they had a higher P/S ratio, making them promising for cultivating high-yield C. oleifera plants. These findings will provide a theoretical basis for the selection of breeding materials, will promote purposeful breeding of healthy cultivars and will decrease breeding costs.
Supporting information S1