Identification of the cold-responsive TPP gene of A. arguta and its subcellular localization
From the transcriptome data of A. arguta cv. Kuilv (KL), we found 39 trehalose-metabolism-related genes including 21 TPSs, 17 TPPs, and one TRE. Comparative analysis of their expression levels between normal (CK) and low temperature (LT) conditions showed that among these genes, TPP#3, #5, #8, and #9 expressed differentially in response to cold stress (Fig. 1a). Since A. arguta (KL) is a tetraploid plant (2n = 4×=116) (Liu et al. 2011), these cold-responsive TPPs might be alleles. Therefore, we cloned all cDNAs of TPP#3, #5, #8, and #9, and found more than 22 different sequences (Fig. S1a). Multiple alignment and phylogenetic tree of these sequences showed that the cloned TPP genes were largely divided into three groups regardless of the primers used for PCR amplification (Fig. S1b). The identities between groups were over 84.6% (Fig. S1c), suggesting that these cold-responsive TPP genes of each group were alleles or duplicate genes, as expected. Among these genes, we selected the Group I TPP gene as the purpose cold-responsive TPP gene, which has the highest expression difference between CK and LT conditions (TPP#5; Appendix S1). This cold-responsive TPP gene of A. arguta had the highest homology with TPPA (trehalose-6-phosphate phosphatase A) of Populus trichocarpa and Arabidopsis thaliana (Fig. 1b). From this, we termed the cold-responsive A. arguta TPP gene as AaTPPA. The similarity of AaTPPA with AtTPPA, PtrTPPA, and OtsB (the E. coli TPP gene) was 66.8%, 73.3%, and 21.4%, respectively (Fig. S2).
The TPP gene family belongs to the HAD (L-2-haloacid dehalogenase) superfamily of phosphatases and other hydrolases, and the proteins of this superfamily are characterized by having three conserved motifs (HAD motif I, II, and III) that form the catalytic active site (Fieulaine et al. 2005). AaTPPA also had all the HAD motifs, considering that it could have the enzymatic activity of trehalose-6-phosphate phosphatase (Fig. S2). To determine its enzymatic activity, the His6-tag-linked AaTPPA protein was expressed in E. coli, and purified by Ni-affinity chromatography (Fig. S3). The His6-tag-linked OtsB protein was used as a positive control. The recombinant AaTPPA protein was mainly expressed in the form of an inclusion body in the E. coli host cells, while the recombinant OtsB protein was present in the cytoplasmic soluble fraction, suggesting that AaTPPA is largely different from OtsB in physiochemical properties. The enzymatic activity of AaTPPA was more than twice that of OtsB at 9–25°C and decreased largely above 35°C (Fig. 1d).
AaTPPA was predicted to have a monopartite nuclear localization signal (NLS) of RKKAK(73–77) (Fig. S2a). In addition, the predicted 3D structure of AaTPPA showed that it has a catalytic domain similar to OtsB, but AaTPPA has an additional N-terminal domain of an unknown function which is common in plant TPP genes (Fichtner and Lunn 2021) (Fig. 1c). The analysis of subcellular localization in onion epidermal cells and tobacco leaf epidermal cells showed that both AaTPPA and OtsB were present in the cell nucleus and the whole cytoplasm (Fig. 1e, f). This result is consistent with a previous report that AtTPPA, AtTPPB, AtTPPC, AtTPPF, and AtTPPH are present in the cytosol and nucleus (Krasensky et al. 2014).
AaTPPA expression is induced in response to mild low-temperature
A. arguta acquires cold acclimatization as the weather temperature gradually decreases in autumn (Van Labeke et al. 2015). Thus, it was necessary to determine at what temperature the expression of the cold-responsive AaTPPA was induced. As a result, the expression level of AaTPPA in mature leaves gradually increased under the mild low-temperature of 15°C, and it was 3.5 times higher at 6 h than before the low-temperature treatment. AaTPPA expression decreased gradually over time at 4°C and 9°C (Fig. 2a).
A. arguta has a strong overwintering ability, but the developing young tissues are very susceptible to frost damage in spring (Marosz 2009). So, we investigated the low-temperature responsiveness of AaTPPA at different leaf growth stages. At the sprouting bud stage, AaTPPA did not respond to mild low temperature, and at the young leaf stage, its expression level increased about twice only after 6 h at 15°C (Fig. 2b). In addition, AaTPPA expression was not induced at all in the sprouting buds and the young leaves under 4°C and 9°C treatments.
Based on the above results, we investigated changes in AaTPPA expression levels in mature leaves of A. arguta grown in a garden as the weather temperature dropped in autumn. Consistent with the indoor experiment, the expression levels of AaTPPA in mature leaves gradually increased until August 23, 2022 (Harbin, China), when the daily low-temperature elapsed 13°C, and then decreased as the temperature dropped below 10°C (Fig. 2d). The weather temperatures below 10oC suppressed the growth of shoot apical meristems, and accelerated leaf senescence (Fig. 2c). Leaves began to fall on September 25, the deadline for sampling. This result indicated that AaTPPA was transiently induced in mature leaves in response to temporary mild low temperature, and its expression was not maintained at the elevated level during the entire autumn, suggesting that temporary expression of AaTPPA in leaf tissues might not make a critical contribution to the acquisition of overwintering ability of A. arguta.
Therefore, to determine the main acting sites of AaTPPA for acquiring cold acclimatization of A. arguta, the plant's roots, stem cambium, lateral buds, and mature leaves were collected on September 25 in 2022, and AaTPPA expression levels were analyzed (Fig. 2e). The result showed that AaTPPA expression levels in lateral buds, roots, and stem cambium were much higher than that in leaves. Since lateral buds, roots, and stem cambium are dormant tissues in winter, the high expression levels of AaTPPA in these tissues implied that this gene might play an essential role in the cold acclimatization of A. arguta.
Developmental and morphological characteristics of AaTPPA-overexpressing Arabidopsis
Among the transgenic A. thaliana plants, #4 and #7 lines with high expression levels of AaTPPA were selected for AaTPPA-overexpressing (OE) lines (Fig. S4a). We also constructed A. thaliana transformed with OtsB which only has a catalytic domain, unlike the plant TPP enzymes, and OtsB-OE #6 and #7 lines were selected as positive controls (Fig. S4b). Seeds of the selected OE lines were all hygromycin-resistant at T3 generation, confirming that these lines were homozygous.
First, we measured trehalose and Tre6P levels in AaTPPA and OtsB-OE A. thaliana lines. Surprisingly, trehalose levels in AaTPPA-OE lines were 8 ~ 11 times higher than that of the wild type (WT), whereas no significant changes were observed in OtsB-OE lines (Fig. 3a). T6P levels in A. thaliana were too low to distinguish significant differences between WT and the OE lines, but the result showed a tendency that T6P levels of OtsB-OE lines were lower than those of WT and AaTPPA-OE lines (Fig. S4c). To assay changes in Tre6P levels, we further investigated the expression levels of several SnRK1 marker genes (Fig. S4e). As a result, changes in the expression levels of SnRK1 marker genes in OtsB-OE lines were consistent with the previous report (Zhang et al. 2009), but in AaTPPA-OE lines, there was no significant difference from WT except for the expression of AtDXP. This suggested that Tre6P levels of AaTPPA-OE A. thaliana did not change significantly enough to affect the SnRK1 activity, unlike that of OtsB-OE lines.
Next, we compared the developmental and morphological characteristics between AaTPPA and OtsB-OE A. thaliana. The seeds of all the OE lines were smaller and misshapen compared to WT (Fig. 3b, c). Especially, the seeds of AaTPPA-OE lines were smaller even than those of OtsB-OE lines. However, the small seed size did not affect the germination rate, suggesting that overexpression of the TPP genes only affected seed filling but not embryogenesis (Fig. S5a). Seedlings of the OE lines grown in MS media were also smaller in size and had short primary root lengths at first (Fig. S5b, c; Fig. 3d, e), but afterward, there was no significant difference in growth compared to WT (Fig. S5d), demonstrating that the initial growth disturbance was only due to a lack of nutrients in the seeds of the TPP-OE lines. AaTPPA overexpression rather seemed to enhance vegetative growth in the later developmental stages (Fig. 3f).
AaTPPA-OE lines had elongated rosette leaves like WT, but OtsB-OE lines had broad-shaped leaves (Fig. S5d; Fig. 3g). In addition, chlorophyll contents of rosette leaves were significantly lower in OtsB-OE lines, but there was no difference between AaTPPA-OE lines and WT (Fig. S5e). As in previous studies, bolting periods of all the OE lines were delayed than WT (Fig. 3f), and even at the time when all the rosette leaves of WT had wilted, the rosette leaves of the OE lines remained vivid (Fig. 3h). Interestingly, AaTPPA-OE lines had more branches and siliques than WT, but OtsB-OE lines had significantly fewer branches and siliques (Fig. 3i; Fig. S5g). No significant difference was observed in heights between the OE lines and WT (Fig. S5f).
Overall, AaTPPA-OE A. thaliana had much higher trehalose levels and showed distinct developmental and morphological characteristics. Like OtsB-OE lines, AaTPPA-OE lines had misshapen seeds and delayed reproductive growth and senescence. However, unlike OtsB, AaTPPA overexpression did not affect the leaf shape and chlorophyll content. Furthermore, AaTPPA-OE lines had more branches and siliques compared to WT and OtsB-OE lines.
Heterogeneous overexpression of AaTPPA enhances the freezing tolerance of Arabidopsis
First, we investigated the freezing tolerance of the OE A. thaliana lines at the seedling stage. The survival rate of AaTPPA-OE #4 line was significantly higher than that of WT after freezing treatment at -10°C for 0.5 and 1 h. On the contrary, OtsB-OE lines showed a relatively low survival rate compared to WT (Fig. S6a, b). Meanwhile, the ion leakages of both AaTPPA-OE and OtsB-OE seedlings were lower than that of WT after 1 h of freezing treatment, suggesting that the TPP overexpression had a protective effect on biological membranes (Fig. S6c).
Next, A. thaliana plants grown in pots were subjected to freezing stress at -10°C. The leaf wilting rates of AaTPPA-OE lines of 4 h treatment group were significantly lower than that of other lines (Fig. 4a, b). Surprisingly, after normal cultivation of this group for 7 d, the reviving rate of AaTPPA-OE #4 line was twice that of WT, while the reviving rates of OtsB-OE lines were significantly lower than that of WT (Fig. 4c).
Freezing stress-treated A. thaliana lines were stained with NBT to assay the levels of superoxide anion (O2−) (Fig. S6d, e). Under normal conditions, O2− levels of all the OE lines were significantly lower than that of WT, suggesting that the TPP overexpression increased ROS scavenging capacity. Under freezing conditions, O2− levels of all the lines decreased as freezing time increased, and its levels were maintained higher in all the OE lines than in WT. This result could be interpreted that photosystem I (PSI), the generator of O2− in chloroplasts, was destructed under extreme freezing conditions, but the TPP overexpression prevented PSI from freezing destruction. The fact that the plant photosystems were protected in the OE lines was also proved by chlorophyll fluorescence assay. After cold stress at -2°C for 2h, Fv/Fm values remained the highest in AaTPPA-OE lines, and the values of OtsB-OE lines were also higher than that of WT (Fig. 4d, e).
We also evaluated the degree of hydrogen peroxide (H2O2) accumulation using DAB staining (Fig. S6f, g). AaTPPA-OE lines tended to accumulate less H2O2 than WT and OtsB-OE lines under normal and freezing stress, but no significant difference was recognized. Moreover, the activities of SOD, POD, and CAT increased more in AaTPPA-OE lines than in WT and OtsB-OE lines under normal and cold stress conditions (Fig. S6h, i, j), and the MDA contents also increased more in AaTPPA-OE lines after cold stress (Fig. 4f), indicating that AaTPPA overexpression plays a positive role in activating the plant antioxidant system.
The above results implied that AaTPPA overexpression could have a protective effect on biological membranes and biomolecules by accumulating osmotic protectants. For proof of concept, the contents of sucrose, soluble sugars, and proline, which are the main osmotic protectants of plants, were measured, showing that under normal and cold stress conditions more sucrose and soluble sugars were accumulated in both AaTPPA and OtsB-OE lines than those in WT, and proline contents were higher only in AaTPPA-OE lines (Fig. 4g, h, i).
In summary, AaTPPA overexpression enhanced the cold tolerance of A. thaliana, but OtsB overexpression did not affect the cold tolerance obviously and was rather negative in some aspects.
Transient overexpression of AaTPPA in A. arguta leaves improves cold tolerance indices
Through the Agrobacterium tumefaciens-mediated transformation, target genes can be transiently overexpressed in plant leaf tissues. However, this method has low reproducibility, because the expression efficiency of the target gene is different each time depending on the growth state of the infiltrated leaves and the experimenter's technique. To improve the reproducibility of this transient transformation method, we created a new plant expression vector (GVloxP/DsRed) which contains a red fluorescence protein (DsRed) as a marker, allowing to monitor the degree of transfection in the infiltrated leaves by a non-invasive manner (Fig. S7). The overexpression vector, GVloxP/DsRed-AaTPPA, was prepared by replacing the CDS of the hygromycin resistance gene (HygR) of this vector with AaTPPA CDS (Fig. 5a).
Young leaves of A. arguta cv. Longcheng No.2 (L2) plantlets grown in pots were infiltrated with Agrobacterium containing GVloxP/DsRed-AaTPPA vector, and after 3 d, they were observed under a fluorescence stereomicroscope, and leaves with obvious red fluorescence were collected to determine the gene expression level (Fig. 5b, c). The result demonstrated that the AaTPPA expression level was incomparably higher in AaTPPA-OE group than in EV (the mock group only infiltrated with GVloxP/DsRed empty vector) and CK (the control group not transfected).
Next, A. arguta (L2) plantlets with red-fluorescent leaves that were transformed as described above were subjected to cold stress at -4°C for 2 h, and then several cold tolerance-related indices, such as ROS accumulation, ion leakage, MDA, soluble sugars, and proline contents, were analyzed. Interestingly, DAB staining showed that H2O2 accumulation in EV group was particularly higher, explaining that the leaf tissues might be damaged by the transfection process, or that the overexpression of marker genes, such as DsRed and HygR, might be harmful to the physiology of the plant (Fig. 5d, e). Nevertheless, H2O2 level of OE group was significantly lower than that of CK group after cold stress. In addition, after cold stress, the ion leakage and MDA contents of AaTPPA-OE leaves were significantly lower than those of CK and EV groups (Fig. 5f). Soluble sugars and proline contents were higher in AaTPPA-OE group even under normal conditions. Cold stress increased the contents of these osmotic protectants in all groups, but AaTPPA-OE group still maintained the highest level (Fig. 5f).
Cloning of the AaTPPA promoter and verification of its transcriptional activity
Since the genome sequence of A. arguta has not been published yet, the promoter sequence of AaTPPA had to be cloned to seek upstream transcription factors that affect the AaTPPA expression in response to low temperature. According to the genome data of A. chinensis cv. Red5, the homologous TPPA gene had an exceptionally long 5'-UTR (3,162 bp) upstream of the ATG translational start codon (Fig. S8h). Therefore, to obtain the promoter sequence of AaTPPA, we had to clone a region of at least 4,000 bp upstream of the ATG codon. By combining the conventional PCR and TAIL-PCR (thermal asymmetric interlaced PCR), a total region of 4,138 bp upstream of the ATG codon was cloned and sequenced (Fig. S8, Appendix S2). Alignment with the TPPA gene sequence of A. chinensis cv. Red5 showed that among the whole cloned sequence, the promoter region was approximately 1185 bp, and the remaining section (2953 bp) belonged to 5'-UTR (Fig. S9a).
In the previous study on the promoters of TPP genes, ~ 2000 bp sequence upstream to CDS was considered to be the promoter region (Wang et al. 2022). So, to determine whether the long 5'-UTR of AaTPPA affects the gene expression, the CaMV 35S promoter of the GUS gene in pCAMBIA1303 vector was substituted by four sections of the cloned 5'-UTR plus promoter region to construct the promoter activity checking vector (Fig. S9a, b). The AaTPPApro2 section contained a 1685 bp sequence from the ATG codon to the middle of intron 1, the AaTPPApro3 section contained a 2688 bp sequence including the almost intron 1, and the AaTPPApro9 included the whole 5'-UTR and promoter region of AaTPPA (4138 bp). The AaTPPApro8 section contained only the predicted promoter region (Fig. S9a). The GUS expression cassettes driven by the above four different sections were transformed into A. thaliana, and the screened transgenic T1 plants were verified by PCR test on the corresponding sections (Fig. S9c-j).
GUS gene driven by AaTPPApro9 was expressed in almost all tissues including roots, vascular tissues, mesophyll, stems and branches, sepals and petals, siliques, and developing seeds (Fig. 6a). In comparison, GUS driven by AaTPPApro8 had the same expression pattern as AaTPPApro9, except for young leaves and developing seeds. On the contrary, AaTPPApro2 and AaTPPApro3-driven GUS were not expressed in roots and were mainly expressed in vascular tissues. AaTPPApro3::GUS was also specifically expressed in developing seeds. The above results showed that the AaTPPApro8 sequence had sufficient promoter activity in almost all tissues except for developing seeds and young leaves.
Next, to determine whether the AaTPPApro8 section has cold-stress responsive characteristics, we analyzed the expression levels of the AaTPPApro8-driven GUS gene under normal and low temperature (9°C) conditions (Fig. 6b, c). The result showed that GUS staining was more intense in the low-temperature treatment group than in the normal group, indicating that the AaTPPApro8 sequence has a low-temperature responsive cis-acting element.
Taken together, it can be concluded that AaTPPA has a long 5'-UTR of about 3000 bp and that the AaTPPApro8 sequence of 1274 bp has a low-temperature responsive characteristic as well as constitutive promoter activity.
AaERF64 binds to the AaTPPA promoter and activates its transcription
We searched for cis-acting regulatory elements in the promoter region of 1185 bp upstream from the predicted transcriptional start site using the online software PlantCARE (Fig. 7a, Table S1). Based on the motif analysis, we adopted three sections which are clustered with functional motifs, and each fragment was inserted into pAbAi vector, respectively, to prepare bait vectors. Among the Y1HGold strains transformed with these bait vectors, only the AaTPPAproS3 introduced strain was completely suppressed on 300 ng/mL aureobasidin A (AbA) containing medium, and the others were insensitive to AbA because of their constitutive promoter activity. Therefore, proteins interacting with the AaTPPAproS3 section were screened from the A. arguta cDNA library using the yeast-one-hybridization (Y1H) method. As a result, three ethylene-responsive transcription factors (ERF) were identified as transcription factors. From the transcriptome data of A. arguta, the ERF gene (gene28680) which was responsive to low temperature was selected for further study (Table S2). Its interaction with the AaTPPA promoter was reconfirmed by using the Y1H assay (Fig. 7b). This ERF gene had the highest homology with ERF64 of Actinidia deliciosa (AdERF64) in NCBI's nucleotide database (Fig. S10a), so we called it AaERF64 (Appendix S3).
In plants, ERF genes constitute a large family of transcription factors which is a part of the AP2/ERF (APETALA2/ETHYLENE RESPONSE FACTOR) superfamily (Riechmann et al. 2000). Generally, proteins of AP2/ERF superfamily have an AP2 DNA-binding domain of approximately 60–70 amino acids (Sakuma et al. 2002). AP2/ERF superfamily is divided into 12 phylogenetic groups in A. thaliana, and some of them are further divided into several subgroups (Nakano et al. 2006). The phylogenetic tree with 122 ERF family members of A. thaliana showed that AaERF64 belongs to the V group of the ERF family, and has the closest relation to ERF#003 (At5g25190.1) of A. thaliana (Fig. S10a). This group is characterized to have one CMV-2 motif (Nakano et al. 2006). Bioinformatic analysis showed that both the AP2 DNA-binding domain and CMV-2 motif are well conserved in the AaERF64 protein (Fig. S10c). In addition, the N-terminal contains NLS peptide, indicating that this protein could be localized in the cell nucleus. The predicted 3D structure showed that the AaERF64 protein is composed of two domains including the AP2 DNA-binding domain (4–69 aa) and a poorly predicted C-terminal domain (70–169 aa). The C-terminal domain was thought to be a transcriptional activation domain due to containing the CMV-2 motif (Fig. S10b).
Subcellular localization assay of AaERF64 in tobacco leaf epidermal cells showed that AaERF64 was localized in the cell nucleus, satisfying the first requirement for acting as a transcription factor (Fig. 7c). Its transcriptional activation ability was examined by the Y2H self-transactivation test (Fig. 7d). The result demonstrated that AaERF64 had a transcriptional activation function and that the C-terminal domain (70–169 aa) was a transcriptional activation domain, as expected. In addition, the dual luciferase assay also demonstrated that AaERF64 could interact with an element of the AaTPPA promoter to activate the transcription of the downstream gene in plant cells (Fig. 9e). Furthermore, we transiently transfected A. arguta leaves with the GVloxP/DsRed-AaERF64 vector and proved that AaERF64 overexpression upregulates the expression of the endogenous AaTPPA in A. arguta genome (Fig. 8a). These results revealed that AaERF64 acts as a transcriptional activator of AaTPPA in A. arguta.
Next, we searched for the AaERF64-binding cis-element within the AaTPPA promoter sequence. All the ERF family proteins have the conservative AP2 DNA-binding domain, which has been known to specifically bind to GCC-box (Ohme-Takagi and Shinshi 1995). We found three GCC-box motifs in the AaTPPAproS3 section, which was used as a bait sequence for the AaERF64 screen (Fig. 7f). Three DNA fragments (EM1 ~ 3) containing each GCC-box motif were introduced into the pAbAi vector, respectively, and their binding potential with AaERF64 was examined by employing the Y1H assay. The result revealed that the EM2 fragment (22 bp) interacted with AaERF64. Furthermore, electrophoretic mobility shift assay (EMSA) proved that AaERF64 directly binds to the GCC-box motif (GCCGAC) of the EM2 sequence (Fig. 7g).
AaERF64 expression is induced through the ethylene signal transduction pathway
As shown above, AaTPPA expression levels were particularly high in overwintering tissues of A. arguta in autumn (Fig. 2e). This time, we investigated changes in the expression levels of AaTPPA and AaERF64 in these tissues according to the temperature decrease in autumn (Fig. 8b). Based on the previous experimental data, we selected three sampling times; the first time was the day when the weather temperature was the highest in summer, 2023, the second time was the day when the daily low-temperature dropped below 15°C, and the last time was the day when the daily low-temperature dropped below 10°C. At the last time, the shoot apical buds began to degenerate.
In mature leaves, stem cambium, and lateral buds, the expression levels of AaTPPA were positively correlated to those of AaERF64, supporting the prior result that AaERF64 overexpression induced upregulation of AaTPPA expression (Fig. 10a). Unexpectedly, in roots, no significant changes were observed in the expression levels of AaERF64 depending on the period, but the expression levels of AaTPPA gradually increased as the temperature became colder. This suggested that AaERF64 has no action in roots and that AaTPPA has a high expression level specifically in roots, and may be regulated by other pathways.
In mature leaves, consistent with the results in 2022 (Fig. 2d), the expression levels of AaERF64 and AaTPPA were the highest at the mild low temperature (daily high-temperature 24°C, daily low-temperature 14°C), and afterward, the expression of both genes decreased. However, AaERF64 expression in lateral buds and stem cambium continued to increase until the last day (September 15, 2023), and so did AaTPPA expression. This result suggested that under excessively low temperature conditions, leaves become withered, leading to the reduction of overall gene expression, whereas in lateral buds and stem cambium, expression of genes, especially the genes related to cold stress, continue to be activated even in late autumn.
Next, we conducted research to identify the signal transduction pathway that induces the expression of the AaERF64 transcription factor. According to previous studies, plant TPP genes were induced by several phytohormones, such as abscisic acid (ABA), jasmonic acid methyl ester (MeJA), and gibberellin (GA) (Ge et al. 2008; Lin et al. 2020; Huang et al. 2022b), and the cis-acting element analysis of the AaTPPA promoter predicted that there exists MeJA and ABA-responsive motifs (Fig. 7a, Table S1). For AaERF64 as an ethylene-responsive factor, we selected ethephon (the precursor of ethylene), MeJA, and ABA as the putative inducible phytohormones (Fig. S12c). In A. arguta leaves, AaTPPA was induced by ethylene and ABA, and AaERF64 was induced only by ethylene treatment. Induction of AaERF64 and AaTPPA expression by ethylene occurred at 1–2 d of treatment, whereas induction of AaTPPA expression by ABA peaked at 12 h and then decreased. Interestingly, ABA treatment rather suppressed AaERF64 expression. This suggested that the induction of AaTPPA expression by ABA occurs through a pathway different from AaERF64 and that the induction by ABA occurs transiently in leaves, which is reversible when the ABA signal disappears.
We further investigated changes in the expression levels of AaTPPA and AaER64 in the overwintering tissues of A. arguta by ethylene treatment (Fig. 8c). In the lateral buds and stem cambium, the expression levels of AaTPPA and AaER64 continued to increase until 15 d of treatment, reaching over 10-fold higher than before treatment. However, their expression levels in leaves did not increase further after 5 d of treatment. Soluble sugars and proline are the main osmotic protectants in plants (Ding et al. 2020). Undoubtedly, lateral buds, stem cambium, and roots of A. arguta also contained much soluble sugars and proline in late autumn (Fig. S12a, b). Ethylene treatment also increased the contents of soluble sugars and proline in the lateral buds and stem cambium more than in the leaves (Fig. 8c). However, in roots, AaTPPA and AaER64 expression did not respond to ethylene treatment, and there were no significant changes in the contents of soluble sugars and proline. This suggested that the increase in AaTPPA expression and the accumulation of soluble sugars and proline in roots in autumn might occur depending on other signaling pathways, maybe by ABA (Fig. 8b).