Drought Hardening Contributes to the Maintenance of Proportions of Non-Embolized Xylem and Cambium Status during Consecutive Dry Treatment in Container-Grown Seedling of Japanese Cedar ( Cryptomeria japonica

: Climate models in Japan predict that the annual mean air temperature and number of consecutive dry days will increase in the future, leading to high seedling mortality rates. Maintaining high survival rates of Cryptomeria japonica seedlings, a commercially important tree species, is therefore, important in terms of appropriate forest management under climate change. Although drought hardening, in which seedlings are acclimated to dry conditions in the nursery prior to planting, contributes to increased survival under drought conditions, little is known about the e ﬀ ective irrigation frequency of drought hardening in C. japonica seedlings. In this study, we therefore, examine the e ﬀ ectiveness of di ﬀ erent drought-hardening treatments in C. japonica . We ﬁrst clarify the e ﬀ ects on physiological and morphological traits by comparing three drought-hardening treatments [control ( C ): Irrigation once daily; mild ( M ): irrigation once every three days; and severe ( S ): irrigation once every ﬁve days] for one month. Next, to conﬁrm the e ﬀ ects during consecutive dry treatment, we stopped irrigation for 13 days and once again compared the physiological traits between the three drought-hardening treatments. Drought hardening reduced whole-plant transpiration ( E whole ), resulting in conserved water use, and this tendency was particularly evident under the S treatment. Moreover, during consecutive dry treatment, the E whole , proportions of non-embolized xylem, and cambium status of basal stem regions were maintained for the longest duration under the S treatment, followed by the M treatment. Our ﬁndings suggest that the e ﬃ ciency of drought hardening increased with drought severity. Furthermore, one month of drought hardening contributed to both water conservation and the maintenance of cell di ﬀ erentiation under consecutive dry treatment, likely increasing the tolerance and survival of C. japonica seedlings under prolonged drought.


Introduction
Cryptomeria japonica D. Don (Japanese cedar) is a commercially important tree species in Japan, covering 36.4% of the total forested area [1]. According to afforestation data from 2006 to 2014 for compensation of tree die-off due to natural disasters, the area of afforestation resulting from drought is estimated at 1919 ha (14.9% of the total area of afforestation) [2]. In Japan, drought-induced annual mortality rates are particularly high in 1-to 5-year-old seedlings, at approximately 27% [3]. As approximately 85% of the afforestation costs for C. japonica are spent on seedlings [4], the high rate

Plant Materials
Three-year-old container-grown C. japonica seedlings (Ringyo-Shubyo Association, Ibaraki, Japan) were used. The seedlings were grown in the field for 1 year and transplanted to 150-cc multi-cavity containers (OY-150; Zenbyoren, Tokyo, Japan) filled with coir peat. They were then periodically fertilized using solid and liquid fertilizer and were irrigated daily in the field for 1 year. After the development of leaves and cessation of xylem development (start of September 2018), the seedlings were transplanted into 200 cc single-cavity pots filled with coir peat and placed in a phytotron chamber under an artificial illumination system (Koitotron KG; Koito Industries, Yokohama, Japan) at the Forestry and Forest Products Research Institute (Tsukuba, Japan; 36 • 01 N, 140 • 08 E) for 1 month. The seedlings were maintained at day and night temperatures of 28 • C and 25 • C, respectively, with a linear increase from 25 • C to 28 • C for 1 h followed by 28 • C for 10 h, a linear decrease to 25 • C for 1 h, and then 25 • C for 12 h. Humidity was maintained at 60% in the day and 70% at night, respectively, throughout the experiment, with a linear decrease from 70% to 60% for 1 h followed by 60% for 10 h, a linear increase to 70% for 1 h, and 70% for 12 h. Day and night PPFD was maintained at 1200 and 0 µmo m −2 s −1 , with a linear increase from 0 to 1200 µmo m −2 s −1 for 1 h followed by 1200 µmo m −2 s −1 for 10 h, a linear decrease to 0 µmo m −2 s −1 for 1 h, and 0 µmo m −2 s −1 for 12 h. Before drought hardening, the seedlings were fully irrigated and fertilized with liquid fertilizer in the form of 0.5 L 1000× Hyponex solution (22 mg N L −1 , 36 mg PO 4 3− L −1 , and 18 mg K L −1 ; 6:10:5, N:P:K; Hyponex Japan Corp., Ltd., Osaka, Japan) for 1 month to minimize the effects of transplanting and transfer on seedling growth. The position of the seedlings was randomly changed approximately once a week to minimize the effect of micro-environmental factors within the phytotron chamber.

Drought Hardening and Consecutive Dry Treatment
A total of 36 seedlings [65.6 ± 6.6 (mean ± 1SD) cm in height with a basal stem diameter of 7.6 ± 1.0 mm (mean ± 1SD)] were transferred to the phytotron chamber. All seedlings satisfied the standard height (30 cm height) and basal stem diameter (3.5 mm diameter) for afforestation in Japan. Twelve seedlings were subjected to each of three drought-hardening treatments for 1 month: irrigation once daily [control (C) treatment], irrigation once every 3 days [mild drought (M) treatment], and irrigation once every 5 days [severe drought (S) treatment] ("drought hardening"). To minimize the effects of water distribution heterogeneity in the soil, irrigation was carried out by submerging the pots in water for 45 min. The volumetric soil water content (SWC) was then measured in four pots per treatment using a soil moisture sensor (S-SMC-M005; Onset, Kyoto, Japan). The measurement interval for SWC was 15 min. After drought hardening, irrigation was ceased in all treatments to verify the efficacy of drought hardening on physiological traits during consecutive dry days ("consecutive dry treatment"). Additionally, three seedlings were irrigated on day 8 of the consecutive dry treatment to verify the refilling of embolized tracheids in the basal stem xylem (see below for details).

Leaf Gas Exchange Traits during Drought Hardening and Consecutive Dry Treatment
Leaf gas exchange traits were measured in the morning on days 1, 16, and 31 of drought hardening after at least 30 min when the seedlings were exposed to full irradiation. The following traits were measured using a portable open gas exchange system (LI-6400; Li-Cor Inc., Lincoln, NE) with a leaf chamber fluorometer (LI-6400-40; Li-Cor Inc.): Net assimilation (A), water vapor stomatal conductance (g s ), and the transpiration rate (E). All measurements were conducted under a photosynthetic photon flux (PPF) of 1200 µmo m −2 s −1 using red-blue LED lamps, at a leaf temperature of 24.8 • C to 27.1 • C, and a CO 2 concentration of 400 µmol mol −1 in the inlet gas stream of the leaf chamber. To measure the response of g s to leaf-to-air vapor pressure deficit (VPD), VPD inside the leaf chamber was adjusted four points from 0.8 to 3.3 kPa, by gradually increasing the volume of dry air and by adding Drierite (97% CaSO 4 , 3% CaCl 2 ) to the inlet gas stream. After at least 15 min from the adjustment of VPD, g s was measured. As g s decreased exponentially with increasing leaf-to-air VPD, an analysis of stomatal sensitivity was carried out according to Oren et al. [24], where m is the stomatal sensitivity to VPD (mmol m −2 s −1 ln (kPa) −1 ) and b is the reference stomatal conductance at a leaf-to-air VPD of 1 kPa. The value of m was determined by extrapolation and curve fitting using least squares regression analysis. Before turning on the lights, all the pots, excluding the above-ground seedling, were covered in a plastic bag to prevent evaporation from the soil surface and the weight of each pot was measured (W start ). After turning off the lights (following exposure to light for 12 h), the pots were re-measured (W end ) and the plastic bags were removed to minimize the effect on root growth. Whole-plant transpiration(E whole ) was then calculated as follows, where W shoot is the shoot dry weight.

Growth and Biomass Allocation during Drought Hardening
Tree height and the diameter of the basal stem region were measured at both the start and end of drought hardening to determine the growth rate. Three seedlings per drought-hardening treatment were also sampled on days 1, 8, and 13 of consecutive dry treatment. Seedlings were then divided into leaf, stem, coarse-root (>2 mm in diameter), and fine-root (<2 mm in diameter) samples and dried at 70 • C for at least 3 days to determine the dry mass. Shoot and fine-root biomass, the ratio of shoot to whole-plant dry weight [shoot weight ratio (SWR)] and the ratio of fine-root to shoot dry weight (fine-root: shoot ratio) at the end day of drought hardening (start of consecutive dry treatment) were then determined.

Cryo-SEM Observations of Xylem and Cambium in the Basal Stem Region
To visualize the water distribution in the xylem tracheids and cambium under progressive drought during consecutive dry treatment, cryo-SEM observations of the xylem and cambium in the basal stem region were conducted on days 1, 8, and 13 of the consecutive dry treatment with one frozen basal stem per drought-hardening treatment per sampling day. Additionally, the seedlings that were irrigated on day 8 of consecutive dry treatment were also sampled on day 13 to verify recovery of xylem embolism. All seedlings were placed in a black plastic bag for at least 2 h to prevent transpiration [25] before the top shoot region was removed under running water. The top shoot region was submerged in water for at least 30 min to release negative pressure in the xylem conduits, since non-negligible artifacts have been reported under negative pressure [26,27]. Immediately after sampling, the main stems were submerged in liquid nitrogen for at least 1 min until the stem and water froze. They were stored in a deep freezer at −80 • C until cryo-SEM observation. The transverse face of the frozen sample was cleanly cut using a cryostat microtome (Cryo Star NX70, Thermo Fisher, Shanghai, China) then transferred to the cold stage of the cryo-SEM (JSM-6510 attached to a cryo-SEM unit; JEOL Ltd., Tokyo, Japan) and subjected to etching at −90 • C for 3-5 min. The specimens were then observed using cryo-SEM under an accelerating voltage of 3 kV and a cold stage temperature of approximately −130 • C.

Statistical Analysis
For statistical analysis of the relationships between the different drought-hardening treatments and the physiological and morphological traits in the C. japonica seedlings, generalized linear mixed models (GLMM) were developed and the best fit model was selected based on Akaike Information Criterion (AIC) values. All models were developed with a GLMM, assuming Gaussian error distribution with identity link function. The response variables were the physiological and morphological traits (numeric class), the explanatory variable was the different irrigation frequency during drought hardening (factor class), and the random effect was the holder of 200-cc single-cavity pots used in the phytotron chamber. Physiological traits were the leaf gas exchange parameters (A, g s , E, and E whole ) and morphological traits were the growth traits (growth rates of height and basal stem diameter during 1 month of drought hardening, SWR, and fine-root: shoot ratio). Five models were constructed using dummy explanatory variables as follows: Model 1 [y = 1 (null model)] indicated no difference between irrigation frequencies during drought hardening in any given response variable. Model 2 [y = (C treatment + M treatment) + S treatment − 1] indicated a difference in the response variable under the S treatment compared with the C and M treatments. In this model, (C treatment + M treatment) and (S treatment) were used as the dummy variables "ab" and "c", respectively. Model 3 [y = (C treatment + S treatment) + (M treatment) − 1] indicated a difference in the response variable under the M treatment compared with the C and S treatments, with (C treatment + S treatment) and (M treatment) used as the dummy variables "ac", and "b", respectively. Model 4 [y = (C treatment) + (M treatment + S treatment) − 1] indicated a difference in the response variable under the C treatment compared with the M and S treatments, with (C treatment) and (M treatment + S treatment) used as the dummy variables "a" and "bc", respectively. Model 5 [y = (C treatment) + (M treatment) + (S treatment) − 1] indicated a difference in the response variable between drought hardening treatments. In this model, (C treatment), (M treatment), and (S treatment) were used as the dummy variables "a", "b", and "c", respectively.
During consecutive dry treatment, mean E whole was also compared between measurement days in the three drought-hardening treatments using one-way repeated measures ANOVA followed by Dunnett's test at a significance level of α = 0.05, for multiple comparisons between the first day and subsequent measurement days. All statistical analyses were conducted using 'R' (version 2.15.1; Development Core Team), with the "lme4" and "multcomp" packages for model selection and multiple comparisons, respectively.  During consecutive dry treatment, mean Ewhole was also compared between measurement days in the three drought-hardening treatments using one-way repeated measures ANOVA followed by Dunnett's test at a significance level of α = 0.05, for multiple comparisons between the first day and subsequent measurement days. All statistical analyses were conducted using 'R' (version 2.15.1; Development Core Team), with the "lme4" and "multcomp" packages for model selection and multiple comparisons, respectively.

Volumetric Soil Water Content during Drought Hardening and Consecutive Dry Treatment
During drought hardening, the mean SWC was lower under S treatment [0.66 ± 0.22 (mean ± 1 standard deviation (SD)) m 3

Leaf Gas Exchange Traits during Drought Hardening and Consecutive Dry Treatment
The results of model selection in leaf gas exchange traits during drought hardening are shown in Table 1 and Appendix Table A1. There were no significant differences in A and gs on the first day of drought hardening between drought-hardening treatments. Although, a decrease was observed

Leaf Gas Exchange Traits during Drought Hardening and Consecutive Dry Treatment
The results of model selection in leaf gas exchange traits during drought hardening are shown in Tables 1 and A1. There were no significant differences in A and g s on the first day of drought hardening between drought-hardening treatments. Although, a decrease was observed under the S treatment compared with the C and M treatment on day 16 (the mid-point of drought hardening) and day 31 (the end day of drought hardening) (Figure 2a,b). In contrast, E and m were lower under the M treatment than the C and S treatments on the first day of drought hardening ( Figure 2c,d), while at the end of drought hardening, E was highest under the C treatment followed by the M treatment ( Figure 2c). There were no significant differences in m between treatments on the final day of drought hardening ( Figure 2d). E whole was higher under the S treatment than the C and M treatment on the first day of drought hardening, while at the mid-point to the end of drought hardening, E whole was the lowest under the S treatment followed by the M treatment (Figure 2e). under the S treatment compared with the C and M treatment on day 16 (the mid-point of drought hardening) and day 31 (the end day of drought hardening) (Figure 2a,b). In contrast, E and m were lower under the M treatment than the C and S treatments on the first day of drought hardening ( Figure 2c,d), while at the end of drought hardening, E was highest under the C treatment followed by the M treatment ( Figure 2c). There were no significant differences in m between treatments on the final day of drought hardening ( Figure 2d). Ewhole was higher under the S treatment than the C and M treatment on the first day of drought hardening, while at the mid-point to the end of drought hardening, Ewhole was the lowest under the S treatment followed by the M treatment ( Figure 2e). During consecutive dry treatment, the effect of the day of treatment on Ewhole was significant for all drought-hardening treatments (Appendix Table A2). Ewhole remained high for 4 days under the S treatment, while Ewhole remained high for 3 days under the M treatment, and under the C treatment, Ewhole remained high for only 2 days (Figure 3). On day 13 of consecutive dry treatment, Ewhole decreased to a low value in all drought-hardening treatments.
. were compared between drought-hardening treatments [control (C), mild drought (M), and severe drought (S)] during drought hardening. Measurements of leaf gas exchange traits were conducted in the morning on days 1, 16, and 31 of drought hardening after full irrigation. Boxplots show the median (bar), interquartile range (box), 5th and 95th percentiles (whiskers), and outliers (points) (n of A, gs, m, and E were 4, n of Ewhole was 12). For statistical analysis in the relationships between different drought-hardening treatments and the physiological traits in C. japonica seedlings, we developed generalized liner mixed models (GLMM) and then we selected the best fit model between these models based on AIC (Akaike Information Criterion) values. Different letters indicate statistically significant differences (see methods for derails of statistical analysis). "n.s." means no statistically significant differences.  , and outliers (points) (n of A, g s , m, and E were 4, n of E whole was 12). For statistical analysis in the relationships between different drought-hardening treatments and the physiological traits in C. japonica seedlings, we developed generalized liner mixed models (GLMM) and then we selected the best fit model between these models based on AIC (Akaike Information Criterion) values. Different letters indicate statistically significant differences (see methods for derails of statistical analysis). "n.s." means no statistically significant differences. Table 1. Estimated coefficients, standard error, and AIC values of the selected models in gas exchange traits in shoot. We developed generalized liner mixed models (GLMM) and then we selected the best fit model between these models based on AIC (Akaike Information Criterion) values. In these models, the response variables are leaf gas exchange (A, g s , E, m, and E whole ) (numeric class), the explanatory variable is the different irrigation frequency during drought hardening (factor class) and random effect is the holder of 200-cc single-cavity pots used in the phytotron chamber. In the column of selected model and response variables, C, M, and S means C treatment, M treatment, and S treatment during drought hardening. During consecutive dry treatment, the effect of the day of treatment on E whole was significant for all drought-hardening treatments (Table A2). E whole remained high for 4 days under the S treatment, while E whole remained high for 3 days under the M treatment, and under the C treatment, E whole remained high for only 2 days (Figure 3). On day 13 of consecutive dry treatment, E whole decreased to a low value in all drought-hardening treatments.  , and outliers (points) (n of Ewhole was nine from one to eight day and three from eight to thirteen day during consecutive dry treatment). The means of Ewhole were compared for each measurement days in all drought-hardening treatments using one-way repeated measures ANOVA followed by Dunnett's test at a significance level of α = 0.05, for multiple comparisons between the first day and subsequent measurement days. Asterisk indicate statistically significant differences from 1st day.

Growth and Biomass Allocation during Drought Hardening
The results of model selection in terms of growth and biomass allocation during drought hardening are shown in Table 2 and Appendix Table A3. During 1 month of drought hardening, the height growth rate, shoot dry weight, and SWR remained unchanged between drought-hardening treatments (Figure 4a,c,e), whereas diameter growth of the basal stem region was lower under the S treatment than the C and M treatments (Figure 4b). In addition, fine-root dry weight and the fineroot:shoot ratio were higher under the S treatment than the C and M treatments (Figure 4d,f). These findings suggest that growth and biomass allocation changed from diameter growth of the basal stem region to growth of the fine-roots during severe drought hardening. , and outliers (points) (n of E whole was nine from one to eight day and three from eight to thirteen day during consecutive dry treatment). The means of E whole were compared for each measurement days in all drought-hardening treatments using one-way repeated measures ANOVA followed by Dunnett's test at a significance level of α = 0.05, for multiple comparisons between the first day and subsequent measurement days. Asterisk indicate statistically significant differences from 1st day.

Growth and Biomass Allocation during Drought Hardening
The results of model selection in terms of growth and biomass allocation during drought hardening are shown in Tables 2 and A3. During 1 month of drought hardening, the height growth rate, shoot dry weight, and SWR remained unchanged between drought-hardening treatments (Figure 4a,c,e), whereas diameter growth of the basal stem region was lower under the S treatment than the C and M treatments (Figure 4b). In addition, fine-root dry weight and the fine-root:shoot ratio were higher under the S treatment than the C and M treatments (Figure 4d,f). These findings suggest that growth and biomass allocation changed from diameter growth of the basal stem region to growth of the fine-roots during severe drought hardening.

Cryo-SEM Observations of Xylem and Cambium in Basal Stem Region during Consecutive Dry Treatment
No embolism of xylem tracheids was observed under any of the drought hardening treatments on day 1 of consecutive dry treatment (Figure 5a). However, on day 8, embolism was shown under the C treatment, even in current-year xylem, while proportions of embolized xylem was small under the M and S treatments (Figure 5b). On day 13, most tracheids had lost water under the C and M treatments, while approximately half of tracheids in current-year xylem were filled with water under the S treatment (Figure 5c). In addition, cambium cells remained unchanged under the S treatment on day 13 of consecutive dry treatments, whereas cambium cell deformation was observed under the C and M treatments ( Figure 6).
Under the S treatment, comparisons of proportions of non-embolized xylem in the basal stem region before and after irrigation treatment on day 8 of consecutive dry treatment revealed no clear recovery of the embolism area even after re-watering ( Figure 7). This trend was also observed under the C and M treatments (data not shown). Table 2. Estimated coefficients, standard error, and AIC values of the selected models in growth, shoot and fine-root dry weight, the ratio of shoot dry weight to whole-plant dry weight (SWR), and the ratio of fine-root dry weight to shoot dry weight (fine-root:shoot ratio). We developed generalized liner mixed models (GLMM) and then selected the best fit between these models based on AIC (Akaike Information Criterion) values. In these models, the response variables are growth traits (the rate of height and basal diameter, shoot and fine-root dry weight, SWR, and fine-root:shoot ratio) (numeric class), the explanatory variable is the different irrigation frequency during drought hardening (factor class) and random effect is the holder of 200-cc single-cavity pots used in the phytotron chamber. In the column of selected model and response variables, C, M, and S means C treatment, M treatment, and S treatment during drought-hardening.   ; e), and the ratio of fine-root dry weight to shoot dry weight (fine-root:shoot ratio; f) were compared between drought-hardening treatment [control (C), mild drought (M), and severe drought (S)] during drought hardening. Boxplots show the median (bar), interquartile range (box), 5th and 95th percentiles (whiskers), and outliers (points) (n of height and diameter growth were 12. n of shoot and fine-root biomass,SWR, and fine-root:shoot ratio were 3). For statistical analysis in the relationships between different drought-hardening treatments and the morphological traits in C. japonica seedlings, we developed generalized liner mixed models (GLMM) and then we selected the best fit model between these models based on AIC (Akaike Information Criterion) values. Different letters indicate statistically significant differences (see methods for derails of statistical analysis). "n.s." means no statistically significant differences.

Cryo-SEM Observations of Xylem and Cambium in Basal Stem Region during Consecutive Dry Treatment
No embolism of xylem tracheids was observed under any of the drought hardening treatments  . Height (a) and diameter growth (b) in basal stem region, shoot (c) and fine-root biomass (d), the ratio of shoot dry weight to whole-plant dry weight (SWR; e), and the ratio of fine-root dry weight to shoot dry weight (fine-root:shoot ratio; f) were compared between drought-hardening treatment [control (C), mild drought (M), and severe drought (S)] during drought hardening. Boxplots show the median (bar), interquartile range (box), 5th and 95th percentiles (whiskers), and outliers (points) (n of height and diameter growth were 12. n of shoot and fine-root biomass, SWR, and fine-root:shoot ratio were 3). For statistical analysis in the relationships between different drought-hardening treatments and the morphological traits in C. japonica seedlings, we developed generalized liner mixed models (GLMM) and then we selected the best fit model between these models based on AIC (Akaike Information Criterion) values. Different letters indicate statistically significant differences (see methods for derails of statistical analysis). "n.s." means no statistically significant differences.     Prior to the consecutive dry treatment, the seedlings were exposed to the severe drought-hardening treatment (S). Notice embolized tracheids are characterized by dark lumina. Scale bars indicate 100 µm.

Discussion
One month of drought hardening enhanced water conservation under consecutive dry treatment without inducing leaf shedding in C. japonica seedlings. Mild and severe drought hardening (M and S treatment, respectively) promoted conserved water use through a decrease in whole-plant transpiration (E whole ) (physiological response) even after full irrigation during drought -hardening treatment (Figure 2e). This tendency was particularly evident under S treatment, whereas E whole for the C and M treatments on days 1 and 16 of the drought-hardening treatment could be underestimated because shoot biomass may increase, considering the results for diameter growth (Figure 4b). Therefore, our methodological approach needs to be interpreted with attention regarding the normalization of E whole . In addition, S treatment promoted a shift in carbon allocation from diameter growth of the basal stem region to fine-root growth (morphological response) (Figure 4b,d,f). These physiological and morphological responses to drought hardening probably contributed to the maintenance of both E whole ( Figure 3) and proportions of non-embolized xylem (Figure 5b,c) during consecutive dry treatment. E whole and proportions of non-embolized xylem were maintained longer under S treatment than M treatment (Figures 3 and 5). Considering that the seedlings were subjected to the drought-hardening treatment after the cessation of xylem development, the maintenance of E whole and proportions of non-embolized xylem during consecutive dry treatment may reflect conserved water use rather than the properties of xylem such as pit ultrastructure. Therefore, our results suggest that the effect of drought hardening on physiological and morphological traits depends on the frequency of irrigation of drought hardening.
The decrease in E whole as a result of drought hardening was inconsistent with the results of a previous study of acclimation to drought in 1.5-month-old C. japonica seedlings [21]. Nagakura et al. [21] found that the whole-plant transpiration rate remains unchanged, compared with control seedlings, following acclimation to dry conditions (drought hardening) for 12 weeks. In this previous study, watering was performed when the soil water content reached 0.30 m 3 m −3 , similar to the drought hardening M treatment in our study (Figure 1). We used three-year-old container-grown seedlings, and E whole under the M treatment was already lower than that under the C treatment on day 16 of drought hardening, even after full irrigation (Figure 2e). Therefore, the effect of drought hardening on conserved water showed a decrease in E whole , which may depend on the age and/or size of the seedling.
Our findings also suggest that our drought-hardening treatment has the potential to shorten the duration of drought hardening prior to planting for afforestation. Nagakura et al. [20] found that 19.5 weeks of drought hardening, which was 15.5 weeks longer than that in our study, decreased E whole in C. japonica seedlings. Additionally, for seedlings of Douglas fir, lodgepole pine, and white spruce, the duration of drought hardening is about 10 weeks [7], which is six weeks longer than in our study. The methodological differences between these previous studies and our study, included differences in the tree species and tree body size, pot size, irrigation interval, and soil type, which are related to the soil water holding capacity. Although, the age and size of C. japonica seedlings, soil types, and tree species considered in our study were limited, our drought-hardening treatment is probably applicable to seedlings up to the standard for height and basal stem diameter in C. japonica for afforestation in Japan.
Maintenance of xylem transport capacity is important for survival under consecutive drought in seedlings [28]. Embolized tracheids in the xylem are related to hydraulic failure, which leads to rapid or gradual dysfunction of both physiological and metabolic traits, including decreases in photosynthesis and respiration [29][30][31], which in turn leads to death if xylem embolism is not repaired. The repair of xylem embolism has been explained by both refilling of the xylem, which results in short-term repair (within hours or days) [32], and regrowth of new xylem, which provides long-term repair [33]. We did not detect the refilling of embolized tracheids five days after full irrigation during consecutive dry treatment in C. japonica seedlings (Figure 7), which is consistent with a previous study in the Cupressaceae family [34]. These findings suggest that species in the Cupressaceae family, including C. japonica, need to produce new xylem to recover xylem transport capacity after drought-induced hydraulic failure, due to an inability to refill their embolized tracheids. Moreover, Abe et al. [34] found that both the number of tracheids and mean radial diameter of tracheids decreased due to cambium cell distortion under water deficits. This report suggests that maintaining the cambium status is important for survival under drought conditions [33]. Our data showed that drought hardening contributes to the maintenance of the cambium status (Figure 6), and this could be caused by a reduction in E whole (Figure 3). Therefore, drought hardening contributes to the maintenance of cell differentiation related to the production of hydraulically functional xylem, thereby increasing drought tolerance under prolonged drought.

Conclusions
Drought hardening promoted water use conservation, which is related to the maintenance of xylem transport capacity and the cambium status during consecutive dry treatment. This is expected to increase the tolerance and survival of C. japonica seedlings under prolonged drought. Our results suggest that the effect of drought hardening on water use depends on stomatal regulation under drought conditions, but we cannot exclude impacts of other factors, such as modifications in root morphology (fine roots and root hairs), aquaporin expressions (both in roots and leaves), osmotic adjustment, etc. As our findings are based only on C. japonica seedlings, and are focused on stomatal regulation, we will examine the efficacy of drought hardening with a focus on variation in stomatal control and other factors, such as aquaporin expressions, under drought conditions in a wider range of tree species (including isohydric and anisohydric tree species) [35] to validate the generalizability of our findings.
Due to the high genetic divergence between elite trees of C. japonica in Japan, differences in the physiological and morphological responses to environmental conditions associated with genetic variation may limit the scope of our findings. Hiraoka et al. [36] reported that annual height increments differ among 12 genotypes of C. japonica in a common garden, suggesting that the efficacy of drought hardening is influenced by genotype. Further studies, focusing on the roles of genetic factors, are needed to establish effective methods for drought hardening in different genotypes in C. japonica.
Oguchi and Fujibe [6] projected that consecutive dry days increase by around 1.0 day per year over the wide area of Japan. Moreover, the average length of consecutive dry days is around 20 days, with annual fluctuations [6]. These projections indicate that seedlings for afforestation are expected to be exposed to unexperienced severe drought in the future. We showed that the proportions of non-embolized xylem and the phloem status, under the S treatment, were maintained for 13 days without irrigation, which is at least five days longer the duration for the C and M treatments ( Figures 5  and 6, respectively). However, we only examined the effectiveness of different levels of drought hardening before planting in C. japonica seedlings. Therefore, we intend to examine the effect of drought hardening on drought tolerance using planted seedlings after drought hardening for validation of our findings. Table A1. AIC (Akaike Information Criterion) values of all generalized liner mixed models (GLMM) in gas exchange traits in shoot. In these models, the response variables are leaf gas exchange (A, g s , E, m, and E whole ) (numeric class), the explanatory variable is the different irrigation frequency during drought hardening (factor class) and random effect is the holder of 200-cc single-cavity pots used in the phytotron chamber. In the column of selected model and response variables, C, M, and S means C treatment, M treatment, and S treatment during drought hardening. Gray colored cells indicate the best fit models based on AIC values.  Table A1. Cont.   Table A3. AIC (Akaike Information Criterion) values of all generalized liner mixed models (GLMM) in growth, shoot and fine-root dry weight, the ratio of shoot dry weight to whole-plant dry weight (SWR), and the ratio of fine-root dry weight to shoot dry weight (fine-root:shoot ratio). In these models, the response variables are growth traits (the rate of height and basal diameter, shoot and fine-root dry weight, SWR, and fine-root:shoot ratio) (numeric class), the explanatory variable is the different irrigation frequency during drought hardening (factor class) and the random effect is the holder of 200-cc single-cavity pots used in the phytotron chamber. In the column of selected model and response variables, C, M, and S means C treatment, M treatment, and S treatment during drought-hardening treatments. Gray colored cells indicate the best fit models based on AIC values.  Table A3. Cont.