Individual Cryptomeria fortunei Hooibrenk Clones Show Varying Degrees of Chilling Stress Resistance

: Low temperature (LT) is an important abiotic factor a ﬀ ecting plant survival, growth and distribution. The response of Cryptomeria fortunei (Chinese cedar) to LT is not well known, limiting its application in production and ornamental value. In this study, we ﬁrst screened ten clones (#3, #25, #32, #42, #54, #57, #68, #66, #74, #X1), originating from ﬁve di ﬀ erent locations in China, for their degrees of cold resistance. We then selected the two showing the highest (#32) and lowest (#42) cold resistance to see the physiological and morphological response of di ﬀ erent cold-resistant C. fortunei clones to LT. We found that the electrolyte leakage of all ten clones increased strongly between 0 and − 8 ◦ C, while below − 8 or between 4 and 0 ◦ C did not yield additional increases. Under cold stress, clones #32 and #42 showed di ﬀ erent degrees of needle browning. From 25 to − 20 ◦ C, maximum and e ﬀ ective quantum yields of photosystem II ( F v / F m and YII) and photochemical and non-photochemical quenching (qP and NPQ) decreased continuously in two clones with decreasing temperature, where #42 was more strongly a ﬀ ected compared with #32. The chlorophyll content ﬁrst decreased signiﬁcantly to the lowest from 25 to − 12 ◦ C, then increased signiﬁcantly at − 16 / − 20 ◦ C compared with − 12 ◦ C. We observed changes in needle cellular ultrastructure at − 8 ◦ C, with chloroplasts of #32 swelling, while those of #42 were destroyed. Correlation analysis indicated that needle browning and chlorophyll ﬂuorescence were closely related to temperature, and cellular ultrastructure changed notably around semi-lethal temperature (LT 50 ), which can be used as physiological indicators for the identiﬁcation of cold resistance. We found a clear di ﬀ erence in cold tolerance between clones of #32 and #42, with #32 being more tolerant, which can be exploited in breeding programs. We conclude that strongly cold-resistant clones have more stable physiological states and a wider adaptability to LT compared with weak ones.


Introduction
Low temperature (LT) is one of the most important abiotic stress factors affecting plant survival and growth, and it is also the main factor that determines the geographical distribution and growth season of plants [1]. Cold resistance is a complex feature and can be a major influence on plant survival [2]. During their life cycles, most terrestrial plants that are grown in temperate or cold regions are affected

LT Treatment
In June 2019, one-year semi-lignified branches with no visible damage from 10 clones of C. fortunei (5-year-old cuttage seedlings, at heights of 100-120 cm) were collected and transferred to our laboratory in Nanjing (118 • 50' E, 32 • 05' N) at ambient temperature in the dark, a process that took around 2-3 h. Branches were rinsed with deionized water and surface-dried with clean filter papers. Material from each clone was divided into 8 groups for different treatments. Each treatment was performed in triplicate (10 × 3 branches) and packed in a 0.03 mm thick polyvinyl chloride (PVC) storage bag. One group, placed in the dark at room temperature (25 • C), served as control (CK), the other groups were placed in a program-controlled refrigerator (BCD-405WDGQU1, Haier, Qingdao, Shandong Province, China) for LT treatment (4, 0, −4, −8, −12, −16 and −20 • C) for 12 h. The temperature was decreased at a rate of 4 • C/h, starting from 4 • C down to the aimed temperature [36]. After 12 h treatment, samples above 0 • C (25 and 4 • C) could be directly used for the measurement of physiological indexes, whereas samples treated at and below 0 • C were raised to 4 • C at a rate of 4 • C/h, then the physiological indexes were analyzed. The same treatments were used for #32 and #42 clones for different time points (2 and 24 h).

Electrolyte Leakage (EL) and Semi-Lethal Temperature (LT 50 )
EL was measured according to the method of Tsarouhas et al. [37] and Ma et al. [38] with some modifications. Briefly, the needles of LT-treated samples were collected and incubated in 20 mL of distilled water at room temperature. After 12 h incubation, 0.2 g of needles was used, and its conductivity was measured as E 1 with a conductivity meter (DDS-307, Leici Instruments Co., Shanghai, China), setting the conductivity of distilled water as E 0 . The samples were boiled for 20 min. After cooling to room temperature, the ultimate conductivity was measured as E 2 . EL was calculated using the following formula: A series of EL was combined with the exponential equation y = k / (1 + ae −bt ) to obtain intermediate parameters for the plant LT 50 . Within this formula, y is EL, k is the maximum EL, t is temperature, e (about 2.718) is constant, and a and b are values calculated by the equation. Let Y = ln (k − y) / y, the equation is transformed to Y = ln a − bx, then x = (ln a) / b and x is the LT 50 [39,40].

Needle Phenotype Observations
Fifteen branches of #32 and #42 clones were treated in the same way as above for 2, 12 and 24 h, respectively, and then placed in a 4 • C refrigerator (BCD-405WDGQU1, Haier, Qingdao, Shandong Province, China) for 12 h. Branches placed at 25 • C were used as controls (CK). The needles were phenotyped. The browning degrees were classified according to the depth and areas of browning, then the browning index was calculated [39]. The browning degrees were divided into 5 grades: grade 1, nonbrowning, green; grade 2, slight browning, 0-1/4 area of browning on the surface; grade 3, 1/4-1/2 area of browning on the surface; grade 4, 1/2-3/4 area of browning on the surface; grade 5, greater than 3/4 area of browning on the surface, dark brown. The browning index was calculated by browning grade and browning area, and the equation was Browning index (%) = (each browning grade × percentage of browning area of each grade) (the highest grade of browning × total number of branches)

Chlorophyll Content
The needles of #32 and #42 clones were collected after LT treatment as described above for 2, 12 and 24 h. Chlorophyll content was measured using the alcohol extraction method described by Lichtenthaler and Wellburn [41]. The chlorophyll was extracted by grinding each sample (about 0.1 g) with calcium carbonate, quartz sand and 96% alcohol and then incubating for several minutes until tissues were blanched. After centrifuging (Allegra X-22R, Beckman Coulter, Inc, Breabbe, California, USA) 3 times at 4000 rpm for 20 min each time at room temperature, the OD 665 (optical density at 665 nm) and OD 649 (optical density at 649 nm) of the supernatant were determined using a spectrophotometer (V-1200, Shanghai Mapada Instruments Co., Ltd., Shanghai, China). The concentration of chlorophyll (mg g −1 fresh weight) was calculated using Lichtenthaler and Wellburn equations [41].

Chlorophyll Fluorescence Characteristics
The needles of #32 and #42 clones were collected after LT treatment as described above for 2, 12 and 24 h, then the chlorophyll fluorescence transient was determined using a fluorimeter (Dual-PAM/F-100, Walz, Effeltrich, Baden-Wuertenberg, Germany) [42]. The third to the fifth mature needles from top to bottom were placed in the dark for 30 min before measurements. A saturating light (3500 µmol m −2 s −1 ) along with actinic light (300 µmol m −2 s −1 ) was applied, and the dual-channel mode (fluo + P700) was adopted to measure the induction curve. The fluorescence signal was stable for 5 min under photochemical conditions, and the following fluorescence parameters were measured as indicators of the functionality of the photosynthetic apparatus: initial fluorescence (F o ), maximum fluorescence (F m ), , the effective quantum yield of PSII (YII), photochemical quenching (qP) and non-photochemical quenching (NPQ).

Ultrastructural Observations
Branches of #32 and #42 treated at −4, −8 and −20 • C for 24 h were collected, and samples at 25 • C were used as a control (CK). For each treatment, three needles (the third to the fifth mature needles from top to bottom) were sampled and fixed in 3% glutaraldehyde solution in phosphate buffer with a pH of 7.0 for 24 h at room temperature, followed by secondary fixation in 2.5% osmium tetroxide solution for 2 h. The specimens were rinsed and dehydrated in a graded series of alcohols and acetone of increasing strength until 100% and embedded in turn in a gradient of absolute ethanol and Poly Bed 812 epoxy resin (3:1, 1:1, 1:3, 0:3). The samples were stored dry for 8 h. Ultrathin sections (80−100 nm) were obtained using an ultratome III ultramicrotome (LKB, Bromma, Stockholm Province, Sweden) and observed using transmission electron microscopy (TEM) (JEM-2100, JEOL Ltd., Tokyo, Japan) after being stained in uranyl acetate and lead citrate [43]. From the electronic microscope images, 10 chloroplasts were selected for observation and measurement. The chloroplasts were approximated as ellipses, then the lengths of the semi major and semi minor axes of the ellipse were measured, respectively, and the approximate areas of the transverse sections of the chloroplasts in mesophyll cells were calculated. The elliptical area was calculated by the following equation: where S is the elliptical area (µm 2 ); π is pi; and a and b are the length of semi major axis (µm) and semi minor axis (µm), respectively.

Statistical Analysis
The experimental data were first collected, then statistically analyzed and plotted using Excel 2003 (Microsoft Co., Redmond, Washington, USA). All data shown in this study are mean values and standard errors (SEs). The data checked for the normal distribution (p-p chart) were statistically analyzed using one-way analysis (temperature as a factor) of variance in SPSS v25.0 (SPSS, Inc, an IBM Co., Chicago, Illinois, USA), and Duncan's multiple range test at a confidence level of 95% (p < 0.05) was adopted to define the statistical significance of the difference [44].

Individual C. fortunei Clones Show Varying Degrees of Cold Resistance
We aimed to preselect the most and least cold resistant out of 10 C. fortunei clones (isolated from 5 places in China) by analyzing their electrolyte leakage (EL) and semilethal temperature (LT 50 ) values in response to cold stress. As temperature decreased, all clones showed a similar S-curve response of increasing EL values ( Figure 1). The EL of needles of most C. fortunei clones increased slowly from 4 to 0 • C, and only #42 significantly increased by 121.05%. Whereas, the EL of each clone rapidly increased between 0 and −8 • C to 61%-73%. Any further reduction in temperature below −8 • C did not yield additional increase in EL values ( Figure 1). We used an exponential equation to derive a response curve function for each of the clones, which gave a significant fitting degree (R 2 ) in all cases ( Table 2). From these curve functions we could derive that clone #32 showed the lowest LT 50 (−7.8 • C), while #42 showed the highest (−4.7 • C) ( Table 2). We then selected these two clones for further experiments in order to study the widest possible range of cold-resistant responses within C. fortunei [45][46][47].

Low-Temperature Stress (LTS) Causes Needle Damage in C. fortunei Clones
Cold acclimation often leads to changes in leaf color and shape [48,49]. To test how and whether needle phenotypes related to cold resistance, we calculated the browning indexes of needles of clones #32 and #42 under LTS at a series of decreasing temperatures. After 2h LT treatment, needles of both #32 and #42 did not show significant differences in both needle shape and color at temperatures ≥−4 °C, compared to the control (25 °C), yet at ≤−8 °C they showed varying degrees of chlorosis (Figure 2a,b, Table 3). After treating for 12 or 24 h, only needles kept at ≥0 °C showed no obvious changes in needle color, while chlorosis started to appear at ≤−4 °C (Table  3). Furthermore, this prolonged treatment below −12 °C also caused top needles to wilt and droop, making them thicker and easier to detach (Figure 2c-

Low-Temperature Stress (LTS) Causes Needle Damage in C. fortunei Clones
Cold acclimation often leads to changes in leaf color and shape [48,49]. To test how and whether needle phenotypes related to cold resistance, we calculated the browning indexes of needles of clones #32 and #42 under LTS at a series of decreasing temperatures. After 2h LT treatment, needles of both #32 and #42 did not show significant differences in both needle shape and color at temperatures ≥−4 • C, compared to the control (25 • C), yet at ≤−8 • C they showed varying degrees of chlorosis ( Figure 2a,b, Table 3). After treating for 12 or 24 h, only needles kept at ≥0 • C showed no obvious changes in needle color, while chlorosis started to appear at ≤−4 • C (Table 3). Furthermore, this prolonged treatment below −12 • C also caused top needles to wilt and droop, making them thicker and easier to detach (Figure 2c-f). With a decrease in temperature and prolongation of LT treatment, the degree of browning damage increased significantly compared with the control (Figure 2, Table 3). In addition, any phenotypic changes in response to cold stress were less obvious in clones #32 than in #42, consistent with the fact that #32 showed a lower LT 50 than that of #42 (Table 2).

Chlorophyll Contents of the C. fortunei Clones Change under Different LTS
Since we found that increasing LTS caused a greater degree of needle browning in C. fortunei, we sought to determine whether the chlorosis we observed was related to a decrease in chlorophyll content and low cold resistance. We found that with decreasing temperature, chlorophyll a, b and a + b content of #32 and #42 needles would initially decline steadily from 25 to −12 • C (p < 0.05), only to rise again at −16 and/or −20 • C (Figure 3). With the prolongation of LT treatment time (2, 12, 24 h), the chlorophyll content of clone #32 was rather stable, while that of #42 needles decreased further (Figure 3). Chlorophyll a, b and a + b content of clone #42 needles decreased significantly faster than those of clone #32. These results suggest that the variations of clone #32 were smaller and the change processes were slower at LT for no more than 24 h, indicating that clone #32 was more stable to cold stress than clone #42.
we sought to determine whether the chlorosis we observed was related to a decrease in chlorophyll content and low cold resistance. We found that with decreasing temperature, chlorophyll a, b and a + b content of #32 and #42 needles would initially decline steadily from 25 to −12 °C (p < 0.05), only to rise again at −16 and/or −20 °C (Figure 3). With the prolongation of LT treatment time (2, 12, 24 h), the chlorophyll content of clone #32 was rather stable, while that of #42 needles decreased further (Figure 3). Chlorophyll a, b and a + b content of clone #42 needles decreased significantly faster than those of clone #32. These results suggest that the variations of clone #32 were smaller and the change processes were slower at LT for no more than 24 h, indicating that clone #32 was more stable to cold stress than clone #42. (c) chlorophyll a + b content. Each value in the graph is the mean ± standard error (n = 6). Different patterns indicate different treatment temperatures, and lowercase letters above each bar indicate significant differences (p < 0.05).

Chlorophyll Fluorescence of C. fortunei is Negatively Affected by LT
We next tested whether the needle browning we observed in C. fortunei clones under cold stress indeed resulted in lowered chlorophyll fluorescence by analyzing their Fv/Fm, YII, qP and NPQ in detail. We found that Fv/Fm ratios of #32 and #42 needles were reduced with a decrease in temperature and prolongation of LTS time (Figure 4a

Chlorophyll Fluorescence of C. fortunei is Negatively Affected by LT
We next tested whether the needle browning we observed in C. fortunei clones under cold stress indeed resulted in lowered chlorophyll fluorescence by analyzing their F v /F m , YII, qP and NPQ in detail. We found that F v /F m ratios of #32 and #42 needles were reduced with a decrease in temperature and prolongation of LTS time (Figure 4a). For clone #32, the F v /F m of the needles decreased rapidly below −12 • C when treating for 2 h or between −8 and -20 • C when treating for 12 or 24 h. For clone #42, under cold treatment for 2 h, F v /F m of the needles decreased sharply between −8 and −12 • C, while at 12 and 24 h treatment times, this effect was seen between −4 and − 8 • C (Figure 4a). The F v /F m value of clone #32 at −20 • C was similar to or higher than that of clone #42 at −12 • C. The trends of YII, qP and NPQ were almost the same as that of F v /F m , showing a stable trend (from 25 to −4 • C), then rapid decline (from −4 to −8/−12/−16 • C) and finally a stable trend. For temperatures between −4 and −8/−12 • C, the decreasing range of #32 was less than that of #42 and the value of #32 was higher than that of #42 (Figure 4b All these results suggest that clone #32 may cope with cold stress better than clone #42. 12 or 24 h. For clone #42, under cold treatment for 2 h, Fv/Fm of the needles decreased sharply between −8 and −12 °C, while at 12 and 24 h treatment times, this effect was seen between −4 and − 8°C (Figure 4a). The Fv/Fm value of clone #32 at −20 °C was similar to or higher than that of clone #42 at −12 °C. The trends of YII, qP and NPQ were almost the same as that of Fv/Fm, showing a stable trend (from 25 to −4 °C), then rapid decline (from −4 to −8/−12/−16 °C) and finally a stable trend. For temperatures between −4 and −8/−12 °C, the decreasing range of #32 was less than that of #42 and the value of #32 was higher than that of #42 (Figure 4b-d

Mesophyll Cells of C. fortunei Needles Undergo Ultrastructural Changes under LTS
Since we observed eventual structural collapse of needles at very low temperature, we sought to test whether the cellular structure of needles was affected by LTS using electron microscopy. Indeed, we found that the ultrastructure of C. fortunei mesophyll cell showed significant changes under decreasing LTS ( Figure 5). At 25 • C, the cell nucleus, cell walls, chloroplasts and mitochondria had a normal structure with clearly visible membranes and homogeneous karyoplasm, and distinct nucleoli could also be observed (Figure 5a-f). Large central vacuoles occupied most of the cellular space, squeezing other organelles to the edge (Figure 5a,d,e). A large number of well-developed chloroplasts were distributed along the long axis on the inner side of the plasma membrane, showing dense stroma and regularly stacked grana thylakoids. A few osmiophilic plastoglobules and starch granules appeared in the chloroplasts, which were themselves densely surrounded by mitochondria (Figure 5c,f).  When the temperature decreased to −4 • C, the majority of cell organellar structures remained clear and complete (Figure 5g-l). No obvious difference in ultrastructure was observed in needles derived from clone #32 compared to the control condition (25 • C). However, we did observe that most chloroplasts became swollen. In the chloroplasts, the number of plastoglobules markedly increased, whereas starch granules reduced significantly (Figure 5g-i). Additionally, we observed dilated thylakoid lumina. Similar changes, yet to a greater extent, were observed in clone #42 than in #32: the cytoplasmic membrane was partly invaginated (Figure 5j), multiple vacuoles were found ( Figure 5k) and their inclusion increased sharply, most chloroplasts were swollen into the hemisphere and some moved towards the center of cells (Figure 5j,l) and vacuolizing mitochondria were also visible ( Figure 5l).
Obvious changes in ultrastructures were observed as the temperature was further reduced to −8 • C, such as unclear thylakoids, the separation of chloroplasts from the cell wall and serious plasmolysis (Figure 5m-r). In the cells of #32, the nucleus was slightly blurred, and vacuole inclusions increased significantly (Figure 5n). The chloroplasts were further swollen and deformed (Figure 5o). Intervals between stroma lamellae, a small amount of plastoglobules as well as vacuolizing mitochondria were visible (Figure 5n,o). The cells of #42 displayed more serious damage. Most of the mesophyll cells showed a greater degree of plasmolysis (Figure 5p). The cellular structure was scattered, the nucleus blurred, the cell wall folded (Figure 5p) and vacuole membrane ruptured, leading to contents flowing out. Some chloroplasts were swollen, distorted, and even their bilayer membranes were damaged (Figure 5r). Starch granules were not observed in cells, and the number of mitochondria decreased dramatically (Figure 5q,r).
At −20 • C, the ultrastructures of the two clones had the most severe defects (Figure 5s-x). Most organelles were collapsed and degraded. The intact structure of chloroplasts was lost. The stroma lamellae were broken, discrete, elongated and thinned, showing a severe filamentous layer structure. Mitochondria were also decomposed and disappeared. Based on these results, we concluded that the ultrastructure of cells indeed responded to cold stress, and the clone (#32) more resistant to cold stress appeared less affected than the sensitive one (#42).

Chloroplast Size of C. fortunei is Affected by LTS
Plant photosynthesis depends not only on the increase in the number of chloroplasts in cells, but also on the expansion of their cell surface area [50]. To determine whether the sizes of chloroplasts we observed in mesophyll cells was changed by an increasing degree of LTS, we turned to directly measuring cross-sections of chloroplasts. Under LTS for 24 h, with decreasing temperature, no significant differences in the lengths of the long or short axes of the chloroplasts within each clone were found ( Table 4). The cross-section areas of chloroplasts of #32 and #42 first ascended and then descended, and the minimum values were 14.55 and 13.25 µm 2 at −8 and −4 • C (the temperature with the minimum cross-section area was consistent with the LT 50 ), respectively (Table 4). Thus, we speculated that as the temperature (>LT 50 ) decreased, the cross-section areas of chloroplasts reduced; at LT 50 or lower, the integrity of the chloroplast will be destroyed, and the cross-section area of the chloroplast will be further expanded. The cross-section areas of chloroplasts from clone #32 were not significantly affected by increasing intensity of LTS, while chloroplasts from #42 did show a significant response to LT (p < 0.05) ( Table 4). We therefore conclude that, in general, clone #32, which had stronger cold resistance, had smaller changes to chloroplast size. Values are mean ± standard error, n = 10, lowercase letters in the same column indicate significant differences (p < 0.05).

Discussion
Increased tolerance to biotic and abiotic stresses in plants is necessary to increase productivity under conditions with LT, high salinity, limited water supply and so forth [51]. In this study, the physiological responses to cold stress were studied to identify key elements that may be responsible for LT tolerance in C. fortunei.
Plant cell membranes play a crucial role in maintaining a stable intracellular environment [52]. Under LTS, the integrity of the cell membrane is destroyed, causing leakage of intracellular soluble substances and an increase in EL. There is a positive correlation between the rate of lipid peroxidation and EL [53]: the greater the LTS, the greater the relative permeability of plasma membrane and EL. In this study, the EL of C. fortunei clones showed a typical rising S-curve in response to decreasing temperatures, with a slow increase in EL between 4 and 0 • C, similar to the results obtained using Salix eriocephala [37]. During initial or mild injury (sublethal dose) and adverse stress of C. fortunei, there was a stage in which EL occurred with a slight decrease or no significant change at 0 • C, compared to 4 • C, because of the self-protection mechanism of plants. During the frosting process (from 0 to −8 • C) EL was found to increase sharply, which is consistent with previous reports on Forsythia × intermedia ('Lynwood Gold' forsythia) cultivars [54]. It might be because at this temperature, the biofilm system is seriously damaged, and its permeability is sharply increased, resulting in a large amount of extravasation of electrolytes [55]. As the temperature was further reduced, we found that EL did not change significantly, which was consistent with research on cherry rootstocks [56]. It is possible that most cell damage had reached the saturation state and caused irreversible damage. Reports in tea leaves [57] and evergreen alpine shrubs [58] showed that the LT 50 (at LT) can reflect cold resistance accurately: the lower the value, the stronger the cold resistance of plants. Therefore, according to the LT 50 value, we designated #32 as having the highest cold resistance and #42 as having the lowest cold resistance.
Our results show that needle shapes and colors changed significantly in response to cold stress, and they were extremely significantly negative (p < 0.01) related to cold resistance (temperature) of C. fortunei (Appendix A, Table A1). This was similar to results obtained from previous studies on Fagus sylvatica (European beech) [59], in which leaves did not remain freshly green as a result of freezing stress. We found that, with a further decrease of LT (≥−12 • C) and the prolongation of LT treatment, the chlorophyll contents of our two C. fortunei clones reduced by varying degrees. This was similar to studies performed on Eupatorium adenophorum (crofton weed) [60], rice [61] and Dactylis glomerata (orchardgrass) [62], supporting that, under freezing stress, chlorophyll content decreases as temperature decreases. We also found that within a certain temperature range (from 25 to −12 • C), the changes of needle color were related to chlorophyll content (Appendix A, Table A2). It may be that free radicals cause damage to the plant chloroplasts' membrane [63] and inhibit the synthesis of light-harvesting complexes in chloroplast lamellae as well as protochlorophyllide oxidoreductase [64], leading to a blockade in chlorophyll synthesis. On the other hand, LT can enhance chlorophyll degradation and further reduce the chlorophyll content [65,66], resulting in a high degree of plant needle tissue necrosis. Interestingly, we found that chlorophyll content showed a resurgence at −16 and/or −20 • C. We think this may be caused by needle wilting and excessive water loss, resulting in a high content of chlorophyll per gram of fresh weight, and might be due to a too low temperature, causing the needle to freeze quickly and leaving no time for the chemical reaction. In general, the LT change mechanism of the chlorophyll content is more complex [67].
We found that F v /F m continuously decreased as stress increased. This agrees with much of the research reported in Eupatorium adenophorum, sunflower and apricot fruit, which states that LTS can cause damage to the PSII reaction center in plants, resulting in a continuous decline in F v /F m [60,68,69]. A series of fluorescence parameters, such as yield, qP and NPQ, related to PSII electron transfer ability also decreased in the leaves of C. fortunei under LTS. This shows that LT not only causes changes to the structure of the light harvesting antenna and reaction center, but it also reduces the chemical conversion efficiency of PSII [70] and inhibits the electron transfer of PSII. We speculated that the basic reason may be that LT causes the changes to the structure or conformation of the reaction center, light harvesting antenna, electron transfer body and other proteins on the thylakoid membrane, and finally causes the omnidirectional damage and suppression of the PSII reaction center [20]. Fluorescence parameters were closely correlated to the cold resistance of C. fortunei (Appendix A, Table A1). The change of F v /F m is widely used and relatively clear, making it a suitable indicator of cold tolerance. At the same time, combined with the temperature when F v /F m > 0.5, we speculated that #32 can tolerate −12 • C for 2 h or −8 • C for a long time (12, 24 h). In contrast, 42# can tolerate −8 • C for 2 h or −4 • C for a long time. We also think the tolerable temperature may be close to LT 50 . Future experiments could include different seasons for a more accurate and detailed record of cold resistance. At the same time, different tissues of C. fortunei clones can be sampled to study whether physiological changes might differ across organs.
We found that with increasing cold stress, the cross-section area of chloroplasts decreased first, then increased. Chloroplast enlargement appears to occur in both cold-resistant plants such as Arabidopsis thaliana (thale cress) and cold-sensitive plants such as mung bean [71]. Chloroplasts are believed to be the most sensitive organelles [72]. An osmotic imbalance between stroma and cytoplasm may be the cause of chloroplast swelling [28,73]. Slight swelling of the chloroplasts may reduce the activity of starch degradation enzymes and soluble sugar content, thereby reducing the amount of photosynthetic product in the chloroplast as well as the matrix water potential [28]. At LT 50 or lower, the chloroplast further expands, resulting in the loss of its intact structure and disorganization of its outer membrane, resulting in an increase in its area. Our study showed that LTS caused plasmolysis as well as numerous invaginations of the plasma membrane, vacuolation of plastids and mitochondria, condensation of chromatin in the nucleus, multiple vacuoles and an increase in their inclusion in C. fortunei needles. Similar alterations in cellular ultrastructure have been observed by cold stress in leaves of Camellia sinensis (tea) [74], Arabidopsis thaliana [75], winter wheat cv. Volodarka [76], as well as one-year-old young grape plants (Vitis vinifera) [77]. Ultrastructural cold damage to plant cells is significantly similar to the ultrastructural changes in the programmed death process, and the physiological response of plant cells exposed to LT (above-freezing temperatures) also resembles that seen in putative programmed cell death [28]. The reasons for mitochondrial damage might be that mitochondria accumulate a large amount of reactive oxygen species, which could damage mitochondrial cristae, decreasing the mitochondrial energy production rate and the metabolism of mesophyll cells, leaf tissue or the whole plant [78]. The phenomena of multiple vacuoles and their increased inclusion may be related to vacuole-clearing, senescent organelles and cells under LT conditions. Alternatively, it may be a general plant stress response. The nucleus is the regulatory center of physiological and biochemical processes for cells, and damage to the nucleus results in cells being unable to maintain their normal metabolic processes, ultimately leading to cell death [79]. We did find that the majority of cell organellar structures of #32 and #42 remained clear and complete at −4 • C, so we think the needles of both clones were not damaged irreparably and can tolerate −4 • C for 24 h. Whereas, the cellular ultrastructure changed obviously at −8 • C, and for #32, chloroplasts were swollen, the nucleus was slightly blurred and vacuole inclusions increased (the structures were complete), while those structures (most mesophyll cells, cellular structure, nucleus, cell wall, vacuole membrane, chloroplast) of #42 were damaged seriously or even destroyed. Therefore, we speculated that #32 was able to endure at −8 • C for 24 h, while #42 can tolerate −4 • C for 24 h, and we think the tolerable temperature may be also close to LT 50 . In addition, metabolic disruptions induced by environmental factors are manifested in the ultrastructure of cell organelles, and cellular components have varied tolerance to LT [80]. Similarly, under LTS, damage inflicted on the same types of organelles was different between cells, which ensured that cells can maintain normal physiological function for a period of time under external stress, and it may be beneficial to the survival of C. fortunei at moderate levels of LTS. At temperatures lower than the LT 50 , the specific effects of different LT with small gradients (such as 1 • C, 2 • C) on the ultrastructural changes of plants need further study.

Conclusions
We found the LT 50 of 10 tested clones of C. fortunei to be between −4.7 and −7.8 • C, with #32 having the lowest value at −7.8 • C and #42 the highest at −4.7 • C, suggesting that #32 and #42 have the strongest and weakest cold resistances, respectively. Furthermore, under LTS, the degree of detrimental effects on needle browning and chlorophyll fluorescence parameters are closely correlated to the amount of cold resistance that a C. fortunei clone has. In a certain temperature range (from 25 to −12 • C), chlorophyll content is closely related to cold resistance. Needle browning is the most sensitive indicator of LTS or of the difference between #32 and #42, as needles begin to turn reddish brown above LT 50 . The changes in chlorophyll are more complex compared to needle browning or other parameters, and chlorophyll contents reduced to varying degrees as the temperature decreased (≥−12 • C) and then recovered at −16 and/or −20 • C. Meanwhile, we found that the parameters most closely correlated to the LT 50 are F v /F m and ultrastructure changes, suggesting that F v /F m and ultrastructure changes are more suitable as physiological indicators for cold resistance. Clone #32 can tolerate −8 • C for a long time (12 h, 24 h) or −12 • C for a short time (2 h). Similarly, #42 can tolerate −4 • C for a long time or −8 • C for a short time. The physiological states of #32 were more stable and have wider adaptability to LTS. It may be useful for breeders to breed a cultivar that is more resistant to chilling conditions by inbreeding clone #32 to pursue and isolate LT-dominant plants or doing more molecular studies between #42 and #32.