Skip to main content

Modulation of early gene expression responses to water deprivation stress by the E3 ubiquitin ligase ATL80: implications for retrograde signaling interplay

Abstract

Background

Primary response genes play a pivotal role in translating short-lived stress signals into sustained adaptive responses. In this study, we investigated the involvement of ATL80, an E3 ubiquitin ligase, in the dynamics of gene expression following water deprivation stress. We observed that ATL80 is rapidly activated within minutes of water deprivation stress perception, reaching peak expression around 60 min before gradually declining. ATL80, despite its post-translational regulation role, emerged as a key player in modulating early gene expression responses to water deprivation stress.

Results

The impact of ATL80 on gene expression was assessed using a time-course microarray analysis (0, 15, 30, 60, and 120 min), revealing a burst of differentially expressed genes, many of which were associated with various stress responses. In addition, the diversity of early modulation of gene expression in response to water deprivation stress was significantly abolished in the atl80 mutant compared to wild-type plants. A subset of 73 genes that exhibited a similar expression pattern to ATL80 was identified. Among them, several are linked to stress responses, including ERF/AP2 and WRKY transcription factors, calcium signaling genes, MAP kinases, and signaling peptides. Promoter analysis predicts enrichment of binding sites for CAMTA1 and CAMTA5, which are known regulators of rapid stress responses. Furthermore, we have identified a group of differentially expressed ERF/AP2 transcription factors, proteins associated with folding and refolding, as well as pinpointed core module genes which are known to play roles in retrograde signaling pathways that cross-referenced with the early ATL80 transcriptome.

Conclusions

Based on these findings, we propose that ATL80 may target one or more components within the retrograde signaling pathways for degradation. In essence, ATL80 serves as a bridge connecting these signaling pathways and effectively functions as an alarm signal.

Peer Review reports

Background

Cells respond to signals activating primary and secondary transcriptional programs. During primary responses, expression of an array of genes is induced within minutes, independently of de novo protein production; these genes are known as primary response genes (PRG)s. Secondary response genes are subsequently expressed that require de novo protein synthesis, and these genes are more numerous than primary response genes, which serve as master switches from brief and transient signals to extended responses [1]. Mammalian proto-oncogenes c-fos and c-jun are prototypes of primary response genes which are implicated in cell proliferation and differentiation in response to diverse stimuli (i.e., developmental and stress signals, growth factors, cytokines, xenobiotics) [2]. c-fos is among various primary response genes that are transiently expressed in neurons within minutes after synaptic activation [3]. Conversely, c-fos expression is constitutive and not regulated in many cancer cells [4].

In plants, many developmental and cellular processes and responses to physical injury and abiotic and biotic stresses are triggered by primary response genes. Auxin action guides many growth and development processes through the transcriptional regulation of a short signaling pathway. Auxin induces the expression of primary response genes within minutes, including members of the Auxin/Indole-3-Acetic Acid (Aux/IAA) and the auxin response factor (ARF) families. Members of the Aux/IAA and ARF families interact between them, generating a group of diverse dimers that function as transcriptional repressors routing many auxin responses [5]. The primary response cytokinin response genes comprise the type-B Response Regulator (RRBs) transcriptional factors. Transcription factor (TF) complexes originated by homo- and heterodimerization of RRBs and by interaction with other TF (i.e., DELLA proteins, EIN3, HD-ZIII, BZIP63, SPL9) mediate the response to this phytohormone. Both plants and animals have preserved an innate immune system during evolution to respond to tissue injury and pathogen attack. Cell surface receptors interact with foreign molecules, including endogenous damage-associated molecular patterns (DAMPs) and microbe- or pathogen-associated molecular patterns induce (MAMPs/PAMPs) to trigger strong defense responses. Primary response genes are essential components of these responses, with the perception of DAMPs and MAMPs/PAMPs leading to the activation of many primary immune response TFs and signaling components. These pathways of ten converge on regulatory mitogen-activated protein kinase (MAPK) cascades [6].

Plant-specific A-Type RING Ligases (ATLs), formerly known as Arabidopsis Tóxicos en Levadura, make up a large family of RING-type ubiquitin-ligases. They emerged around 470 million years ago when plants colonized land. Arabidopsis thaliana, for example, has 100 ATLs [7]. These ATLs are key players in the ubiquitination pathway, interacting with ubiquitin and ubiquitin-conjugating enzymes (E2s) to mark target proteins for degradation. A typical ATL comprises a RING-H2 finger for interaction with E2 conjugases, transmembrane helices at the N-terminus for directing ATLs to specific cellular compartments, and a conserved 12–16 amino acid motif known as GLD, located between the hydrophobic domain and the RING-H2 [7].

ATLs are frequently involved in plant immunity and responses to environmental stresses. For instance, ATL2 is rapidly and transiently activated in response to pathogen-associated molecular patterns (PAMPs), ATL9 is involved in the expression of defense-associated genes, cell death, and callose deposition, and ATL12 contributes to resistance against fungal pathogens [8, 9]. A pair of closely related ATL paralogs, ATL6 and ATL31, and ATL44/RHA3A and ATL45/RHA3B mediate the homeostasis of BOTRYTIS-INDUCED KINASE 1 (BIK1), a central component regulating immune responses. ATL6 and ATL31 facilitate BIK1 stability by targeting CALCIUM-DEPENDENT PROTEIN KINASE28 (CPK28), a regulator of BIK1 turnover, and ATL44/RHA3A and ATL45/RHA3B mediate the monoubiquitination of BIK1 which is essential to boost immune responses [10, 11]. ATLs have been shown to grant tolerance to environmental stresses. For instance, ectopic expression of ATL61 enhances drought tolerance in Arabidopsis, and the sweet potato IbATL38 improves salt stress tolerance [12, 13]. In rice, decreased expression of OsATL38 enhances tolerance to cold stress through mono-ubiquitination of the 14-3-3 protein OsGF14d [14].

ATL78 orthologs are implicated in tolerance to drought stress in A. thaliana and tomato. In A. thaliana, an atl78 mutant shows reduced tolerance to drought stress; consistently, in tomato, the ectopic expression of ShATL78 enhances tolerance to drought [15, 16]. Interestingly, ATL78 was linked to an adaptation to drought tolerance during the evolution of Brassicaceae. As a result of a segmental duplication event, ATL81 with diminished expression generated ATL78, a paralog with an enhanced expression that conferred a tolerance phenotype [17]. Among ATLs whose expression is induced in drought stress, we identified ATL80, which shows features of a primary response gene that its expression following water deprivation stress (WDS) [18]. To understand the impact of ATL80 on differential gene expression, a time-course microarray experiment was conducted to compare responses to WDS between wild-type and atl80 mutant plants. The results revealed that ATL80 significantly influenced gene expression, leading to a severe decrease in differentially expressed genes (DEGs) in the atl80 mutant compared to the wild type. Additionally, among the DEGs, components common to the retrograde signaling pathways were identified, suggesting the existence of converging signaling pathways.

Methods

Arabidopsis thaliana growth conditions and water deprivation treatments

We utilized WT Arabidopsis thaliana Col-0 and the mutant atl80 (bearing a T-DNA insertion SALK_046204) for our research. The SALK_046204 mutant line was acquired from ABRC, genotyped using the SIGnAL iSect Toolbox (http://signal.salk.edu/isects.html), and underwent backcrossing prior to analysis. The T-DNA insertion in SALK_046204 disrupts the coding sequence of the intronless ATL80 gene, and this insertion is expected to result in a severe loss-of-function phenotype. The plants were cultivated under controlled environmental conditions with a 16-hour light/8-hour dark cycle at 23 °C. In vitro growth conditions for both WT Col-0 and the atl80 mutant were conducted following previous protocols. Approximately 25 seeds from each Arabidopsis WT Col-0 and atl80 were surface-sterilized and placed on homemade rafts made of fabric (polyester and rayon) positioned on a polypropylene mesh supported by a sponge. These rafts were then positioned inside Magenta boxes containing 30 ml of liquid MS medium and 0.5% sucrose. The seeds were subjected to stratification for four days in darkness within a cold room and were subsequently incubated for 13 days in a growth chamber with a 16-hour light/8-hour dark cycle at 23 °C. At the conclusion of the incubation period, the seeds were transferred to new Magenta boxes containing MS medium without sucrose and incubated for an additional five days. On Day 5, the rafts were removed, placed on filter paper, and exposed to continuous air dehydration for various time intervals in a laminar flow hood, including 0, 15, 30, 60, and 120 min. For experiments involving mature plants, the plants were watered every 4 days for 20 days following germination (hydrated condition). This was followed by a suspension of irrigation for 15 days (dehydration condition), and then rewatering for 1 day (rehydrated condition). Cycloheximide (SIGMA) treatment was administered to eighteen-day-old WT seedlings submerged in MS medium supplemented with 70 μm cycloheximide for 0, 15, 30, 60, and 120 min.

Determination of relative fresh weight, chlorophyll content, and reactive oxygen species following water deprivation

Physiological experiments were performed in WT and atl80 mutant plants grown in vitro in a growth chamber for 13 days and then for five more days in media without sucrose, as described above. For each dehydration time point (0, 15, 30, 60, and 120 min), 100 mg of tissue from the seedlings was collected and weighed. Excess liquid media was briefly blotted on filter paper, and samples were weighed immediately. Two replicas were collected for each time point. Fresh weight (FW) was obtained by collecting and immediately weighing the samples after removing excess water.

To quantify total chlorophyll (TCh) content, WT and atl80 seedlings for each time point were suspended in 96% ethanol and kept in the dark for 1 h before they were ground and homogenized with 2.5 ml 80% acetone in 5 mM sodium phosphate buffer (pH 7.5). The supernatant was collected and placed in a glass cuvette that was placed in a spectrophotometer to measure the absorbance at 663 and 645 nm (A645 and A663, respectively). The readings were used to determine the TCh concentration with the following equation: TCh (µg ml− 1) = 5.24A664.2 + 22.24A648.6; as previously described. To analyze reactive oxygen species, a hydrogen peroxide kit (National Diagnostics, Atlanta, GA, USA) was used. WT and atl80 mutant seedlings were ground into a powder that was dissolved in 500 µl Milli-Q water. The mixture was stirred for 20 min and then centrifuged for 15 min at 12,000 g. Aliquots (10 µl) of each sample were mixed with 90 µl of the kit solution (ferrosol iron + xylenol orange) and incubated at room temperature for 30 min. Then, the absorbance was read at 560 nm and the H2O2 concentration was calculated from a standard curve.

The data were analyzed using the STATISTICA program version 12.5. One-way and two-way ANOVA analyses were conducted, followed by a post hoc test (Tukey’s HSD) to examine differences between groups. Significance levels were determined at P < 0.05, and distinct letters were employed to indicate statistically significant differences (Tukey’s HSD, P < 0.05).

Drought transcriptome profiling

WT and atl80 seedlings from each time point were collected for control and dehydrated conditions. Two biological replicas were obtained from pools of seedlings grown under the same conditions. RNA was isolated from frozen seedlings using a RNeasy plant mini kit (Qiagen; Hilden, Germany). Labeling of total RNA, microarray processing data analysis, and normalization was performed by Oaklabs GmbH (Hennigsdorf, Germany). An ArrayXS Arabidopsis v2 (XS-5010) microarray in the Agilent 8 × 60 K format representing 30,541 Arabidopsis genes was used. TAIR9 release was used to annotate the genes, as provided by Oaklabs. Subsequently, TAIR10 and/or Araport111 were used to annotate all differentially expressed genes as defined by comparison of the atl80 mutant transcriptome with that of WT under each of the three conditions. P-values and log2 fold-changes were calculated among the groups of replicate samples. To calculate P-values, a two-tailed Student’s t-test was performed. The P-value cutoff was generally at 5% (calculated by two-tailed Student’s t-test) with an absolute log2 fold-change ≥ 1.7.

qRT‑PCR analysis

RNA was extracted from frozen seedlings using the RNeasy Plant Mini Kit (Qiagen, Hilden, Germany). cDNA synthesis was carried out with 3 µg of total RNA using the SuperScript III First-Strand Synthesis SuperMix for qRT-PCR (Invitrogen). Quantitative RT-PCR was conducted on a CFX96 Real-time System (Bio-Rad, Hercules, CA, USA) with Maxima SYBR Green qPCR Master Mix (2X) and analyzed using CFX Manager software version 3.1 (Bio-Rad, Hercules, CA, USA). The PCR conditions consisted of an initial denaturation at 95 °C for 3 min, followed by 40 cycles of 95 °C for 15 s and 63 °C for 30 s. All qRT-PCR experiments were performed with two biological replicates and three technical replicates. To normalize the data, the expression of the reference genes AT3G18780 (ACTIN2) and AT1G13440 (GAPC2) was used.

Gene Ontology (GO) term enrichment

For a comprehensive appraisal of GO term enrichment, we introduced more suitable GO classes. We categorized differentially expressed into five classes and into a class termed “other components” (a class that included non-categorized genes). The five classes were the 43 categories to assist the identification of archetypal genes from each time point are as follow: The six classes and 43 categories (in parentheses) are: gene expression (basic transcriptional machinery, maintenance of the genome/chromatin organization and remodeling, non-coding RNA, posttranscriptional regulation, TFs, tRNA and translational machinery), regulation of signaling (calcium signaling and binding, dephosphorylation, ion/biomolecule binding, lectin gene family, leucine-rich repeat, peptides, phosphorylation, protein-protein interaction, receptor-like protein kinase, ubiquitination), growth and development (senescence, cell cycle and cytoskeleton, circadian clock, hormones, modulator of growth, reproductive development, senescence, tropism/hydrotropism), metabolism (glucosinolate metabolism, lipid metabolism, nutrition, oxidation-reduction process, photosynthesis/light response/chloroplast function, primary metabolic process, secondary metabolism, endocytic pathway, transport), protection (cell wall and membrane dynamics, dehydration/hypoxia/freezing tolerance, detoxification, lipid transfer, proteases, protein folding/refolding, response to stress factors, response to stress factors-abiotic, response to stress factors-biotic, sugar and osmolytes, water transport).

cis-regulatory element enrichment on promoters of PRGs

The cis-element enrichment analysis was performed using the CisCross web service (https://plamorph.sysbio.ru/ciscross/). The analysis aimed to predict upstream regulators, and the following parameters were employed: promoter length, 1500 bp; background, ARAPORT11; DAP-seq collection, Plant Cistrome, CisCross MACS2; multiple comparisons, Benjamini-Hochberg; adjusted p-value threshold, 0.001 [19, 20].

Results

Influence of ATL80 on the early response to water deprivation stress

Stress induced by water deficit significantly hampers plant growth and development. In response to this stress, a multitude of changes in gene expression takes place. Most investigations focusing on alterations in gene expression have historically involved prolonged exposure to drought stress [21]. In our examination of ATL-type E3 ligases involved in drought stress within community databases for Arabidopsis, we identified an interesting observation. In the AtGenExpress abiotic stress microarray analysis, ATL80 exhibited a rapid and transient response to water loss, suggesting that ATL80 acts as a primary response gene [18]. The drought experiments conducted as part of the AtGenExpress project aimed to assess the impact of a 15-minute water deficit stress (WDS) exposure on gene expression. Seedlings grown in liquid media within Magenta boxes were removed and exposed to an air stream for 15 min. Subsequently, they were returned to the Magenta boxes, incubated in the growth chamber, and explants were collected after 15, 30, 60, 180, and 360 min [18]. In a bid to investigate the effect of continuous WDS on early gene expression under similar conditions, we designed an experiment (Fig. 1A). Our findings revealed a swift and transient transcriptional response of ATL80 when seedlings were subjected either to WDS for 15 min or continuous WDS. ATL80 expression peaked after 30 min of WDS and subsequently decreased after 90 min (Fig. 1B). Remarkably, no other biotic or abiotic stimuli tested elicited such a transient response (Kilian et al., 2007). Notably, ATL80 expression exhibited a substantial increase even in the absence of new protein synthesis, which aligns with the characteristics of primary response genes. This induction of ATL80 expression was notably pronounced in the presence of the translational inhibitor cycloheximide (Fig. 1C).

To contextualize the expression of ATL80, we compared it to other ATLs that exhibited early expression in response to WDS. ATL80 displayed an expression pattern reminiscent of a primary response gene. ATL15 exhibited an initial increase at 15 min but showed a minor decline thereafter. In contrast, ATL17, ATL31, ATL41, and ATL56 showed a peak of expression after 60 min, and, with the exception of ATL17, their expression initially decreased in the first 15 min after WDS before peaking (Fig. 1D). Based on these observations, we postulated that ATL80 might play a role in modulating early transcriptional responses of plants to dehydration.

Fig. 1
figure 1

Influence of ATL80 on the early response to WDS. (A) Experimental design to inspect the impact of ATL80 on the initial response to continuous WDS. Eighteen-day-old WT and atl80 seedlings, cultivated in vitro, were exposed to continuous air dehydration in a laminar flow hood. Explants were collected at 0, 15, 30, 60, and 120 min. The expression of ATL80 was also evaluated in the atl80 mutant insertional line (dashed line). The details of the experimental setup can be found in the Methods section. (B) Transcriptional response of ATL80 to WDS. For the 15 min response to WDS, data was obtained from the AtGenExpress project (dark line) [18]. The continuous WDS response was conducted as described in (A) and quantified using real-time quantitative PCR (pPCR) (gray line). (C) Effect of the translational inhibitor cycloheximide on ATL80 expression. Eighteen-day-old seedlings WT were incubated with 70 μm of cycloheximide for 0, 15, 30, 60, and 120 min and relative expression was evaluated by qPCR. (D) Evaluation of the expression levels of early responsive ATLs to WDS among the 100 A. thaliana ATLs. Genes displaying enhanced early expression in response to WDS are shown

The atl80 mutant did not exhibit noticeable disruptions in basic physiological processes associated with WDS. This was evident from the similar response observed in both wild-type (WT) and mutant lines, which manifested as a progressive decrease in fresh weight, chlorophyll content, and H2O2 production (Fig. 2A). In an attempt to identify a water stress-related phenotype in atl80, adult plants were subjected to a dehydration/rehydration experiment. Seedlings were watered every 4 days for 20 days after germination, followed by a suspension of irrigation for 15 days, and then rewatered for one day. After rehydration, it was observed that 84% of WT plants and 23% of the control atl78 survived, whereas 96% of the atl80 mutant did, suggesting a moderate effect of ATL80 on the survival of adult plants after rehydration (Fig. 2B). We also determined the number leaves in the terminal stage of senescence and H2O2 production. The number of chlorotic leaves is lower in atl80 compared to WT and atl78 plants (Fig. 2C, left panel). During dehydration, there was a progressive increase in H2O2 production observed in all three lines. However, upon rehydration, H2O2 production significantly decreased in the WT and in the atl78 mutant, and to a somewhat lesser extent in the atl80 mutant. Additionally, the atl80 mutant displayed an accelerated flowering time, indicative of a higher growth rate [22].

Fig. 2
figure 2

Impact of ATL80 on standard responses to WDS. (A) Comparison of WT (dark line) and the atl80 (gray line) mutant on relative water content, total chlorophyll content, and H2O2 content during the time course experiment in seedlings. (B) Survival rates of WT, atl78, atl80 lines after WDS. Plants were watered every 5 days for 20 days after germination, followed by suspension of irrigation for 11 days, and then rehydrated for one day. Survival percentages were recorded one day after rehydration. (C) Number of leaves in the terminal stage of senescence (chlorotic leaves) recorded after rehydration (left panel), and H2O2 content (right panel) obtained from the experiment in (B)

Time course microarray analysis of primary response genes in wild-type and the atl80 mutant after WDS

The expression of ATL80 exhibits remarkable dynamics in response to WDS, undergoing rapid changes within minutes. We reasoned that gene expression changes dependent on ATL80 might occur during these brief time intervals, and that these changes could be readily detected than alterations in protein stability regulated by the E3 ubiquitin ligase activity of ATL80. Accordingly, as an initial step in analyzing ATL80 function we conducted a time course microarray analysis experiment spanning 0, 15, 30, 60, and 120 min to assess the influence of ATL80 on the early transcriptome in response to WDS. Its impact on early gene expression in response to WDS. We employed an ArrayXS Arabidopsis chip containing 32,073 genes for this purpose. The transcriptome time-course after exposure to WDS was profiled using “in vitro” grown seedlings, and two biological replicates were included in the analysis (Fig. 1A and Methods). Principal component analysis (PCA) conducted on the raw data from all samples shows the biological variation among sample replicates and distinguishes differences between treatments and genotypes. The PCA analysis revealed that the transcriptional profiles of the replicates were highly similar, whereas significant differences emerged at various time points, primarily at 15 and 30 min after dehydration in the WT. Interestingly, it becomes evident that the atl80 mutant fails to exhibit a response to dehydration, given the similarity observed across different time points. This dissimilarity between the WT and the mutant plant underscores the important role of ATL80 in the initial response to dehydration (Fig. 3A). Differentially expressed genes (DEGs) were identified for both WT and the atl80 mutant at various time points. The results indicated that ATL80 had a significant impact on gene expression, with a striking 95% reduction in DEGs observed in the atl80 mutant compared to the WT (Fig. 3B). Notably, nearly 50% of these DEGs were detected in WT plants within the first 15 min of exposure to WDS. During this initial 15-30-minute interval, over 95% of DEGs exhibited upregulation.

Fig. 3
figure 3

Influence of ATL80 on the transcriptome in response to WDS. (A) Principal Component Analysis (PCA) of Raw Expression Data. PCA was performed on the raw expression data of all samples. Dots represent WT samples, color-coded for easy distinction, while triangles represent atl80 mutant samples. The axes represent the principal components. (B) Comparison of WT and atl80 transcriptomes between time points. The number of DEGs in WT and atl80 are presented in the left and right panels, respectively

Dynamic bursts of differentially expressed genes between WT and atl80 transcriptomes

Next, DEGs were identified through a comparative analysis of the WT and atl80 transcriptomes across various time points (Fig. 4A). Fold change values were calculated using log2-normalized data from replicate samples, with a significance threshold of P-value < 0.05 and an absolute log2 fold-change cutoff ≥ 1.7. The number and proportion of upregulated and downregulated genes varied across different time points. Specifically, higher numbers of DEGs were observed at 0 and 120 min, and distinct patterns of upregulation and downregulation were evident within 0 and 15 min, indicating that ATL80 had a significant impact on gene expression. A number of differences between WT and atl80 were detected before WDS (Fig. 4B). Furthermore, several differences between the WT and atl80 transcriptomes were detected even before the onset of WDS (Fig. 4B). Genes related to the signaling class, including receptor-like proteins and leucine-rich repeat proteins, exhibited a general trend of upregulation both before and after WDS. In contrast, genes involved in phosphorylation, protein-protein interactions, ubiquitination, and peptides showed an overall trend of upregulated expression before WDS and at 15 min, followed by downregulation at 30, 60, and 120 min post-WDS. Similarly, genes related to metabolism and the protection class displayed a similar pattern of upregulated expression before WDS and after 15 min (Fig. 4B).

Although TFs were similarly regulated before and after WDS, they exhibited upregulation before WDS and significant downregulation after WDS (Fig. 4C). In contrast, WRKY, NAC DOMAIN CONTAINING PROTEIN (NAC), and zinc finger factors generally showed upregulation, while Ethylene-Responsive Factor/APETALA2 (ERF/AP2), MADS-box (MADS), and homeobox factors were predominantly downregulated. Basic helix-loop-helix (bHLH), Basic leucine-zipper (bZIP), and MYB domain proteins (MYBs) exhibited upregulation preceding WDS but were largely downregulated after the stress (Fig. 4C). NAC, MADS, homeobox, bHLH, and bZIP transcription factors are known to play diverse roles in growth and developmental processes (Supplementary Table 1). Conversely, other types of transcription factors primarily function in regulating stress responses. Several WRKY transcription factors, such as WRKY17, WRKY18, WRKY48, WRKY50, WRKY51, and WRKY70, are involved in plant immunity) [23]. Zinc finger factors like ZAT12 and ZAT18 are associated with responses to oxidative stress and drought stress tolerance [24, 25]. ERF/AP2 TFs play crucial roles in responding to abiotic stresses. For example, DREB2A and DREB19 respond to WDS, DDF1 and DDF2 are involved in gibberellic acid-mediated stress responses, and ERF96, ORA47, DREB1A/CBF3 participate in ABA-mediated responses to various stresses [26,27,28,29,30,31]. Furthermore, eleven ERF/AP2 TFs, including ERF022, ERF025, ERF017, ERF013, ERF011, ERF054, ERF6, ERF5, ERF105, ERF104, and ERF109/RRTF1, respond rapidly to high (H)-light and are involved in a retrograde signaling pathway [32]. In addition, MYC2 plays a significant role in JA-regulated signaling for stress responses and plant growth and development, while AITR2 is involved in the feedback regulation of ABA signaling [33, 34].

Before WDS, genes in the maintenance and chromatin organization category showed upregulation. This included histone superfamily genes, DNA glycosylases, origin of replication complex protein 1B, and condensing complex subunit, among others (Fig. 4D). Notably, the expression of ARGONAUTE9 (AGO9), known for its role in silencing transposable elements, as well as several stress response genes, were upregulated at both 15 and 60 min. Additionally, CAF1D, a component of the deadenylase complex, exhibited downregulation at 15 min [35, 36].

Genes associated with growth and development also displayed upregulation prior to WDS. Two phosphatidylethanolamine-binding protein (PEBP) genes involved in reproductive development, namely FLOWERING LOCUS T (FT) at 15 and 30 min and TWIN SISTER OF FT (TSF) at 120 min, were found to be downregulated (Fig. 4D,E) [37]. Moreover, genes related to the cell cycle and cytoskeleton exhibited upregulation before WDS. The circadian clock master regulators CCA1 and LHY1 were upregulated at all time points except at 60 min; after 15 min of WDS, other regulators and clock-regulated genes, including RVE1, RVE8, PRR9, LNK2, and LNK3, displayed upregulation (Fig. 4E) [38, 39]. Additionally, genes regulated by auxin were upregulated prior to WDS, while those related to jasmonic acid (JA) showed downregulation after 30 min (Fig. 4F). Conversely, genes associated with protein folding/refolding were downregulated prior to WDS (Fig. 4G). It’s worth noting that cytosolic chaperones from the HSP90 and HSP70 families are known to regulate retrograde signaling [40].

Fig. 4
figure 4

Comparison of WT and atl80 transcriptomes in response to WDS. (A) Heat map displaying 1,938 DEGs at time points 0, 15, 30, 60, and 120 min following WDS. The number of upregulated and downregulated genes at each time point is indicated to the right. (B) GO visualization of relevant categories across five classes. Gene clustering was performed using the K-means method with the Complex Heatmap package in Rstudio. Additional details for specific terms are shown in: (C) TFs with a further breakdown of ERF/AP2s, *involved in retrograde signaling [32]; (D) posttranscriptional regulation, maintenance, and chromatin organization; growth and development category; (E) cell cycle and cytoskeleton and circadian clock master regulators; (F) hormones; (G) the protein folding/refolding. Supplementary Table 1 display a comprehensive list of the genes included in each class and category

Effect of ATL80 on transcriptional responses induced by WDS

Previous studies exploring early responses to abiotic stresses have revealed a wide array of distinct gene expression patterns in reaction to stress exposure. In our investigation, we assumed the existence of multiple transcriptional waves triggered by WDS. To examine this, we utilized fold-change values to group genes based on their fold-differences, ultimately identifying 2,156 genes that displayed differential expression in at least one of the comparisons between consecutive time points during our time-course experiment. Subsequently, we organized these WDS-responsive genes into 12 distinct clusters. The expression values of each cluster relative to the atl80 mutant were extracted from the transcriptomic data and clustered accordingly (Fig. 5A, Supplementary Table 2). Our primary objective was to ascertain whether ATL80 exerts an influence on the overall expression of WDS-responsive genes. Remarkably, the comparison of these twelve sets of clusters demonstrated that the expression patterns of WDS-responsive genes were predominantly reliant on ATL80, as the clustering patterns were disrupted in the atl80 mutant. Additionally, qPCR analysis of PRGs and other genes responsive to WDS highlighted that ATL80 played a pivotal role in regulating these genes, as evidenced by the disruption of expression in the atl80 mutant (Fig. 5B).

Fig. 5
figure 5

Patterns of rapidly induced transcription in response to WDS are abolished in atl80. (A) K-means cluster analysis of 2156 (DEGs) between time points in WT transcriptomes. DEGs log2FC values were extracted from the atl80 mutant transcriptomic data and clustered by the K-means method. Cluster 11 includes ATL80. (B) Comparison of relative expression between WT and atl80 in response to WDS, assessed by qRT-PCR analysis. CAMBP25, AP2C1, WRKY40 and ERF13 are representative genes selected from cluster 11, which contains 77 PRGs. GGL19, Histone3, FT, GSTU20, CASPL1D1, NAC102, DJC66 and HAP17.6 A represent genes from other clusters

Primary response genes among the genes responsive to WDS

Among the genes responsive to WDS, we have identified a subset of PRGs. These PRGs exhibit a transient induction of expression, peaking at 30 min following WDS. To gain insights into the regulation and function of ATL80, we focused on the subset of 73 PRGs that closely match the expression pattern of ATL80 among PRGs with this transient induction. Within this group of 73 WDS-responsive PRGs (as depicted in Fig. 6A), we performed Gene Ontology (GO) term enrichment analysis (Fig. 6B). Notably, we observed the presence of ERF/AP2 and WRKY transcription factors (TFs), with 10 members belonging to the ERF/AP2 family and 4 to the WRKY family. Among the ERF/AP2 TFs, OCTADECANOID-RESPONSIVE AP2/ERF 47 (ORA47) is known for facilitating a rapid response to environmental signals. C-REPEAT/DRE BINDING FACTOR1 (CBF1) plays a critical role in cold acclimation, while REDOX RESPONSIVE TRANSCRIPTION FACTOR1 (ERF109/RRTF1) is involve in several stress responses and integrates age and wound signals for root regeneration [41,42,43]. On the other hand, WRKY40 is primarily involved in perceiving and responding to biotic stress, WRKY46 is associated with drought stress tolerance, and WRKY53, along with WRKY18, activates early expression of several sugar-responsive genes [44, 45]. In the posttranscriptional regulation category, we found CCR4-ASSOCIATED FACTOR1 (CAF1A), CAF1D, and CAF1E, which are orthologs of the yeast deadenylase component CAF1 [36].

In the regulation of signaling category, two key genes are CALMODULIN LIKE 37 (CML37), which plays a crucial role in defense responses under stress conditions and CAMBP25, which acts as a negative effector of osmotic stress tolerance [46,47,48]. Additionally, CML39 is involved in various developmental processes and has a role in light-promoted seedling development [49, 50].

In the phosphorylation category, various genes are highlighted. LYSM-CONTAINING RECEPTOR-LIKE KINASE 5 (LYK5) is crucial for chitin perception [10, 51, 52]. AP2C1, which is closely related to PROTEIN PHOSPHATASE 2C5 and interacts with stress-induced Mitogen-Activated Protein Kinases (MAPKs) MPK3, MPK4, and MPK6. AP2C1 functions as a MAPK phosphatase, regulating ABA-induced gene activation [53]. MPK3 and MPK6 are PRGs that enhance ethylene biosynthesis in response to wounding, and MAP3K14, whose expression is rapidly induced in response to wounding, is part of two modules that rapidly respond to wounding [54].

In the hormone category, 1-AMINOCYCLOPROPANE-1-CARBOXYLIC ACID SYNTHASE 6 (ACS6) is a member of the ACC synthase family of ethylene biosynthetic enzymes. Its activation is involved in regulating stomatal density on the leaf epidermis in response to drought [55]. Additionally, F-BOX STRESS INDUCED 1 (FBS1) influences the expression of genes regulated by JA and abscisic acid [56,57,58]. Several genes involved in biotic and abiotic stresses were also share a primary response genes pattern. BYPASS1-LIKE (B1L) regulates growth and cold tolerance [59,60,61,62,63]. Besides, FLOTILLIN1 (FLOT1) and FLOT2 are genes that form membrane nanodomains, potentially assisting in communication between the plasma membrane and the extracellular environment [64].

Convergence among retrograde signaling pathways and the influence of ATL80 on gene expression were deduced in our previous analysis (Fig. 4C and G). To provide support to such assumption, we cross-referenced data associated with retrograde signaling pathways against the early ATL80 transcriptome. We employ a set of 39 genes that form a core response module for retrograde signaling, along with a set of 64 reference genes analyzed in an experiment on conditional acclimation following six hours of H-light treatment [65, 66]. We investigated whether among them they were WDS-responsive genes. Thirteen genes were identified within the 39 core response genes, while 29 were detected in the 64 genes responsive to H-light (Fig. 6C and D) (Supplementary Table 3). The essential regulatory role of ATL80 in early expression of several these genes is evident, as demonstrated by the observed alteration in expression in the atl80 mutant (Fig. 6C and D; left graphs).

Among the genes involved in retrograde signaling pathways, six TFs exhibited increased expression early in response to WDS. Notably, ERD109/RRTF1 has roles in both abiotic and biotic stress responses (Fig. 6D, lower graphs) [43]. Additionally, MYC2 serves as a key regulator of the JA signaling, impacting stress responses and plant growth and development (Fig. 6D, middle graphs) [67]. ZAT10 and ZAT12, are implicated in diverse stress responses, including drought, H-light, salt, cold, oxidative stress and osmotic stress (Fig. 6D, middle graphs) [68]. In the chromatin organization category, HIS1-3 encode a histone variant whose expression is induced by drought, salt and ABA (Fig. 6C, upper graphs) [69].

Within the calcium signaling and binding category, two pertinent genes are CML24 whose expression is enhanced in response to environmental changes, hydrogen peroxide, ABA, and IAA (Fig. 6C, upper graphs) [70]. Furthermore, BAP1, in conjunction with BON1 control plant growth homeostasis. BON1 proteins regulate global stress responses, including calcium signaling (Fig. 6D, lower graphs) [71]. Remarkably, COR47 and ATDI8, members of the dehydrin family that share highly conserved structural features (Fig. 6D, upper and lower graphs respectively), along with the RING-H2 E3 ubiquitin ligase ATL78 that has a role in adaptation to drought tolerance and FLOT3 found in membrane nanodomains, are also present among the cross-referenced genes (Fig. 6C, middle graphs; Fig. 6D, second graphs respectively) [17, 64, 72].

Fig. 6
figure 6

Primary response genes pattern in response to WDS. (A) Heat map displays the expression data of 73 PRGs in response to WDS, along with their respective expression values extracted from transcriptomic data of the atl80 mutant. The 73 genes were specifically chosen from Cluster 11, which comprises a total of 129 genes, due to their pronounced similarity to the expression pattern observed in ATL80. (B) GO analysis of the 73 PRGs. (C) Cross-reference of transcriptomic data related to core response module of retrograde signaling pathway genes with the ATL80 transcriptome. The mean expression of 13 genes is illustrated in three graphs based on expression intensity and is compared to their respective expression levels in the atl80 mutant (Supplementary Table 3). (D) Cross-reference of transcriptomic data related to genes responding to H-light acclimation. The mean expression of 29 genes is depicted in five graphs based on expression intensity and is compared to their corresponding expression in the atl80 mutant (Supplementary Table 3)

Enrichment for calmodulin-binding transcription activator (CAMTA) TFs binding sites in the PRG promoters

TFs play a pivotal role in the initial control of gene expression by binding to specific DNA sequences, typically located in gene promoter regions. These DNA sequences in the promoters are specific cis-regulatory elements (CREs), which, on a genome-wide scale, collectively form the cistrome. To examine the potential association between TFs to cis-regulatory elements using the CisCross tool which predicts upstream regulators [19]. We evaluated the enrichment of cis-elements on the 73 set of PRGs and compared this the enrichment with cluster 11, the originally cluster that included 129 genes, and to 2027 WDS-responding genes out of this cluster (Fig. 7).

We then compare the enrichment of TFs to cis-regulatory elements between PRGs and those that do not exhibit a primary response gene pattern (Fig. 7). Our analysis revealed the presence of TFs exclusively enriched in PRGs, indicating their potential involvement in the regulation of these genes. Particularly noteworthy was the exclusive enrichment of binding sites for CALMODULIN-BINDING TRANSCRIPTION ACTIVATOR (CAMTA) TFs. CAMTA TFs are known to play a crucial role in the rapid response to a wide range of stresses. Specifically, CAMTA1 has been implicated as a regulator of drought stress, while CAMTA5 is rapidly activated in response to an immediate decrease in temperature, where it activates the expression of DREB1 genes [73,74,75]. Our analysis also revealed the enrichment of several WRKY TFs when all genes from cluster 11 were included. Furthermore, we observed a significant assortment and enrichment of TFs among other WDS-responsive genes, highlighting a diversity of TFs enriched among the genes in the other clusters (bZIP, HD, MYB, NAC) (Fig. 7).

Fig. 7
figure 7

Enrichment analysis of cis-regulatory elements on PRGs. Heat map displaying the results of the cis-regulatory elements enrichment analysis of PRGs. This analysis was performed using the CisCross web service with default parameters. TF types are listed along with the family and representative sequence features for each type. The 73 PRGs were analyzed independently or together with the 129 genes in Cluster 11. “Other” refers to genes from the rest of the clusters. The respective logos were downloaded from the website http://neomorph.salk.edu/dap_web/pages/index.php [20]. The analysis was conducted independently for the 73 PRGs, collectively for the 129 genes in Cluster 11 that include the 73 PRGs, and as the group for genes from the rest of the clusters labeled as “Other”

Discussion

Primary response genes serve as key regulators in the transition of short-lived signals into prolonged responses. Their expression is rapidly triggered within minutes of perceiving stimuli, reaching its peak at around 60 min before gradually declining. The swift transcriptional responses to water deficit stress at brief time intervals imply that these responses may play a pivotal role in enabling plants to adapt to such challenging conditions. Despite ATL80 being characterized as an E3 ubiquitin ligase, typically associated with post-translational regulation, our hypothesis posits that, as a PRG, it might have a significant influence on gene expression. Our findings provide evidence that ATL80 indeed plays a critical role in modulating early responses to WDS, as demonstrated by the observation that gene expression within the 0-120 min time frame was profoundly disrupted in atl80 mutant plants. This assertion is further supported by the presence of numerous DEGs between the WT and atl80 mutant, many of which are known to have substantial roles in stress responses or response to WDS. Notable examples include several WRKY genes, ZAT12, ZAT18, DREB2A, DREB19, DDF1, DDF2, ERF96, ORA47, and DREB1A/CBF3 (Fig. 4). For instance, WRKY46 has been implicated in drought tolerance, as demonstrated by loss-of-function mutants that revealed its importance in this context [44, 45]. Additionally, the ectopic expression of ERF019, DREB2A, or DDF1 has been shown to confer significant drought stress tolerance and activate various genes involved in the response to WDS [26,27,28] (Fig. 4).

The analysis of transcriptome data at various time points following WDS revealed a significant reduction (up to 95%) in the number of DEGs in the atl80 mutant compared to the WT, indicating that ATL80 plays a pivotal role in modulating gene expression during the initial two hours after WDS (see Fig. 3b). It is conceivable that ATL80 may also exhibit peaks or additional expression patterns at extended time intervals in response to WDS (beyond several hours). For example, ZAT12 and WRKY40 transcription factors have been shown to display three distinct expression patterns at seconds, minutes, and hours in response to stress [76]. The presence of such expression patterns suggests that bursts of regulatory gene expression exert a significant influence on the transcriptomic landscape, contributing to the shaping and fine-tuning of the stress response.

In the absence of stress, ATL80 exerts an influence on the expression of several classes of genes. These genes span various functional categories, with upregulated genes related to the maintenance of genome/chromatin organization, remodeling, growth and development, cell cycle, cytoskeleton, and genes involved in auxin, brassinosteroids, cytokinin, gibberellic acid, and ethylene metabolism. Conversely, genes categorized under protein folding/refolding showed a downregulation, including notable chaperones like HSP90.1, DNAJ-domain proteins, HSP70, HSP17, and HSP20-like chaperones. These molecular chaperones play pivotal roles in diverse growth and plant development processes, facilitating the proper folding of polypeptide chains to prevent protein aggregation during stress responses [77,78,79]. Following exposure to WDS, the expression of these genes was no longer differential, suggesting a rapid transcriptome reprogramming in response to stress. Similar swift responses have been previously documented, where various genes displayed altered expression patterns within minutes or seconds in response to abiotic stresses [80, 81].

The rapid response to WDS results in the activation of calcium signaling. Calcium has also been associated with rapid, widespread signaling activity across the plant [82]. Among primary response genes to WDS included in the calcium signaling and calcium binding category genes, are regulators of the ABA pathway. Additionally, CML37 is a calmodulin-like gene that functions as a positive regulator of ABA [47]. A search for general stress response regulators identified ORA47, a jasmonic acid-induced gene, and JAZ1 which acts as a repressor may serve as a transcriptional regulator for a central hub controlling the general stress response [30]. The finding that several genes encoding JA precursor of biosynthetic proteins and signaling cascade components that are among the primary response genes identified here support such a mechanism that entails the interaction between the JA metabolism and the reprograming of stress responses.

In the subset of 73 PRGs, there are several regulatory genes implicated in stress responses, such as ERF/AP2 and WRKY TFs, calcium signaling genes, MAP kinases, and various signaling peptides associated with stress (Fig. 6). The expression of these specific genes may share common cis-regulatory elements, as indicated by the analysis of the promoters of these 73 PRGs, which predicts an enrichment of binding sites for CAMTA1 and CAMTA5. CAMTA TFs play a crucial role in rapidly responding to various stresses. Notably, CAMTA family genes have been studied in various plant species. CAMTA1 has been identified as a key regulator of drought stress in A. thaliana [73, 74]. Transcriptome data from Heimia myrtifolia revealed three upregulated CAMTA genes in response to WDS [83]. Similarly, Phyllostachys edulis encode CAMTA genes associated with drought stress, with promoter analysis indicating the presence of stress-related cis-elements [84]. Camellia sinensis displayed distinctive expression patterns of CAMTA genes across organs and diverse stress conditions [85]. Notably, studies in Quinoa have identified the role of CqCAMTA03 in enhancing drought stress tolerance [86]. In wheat, TaCAMTA1b-B.1 emerged as a participant in the drought stress response during the seedling stage [87]. Additionally, specific CAMTA genes in soybeans were identified as negative regulators influencing both development and responses to drought stress [88]. Additionally, it is worth noting that among the 73 PRGs are CAF1A, CAF1D, and CAF1E, which may be components of the deadenylase complex. These findings suggest a potential role in mediating the transient half-life of rapidly induced mRNAs, including PRGs, in response to WDS [36].

Among DEGs, 52 ERF/AP2 TFs were identified, with 44 of them showing gene expression fluctuations within 120 minutes after subjecting the plants to a H-light intensity. Previously, 19 ERF/AP2 TFs that responded rapidly when transitioning low-light-acclimated A. thaliana plants to H-light conditions had been identified. Notably, 11 of these ERF/AP2 TFs are among the DEGs. Upon exposure to H-light conditions, there is a rapid communication pathway from the chloroplast to the nucleus, known as retrograde signaling. This pathway swiftly conveys information that impacts gene expression in response to developmental signals and environmental stresses [32]. The rapid induction of a similar set of ERF/AP2s in response to both WDS and H-light suggests a convergence of signaling pathways (Fig. 8). GENOMES UNCOUPLED1 (GUN1), a protein localized in the chloroplast with a central role in retrograde signaling, interacts with ClpC1, a chloroplast chaperone associated with the protein import machinery. In a double gun1clpc1 mutant, the expression of members from the HSP90 and HSP70 families is significantly increased (Fig. 8). This finding reveals a role for cytosolic HSP90-HSP70 chaperone complexes in retrograde signaling [40]. Genes encoding HSP90.1 and HSP70 chaperones exhibited downregulation in the atl80 mutant prior to WD, suggesting a convergence of both signaling pathways. The early expression of several core module genes, identified through the analysis of transcriptional responses in mutant lines compared to WT lines across six different experiments impacting retrograde signaling, as well as a set of reference signaling genes derived from an independent experiment on the initiation of H-light treatment, substantiates the convergence of these pathways in response to WDS (Fig. 6C and D) [65, 66]. Another pathway with a role in retrograde signaling, responding to drought and H-light stress, is the SAL1-nucleotide 3’-phosphoadenosine 5’-phosphate (PAP) signaling pathway. Drought and H-light stress inhibit SAL1 activity, leading to an increase in PAP accumulation within the chloroplasts. PAP is proposed to serve as a retrograde signal that induces the expression of nuclear genes through a yet-unknown mechanism [89]. These observations highlight ATL80 as a potential candidate for bridging various retrograde signaling pathways (Fig. 8).

Our results suggest that ATL80 has a rapid and significant impact on gene expression, as evidenced by a remarkable reduction in DEGs between the time points 0 to 120 min during the time course experiment when comparing the atl80 mutant to the WT. A similar rapid impact on gene expression is observed in response to H-light intensity. The expression of ERF/AP2 TFs and MITOGEN-ACTIVATED PROTEIN KINASE6 (MPK6) genes is highly induced within minutes. Notably, this rapid response is disrupted in the triose phosphate/phosphate translocator 1 (tpt1) and tpt2 mutant lines, implying that TPT plays a crucial role in early transcriptional responses to H-light. TPT facilitates the translocation of triose phosphate from the chloroplast, constituting a retrograde signaling pathway [32]. Furthermore, our findings suggest that components of the retrograde signaling pathway, which regulate transcriptional responses within minutes, consist of proteins other than TFs. This hints at the possibility that the E3 ubiquitin ligase ATL80 may target one or more of these components.

Fig. 8
figure 8

Schematic of the early ATL80-dependent response to WDS and its presumed function in retrograde signaling pathways. The ATL80 response to WDS is depicted in two columns: one representing the condition before stress (0 min, light colors), and the other showing the response after 15 to 120 min following stress (darker colors). Upregulation and downregulation events are indicated by yellow and violet arrows, respectively. GO classes, categories, or individual genes are denoted. Below this representation, we illustrate the expression patterns of ATL80 along with other PRGs. We postulate that CAMTA TFs are involved in regulating the expression of this set of PRGs. On the right side of the figure, we outline the basic retrograde signaling pathways. The pathway triggered by H-light exposure involves the translocation of triose phosphate, a product of the Calvin-Benson cycle, from the chloroplast through a triose phosphate translocator (TPT). This process activates the kinase MPK6 and ERF/AP2 TFs. The GUN1–HSC70 complex initiates signaling, which is mediated by cytosolic chaperones HSP90 and HSP70. The PAP-SAL1 signaling pathway responds to drought and H-light stress [90]. We have encircled components that were both detected as DEGs in our ATL80 transcriptome and are involved in retrograde signaling with a cloud shape (ERF/AP2 and proteins associated with folding/refolding/HSP90/HSP70). We hypothesize that one or more components of the retrograde signaling pathways (highlighted in light gray) serve as targets for degradation by ATL80, indicated by the arrow pointing toward the shaded region, functioning as an alarm signal

Conclusion

To maintain cellular homeostasis under stressful conditions, transcriptional reprogramming plays a critical role. Reprogramming of the transcriptome can happen within seconds following a cellular stressor (i.e., WDS), initiating a rapid response that can activate both well-known and previously unrecognized proteins and pathways. Chloroplasts are capable of sensing drought stress and generating retrograde signals that subsequently regulate the cellular stress response. Several components of retrograde signaling serve the function of modulating transcription, such as TPT. ATL80 is an E3 ligase primarily localized to the plasma membrane [22]. Our results suggest that ATL80 has a positive impact on gene expression when responding to WDS. We propose that ATL80 plays a role in the early response to WDS by reducing the half-life of components within retrograde signaling pathways, acting as an E3 ligase. Furthermore, it is possible that ATL80 might exhibit varying peaks or show additional expression patterns at extended time intervals following exposure to WDS (lasting for more than a few hours). For instance, certain PRG TFs exhibit three distinct expression patterns, which occur within seconds, minutes, and hours in response to stress [76]. The presence of this expression pattern indicates that bursts in regulatory gene activity impact the transcriptomic landscape of the plant, shaping and fine-tuning the stress response. This, in turn, modulates distinct retrograde signaling pathways, serving as an early alarm signal to mitigate damage during the initial stages of the stress response (Fig. 8).

Data availability

The authors declare that all data supporting the findings in this study are available within the paper and Supplementary information. All data are available upon request to the corresponding author.

Data availability

The authors declare that all data supporting the findings in this study are available within the paper and Supplementary information. All data are available upon request to the corresponding author.

Abbreviations

PRG:

Primary response gene

WDS:

Water deprivation stress

TF:

Transcription factor

References

  1. Fowler T, Sen R, Roy AL. Regulation of primary response genes. Mol Cell. 2011;44(3):348–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Jariel-Encontre I, Salvat C, Steff A-M, Pariat M, Acquaviva C, Furstoss O, Piechaczyk M. Complex mechanisms for c-fos and c-jun degradation. Mol Biol Rep. 1997;24:51–6.

    Article  CAS  PubMed  Google Scholar 

  3. Sheng M, Greenberg ME. The regulation and function of c-fos and other immediate early genes in the nervous system. Neuron. 1990;4(4):477–85.

    Article  CAS  PubMed  Google Scholar 

  4. Van Riggelen J, Buchwalter G, Soto U, Arce JD-C, Zur Hausen H, Wasylyk B, Rösl F. Loss of net as repressor leads to constitutive increased c-fos transcription in cervical cancer cells. J Biol Chem. 2005;280(5):3286–94.

    Article  PubMed  Google Scholar 

  5. Ang ACH, Østergaard L. Save your TIRs–more to auxin than meets the eye. New Phytol. 2023;238:971–6.

    Article  PubMed  Google Scholar 

  6. Li B, Meng X, Shan L, He P. Transcriptional regulation of pattern-triggered immunity in plants. Cell Host Microbe. 2016;19(5):641–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Jiménez-López D, Muñóz-Belman F, González-Prieto JM, Aguilar-Hernández V, Guzmán P. Repertoire of plant RING E3 ubiquitin ligases revisited: new groups counting gene families and single genes. PLoS ONE. 2018;13(8):e0203442.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Guo T, Kong F, Burton C, Scaglione S, Beagles B, Ray J, Ramonell KM. The E3 ubiquitin ligase ATL9 affects expression of Defense related genes, cell death and Callose Deposition in response to fungal infection. Pathogens. 2022;11(1):68.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Kong F, Guo T, Ramonell KM. Arabidopsis Toxicos en Levadura 12 (ATL12): a gene involved in Chitin-Induced, hormone-related and NADPH oxidase-mediated defense responses. J Fungi. 2021;7(10):883.

    Article  CAS  Google Scholar 

  10. Liu X, Zhou Y, Du M, Liang X, Fan F, Huang G, Zou Y, Bai J, Lu D. The calcium-dependent protein kinase CPK28 is targeted by the ubiquitin ligases ATL31 and ATL6 for proteasome-mediated degradation to fine-tune immune signaling in Arabidopsis. Plant Cell. 2022;34(1):679–97.

    Article  PubMed  Google Scholar 

  11. Ma X, Claus LA, Leslie ME, Tao K, Wu Z, Liu J, Yu X, Li B, Zhou J, Savatin DV. Ligand-induced monoubiquitination of BIK1 regulates plant immunity. Nature. 2020;581(7807):199–203.

    Article  CAS  PubMed  PubMed Central  ADS  Google Scholar 

  12. Yang R, Wang T, Shi W, Li S, Liu Z, Wang J, Yang Y. E3 ubiquitin ligase ATL61 acts as a positive regulator in abscisic acid mediated drought response in Arabidopsis. Biochem Biophys Res Commun. 2020;528(2):292–8.

    Article  CAS  PubMed  Google Scholar 

  13. Du B, Nie N, Sun S, Hu Y, Bai Y, He S, Zhao N, Liu Q, Zhai H. A novel sweetpotato RING-H2 type E3 ubiquitin ligase gene IbATL38 enhances salt tolerance in transgenic Arabidopsis. Plant Sci. 2021;304:110802.

    Article  CAS  PubMed  Google Scholar 

  14. Cui LH, Min HJ, Yu SG, Byun MY, Oh TR, Lee A, Yang HW, Kim WT. OsATL38 mediates mono-ubiquitination of the 14-3-3 protein OsGF14d and negatively regulates the cold stress response in rice. J Exp Bot. 2022;73(1):307–23.

    Article  CAS  PubMed  Google Scholar 

  15. Kim SJ, Kim WT. Suppression of Arabidopsis RING E3 ubiquitin ligase AtATL78 increases tolerance to cold stress and decreases tolerance to drought stress. FEBS Lett. 2013;587(16):2584–90.

    Article  CAS  PubMed  Google Scholar 

  16. Song J, Xing Y, Munir S, Yu C, Song L, Li H, Wang T, Ye Z. An ATL78-Like RING-H2 finger protein confers abiotic stress tolerance through interacting with RAV2 and CSN5B in tomato. Front Plant Sci. 2016;7:1305.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Jiménez-Morales E, Aguilar‐Hernández V, Aguilar‐Henonin L, Guzmán P. Molecular basis for neofunctionalization of duplicated E3 ubiquitin ligases underlying adaptation to drought tolerance in Arabidopsis thaliana. Plant J. 2020;104(2):474–92.

    Article  PubMed  Google Scholar 

  18. Kilian J, Whitehead D, Horak J, Wanke D, Weinl S, Batistic O, D’Angelo C, Bornberg-Bauer E, Kudla J, Harter K. The AtGenExpress global stress expression data set: protocols, evaluation and model data analysis of UV‐B light, drought and cold stress responses. Plant J. 2007;50(2):347–63.

    Article  CAS  PubMed  Google Scholar 

  19. Lavrekha VV, Levitsky VG, Tsukanov AV, Bogomolov AG, Grigorovich DA, Omelyanchuk N, Ubogoeva EV, Zemlyanskaya EV, Mironova V. CisCross: a gene list enrichment analysis to predict upstream regulators in Arabidopsis thaliana. Front Plant Sci. 2022;13:942710.

    Article  PubMed  PubMed Central  Google Scholar 

  20. O’Malley RC, Huang S-sC, Song L, Lewsey MG, Bartlett A, Nery JR, Galli M, Gallavotti A, Ecker JR. Cistrome and epicistrome features shape the regulatory DNA landscape. Cell. 2016;165(5):1280–92.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Dubois M, Inzé D. Plant growth under suboptimal water conditions: early responses and methods to study them. J Exp Bot. 2020;71(5):1706–22.

    Article  CAS  PubMed  Google Scholar 

  22. Suh JY, Kim WT. Arabidopsis RING E3 ubiquitin ligase AtATL80 is negatively involved in phosphate mobilization and cold stress response in sufficient phosphate growth conditions. Biochem Biophys Res Commun. 2015;463(4):793–9.

    Article  CAS  PubMed  Google Scholar 

  23. Wani SH, Anand S, Singh B, Bohra A, Joshi R. WRKY transcription factors and plant defense responses: latest discoveries and future prospects. Plant Cell Rep. 2021;40:1071–85.

    Article  CAS  PubMed  Google Scholar 

  24. Yin M, Wang Y, Zhang L, Li J, Quan W, Yang L, Wang Q, Chan Z. The Arabidopsis Cys2/His2 zinc finger transcription factor ZAT18 is a positive regulator of plant tolerance to drought stress. J Exp Bot. 2017;68(11):2991–3005.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Davletova S, Rizhsky L, Liang H, Shengqiang Z, Oliver DJ, Coutu J, Shulaev V, Schlauch K, Mittler R. Cytosolic ascorbate peroxidase 1 is a central component of the reactive oxygen gene network of Arabidopsis. Plant Cell. 2005;17(1):268–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Sakuma Y, Maruyama K, Osakabe Y, Qin F, Seki M, Shinozaki K, Yamaguchi-Shinozaki K. Functional analysis of an Arabidopsis transcription factor, DREB2A, involved in drought-responsive gene expression. Plant Cell. 2006;18(5):1292–309.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Scarpeci TE, Frea VS, Zanor MI, Valle EM. Overexpression of AtERF019 delays plant growth and senescence, and improves drought tolerance in Arabidopsis. J Exp Bot. 2017;68(3):673–85.

    CAS  PubMed  Google Scholar 

  28. Kang H-G, Kim J, Kim B, Jeong H, Choi SH, Kim EK, Lee H-Y, Lim PO. Overexpression of FTL1/DDF1, an AP2 transcription factor, enhances tolerance to cold, drought, and heat stresses in Arabidopsis thaliana. Plant Sci. 2011;180(4):634–41.

    Article  CAS  PubMed  Google Scholar 

  29. Wang X, Liu S, Tian H, Wang S, Chen J-G. The small ethylene response factor ERF96 is involved in the regulation of the abscisic acid response in Arabidopsis. Front Plant Sci. 2015;6:1064.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Zeng L, Chen H, Wang Y, Hicks D, Ke H, Pruneda-Paz J, Dehesh K. ORA47 is a transcriptional regulator of a general stress response hub. Plant J. 2022;110(2):562–71.

    Article  CAS  PubMed  Google Scholar 

  31. Ma T, Wang S, Sun C, Tian J, Guo H, Cui S, Zhao H, Arabidopsis LFR. A SWI/SNF complex component, interacts with ICE1 and activates ICE1 and CBF3 expression in cold acclimation. Front Plant Sci. 2023;14:1097158.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Vogel MO, Moore M, König K, Pecher P, Alsharafa K, Lee J, Dietz K-J. Fast retrograde signaling in response to high light involves metabolite export, MITOGEN-ACTIVATED PROTEIN KINASE6, and AP2/ERF transcription factors in Arabidopsis. Plant Cell. 2014;26(3):1151–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Dombrecht B, Xue GP, Sprague SJ, Kirkegaard JA, Ross JJ, Reid JB, Fitt GP, Sewelam N, Schenk PM, Manners JM. MYC2 differentially modulates diverse jasmonate-dependent functions in Arabidopsis. Plant Cell. 2007;19(7):2225–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Tian H, Chen S, Yang W, Wang T, Zheng K, Wang Y, Cheng Y, Zhang N, Liu S, Li D. A novel family of transcription factors conserved in angiosperms is required for ABA signalling. Plant Cell Environ. 2017;40(12):2958–71.

    Article  CAS  PubMed  Google Scholar 

  35. Olmedo-Monfil V, Durán-Figueroa N, Arteaga-Vázquez M, Demesa-Arévalo E, Autran D, Grimanelli D, Slotkin RK, Martienssen RA, Vielle-Calzada J-P. Control of female gamete formation by a small RNA pathway in Arabidopsis. Nature. 2010;464(7288):628–32.

    Article  CAS  PubMed  PubMed Central  ADS  Google Scholar 

  36. Walley JW, Kelley DR, Nestorova G, Hirschberg DL, Dehesh K. Arabidopsis deadenylases AtCAF1a and AtCAF1b play overlapping and distinct roles in mediating environmental stress responses. Plant Physiol. 2010;152(2):866–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Karlgren A, Gyllenstrand N, Källman T, Sundström JF, Moore D, Lascoux M, Lagercrantz U. Evolution of the PEBP gene family in plants: functional diversification in seed plant evolution. Plant Physiol. 2011;156(4):1967–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Mizoguchi T, Wheatley K, Hanzawa Y, Wright L, Mizoguchi M, Song H-R, Carré IA, Coupland G. LHY and CCA1 are partially redundant genes required to maintain circadian rhythms in Arabidopsis. Dev Cell. 2002;2(5):629–41.

    Article  CAS  PubMed  Google Scholar 

  39. Xie Q, Wang P, Liu X, Yuan L, Wang L, Zhang C, Li Y, Xing H, Zhi L, Yue Z. LNK1 and LNK2 are transcriptional coactivators in the Arabidopsis circadian oscillator. Plant Cell. 2014;26(7):2843–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Wu G-Z, Meyer EH, Richter AS, Schuster M, Ling Q, Schöttler MA, Walther D, Zoschke R, Grimm B, Jarvis RP. Control of retrograde signalling by protein import and cytosolic folding stress. Nat Plants. 2019;5(5):525–38.

    Article  CAS  PubMed  Google Scholar 

  41. Dong X, Yan Y, Jiang B, Shi Y, Jia Y, Cheng J, Shi Y, Kang J, Li H, Zhang D. The cold response regulator CBF1 promotes Arabidopsis hypocotyl growth at ambient temperatures. EMBO J. 2020;39(13):e103630.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Ye B-B, Shang G-D, Pan Y, Xu Z-G, Zhou C-M, Mao Y-B, Bao N, Sun L, Xu T, Wang J-W. AP2/ERF transcription factors integrate age and wound signals for root regeneration. Plant Cell. 2020;32(1):226–41.

    Article  CAS  PubMed  Google Scholar 

  43. Yang C-L, Huang Y-T, Schmidt W, Klein P, Chan M-T, Pan I. Ethylene response Factor109 attunes immunity, photosynthesis, and iron homeostasis in arabidopsis leaves. Front Plant Sci. 2022;13:492.

    Google Scholar 

  44. Chen J, Nolan TM, Ye H, Zhang M, Tong H, Xin P, Chu J, Chu C, Li Z, Yin Y. Arabidopsis WRKY46, WRKY54, and WRKY70 transcription factors are involved in brassinosteroid-regulated plant growth and drought responses. Plant Cell. 2017;29(6):1425–39.

    CAS  PubMed  PubMed Central  Google Scholar 

  45. Chen Q, Xu X, Xu D, Zhang H, Zhang C, Li G. WRKY18 and WRKY53 coordinate with HISTONE ACETYLTRANSFERASE1 to regulate rapid responses to sugar. Plant Physiol. 2019;180(4):2212–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Choudhary A, Senthil-Kumar M. Drought attenuates plant defence against bacterial pathogens by suppressing the expression of CBP60g/SARD1 during combined stress. Plant Cell Environ. 2022;45(4):1127–45.

    Article  CAS  PubMed  Google Scholar 

  47. Scholz SS, Reichelt M, Vadassery J, Mithöfer A. Calmodulin-like protein CML37 is a positive regulator of ABA during drought stress in Arabidopsis. Plant Signal Behav. 2015;10(6):e1011951.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Perruc E, Charpenteau M, Ramirez BC, Jauneau A, Galaud JP, Ranjeva R, Ranty B. A novel calmodulin-binding protein functions as a negative regulator of osmotic stress tolerance in Arabidopsis thaliana seedlings. Plant J. 2004;38(3):410–20.

    Article  CAS  PubMed  Google Scholar 

  49. Bender KW, Rosenbaum DM, Vanderbeld B, Ubaid M, Snedden WA. The A rabidopsis calmodulin-like protein, CML 39, functions during early seedling establishment. Plant J. 2013;76(4):634–47.

    Article  CAS  PubMed  Google Scholar 

  50. Zhou YP, Wu JH, Xiao WH, Chen W, Chen QH, Fan T, Xie CP, Tian C-E. Arabidopsis IQM4, a novel calmodulin-binding protein, is involved with seed dormancy and germination in Arabidopsis. Front Plant Sci. 2018;9:721.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Huang C, Yan Y, Zhao H, Ye Y, Cao Y. Arabidopsis CPK5 phosphorylates the chitin receptor LYK5 to regulate plant innate immunity. Front Plant Sci. 2020;11:702.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Beaugelin I, Chevalier A, D’alessandro S, Ksas B, Novák O, Strnad M, Forzani C, Hirt H, Havaux M, Monnet F. OXI1 and DAD regulate light-induced cell death antagonistically through jasmonate and salicylate levels. Plant Physiol. 2019;180(3):1691–708.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Brock AK, Willmann R, Kolb D, Grefen L, Lajunen HM, Bethke G, Lee J, Nürnberger T, Gust AA. The Arabidopsis mitogen-activated protein kinase phosphatase PP2C5 affects seed germination, stomatal aperture, and abscisic acid-inducible gene expression. Plant Physiol. 2010;153(3):1098–111.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Sözen C, Schenk ST, Boudsocq M, Chardin C, Almeida-Trapp M, Krapp A, Hirt H, Mithöfer A, Colcombet J. Wounding and insect feeding trigger two independent MAPK pathways with distinct regulation and kinetics. Plant Cell. 2020;32(6):1988–2003.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Jia M-z, Liu L-y, Geng C, Jiang J. Activation of 1-aminocyclopropane-1-carboxylic acid synthases sets stomatal density and clustered ratio on leaf epidermis of Arabidopsis in response to drought. Front Plant Sci. 2021;12:758785.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Hasegawa Y, Huarancca Reyes T, Uemura T, Baral A, Fujimaki A, Luo Y, Morita Y, Saeki Y, Maekawa S, Yasuda S. The TGN/EE SNARE protein SYP61 and the ubiquitin ligase ATL31 cooperatively regulate plant responses to carbon/nitrogen conditions in Arabidopsis. Plant Cell. 2022;34(4):1354–74.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Cho SK, Ryu MY, Song C, Kwak JM, Kim WT. Arabidopsis PUB22 and PUB23 are homologous U-Box E3 ubiquitin ligases that play combinatory roles in response to drought stress. Plant Cell. 2008;20(7):1899–914.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Gonzalez LE, Keller K, Chan KX, Gessel MM, Thines BC. Transcriptome analysis uncovers Arabidopsis F-BOX STRESS INDUCED 1 as a regulator of jasmonic acid and abscisic acid stress gene expression. BMC Genomics. 2017;18(1):1–15.

    Article  Google Scholar 

  59. Dutta A, Choudhary P, Gupta-Bouder P, Chatterjee S, Liu P-P, Klessig DF, Raina R. Arabidopsis small defense-associated protein 1 modulates pathogen defense and tolerance to oxidative stress. Front Plant Sci. 2020;11:703.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Chen T, Zhang W, Yang G, Chen J-H, Chen B-X, Sun R, Zhang H, An L-Z. TRANSTHYRETIN-LIKE and BYPASS1-LIKE co-regulate growth and cold tolerance in Arabidopsis. BMC Plant Biol. 2020;20(1):1–11.

    Article  CAS  Google Scholar 

  61. Cai C, Wang W, Ye S, Zhang Z, Ding W, Xiang M, Wu C, Zhu Q. Overexpression of a novel Arabidopsis gene SUPA leads to various morphological and abiotic stress tolerance alternations in Arabidopsis and Poplar. Front Plant Sci. 2020;11:560985.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Ge X, Li G-J, Wang S-B, Zhu H, Zhu T, Wang X, Xia Y. AtNUDT7, a negative regulator of basal immunity in Arabidopsis, modulates two distinct defense response pathways and is involved in maintaining redox homeostasis. Plant Physiol. 2007;145(1):204–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Pereira Mendes M, Hickman R, Van Verk MC, Nieuwendijk NM, Reinstädler A, Panstruga R, Pieterse CM, Van Wees SC. A family of pathogen-induced cysteine-rich transmembrane proteins is involved in plant disease resistance. Planta. 2021;253:1–17.

    Article  Google Scholar 

  64. Daněk M, Angelini J, Malínská K, Andrejch J, Amlerová Z, Kocourková D, Brouzdová J, Valentová O, Martinec J, Petrášek J. Cell wall contributes to the stability of plasma membrane nanodomain organization of Arabidopsis thaliana FLOTILLIN2 and HYPERSENSITIVE INDUCED REACTION1 proteins. Plant J. 2020;101(3):619–36.

    Article  PubMed  Google Scholar 

  65. Gläßer C, Haberer G, Finkemeier I, Pfannschmidt T, Kleine T, Leister D, Dietz K-J, Häusler RE, Grimm B, Mayer KFX. Meta-analysis of retrograde signaling in Arabidopsis thaliana reveals a core module of genes embedded in complex cellular signaling networks. Mol Plant. 2014;7(7):1167–90.

    Article  PubMed  Google Scholar 

  66. Alsharafa K, Vogel MO, Oelze M-L, Moore M, Stingl N, König K, Friedman H, Mueller MJ, Dietz K-J. Kinetics of retrograde signalling initiation in the high light response of Arabidopsis thaliana. Philosophical Trans Royal Soc B: Biol Sci. 2014;369(1640):20130424.

    Article  Google Scholar 

  67. Kazan K, Manners JM. MYC2: the master in action. Mol Plant. 2013;6(3):686–703.

    Article  CAS  PubMed  Google Scholar 

  68. Liu Y, Khan AR, Gan Y. C2H2 zinc finger proteins response to abiotic stress in plants. Int J Mol Sci. 2022;23(5):2730.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Han S-K, Wagner D. Role of chromatin in water stress responses in plants. J Exp Bot. 2014;65(10):2785–99.

    Article  CAS  PubMed  Google Scholar 

  70. Delk NA, Johnson KA, Chowdhury NI, Braam J. CML24, regulated in expression by diverse stimuli, encodes a potential Ca2 + sensor that functions in responses to abscisic acid, daylength, and ion stress. Plant Physiol. 2005;139(1):240–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Hua J, Grisafi P, Cheng S-H, Fink GR. Plant growth homeostasis is controlled by the Arabidopsis BON1 and BAP1 genes. Genes Dev. 2001;15(17):2263–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Melgar AE, Zelada AM. Evolutionary analysis of angiosperm dehydrin gene family reveals three orthologues groups associated to specific protein domains. Sci Rep. 2021;11(1):23869.

    Article  CAS  PubMed  PubMed Central  ADS  Google Scholar 

  73. Bjornson M, Pimprikar P, Nürnberger T, Zipfel C. The transcriptional landscape of Arabidopsis thaliana pattern-triggered immunity. Nat Plants. 2021;7(5):579–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Pandey N, Ranjan A, Pant P, Tripathi RK, Ateek F, Pandey HP, Patre UV, Sawant SV. CAMTA 1 regulates drought responses in Arabidopsis thaliana. BMC Genomics. 2013;14(1):1–23.

    Article  Google Scholar 

  75. Kidokoro S, Yoneda K, Takasaki H, Takahashi F, Shinozaki K, Yamaguchi-Shinozaki K. Different cold-signaling pathways function in the responses to rapid and gradual decreases in temperature. Plant Cell. 2017;29(4):760–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Kollist H, Zandalinas SI, Sengupta S, Nuhkat M, Kangasjärvi J, Mittler R. Rapid responses to abiotic stress: priming the landscape for the signal transduction network. Trends Plant Sci. 2019;24(1):25–37.

    Article  CAS  PubMed  Google Scholar 

  77. Pulido P, Leister D. Novel DNAJ-related proteins in Arabidopsis thaliana. New Phytol. 2018;217(2):480–90.

    Article  CAS  PubMed  Google Scholar 

  78. Bourgine B, Guihur A. Heat shock signaling in land plants: from plasma membrane sensing to the transcription of small heat shock proteins. Front Plant Sci. 2021;12:710801.

    Article  PubMed  PubMed Central  Google Scholar 

  79. Jacob P, Hirt H, Bendahmane A. The heat-shock protein/chaperone network and multiple stress resistance. Plant Biotechnol J. 2017;15(4):405–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Hahn A, Kilian J, Mohrholz A, Ladwig F, Peschke F, Dautel R, Harter K, Berendzen KW, Wanke D. Plant core environmental stress response genes are systemically coordinated during abiotic stresses. Int J Mol Sci. 2013;14(4):7617–41.

    Article  PubMed  PubMed Central  Google Scholar 

  81. Suzuki N, Devireddy AR, Inupakutika MA, Baxter A, Miller G, Song L, Shulaev E, Azad RK, Shulaev V, Mittler R. Ultra-fast alterations in mRNA levels uncover multiple players in light stress acclimation in plants. Plant J. 2015;84(4):760–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Gilroy S, Białasek M, Suzuki N, Górecka M, Devireddy AR, Karpiński S, Mittler R. ROS, calcium, and electric signals: key mediators of rapid systemic signaling in plants. Plant Physiol. 2016;171(3):1606–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. Yang L, Zhao Y, Zhang G, Shang L, Wang Q, Hong S, Ma Q, Gu C. Identification of CAMTA Gene Family in Heimia myrtifolia and expression analysis under Drought stress. Plants. 2022;11(22):3031.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  84. Liu C, Tang D. Comprehensive identification and expression analysis of CAMTA gene family in Phyllostachys edulis under abiotic stress. PeerJ. 2023;11:e15358.

    Article  PubMed  PubMed Central  Google Scholar 

  85. Zhou Q, Zhao M, Xing F, Mao G, Wang Y, Dai Y, Niu M, Yuan H. Identification and expression analysis of CAMTA genes in tea plant reveal their complex regulatory role in stress responses. Front Plant Sci. 2022;13:910768.

    Article  PubMed  PubMed Central  Google Scholar 

  86. Zhu X, Wang B, Wei X, Du X. Characterization of the CqCAMTA gene family reveals the role of CqCAMTA03 in drought tolerance. BMC Plant Biol. 2022;22(1):428.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  87. Wang D, Wu X, Gao S, Zhang S, Wang W, Fang Z, Liu S, Wang X, Zhao C, Tang Y. Systematic analysis and identification of drought-responsive genes of the CAMTA gene family in wheat (Triticum aestivum L). Int J Mol Sci. 2022;23(9):4542.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  88. Baek D, Cho HM, Cha YJ, Jin BJ, Lee SH, Park MS, Chun HJ, Kim MC. Soybean calmodulin-binding transcription activators, GmCAMTA2 and GmCAMTA8, coordinate the circadian regulation of developmental processes and Drought stress responses. Int J Mol Sci. 2023;24(14):11477.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  89. Estavillo GM, Crisp PA, Pornsiriwong W, Wirtz M, Collinge D, Carrie C, Giraud E, Whelan J, David P, Javot H. Evidence for a SAL1-PAP chloroplast retrograde pathway that functions in drought and high light signaling in Arabidopsis. Plant Cell. 2011;23(11):3992–4012.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  90. Cejudo FJ, González M-C, Pérez-Ruiz JM. Redox regulation of chloroplast metabolism. Plant Physiol. 2021;186(1):9–21.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors also thank the Salk Institute Genomic Analysis Laboratory for providing the T-DNA insertion mutants SALK_046204 and the Arabidopsis Biological Resource Center for distributing the seeds. We are grateful for the technical contributions of Laura Aguilar-Henonin in the initial stages of this work.

Funding

This work was supported by the Ciencia de Frontera grant 2014 from the Consejo Nacional de Humanidades Ciencia y Tecnología, Mexico, to PG. D. S-C. and M.M-G are recipients of a postdoctoral scholarship from the program “Estancias Postdoctorales por México” by Conahcyt.

Author information

Authors and Affiliations

Authors

Contributions

P.G., E.J-M., D. S-C., and M.M-G. conceived and designed the project and conducted experiments. E.J-M., D. S-C. and M.M-G. performed the experiments and analyzed the data. P.G. and M.M-G. wrote the manuscript. P.G. obtained the funding. All authors have read and approved the manuscript.

Corresponding author

Correspondence to Plinio Guzmán.

Ethics declarations

Ethics approval and consent to participate

The seeds employed in this study were obtained from the ARABIDOPSIS BIOLOGICAL RESOURCE CENTER (ABRC), obviating the requirement for field permissions to gather samples. The authors declare that the experimental research conducted on the plants detailed in this paper conforms to institutional, national, and international guidelines.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

12870_2024_4872_MOESM1_ESM.xlsx

Supplementary Material 1: Table 1. Differentially expressed genes between the atl80 mutant and the WT transcriptomes at different time points.

12870_2024_4872_MOESM2_ESM.xlsx

Supplementary Material 2: Table 2. Gene expression rapidly induced in response to water deprivation stress compiled in twelve clusters.

12870_2024_4872_MOESM3_ESM.xlsx

Supplementary Material 3: Table 3. Cross-reference of transcriptomic data related to retrograde signaling pathways with the ATL80 transcriptome.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Méndez-Gómez, M., Sierra-Cacho, D., Jiménez-Morales, E. et al. Modulation of early gene expression responses to water deprivation stress by the E3 ubiquitin ligase ATL80: implications for retrograde signaling interplay. BMC Plant Biol 24, 180 (2024). https://doi.org/10.1186/s12870-024-04872-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-024-04872-5

Keywords