Tandem stretches of guanines can associate in hydrogen-bonded arrays to form G-quadruplexes, which are stabilized by K+ ions. Using computational methods, we searched for G-Quadruplex Sequence (GQS) patterns in the model plant species Arabidopsis thaliana. We found ∼1200 GQS with a G3 repeat sequence motif, most of which are located in the intergenic region. Using a Markov modeled genome, we determined that GQS are significantly underrepresented in the genome. Additionally, we found ∼43 000 GQS with a G2 repeat sequence motif; notably, 80% of these were located in genic regions, suggesting that these sequences may fold at the RNA level. Gene Ontology functional analysis revealed that GQS are overrepresented in genes encoding proteins of certain functional categories, including enzyme activity. Conversely, GQS are underrepresented in other categories of genes, notably those for non-coding RNAs such as tRNAs and rRNAs. We also find that genes that are differentially regulated by drought are significantly more likely to contain a GQS. CD-detected K+ titrations performed on representative RNAs verified formation of quadruplexes at physiological K+ concentrations. Overall, this study indicates that GQS are present at unique locations in Arabidopsis and that folding of RNA GQS may play important roles in regulating gene expression.

Changes in RNA conformation can alter gene expression. The guanine quadruplex sequence (GQS) is an RNA motif that folds in the presence of K(+) ions. Changes in the conformation of this motif could be especially important in regulating gene expression in plants because intracellular K(+) concentrations often increase during drought stress. Little is known about the folding thermodynamics of RNA GQS. We show here that RNA GQS with tracts containing three G's [e.g., (GGGxx)(4)] have a modest dependence on the K(+) concentration, folding with no or even negative cooperativity (Hill coefficients ≤1), and are associated with populated folding intermediates. In contrast, GQS with tracts containing just two G's [e.g., (GGxx)(4)] have a steep dependence on the K(+) concentration and fold with positive cooperativity (Hill coefficients of 1.7-2.7) without significantly populating intermediate states. We postulate that in plants, the more stable G3 sequences are largely folded even under unstressed conditions, while the less stable G2 sequences fold only at the higher K(+) concentrations associated with cellular stress, wherein they respond sharply to changing K(+) concentrations. Given the binary nature of their folding, G2 sequences may find application in computation with DNA and in engineering of genetic circuits.

No abstract