Saintpaulia (Saintpaulia ionantha), a popular indoor ornamental potted plant, is native to the highlands of Kenya and Tanzania where temperatures rarely fall below 4°C. Chilling injury caused during cultivation and transportation is a major commercial problem in Saintpaulia. In this study, we investigated chilling acclimation in Saintpaulia ‘Kilauea’. Plants grown at 20°C (14 h light/10 h dark) displayed rapid and severe chilling injury after 24-hour exposure to 4°C. However, chilling injury at 4°C could be dramatically reduced by pre-treating the plants at 10°C but not at 6°C. From whole genome analysis, 161 Ethylene-Responsive Factors (ERFs) were identified and classified into 12 clades according to existing reports. Among these ERFs, 43, 8, and 4 ERFs were upregulated at 12, 24, and 48 h after 10°C treatment, respectively. Most of these ERFs had GCC box and/or DRE/CRT core motifs-like sequences in their upstream regions. Finally, we compared the expression of ERFs between the treatments for 24 h at 10°C, an effective temperature for chilling acclimation, and at 6°C, an ineffective temperature. The results showed that the expressions of all six ERFs we investigated were increased by the 10°C treatment, but not or only barely increased by the 6°C treatment. This study suggests that Saintpaulia, a subtropical plant, also has the ability to acclimate to low temperatures and that ERF upregulation is involved in chilling acclimation.
Research Article
Comprehensive expression analysis of ERF transcription factors during chilling acclimation in Saintpaulia
https://doi.org/10.21203/rs.3.rs-3960025/v1
This work is licensed under a CC BY 4.0 License
You are reading this latest preprint version
chilling stress
chilling acclimation
Saintpaulia
ERF
CBF/DREB
Saintpaulia or African violet (Saintpaulia ionantha) is an ornamental plant that belongs to the Gesneriaceae family. It is popular as a potted flower for indoor use because of its high shade tolerance and year-round flowering under low light intensity. Saintpaulia is native to Kenya and Tanzania, where temperatures are relatively stable; therefore, it is sensitive to temperature change (Johansson 1978; Baatvik 1993). When exposed to low temperatures, Saintpaulia leaves turn brown and visible damage occurs even within a short period (approximately 12 h) (Larcher and Bodner 1980). For chilling sensitive plants, chilling is a significant cause of damage not only in production but also in transportation and interior design (Chen et al. 2004). Therefore, the cultivation and shipping of Saintpaulia under low temperatures during the winter season require extreme care.
Chilling injuries are caused by non-freezing temperatures (Levitt 1980). Most studies on cold acclimation have investigated responses to freezing stress. Cold acclimation occurs when plants are exposed to suboptimal low temperatures, which do not cause plant death, but enhance their cold tolerance (Levitt 1980; Gilmour et al. 1988; Thomashow 1999). Research has been conducted using model plants native to temperate, such as Arabidopsis thaliana, to investigate cold tolerance and acclimation by exposure to sublethal temperatures. A. thaliana dies at temperatures near freezing, but exposure to 4°C for some time improves its survival rate at − 9°C (Gilmour et al. 1988). Many plants native to tropical and subtropical regions, including Saintpaulia, suffer from chilling injury or die at low temperatures (~ 10°C). If chilling acclimation is also possible in Saintpaulia, it would allow for safe cultivation and commercial transportation during the winter season.
Genetic studies of cold acclimation have been conducted extensively in A. thaliana (Thomashow 1990; Xin and Browse 2000; Yamaguchi-Shinozaki and Shinozaki 2006). C-repeat binding factor /Dehydration responsive element binding 1 (CBF/DREB1) transcription factors, which belong to the APETALA2/Ethylene-Responsive Factor (AP2/ERF) family, are the most well-known genes that are upregulated during cold acclimation and are known to improve cold tolerance. AtCBF1, AtCBF2 and AtCBF3 of A. thaliana (also known as AtDREB1B, AtDREB1C, and AtDREB1A, respectively) bind to the CRT/DRE element in the promoter region of downstream genes, such as the cold-regulated (COR) genes, and contribute to the improvement of cold tolerance by regulating their transcription (Jaglo-Ottosen et al. 1998; Liu et al. 1998; Gilmour et al. 2000, 2004; Yamaguchi-Shinozaki and Shinozaki 2006). AtCBF2 inhibits AtCBF1/3, and cold tolerance is improved in cbf2 mutants of A. thaliana, indicating a complex regulatory mechanism (Novillo et al. 2004). CBF/DREB1 improves cold tolerance not only in A. thaliana but also in various other plants, such as rice (Oryza sativa), maize (Zea mays), and tomato (Solanum lycopersicum) (Dubouzet et al. 2003; Qin et al. 2004; Zhang et al. 2004). In addition, several ERFs other than CBF/DREB1 have been reported to contribute to the improvement of cold tolerance in diverse plant species, such as AtERF102/103/104/105 (A. thaliana), PtrERF9/108/109 (Poncirus trifoliata), MdERF1B, MdABI4 (Malus × domestica), and VaERF057, 092 (Vitis amurensis) (Sun et al. 2016; Bolt et al. 2017; Sun et al. 2019; Wang et al. 2019; Illgen et al. 2020; Khan et al. 2021; Wang et al. 2021; An et al. 2022; Zhang et al. 2022). Although various transcription factors other than ERFs are involved in cold acclimation, ERF transcription factors have a particularly central role (Guo et al. 2022).
On the other hand, few studies have been done on chilling acclimation of plants native to tropical and subtropical regions such as Saintpaulia. And, there are also few reports on whether ERF is involved in the chilling tolerance of subtropical and tropical native plants. In this study, we investigated two inquiries. First, we investigated whether chilling acclimation occurs in Saintpaulia. Second, we conducted the comprehensive expression analysis of ERFs based on the whole genome sequence and transcriptome analysis during sublethal temperature exposure in Saintpaulia. Finally, we discussed the involvement of ERFs in the chilling acclimation of Saintpaulia.
Plant materials
To eliminate the influence of genetic background on chilling tolerance, we used vegetatively propagated adventitious shoots of S. ionantha ‘Kilauea’ obtained by cutting leaves from the same plant. The plants were grown under conventional conditions at 20°C, 14 h light/10 h dark, using plant growth fluorescent lamps (BioLux, 20W, HotaluX, PPFD 82.5 µmol cm− 2 s− 1). Plants were potted in 220-ml plastic pots using Metro-Mix 360 soil (Sun Gro Horticulture, Massachusetts, United States). The temperature data in this study were consistently measured using a data logger (Thermo Recorder Ondotori Jr. TR-50U2, T&D, Japan) and carefully kept within 20 ± 1°C.
Identification of chilling injury-causing temperatures
The temperature at which Saintpaulia exhibits chilling injury slightly varies depending on the environment and plant cultivar (Yang et al. 2003). To identify that temperature in this study, plants grown at 20°C were transferred to an incubator with temperatures varying between 4°C and 20°C. Plants at their flowering stage with five or more leaves were used in all treatments. Each temperature treatment was conducted for 24 h in an incubator (MIR-553, SANYO, Japan) which was equipped with a circulator to stabilize the temperature. The following conditions were tested: 4, 6, 8, 10, 12, 16, and 20°C under a 12 h light/12 h dark cycle. After that, electrolyte leakage was measured as an indicator of chilling injury as described below.
Fifth and sixth leaves (5–6 cm diameter) from the youngest expanded leaves were collected from each plant, and three leaf disks were obtained per leaf using a biopsy punch (Kai Industries, Gifu, Japan. BP-60F. φ6mm). The leaf disks were washed with distilled water for about 10 seconds to remove electrolytes from the cut surface, soaked in 3 ml of distilled water, and shaken at 80 rpm at 20°C under dark conditions for 24 h. Next, the electrical conductivity (EC) of the water was measured using a compact conductivity meter, LAQUAtwin (HORIBA, Kyoto, Japan). Samples were then autoclaved to destroy all cells and let the electrolytes leak into the solution, and the EC was measured again. The electrolyte leakage was calculated as the ratio EC before autoclave to EC after autoclave and expressed as a percentage. Thereafter, electrolyte leakage rates were measured in the same manner.
Next, plants were placed at 4°C, which is a stable temperature for chilling injury to occur, in the dark to mimic transportation conditions. Plants were left at 4°C for 1, 6, 12, 18, and 24 h to identify the best timing to induce chilling injury and electrolyte leakage was measured.
Investigation of chilling acclimation conditions
To investigate the suitable temperature for chilling acclimation in S. ionantha ‘Kilauea’, plants were placed at 6, 8, 10, 12, 14, 16, 18, and 20°C for 24 h followed by a 24 h treatment at 4°C. After the treatments, electrolyte leakage was measured from leaves, and the best acclimation-inducing temperature was identified.
From the result of the experiment, 10°C was identified as the acclimation-inducing temperature. To investigate the optimal time for chilling acclimation at 10°C, plants were moved to a 10°C incubator for different durations. After 12, 24, 48, 72, and 144 h, plants were exposed to a chilling treatment at 4°C for 24 h, and the electrolyte leakage of leaves was measured. Plant images were captured before/after the chilling treatment at 4°C. The survival rate was measured 20 days after the chilling treatment with or without pretreatment at 10°C for 144 h.
Genome sequencing, assembly, gene prediction and functional annotation
Genomic DNA was extracted from young leaves of S. ionantha ‘Kilauea’ with the CTAB (Cetyl trimethyl ammonium bromide) method (Dempster et al. 1999). A long-read sequence library was constructed using a SMRTbell Express Template Prep Kit 1.0 (PacBio, Menlo Park, CA, USA) and sequenced on SMRT cells (1 M v3) in a PacBio Sequel system. A short-read library was prepared with the TruSeq Nano DNA Sample Prep Kit (Illumina, San Diego, CA, USA) and sequenced on a MiSeq (Illumina) instrument in a paired-end, 301-bp mode.
The Sequel long-reads were assembled by FALCON v1.2.5 (Chin et al. 2016), and the resultant primary contigs and associated contigs were applied to FALCON-Unzip v1.1.5 (Chin et al. 2016) for phasing. The primary contigs were polished by Arrow (https://github.com/PacificBiosciences/GenomicConsensus) using Sequel long-reads, and further polished by Pilon (Walker et al. 2014) using Illumina paired-end reads.
Ab initio gene prediction has been conducted on the polished primary contigs using BRAKER2 v2.1.5 (Hoff et al. 2019). The quality trimming of RNA-seq reads sampled from petal (1_blue; DRR462688, Kira_P-petal; DRR462690, Kira_Str-Pin; DRR462691, Kira_Str-Pout; DRR462692) and leaf (Kira_P-leaf; DRR462689) samples was performed by PRINSEQ v0.20.4 (Schmieder and Edwards 2011), and adaptor trimming was further performed by FASTX-toolkit (http://hannonlab.cshl.edu/fastx_toolkit). The trimmed reads were mapped onto the polished primary contigs by HISAT2 v2.2.0 (Kim et al. 2019), and the bam files were used for BRAKER2 as a training set. Similarity searches against the predicted genes were conducted against UniProtKB (The UniProt Consortium, 2023), NCBI’s NR (https://www.ncbi.nlm.nih.gov/refseq/about/nonredundantproteins/) using DIAMOND (Buchfink et al. 2021) in a more sensitive mode, and the protein sequences in Araport 11 (Cheng et al. 2017) using BLASTP (Camacho et al. 2009). Functional annotation based on orthology relationships against EggNOG was conducted by eggNOG-mappar (http://eggnog-mapper.embl.de/, Huerta-Cepas et al. 2017).
Determination and classification of candidate ERFs in Saintpaulia genome
Based on the annotation information by the eggNOG-mapper and the BLASTP results, 161 Saintpaulia ERFs were selected. A molecular phylogenetic tree was constructed using the amino acid sequences of 161 Saintpaulia ERFs, 121 A. thaliana ERFs amino acid sequences obtained from GenBank, and 137 S. lycopersicum ERFs amino acid sequences obtained from PlantTFDB (Guo et al. 2007, http://planttfdb.gao-lab.org/). The amino acid sequences were filtered using PREQUAL v1.0.2 (Whelan et al. 2018). Subsequently, multiple sequence alignment was performed using MUSCLE v5.1 (Edgar 2022) with super5 default settings. The aligned sequences were further trimmed using TrimAl v1.4.1 (Capella-Gutierrez et al. 2009). Using trimmed alignment sequences, a maximum likelihood (ML) method-based molecular phylogenetic tree was constructed by IQ-TREE2 v2.2.0.3 (Minh et al. 2020). When constructing, the Ultrafast bootstrap method (UF-boot) and SH-aLRT methods were used 1,000 times, and branches with confidence levels of ≥ 95% and ≥ 80% respectively were considered robust (Anisimova et al. 2011; Hoang et al. 2018). The clade to which Saintpaulia ERFs belong was estimated according to Nakano et al. (2006).
RNA-seq analysis
To determine which genes are expressed during chilling acclimation in Saintpaulia, RNA-seq analysis was performed. Total RNA was extracted from the youngest fully expanded leaf at 0, 12, 24, 48, and 144 h after treatment at 10°C using Sepasol®-RNA I Super G (Nacalai tesque, Kyoto, Japan), according to the manufacturer’s protocol. RNA-seq was performed by Rhelixa (Kanagawa, Japan). RNA-seq libraries were prepared using NEBNext® Poly(A) mRNA Magnetic Isolation Module and NEBNext® UltraTM II Directional RNA Library Prep Kit (New England Biolabs). Paired-end sequencing was performed using the Illumina NovaSeq 6000 with a read length of 150 bp. Adapters and low-quality reads were removed using Trimmomatic v0.39 (Bolger et al. 2014). The SIO_r2.0 genome data was used as a reference sequence and mapping was performed using HISAT2 v2.2.1 (Kim et al. 2019). Transcripts per million (TPM) values of the predicted genes were calculated by Stringtie v2.2.1 (Pertea et al. 2015).
GO enrichment analysis was performed using GO terms of genes whose expression was increased more than 2-fold by 10°C treatment. The R package topGO v2.41.0 was used for GO enrichment analysis (with settings; algorithm = “elim”, statistic = “fisher”) (Alexa and Rahnenfuhrer 2023). The top 10 significant GO terms enriched at each treatment time after 10°C are shown in the heatmap.
The TPM expression data of 161 Saintpaulia ERFs were used to create a heatmap using the R package ComplexHeatmap v2.14.0 (Gu et al. 2016, Gu 2022) (with cluster method settings = ward.D2).
Quantitative RT-PCR
The expression levels of some ERFs were confirmed by quantitative RT-PCR (qRT-PCR). Total RNA was extracted from Saintpaulia leaves previously treated at 10°C for 0 (no treatment as control), 3, 6, 12, 24, 48, 72 and 144 h using Sepasol®-RNA I Super G (Nacalai tesque). First-strand cDNA was synthesized from the extracted RNA using ReverTra Ace® (TOYOBO, Osaka, Japan), and qRT-PCR was performed using the cDNA and SYBR THUNDER BIRD® qPCR Mix (TOYOBO) with the CFX ConnectTM Real-Time PCR Detection System (Bio-Rad, California, United States). Primers were designed to amplify Saintpaulia ERFs that may be associated with chilling tolerance and were used for qRT-PCR (Supplementary Table S1). Relative quantification was performed using the standard curve method. A standard curve with a correlation coefficient of r ≥ 0.99 was used. SiActin (Sio_r2.0_p0205.1.g00590.1) was used as a housekeeping gene.
No injury occurred in the 24 h treatment at 6°C, but neither did chilling acclimation (described in results section). Therefore, we extracted RNA from leaves treated at 6°C for 24 h and compared the expression of the ERFs between no treatment (0 h, as control), 6°C 24 h treatment, and 10°C 24 h treatment by qRT-PCR. The procedure is the same as described above.
Promoter analysis
We obtained 2 kb upstream from the start codon of Saintpaulia ERFs whose expression was increased by low temperature and performed promoter analysis on MEME Suite (Bailey et al. 2015). Motif enrichment analysis was performed using MEME (v.5.5.3) with “-mod zoops -minw 5 -maxw 20 -markov_order 0 -objfun classic” settings (Bailey et al. 1994). A motif with an E-value < 0.05 was considered a significant motif. Enriched motifs were used for Tomtom (Gupta et al. 2007). Enriched motifs were used as queries to search for similar motifs in the A. thaliana DAP-seq motif database (E-value < 0.05, O’ Malley et al. 2016).
Statistical analysis
All statistical analyses were performed using R software (ver. 4.2.1). Tukey’s range test was performed using the stats package, and Dunnett’s test was performed using the multcomp package (Hothorn et al. 2008).
Chilling injury in S. ionantha ‘Kilauea’
We investigated the chilling injury in S. ionantha ‘Kilauea’. We exposed the plants to various temperature conditions for 24 h and calculated electrolyte leakage as an indicator of chilling injury. The electrolyte leakage ratio was approximately 20% and lower at temperatures higher than 6°C (Fig. 1A). However, when exposed to 4°C for 24 h, the electrolyte leakage ratio was approximately 60% (Fig. 1A). We then investigated the relationship between the exposure time at 4°C and the occurrence of chilling injury using plants grown under the same conditions. The electrolyte leakage was similar to non-treated conditions for up to 12 h after exposure to 4°C (Fig. 1B), but it increased rapidly after 12 h (Fig. 1B). Plants exposed to 4°C for 24 h exhibited withered young leaves near their shoot tip and subsequently died (Fig. 2C, D). Based on these results, it was determined that a chilling treatment at 4°C for 24 h would be used as a chilling treatment in subsequent experiments.
Chilling acclimation conditions for S. ionantha ‘Kilauea’
To investigate the optimal low temperature for chilling acclimation, plants were exposed to different temperatures for 24 h before undergoing the chilling treatment at 4°C. When plants were exposed to temperatures below 8°C for 24 h before the chilling treatment, electrolyte leakage ratio was around 80%, indicating severe chilling injury (Fig. 2A). However, when plants were exposed to 10 − 16°C for 24 h, electrolyte leakage was lower than 6°C or 8°C treatment, indicating reduced chilling injury (Fig. 2A). Especially, results of comparison the normal growing conditions (20°C) showed that treatment at both 10 and 12°C was most effective. Although both temperatures are effective for chilling acclimation, in this study, the 10°C treatment was chosen as the treatment temperature for chilling acclimation.
To investigate the optimal duration of chilling acclimation, plants were exposed to 10°C for varying periods before the chilling treatment. Electrolyte leakage gradually decreased with increasing exposure duration up to 72 h, and slightly increased up to 144 h at 10°C (Fig. 2B). When chilling acclimated plants were subjected to chilling treatment, mid- to outer leaves (older than eighth leaves from the youngest expanded leaf) turned pale but did not dry up. The shoot tip stayed alive without withering even after 20 days of returning to recovery temperature, unlike non-acclimated plants (Fig. 2C). The acclimation treatment avoided leaf and plant death, but leaf bleaching on outer leaves was observed (Fig. 2C). After 144 h of 10°C pretreatment, all plants survived (Fig. 2D).
Genome assembly, gene prediction, and determination of candidate Saintpaulia ERF
The statistics of the assembled genome sequences are summarized in Table 1, and Table S2. In total, 71.6 Gb of continuous long reads were obtained from 12 SMRT cells. The 12 cells of the Sequel long-reads were assembled by FALCON, and the 2,967 primary contigs and 1,266 associate contigs were constructed. The primary contigs and associated contigs were phased by FALCON-Unzip, and the 2,200 primary contigs and 7,658 haplotigs were constructed. Next, the primary contigs and haplotigs were polished by Arrow and Pilon. The resultant 2,200 primary contigs and 7,658 haplotigs were designated SIO_r2.0. The BUSCO completeness of the primary contigs was 97.9% (Table 1 and Table S2).
Table 1 Statistics of the assembled genome sequences and predicted genes (SIO_r2.0)
The statistics of the predicted genes are summarized in Table 1 and Table S2. The gene prediction was performed against the 2,200 primary contigs by BRAKER2, and 152,177 genes were predicted. Next duplicated genes in splicing variants were removed, and then 145,403 genes were defined as “best genes”. The “best genes” were searched against UniProtKB and NCBI’s NR by DIAMOND in the more sensitive mode with an e-value of 1E-100, identity ≥ 75%, and length coverage ≥ 50%, and those were searched against the proteins encoded in the genome sequence of Dorcoceras hygrometricum XS01 (Accession: GCA_001598015.1) by BLASTP with E-value of 1E-100, identity ≥ 60%, length coverage ≥ 50% and those in A. thaliana provided from Araport 11 (Cheng et al. 2017) with an e-value of 1E-100, and identity ≥ 40%, and length coverage ≥ 50%. The “best genes” were also searched against the protein domain database Pfam31.1 by HMMER v3.3 (Eddy 2011) with an e-value of 1e-20. The expressed genes were searched by mapping of the RNA-seq data applied to the gene prediction by Salmon v1.2.1 (Roberts and Pachter 2013) with a TPM value ≥ 0.4. The 46,536 displaying similarities against the proteins of UniProtKB and NR related to transposable elements were classified as Transposable Element (TE) genes. The 37,682 “best genes” displaying similarities and expressions in all of the analyses described above were classified as high-confidence (HC) genes, and the remaining 61,185 genes were classified as low-confidence (LC) genes. The BUSCO completeness of the genes classified into HC was 93.8%. The 37,682 genes classified in HC were used for further analyses as intrinsic genes.
In the S. ionantha genome database SIO_r2.0, 161 sequences were annotated as ERF transcription factors. Based on previous reports, the ERFs of Saintpaulia were classified into 12 clades (Nakano et al. 2006). The results were as follows: 10 ERFs in clade I, 19 in clade II, 34 in clade III, 10 in clade IV, 9 in clade V, 11 in clade VI, 4 in clade VI-L, 5 in clade VII, 18 in clade VIII, 28 in clade IX, 11 in clade X, 2 in clade Xb-L (Fig. 3). ERFs of III clade were the most abundant, followed by IX clade.
Expression analysis during chilling acclimation and promoter analysis
To confirm the response of the plants during chilling acclimation, GO enrichment analysis was performed on the genes whose expression increased more than 2-fold after the 10°C treatment. In Figure S1, the top 10 most significantly enriched GO terms at each treatment time are shown. GO terms at each treatment time included GO terms related to stress response and phytohormone.
The 161 Saintpaulia ERFs found in SIO_r2.0 were clustered based on their expression data during chilling acclimation (Fig. 4A). As a result, some ERFs formed a cluster (Fig. 4A). We focused on this cluster and designated it as "cluster of interest (COI)" consisting of 61 ERFs (Fig. 4A). The TPM expression levels of ERFs belonging to the COI were relatively higher than those of ERFs not belonging to the COI after the 10°C treatment (Fig. 4B). Among them, 43, 8, and 4 ERFs (a total of 55 ERFs) showed more than 2-fold increase in expression after 12 h, 24 h, and 48 h of treatment at 10°C, respectively. From the 43 ERFs whose expression increased at 12 h, we compared the expressions of 7 ERFs belonging to different clades that are reported to be involved in cold acclimation by qRT-PCR. In addition, we examined the expression levels of these 7 ERFs at 3 h and 6 h after the 10°C treatment. Except for Sio_r2.0_p0024.1.g05000.1, the results of qRT-PCR and RNA-seq were almost identical (Fig. 5, Fig. S2). The expression levels of some ERFs were already increased after 6 h (Fig. 5). Expression of ERFs homologous to CBF/DREB1 was not increased by 10°C treatment (Fig. S3).
The Saintpaulia ERFs whose expression increased greatly at an early stage after the 10°C treatment are probably involved in the chilling acclimation. Therefore, we analyzed the promoters of 55 ERFs from the COI, whose expression increased more than 2-fold after 12, 24, and 48 hours. As a result, 9 motifs were enriched from ERFs upregulated at 12 h and 2 motifs from ERFs upregulated at 24 h (Tables S5 and S6). No significant motifs were obtained from the ERFs upregulated at 48 h. The 11 enriched motifs were used to investigate the homology with the motifs of A. thaliana by Tomtom (Gupta et al. 2007). Except for Motif 2/3, no similar motifs were found. When Motif 2 was used as a query, several motifs of GCC box-like motif, to which ERF transcription factors bind, were hit (with E-value < 0.05, Table S7). When Motif 3 was used as a query, several DRE/CRT core sequences like motif, to which DREB subfamily members bind, were hit (Table S7).
Finally, we compared the expression of ERFs between the treatments for 24 h at 10°C, an effective temperature for chilling acclimation, and at 6°C, an ineffective temperature. The results showed that the expressions of all six ERFs we investigated were not or only barely increased by the 6°C treatment (Fig. 6).
Saintpaulia, a subtropical plant, is capable of chilling acclimation
Saintpaulia is extremely sensitive to low temperatures, and once chilling injury occurs, the plant does not re-green (Kratsch and Wise 2000). In this study, the electrolyte leakage of S. ionantha 'Kilauea' grown at a stable temperature of 20°C increased rapidly after 24 h of exposure to 4°C, causing the plant to die (Fig. 1A, B and Fig. 2C). In contrast to the 4°C treatment, electrolyte leakage did not increase after 24 h of treatment above 6°C (Fig. 1A). A pretreatment at 10°C before the chilling treatment at 4°C suppressed the increase in electrolyte leakage, resulting in chilling acclimation (Fig. 2). Even in extremely low temperature-sensitive plants such as mung bean (Vigna radiata) and upland cotton (Gossypium hirsutum), chilling acclimation can develop chilling tolerance (Kratsch and Wise 2000; Chang et al. 2001; Kargiotidou et al. 2010). Our results indicate Saintpaulia can also undergo chilling acclimation and acquire chilling tolerance.
It has been reported that cold water applied to the leaves of Saintpaulia results in the death of palisade cells. This disorder is reported to be caused by a sudden drop in leaf temperature. For example, a sudden drop in leaf temperature from 35°C to 25°C also injured palisade cells (Elliot 1946; Maekawa et al. 1987; Yang et al. 2003), which is said to be different from cell death caused by chilling stress. In our study, some leaves of plants that survived the 4°C treatment by chilling acclimation (at 10°C for 72–144 h) showed partial bleaching, similar to the injury caused by a sudden temperature drop (Fig. 2C). This bleaching was not observed in plants at 4°C without acclimation treatment., suggesting that bleaching phenomenon was not caused by a sudden temperature drop in 4°C treatment. However, the similarity in appearance with the bleaching caused by a sudden temperature drop suggests that the injury caused by a sudden temperature drop and the low-temperature injury have something in common.
GO enrichment analysis showed that the enriched terms included not only 'response to cold' but also 'response to abscisic acid' and 'trehalose metabolism in response to stress' (Fig. S1). Abscisic acid is a well-known stress hormone and is synthesized in response to various stresses, including low-temperature stress (Swamy et al. 1999). Trehalose is also known to respond to various stresses and act as an osmolyte and a protectant of cell membranes and proteins (Iordachescu et al. 2008). The fact that these GO terms were observed in Saintpaulia after 10°C treatment suggests that abscisic acid and trehalose may contribute to the improvement of low-temperature tolerance, like other plants.
The expression of multiple ERFs is induced during chilling acclimation in Saintpaulia
According to the whole genome sequence (Table 1, Table S2), 161 ERFs were determined in Saintpaulia, with the largest number of ERFs belonging to clade III, followed by clade IX (34 in clade III and 28 in clade IX) (Fig. 3). The increase in AP2/ERF family transcription factors through whole genome and tandem duplication may have contributed to the adaptation of plant ancestors to low-temperature stress caused by climate change (Guo et al. 2022). In addition, it has been suggested that the number of ERFs belonging to clade III and IX was the most expanded by whole genome duplication or tandem duplication (Guo et al. 2022). Also in Saintpaulia, it is possible that whole genome duplication or tandem duplication is responsible for the increase in the number of ERFs in the various clades, especially in clades III and IX. Clade III includes the CBF/DREB1 transcription factors, which have been extensively studied for their role in cold tolerance in many plants, including A. thaliana (Jaglo-Ottosen et al. 1998; Liu et al. 1998; Gilmour et al. 2000; Dubouzet et al. 2003; Gilmour et al. 2004; Qin et al. 2004; Zhang et al. 2004). In Saintpaulia, ERFs similar to A. thaliana and S. lycopersicum CBF/DREB1 were found (Fig. 3, Fig. S3). However, these ERFs did not show an increased expression under low temperature; instead, their expression levels decreased (Fig. 4A, Fig. 2). These results suggest that Saintpaulia is unable to induce the expression of CBF/DREB1 in response to low temperatures. Taiwan banana (Musa acuminata), native to the Malay Peninsula, and Sea Island cotton (Gossypium barbadense), native to the Caribbean Sea, are tropical and subtropical plants growing in areas where the temperature remains stable, around 26°C. In these plants, the expression of CBF/DREB1 was upregulated by low temperatures (MaDREB1F, GbCBF1), and the genetic transformation of these CBF/DREB1 improved low-temperature tolerance (Guo et al. 2011; Xu et al. 2023). Thus, CBF/DREB1 is upregulated in response to low temperatures even in tropical plants that do not usually encounter low temperatures, and enhances chilling tolerance by upregulation. As the cbf2 mutant in A. thaliana has enhanced cold tolerance (Novillo et al. 2004), the downregulated expression of CBF/DREB1 like may contribute to the enhanced low temperature tolerance in Saintpaulia. Unlike ERFs homologous to CBF/DREB1 in Saintpaulia, other ERFs from various clades exhibited an increased expression under low temperatures (Fig. 4, 6). AtDDF1 (clade III) and AtCRF4 (clade VI) in A. thaliana are known to improve cold tolerance (Kang et al. 2011; Zwack et al. 2016). Additionally, VaERF057 (clade VII) in V. amurensis, PtERF9 (clade VIII) and PtERF108/109 (clade X) in P. trifoliata and MdABI4 (clade IV) in M. domestica are also involved in enhancing cold tolerance (An et al. 2022; Sun et al. 2016; Wang et al. 2019; Khan et al. 2021; Zhang et al. 2022). In many plants, the increased expression of ERFs belonging to various clades has been reported to be involved in the improvement of low-temperature tolerance. Therefore, the increased expression of various clades may also be involved in chilling acclimation in Saintpaulia. When the pretreatment duration at 10°C was 24 − 48 h, the increase in electrolyte leakage caused by the subsequent 4°C treatment was suppressed. However, the plants were unable to recover after the 4°C treatment and eventually withered (Fig. 2C). When pre-treated at 10°C for 72 h or more, the plants survived after the subsequent 4°C treatment without withering (Fig. 2C, D). The increased expression of ERFs after 12 h of treatment at 10°C suggests a time lag between that increase expression of ERFs and the acquisition of chilling tolerance by acclimation. Although the expression of several ERFs, which are expected to be associated with chilling acclimation, was increased at 10°C (Fig. 5), when the plants were exposed to 6°C, expression levels of these ERFs are not upregulated (Fig. 6). This may suggest that ERF expression is important for Saintpaulia chilling acclimation.
The promoter regions of some Saintpaulia ERFs whose expression increased during cold acclimation had cis-elements to which AP2/ERF transcription factors could bind. Thus, the ERFs whose expression is increased by low temperature may regulate the expression of other ERFs by binding to the promoter region of their own or other ERFs. However, since GCC boxes and DRE/CRT core motifs were also found in the promoter regions of ERFs whose expression was not increased by low temperature, the increased expression of multiple ERFs may require other factors. Also, it is reported that the ability of ERFs to bind to cis-elements varies even in the same clade of ERFs (Shoji et al. 2013). Such ability variation is thought to be caused by the difference in the amino acid sequence of the AP2 domain, which is the DNA-binding region. Among the Saintpaulia 55 ERFs induced by low temperature, the sequence of each AP2 domain is different (Fig. S4). This suggests that each ERF regulates different downstream genes, and the expression of various ERFs may induce more downstream genes to improve chilling tolerance.
This study revealed that Saintpaulia could be chilling acclimatized by exposure to 10°C for more than 72 h, and can survive at 4°C for 24 h, which would normally kill the plants. In addition, we constructed the SIO_r2.0 genome data and found that 161 ERF transcription factors, including CBF/DREB1, were present in Saintpaulia. The expression of multiple ERFs from various clades was upregulated in response to acclimation at 10°C.
Chilling injury is a major problem in the production and transportation of Saintpaulia in temperate regions. Near perspectives include confirming whether the chilling acclimation of Saintpaulia can have a practical effect during cultivation and transportation and optimizing the acclimation conditions by testing various conditions and cultivars. In addition, if the genes associated with chilling acclimation can be identified, the genetic markers or the expression of these genes could be used as indicators to evaluate chilling tolerance in Saintpaulia and be helpful with transportation challenges.
Funding
This work was supported by JSPS KAKENHI Grant Number 16H06279 (PAGS), 17H03769, 22H05172, 22H05181, the Kazusa DNA Research Institute Foundation and Kindai University foundation (A Grant for Scientific Research from the Faculty of Agriculture, Kindai University).
Acknowledgements
We thank S. Nakayama, C. Minami, S. Sasamoto, H. Tsuruoka, and A. Watanabe (Kazusa DNA Research Institute) and for technical assistance. We also thank T. Tsuzaki and R. Kunimune (Graduate School of Agriculture, Kindai University) for their help in cultivation management.
Conflict of interest The authors declare that they have no conflict of interest.
Author contributions
Daichi Kurata: Data curation, Investigation, Resources, Visualization, Writhing-Original draft. Kento Fukutomi: Investigation, Resources. Kanae Kubo: Investigation, Resources. Kenta Shirasawa: Funding acquisition, Data curation, Formal analysis, Writing-Original draft. Hideki Hirakawa: Funding acquisition, Data curation, Formal analysis, Writing-Original draft. Munetaka Hosokawa: Funding acquisition, Supervision, Project administration, Conceptualization, Writing-review & editing.
Data availability
The genome sequence and raw data are available under the BioProject PRJDB15710. The accession numbers of the raw data are summarized into Table S3. The accession numbers of the amino acid of A. thaliana ERFs used in phylogenetic tree construct are summarized into Table S4. The cds and pep files of the predicted genes (HC, LC, and TE) and gff3 files were provided from FigShare (DOI: 10.6084/m9.figshare.24548854). (for reviewer; https://figshare.com/account/items/24548854/edit.
- Alexa A, Rahnenfuhrer J (2023) topGO: enrichment analysis for gene ontology
- An JP, Xu RR, Liu X, Su L, Yang K, Wang XF, Wang GL, You CX (2022) Abscisic acid insensitive 4 interacts with ICE1 and JAZ proteins to regulate ABA signaling-mediated cold tolerance in apple. J Exp Bot 73:980–997. https://doi.org/10.1093/jxb/erab433
- Anisimova M, Gil M, Dufayard JF, Dessimoz C, Gascuel O (2011) Survey of branch support methods demonstrates accuracy, power, and robustness of fast likelihood-based approximation schemes. Syst Biol 60:685–699. https://doi.org/10.1093/sysbio/syr041
- Baatvik ST (1993) The genus Saintpaulia (Gesneriaceae) 100 years: history, taxonomy, ecology, distribution and conservation. Fragm Flor Geobot Suppl 2:97–112
- Bailey TL, Elkan C (1994) Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol 2:28–36
- Bailey TL, Johnson J, Grant CE, Noble WS (2015) The MEME suite. Nucleic Acids Res 43(W1):W39–49. https://doi.org/10.1093/nar/gkv416
- Bolger AM, Lohse M, Usadel B (2014) Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30:2114–2120. https://doi.org/10.1093/bioinformatics/btu170
- Bolt S, Zuther E, Zintl S, Hincha DK, Schmülling T (2017) ERF105 is a transcription factor gene of Arabidopsis thaliana required for freezing tolerance and cold acclimation. Plant Cell Environ 40:108–120. https://doi.org/10.1111/pce.12838
- Buchfink B, Reuter K, Drost HG (2021) Sensitive protein alignments at tree-of-life scale using DIAMOND. Nat Methods 18:366–368. https://doi.org/10.1038/s41592-021-01101-x
- Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL (2009) BLAST+: architecture and applications. BMC Bioinform 10:421. https://doi.org/10.1186/1471-2105-10-421
- Capella-Gutierrez S, Silla-Martinez JM, Gabaldon T (2009) trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25:1972–1973. https://doi.org/10.1093/bioinformatics/btp348
- Chang MY, Chen SL, Lee CF, Chen YM (2001) Cold-acclimation and root temperature protection from chilling injury in chilling-sensitive mungbean (Vigna radiata L) seedlings. Bot Bull Acad Sin 42:53–60
- Chen J, McConnell DB, Henny RJ, Norman DJ (2004) The foliage plant industry. Hort Rev (Amer Soc Hort Sci) 31:47–112. https://doi.org/10.1002/9780470650882.ch2
- Cheng C, Krishnakumar V, Chan AP, Thibaud-Nissen F, Schobel S, Town CD (2017) Araport11: a complete reannotation of the Arabidopsis thaliana reference genome. Plant J 89:789–804. https://doi.org/10.1111/tpj.13415
- Chin CS, Peluso P, Sedlazeck FJ, Nattestad M, Concepcion GT, Clum A, Dunn C, O’Malley R, Figueroa-Balderas R, Morales-Cruz A, Cramer GR, Delledonne M, Luo C, Ecker JR, Cantu D, Rank DR, Schatz MC (2016) Phased diploid genome assembly with single-molecule real-time sequencing. Nat Methods 13:1050–1054. https://doi.org/10.1038/nmeth.4035
- Dempster EL, Pryor KV, Francis D, Young JE, Rogers HJ (1999) Rapid DNA extraction from ferns for PCR–based analyses. Biotechniques 27:66–68. https://doi.org/10.2144/99271bm13
- Dubouzet JG, Sakuma Y, Ito Y, Kasuga M, Dubouzet EG, Miura S, Seki M, Shinozaki K, Yamaguchi-Shinozaki K (2003) OsDREB genes in rice, Oryza sativa L, encode transcription activators that function in drought-, high-salt- and cold-responsive gene expression: DREB transcription activators in rice. Plant J 33:751–763. https://doi.org/10.1046/j.1365-313X.2003.01661.x
- Eddy SR (2011) Accelerated profile HMM searches. PLoS Comput Biol 7:e1002195. https://doi.org/10.1371/journal.pcbi.1002195
- Edgar RC (2022) Muscle5: high-accuracy alignment ensembles enable unbiased assessments of sequence homology and phylogeny. Nat Commun 13:6968. https://doi.org/10.1038/s41467-022-34630-w
- Elliot FH (1946) Saintpaulia leaf spot and temperature differential. Proc Amer Soc Hort Sci 47:511–514
- Gilmour SJ, Fowler SG, Thomashow MF (2004) Arabidopsis transcriptional activators CBF1, CBF2, and CBF3 have matching functional activities. Plant Mol Biol 54:767–781. https://doi.org/10.1023/B:PLAN.0000040902.06881.d4
- Gilmour SJ, Hajela RK, Thomashow MF (1988) Cold acclimation in Arabidopsis thaliana. Plant Physiol 87:745–750. https://doi.org/10.1104/pp.87.3.745
- Gilmour SJ, Sebolt AM, Salazar MP, Everard JD, Thomashow MF (2000) Overexpression of the Arabidopsis CBF3 transcriptional activator mimics multiple biochemical changes associated with cold acclimation. Plant Physiol 124:1854–1865. https://doi.org/10.1104/pp.124.4.1854
- Gu Z, Eils R, Schlesner M (2016) Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics 32:2847–2849. https://doi.org/10.1093/bioinformatics/btw313
- Gu Z (2022) Complex heatmap visualization. IMeta 1:e43. https://doi.org/10.1002/imt2.43
- Guo AY, Chen X, Gao G, Zhang H, Zhu QH, Liu XC, Zhong YF, Gu X, He K, Luo J (2007) PlantTFDB: a comprehensive plant transcription factor database. Nucleic Acids Res 36(suppl1):D966–D969. https://doi.org/10.1093/nar/gkm841
- Guo HM, Li ZC, Zhang H, Xin YZ, Cheng HM (2011) Cloning of cotton CBF gene for cold tolerance and its expression in transgenic tobacco. Acta Agron Sin 37:286–293. https://doi.org/10.1016/S1875-2780(11)60009-6
- Guo L, Wang S, Nie Y, Shen Y, Ye X, Wu W (2022) Convergent evolution of AP2/ERF III and IX subfamilies through recurrent polyploidization and tandem duplication during eudicot adaptation to paleoenvironmental changes. Plant Commun 3:100420. https://doi.org/10.1016/j.xplc.2022.100420
- Gupta S, Stamatoyannopoulos JA, Bailey TL, Noble WS (2007) Quantifying similarity between motifs. Genome Biol 8:R24. https://doi.org/10.1186/gb-2007-8-2-r24
- Hoang DT, Chernomor O, von Haeseler A, Minh BQ, Vinh LS (2018) UFBoot2: improving the ultrafast bootstrap approximation. Mol Biol Evol 35:518–522. https://doi.org/10.1093/molbev/msx281
- Hoff KJ, Lomsadze A, Borodovsky M, Stanke M (2019) Whole-genome annotation with BRAKER in: Kollmar M (eds) Gene Prediction Methods Mol Biol 1962:65–95. https://doi.org/10.1007/978-1-4939-9173-0_5
- Hothorn T, Bretz F, Westfall P (2008) Simultaneous inference in general parametric models. Biom J 50:346–363. https://doi.org/10.1002/bimj.200810425
- Huerta-Cepas J, Forslund K, Coelho LP, Szklarczyk D, Jensen LJ, von Mering C, Bork P (2017) Fast genome-wide functional annotation through orthology assignment by eggNOG-mapper. Mol Biol Evol 34:2115–2122. https://doi.org/10.1093/molbev/msx148
- Illgen S, Zintl S, Zuther E, Hincha DK, Schmülling T (2020) Characterisation of the ERF102 to ERF105 genes of Arabidopsis thaliana and their role in the response to cold stress. Plant Mol Biol 103:303–320. https://doi.org/10.1007/s11103-020-00993-1
- Iordachescu M, Imai R (2008) Trehalose biosynthesis in response to abiotic stresses. J Integr Plant Biol 50:1223–1229
- Jaglo-Ottosen KR, Gilmour SJ, Zarka DG, Schabenberger O, Thomashow MF (1998) Arabidopsis CBF1 overexpression Induces COR genes and enhances freezing tolerance. Science 280:104–106. https://doi.org/10.1126/science.280.5360.104
- Johansson DR (1978) Saintpaulias in their natural environment with notes on their present status in Tanzania and Kenya. Biol Conserv 14:45–62. https://doi.org/10.1016/0006-3207(78)90005-8
- Kang HG, Kim J, Kim B, Jeong H, Choi SH, Kim EK, Lee HY, Lim PO (2011) Overexpression of FTL1/DDF1, an AP2 transcription factor, enhances tolerance to cold, drought, and heat stresses in Arabidopsis thaliana. Plant Sci 180:634–641. https://doi.org/10.1016/j.plantsci.2011.01.002
- Kargiotidou A, Kappas I, Tsaftaris A, Galanopoulou D, Farmaki T (2010) Cold acclimation and low temperature resistance in cotton: Gossypium hirsutum phospholipase Dα isoforms are differentially regulated by temperature and light. J Exp Bot 61:2991–3002. https://doi.org/10.1093/jxb/erq124
- Khan M, Hu J, Dahro B, Ming R, Zhang Y, Wang Y, Alhag A, Li C, Liu J (2021) ERF108 from Poncirus trifoliata (L) Raf functions in cold tolerance by modulating raffinose synthesis through transcriptional regulation of PtrRafS. Plant J 108:705–724. https://doi.org/10.1111/tpj.15465
- Kim D, Paggi JM, Park C, Bennett C, Salzberg SL (2019) Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol 37:907–915. https://doi.org/10.1038/s41587-019-0201-4
- Kratsch HA, Wise RR (2000) The ultrastructure of chilling stress. Plant Cell Environ 23:337–350. https://doi.org/10.1046/j.1365-3040.2000.00560.x
- Larcher W, Bodner M (1980) Dosisletalitat-nomogramm zur charakterisierung der erkaltungsempfindlichkeit tropischer pflanzen (in German). Angew Bot 54:273–278
- Levitt J (1980) Responses of plants to environmental stresses, 2nd edn. Academic, New York and London
- Liu Q, Kasuga M, Sakuma Y, Abe H, Miura S, Yamaguchi-Shinozaki K, Shinozaki K (1998) Two transcription factors, DREB1 and DREB2, with an EREBP/AP2 DNA binding domain separate two cellular signal transduction pathways in drought- and low- temperature-responsive gene expression, respectively, in Arabidopsis. Plant Cell 10:1391–1406. https://doi.org/10.1105/tpc.10.8.1391
- Maekawa S, Torisu Y, Inagaki N, Terabun M (1987) Leaf injury caused by drop in leaf temperature of Saintpaulia ionantha (in Japanese). J Jpn Soc Hortic Sci 55:484–489
- Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, Lanfear R (2020) IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol 37:1530–1534. https://doi.org/10.1093/molbev/msaa015
- Nakano T, Suzuki K, Fujimura T, Shinshi H (2006) Genome-wide analysis of the ERF gene family in Arabidopsis and Rice. Plant Physiol 140:411–432. https://doi.org/10.1104/pp.105.073783
- Novillo F, Alonso JM, Ecker JR, Salinas J (2004) CBF2/DREB1C is a negative regulator of CBF1/DREB1B and CBF3/DREB1A expression and plays a central role in stress tolerance in Arabidopsis. Proc Natl Acad Sci USA 101:3985–3990. https://doi.org/10.1073/pnas.0303029101
- O'Malley RC, Huang SC, Song L, Lewsey MG, Bartlett A, Nery JR, Galli M, Gallavotti A, Ecker JR (2016) Cistrome and epicistrome features shape the regulatory DNA landscape. Cell 165:1280–1292. https://doi.org/10.1016/j.cell.2016.04.038
- Pertea M, Pertea GM, Antonescu CM, Chang TC, Mendell JT, Salzberg SL (2015) StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol 33:290–295. https://doi.org/10.1038/nbt.3122
- Qin F, Sakuma Y, Li J, Liu Q, Li YQ, Shinozaki K, Yamaguchi-Shinozaki K (2004) Cloning and functional analysis of a novel DREB1/CBF transcription factor involved in cold-responsive gene expression in Zea mays L. Plant Cell Physiol 45:1042–1052. https://doi.org/10.1093/pcp/pch118
- Roberts A, Pachter L (2013) Streaming fragment assignment for real-time analysis of sequencing experiments. Nat Methods 10:71–73. https://doi.org/10.1038/nmeth.2251
- Schmieder R, Edwards R (2011) Quality control and preprocessing of metagenomic datasets. Bioinformatics 27:863–864. https://doi.org/10.1093/bioinformatics/btr026
- Shoji T, Mishima M, Hashimoto T (2013) Divergent DNA-binding specificities of group of ETHYLENE RESPONSE FACTOR transcription factors involved in plant defense. Plant Physiol 162:977–990. https://doi.org/10.1104/pp.113.217455
- Sun X, Zhang L, Wong DCJ, Wang Y, Zhu Z, Xu G, Wang Q, Li S, Liang Z, Xin H (2019) The ethylene response factor VaERF092 from Amur grape regulates the transcription factor VaWRKY33, improving cold tolerance. Plant J 99:988–1002. https://doi.org/10.1111/tpj.14378
- Sun X, Zhao T, Gan S, Ren X, Fang L, Karungo SK, Wang Y, Chen L, Li S, Xin H (2016) Ethylene positively regulates cold tolerance in grapevine by modulating the expression of ETHYLENE RESPONSE FACTOR 057. Sci Rep 6:24066. https://doi.org/10.1038/srep24066
- Swamy PM, Smith BN (1999) Role of abscisic acid in plant stress tolerance. Curr Sci 76:1220–1227. http://www.jstor.org/stable/24101946
- The UniProt Consortium (2023) UniProt: the universal protein knowledgebase in 2023. Nucleic Acids Res 51(D1):D523–D531. https://doi.org/10.1093/nar/gkac1052
- Thomashow MF (1990) Molecular genetics of cold acclimation in higher plants. Adv Genet 28:99–131. https://doi.org/10.1016/S0065-2660(08)60525-8
- Thomashow MF (1999) Plant cold acclimation: freezing tolerance genes and regulatory mechanisms. Annu Rev Plant Physiol Plant Mol Biol 50:571–599. https://doi.org/10.1146/annurev.arplant.50.1.571
- Walker BJ, Abeel T, Shea T, Priest M, Abouelliel A, Sakthikumar S, Cuomo CA, Zeng Q, Wortman J, Young SK, Earl AM (2014) Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS ONE 9:e112963. https://doi.org/10.1371/journal.pone.0112963
- Wang M, Dai W, Du J, Ming R, Dahro B, Liu J (2019) ERF 109 of trifoliate orange (Poncirus trifoliata (L) Raf) contributes to cold tolerance by directly regulating expression of Prx1 involved in antioxidative process. Plant Biotechnol J 1:1316–1332. https://doi.org/10.1111/pbi.13056
- Wang Y, Jiang H, Mao Z, Liu W, Jiang S, Xu H, Su M, Zhang J, Wang N, Zhang Z, Chen X (2021) Ethylene increases the cold tolerance of apple via the MdERF1B–MdCIbHLH1 regulatory module. Plant J 106:379–393. https://doi.org/10.1111/tpj.15170
- Whelan S, Irisarri I, Burki F (2018) PREQUAL: detecting non-homologous characters in sets of unaligned homologous sequences. Bioinformatics 34:3929–3930. https://doi.org/10.1093/bioinformatics/bty448
- Xin Z, Browse J (2000) Cold comfort farm: the acclimation of plants to freezing temperatures. Plant Cell Environ 23:893–902. https://doi.org/10.1046/j.1365-3040.2000.00611.x
- Xu Y, Hu W, Song S, Ye X, Ding Z, Liu J, Wang Z, Li J, Hou X, Xu B, Jin Z (2023) MaDREB1F confers cold and drought stress resistance through common regulation of hormone synthesis and protectant metabolite contents in banana. Hort Res 10:uhac275. https://doi.org/10.1093/hr/uhac275
- Yamaguchi-Shinozaki K, Shinozaki K (2006) Transcriptional regulatory networks in cellular responses and tolerance to dehydration and cold stresses. Annu Rev Plant Biol 57:781–803. https://doi.org/10.1146/annurev.arplant.57.032905.105444
- Yang SJ, Hayashi T, Hosokawa M, Yazawa S (2003) Leaf temperature drop measured by thermography and occurrence of leaf browning injury in Saintpaulia (in Japanese). Environ Control Biol 41:265–270. https://doi.org/10.2525/ecb1963.41.265
- Zhang X, Fowler SG, Cheng H, Lou Y, Rhee SY, Stockinger EJ, Thomashow MF (2004) Freezing-sensitive tomato has a functional CBF cold response pathway, but a CBF regulon that differs from that of freezing-tolerant Arabidopsis. Plant J 39:905–919. https://doi.org/10.1111/j.1365-313X.2004.02176.x
- Zhang Y, Ming R, Khan M, Wang Y, Dahro B, Xiao W, Li C, Liu J (2022) ERF9 of Poncirus trifoliata (L) Raf undergoes feedback regulation by ethylene and modulates cold tolerance via regulating a glutathione S-transferase U17 gene. Plant Biotechnol J 20:183–200. https://doi.org/10.1111/pbi.13705
- Zwack PJ, Compton MA, Adams CI, Rashotte AM (2016) Cytokinin response factor 4 (CRF4) is induced by cold and involved in freezing tolerance. Plant Cell Rep 35:573–584. https://doi.org/10.1007/s00299-015-1904-8