Ruth Ndathe
Renee Dale
Naohiro Kato- 1Department of Biological Sciences, Louisiana State University, Baton Rouge, LA, United States
- 2Donald Danforth Plant Science Center, St. Louis, MO, United States
The abscisic acid (ABA) signaling pathway is the key defense mechanism against drought stress in plants. In the pathway, signal transduction among four core proteins, pyrabactin resistance (PYR), protein phosphatase 2C (PP2C), sucrose-non-fermenting-1-related protein kinase 2 (SnRK2), and ABRE binding factor (ABF) leads to altered gene expression kinetics that is driven by an ABA-responsive element (ABRE). A most recent and comprehensive study provided data suggesting that ABA alters the expression kinetics in over 6,500 genes through the ABF-ABRE associations in Arabidopsis. Of these genes, termed ABA gene regulatory network (GRN), over 50% contain a single ABRE within 4 kb of the gene body, despite previous findings suggesting that a single copy of ABRE is not sufficient to drive the gene expression. To understand the expression system of the ABA GRN by the single ABRE, a dynamic model of the gene expression for the desiccation 29A (RD29A) gene was constructed with ordinary differential equations. Parameter values of molecular-molecular interactions and enzymatic reactions in the model were implemented from the data obtained by previously conducted in vitro experiments. On the other hand, parameter values of gene expression and translation were determined by comparing the kinetics of gene expression in the model to the expression kinetics of RD29A in real plants. The optimized model recapitulated the trend of gene expression kinetics of RD29A in ABA dose–response that were previously investigated. Further analysis of the model suggested that a single ABRE controls the time scale and dynamic range of the ABA-dependent gene expression through the PP2C feedback regulation even though an additional cis-element is required to drive the expression. The model construed in this study underpins the importance of a single ABRE in the ABA GRN.
Introduction
Plants possess defense mechanisms against abiotic stresses (Basu et al., 2016; Kumar et al., 2018; Takahashi et al., 2020). One of the primary mechanisms is the abscisic acid (ABA) signaling pathway. ABA is a phytohormone that is produced under abiotic stresses such as drought conditions (Zeevaart and Creelman, 1988; Sauter et al., 2001; Ikegami et al., 2008). The ABA signaling pathway has been well-characterized, leading to downstream ABA responses such as stomatal closure and gene expression that help the plant acquire drought stress resistance (Steuer et al., 1988; Fujii et al., 2009; Umezawa et al., 2009). The most upstream of the core components in the ABA signaling pathway are ABA receptors named pyrabactin resistance/pyr1-like/ regulatory components of ABA receptors (PYR/PYL/RCAR) that bind ABA and, in turn, interact with different protein phosphatase 2Cs (PP2Cs), namely aba insensitive1/2 (ABI1/ABI2), hypersensitive to aba1/2 (HAB1/HAB2), aba-hypersensitive germination 3 (AHG3/PP2CA), and highly aba induced 1/2/3 (HA1/2/3). The PP2Cs inhibit SNF1-related protein kinase 2 s (SnRK2s), including SnRK2.2, SnRK2.3, and SnRK2.6, when they are not interacting with PYR (Rodriguez et al., 1998; Gosti et al., 1999; Merlot et al., 2001; Saez et al., 2004; Ma et al., 2009; Melcher et al., 2009; Nishimura et al., 2009; Park et al., 2009; Santiago et al., 2009; Yin et al., 2009; Soon et al., 2012). For gene expression, activated SnRK2s phosphorylate ABA-responsive elements (ABRE) binding factors 1/2/3/4 (ABF1/2/3/4). These phosphorylated transcription factors bind ABRE, a regulatory region of ABA-induced genes (Choi et al., 2000; Uno et al., 2000; Yoshida et al., 2015). For stomatal closure, the activated SnRK2, namely SnRK2.6 kinase, phosphorylate the slow-anion channels (SLAC1), leading to an anion and K+ efflux and eventual solute loss from the guard cells (Schroeder et al., 1984; Geiger et al., 2009; Lee et al., 2009; Albert et al., 2017).
A relationship between ABF-ABRE associations and kinetics of ABA-dependent gene expression has been studied extensively to understand the downstream response that changes growth and physiology with a function of time in plants. The current understanding is that a single copy of ABRE is not sufficient for the ABA-dependent expression (Shen et al., 1996; Hobo et al., 1999). Genes containing several ABREs in the promoter are mainly regulated by the ABA signaling pathway (Yamaguchi-Shinozaki and Shinozaki, 1994; Fujii et al., 2009; Ma et al., 2009). Namely, a pair of ABREs has been shown to be overrepresented in the promoter region of ABA-inducible genes in Arabidopsis (Zhang et al., 2005; Gómez-Porras et al., 2007). However, the data obtained in the most recent study with DAP-Seq (DNA affinity purification and sequencing) and RNA-seq (RNA sequencing) suggested that over 50% of the genes regulated by the ABF-ABRE associations contain a single ABRE within 4 kb of the gene body (Sun et al., 2022). This new data set suggests that genes containing a single ABRE are the primary target of the ABA signaling pathway. Although the requirement of the second cis-element for the expression of these genes has been shown (Narusaka et al., 2003), how the ABRE and second cis-elements co-operatively change the kinetics of ABA-dependent gene expression is little understood.
Dynamic modeling is a powerful tool that integrates extensive experimental data of pathway components, improving our understanding of the signaling pathway dynamics and making novel hypotheses and predictions (Poolman et al., 2004; Aldridge et al., 2006; Janes and Yaffe, 2006; Thakar et al., 2007). The network connectivity of the core components of the ABA signaling pathway has been revealed. Furthermore, in vitro parameters for many of the interactions of the core components have been experimentally determined, making this a good candidate for modeling.
This study aims to build a dynamic model of ABA-responsive gene expression with a single ABRE element, namely the RD29A gene that has been used as the marker for the ABA signaling pathway in Arabidopsis. The RD29A gene requires DRE (dehydration responsive element) to which DREB2A (DRE binding protein 2A) binds (Yamaguchi-Shinozaki and Shinozaki, 1994; Liu et al., 1998) in addition to ABRE for its expression (Yamaguchi-Shinozaki et al., 1995; Narusaka et al., 2003). Because the network structure of the RD29A gene was well established and the parameter values were already obtained experimentally, testing the effects of network structures on the model performance was not focused on in this study. Approximate curve fitting of the model output to actual plant data was conducted by optimizing parameter values of transcription and translation, which were not determined previously. This report describes how we built, optimized, and validated the model. The resulting model led us to form a new hypothesis that a single ABRE is not sufficient to derive the gene expression yet controls the time scale and the dynamic range of the ABA-dependent gene expression kinetics.
Description
Identification of ABA-induced genes containing a single ABRE in Arabidopsis
Supplementary Data 2 in the publication by Sun et al. (2022) was used to extract the data on genes that are differentially expressed by ABA and the number of ABRE in the genes in Arabidopsis. In the publication, gene expression was determined with RNA-seq. DEG (differentially expressed gene) was defined by DESeq2 as a gene whose value of p for the significance of differential gene expression is less than 0.05 in 10 μM ABA-treated samples compared to controls. Seven days old seedlings (shoot and root tissues) and two-time points (3 h and 24 h after the treatment) were compared. The number of ABREs in a gene was determined with DAP-seq. The ABRE element was defined as an ACGT sequence onto which ABF1, 2, 3, or 4 binds within the 5’ 2Kb upstream and 3’ 2Kb downstream in a gene’s protein-coding sequence (gene body).
Construction of the dynamic model of the RD29A gene expression
A previous study defined a minimal set of core components and the signal transduction that led to ABA-induced gene expression (Yoshida et al., 2014; Singh and Laxmi, 2015; Wang et al., 2019). The components are ABA, PYR, PP2C, SnRK2, ABF, DREB2A, ABRE, and DRE (Figure 1).
Figure 1. Schematic drawing of the ABA signaling pathway for the RD29A expression. ABA signaling core components and their interactions are shown with and without ABA. In the presence of ABA, downstream regulators ABF and DREB2A bind ABRE and DRE, the cis-elements on the RD29A gene, leading to its transient expression. The drawing was modified from the figure in (Wang et al., 2019). P with a circle represents the phosphorylation of a protein. A different number of arrowheads indicates a different level of gene expression. X indicates no gene expression.
We included other components necessary to connect each core component functionally to make a dynamic model. Previous studies have determined that the PP2C phosphatases dephosphorylate phosphorylated ABFs (Antoni et al., 2012; Lynch et al., 2012). The components were hence included in the model. In addition, another study identified that SnRK2s are enzymatically phosphorylated by MAP3Ks, RAF-like kinases, although the regulation of MAP3Ks by ABA has not yet been revealed (Katsuta et al., 2020; Lin et al., 2020; Takahashi et al., 2020). It is also known that BIN2 (BRASSINOSTEROID INSENSITIVE 2 kinase) phosphorylates SnRK2s, although its regulation by ABA is unknown (Cai et al., 2014). To this end, we added components that phosphorylate SnRK2s independently from the ABA regulation (presented as MAP3K in the model). DREB2A is subjected to 26 s proteasome proteolysis in a normal condition, but stress conditions block the proteolysis of DREB2A through a yet unknown mechanism (Qin et al., 2008). We added two interaction reactions to implement the finding in the model. One is between DREB2A and 26 s proteasome, the complex of which leads to degradation of DREB2A. The other is between 26 s proteasome and ABA, a complex which deactivates the degradation. We also included the feedback regulation in which the expression of PP2C, ABF, and DREB2A genes is upregulated by the ABRE promoter activity (Kim et al., 2011; Wang et al., 2019). In the dynamic model, a set of 33 variables and 63 parameters representing biochemical reactions of each component were constructed based on the law of mass action. In the model, a single protein in a homologous protein family that redundantly function in the cells regulates the system. Values of parameters in the equations were obtained from the literature (Table 1). The equations, initial conditions (concentrations), and parameter values were compiled and analyzed numerically with default settings using MATLAB SimBiology (MathWorks) (Supplementary File S1).
Table 1. Curated values from literature and the values chosen as parameters for the model.
In the model, we assumed:
• ABA signal transduction occurs through the system composed of ABA, PYR, PP2C, SnRK2, ABF, DREB2A, ABRE, and DRE, in which RD29A is expressed.
• Enzymatic reactions follow Michaelis–Menten kinetics.
• All molecules freely diffuse in the cell. This is not true for all components, but the assumption is necessary for modeling.
• The cell volume is 50 μm3.
• The Michaelis constant is = , where is the dissociation rate constant, is the catalytic rate constant, and is the association rate constant.
• A molecule associates with another molecule at a rate constant of = 1,000 μM−1 s−1 (Milo and Phillips, 2015).
• Proteins are generated by reactions of gene expression and protein translation, then subject to degradation.
• The protein concentration in a cell remains at 0.1 μM at a steady state without ABA activation and feedback regulation.
• A gene (mRNA) is expressed from a pair of gene loci with a constitutively active promoter and then subjected to degradation.
• Genes with feedback regulation (ABF, PP2C, DREB2A) have ABREs and a constitutively active promoter. Hence, upon binding the activated ABF, expression is increased compared to the level of constitutive expression.
• The binding of ABF and DREB2A to ABRE and DRE, respectively, is required for RD29A expression.
• Initial values of variables (components) are 0 (zero) except for the gene in the genome.
The model was first to run for 300 equivalent hours in numerical analysis with the variable ABA (representing intracellular ABA) set at 0 μM. This allows the system to reach a quasi-steady state. After the 300 equivalent hours, the variable ABA was set to 100 μM. Changes in all variables in the model from the quasi-steady state were then monitored for another 300 equivalent hours. This report presents the time when the variable ABA is changed to time zero.
Optimization of parameters, validation of the model, and analyzing identifiability of model parameters
We approximately curve fit model output to experimental data to optimize selected model parameters. We focused on changes in the variable RD29A, representing accumulated mRNA expressed. Three parameters, 1. transcription of feedback ABRE containing genes, 2. transcription of RD29A, 3. translation of PP2C, ABF, and DREB2A affected by feedback, were manually changed to obtain a qualitatively good fit to experimental data. The remaining model parameters were unchanged (fixed). To validate the model, we quantitatively evaluated changes of the variable RD29A. Fold changes calculated by the model were compared to previously published data or newly obtained in this study. Our new experiments used the transgenic Arabidopsis thaliana plant (Supplementary Methods). The transgenic plant carries the RD29A::LUC gene expression cassette that drives the expression of luciferase (LUC) from the RD29A promoter in the genome containing ABRE and DRE cis-elements (Zhan et al., 2012). Because the half-life of luciferase is shorter than that of the mRNA, the activity of luciferase that is detected as the emission of light (luminescence) can be used to track the accumulation of the mRNA in near real-time. To analyze the identifiability of the variable RD29A dynamics, we conducted a sensitivity analysis using Calculate Sensitivity in Model Analyzer in SimBiology with default settings.
Results
Over 50% of the genes in the ABA gene regulatory network, including RD29A, contain a single ABRE in the gene body
A most recent and comprehensive study conducted by Sun et al. identified genes that are differentially expressed by ABA and the number of ABRE in the genes in Arabidopsis (Sun et al., 2022). A set of genes whose expression is significantly altered by 10 μM ABA through binding of activated ABF on ABRE in the root or shoot tissue in 7 days old seedlings were defined as an ABA GRN (gene regulatory network) (Sun et al., 2022). We found that 53% of the GRN carry a single ABRE within the gene body (Figure 2A and Supplementary Table S1). The range of a fold-change in the gene expression in the genes carrying a single ABRE in their gene body is as wide as those carrying multiple ABREs in their gene body (Figure 2B and Supplementary Table S1). This suggests that the genes carrying a single ABRE are the major contributors to the ABA GRN. RD29A, one of the most studied genes as the ABA-induced gene marker, also carries a single ABRE in the gene body (Sun et al., 2022).
Figure 2. Genes regulated through ABRE-ABF binding in Arabidopsis. (A) Genes in the ABA GRN (Sun et al., 2022) are categorized by the number of ABRE in their gene body. The number of genes belonging to each category and percentage in the ABA GRN is shown with a separate marker “|.” The bars in the figure indicate the percentages. (B) Box plots of genes in the ABA GRN. X-axis variables are the number of ABRE in a gene. Y-axis variables are changes in gene expression (lg2) by ABA in root tissue (3 h after the ABA treatment). Black dots indicate points of individual genes. All genes in the ABA GRN are listed, together with gene expression in root and shoot tissues at different time points, in Supplementary Table S1.
Parameter values of the RD29A gene expression were obtained by literature curation
To understand the alteration of gene expression kinetics by the single ABRE, a dynamic model of the RD29A gene expression was constructed with ordinary differential equations. We curated previously published data to define parameters in the model of the ABA signaling pathway that activates the ABF, resulting in the activation of the gene promoter containing the ABRE and DRE. The summary of our curation is shown below (Table 1).
While parameter values for protein–protein interactions and enzymatic reactions were characterized in vitro studies using recombinant proteins, no studies related to parameter values of DNA-protein binding, gene expression, protein translation, and degradation were found for the ABA signaling pathway. To this end, we implemented parameter values from studies using non-plant eukaryotic organisms. These parameters had a wide range to select from 1. equilibrium dissociation constant between the transcription factors and respective cis-elements, (from 2 nM to 2 μM) (Geertz et al., 2012), 2. translation rate of protein from mRNA expressed (less than 10,000 h−1) (Hausser et al., 2019), 3. transcription rates (slower than the translation rate) (Hausser et al., 2019). We selected the translation and transcription rates for genes at 4.5 h−1 and 1 h−1, respectively, and 2 nM for transcription factor-cis-element binding. This is because the average rate of gene transcription in multicellular eukaryotes is 1 h−1 (Hausser et al., 2019), while the average concentration of proteins involved in signal transduction is 0.1 μM (Milo and Phillips, 2015). Setting translation rate at 4.5 h−1 and transcription rate at 1 h−1 makes the concentration of a protein at a quasi-steady state to 0.1 μM without ABA and feedback regulation in our model. The affinity of transcription factor-cis-element binding was set at 2 nM to fit the kinetics of the variable RD29A with actual gene expression (Figure 3). Protein degradation was set at 0.05 h−1 (Hausser et al., 2019). The equilibrium dissociation constant between SnRK2 (non-phosphorylated SnRK2) and PP2C was set at 100 pM, representing complete inhibition of SnRK2 kinase activity by PP2C at an equal molar concentration (Soon et al., 2012).
Figure 3. The dynamic model agrees with ABA-induced gene expression in real plants after optimization. (A) Kinetics of luciferase activity in the RD29A::LUC plant after exposure to 200 μM ABA (+ABA) or DMSO for control (-ABA). The graph shows the mean of three independent experiments. Error bars represent standard error from the mean. (B) Kinetics of RD29A gene accumulation in the previously published data with 10 μM ABA in Arabidopsis (Song et al., 2016). (C) Model output with showing transient expression after optimization of parameters kf26, kf27, and kf28.
The rates of transcription and translation in RD29A and the feedbacked genes, ABFs, PP2Cs, and DREB2A, were optimized in the model to capture observed dynamics in experimental data
To understand the connectivity of the components, we compared the kinetics of gene expression in the model and experimental data in actual plants. We compared the simulation data of the variable RD29A, to two independent data sets that were experimentally obtained using actual plants. One set of data was obtained by our new experiments using transgenic Arabidopsis thaliana carrying the RD29A::LUC (Zhan et al., 2012). The other set was obtained from previously published data that show a change in the RD29A gene expressed in Arabidopsis thaliana (Song et al., 2016). Kinetics of the gene expression in the plants and the variable RD29A were compared within the first 24 h (Figure 3).
Experimental data from the transgenic RD29A::LUC plants showed transient activation of its promoter with an initial increase and then a decrease after 5 h (Figure 3A). Similar transient expressions of the RD29A gene were observed in non-transgenic Arabidopsis plants, (Song et al., 2016) (Figure 3B). Despite the RD29A transient expression pattern being consistently observed in both transgenic and wild-type plants, the mechanism that explains the transient expression was unknown. Therefore, we investigated the mechanism using the dynamic model we constructed. When we simulated the kinetics of the variable RD29A, the kinetics were logarithmic upon adding ABA. We therefore optimized the parameters so that the model’s kinetics in the gene expression qualitatively agree with that in actual plants (Figure 3C). We altered three parameters, the feedback transcription rate constant of the ABRE promoter (parameter kf26), the transcription rate constant of RD29A with both DRE and ABRE elements (parameter kf27), and the feedback translation rate constants of ABF, PP2C, and DREB2A (parameter kf28). These three parameters had not been determined previously, and studies in other eukaryotic cells indicate wide ranges of reasonable values (Table 1). Hence, we manually altered the values within the ranges of a previous biological study (Hausser et al., 2019) so that the kinetics of the variable RD29A resembles the actual plant data. Changes on these parameters most affected the aspect of transient increase of the variable RD29A. The values 5 < kf26 < 10 h−1, 5 < kf27 < 10 h−1, and 150 < kf28 < 250 h−1 recapitulated the trend of gene expression (Figure 3).
Approximation of the model was validated by determining model responses to different doses of ABA or a set of gene null-mutations
To validate the model, we first compared the ABA-dose-dependent response in actual plants to the dynamics of the variable RD29A (Figure 4). In the model, changes of the variable RD29A increased in an ABA-dose-dependent manner from 0 to 200 μM (Figure 4A). With the RD29A::LUC transgenic plants, changes in luminescence increased in an ABA-dose-dependent way in the range from 0 to 200 μM (Figure 4B). The results would depend very much on time after treatment with ABA. Hence, we sampled the actual plants 5 h after exposure to ABA when the largest dose–response would be expected. The comparison of the outcomes suggested that the model is approximated to actual plants concerning ABA sensitivity. However, the response in the model seems to have narrower sensitivity against the ABA concentration (i.e., from 0 to 50 μM) compared to that in the actual plants (i.e., from 0 to 200 μM) (Figure 4B) (Gampala et al., 2001; Lee et al., 2016).
Figure 4. RD29A expression increases with a function of ABA concentration in the model as it is observed in actual plants. (A) Model output of the variable RD29A with different values of the variable ABA. (B) The relative luminescence unit in 25-day-old RD29A::LUC plants was determined at 5 h after spraying different concentrations of ABA. The bars represent the mean relative luminescence of three replicates, with error bars representing standard error from the mean (15 seedlings).
We also validated changes of the variable RD29A in gene-knockout simulations. Namely, we simulated the expression RD29A gene in gene null-mutations of pyr, pp2c, snrk2, and abf, which were previously studied (Fujita et al., 2009; Rubio et al., 2009; Nishimura et al., 2010; Yoshida et al., 2015). In these studies, different stages of seedlings were treated differently with ABA. Hence, we evaluate the expression of RD29A qualitatively (increased or decreased) but not quantitatively (i.e., comparing changes in the expression kinetics). We simulated knockout mutations by setting the translation rate constant (kf30) to zero for the variable PYR, PP2C, SnRK2, and ABF. In addition, we also set the translation rates of the feedback regulations kf28 to zero for ABF and PP2C, respectively. The mimicked null-mutant in pyr, snrk2, and abf, all showed reduced levels of the variable RD29A, while the mimicked null-mutant in pp2c showed elevated levels (Table 2).
Table 2. Mutant simulations in the model show qualitative similarity to actual mutant plants, concerning the RD29A expression.
Experimental data in actual plants shows that pyr null-mutants are impaired in ABA-induced gene expression (Park et al., 2009; Nishimura et al., 2010; Gonzalez-Guzman et al., 2012). Similarly, experimental data on snrk2.2/ snrk2.3/ snrk2.6 triple knockout mutants showed that the expression of ABA-induced genes was impaired (Fujii and Zhu, 2009; Fujita et al., 2009; Thalmann et al., 2016). Triple areb/abf mutants were found to have reduced ABA-induced gene expression (Yoshida et al., 2015; Thalmann et al., 2016). On the other hand, null mutants of pp2cs in actual plants show a higher and constitutive ABA response (Rubio et al., 2009; Antoni et al., 2012). Based on the two validations described above, we concluded that the model constructed, and parameters implemented in the model are approximated to actual plants.
The dynamic model predicts that a single ABRE makes the RD29A expression transient through the PP2C feedback loop during ABA exposure
To understand which parameters are sensitive for the RD29A gene expression, we conducted a sensitivity analysis on fourteen parameters that determine protein binding, enzymatic activity, and protein degradation of the key compartments were selected (Figure 5). The analysis showed that the parameter related to the ABF-P binding to ABRE for the feedbacked regulation genes (ABF, PP2C, ABF, DREB2A) (kr10) and the parameter related to the ABF-P biding to ABRE for the RD29A expression (kr14) were the most sensitive. Parameters associated with DREB2A binding of the DRE were not significantly sensitive as those related to ABF-P binding to ABRE (Figure 5). However, the binding affinity of DREB2A-DRE certainly affected the expression levels of the RD29A mRNA (Supplementary Figure S1). Interaction of DREB2A and 26 s proteasome (kr16) was most sensitive after those related to ABF-P binding to ABRE.
Figure 5. Sensitivity analysis of different parameters on RD29A gene expression. A sensitivity analysis conducted against the variable RD29A gene determined that the parameters for ABF-P-ABRE binding are most sensitive to the kinetics of RD29A gene expression. The y-axis shows the level of influence of the respective parameter on the output or RD29A gene expression. The higher the value, the higher the effect of the parameter on the gene expression.
Because the parameter representing feedback loops of the ABA signaling pathway was identified as the most sensitive for the RD29A expression, we individually examined the effect of the ABF, DREB2A, and PP2C feedback loops. On removal of the feedback loop on ABF or DREB2A, the expression dynamics did not change. Still removing the PP2C feedback loop created logarithmic gene expression (Figure 6). Also, the concentration of RD29A mRNA increased about 100-fold higher without the PP2C feedback loop. These indicate that the feedback loop of the PP2C gene expression makes the RD29A gene expression transient and reduces the dynamic range.
Figure 6. Effect of the different feedback loops on RD29A gene expression dynamics. (A) Without the ABF feedback loop. (B) Without the DREB2A feedback loop. (C) Without the PP2C feedback loop.
Discussion
The ABRE is the main element involved in ABA-mediated gene expression. However, previous studies have shown that a single ABRE copy is not sufficient for gene expression, and other elements like DRE, CE1, or CE3 are required as coupling elements. (Shen et al., 1996; Abe et al., 1997; Hobo et al., 1999; Uno et al., 2000; Zhu, 2002). Yet, the most comprehensive genome-wide study indicated that genes carrying a single ABRE are the major genes in the ABA GRN (Sun et al., 2022) (Figure 2).
The role of the single ABRE on the gene expression dynamics in the ABA GRN was unknown (Narusaka et al., 2003; Sun et al., 2022). Here, we presented a model of the ABA signaling pathway describing the expression dynamics of RD29A under the control of the ABRE and DRE (Figure 3). The model was built with fixed parameter values of protein–protein interactions and enzymatic kinetics obtained by in vitro experiments from the literature. The model was further validated by comparing RD29A gene expression from the model’s output to actual plant data (Figure 4 and Table 2). The model suggests that a single ABRE controls the transient nature of the gene expression through the feedback loop of the PP2C gene expression (Figure 6). In addition, the model suggests that PP2C largely suppressed the expression level of RD29A (Figure 6). On the other hand, the feedback loop of DREB2A that binds to DRE does not affect the kinetics of the RD29A expression (Figure 6), despite a fact that DRE also regulates the level of the ABA-dependent RD29A gene expression (Supplementary Figure S1). Based on the model analyses, we hypothesize that a single ABRE regulates the time scale and dynamic range of the expression in the ABA GRN although the single ABRE itself is not sufficient to drive the expression.
In this study, a direct quantitative comparison between the model output and actual data was not conducted because we modeled the system as a single cell, whereas the actual data obtained from the multicellular system. Furthermore, each homologous protein that redundantly functions in the signaling pathway may have a different parameter. For instance, the homologous protein family of 14 PYR and 9 PP2Cs were shown to have different affinities in the ABA-responsive gene expression (Tischer et al., 2017). Accordingly, our model cannot simulate the plant response exactly. For instance, when an ABA-concentration-dependent response of the ABRE promoter was determined, the response range was narrower in the model than in existing plants (Figure 4). Optimization of parameter values fixed in this study may be required to improve model performance.
Nevertheless, our model successfully builds off existing work to represent the relationship between the ABA signaling pathway and the gene expression regulation by a single ABRE. The model construed in this study underpins the importance of a single ABRE in the ABA GRN. Revealing the interaction between the single ABRE with other cis-elements in the regulatory region of each gene would be the next frontier for understanding the ABA GRN.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/genbank/, AY142623, BT008860, BT015409, BT025246, BT002082, AY081538, BT026443, AY081467, and AY091298.
Author contributions
NK: conceptualization, methodology, and funding acquisition. RN and RD: validation. RN: experiments. RN and NK: formal analysis and writing—original draft preparation. RN, RD, and NK: writing—review and editing. All authors contributed to the article and approved the submitted version.
Funding
This study is partly supported by Economic Development Assistantships from Louisiana State.
Acknowledgments
NK thanks Dr. Dong-Ha Oh for his assistance.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.928718/full#supplementary-material
References
Abe, H., Yamaguchi-Shinozaki, K., Urao, T., Iwasaki, T., Hosokawa, D., and Shinozaki, K. (1997). Role of arabidopsis MYC and MYB homologs in drought-and abscisic acid-regulated gene expression. Plant Cell 9, 1859–1868. doi: 10.1105/tpc.9.10.1859
Albert, R., Acharya, B. R., Jeon, B. W., Zañudo, J. G. T., Zhu, M., Osman, K., et al. (2017). A new discrete dynamic model of ABA-induced stomatal closure predicts key feedback loops. PLoS Biol. 15:e2003451. doi: 10.1371/journal.pbio.2003451
Aldridge, B. B., Burke, J. M., Lauffenburger, D. A., and Sorger, P. K. (2006). Physicochemical modelling of cell signalling pathways. Nat. Cell Biol. 8, 1195–1203. doi: 10.1038/ncb1497
Antoni, R., Gonzalez-Guzman, M., Rodriguez, L., Rodrigues, A., Pizzio, G. A., and Rodriguez, P. L. (2012). Selective inhibition of clade A phosphatases type 2C by PYR/PYL/RCAR Abscisic acid receptors 1[C][W]. Plant Physiol. 158, 970–980. doi: 10.1104/pp.111.188623
Bar-Even, A., Noor, E., Savir, Y., Liebermeister, W., Davidi, D., Tawfik, D. S., et al. (2011). The moderately efficient enzyme: evolutionary and physicochemical trends shaping enzyme parameters. Biochemistry 50, 4402–4410. doi: 10.1021/bi2002289
Basu, S., Ramegowda, V., Kumar, A., and Pereira, A. (2016). Plant adaptation to drought stress. F1000Research 5, F1000 faculty Rev-1554 5. doi: 10.12688/f1000research.7678.1
Cai, Z., Liu, J., Wang, H., Yang, C., Chen, Y., Li, Y., et al. (2014). GSK3-like kinases positively modulate abscisic acid signaling through phosphorylating subgroup III SnRK2s in Arabidopsis. Proc. Natl. Acad. Sci. U. S. A. 111, 9651–9656. doi: 10.1073/pnas.1316717111
Choi, H., Hong, J., Ha, J., Kang, J., and Kim, S. Y. (2000). ABFs, a family of ABA-responsive element binding factors *. J. Biol. Chem. 275, 1723–1730. doi: 10.1074/jbc.275.3.1723
Dupeux, F., Santiago, J., Betz, K., Twycross, J., Park, S.-Y., Rodriguez, L., et al. (2011). A thermodynamic switch modulates abscisic acid receptor sensitivity. EMBO J. 30, 4171–4184. doi: 10.1038/emboj.2011.294
Fujii, H., Chinnusamy, V., Rodrigues, A., Rubio, S., Antoni, R., Park, S.-Y., et al. (2009). In vitro reconstitution of an ABA signaling pathway. Nature 462, 660–664. doi: 10.1038/nature08599
Fujii, H., and Zhu, J.-K. (2009). Arabidopsis mutant deficient in 3 abscisic acid-activated protein kinases reveals critical roles in growth, reproduction, and stress. Proc. Natl. Acad. Sci. 106, 8380–8385. doi: 10.1073/pnas.0903144106
Fujita, Y., Nakashima, K., Yoshida, T., Katagiri, T., Kidokoro, S., Kanamori, N., et al. (2009). Three SnRK2 protein kinases are the Main positive regulators of Abscisic acid signaling in response to water stress in Arabidopsis. Plant Cell Physiol. 50, 2123–2132. doi: 10.1093/pcp/pcp147
Gampala, S. S. L., Hagenbeek, D., and Rock, C. D. (2001). Functional interactions of lanthanum and phospholipase D with the Abscisic acid signaling effectors VP1 and ABI1-1 in Rice protoplasts *. J. Biol. Chem. 276, 9855–9860. doi: 10.1074/jbc.M009168200
Geertz, M., Shore, D., and Maerkl, S. J. (2012). Massively parallel measurements of molecular interaction kinetics on a microfluidic platform. Proc. Natl. Acad. Sci. 109, 16540–16545. doi: 10.1073/pnas.1206011109
Geiger, D., Scherzer, S., Mumm, P., Stange, A., Marten, I., Bauer, H., et al. (2009). Activity of guard cell anion channel SLAC1 is controlled by drought-stress signaling kinase-phosphatase pair. Proc. Natl. Acad. Sci. 106, 21425–21430. doi: 10.1073/pnas.0912021106
Ghose, R. (2019). Nature of the pre-chemistry ensemble in mitogen-activated protein Kinases. J. Mol. Biol. 431, 145–157. doi: 10.1016/j.jmb.2018.12.007
Gómez-Porras, J. L., Riaño-Pachón, D. M., Dreyer, I., Mayer, J. E., and Mueller-Roeber, B. (2007). Genome-wide analysis of ABA-responsive elements ABRE and CE3 reveals divergent patterns in Arabidopsis and rice. BMC Genomics 8:260. doi: 10.1186/1471-2164-8-260
Gonzalez-Guzman, M., Pizzio, G. A., Antoni, R., Vera-Sirera, F., Merilo, E., Bassel, G. W., et al. (2012). Arabidopsis PYR/PYL/RCAR receptors play a major role in quantitative regulation of Stomatal aperture and transcriptional response to Abscisic acid. Plant Cell 24, 2483–2496. doi: 10.1105/tpc.112.098574
Gosti, F., Beaudoin, N., Serizet, C., Webb, A. A. R., Vartanian, N., and Giraudat, J. (1999). ABI1 protein phosphatase 2C is a negative regulator of Abscisic acid signaling. Plant Cell 11, 1897–1909. doi: 10.1105/tpc.11.10.1897
Hausser, J., Mayo, A., Keren, L., and Alon, U. (2019). Central dogma rates and the trade-off between precision and economy in gene expression. Nat. Commun. 10:68. doi: 10.1038/s41467-018-07391-8
Hobo, T., Asada, M., Kowyama, Y., and Hattori, T. (1999). ACGT-containing abscisic acid response element (ABRE) and coupling element 3 (CE3) are functionally equivalent. Plant J. 19, 679–689. doi: 10.1046/j.1365-313x.1999.00565.x
Ikegami, K., Okamoto, M., Seo, M., and Koshiba, T. (2008). Activation of abscisic acid biosynthesis in the leaves of Arabidopsis thaliana in response to water deficit. J. Plant Res. 122, 235–243. doi: 10.1007/s10265-008-0201-9
Janes, K. A., and Yaffe, M. B. (2006). Data-driven modelling of signal-transduction networks. Nat. Rev. Mol. Cell Biol. 7, 820–828. doi: 10.1038/nrm2041
Katsuta, S., Masuda, G., Bak, H., Shinozawa, A., Kamiyama, Y., Umezawa, T., et al. (2020). Arabidopsis Raf-like kinases act as positive regulators of subclass III SnRK2 in osmostress signaling. Plant J. 103, 634–644. doi: 10.1111/tpj.14756
Kim, J.-S., Mizoi, J., Yoshida, T., Fujita, Y., Nakajima, J., Ohori, T., et al. (2011). An ABRE promoter sequence is involved in osmotic stress-responsive expression of the DREB2A gene, which encodes a transcription factor regulating drought-inducible genes in Arabidopsis. Plant Cell Physiol. 52, 2136–2146. doi: 10.1093/pcp/pcr143
Kumar, S., Sachdeva, S., Bhat, K. V., and Vats, S. (2018). “Plant responses to drought stress: physiological, biochemical and molecular basis,” in Biotic and Abiotic Stress Tolerance in Plants. ed. S. Vats (Singapore: Springer), 1–25.
Lee, S. C., Lan, W., Buchanan, B. B., and Luan, S. (2009). A protein kinase-phosphatase pair interacts with an ion channel to regulate ABA signaling in plant guard cells. Proc. Natl. Acad. Sci. U. S. A. 106, 21419–21424. doi: 10.1073/pnas.0910601106
Lee, S. Y., Boon, N. J., Webb, A. A. R., and Tanaka, R. J. (2016). Synergistic activation of RD29A via integration of salinity stress and Abscisic acid in Arabidopsis thaliana. Plant Cell Physiol. 57, 2147–2160. doi: 10.1093/pcp/pcw132
Lin, Z., Li, Y., Zhang, Z., Liu, X., Hsu, C.-C., Du, Y., et al. (2020). A RAF-SnRK2 kinase cascade mediates early osmotic stress signaling in higher plants. Nat. Commun. 11:613. doi: 10.1038/s41467-020-14477-9
Liu, Q., Kasuga, M., Sakuma, Y., Abe, H., Miura, S., Yamaguchi-Shinozaki, K., et al. (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. doi: 10.1105/tpc.10.8.1391
Lynch, T., Erickson, B. J., and Finkelstein, R. R. (2012). Direct interactions of ABA-insensitive (ABI)-clade protein phosphatase (PP)2Cs with calcium-dependent protein kinases and ABA response element-binding bZIPs may contribute to turning off ABA response. Plant Mol. Biol. 80, 647–658. doi: 10.1007/s11103-012-9973-3
Ma, Y., Szostkiewicz, I., Korte, A., Moes, D., Yang, Y., Christmann, A., et al. (2009). Regulators of PP2C phosphatase activity function as Abscisic acid sensors. Science 324, 1064–1068. doi: 10.1126/science.1172408
Melcher, K., Ng, L.-M., Zhou, X. E., Soon, F.-F., Xu, Y., Suino-Powell, K. M., et al. (2009). A gate–latch–lock mechanism for hormone signalling by abscisic acid receptors. Nature 462, 602–608. doi: 10.1038/nature08613
Merlot, S., Gosti, F., Guerrier, D., Vavasseur, A., and Giraudat, J. (2001). The ABI1 and ABI2 protein phosphatases 2C act in a negative feedback regulatory loop of the abscisic acid signalling pathway. Plant J. Cell Mol. Biol. 25, 295–303. doi: 10.1046/j.1365-313x.2001.00965.x
Narusaka, Y., Nakashima, K., Shinwari, Z. K., Sakuma, Y., Furihata, T., Abe, H., et al. (2003). Interaction between two cis-acting elements, ABRE and DRE, in ABA-dependent expression of Arabidopsis rd29A gene in response to dehydration and high-salinity stresses. Plant J. 34, 137–148. doi: 10.1046/j.1365-313X.2003.01708.x
Nishimura, N., Hitomi, K., Arvai, A. S., Rambo, R. P., Hitomi, C., Cutler, S. R., et al. (2009). Structural mechanism of Abscisic acid binding and signaling by Dimeric PYR1. Science 326, 1373–1379. doi: 10.1126/science.1181829
Nishimura, N., Sarkeshik, A., Nito, K., Park, S.-Y., Wang, A., Carvalho, P. C., et al. (2010). PYR/PYL/RCAR family members are major in-vivo ABI1 protein phosphatase 2C-interacting proteins in Arabidopsis. Plant J. Cell Mol. Biol. 61, 290–299. doi: 10.1111/j.1365-313X.2009.04054.x
Pan, C., Tang, J., Xu, Y., Xiao, P., Liu, H., Wang, H., et al. (2015). The catalytic role of the M2 metal ion in PP2Cα. Sci. Rep. 5, 1–9. doi: 10.1038/srep08560
Park, S.-Y., Fung, P., Nishimura, N., Jensen, D. R., Fujii, H., Zhao, Y., et al. (2009). Abscisic acid inhibits type 2C protein phosphatases via the PYR/PYL family of START proteins. Science 324, 1068–1071. doi: 10.1126/science.1173041
Poolman, M. G., Assmus, H. E., and Fell, D. A. (2004). Applications of metabolic modelling to plant metabolism. J. Exp. Bot. 55, 1177–1186. doi: 10.1093/jxb/erh090
Qin, F., Sakuma, Y., Tran, L.-S. P., Maruyama, K., Kidokoro, S., Fujita, Y., et al. (2008). Arabidopsis DREB2A-interacting proteins function as RING E3 ligases and negatively regulate plant drought stress–responsive gene expression. Plant Cell 20, 1693–1707. doi: 10.1105/tpc.107.057380
Rodriguez, P. L., Leube, M. P., and Grill, E. (1998). Molecular cloning in Arabidopsis thaliana of a new protein phosphatase 2C (PP2C) with homology to ABI1 and ABI2. Plant Mol. Biol. 38, 879–883. doi: 10.1023/a:1006012218704
Rubio, S., Rodrigues, A., Saez, A., Dizon, M. B., Galle, A., Kim, T.-H., et al. (2009). Triple loss of function of protein phosphatases type 2C leads to partial constitutive response to endogenous Abscisic acid. Plant Physiol. 150, 1345–1355. doi: 10.1104/pp.109.137174
Saez, A., Apostolova, N., Gonzalez-Guzman, M., Gonzalez-Garcia, M. P., Nicolas, C., Lorenzo, O., et al. (2004). Gain-of-function and loss-of-function phenotypes of the protein phosphatase 2C HAB1 reveal its role as a negative regulator of abscisic acid signalling. Plant J. Cell Mol. Biol. 37, 354–369. doi: 10.1046/j.1365-313x.2003.01966.x
Santiago, J., Rodrigues, A., Saez, A., Rubio, S., Antoni, R., Dupeux, F., et al. (2009). Modulation of drought resistance by the abscisic acid receptor PYL5 through inhibition of clade A PP2Cs. Plant J. Cell Mol. Biol. 60, 575–588. doi: 10.1111/j.1365-313X.2009.03981.x
Sauter, A., Davies, W. J., and Hartung, W. (2001). The long-distance abscisic acid signal in the droughted plant: the fate of the hormone on its way from root to shoot. J. Exp. Bot. 52, 1991–1997. doi: 10.1093/jexbot/52.363.1991
Schroeder, J. I., Hedrich, R., and Fernandez, J. M. (1984). Potassium-selective single channels in guard cell protoplasts of Vicia faba. Nature 312, 361–362. doi: 10.1038/312361a0
Shen, Q., Zhang, P., and Ho, T. H. (1996). Modular nature of abscisic acid (ABA) response complexes: Composite promoter units that are necessary and sufficient for ABA induction of gene expression in barley. Plant Cell 8, 1107–1119. doi: 10.1105/tpc.8.7.1107
Singh, D., and Laxmi, A. (2015). Transcriptional regulation of drought response: A tortuous network of transcriptional factors. Front. Plant Sci. 6:895. doi: 10.3389/fpls.2015.00895
Song, L., Huang, S. C., Wise, A., Castanon, R., Nery, J. R., Chen, H., et al. (2016). A transcription factor hierarchy defines an environmental stress response network. Science, 354. Available at: http://www.science.org/doi/abs/10.1126/science.aag1550 (Accessed September 8, 2021).
Soon, F.-F., Ng, L.-M., Zhou, X. E., West, G. M., Kovach, A., Tan, M. H. E., et al. (2012). Molecular mimicry regulates ABA signaling by SnRK2 kinases and PP2C phosphatases. Science 335, 85–88. doi: 10.1126/science.1215106
Steuer, B., Stuhlfauth, T., and Fock, H. P. (1988). The efficiency of water use in water stressed plants is increased due to ABA induced stomatal closure. Photosynth. Res. 18, 327–336. doi: 10.1007/BF00034837
Sun, Y., Oh, D.-H., Duan, L., Ramachandran, P., Ramirez, A., Bartlett, A., et al. (2022). Divergence in the ABA gene regulatory network underlies differential growth control. Nat. Plants 8, 549–560. doi: 10.1038/s41477-022-01139-5
Takahashi, Y., Zhang, J., Hsu, P.-K., Ceciliato, P. H. O., Zhang, L., Dubeaux, G., et al. (2020). MAP3Kinase-dependent SnRK2-kinase activation is required for abscisic acid signal transduction and rapid osmotic stress response. Nat. Commun. 11:12. doi: 10.1038/s41467-019-13875-y
Thakar, J., Pilione, M., Kirimanjeswara, G., Harvill, E. T., and Albert, R. (2007). Modeling systems-level regulation of host immune responses. PLoS Comput. Biol. 3:e109. doi: 10.1371/journal.pcbi.0030109
Thalmann, M., Pazmino, D., Seung, D., Horrer, D., Nigro, A., Meier, T., et al. (2016). Regulation of leaf starch degradation by Abscisic acid is important for osmotic stress tolerance in plants [OPEN]. Plant Cell 28, 1860–1878. doi: 10.1105/tpc.16.00143
Tischer, S. V., Wunschel, C., Papacek, M., Kleigrewe, K., Hofmann, T., Christmann, A., et al. (2017). Combinatorial interaction network of abscisic acid receptors and coreceptors from Arabidopsis thaliana. Proc. Natl. Acad. Sci. 114, 10280–10285. doi: 10.1073/pnas.1706593114
Umezawa, T., Sugiyama, N., Mizoguchi, M., Hayashi, S., Myouga, F., Yamaguchi-Shinozaki, K., et al. (2009). Type 2C protein phosphatases directly regulate abscisic acid-activated protein kinases in Arabidopsis. Proc. Natl. Acad. Sci. 106, 17588–17593. doi: 10.1073/pnas.0907095106
Uno, Y., Furihata, T., Abe, H., Yoshida, R., Shinozaki, K., and Yamaguchi-Shinozaki, K. (2000). Arabidopsis basic leucine zipper transcription factors involved in an abscisic acid-dependent signal transduction pathway under drought and high-salinity conditions. Proc. Natl. Acad. Sci. 97, 11632–11637. doi: 10.1073/pnas.190309197
Wang, X., Guo, C., Peng, J., Li, C., Wan, F., Zhang, S., et al. (2019). ABRE-BINDING FACTORS play a role in the feedback regulation of ABA signaling by mediating rapid ABA induction of ABA co-receptor genes. New Phytol. 221, 341–355. doi: 10.1111/nph.15345
Xie, T., Ren, R., Zhang, Y., Pang, Y., Yan, C., Gong, X., et al. (2012). Molecular mechanism for inhibition of a critical component in the Arabidopsis thaliana Abscisic acid signal transduction pathways, SnRK2.6, by protein phosphatase ABI1. J. Biol. Chem. 287, 794–802. doi: 10.1074/jbc.M111.313106
Yamaguchi-Shinozaki, K., and Shinozaki, K. (1994). A novel cis-acting element in an Arabidopsis gene is involved in responsiveness to drought, low-temperature, or high-salt stress. Plant Cell 6, 251–264. doi: 10.1105/tpc.6.2.251
Yamaguchi-Shinozaki, K., Urao, T., and Shinozaki, K. (1995). Regulation of genes that are induced by drought stress in Arabidopsis thaliana. J. Plant Res. 108, 127–136. doi: 10.1007/BF02344316
Yin, P., Fan, H., Hao, Q., Yuan, X., Wu, D., Pang, Y., et al. (2009). Structural insights into the mechanism of abscisic acid signaling by PYL proteins. Nat. Struct. Mol. Biol. 16, 1230–1236. doi: 10.1038/nsmb.1730
Yoshida, T., Fujita, Y., Maruyama, K., Mogami, J., Todaka, D., Shinozaki, K., et al. (2015). Four Arabidopsis AREB/ABF transcription factors function predominantly in gene expression downstream of SnRK2 kinases in abscisic acid signalling in response to osmotic stress. Plant Cell Environ. 38, 35–49. doi: 10.1111/pce.12351
Yoshida, T., Mogami, J., and Yamaguchi-Shinozaki, K. (2014). ABA-dependent and ABA-independent signaling in response to osmotic stress in plants. Curr. Opin. Plant Biol. 21, 133–139. doi: 10.1016/j.pbi.2014.07.009
Zeevaart, J. A. D., and Creelman, R. A. (1988). Metabolism and physiology of Abscisic acid. Annu. Rev. Plant Physiol. Plant Mol. Biol. 39, 439–473. doi: 10.1146/annurev.pp.39.060188.002255
Zhan, X., Wang, B., Li, H., Liu, R., Kalia, R. K., Zhu, J.-K., et al. (2012). Arabidopsis proline-rich protein important for development and abiotic stress tolerance is involved in micro RNA biogenesis. Proc. Natl. Acad. Sci. 109, 18198–18203. doi: 10.1073/pnas.1216199109
Zhang, W., Ruan, J., Ho, T.-H. D., You, Y., Yu, T., and Quatrano, R. S. (2005). Cis-regulatory element based targeted gene finding: Genome-wide identification of abscisic acid-and abiotic stress-responsive genes in Arabidopsis thaliana. Bioinforma. Oxf. Engl. 21, 3074–3081. doi: 10.1093/bioinformatics/bti490
Keywords: ODEs, ABA, gene regulatory network, ABRE, DRE, RD29A, mathematical model, feedback regulation
Citation: Ndathe R, Dale R and Kato N (2022) Dynamic modeling of ABA-dependent expression of the Arabidopsis RD29A gene. Front. Plant Sci. 13:928718. doi: 10.3389/fpls.2022.928718
Edited by:
Pil Joon Seo, Seoul National University, South KoreaReviewed by:
Takuya Yoshida, Technical University of Munich, GermanyRyosuke Mega, Yamaguchi University, Japan
Copyright © 2022 Ndathe, Dale and Kato. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Naohiro Kato, a2F0b0Bsc3UuZWR1