Boolean modeling of transcriptome data reveals novel modes of heterotrimeric G-protein action

Preparation of isolated guard cells, containing cell walls, from Vicia faba

Stomatal movement, initiated by specialized epidermal cells known as guard cells (GCs), plays a pivotal role in plant gas exchange and water use efficiency. Despite protocols existing for isolating GCs through proplasting for carrying out biochemical, physiological, and molecular studies, protocals for isolating GCs with their cell walls still intact have been lacking in the literature. In this paper, we introduce a method for the isolation of complete GCs from Vicia faba and show their membrane to remain impermeable through propidium iodide staining. This methodology enables further in-depth analyses into the cell wall composition of GCs, facilitating our understanding of structure-function relationship governing reversible actuation within cells.

Moving beyond the arabidopsis-centric view of G-protein signaling in plants

Heterotrimeric G-protein-mediated signaling is a key mechanism to transduce a multitude of endogenous and environmental signals in diverse organisms. The scope and expectations of plant G-protein research were set by pioneering work in metazoans. Given the similarity of the core constituents, G-protein-signaling mechanisms were presumed to be universally conserved. However, because of the enormous diversity of survival strategies and endless forms among eukaryotes, the signal, its interpretation, and responses vary even among different plant groups. Earlier G-protein research in arabidopsis (Arabidopsis thaliana) has emphasized its divergence from Metazoa. Here, we compare recent evidence from diverse plant lineages with the available arabidopsis G-protein model and discuss the conserved and novel protein components, signaling mechanisms, and response regulation.

PLC1 mediated Cycloastragenol-induced stomatal movement by regulating the production of NO in Arabidopsis thaliana

BACKGROUND

Astragalus grows mainly in drought areas. Cycloastragenol (CAG) is a tetracyclic triterpenoid allelochemical extracted from traditional Chinese medicine Astragalus root. Phospholipase C (PLC) and Gα-submit of the heterotrimeric G-protein (GPA1) are involved in many biotic or abiotic stresses. Nitric oxide (NO) is a crucial gas signal molecule in plants.

RESULTS

In this study, using the seedlings of Arabidopsis thaliana (A. thaliana), the results showed that low concentrations of CAG induced stomatal closure, and high concentrations inhibited stomatal closure. 30 µmol·L-1 CAG significantly increased the relative expression levels of PLC1 and GPA1 and the activities of PLC and GTP hydrolysis. The stomatal aperture of plc1, gpa1, and plc1/gpa1 was higher than that of WT under CAG treatment. CAG increased the fluorescence intensity of NO in guard cells. Exogenous application of c-PTIO to WT significantly induced stomatal aperture under CAG treatment. CAG significantly increased the relative expression levels of NIA1 and NOA1. Mutants of noa1, nia1, and nia2 showed that NO production was mainly from NOA1 and NIA1 by CAG treatment. The fluorescence intensity of NO in guard cells of plc1, gpa1, and plc1/gpa1 was lower than WT, indicating that PLC1 and GPA1 were involved in the NO production in guard cells. There was no significant difference in the gene expression of PLC1 in WT, nia1, and noa1 under CAG treatment. The gene expression levels of NIA1 and NOA1 in plc1, gpa1, and plc1/gpa1 were significantly lower than WT, indicating that PLC1 and GPA1 were positively regulating NO production by regulating the expression of NIA1 and NOA1 under CAG treatment.

CONCLUSIONS

These results suggested that the NO accumulation was essential to induce stomatal closure under CAG treatment, and GPA1 and PLC1 acted upstream of NO.

FAR-RED INSENSITIVE 219 and phytochrome B corepress shade avoidance via modulating nuclear speckle formation

Plants can sense the shade from neighboring plants by detecting a reduction of the red:far-red light (R:FR) ratio. Phytochrome B (phyB) is the primary photoreceptor that perceives shade light and regulates jasmonic acid (JA) signaling. However, the molecular mechanisms underlying phyB and JA signaling integration in shade responses remain largely unknown. Here, we show the interaction of phyB and FAR-RED INSENSITIVE 219 (FIN219)/JASMONATE RESISTANT1 (JAR1) in a functional demand manner in Arabidopsis (Arabidopsis thaliana) seedling development. Genetic evidence and interaction studies indicated that phyB and FIN219 synergistically and negatively regulate shade-induced hypocotyl elongation. Moreover, phyB interacted with various isoforms of FIN219 under high and low R:FR light. Methyl jasmonate (MeJA) treatment, FIN219 mutation, and PHYBOE digalactosyldiacylglycerol synthase1-1 (dgd1-1) plants, which show increased levels of JA, altered the patterns of phyB-associated nuclear speckles under the same conditions. Surprisingly, PHYBOE dgd1-1 showed a shorter hypocotyl phenotype than its parental mutants under shade conditions. Microarray assays using PHYBOE and PHYBOE fin219-2 indicated that PHYB overexpression substantially affects defense response-related genes under shade light and coregulates expression of auxin-responsive genes with FIN219. Thus, our findings reveal that phyB substantially crosstalks with JA signaling through FIN219 to modulate seedling development under shade light.

Flavonols modulate plant development, signaling, and stress responses

Flavonols are plant-specialized metabolites with important functions in plant growth and development. Isolation and characterization of mutants with reduced flavonol levels, especially the transparent testa mutants in Arabidopsis thaliana, have contributed to our understanding of the flavonol biosynthetic pathway. These mutants have also uncovered the roles of flavonols in controlling development in above- and below-ground tissues, notably in the regulation of root architecture, guard cell signaling, and pollen development. In this review, we present recent progress made towards a mechanistic understanding of flavonol function in plant growth and development. Specifically, we highlight findings that flavonols act as reactive oxygen species (ROS) scavengers and inhibitors of auxin transport in diverse tissues and cell types to modulate plant growth and development and responses to abiotic stresses.

Mathematical model of the cell signaling pathway based on the extended Boolean network model with a stochastic process

BACKGROUND

In cell signaling pathways, proteins interact with each other to determine cell fate in response to either cell-extrinsic (micro-environmental) or intrinsic cues. One of the well-studied pathways, the mitogen-activated protein kinase (MAPK) signaling pathway, regulates cell processes such as differentiation, proliferation, apoptosis, and survival in response to various micro-environmental stimuli in eukaryotes. Upon micro-environmental stimulus, receptors on the cell membrane become activated. Activated receptors initiate a cascade of protein activation in the MAPK pathway. This activation involves protein binding, creating scaffold proteins, which are known to facilitate effective MAPK signaling transduction.

RESULTS

This paper presents a novel mathematical model of a cell signaling pathway coordinated by protein scaffolding. The model is based on the extended Boolean network approach with stochastic processes. Protein production or decay in a cell was modeled considering the stochastic process, whereas the protein-protein interactions were modeled based on the extended Boolean network approach. Our model fills a gap in the binary set applied to previous models. The model simultaneously considers the stochastic process directly. Using the model, we simulated a simplified mitogen-activated protein kinase (MAPK) signaling pathway upon stimulation of both a single receptor at the initial time and multiple receptors at several time points. Our simulations showed that the signal is amplified as it travels down to the pathway from the receptor, generating substantially amplified downstream ERK activity. The noise generated by the stochastic process of protein self-activity in the model was also amplified as the signaling propagated through the pathway.

CONCLUSIONS

The signaling transduction in a simplified MAPK signaling pathway could be explained by a mathematical model based on the extended Boolean network model with a stochastic process. The model simulations demonstrated signaling amplifications when it travels downstream, which was already observed in experimental settings. We also highlight the importance of stochastic activity in regulating protein inactivation.

Relative importance of composition structures and biologically meaningful logics in bipartite Boolean models of gene regulation

Boolean networks have been widely used to model gene networks. However, such models are coarse-grained to an extent that they abstract away molecular specificities of gene regulation. Alternatively, bipartite Boolean network models of gene regulation explicitly distinguish genes from transcription factors (TFs). In such bipartite models, multiple TFs may simultaneously contribute to gene regulation by forming heteromeric complexes, thus giving rise to composition structures. Since bipartite Boolean models are relatively recent, an empirical investigation of their biological plausibility is lacking. Here, we estimate the prevalence of composition structures arising through heteromeric complexes. Moreover, we present an additional mechanism where composition structures may arise as a result of multiple TFs binding to cis-regulatory regions and provide empirical support for this mechanism. Next, we compare the restriction in BFs imposed by composition structures and by biologically meaningful properties. We find that though composition structures can severely restrict the number of Boolean functions (BFs) driving a gene, the two types of minimally complex BFs, namely nested canalyzing functions (NCFs) and read-once functions (RoFs), are comparatively more restrictive. Finally, we find that composition structures are highly enriched in real networks, but this enrichment most likely comes from NCFs and RoFs.

Interplay between ARABIDOPSIS Gβ and WRKY transcription factors differentiates environmental stress responses

Different environmental stresses often evoke similar physiological disorders such as growth retardation; however, specific consequences reported among individual stresses indicate potential mechanisms to distinguish different stress types in plants. Here, we examined mechanisms to differentiate between stress types in Arabidopsis (Arabidopsis thaliana). Gene expression patterns recapitulating several abiotic stress responses suggested abscisic acid (ABA) as a mediator of the common stress response, while stress type-specific responses were related to metabolic adaptations. Transcriptome and metabolome analyses identified Arabidopsis Gβ (AGB1) mediating the common stress-responsive genes and primary metabolisms under nitrogen excess. AGB1 regulated the expressions of multiple WRKY transcription factors. Gene Ontology and mutant analyses revealed different roles among WRKYs: WRKY40 is involved in ABA and common stress responses, while WRKY75 regulates metabolic processes. The AGB1-WRKY signaling module controlled developmental plasticity in roots under nitrogen excess. Signal transmission from AGB1 to a selective set of WRKYs would be essential to evoke unique responses to different types of stresses.

Inference of a Boolean Network From Causal Logic Implications

Biological systems contain a large number of molecules that have diverse interactions. A fruitful path to understanding these systems is to represent them with interaction networks, and then describe flow processes in the network with a dynamic model. Boolean modeling, the simplest discrete dynamic modeling framework for biological networks, has proven its value in recapitulating experimental results and making predictions. A first step and major roadblock to the widespread use of Boolean networks in biology is the laborious network inference and construction process. Here we present a streamlined network inference method that combines the discovery of a parsimonious network structure and the identification of Boolean functions that determine the dynamics of the system. This inference method is based on a causal logic analysis method that associates a logic type (sufficient or necessary) to node-pair relationships (whether promoting or inhibitory). We use the causal logic framework to assimilate indirect information obtained from perturbation experiments and infer relationships that have not yet been documented experimentally. We apply this inference method to a well-studied process of hormone signaling in plants, the signaling underlying abscisic acid (ABA)—induced stomatal closure. Applying the causal logic inference method significantly reduces the manual work typically required for network and Boolean model construction. The inferred model agrees with the manually curated model. We also test this method by re-inferring a network representing epithelial to mesenchymal transition based on a subset of the information that was initially used to construct the model. We find that the inference method performs well for various likely scenarios of inference input information. We conclude that our method is an effective approach toward inference of biological networks and can become an efficient step in the iterative process between experiments and computations.

CBF4/DREB1D represses XERICO to attenuate ABA, osmotic and drought stress responses in Arabidopsis

Water stress can severely impact plant growth, productivity and yield. Consequently, plants have evolved various strategies through which they can respond and adapt to their environment. XERICO (XER) is a stress-responsive RING E3 ubiquitin ligase that modulates abscisic acid (ABA) levels and promotes drought tolerance when overexpressed. To better understand the biological role of XER in stress responses, we characterized a xer-1 hypomorphic mutant and a CRISPR/Cas9-induced xer-2 null mutant in Arabidopsis. Both xer mutant alleles exhibited increased drought sensitivity, supporting the results from overexpression studies. Furthermore, we discovered that both xer mutants have greater stomatal indices and that XER is expressed in epidermal cells, indicating that XER functions in the epidermis to repress stomatal development. To explore XER spatiotemporal and stress-dependent regulation, we conducted a yeast one-hybrid screen and found that CBF4/DREB1D associates with the XER 5' untranslated region (5'-UTR). We generated three cbf4 null mutants with CRISPR/Cas9 and showed that CBF4 negatively regulates ABA responses, promotes stomatal development and reduces drought tolerance, in contrast to the roles shown for XER. CBF4 is induced by ABA and osmotic stress, and localizes to the nucleus where it downregulates XER expression via the DRE element in its 5'-UTR. Lastly, genetic interaction studies confirmed that xer is epistatic to cbf4 in stomatal development and in ABA, osmotic and drought stress responses. We propose that the repression of XER by CBF4 functions to attenuate ABA signaling and stress responses to maintain a balance between plant growth and survival under adverse environmental conditions.

SymRK-dependent phosphorylation of Gα protein and its role in signaling during soybean (Glycine max) nodulation

Heterotrimeric G proteins, comprised of Gα, Gβ and Gγ subunits, influence signaling in most eukaryotes. In metazoans, G proteins are activated by G protein-coupled receptor (GPCR)-mediated GDP to GTP exchange on Gα; however, the role(s) of GPCRs in regulating plant G-protein signaling remains equivocal. Mounting evidence suggests the involvement of receptor-like kinases (RLKs) in regulating plant G-protein signaling, but their mechanistic details remain scarce. We have previously shown that during Glycine max (soybean) nodulation, the nod factor receptor 1 (NFR1) interacts with G-protein components and indirectly affects signaling. We explored the direct regulation of G-protein signaling by RLKs using protein-protein interactions, receptor-mediated in vitro phosphorylations and the effects of such phosphorylations on soybean nodule formation. Results presented in this study demonstrate a direct, phosphorylation-based regulation of Gα by symbiosis receptor kinase (SymRK). SymRKs interact with and phosphorylate Gα at multiple residues in vitro, including two in its active site, which abolishes GTP binding. Additionally, phospho-mimetic Gα fails to interact with Gβγ, potentially allowing for constitutive signaling by the freed Gβγ. These results uncover an unusual mechanism of G-protein cycle regulation in plants where the receptor-mediated phosphorylation of Gα not only affects its activity but also influences the availability of its signaling partners, thereby exerting a two-pronged check on signaling.

Phosphatidic Acid in Plant Hormonal Signaling: From Target Proteins to Membrane Conformations

Cells sense a variety of extracellular signals balancing their metabolism and physiology according to changing growth conditions. Plasma membranes are the outermost informational barriers that render cells sensitive to regulatory inputs. Membranes are composed of different types of lipids that play not only structural but also informational roles. Hormones and other regulators are sensed by specific receptors leading to the activation of lipid metabolizing enzymes. These enzymes generate lipid second messengers. Among them, phosphatidic acid (PA) is a well-known intracellular messenger that regulates various cellular processes. This lipid affects the functional properties of cell membranes and binds to specific target proteins leading to either genomic (affecting transcriptome) or non-genomic responses. The subsequent biochemical, cellular and physiological reactions regulate plant growth, development and stress tolerance. In the present review, we focus on primary (genome-independent) signaling events triggered by rapid PA accumulation in plant cells and describe the functional role of PA in mediating response to hormones and hormone-like regulators. The contributions of individual lipid signaling enzymes to the formation of PA by specific stimuli are also discussed. We provide an overview of the current state of knowledge and future perspectives needed to decipher the mode of action of PA in the regulation of cell functions.

Transcriptional regulation of KCS gene by bZIP29 and MYB70 transcription factors during ABA-stimulated wound suberization of kiwifruit (Actinidia deliciosa)

BACKGROUND

Our previous study has demonstrated that the transcription of AchnKCS involved in suberin biosynthesis was up-regulated by exogenous abscisic acid (ABA) during the wound suberization of kiwifruit, but the regulatory mechanism has not been fully elucidated.

RESULTS

Through subcellular localization analysis in this work, AchnbZIP29 and AchnMYB70 transcription factors were observed to be localized in the nucleus. Yeast one-hybrid and dual-luciferase assay proved the transcriptional activation of AchnMYB70 and transcriptional suppression of AchnbZIP29 on AchnKCS promoter. Furthermore, the transcription level of AchnMYB70 was enhanced by ABA during wound suberization of kiwifruit, but AchnbZIP29 transcription was reduced by ABA.

CONCLUSIONS

Therefore, it was believed that ABA enhanced the transcriptional activation of AchnMYB70 on AchnKCS by increasing AchnMYB70 expression. On the contrary, ABA relieved the inhibitory effect of AchnbZIP29 on transcription of AchnKCS by inhibiting AchnbZIP29 expression. These results gave further insight into the molecular regulatory network of ABA in wound suberization of kiwifruit.

Raf-like kinases and receptor-like (pseudo)kinase GHR1 are required for stomatal vapor pressure difference response

Stomatal pores close rapidly in response to low-air-humidity-induced leaf-to-air vapor pressure difference (VPD) increases, thereby reducing excessive water loss. The hydroactive signal-transduction mechanisms mediating high VPD-induced stomatal closure remain largely unknown. The kinetics of stomatal high-VPD responses were investigated by using time-resolved gas-exchange analyses of higher-order mutants in guard-cell signal-transduction branches. We show that the slow-type anion channel SLAC1 plays a relatively more substantial role than the rapid-type anion channel ALMT12/QUAC1 in stomatal VPD signaling. VPD-induced stomatal closure is not affected in mpk12/mpk4GC double mutants that completely disrupt stomatal CO2 signaling, indicating that VPD signaling is independent of the early CO2 signal-transduction pathway. Calcium imaging shows that osmotic stress causes cytoplasmic Ca2+ transients in guard cells. Nevertheless, osca1-2/1.3/2.2/2.3/3.1 Ca2+-permeable channel quintuple, osca1.3/1.7-channel double, cngc5/6-channel double, cngc20-channel single, cngc19/20crispr-channel double, glr3.2/3.3-channel double, cpk-kinase quintuple, cbl1/4/5/8/9 quintuple, and cbl2/3rf double mutants showed wild-type-like stomatal VPD responses. A B3-family Raf-like mitogen-activated protein (MAP)-kinase kinase kinase, M3Kδ5/RAF6, activates the OST1/SnRK2.6 kinase in plant cells. Interestingly, B3 Raf-kinase m3kδ5 and m3kδ1/δ5/δ6/δ7 (raf3/6/5/4) quadruple mutants, but not a 14-gene raf-kinase mutant including osmotic stress-linked B4-family Raf-kinases, exhibited slowed high-VPD responses, suggesting that B3-family Raf-kinases play an important role in stomatal VPD signaling. Moreover, high VPD-induced stomatal closure was impaired in receptor-like pseudokinase GUARD CELL HYDROGEN PEROXIDE-RESISTANT1 (GHR1) mutant alleles. Notably, the classical transient "wrong-way" VPD response was absent in ghr1 mutant alleles. These findings reveal genes and signaling mechanisms in the elusive high VPD-induced stomatal closing response pathway.

Single-cell-type transcriptomic analysis reveals distinct gene expression profiles in wheat guard cells in response to abscisic acid

Stomatal closure, driven by shrinking guard cells in response to the accumulation of abscisic acid (ABA) under drought stress, has a great impact on plant growth and environmental acclimation. However, the molecular regulatory mechanism underlying the turgor alteration of guard cells remains elusive, especially in cereal grasses. Here, we develop a modified enzyme digestion-based approach for the isolation of wheat (Triticum aestivum L.) guard cells. With this approach, we can remove mesophyll, pavement cells and subsidiary cells successively from the epidermis of the trichomeless coleoptile in wheat and preserve guard cells on the cuticle layers in an intact and physiologically active conditions. Using a robust single-cell-type RNA sequencing analysis, we discovered 9829 differentially expressed genes (DEGs) as significantly up- or down-regulated in guard cells in response to ABA treatment. Transcriptome analysis revealed a large percent of DEGs encoding multiple phytohormone signalling pathways, transporters, calcium signalling components, protein kinases and other ABA signalling-related proteins, which are primarily involved in key signalling pathways in ABA-regulated stomatal control and stress response. Our findings provide valuable resource for investigating the transcriptional regulatory mechanism underlying wheat guard cells in response to ABA.

Lipid Signaling through G Proteins

N-Acylethanolamine (NAE) signaling has received considerable attention in vertebrates as part of the endocannabinoid signaling system, where anandamide acts as a ligand for G protein-coupled cannabinoid receptors. Recent studies indicate that G proteins also are required for some types of NAE signaling in plants. The genetic ablation of the Gβγ dimer or loss of the full set of extra-large G proteins strongly attenuated NAE-induced chloroplast responses in seedlings. Intriguing parallels and distinct differences have emerged between plants and animals in NAE signaling, despite the conserved use of these lipid mediators to modulate cellular processes. Here we compare similarities and differences and identify open questions in a fundamental lipid signaling pathway in eukaryotes with components that are both conserved and diverged in plants.

Seedling Chloroplast Responses Induced by N-Linolenoylethanolamine Require Intact G-Protein Complexes

In animals, several long-chain N-acylethanolamines (NAEs) have been identified as endocannabinoids and are autocrine signals that operate through cell surface G-protein-coupled cannabinoid receptors. Despite the occurrence of NAEs in land plants, including nonvascular plants, their precise signaling properties and molecular targets are not well defined. Here we show that the activity of N-linolenoylethanolamine (NAE 18:3) requires an intact G-protein complex. Specifically, genetic ablation of the Gβγ dimer or loss of the full set of atypical Gα subunits strongly attenuates an NAE-18:3-induced degreening of cotyledons in Arabidopsis (Arabidopsis thaliana) seedlings. This effect involves, at least in part, transcriptional regulation of chlorophyll biosynthesis and catabolism genes. In addition, there is feedforward transcriptional control of G-protein signaling components and G-protein interactors. These results are consistent with NAE 18:3 being a lipid signaling molecule in plants with a requirement for G-proteins to mediate signal transduction, a situation similar, but not identical, to the action of NAE endocannabinoids in animal systems.

Ureide metabolism in Arabidopsis thaliana is modulated by C:N balance

Plants can respond and adapt to changes in the internal content of carbon and nitrogen by using organic compounds that widely differ in their carbon/nitrogen ratio. Among them, the amides asparagine and glutamine are believed to be preferred by most plants, including Arabidopsis. However, increases in the ureides allantoin and/or allantoate concentrations have been observed in different plant species under several environmental conditions. In this work, changes in the ratio between carbon skeletons and reduced nitrogen were investigated by varying the concentrations of nitrogen and sucrose in the growth media. Allantoin accumulation was observed when plants were grown in media with high ammonia concentrations. This increase was reverted by adding sucrose as additional carbon source. Moreover, mutant plants with a decreased capability to degrade allantoin showed a compromised growth compared to WT in ammonia supplemented media. Together, our results indicate that allantoin accumulation is induced by low carbon/nitrogen ratio and suggest that its degradation is critical for proper plant growth and development.

Molecular changes in Mesembryanthemum crystallinum guard cells underlying the C3 to CAM transition

KEY MESSAGE: The timing and transcriptomic changes during the C₃ to CAM transition of common ice plant support the notion that guard cells themselves can shift from C₃ to CAM. Crassulacean acid metabolism (CAM) is a specialized type of photosynthesis: stomata close during the day, enhancing water conservation, and open at night, allowing CO₂ uptake. Mesembryanthemum crystallinum (common ice plant) is a facultative CAM species that can shift from C₃ photosynthesis to CAM under salt or drought stresses. However, the molecular mechanisms underlying the stress induced transition from C₃ to CAM remain unknown. Here we determined the transition time from C₃ to CAM in M. crystallinum under salt stress. In parallel, single-cell-type transcriptomic profiling by 3'-mRNA sequencing was conducted in isolated stomatal guard cells to determine the molecular changes in this key cell type during the transition. In total, 495 transcripts showed differential expression between control and salt-treated samples during the transition, including 285 known guard cell genes, seven CAM-related genes, 18 transcription factors, and 185 other genes previously not found to be expressed in guard cells. PEPC1 and PPCK1, which encode key enzymes of CAM photosynthesis, were up-regulated in guard cells after seven days of salt treatment, indicating that guard cells themselves can shift from C₃ to CAM. This study provides important information towards introducing CAM stomatal behavior into C₃ crops to enhance water use efficiency.

Convergent Loss of an EDS1/PAD4 Signaling Pathway in Several Plant Lineages Reveals Coevolved Components of Plant Immunity and Drought Response

Plant innate immunity relies on nucleotide binding leucine-rich repeat receptors (NLRs) that recognize pathogen-derived molecules and activate downstream signaling pathways. We analyzed the variation in NLR gene copy number and identified plants with a low number of NLR genes relative to sister species. We specifically focused on four plants from two distinct lineages, one monocot lineage (Alismatales) and one eudicot lineage (Lentibulariaceae). In these lineages, the loss of NLR genes coincides with loss of the well-known downstream immune signaling complex ENHANCED DISEASE SUSCEPTIBILITY 1 (EDS1)/PHYTOALEXIN DEFICIENT 4 (PAD4). We expanded our analysis across whole proteomes and found that other characterized immune genes were absent only in Lentibulariaceae and Alismatales. Additionally, we identified genes of unknown function that were convergently lost together with EDS1/PAD4 in five plant species. Gene expression analyses in Arabidopsis (Arabidopsis thaliana) and Oryza sativa revealed that several homologs of the candidates are differentially expressed during pathogen infection, drought, and abscisic acid treatment. Our analysis provides evolutionary evidence for the rewiring of plant immunity in some plant lineages, as well as the coevolution of the EDS1/PAD4 pathway and drought responses.

Overexpression of Isoprene Synthase Affects ABA- and Drought-Related Gene Expression and Enhances Tolerance to Abiotic Stress

Isoprene is the most abundant single biogenic volatile compound emitted by plants. Despite the relevance of this molecule to plant abiotic resistance and its impact on global atmospheric chemistry, little is known about the details of its mechanism of action. Here, we characterized through both physiological and molecular methods the mechanisms of action of isoprene using model transgenic arabidopsis lines overexpressing a monocot isoprene synthase gene. Our results demonstrated the effect that isoprene had on ABA signaling at different tissue-specific, spatial, and temporal scales. In particular, we found that isoprene enhanced stomatal sensitivity to ABA through upregulation of RD29B signaling gene. By contrast, isoprene decreased sensitivity to ABA in germinating seeds and roots, suggesting tissue-specific mechanisms of action. In leaves, isoprene caused the downregulation of COR15A and P5CS genes, suggesting that the enhanced tolerance to water-deprivation stress observed in isoprene-emitting plants may be mediated chiefly by an enhanced membrane integrity and tolerance to osmotic stress.

Crosstalk between heterotrimeric G protein-coupled signaling pathways and WRKY transcription factors modulating plant responses to suboptimal micronutrient conditions

Nutrient stresses induce foliar chlorosis and growth defects. Here we propose heterotrimeric G proteins as signaling mediators of various nutrient stresses, through meta-analyses of >20 transcriptomic data sets associated with nutrient stresses or G protein mutations. Systematic comparison of transcriptomic data yielded 104 genes regulated by G protein subunits under common nutrient stresses: 69 genes under Gβ subunit (AGB1) control and 35 genes under Gα subunit (GPA1) control. Quantitative real-time PCR experiments validate that several transcription factors and metal transporters changed in expression level under suboptimal iron, zinc, and/or copper concentrations, while being misregulated in the Arabidopsis Gβ-null (agb1) mutant. The agb1 mutant had altered metal ion profiles and exhibited severe growth arrest under zinc stress, and aberrant root waving under iron and zinc stresses, while the Gα-null mutation attenuated leaf chlorosis under iron deficiency in both Arabidopsis and rice. Our transcriptional network analysis inferred computationally that WRKY-family transcription factors mediate the AGB1-dependent nutrient responses. As corroborating evidence of our inference, ectopic expression of WRKY25 or WRKY33 rescued the zinc stress phenotypes and the expression of zinc transporters in the agb1-2 background. These results, together with Gene Ontology analyses, suggest two contrasting roles for G protein-coupled signaling pathways in micronutrient stress responses: one enhancing general stress tolerance and the other modulating ion homeostasis through WRKY transcriptional regulatory networks. In addition, tolerance to iron stress in the rice Gα mutant provides an inroad to improve nutrient stress tolerance of agricultural crops by manipulating G protein signaling.

ROS of Distinct Sources and Salicylic Acid Separate Elevated CO2-Mediated Stomatal Movements in Arabidopsis

Elevated CO2 (eCO2) often reduces leaf stomatal aperture and density thus impacts plant physiology and productivity. We have previously demonstrated that the Arabidopsis BIG protein distinguishes between the processes of eCO2-induced stomatal closure and eCO2-inhibited stomatal opening. However, the mechanistic basis of this action is not fully understood. Here we show that eCO2-elicited reactive oxygen species (ROS) production in big mutants was compromised in stomatal closure induction but not in stomatal opening inhibition. Pharmacological and genetic studies show that ROS generated by both NADPH oxidases and cell wall peroxidases contribute to eCO2-induced stomatal closure, whereas inhibition of light-induced stomatal opening by eCO2 may rely on the ROS derived from NADPH oxidases but not from cell wall peroxidases. As with JA and ABA, SA is required for eCO2-induced ROS generation and stomatal closure. In contrast, none of these three signals has a significant role in eCO2-inhibited stomatal opening, unveiling the distinct roles of plant hormonal signaling pathways in the induction of stomatal closure and the inhibition of stomatal opening by eCO2. In conclusion, this study adds SA to a list of plant hormones that together with ROS from distinct sources distinguish two branches of eCO2-mediated stomatal movements.

Dynamic control of plant water use using designed ABA receptor agonists

Drought causes crop losses worldwide, and its impact is expected to increase as the world warms. This has motivated the development of small-molecule tools for mitigating the effects of drought on agriculture. We show here that current leads are limited by poor bioactivity in wheat, a widely grown staple crop, and in tomato. To address this limitation, we combined virtual screening, x-ray crystallography, and structure-guided design to develop opabactin (OP), an abscisic acid (ABA) mimic with up to an approximately sevenfold increase in receptor affinity relative to ABA and up to 10-fold greater activity in vivo. Studies in Arabidopsis thaliana reveal a role of the type III receptor PYRABACTIN RESISTANCE-LIKE 2 for the antitranspirant efficacy of OP. Thus, virtual screening and structure-guided optimization yielded newly discovered agonists for manipulating crop abiotic stress tolerance and water use.

Flexible functional interactions between G-protein subunits contribute to the specificity of plant responses

Plants being sessile integrate information from a variety of endogenous and external cues simultaneously to optimize growth and development. This necessitates the signaling networks in plants to be highly dynamic and flexible. One such network involves heterotrimeric G-proteins comprised of Gα, Gβ, and Gγ subunits, which influence many aspects of growth, development, and stress response pathways. In plants such as Arabidopsis, a relatively simple repertoire of G-proteins comprised of one canonical and three extra-large Gα, one Gβ and three Gγ subunits exists. Because the Gβ and Gγ proteins form obligate dimers, the phenotypes of plants lacking the sole Gβ or all Gγ genes are similar, as expected. However, Gα proteins can exist either as monomers or in a complex with Gβγ, and the details of combinatorial genetic and physiological interactions of different Gα proteins with the sole Gβ remain unexplored. To evaluate such flexible, signal-dependent interactions and their contribution toward eliciting a specific response, we have generated Arabidopsis mutants lacking specific combinations of Gα and Gβ genes, performed extensive phenotypic analysis, and evaluated the results in the context of subunit usage and interaction specificity. Our data show that multiple mechanistic modes, and in some cases complex epistatic relationships, exist depending on the signal-dependent interactions between the Gα and Gβ proteins. This suggests that, despite their limited numbers, the inherent flexibility of plant G-protein networks provides for the adaptability needed to survive under continuously changing environments.

Arabidopsis Lectin EULS3 Is Involved in ABA Signaling in Roots

The Arabidopsis thaliana lectin ArathEULS3 is upregulated in particular stress conditions and upon abscisic acid (ABA) treatment. ABA is a plant hormone important for plant growth and stress responses. During stress ABA is perceived by PYR/PYL/RCAR receptors, inhibiting protein phosphatases PP2Cs thereby enabling SNRK2s kinases to start downstream phosphorylation cascades and signaling. PYL9, one of the ABA receptors was identified as an interacting partner for ArathEULS3. Promoter::GUS activity studies revealed the expression of ArathEULS3 in the central root cylinder and the cells flanking young lateral root primordia, and showed enhanced expression in root tips after ABA treatment. Transcript levels for ArathEULS3 increased after exposure to ABA and osmotic treatments. ArathEULS3 CRISPR KO mutants served as a tool to expand the knowledge on the role of ArathEULS3 in plant development. KO lines revealed a longer root system compared to WT plants, and showed reduced sensitivity to ABA, salt, and osmotic conditions. Additionally it was noted that the KO mutants had more emerged lateral roots when grown in high osmotic conditions. Together these data suggest that ArathEULS3 may be an important player in ABA responses in roots.

The Arabidopsis heterotrimeric G-protein β subunit, AGB1, is required for guard cell calcium sensing and calcium-induced calcium release

Cytosolic calcium concentration ([Ca²⁺ ]_cyt ) and heterotrimeric G-proteins are universal eukaryotic signaling elements. In plant guard cells, extracellular calcium (Ca_o ) is as strong a stimulus for stomatal closure as the phytohormone abscisic acid (ABA), but underlying mechanisms remain elusive. Here, we report that the sole Arabidopsis heterotrimeric Gβ subunit, AGB1, is required for four guard cell Ca_o responses: induction of stomatal closure; inhibition of stomatal opening; [Ca²⁺ ]_cyt oscillation; and inositol 1,4,5-trisphosphate (InsP3) production. Stomata in wild-type Arabidopsis (Col) and in mutants of the canonical Gα subunit, GPA1, showed inhibition of stomatal opening and promotion of stomatal closure by Ca_o . By contrast, stomatal movements of agb1 mutants and agb1/gpa1 double-mutants, as well as those of the agg1agg2 Gγ double-mutant, were insensitive to Ca_o . These behaviors contrast with ABA-regulated stomatal movements, which involve GPA1 and AGB1/AGG3 dimers, illustrating differential partitioning of G-protein subunits among stimuli with similar ultimate impacts, which may facilitate stimulus-specific encoding. AGB1 knockouts retained reactive oxygen species and NO production, but lost YC3.6-detected [Ca²⁺ ]_cyt oscillations in response to Ca_o , initiating only a single [Ca²⁺ ]_cyt spike. Experimentally imposed [Ca²⁺ ]_cyt oscillations restored stomatal closure in agb1. Yeast two-hybrid and bimolecular complementation fluorescence experiments revealed that AGB1 interacts with phospholipase Cs (PLCs), and Ca_o induced InsP3 production in Col but not in agb1. In sum, G-protein signaling via AGB1/AGG1/AGG2 is essential for Ca_o -regulation of stomatal apertures, and stomatal movements in response to Ca_o apparently require Ca²⁺ -induced Ca²⁺ release that is likely dependent on Gβγ interaction with PLCs leading to InsP3 production.

A Boolean network control algorithm guided by forward dynamic programming

Control problem in a biological system is the problem of finding an interventional policy for changing the state of the biological system from an undesirable state, e.g. disease, into a desirable healthy state. Boolean networks are utilized as a mathematical model for gene regulatory networks. This paper provides an algorithm to solve the control problem in Boolean networks. The proposed algorithm is implemented and applied on two biological systems: T-cell receptor network and Drosophila melanogaster network. Results show that the proposed algorithm works faster in solving the control problem over these networks, while having similar accuracy, in comparison to previous exact methods. Source code and a simple web service of the proposed algorithm is available at http://goliaei.ir/net-control/www/.

GCR1 and GPA1 coupling regulates nitrate, cell wall, immunity and light responses in Arabidopsis

G-protein signaling components have been attributed many biological roles in plants, but the extent of involvement of G-protein coupled receptor 1 (GCR1) with the Gα (GPA1) remained unknown. To address this, we have performed transcriptomic analyses on Arabidopsis gpa1-5gcr1-5 double mutant and identified 656 differentially expressed genes (DEGs). MapMan and Gene Ontology analyses revealed global transcriptional changes associated with external stimulus, cell wall organization/biogenesis and secondary metabolite process among others. Comparative transcriptomic analyses using the single and double mutants of gcr1-5 and gpa1-5 identified 194, 139 and 391 exclusive DEGs respectively, whereas 64 DEGs were common to all three mutants. Further, pair wise comparison of DEGs of double mutant with single mutants of gcr1-5 or gpa1-5 showed about one-third and over half common DEGs, respectively. Further analysis of the DEGs exclusive to the double mutant using protein-protein interaction networks revealed molecular complexes associated with nitrate and light signaling and plant-pathogen interactions among others. Physiological and molecular validation of nitrate-response revealed the sensitivity of germination to low N in the double mutant and differential expression of nitrate transporter (and nitrate reductase in all three mutants). Taken together, GCR1 and GPA1 work in partnership as well as independently to regulate different pathways.

The Role of Gβ Protein in Controlling Cell Expansion via Potential Interaction with Lipid Metabolic Pathways

Heterotrimeric G-proteins influence almost all aspects of plant growth, development, and responses to biotic and abiotic stresses in plants, likely via their interaction with specific effectors. However, the identity of such effectors and their mechanism of action are mostly unknown. While investigating the roles of different G-protein subunits in modulating the oil content in Camelina (Camelina sativa), an oil seed crop, we uncovered a role of Gβ proteins in controlling anisotropic cell expansion. Knockdown of Gβ genes causes reduced longitudinal and enhanced transverse expansion, resulting in altered cell, tissue, and organ shapes in transgenic plants during vegetative and reproductive development. These plants also exhibited substantial changes in their fatty acid and phospholipid profiles, which possibly leads to the increased oil content of the transgenic seeds. This increase is potentially caused by the direct interaction of Gβ proteins with a specific patatin-like phospholipase, pPLAIIIδ. Camelina plants with suppressed Gβ expression exhibit higher lipase activity, and show phenotypes similar to plants overexpressing pPLAIIIδ, suggesting that the Gβ proteins are negative regulators of pPLAIIIδ. These results reveal interactions between the G-protein-mediated and lipid signaling/metabolic pathways, where specific phospholipases may act as effectors that control key developmental and environmental responses of plants.

The Complex Fine-Tuning of K⁺ Fluxes in Plants in Relation to Osmotic and Ionic Abiotic Stresses

As the main cation in plant cells, potassium plays an essential role in adaptive responses, especially through its involvement in osmotic pressure and membrane potential adjustments. K⁺ homeostasis must, therefore, be finely controlled. As a result of different abiotic stresses, especially those resulting from global warming, K⁺ fluxes and plant distribution of this ion are disturbed. The hormone abscisic acid (ABA) is a key player in responses to these climate stresses. It triggers signaling cascades that ultimately lead to modulation of the activities of K⁺ channels and transporters. After a brief overview of transcriptional changes induced by abiotic stresses, this review deals with the post-translational molecular mechanisms in different plant organs, in Arabidopsis and species of agronomical interest, triggering changes in K⁺ uptake from the soil, K⁺ transport and accumulation throughout the plant, and stomatal regulation. These modifications involve phosphorylation/dephosphorylation mechanisms, modifications of targeting, and interactions with regulatory partner proteins. Interestingly, many signaling pathways are common to K⁺ and Cl^-/NO3^- counter-ion transport systems. These cross-talks are also addressed.

Insights into ABA-mediated regulation of guard cell primary metabolism revealed by systems biology approaches

Despite the fact that guard cell abscisic acid (ABA) signalling pathway is well documented, our understanding concerning how and to which extent ABA regulates guard cell metabolism remains fragmentary. Here we have adopted different systems approaches to investigate how ABA modulates guard cell central metabolism by providing genes that are possibly ABA-regulated. By using previous published Arabidopsis guard cell transcript profiling data, we carried out an extensive co-expression network analysis using ABA-related genes and those related to the metabolism and transport of sugars, starch and organic acids. Next, we investigated the presence of ABA responsive elements (ABRE) in the promoter of genes that are highly expressed in guard cells, responsive to ABA and co-expressed with ABA-related genes. Together, these analyses indicated that 44 genes are likely regulated by ABA and 8 of them are highly expressed in guard cells in both the presence and absence of ABA, including genes of the tricarboxylic acid cycle and those related to sucrose and hexose transport and metabolism. It seems likely that ABA may modulate both sucrose transport through guard cell plasma membrane and sucrose metabolism within guard cells. In this context, genes associated with sucrose synthase, sucrose phosphate synthase, trehalose-6-phosphate, invertase, UDP-glucose epimerase/pyrophosphorylase and different sugar transporters contain ABRE in their promoter and are thus possibly ABA regulated. Although validation experiments are required, our study highlights the importance of systems biology approaches to drive new hypothesis and to unravel genes and pathways that are regulated by ABA in guard cells.

Genetic and Systematic Approaches Toward G Protein-Coupled Abiotic Stress Signaling in Plants

Heterotrimeric G protein, composed of Gα, Gβ, and Gγ subunits, modulates plant adaptations to environmental stresses such as high salinity, drought, extreme temperatures and high light intensity. Most of these evidence were however derived solely from conventional genetics methods with which stress-associated phenotypes were compared between wild type and various G protein mutant plants. Recent advances in systematic approaches, mainly transcriptome and proteome, have contributed to in-depth understanding of molecular linkages between G proteins and environmental changes. Here, we update our knowledge on the roles of G proteins in abiotic stress responses. Furthermore, we highlight the current whole genome studies and integrated omics approach to better understand the fundamental G protein functions involved in abiotic stress responses. It is our purpose here to bridge the gap between molecular mechanisms in G protein science and stress biology and pave the way toward crop improvement researches in the future.

Emerging themes in heterotrimeric G-protein signaling in plants

Heterotrimeric G-proteins are key signaling components involved during the regulation of a multitude of growth and developmental pathways in all eukaryotes. Although the core proteins (Gα, Gβ, Gγ subunits) and their basic biochemistries are conserved between plants and non-plant systems, seemingly different inherent properties of specific components, altered wirings of G-protein network architectures, and the presence of novel receptors and effector proteins make plant G-protein signaling mechanisms somewhat distinct from the well-established animal paradigm. G-protein research in plants is getting a lot of attention recently due to the emerging roles of these proteins in controlling many agronomically important traits. New findings on both canonical and novel G-protein components and their conserved and unique signaling mechanisms are expected to improve our understanding of this important module in affecting critical plant growth and development pathways and eventually their utilization to produce plants for the future needs. In this review, we briefly summarize what is currently known in plant G-protein research, describe new findings and how they are changing our perceptions of the field, and discuss important issues that still need to be addressed.

The G Protein β-Subunit, AGB1, Interacts with FERONIA in RALF1-Regulated Stomatal Movement

Heterotrimeric guanine nucleotide-binding (G) proteins are composed of Gα, Gβ, and Gγ subunits and function as molecular switches in signal transduction. In Arabidopsis (Arabidopsis thaliana), there are one canonical Gα (GPA1), three extra-large Gα (XLG1, XLG2, and XLG3), one Gβ (AGB1), and three Gγ (AGG1, AGG2, and AGG3) subunits. To elucidate AGB1 molecular signaling, we performed immunoprecipitation using plasma membrane-enriched proteins followed by mass spectrometry to identify the protein interactors of AGB1. After eliminating proteins present in the control immunoprecipitation, commonly identified contaminants, and organellar proteins, a total of 103 candidate AGB1-associated proteins were confidently identified. We identified all of the G protein subunits except XLG1, receptor-like kinases, Ca²⁺ signaling-related proteins, and 14-3-3-like proteins, all of which may couple with or modulate G protein signaling. We confirmed physical interaction between AGB1 and the receptor-like kinase FERONIA (FER) using bimolecular fluorescence complementation. The Rapid Alkalinization Factor (RALF) family of polypeptides have been shown to be ligands of FER. In this study, we demonstrate that RALF1 regulates stomatal apertures and does so in a G protein-dependent manner, inhibiting stomatal opening and promoting stomatal closure in Columbia but not in agb1 mutants. We further show that AGGs and XLGs, but not GPA1, participate in RALF1-mediated stomatal signaling. Our results suggest that FER acts as a G protein-coupled receptor for plant heterotrimeric G proteins.

Preparation and applications of guard cell protoplasts from the leaf epidermis of Solanum lycopersicum

BACKGROUND

Guard cell protoplasts (GCPs) isolated from various plants have proven to be especially useful for studies of signal transduction pathways and plant development. But it is not easy to isolate highly purified preparations of large numbers of GCPs from plants. In this research, our focus is on a method to isolate large numbers of guard cells from tomato leaves. The protocols described yield millions of highly purified, viable GCPs, which are also suitable for studies on guard cell physiology.

RESULTS

We developed an efficient method for isolating GCPs from epidermal fragments of tomato leaves. The protocol requires a two-step digestion to isolate high-quality tomato GCPs. In this procedure, cellulysin (in method L) was replaced by cellulose "Onozuka" RS (in method S) in the first digestion step, which indicated that cellulase RS was more effective than cellulysin. Method S dramatically shortened the time required for obtaining high yields and high-quality GCPs. Moreover, according to the GCP yields, hydroponic plants were more effective than substrate-cultured plants.

CONCLUSIONS

In this paper, protocols for large-scale preparation of GCPs and mesophyll cell protoplasts were described, followed by some success examples of their use in biochemical and molecular approaches such as reverse-transcription polymerase chain reaction, real-time polymerase chain reaction and sodium dodecyl sulfate-polyacrylamide gel electrophoresis. The method was proved to be a more efficient GCP-isolating method, capable of providing high yields with better quality in less time.

Abscisic Acid-Induced Reactive Oxygen Species Are Modulated by Flavonols to Control Stomata Aperture

Abscisic acid (ABA) increases reactive oxygen species (ROS) in guard cells to close Arabidopsis (Arabidopsis thaliana) stomata. In tomato (Solanum lycopersicum), we find that ABA-increased ROS is followed by stomatal closure and that both responses are blocked by inhibitors of ROS-producing respiratory burst oxidase enzymes. ABA-induced ROS sensor fluorescence accumulates in the nucleus, chloroplasts, and endomembranes. The accumulation of flavonol antioxidants in guard cells, but not surrounding pavement cells, was visualized by confocal microscopy using a flavonol-specific fluorescent dye. Decreased flavonols in guard cells in the anthocyanin reduced (are) mutant and elevated levels in the anthocyanin without (aw) mutant were quantified by confocal microscopy and in leaf extracts by mass spectrometry. Consistent with flavonols acting as antioxidants, higher levels of ROS were detected in guard cells of the tomato are mutant and lower levels were detected in aw both at homeostasis and after treatment with ABA. These results demonstrate the inverse relationship between flavonols and ROS. Guard cells of are show greater ABA-induced closure than the wild type, reduced light-dependent guard cell opening, and reduced water loss, with aw having opposite responses. Ethylene treatment of wild-type tomato plants increased flavonol accumulation in guard cells; however, no flavonol increases were observed in Neverripe (Nr), an ethylene receptor mutant. Consistent with lower levels of ROS due to elevated flavonols, ethylene treatments decreased ABA-induced stomatal closure in the wild type, but not Nr, with ethylene responses attenuated in the are mutant. Together, these results are consistent with flavonols dampening the ABA-dependent ROS burst that drives stomatal closure and facilitating stomatal opening to modulate leaf gas exchange.

POLYGALACTURONASE INVOLVED IN EXPANSION3 Functions in Seedling Development, Rosette Growth, and Stomatal Dynamics in Arabidopsis thaliana

Plant cell separation and expansion require pectin degradation by endogenous pectinases such as polygalacturonases, few of which have been functionally characterized. Stomata are a unique system to study both processes because stomatal maturation involves limited separation between sister guard cells and stomatal responses require reversible guard cell elongation and contraction. However, the molecular mechanisms for how stomatal pores form and how guard cell walls facilitate dynamic stomatal responses remain poorly understood. We characterized POLYGALACTURONASE INVOLVED IN EXPANSION3 (PGX3), which is expressed in expanding tissues and guard cells. PGX3-GFP localizes to the cell wall and is enriched at sites of stomatal pore initiation in cotyledons. In seedlings, ablating or overexpressing PGX3 affects both cotyledon shape and the spacing and pore dimensions of developing stomata. In adult plants, PGX3 affects rosette size. Although stomata in true leaves display normal density and morphology when PGX3 expression is altered, loss of PGX3 prevents smooth stomatal closure, and overexpression of PGX3 accelerates stomatal opening. These phenotypes correspond with changes in pectin molecular mass and abundance that can affect wall mechanics. Together, these results demonstrate that PGX3-mediated pectin degradation affects stomatal development in cotyledons, promotes rosette expansion, and modulates guard cell mechanics in adult plants.

Recently duplicated plant heterotrimeric Gα proteins with subtle biochemical differences influence specific outcomes of signal-response coupling

Heterotrimeric G-proteins, comprising Gα, Gβ, and Gγ subunits, regulate key signaling processes in eukaryotes. The Gα subunit determines the status of signaling by switching between inactive GDP-bound and active GTP-bound forms. Unlike animal systems, in which multiple Gα proteins with variable biochemical properties exist, plants have fewer, highly similar Gα subunits that have resulted from recent genome duplications. These proteins exhibit subtle differences in their GTP-binding, GDP/GTP-exchange, and GTP-hydrolysis activities, but the extent to which these differences contribute to affect plant signaling and development remains unknown. To evaluate this, we expressed native and engineered Gα proteins from soybean in an Arabidopsis Gα-null background and studied their effects on modulating a range of developmental and hormonal signaling phenotypes. Our results indicated that inherent biochemical differences in these highly similar Gα proteins are biologically relevant, and some proteins are more flexible than others in influencing the outcomes of specific signals. These observations suggest that alterations in the rate of the G-protein cycle itself may contribute to the specificity of response regulation in plants by affecting the duration of active signaling and/or by the formation of distinct protein-protein complexes. In species such as Arabidopsis having a single canonical Gα, this rate could be affected by regulatory proteins in the presence of specific signals, whereas in plants with multiple Gα proteins, an even more complex regulation may exist, which likely contributes to the specificity of signal-response coupling.

Resolving the central metabolism of Arabidopsis guard cells

Photosynthesis and water use efficiency, key factors affecting plant growth, are directly controlled by microscopic and adjustable pores in the leaf-the stomata. The size of the pores is modulated by the guard cells, which rely on molecular mechanisms to sense and respond to environmental changes. It has been shown that the physiology of mesophyll and guard cells differs substantially. However, the implications of these differences to metabolism at a genome-scale level remain unclear. Here, we used constraint-based modeling to predict the differences in metabolic fluxes between the mesophyll and guard cells of Arabidopsis thaliana by exploring the space of fluxes that are most concordant to cell-type-specific transcript profiles. An independent 13C-labeling experiment using isolated mesophyll and guard cells was conducted and provided support for our predictions about the role of the Calvin-Benson cycle in sucrose synthesis in guard cells. The combination of in silico with in vivo analyses indicated that guard cells have higher anaplerotic CO2 fixation via phosphoenolpyruvate carboxylase, which was demonstrated to be an important source of malate. Beyond highlighting the metabolic differences between mesophyll and guard cells, our findings can be used in future integrated modeling of multi-cellular plant systems and their engineering towards improved growth.

Isolation of Guard-cell Enriched Tissue for RNA Extraction

This is a protocol for isolation of guard cell enriched samples from Arabidopsis thaliana plants for RNA extraction. Leaves are blended in ice-water and filtered through nylon mesh to obtain guard cell enriched fragments. With guard cell enriched samples, gene expression analysis can be done, e.g., comparing different gene expression levels in guard cells versus whole leaf to determine if a gene of interest is predominantly expressed in guard cells. It can also be used to study the effect of treatments or different genetic backgrounds in the regulation of the guard cell expressed genes.

ePlant: Visualizing and Exploring Multiple Levels of Data for Hypothesis Generation in Plant Biology

A big challenge in current systems biology research arises when different types of data must be accessed from separate sources and visualized using separate tools. The high cognitive load required to navigate such a workflow is detrimental to hypothesis generation. Accordingly, there is a need for a robust research platform that incorporates all data and provides integrated search, analysis, and visualization features through a single portal. Here, we present ePlant (http://bar.utoronto.ca/eplant), a visual analytic tool for exploring multiple levels of Arabidopsis thaliana data through a zoomable user interface. ePlant connects to several publicly available web services to download genome, proteome, interactome, transcriptome, and 3D molecular structure data for one or more genes or gene products of interest. Data are displayed with a set of visualization tools that are presented using a conceptual hierarchy from big to small, and many of the tools combine information from more than one data type. We describe the development of ePlant in this article and present several examples illustrating its integrative features for hypothesis generation. We also describe the process of deploying ePlant as an "app" on Araport. Building on readily available web services, the code for ePlant is freely available for any other biological species research.

Phosphatidic acid binding inhibits RGS1 activity to affect specific signaling pathways in Arabidopsis

Modulation of the active versus inactive forms of the Gα protein is critical for the signaling processes mediated by the heterotrimeric G-protein complex. We have recently established that in Arabidopsis, the regulator of G-protein signaling (RGS1) protein and a lipid-hydrolyzing enzyme, phospholipase Dα1 (PLDα1), both act as GTPase-activity accelerating proteins (GAPs) for the Gα protein to attenuate its activity. RGS1 and PLDα1 interact with each other, and RGS1 inhibits the activity of PLDα1 during regulation of a subset of responses. In this study, we present evidence that this regulation is bidirectional. Phosphatidic acid (PA), a second messenger typically derived from the lipid-hydrolyzing activity of PLDα1, is a molecular target of RGS1. PA binds and inhibits the GAP activity of RGS1. A conserved lysine residue in RGS1 (Lys²⁵⁹ ) is directly involved in RGS1-PA binding. Introduction of this RGS1 protein variant in the rgs1 mutant background makes plants hypersensitive to a subset of abscisic acid-mediated responses. Our data point to the existence of negative feedback loops between these two regulatory proteins that precisely modulate the level of active Gα, consequently generating a highly controlled signal-response output.

Cytokinin-Mediated Regulation of Reactive Oxygen Species Homeostasis Modulates Stomatal Immunity in Arabidopsis

Stomata play an important role in preinvasive defense responses by limiting pathogen entry into leaves. Although the stress hormones salicylic acid (SA) and abscisic acid (ABA) are known to regulate stomatal immunity, the role of growth promoting hormones is far from understood. Here, we show that in Arabidopsis thaliana, cytokinins (CKs) function in stomatal defense responses. The cytokinin receptor HISTIDINE KINASE3 (AHK3) and RESPONSE REGULATOR2 (ARR2) promote stomatal closure triggered by pathogen-associated molecular pattern (PAMP) and resistance to Pseudomonas syringae pv tomato bacteria. Importantly, the cytokinin trans-zeatin induces stomatal closure and accumulation of reactive oxygen species (ROS) in guard cells through AHK3 and ARR2 in an SA-dependent and ABA-independent manner. Using pharmacological and reverse genetics approaches, we found that CK-mediated stomatal responses involve the apoplastic peroxidases PRX4, PRX33, PRX34, and PRX71, but not the NADPH oxidases RBOHD and RBOHF. Moreover, ARR2 directly activates the expression of PRX33 and PRX34, which are required for SA- and PAMP-triggered ROS production. Thus, the CK signaling pathway regulates ROS homeostasis in guard cells, which leads to enhanced stomatal immunity and plant resistance to bacteria.

Functional characterization of the Arabidopsis transcription factor bZIP29 reveals its role in leaf and root development

Plant bZIP group I transcription factors have been reported mainly for their role during vascular development and osmosensory responses. Interestingly, bZIP29 has been identified in a cell cycle interactome, indicating additional functions of bZIP29 in plant development. Here, bZIP29 was functionally characterized to study its role during plant development. It is not present in vascular tissue but is specifically expressed in proliferative tissues. Genome-wide mapping of bZIP29 target genes confirmed its role in stress and osmosensory responses, but also identified specific binding to several core cell cycle genes and to genes involved in cell wall organization. bZIP29 protein complex analyses validated interaction with other bZIP group I members and provided insight into regulatory mechanisms acting on bZIP dimers. In agreement with bZIP29 expression in proliferative tissues and with its binding to promoters of cell cycle regulators, dominant-negative repression of bZIP29 altered the cell number in leaves and in the root meristem. A transcriptome analysis on the root meristem, however, indicated that bZIP29 might regulate cell number through control of cell wall organization. Finally, ectopic dominant-negative repression of bZIP29 and redundant factors led to a seedling-lethal phenotype, pointing to essential roles for bZIP group I factors early in plant development.

Allantoin accumulation mediated by allantoinase downregulation and transport by Ureide Permease 5 confers salt stress tolerance to Arabidopsis plants

Allantoin, a metabolite generated in the purine degradation pathway, was primarily considered an intermediate for recycling of the abundant nitrogen assimilated in plant purines. More specifically, tropical legumes utilize allantoin and allantoic acid as major nodule-to-shoot nitrogen transport compounds. In other species, an increase in allantoin content was observed under different stress conditions, but the underlying molecular mechanisms remain poorly understood. In this work, Arabidopsis thaliana was used as a model system to investigate the effects of salt stress on allantoin metabolism and to know whether its accumulation results in plant protection. Plant seedlings treated with NaCl at different concentrations showed higher allantoin and lower allantoic acid contents. Treatments with NaCl favored the expression of genes involved in allantoin synthesis, but strongly repressed the unique gene encoding allantoinase (AtALN). Due to the potential regulatory role of this gene for allantoin accumulation, AtALN promoter activity was studied using a reporter system. GUS mediated coloration was found in specific plant tissues and was diminished with increasing salt concentrations. Phenotypic analysis of knockout, knockdown and stress-inducible mutants for AtALN revealed that allantoin accumulation is essential for salt stress tolerance. In addition, the possible role of allantoin transport was investigated. The Ureide Permease 5 (UPS5) is expressed in the cortex and endodermis of roots and its transcription is enhanced by salt treatment. Ups5 knockout plants under salt stress presented a susceptible phenotype and altered allantoin root-to-shoot content ratios. Possible roles of allantoin as a protectant compound in oxidative events or signaling are discussed.

The role of PLDα1 in providing specificity to signal-response coupling by heterotrimeric G-protein components in Arabidopsis

Heterotrimeric G-proteins comprised of Gα, Gβ and Gγ subunits are important signal transducers in all eukaryotes. In plants, G-proteins affect multiple biotic and abiotic stress responses, as well as many developmental processes, even though their repertoire is significantly limited compared with that in metazoan systems. One canonical and three extra-large Gα, 1 Gβ and 3 Gγ proteins represent the heterotrimeric G-protein complex in Arabidopsis, and a single regulatory protein, RGS1, is one of the few known biochemical regulators of this signaling complex. This quantitative disparity between the number of signaling components and the range of processes they influence is rather intriguing. We now present evidence that the phospholipase Dα1 protein is a key component and modulator of the G-protein complex in affecting a subset of signaling pathways. We also show that the same G-protein subunits and their modulators exhibit distinct physiological and genetic interactions depending on specific signaling and developmental pathways. Such developmental plasticity and interaction specificity likely compensates for the lack of multiplicity of individual subunits, and helps to fine tune the plants' responses to constantly changing environments.

A Specific Transcriptome Signature for Guard Cells from the C4 Plant Gynandropsis gynandra

C4 photosynthesis represents an excellent example of convergent evolution that results in the optimization of both carbon and water usage by plants. In C4 plants, a carbon-concentrating mechanism divided between bundle sheath and mesophyll cells increases photosynthetic efficiency. Compared with C3 leaves, the carbon-concentrating mechanism of C4 plants allows photosynthetic operation at lower stomatal conductance, and as a consequence, transpiration is reduced. Here, we characterize transcriptomes from guard cells in C3 Tareneya hassleriana and C4 Gynandropsis gynandra belonging to the Cleomaceae. While approximately 60% of Gene Ontology terms previously associated with guard cells from the C3 model Arabidopsis (Arabidopsis thaliana) are conserved, there is much less overlap between patterns of individual gene expression. Most ion and CO2 signaling modules appear unchanged at the transcript level in guard cells from C3 and C4 species, but major variations in transcripts associated with carbon-related pathways known to influence stomatal behavior were detected. Genes associated with C4 photosynthesis were more highly expressed in guard cells of C4 compared with C3 leaves. Furthermore, we detected two major patterns of cell-specific C4 gene expression within the C4 leaf. In the first, genes previously associated with preferential expression in the bundle sheath showed continually decreasing expression from bundle sheath to mesophyll to guard cells. In the second, expression was maximal in the mesophyll compared with both guard cells and bundle sheath. These data imply that at least two gene regulatory networks act to coordinate gene expression across the bundle sheath, mesophyll, and guard cells in the C4 leaf.

Preparation of Epidermal Peels and Guard Cell Protoplasts for Cellular, Electrophysiological, and -Omics Assays of Guard Cell Function

Bioassays are commonly used to study stomatal phenotypes. There are multiple options in the choice of plant materials and species used for observation of stomatal and guard cell responses in vivo. Here, detailed procedures for bioassays of stomatal responses to abscisic acid (ABA) in Arabidopsis thaliana are described, including ABA promotion of stomatal closure, ABA inhibition of stomatal opening, and ABA promotion of reaction oxygen species (ROS) production in guard cells. We also include an example of a stomatal bioassay for the guard cell CO2 response using guard cell-enriched epidermal peels from Brassica napus. Highly pure preparations of guard cell protoplasts can be produced, which are also suitable for studies on guard cell signaling, as well as for studies on guard cell ion transport. Small-scale and large-scale guard cell protoplast preparations are commonly used for electrophysiological and -omics studies, respectively. We provide a procedure for small-scale guard cell protoplasting from A. thaliana. Additionally, a general protocol for large-scale preparation of guard cell protoplasts, with specifications for three different species, A. thaliana, B. napus, and Vicia faba is also provided.

Microarray Analysis of Rice d1 (RGA1) Mutant Reveals the Potential Role of G-Protein Alpha Subunit in Regulating Multiple Abiotic Stresses Such as Drought, Salinity, Heat, and Cold

The genome-wide role of heterotrimeric G-proteins in abiotic stress response in rice has not been examined from a functional genomics perspective, despite the availability of mutants and evidences involving individual genes/processes/stresses. Our rice whole transcriptome microarray analysis (GSE 20925 at NCBI GEO) using the G-alpha subunit (RGA1) null mutant (Daikoku 1 or d1) and its corresponding wild type (Oryza sativa Japonica Nipponbare) identified 2270 unique differentially expressed genes (DEGs). Out of them, we mined for all the potentially abiotic stress-responsive genes using Gene Ontology terms, STIFDB2.0 and Rice DB. The first two approaches produced smaller subsets of the 1886 genes found at Rice DB. The GO approach revealed similar regulation of several families of stress-responsive genes in RGA1 mutant. The Genevestigator analysis of the stress-responsive subset of the RGA1-regulated genes from STIFDB revealed cold and drought-responsive clusters. Meta data analysis at Rice DB revealed large stress-response categories such as cold (878 up/810 down), drought (882 up/837 down), heat (913 up/777 down), and salt stress (889 up/841 down). One thousand four hundred ninety-eight of them are common to all the four abiotic stresses, followed by fewer genes common to smaller groups of stresses. The RGA1-regulated genes that uniquely respond to individual stresses include 111 in heat stress, eight each in cold only and drought only stresses, and two genes in salt stress only. The common DEGs (1498) belong to pathways such as the synthesis of polyamine, glycine-betaine, proline, and trehalose. Some of the common DEGs belong to abiotic stress signaling pathways such as calcium-dependent pathway, ABA independent and dependent pathway, and MAP kinase pathway in the RGA1 mutant. Gene ontology of the common stress responsive DEGs revealed 62 unique molecular functions such as transporters, enzyme regulators, transferases, hydrolases, carbon and protein metabolism, binding to nucleotides, carbohydrates, receptors and lipids, morphogenesis, flower development, and cell homeostasis. We also mined 63 miRNAs that bind to the stress responsive transcripts identified in this study, indicating their post-transcriptional regulation. Overall, these results indicate the potentially extensive role of RGA1 in the regulation of multiple abiotic stresses in rice for further validation.