A network-based modeling framework reveals the core signal transduction network underlying high carbon dioxide-induced stomatal closure in guard cells.

Stomata are pores on plant aerial surfaces, each bordered by a pair of guard cells. They control gas exchange vital for plant survival. Understanding how guard cells respond to environmental signals such as atmospheric carbon dioxide (CO2) levels is not only insightful to fundamental biology but also relevant to real-world issues of crop productivity under global climate change. In the past decade, multiple important signaling elements for stomatal closure induced by elevated CO2 have been identified. Yet, there is no comprehensive understanding of high CO2-induced stomatal closure. In this work, we assemble a cellular signaling network underlying high CO2-induced stomatal closure by integrating evidence from a comprehensive literature analysis. We further construct a Boolean dynamic model of the network, which allows in silico simulation of the stomatal closure response to high CO2 in wild-type Arabidopsis thaliana plants and in cases of pharmacological or genetic manipulation of network nodes. Our model has a 91% accuracy in capturing known experimental observations. We perform network-based logical analysis and reveal a feedback core of the network, which dictates cellular decisions in closure response to high CO2. Based on these analyses, we predict and experimentally confirm that applying nitric oxide (NO) induces stomatal closure in ambient CO2 and causes hypersensitivity to elevated CO2. Moreover, we predict a negative regulatory relationship between NO and the protein phosphatase ABI2 and find experimentally that NO inhibits ABI2 phosphatase activity. The experimental validation of these model predictions demonstrates the effectiveness of network-based modeling and highlights the decision-making role of the feedback core of the network in signal transduction. We further explore the model's potential in predicting targets of signaling elements not yet connected to the CO2 network. Our combination of network science, in silico model simulation, and experimental assays demonstrates an effective interdisciplinary approach to understanding system-level biology.

Rock, scissors, paper: How RNA structure informs function.

RNA can fold back on itself to adopt a wide range of structures. These range from relatively simple hairpins to intricate 3D folds and can be accompanied by regulatory interactions with both metabolites and macromolecules. The last 50 yr have witnessed elucidation of an astonishing array of RNA structures including transfer RNAs, ribozymes, riboswitches, the ribosome, the spliceosome, and most recently entire RNA structuromes. These advances in RNA structural biology have deepened insight into fundamental biological processes including gene editing, transcription, translation, and structure-based detection and response to temperature and other environmental signals. These discoveries reveal that RNA can be relatively static, like a rock; that it can have catalytic functions of cutting bonds, like scissors; and that it can adopt myriad functional shapes, like paper. We relate these extraordinary discoveries in the biology of RNA structure to the plant way of life. We trace plant-specific discovery of ribozymes and riboswitches, alternative splicing, organellar ribosomes, thermometers, whole-transcriptome structuromes and pan-structuromes, and conclude that plants have a special set of RNA structures that confer unique types of gene regulation. We finish with a consideration of future directions for the RNA structure-function field.

Genome-wide analysis of the in vivo tRNA structurome reveals RNA structural and modification dynamics under heat stress.

RNA structure plays roles in myriad cellular events including transcription, translation, and RNA processing. Genome-wide analyses of RNA secondary structure in vivo by chemical probing have revealed critical structural features of mRNAs and long ncRNAs. Here, we examine the in vivo secondary structure of a small RNA class, tRNAs. Study of tRNA structure is challenging because tRNAs are heavily modified and strongly structured. We introduce "tRNA structure-seq," a new workflow that accurately determines in vivo secondary structures of tRNA. The workflow combines dimethyl sulfate (DMS) probing, ultra-processive RT, and mutational profiling (MaP), which provides mutations opposite DMS and natural modifications thereby allowing multiple modifications to be identified in a single read. We applied tRNA structure-seq to E. coli under control and stress conditions. A leading folding algorithm predicts E. coli tRNA structures with only ∼80% average accuracy from sequence alone. Strikingly, tRNA structure-seq, by providing experimental restraints, improves structure prediction under in vivo conditions to ∼95% accuracy, with more than 14 tRNAs predicted completely correctly. tRNA structure-seq also quantifies the relative levels of tRNAs and their natural modifications at single nucleotide resolution, as validated by LC-MS/MS. Our application of tRNA structure-seq yields insights into tRNA structure in living cells, revealing that it is not immutable but has dynamics, with partial unfolding of secondary and tertiary tRNA structure under heat stress that is correlated with a loss of tRNA abundance. This method is applicable to other small RNAs, including those with natural modifications and highly structured regions.

Flanking Sequence Cotranscriptionally Regulates Twister Ribozyme Activity.

Small nucleolytic ribozymes are RNAs that cleave their own phosphodiester backbone. While proteinaceous enzymes are regulated by a variety of known mechanisms, methods of regulation for ribozymes remain unclear. Twister is one ribozyme class for which many structural and catalytic properties have been elucidated. However, few studies have analyzed the activity of twister ribozymes in the context of a native flanking sequence, even though ribozymes as transcribed in nature do not exist in isolation. Interactions between the ribozyme and its neighboring sequences can induce conformational changes that inhibit self-cleavage, providing a regulatory mechanism that could naturally determine ribozyme activity in vivo and in synthetic applications. To date, eight twister ribozymes have been identified within the staple crop rice (Oryza sativa). Herein, we select several twister ribozymes from rice and show that they are differentially regulated by their flanking sequence using published RNA-seq data sets, structure probing, and cotranscriptional cleavage assays. We found that the Osa 1-2 ribozyme does not interact with its flanking sequences. However, sequences flanking the Osa 1-3 and Osa 1-8 ribozymes form inactive conformations, referred to here as "ribozymogens", that attenuate ribozyme self-cleavage activity. For the Osa 1-3 ribozyme, we show that activity can be rescued upon addition of a complementary antisense oligonucleotide, suggesting ribozymogens can be controlled via external signals. In all, our data provide a plausible mechanism wherein flanking sequence differentially regulates ribozyme activity in vivo. More broadly, the ability to regulate ribozyme behavior locally has potential applications in control of gene expression and synthetic biology.

Cantil: a previously unreported organ in wild-type Arabidopsis regulated by FT, ERECTA and heterotrimeric G proteins.

We describe a previously unreported macroscopic Arabidopsis organ, the cantil, named for its 'cantilever' function of holding the pedicel at a distance from the stem. Cantil development is strongest at the first nodes after the vegetative to reproductive inflorescence transition; cantil magnitude and frequency decrease acropetally. Cantils develop in wild-type Arabidopsis accessions (e.g. Col-0, Ws and Di-G) as a consequence of delayed flowering in short days; cantil formation is observed in long days when flowering is delayed by null mutation of the floral regulator FLOWERING LOCUS T. The receptor-like kinase ERECTA is a global positive regulator of cantil formation; therefore, cantils never form in the Arabidopsis strain Ler. ERECTA functions genetically upstream of heterotrimeric G proteins. Cantil expressivity is repressed by the specific heterotrimeric complex subunits GPA1, AGB1 and AGG3, which also play independent roles: GPA1 suppresses distal spurs at cantil termini, while AGB1 and AGG3 suppress ectopic epidermal rippling. These G protein mutant traits are recapitulated in long-day flowering gpa1-3 ft-10 plants, demonstrating that cantils, spurs and ectopic rippling occur as a function of delayed phase transition, rather than as a function of photoperiod per se.

Genome-Wide Analysis of RNA Secondary Structure.

Single-stranded RNA molecules fold into extraordinarily complicated secondary and tertiary structures as a result of intramolecular base pairing. In vivo, these RNA structures are not static. Instead, they are remodeled in response to changes in the prevailing physicochemical environment of the cell and as a result of intermolecular base pairing and interactions with RNA-binding proteins. Remarkable technical advances now allow us to probe RNA secondary structure at single-nucleotide resolution and genome-wide, both in vitro and in vivo. These data sets provide new glimpses into the RNA universe. Analyses of RNA structuromes in HIV, yeast, Arabidopsis, and mammalian cells and tissues have revealed regulatory effects of RNA structure on messenger RNA (mRNA) polyadenylation, splicing, translation, and turnover. Application of new methods for genome-wide identification of mRNA modifications, particularly methylation and pseudouridylation, has shown that the RNA "epitranscriptome" both influences and is influenced by RNA structure. In this review, we describe newly developed genome-wide RNA structure-probing methods and synthesize the information emerging from their application.

RNA multimerization as an organizing force for liquid-liquid phase separation.

RNA interactions are exceptionally strong and highly redundant. As such, nearly any two RNAs have the potential to interact with one another over relatively short stretches, especially at high RNA concentrations. This is especially true for pairs of RNAs that do not form strong self-structure. Such phenomena can drive liquid-liquid phase separation, either solely from RNA-RNA interactions in the presence of divalent or organic cations, or in concert with proteins. RNA interactions can drive multimerization of RNA strands via both base-pairing and tertiary interactions. In this article, we explore the tendency of RNA to form stable monomers, dimers, and higher order structures as a function of RNA length and sequence through a focus on the intrinsic thermodynamic, kinetic, and structural properties of RNA. The principles we discuss are independent of any specific type of biomolecular condensate, and thus widely applicable. We also speculate how external conditions experienced by living organisms can influence the formation of nonmembranous compartments, again focusing on the physical and structural properties of RNA. Plants, in particular, are subject to diverse abiotic stresses including extreme temperatures, drought, and salinity. These stresses and the cellular responses to them, including changes in the concentrations of small molecules such as polyamines, salts, and compatible solutes, have the potential to regulate condensate formation by melting or strengthening base-pairing. Reversible condensate formation, perhaps including regulation by circadian rhythms, could impact biological processes in plants, and other organisms.

Two novel GPCR-type G proteins are abscisic acid receptors in Arabidopsis.

In plants, G proteins modulate signaling by the stress hormone, abscisic acid (ABA). We identify and characterize two novel Arabidopsis proteins that show homology to an orphan vertebrate GPCR (GPR89) and interact with the sole Arabidopsis G protein alpha subunit, GPA1, but also have intrinsic GTP-binding and GTPase activity. We have named these proteins GPCR-type G proteins (GTG1 and GTG2). Arabidopsis mutants lacking both GTG1 and GTG2 exhibit ABA hyposensitivity. GTG1 and GTG2 bind ABA specifically. The GDP-bound form of the GTGs exhibits greater ABA binding than the GTP-bound form, the GTPase activity of the GTGs is inhibited by GPA1, and gpa1 null mutants exhibit ABA-hypersensitive phenotypes. These results predict that, unusually, it is the GDP-bound, not the GTP-bound, form of the GTGs that actively relays the signal. We propose that GTG proteins function both as a new type of G protein and as a class of membrane-localized ABA receptors.

Tissue-specific changes in the RNA structurome mediate salinity response in Arabidopsis.

Little is known concerning the effects of abiotic factors on in vivo RNA structures. We applied Structure-seq to assess the in vivo mRNA structuromes of Arabidopsis thaliana under salinity stress, which negatively impacts agriculture. Structure-seq utilizes dimethyl sulfate reactivity to identify As and Cs that lack base-pairing or protection. Salt stress refolded transcripts differentially in root versus shoot, evincing tissue specificity of the structurome. Both tissues exhibited an inverse correlation between salt stress-induced changes in transcript reactivity and changes in abundance, with stress-related mRNAs showing particular structural dynamism. This inverse correlation is more pronounced in mRNAs wherein the mean reactivity of the 5'UTR, CDS, and 3'UTR concertedly change under salinity stress, suggesting increased susceptibility to abundance control mechanisms in transcripts exhibiting this phenomenon, which we name "concordancy." Concordant salinity-induced increases in reactivity were notably observed in photosynthesis genes, thereby implicating mRNA structural loss in the well-known depression of photosynthesis by salt stress. Overall, changes in secondary structure appear to impact mRNA abundance, molding the functional specificity of the transcriptome under stress.

Structure-seq2 probing of RNA structure upon amino acid starvation reveals both known and novel RNA switches in Bacillus subtilis.

RNA structure influences numerous processes in all organisms. In bacteria, these processes include transcription termination and attenuation, small RNA and protein binding, translation initiation, and mRNA stability, and can be regulated via metabolite availability and other stresses. Here we use Structure-seq2 to probe the in vivo RNA structurome of Bacillus subtilis grown in the presence and absence of amino acids. Our results reveal that amino acid starvation results in lower overall dimethyl sulfate (DMS) reactivity of the transcriptome, indicating enhanced protection owing to protein binding or RNA structure. Starvation-induced changes in DMS reactivity correlated inversely with transcript abundance changes. This correlation was particularly pronounced in genes associated with the stringent response and CodY regulons, which are involved in adaptation to nutritional stress, suggesting that RNA structure contributes to transcript abundance change in regulons involved in amino acid metabolism. Structure-seq2 accurately reported on four known amino acid-responsive riboswitches: T-box, SAM, glycine, and lysine riboswitches. Additionally, we discovered a transcription attenuation mechanism that reduces yfmG expression when amino acids are added to the growth medium. We also found that translation of a leader peptide (YfmH) encoded just upstream of yfmG regulates yfmG expression. Our results are consistent with a model in which a slow rate of yfmH translation caused by limitation of the amino acids encoded in YfmH prevents transcription termination in the yfmG leader region by favoring formation of an overlapping antiterminator structure. This novel RNA switch offers a way to simultaneously monitor the levels of multiple amino acids.

GTP binding by Arabidopsis extra-large G protein 2 is not essential for its functions.

The extra-large guanosine-5'-triphosphate (GTP)-binding protein 2, XLG2, is an unconventional Gα subunit of the Arabidopsis (Arabidopsis thaliana) heterotrimeric GTP-binding protein complex with a major role in plant defense. In vitro biochemical analyses and molecular dynamic simulations show that affinity of XLG2 for GTP is two orders of magnitude lower than that of the conventional Gα, AtGPA1. Here we tested the physiological relevance of GTP binding by XLG2. We generated an XLG2(T476N) variant with abolished GTP binding, as confirmed by in vitro GTPγS binding assay. Yeast three-hybrid, bimolecular fluorescence complementation, and split firefly-luciferase complementation assays revealed that the nucleotide-depleted XLG2(T476N) retained wild-type XLG2-like interactions with the Gßγ dimer and defense-related receptor-like kinases. Both wild-type and nucleotide-depleted XLG2(T476N) restored the defense responses against Fusarium oxysporum and Pseudomonas syringae compromised in the xlg2 xlg3 double mutant. Additionally, XLG2(T476N) was fully functional restoring stomatal density, root growth, and sensitivity to NaCl, but failed to complement impaired germination and vernalization-induced flowering. We conclude that XLG2 is able to function in a GTP-independent manner and discuss its possible mechanisms of action.

Metabolomics of red-light-induced stomatal opening in Arabidopsis thaliana: Coupling with abscisic acid and jasmonic acid metabolism.

Environmental stimuli-triggered stomatal movement is a key physiological process that regulates CO2 uptake and water loss in plants. Stomata are defined by pairs of guard cells that perceive and transduce external signals, leading to cellular volume changes and consequent stomatal aperture change. Within the visible light spectrum, red light induces stomatal opening in intact leaves. However, there has been debate regarding the extent to which red-light-induced stomatal opening arises from direct guard cell sensing of red light versus indirect responses as a result of red light influences on mesophyll photosynthesis. Here we identify conditions that result in red-light-stimulated stomatal opening in isolated epidermal peels and enlargement of protoplasts, firmly establishing a direct guard cell response to red light. We then employ metabolomics workflows utilizing gas chromatography mass spectrometry and liquid chromatography mass spectrometry for metabolite profiling and identification of Arabidopsis guard cell metabolic signatures in response to red light in the absence of the mesophyll. We quantified 223 metabolites in Arabidopsis guard cells, with 104 found to be red light responsive. These red-light-modulated metabolites participate in the tricarboxylic acid cycle, carbon balance, phytohormone biosynthesis and redox homeostasis. We next analyzed selected Arabidopsis mutants, and discovered that stomatal opening response to red light is correlated with a decrease in guard cell abscisic acid content and an increase in jasmonic acid content. The red-light-modulated guard cell metabolome reported here provides fundamental information concerning autonomous red light signaling pathways in guard cells.

In vivo RNA structural probing of uracil and guanine base-pairing by 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide (EDC).

Many biological functions performed by RNAs arise from their in vivo structures. The structure of the same RNA can differ in vitro and in vivo owing in part to the influence of molecules ranging from protons to secondary metabolites to proteins. Chemical reagents that modify the Watson-Crick (WC) face of unprotected RNA bases report on the absence of base-pairing and so are of value to determining structures adopted by RNAs. Reagents have thus been sought that can report on the native RNA structures that prevail in living cells. Dimethyl sulfate (DMS) and glyoxal penetrate cell membranes and inform on RNA secondary structure in vivo through modification of adenine (A), cytosine (C), and guanine (G) bases. Uracil (U) bases, however, have thus far eluded characterization in vivo. Herein, we show that the water-soluble carbodiimide 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide (EDC) is capable of modifying the WC face of U and G in vivo, favoring the former nucleobase by a factor of ∼1.5, and doing so in the eukaryote rice, as well as in the Gram-negative bacterium Escherichia coli While both EDC and glyoxal target Gs, EDC reacts with Gs in their typical neutral state, while glyoxal requires Gs to populate the rare anionic state. EDC may thus be more generally useful; however, comparison of the reactivity of EDC and glyoxal may allow the identification of Gs with perturbed pKas in vivo and genome-wide. Overall, use of EDC with DMS allows in vivo probing of the base-pairing status of all four RNA bases.

The α subunit of the heterotrimeric G protein regulates mesophyll CO2 conductance and drought tolerance in rice.

Mesophyll conductance gm determines CO2 diffusion rates from mesophyll intercellular air spaces to the chloroplasts and is an important factor limiting photosynthesis. Increasing gm in cultivated plants is a potential strategy to increase photosynthesis and intrinsic water use efficiency (WUEi ). The anatomy of the leaf and metabolic factors such as aquaporins and carbonic anhydrases have been identified as important determinants of gm . However, genes involved in the regulation and modulation of gm remain largely unknown. In this work, we investigated the role of heterotrimeric G proteins in gm and drought tolerance in rice d1 mutants, which harbor a null mutation in the Gα subunit gene, RGA1. d1 mutants in both cv Nipponbare and cv Taichung 65 exhibited increased gm , fostering improvement in photosynthesis, WUEi , and drought tolerance compared with wild-type. The increased surface area of mesophyll cells and chloroplasts exposed to intercellular airspaces and the reduced cell wall and chloroplast thickness in the d1 mutant are evident contributors to the increase in gm . Our results indicate that manipulation of heterotrimeric G protein signaling has the potential to improve crop WUEi and productivity under drought.

Evolution and Structural Characteristics of Plant Voltage-Gated K+ Channels.

Plant voltage-gated K+ channels have been referred to as "plant Shakers" in reference to animal Shaker channels, the first K+ channels identified. Recent advances in our knowledge of K+ channel evolution and structure have significantly deepened the divide between these plant and animal K+ channels, suggesting that it is time to completely retire the "plant Shaker" designation. Evolutionary genomics reveals that plant voltage-gated K+ channels and metazoan Shakers derive from distinct prokaryotic ancestors. The plant channels belong to a lineage that includes cyclic nucleotide-gated channels and metazoan ether-à-go-go and hyperpolarization-activated, cyclic nucleotide-gated channels. We refer to this lineage as the CNBD channel superfamily, because all these channels share a cytoplasmic gating domain homologous to cyclic nucleotide binding domains. The first structures of CNBD superfamily channels reveal marked differences in coupling between the voltage sensor and ion-conducting pore relative to metazoan Shaker channels. Viewing plant voltage-gated K+ channel function through the lens of CNBD superfamily structures should lead to insights into how these channels are regulated.

Genome-wide RNA structurome reprogramming by acute heat shock globally regulates mRNA abundance.

The heat shock response is crucial for organism survival in natural environments. RNA structure is known to influence numerous processes related to gene expression, but there have been few studies on the global RNA structurome as it prevails in vivo. Moreover, how heat shock rapidly affects RNA structure genome-wide in living systems remains unknown. We report here in vivo heat-regulated RNA structuromes. We applied Structure-seq chemical [dimethyl sulfate (DMS)] structure probing to rice (Oryza sativa L.) seedlings with and without 10 min of 42 °C heat shock and obtained structural data on >14,000 mRNAs. We show that RNA secondary structure broadly regulates gene expression in response to heat shock in this essential crop species. Our results indicate significant heat-induced elevation of DMS reactivity in the global transcriptome, revealing RNA unfolding over this biological temperature range. Our parallel Ribo-seq analysis provides no evidence for a correlation between RNA unfolding and heat-induced changes in translation, in contrast to the paradigm established in prokaryotes, wherein melting of RNA thermometers promotes translation. Instead, we find that heat-induced DMS reactivity increases correlate with significant decreases in transcript abundance, as quantified from an RNA-seq time course, indicating that mRNA unfolding promotes transcript degradation. The mechanistic basis for this outcome appears to be mRNA unfolding at both 5' and 3'-UTRs that facilitates access to the RNA degradation machinery. Our results thus reveal unexpected paradigms governing RNA structural changes and the eukaryotic RNA life cycle.

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

Cytosolic calcium concentration ([Ca2+ ]cyt ) and heterotrimeric G-proteins are universal eukaryotic signaling elements. In plant guard cells, extracellular calcium (Cao ) 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 Cao responses: induction of stomatal closure; inhibition of stomatal opening; [Ca2+ ]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 Cao . 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 Cao . 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 [Ca2+ ]cyt oscillations in response to Cao , initiating only a single [Ca2+ ]cyt spike. Experimentally imposed [Ca2+ ]cyt oscillations restored stomatal closure in agb1. Yeast two-hybrid and bimolecular complementation fluorescence experiments revealed that AGB1 interacts with phospholipase Cs (PLCs), and Cao induced InsP3 production in Col but not in agb1. In sum, G-protein signaling via AGB1/AGG1/AGG2 is essential for Cao -regulation of stomatal apertures, and stomatal movements in response to Cao apparently require Ca2+ -induced Ca2+ release that is likely dependent on Gßγ interaction with PLCs leading to InsP3 production.

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

ABSTARCT: KEY MESSAGE: The timing and transcriptomic changes during the C3 to CAM transition of common ice plant support the notion that guard cells themselves can shift from C3 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 CO2 uptake. Mesembryanthemum crystallinum (common ice plant) is a facultative CAM species that can shift from C3 photosynthesis to CAM under salt or drought stresses. However, the molecular mechanisms underlying the stress induced transition from C3 to CAM remain unknown. Here we determined the transition time from C3 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 C3 to CAM. This study provides important information towards introducing CAM stomatal behavior into C3 crops to enhance water use efficiency.

Glyoxals as in vivo RNA structural probes of guanine base-pairing.

Elucidation of the folded structures that RNA forms in vivo is vital to understanding its functions. Chemical reagents that modify the Watson-Crick (WC) face of unprotected nucleobases are particularly useful in structure elucidation. Dimethyl sulfate penetrates cell membranes and informs on RNA base-pairing and secondary structure but only modifies the WC face of adenines and cytosines. We present glyoxal, methylglyoxal, and phenylglyoxal as potent in vivo reagents that target the WC face of guanines as well as cytosines and adenines. Tests on rice (Oryza sativa) 5.8S rRNA in vitro read out by reverse transcription and gel electrophoresis demonstrate specific modification of almost all guanines in a time- and pH-dependent manner. Subsequent in vivo tests on rice, a eukaryote, and Bacillus subtilis and Escherichia coli, Gram-positive and Gram-negative bacteria, respectively, showed that all three reagents enter living cells without prior membrane permeabilization or pH adjustment of the surrounding media and specifically modify solvent-exposed guanine, cytosine, and adenine residues.

A new discrete dynamic model of ABA-induced stomatal closure predicts key feedback loops.

Stomata, microscopic pores in leaf surfaces through which water loss and carbon dioxide uptake occur, are closed in response to drought by the phytohormone abscisic acid (ABA). This process is vital for drought tolerance and has been the topic of extensive experimental investigation in the last decades. Although a core signaling chain has been elucidated consisting of ABA binding to receptors, which alleviates negative regulation by protein phosphatases 2C (PP2Cs) of the protein kinase OPEN STOMATA 1 (OST1) and ultimately results in activation of anion channels, osmotic water loss, and stomatal closure, over 70 additional components have been identified, yet their relationships with each other and the core components are poorly elucidated. We integrated and processed hundreds of disparate observations regarding ABA signal transduction responses underlying stomatal closure into a network of 84 nodes and 156 edges and, as a result, established those relationships, including identification of a 36-node, strongly connected (feedback-rich) component as well as its in- and out-components. The network's domination by a feedback-rich component may reflect a general feature of rapid signaling events. We developed a discrete dynamic model of this network and elucidated the effects of ABA plus knockout or constitutive activity of 79 nodes on both the outcome of the system (closure) and the status of all internal nodes. The model, with more than 1024 system states, is far from fully determined by the available data, yet model results agree with existing experiments in 82 cases and disagree in only 17 cases, a validation rate of 75%. Our results reveal nodes that could be engineered to impact stomatal closure in a controlled fashion and also provide over 140 novel predictions for which experimental data are currently lacking. Noting the paucity of wet-bench data regarding combinatorial effects of ABA and internal node activation, we experimentally confirmed several predictions of the model with regard to reactive oxygen species, cytosolic Ca2+ (Ca2+c), and heterotrimeric G-protein signaling. We analyzed dynamics-determining positive and negative feedback loops, thereby elucidating the attractor (dynamic behavior) repertoire of the system and the groups of nodes that determine each attractor. Based on this analysis, we predict the likely presence of a previously unrecognized feedback mechanism dependent on Ca2+c. This mechanism would provide model agreement with 10 additional experimental observations, for a validation rate of 85%. Our research underscores the importance of feedback regulation in generating robust and adaptable biological responses. The high validation rate of our model illustrates the advantages of discrete dynamic modeling for complex, nonlinear systems common in biology.