Introduction

Chestnut (Castanea mollissima Blume), belonging to Castanea mollissima, has a long cultivation history and is an economically important fruit tree species with wide distribution, diverse processing products and high comprehensive utilization value1. The fruit industry across the world is facing a major constraint in the form of drought due to climate change. This has led to a significant hindrance in its development2,3. Accurate evaluation of drought tolerance of chestnut is very important for selecting cultivars and improving the yield and quality of chestnut4. Water determines the species and quantity of plants5,6. Changes in the physiological responses of plants to drought tolerance have been extensively studied such as changes in physiological indices, e.g., root morphology, osmotic adjustment and antioxidant enzyme activity7,8,9,10. When plants are subjected to drought stress, leaves are most sensitive to drought. Under drought conditions, water transport is slow or even stops, which hinders photosynthesis and affects plant survival and development. Chlorophyll, leaf relative water content and nutrient content stored are significantly reduced, and finally, the yield is reduced11. Generally, the osmotic substances in leaves significantly increase with the deepening of drought such as the accumulation of soluble sugar, proline (PRO) and betaine. Osmotic regulation can help plants stabilize cell expansion, maintain normal photosynthesis, protect enzyme activity and improve drought tolerance12.

Under drought stress, the antioxidant system in plants also responds to drought stress in a timely manner, improving the activities of related enzymes such as superoxide dismutase (SOD) and peroxidase (POD) to remove harmful substances produced by drought in plants and reducing the impact of stress on the normal life activities of plants13. Abiotic stress likes drought or salinity concurrently causes oxidative stress14 and osmotic stress15 which affect the biochemical and physiological processes of plants, including photosynthesis16, energy metabolism, plant water, nutrient uptake imbalance, protein, etc.17, and ultimately affect the growth and yield of crops18. Osmotic stress in plants changes many enzymatic antioxidants, such as CAT, SOD, GOPX19, non-enzymatic antioxidants, such as carotenoids20, pigments like betalains, and chlorophylls, beta-carotene21, phenolic and flavonoids, etc. having high antiradical capacities22 which can manage the adverse impact of abiotic stresses. Under abiotic stress, the plant itself regulates different pathways to increase the application of these antioxidants to detoxify the ROS and mitigate environmental stress23. Therefore, the changes in physiological and biochemical indices and the antioxidant system of plants under drought stress are important manifestations to measure the differences in drought tolerance of plants. Flavonoids are important secondary metabolites that are distributed in various parts of plants. They are not only auxiliary color substances of plant flowers but also help plants stabilize oxidative homeostasis and resist osmotic stress, and metal ions combine to reduce damage to the cytoplasm to improve plant drought tolerance24,25.

Currently, research on the drought tolerance of chestnut has mainly focused on photosynthesis, membrane permeability, enzyme activity, osmotic substances and other physiological indicators26,27,28. Flavonoids and flavonols, as important antioxidants, are rarely used as drought tolerance indicators for evaluation29. In this paper, flavonoids and flavonols were involved in the ranking of drought tolerance of 6 chestnut varieties (Clones). Using the method of natural drought stress, the changes of 15 physiological and biochemical indexes of 6 chestnut varieties (Clones) under drought stress were studied. According to each index, the drought tolerance of 6 chestnut varieties (Clones) was comprehensively evaluated by principal component analysis and membership function method, which was the theoretical basis for the subsequent drought tolerance research of chestnut.

Results

Determination of the soil water content of different chestnut varieties (Clones) under drought stress

Table 1 showed that when the soil relative water content (RWC) of each test variety (clones) was approximately 30%, the sample leaves were collected. For each variety, the measurement was repeated three times to measure the soil relative water content and soil absolute water content, which were at the level of severe drought. It can be concluded that the absolute soil water content obtained by the aluminum box weighing method was significantly correlated with the relative soil water content measured by the instrument in Table 2.

Table 1 Determination of soil water content under drought stress.
Full size table
Table 2 Correlation between soil water content and soil relative water content.
Full size table

Analysis of enzyme system related indices of different chestnut varieties (clones) under drought stress

SOD scavenges harmful substances in the metabolic activities of fruit trees, scavenges free radicals produced by chestnut under drought stress, and repairs cells. It was observed that there was no significant difference between QX42, YQ, QL45, YSZF and YZ, and there was no significant difference between YSZF, YZ and YL. QX42, YQ, QL45, YSZF and YZ were significantly higher than YL in Table 3. Among them, the SOD activity of QX42 was the highest. Thus, QX42 had strong drought tolerance, followed by YQ and QL45.

Table 3 Analysis of enzyme related indexes of different chestnut varieties ( lines ) under drought stress.
Full size table

POD catalyzes many biochemical reactions in fruit trees and has the effect of eliminating cytotoxicity. Table 3 showed that under drought stress, the POD differences of different chestnut varieties showed that YQ, YZ and YL were significantly higher than QL45 and YSZF, and the latter two varieties were significantly higher than QX42. Under drought stress, the activities of YQ and YZ were higher than those of YL. The POD content of YQ was 2.58 times that of QX42.

CAT can decompose hydrogen peroxide into oxygen and water to prevent excessive damage to the cell membrane of fruit tree cells under drought stress. Table 3 shows that the difference between QX 42 and QL 45 was not significant, QX 42 and QL 45 were significantly higher than YZ, YSZF, YQ and YL and YZ, YZ and YL, QX 42 and QX 42, followed by QL 45, in which YL was the lowest and all other varieties were higher than the control YL.

APX is an important enzyme for fruit trees to resist oxidative stress, which helps remove excessive peroxides when fruit trees are subjected to drought stress. There was no significant difference between QL45, YZ, YQ, YSZF and QX42, and QL45 was significantly higher than YL. The activity of QL45 APX was 68.92 μmol/min/g fresh weight, 1.65 times that of YL, as shown in Table 3.

Analysis of cell regulatory substance indices of different chestnut varieties (clones) under drought stress

The accumulation of betaine was slower in plants under drought stress. There was no significant difference in betaine content among YZ, YQ and YSZF, and there was no significant difference among YSZF, QX42 and QL45. The difference in betaine content between YZ and YQ was significantly higher than that between QX42 and QL45, and the lowest content was observed for YL. The betaine content in YZ was 4.41 times that in YL, as shown in Table 4.

Table 4 Analysis of cell material index under drought stress.
Full size table

Soluble sugar is an important regulator of cell permeability. After drought stress treatment, the soluble sugar content of each chestnut variety (clone) was different. Table4 showed that the soluble sugar content of YQ and YSZF was significantly higher than that of QX42, QL45, YZ and YL. The difference in soluble sugar content of QX42 and QL45 was not significant, and the lowest soluble content was observed in YL.

Soluble proteins can dissolve in water, increase the water retention ability of fruit tree cells, and protect fruit trees from drought stress. Table 4 showed that YQ and YZ were significantly higher than QL45 and YL, and the latter two were significantly higher than YSZF and QX42. YQ, YZ, and QL45 were higher than YL, and YSZF and QX42 were lower than YL.

PRO is a common affinity osmoregulatory substance with no toxic effect. Under drought stress, intracellular regulatory substances actively accumulate, and the greater is the accumulation, the stronger is the drought tolerance. Table 4 showed that YL was significantly higher than the other varieties in PRO, and there was little difference between YSZF and YQ. Among them, the highest PRO content was observed in YL, and YZ was the cultivar with the lowest PRO content. It has a certain drought tolerance ability under drought stress, but the degree is different.

Chlorophyll is the main pigment for photosynthesis and can help plants store nutrients. There was no significant difference in chlorophyll content among YZ, YQ, QX42 and YSZF. The chlorophyll content of YL was significantly higher than that of YZ, YQ, QX42 and YSZF. The chlorophyll content of QL45 and YL was significantly higher than that of YZ, YQ and YSZF, as shown in Table 4. The highest chlorophyll content of YL was 2.40 mg/g, 1.17 times that of YSZF, as shown in Table 4.

Flavonoids have a strong oxidation effect and can remove too many oxygen free radicals. The flavonoid content of QL45 leaves was significantly higher than that of QX42, YSZF, YZ, YQ and YL, as shown in Table 4. The content of flavonoids in QL45 was 337.70 mg/g, which was 1.47 times that of YL. In Table 4, the flavonoid content of QL45 was significantly higher than the other five varieties (clones).

Flavonol was a kind of substance with large quantity and range in flavonoids, which can inhibit stomatal closure of plants. As shown in Table 4, the flavonoid content of YQ leaves was significantly higher than that of QX42, YSZF, YZ, YQ and YL. The content of flavonoids in YQ was 23,596.79ng/g, which was 1.47 times that of YL.

Analysis of harmful substance content in different chestnut varieties (clones) under drought stress

MDA is an important product of cell membrane lipid peroxidation, which can combine with proteins and enzymes on the cell membrane and destroy the cell membrane structure. The level of cell membrane lipid peroxidation can be reflected by the MDA content. Therefore, MDA content can be used as one of the detection indices of fruit trees subjected to stress damage. Table 4 showed that the MDA content of each chestnut variety was different under drought stress. The difference between YSZF, YZ and QL45 was not significant, and the difference between QL45, YQ, QX42 and YL was not significant. The MDA content of YSZF and YZ was significantly higher than that of YQ, QX42 and YL. The MDA contents of the six varieties (clones) of chestnut were YSZF, YZ, QL45, YQ, QX42 and YL from high to low, respectively. The MDA content of YSZF was 1.25 times that of YL (Table 5). Compared with the control group, these values were higher than that of YL, indicating that drought stress caused varying degrees of membrane lipid peroxidation, resulting in MDA accumulation.

Table 5 Analysis of harmful substance content under drought stress.
Full size table

Table 5 showed that compared with the OFR content, QX42 and YQ were significantly higher than YZ, YL and YSZF, and the latter three were significantly higher than QL45. Under normal conditions, the content was low and will not cause damage to cells. Under drought stress, it will accumulate too much and cause toxicity to cells. The OFR contents in leaves of six varieties (clones) from high to low were QX42, YQ, YZ, YL, YSZF and QL45, respectively.

The OFR production rate is too fast, which will cause an imbalance in reactive oxygen species production and removal in fruit tree cells. There was no significant difference in the OFR production rate between QX42 and YQ and no significant difference between YZ, YL, YSZF and QL45. QL45 was lower than other varieties (clones) (Table 5).

The H2O2 content is an important representative of the reactive oxygen species content, and the production and removal are balanced. Under drought stress, fruit tree cells produce a large amount of H2O2 to destroy the cell membrane, which reduces the yield of fruit trees. Table 5 showed that the difference in H2O2 content between QL45, QX42, QL and YZ was not significant, the difference in H2O2 content between YZ, YSZF and YQ was not significant, and the difference in QL45, QX42, QL and YL was significantly higher than that in YZ, YSZF and YQ.

Comparison of drought tolerance-related indices reflecting the drought tolerance of different chestnut varieties (clones)

The correlation coefficient matrix was obtained by analyzing the correlation of drought tolerance indices of six chestnut varieties (clones) (Table 6). The correlation coefficient matrix showed that there is a certain correlation between each drought tolerance index pair, and the information reflected by the data overlaps to a certain extent. The direct use of these indices cannot accurately evaluate the drought tolerance of chestnut. Through principal component analysis, the drought tolerance evaluation indices can be better selected by comparing the drought tolerance ability of each index. The principal component comprehensive evaluation can comprehensively and reasonably analyze the dominant factors in many factors, and the evaluation results are highly reliable. According to the principle that the cumulative contribution rate of principal component analysis was greater than 80%30, the first four principal components were extracted as comprehensive indices for drought tolerance evaluation of chestnut (Table 7). The first principal component contribution rate was 32.36%, the second principal component contribution rate was 24.88%, the third principal component contribution rate was 21.11%, and the fourth principal component contribution rate was 12.88%. The first four principal components reflected 91.24% of the information related to drought tolerance. SOD, betaine, soluble sugars, chlorophyll, and PRO had higher loads in the first principal component, indicating that the first principal component mainly reflected the relationship between these indicators and drought tolerance. CAT activity, APX activity, OFR content, OFR production rate, and flavone content had higher loads on the second principal component, indicating that the second principal component mainly reflected the relationship between these indices and drought tolerance. And MDA content, POD activity, soluble protein content, and H2O2 content had higher loads on the third principal component, indicating that the third principal component mainly reflected the relationship between these indicators and drought tolerance. The flavonol content in the fourth principal component has a high load, indicating that the fourth principal component mainly reflects the relationship between flavonol content and drought tolerance.

Table 6 Correlation coefficient matrix of 15 indicators.
Full size table
Table 7 Principal component coefficient, eigenvalue and contribution rate.
Full size table

The corresponding coefficients of each index in the four principal components are obtained by dividing the corresponding data by the square root of the corresponding eigenvalue of the principal component, and the four principal component expressions F1, F2, F3 and F4 are obtained. Then, the principal component comprehensive model was obtained by taking the proportion of the corresponding eigenvalue of each principal component in the total feature and the proportion of the extracted principal component as the weight31:

$${\text{F}} = \lambda_{1} /(\lambda_{1} + \lambda_{2} + \lambda_{3} + \lambda_{4} )F_{1} + \lambda_{2} /(\lambda_{1} + \lambda_{2} + \lambda_{3} + \lambda_{4} )F_{2} + \lambda_{3} /(\lambda_{1} + \lambda_{2} + \lambda_{3} + \lambda_{4} )F_{3} + \lambda_{4} /(\lambda_{1} + \lambda_{2} + \lambda_{3} + \lambda_{4} )F_{4} .$$

Specifically:

$$\begin{gathered} {\text{F}}_{{1}} = \, 0.{\text{352X}}_{{1}} - \, 0.0{\text{16X}}_{{2}} + \, 0.{\text{192X}}_{{3}} + \, 0.{\text{258X}}_{{4}} + \, 0.{\text{428X}}_{{5}} + \, 0.{\text{399X}}_{{6}} + \, 0.{18}0{\text{X}}_{{7}} \hfill \\ - 0.{\text{396X}}_{{8}} - \, 0.{32}0{\text{X}}_{{9}} + \, 0.0{\text{68X}}_{{{1}0}} + \, 0.0{\text{17X}}_{{{11}}} + \, 0.{\text{194X}}_{{{12}}} + \, 0.{\text{148X}}_{{{13}}} + \, 0.{\text{148X}}_{{{14}}} - \, 0.{23}0{\text{X}}_{{{15}}} \hfill \\ {\text{F}}_{{2}} = \, 0.0{\text{67X}}_{{1}} - \, 0.{\text{186X}}_{{2}} + \, 0.{\text{319X}}_{{3}} + \, 0.{4}00{\text{X}}_{{4}} - \, 0.00{\text{8X}}_{{5}} - \, 0.0{\text{99X}}_{{6}} - \, 0.0{\text{57X}}_{{7}} \hfill \\ - 0.{\text{142X}}_{{8}} + \, 0.{\text{236X}}_{{9}} + \, 0.{\text{425X}}_{{{1}0}} - \, 0.{\text{169X}}_{{{11}}} + \, 0.{\text{174X}}_{{{12}}} - \, 0.{\text{376X}}_{{{13}}} - \, 0.{\text{376X}}_{{{14}}} + \, 0.{\text{295X}}_{{{15}}} \hfill \\ {\text{F}}_{{3}} = \, 0.{3}0{\text{2X}}_{{1}} - \, 0.{43}0{\text{X}}_{{2}} + \, 0.{3}0{\text{6X}}_{{3}} - \, 0.0{6}0{\text{X}}_{{4}} - \, 0.{\text{178X}}_{{5}} + \, 0.0{\text{22X}}_{{6}} - \, 0.{\text{332X}}_{{7}} \hfill \\ - 0.0{\text{76X}}_{{8}} + \, 0.0{\text{13X}}_{{9}} + \, 0.0{5}0{\text{X}}_{{{1}0}} + \, 0.{26}0{\text{X}}_{{{11}}} - \, 0.{\text{372X}}_{{{12}}} + \, 0.{\text{285X}}_{{{13}}} + \, 0.{\text{285X}}_{{{14}}} + \, 0.{\text{326X}}_{{{15}}} \hfill \\ {\text{F}}_{{4}} = \, 0.{\text{185X}}_{{1}} + \, 0.{\text{262X}}_{{2}} - \, 0.{\text{159X}}_{{3}} + \, 0.{\text{159X}}_{{4}} - \, 0.0{\text{44X}}_{{5}} + \, 0.{\text{165X}}_{{6}} + \, 0.{4}00{\text{X}}_{{7}} \hfill \\ + \, 0.{\text{175X}}_{{8}} + \, 0.{\text{346X}}_{{9}} + \, 0.{\text{385X}}_{{{1}0}} + \, 0.{5}0{\text{8X}}_{{{11}}} - \, 0.{\text{311X}}_{{{12}}} - \, 0.00{\text{1X}}_{{{13}}} - \, 0.00{\text{1X}}_{{{14}}} - \, 0.0{\text{66X}}_{{{15}}} \hfill \\ \end{gathered}$$

The principal component synthesis model was as follows:

$$\begin{gathered} {\text{F }} = \, 0.{\text{239X}}_{{1}} - \, 0.{\text{119X}}_{{2}} + \, 0.{2}0{\text{4X}}_{{3}} + \, 0.{2}0{\text{9X}}_{{4}} + \, 0.{1}0{\text{2X}}_{{5}} + \, 0.{\text{143X}}_{{6}} + \, 0.0{\text{28X}}_{{7}} \hfill \\ - 0.{\text{172X}}_{{8}} + \, 0.00{\text{3X}}_{{9}} + \, 0.{2}0{\text{6X}}_{{{1}0}} + \, 0.0{\text{92X}}_{{{11}}} - \, 0.0{\text{13X}}_{{{12}}} + \, 0.0{\text{16X}}_{{{13}}} + \, 0.0{\text{16X}}_{{{14}}} + \, 0.0{\text{65X}}_{{{15}}} \hfill \\ \end{gathered}$$

The greater is the weight of each index, the greater is the importance, and vice versa. Therefore, among the drought tolerance indices of chestnut, the order of drought tolerance ability from high to low was as follows: SOD activity, APX activity, flavone content, CAT activity, PRO content, soluble sugar content, POD activity, betaine content, flavonol content, H2O2 content, protein content, soluble protein content, OFR content, OFR production rate, MDA content, chlorophyll content.

Comprehensive evaluation of drought tolerance (membership function method)

The membership function values of different chestnuts were obtained by the membership function method; then, the drought tolerance measurement values were obtained by combining the weight of each index (Table 8) to sort the drought tolerance. The greater is the drought tolerance measure, the stronger is the drought tolerance. Table 7 shows that the maximum drought tolerance was found in QL45 (1.333), and the minimum was found in YL (0.168). The average membership functions of 14 physiological indices in the leaves of 6 varieties (clones), i.e., QX42, QL45, YSZF, YZ, YQ and YL, were determined 1.090, 1.333, 0.773, 0.780, 0.910, 0.168 (Table 7). Based on this, the drought tolerance of six chestnut varieties (clones) was ranked, and the drought tolerance from high to low was as follows: QL45, QX42, YQ, YZ, YSZF, YL.

Table 8 6 subordinate function values, weights and orders of drought tolerance of chestnuts.
Full size table

Discussion

Drought stress creates reactive oxygen species which hampers the water retention of cells, causes cell swelling and destruction of the cell membrane19,32. Drought tolerance was evaluated by comprehensive comparison of physiological and biochemical indices among different chestnut varieties (clones). As a type of drought-resistant fruit tree, it is of great significance to identify the drought tolerance of chestnut for drought-resistant breeding and breeding varieties suitable for cultivation. There are many studies on chestnut enzyme activity, chlorophyll, membrane permeability, osmotic regulation and MDA under drought stress. Castanea mollissima was subjected to drought stress, and MDA and osmotic adjustment substances were increased, indicating that drought stress led to different varieties (clones) of Castanea mollissima to experience varying degrees of membrane damage as well as water loss stress. In addition, a positive response to drought stress was observed, this is consistent with the results of previous studies on chestnut31, lily33 and apricot34. The study of active oxygen protective enzymes in several chestnut varieties (clones) showed that SOD, POD, CAT and APX increased synchronously under drought stress, which indicated that they had synergistic effect, which was also consistent with the results of previous studies26,35. The soluble sugar content of different chestnut varieties (clones) under drought stress was different, and the content was higher than that of the YL of rootstock varieties.

There are various ways to evaluate the drought response of plants, and the genes that control the physiological response of plants to drought are also multiple genes. Different plants have different drought tolerance mechanisms, and the drought tolerance of the same plant is also different. Therefore, it is insufficient to evaluate the drought of plants by a single indicator, and there are certain limitations10,11,12,13,14,15; thus, multiple indicators should be integrated36. In this experiment, grafting trees of six chestnut varieties (clones) were used as experimental materials, and 15 drought tolerance indices were selected for correlation analysis. Principal component analysis was performed to determine the index weight. The component model shows that the drought tolerance indices of chestnut from high to low are as follows: SOD activity, APX activity, flavone content, CAT activity, PRO content, soluble sugar content, POD activity, betaine content, flavonol content, H2O2 content, protein content, soluble protein content, OFR content, OFR production rate, MDA content, chlorophyll content. Based on the membership function value, 15 physiological and biochemical indices were comprehensively analyzed by the principal component evaluation method. The obtained results showed that the drought tolerance of the 6 different varieties (clones) of chestnut was QL45 > QX42 > YQ > YZ > YSZF > YL. In the study of Guo Yan et al., YSZF and YQ drought tolerance belong to medium drought tolerance varieties. In this study, the drought tolerance of the same two varieties was not much different37. In addition, Wu and Guo et al. used potted seedlings and roots as the main indicators of comprehensive evaluation, and found that the drought tolerance of YL was greater than that of YSZF, which was contrary to the results of this study31. Growth states and index comprehensive evaluation results will be different, and twe will be a follow-up study on potted seedlings of chestnut with drought tolerance.

Multivariate statistics based on multiple indicators can accurately reflect the drought tolerance of different chestnut varieties (clones), providing a theoretical basis for the study of drought tolerance in chestnut varieties (clones).

Conclusion

In this study, principal component analysis and membership function method were used to study the drought tolerance of six chestnut varieties (Clones) based on different physiological indexes and two important secondary metabolites of flavonoids and flavonols as evaluation indexes. The results showed that flavonoids and flavonols ranked top among 15 tolerance -related indicators and could be used as tolerance indicators. The drought tolerance of 6 chestnut varieties (clones) was : QL45 > QX42 > YQ > YZ > YSZF > YL.

Materials and methods

Location

Hebei Normal University of Science and Technology Changli Shigezhuang Chestnut Germplasm Resources Nursery is located in the northeast of Hebei Province, southwest of Qinhuangdao. It is between 39° 22′ and 39° 48′ N and 118°45′-119°20′E and has the warm temperate semihumid continental monsoon climate in the eastern monsoon region of China. The annual average temperature is 11.8 °C, and the annual precipitation is 527.0 mm. The annual sunshine is 2719.5 h, annual sunshine percentage is 61%, frost-free period is 210 days, maximum frozen depth is 53 cm, and icing is 117 days.

Experimental design

The grafting trees of six chestnut varieties (clones) were Qianxi 42 (QX42), Qinglong 45 (QL45), Yanshan Zaofeng (YSZF), Yanzi (YZ), Yanqiu (YQ) and Yanlong (YL). In April 2019, 6 varieties (clones) were grafted onto rootstocks of the 5-year-old chestnut cultivar YL at the Changli Shigezhuang base. Each variety (clones) grafted 5 trees. In May 2021, 6 grafting trees of chestnut varieties (clones) in the resource nursery were selected as experimental materials. Unified conventional cultivation management was implemented, and a randomized block test was performed. Each variety (clones) was used to establish a single plant plot, which was repeated three times. In July 2021, the natural drought stress test was performed by watering one time. When the soil relative water content (RWC) decreased to approximately 30%, leaf samples were collected to determine the physiological and biochemical indices related to drought tolerance (the relative water content of agricultural drought grade standard soil was approximately 30%-40% for severe drought38.

Determination method

Three grafting trees of six chestnut varieties (clones) with basically the same growth management were selected to determine the soil relative water content. When the soil RWC was reduced to approximately 30%, mature leaves with similar sizes on the eastern main branch of chestnut varieties (clones) were collected on July 23, 2021, from 8:00 to 10:00 in the morning. Ten leaves were collected from each tree, and a total of 30 leaves were collected from one variety (clones). The labeled leaves were placed in an ice box and brought back to the laboratory. After liquid nitrogen treatment, they were stored in an ultra-low temperature refrigerator for testing.

Etermination of soil water content

The absolute soil water content was determined by drying and weighing method21, and the soil relative water content was measured by the JK-300F soil moisture meter of Xinghua Youke Instrument Co., Ltd.

Determination of physiological and biochemical indicators

The MDA content was determined with reference to thiobarbituric acid colorimetry; the activity of SOD was determined by reference to the nitroblue tetrazolium photoreduction method; POD activity was determined according to the guaiacol colorimetric method39. the chlorophyll content40, PRO content41, H2O2 content41, APX activity41 and CAT activity39 were determined by spectrophotometry; The soluble sugar content was determined by anthrone colorimetry. Soluble protein content was determined by coomassie brilliant blue method; the OFR content40 and OFR production rate41 was determined by the hydroxyl ammonia oxidation method; The betaine content was determined by chemical colorimetry41. The content of flavonoids and flavonol were determined by HPLC42; A single index cannot effectively evaluate the drought tolerance of chestnut. In this study, 15 physiological and biochemical indices related to drought tolerance were selected to comprehensively evaluate the drought tolerance of chestnut43. All of these indexes were measured by using the kit, which came from Suzhou Keming Biotechnology Co., LTD. in China.

Data analysis

Microsoft Excel 2010 was used for data integration. SPSS 26 and DPS 9.01 software were used for analysis of variance and principal component analysis. The membership function method was used to comprehensively evaluate the drought tolerance of chestnut. The related membership function values of each index and drought tolerance were calculated.

Membership function method

The calculation formula of the membership function value of each index and drought tolerance is as follows. When the index is positively correlated with drought tolerance, the calculation result is expressed as R (Xij) = (Xij-Xmin)/(Xmax-Xmin). When the index is negatively correlated with drought tolerance, the calculation results are expressed as R (Xij) = 1 − (Xij − Xmin)/(Xmax − Xmin), where R (Xij) is the drought tolerance membership function value of i variety j, Xij is the measured value of j index of i variety, and xmax and Xmin are the maximum and minimum values of j index. According to the membership function value and the weight of each index, the drought tolerance measurement value was calculated according to the membership function value and the weight of each index. The greater is the drought tolerance measure, the stronger is the drought tolerance ability. The calculation formula of drought tolerance (D) is as follows31:

$$\text{D}={\sum }_{j=1}^{m}(R({X}_{ij}){*W}_{j})$$

In the formula, Wj was the index weight.

Ethics statement

The experiment is in line with relevant institutions, national and international norms and legislation.