Cis-regulatory element based targeted gene finding: genome-wide identification of abscisic acid- and abiotic stress-responsive genes in Arabidopsis thaliana

The harvest-responsive region of the Asparagus officinalis sparagine synthetase promoter reveals complexity in the regulation of the harvest response

The activity of a 1915-bp asparagine synthetase (AS) promoter of Asparagus officinalis L. was induced in mature leaves of transgenic Arabidopsis thaliana (L.) Heynh. plants when the leaves were detached and held in water for 24 h. To understand this induction by harvest, variants of the AS promoter were linked to the β-glucuronidase GUS reporter gene. Harvest induction in the leaves required detachment and was not simply a wound response. Two regions in the AS promoter (Region A, -640 to -523; Region B, -524 to -383) were independently able to confer harvest response to the otherwise unresponsive -383AS (minimal) promoter. Region A was studied in further detail. Various truncations, deletions, or nucleotide substitutions of Region A affected activity and fold induction of the minimal promoter. However, no harvest-inducible cis-acting element within Region A was identified. Although the minimal promoter contained a dehydration-responsive element and ACGT elements similar to ABA-responsive regulatory motifs these were not needed by the upstream regulatory regions for directing harvest response. When four copies of Region A were linked to the minimal promoter it became highly active in leaves before harvest. Deletions within Region A showed that it required its complete 117 bp for driving harvest response, yet the region cannot simply be thought of as a harvest-responsive module, since its concatemerisation led to constitutive expression.

Arabidopsis thaliana regulatory element analyzer

In the Arabidopsis thaliana regulatory element analyzer (AtREA) server, we have integrated sequence data, genome-wide expression data and functional annotation data in three application modules which will be useful to identify major regulatory targets of a user-provided cis-regulatory element (CRE), study different features of CRE distribution and evaluate the role of a set of CREs in the regulation of gene expression--independently as well as in combination with other user-provided CREs.

AVAILABILITY

AtREA is freely available at http://www.bioinformatics.org/grn/atrea.html.

Genome-wide identification of osmotic stress response gene in Arabidopsis thaliana

In this paper, we present a cis-regulatory element based computational approach to genome-wide identification of genes putatively responding to various osmotic stresses in Arabidopsis thaliana. The rationale of our method is that gene expression is largely controlled at the transcriptional level through the interactions between transcription factors and cis-regulatory elements. Using cis-regulatory motifs known to regulate osmotic stress response, we therefore built an artificial neural network model to identify other functionally relevant genes involved in the same process. We performed Gene Ontology enrichment analysis on the 500 top-scoring predictions and found that, except for un-annotated ORFs ( approximately 40%), 91.3% of the enriched GO classification was related to stress response and ABA response. Publicly available gene expression profiling data of Arabidopsis under various stresses were used for cross validation. We also conducted RT-PCR analysis to experimentally verify selected predictions. According to our results, transcript levels of 27 out of 41 top-ranked genes (65.8%) altered under various osmotic stress treatments. We believe that a similar approach could be extensively adopted elsewhere to infer gene function in various cellular processes from different species.

Two ABREs, two redundant root-specific and one W-box cis-acting elements are functional in the sunflower HAHB4 promoter

HAHB4 is a sunflower gene encoding a homeodomain-leucine zipper (HD-Zip) transcription factor. It was previously demonstrated that this gene is regulated at the transcriptional level by several abiotic factors and hormones. A previous analysis in the PLACE database revealed the presence of four putative ABREs. In this work these four elements and also one W-box and two root-specific expression elements were characterized as functional. Site-directed mutagenesis on the promoter, stable transformation of Arabidopis plants as well as transient transformation of sunflower leaves, were performed. The analysis of the transformants was carried out by histochemistry and real time RT-PCR. The results indicate that just one ABRE out of the four is responsible for ABA, NaCl and drought regulation. However, NaCl induction occurs also by an additional ABA-independent way involving another two overlapped ABREs. On the other hand, it was determined that the W-box located 5' upstream is responsive to ethylene and only two root-specific expression elements, among the several detected, are functional but redundant. Conservation of molecular mechanisms between sunflower and Arabidopsis is strongly supported by this experimental work.

STIF: Identification of stress-upregulated transcription factor binding sites in Arabidopsis thaliana

The expressions of proteins in the cell are carefully regulated by transcription factors that interact with their downstream targets in specific signal transduction cascades. Our understanding of the regulation of functional genes responsive to stress signals is still nascent. Plants like Arabidopsis thaliana, are convenient model systems to study fundamental questions related to regulation of the stress transcriptome in response to stress challenges. Microarray results of the Arabidopsis transcriptome indicate that several genes could be upregulated during multiple stresses, such as cold, salinity, drought etc. Experimental biochemical validations have proved the involvement of several transcription factors could be involved in the upregulation of these stress responsive genes. In order to follow the intricate and complicated networks of transcription factors and genes that respond to stress situations in plants, we have developed a computer algorithm that can identify key transcription factor binding sites upstream of a gene of interest. Hidden Markov models of the transcription factor binding sites enable the identification of predicted sites upstream of plant stress genes. The search algorithm, STIF, performs very well, with more than 90% sensitivity, when tested on experimentally validated positions of transcription factor binding sites on a dataset of 60 stress upregulated genes.

AVAILABILITY

Supplementary data is available at http://caps.ncbs.res.in/download/stif.

A third phytoene synthase is devoted to abiotic stress-induced abscisic acid formation in rice and defines functional diversification of phytoene synthase genes

We here report on the characterization of a novel third phytoene synthase gene (PSY) in rice (Oryza sativa), OsPSY3, and on the differences among all three PSY genes with respect to the tissue-specific expression and regulation upon various environmental stimuli. The two already known PSYs are under phytochrome control and involved in carotenoid biosynthesis in photosynthetically active tissues and exhibit different expression patterns during chloroplast development. In contrast, OsPSY3 transcript levels are not affected by light and show almost no tissue-specific differences. Rather, OsPSY3 transcripts are up-regulated during increased abscisic acid (ABA) formation upon salt treatment and drought, especially in roots. The simultaneous induction of genes encoding 9-cis-epoxycarotenoid dioxygenases (NCEDs), involved in the initial steps of ABA biosynthesis, indicate that decreased xanthophyll levels are compensated by the induction of the third PSY gene. Furthermore, OsPSY3 and the OsNCEDs investigated were also induced by the application of ABA, indicating positive feedback regulation. The regulatory differences are mirrored by cis-acting elements in the corresponding promoter regions, with light-responsive elements for OsPSY1 and OsPSY2 and an ABA-response element as well as a coupling element for OsPSY3. The investigation of the gene structures and 5' untranslated regions revealed that OsPSY1 represents a descendant of an ancient PSY gene present in the common ancestor of monocots and dicots. Since the genomic structures of OsPSY2 and OsPSY3 are comparable, we conclude that they originated from the most recent common ancestor, OsPSY1.

Microarray-based analysis of stress-regulated microRNAs in Arabidopsis thaliana

High-salinity, drought, and low temperature are three common environmental stress factors that seriously influence plant growth and development worldwide. Recently, microRNAs (miRNAs) have emerged as a class of gene expression regulators that have also been linked to stress responses. However, the relationship between miRNA expression and stress responses is just beginning to be explored. Here, we identified 14 stress-inducible miRNAs using microarray data in which the effects of three abiotic stresses were surveyed in Arabidopsis thaliana. Among them, 10 high-salinity-, four drought-, and 10 cold-regulated miRNAs were detected, respectively. miR168, miR171, and miR396 responded to all of the stresses. Expression profiling by RT-PCR analysis showed great cross-talk among the high-salinity, drought, and cold stress signaling pathways. The existence of stress-related elements in miRNA promoter regions provided further evidence supporting our results. These findings extend the current view about miRNA as ubiquitous regulators under stress conditions.

N-Acylethanolamine metabolism interacts with abscisic acid signaling in Arabidopsis thaliana seedlings

N-Acylethanolamines (NAEs) are bioactive acylamides that are present in a wide range of organisms. In plants, NAEs are generally elevated in desiccated seeds, suggesting that they may play a role in seed physiology. NAE and abscisic acid (ABA) levels were depleted during seed germination, and both metabolites inhibited the growth of Arabidopsis thaliana seedlings within a similar developmental window. Combined application of low levels of ABA and NAE produced a more dramatic reduction in germination and growth than either compound alone. Transcript profiling and gene expression studies in NAE-treated seedlings revealed elevated transcripts for a number of ABA-responsive genes and genes typically enriched in desiccated seeds. The levels of ABI3 transcripts were inversely associated with NAE-modulated growth. Overexpression of the Arabidopsis NAE degrading enzyme fatty acid amide hydrolase resulted in seedlings that were hypersensitive to ABA, whereas the ABA-insensitive mutants, abi1-1, abi2-1, and abi3-1, exhibited reduced sensitivity to NAE. Collectively, our data indicate that an intact ABA signaling pathway is required for NAE action and that NAE may intersect the ABA pathway downstream from ABA. We propose that NAE metabolism interacts with ABA in the negative regulation of seedling development and that normal seedling establishment depends on the reduction of the endogenous levels of both metabolites.

Genome-wide analysis of ABA-responsive elements ABRE and CE3 reveals divergent patterns in Arabidopsis and rice

Background

In plants, complex regulatory mechanisms are at the core of physiological and developmental processes. The phytohormone abscisic acid (ABA) is involved in the regulation of various such processes, including stomatal closure, seed and bud dormancy, and physiological responses to cold, drought and salinity stress. The underlying tissue or plant-wide control circuits often include combinatorial gene regulatory mechanisms and networks that we are only beginning to unravel with the help of new molecular tools. The increasing availability of genomic sequences and gene expression data enables us to dissect ABA regulatory mechanisms at the individual gene expression level. In this paper we used an in-silico-based approach directed towards genome-wide prediction and identification of specific features of ABA-responsive elements. In particular we analysed the genome-wide occurrence and positional arrangements of two well-described ABA-responsive cis-regulatory elements (CREs), ABRE and CE3, in thale cress (Arabidopsis thaliana) and rice (Oryza sativa).

Results

Our results show that Arabidopsis and rice use the ABA-responsive elements ABRE and CE3 distinctively. Earlier reports for various monocots have identified CE3 as a coupling element (CE) associated with ABRE. Surprisingly, we found that while ABRE is equally abundant in both species, CE3 is practically absent in Arabidopsis. ABRE-ABRE pairs are common in both genomes, suggesting that these can form functional ABA-responsive complexes (ABRCs) in Arabidopsis and rice. Furthermore, we detected distinct combinations, orientation patterns and DNA strand preferences of ABRE and CE3 motifs in rice gene promoters.

Conclusion

Our computational analyses revealed distinct recruitment patterns of ABA-responsive CREs in upstream sequences of Arabidopsis and rice. The apparent absence of CE3s in Arabidopsis suggests that another CE pairs with ABRE to establish a functional ABRC capable of interacting with transcription factors. Further studies will be needed to test whether the observed differences are extrapolatable to monocots and dicots in general, and to understand how they contribute to the fine-tuning of the hormonal response. The outcome of our investigation can now be used to direct future experimentation designed to further dissect the ABA-dependent regulatory networks.

Transcriptomic adaptations in rice suspension cells under sucrose starvation

Sugar is an important resource for energy generation and developmental regulation in plants, and sucrose starvation causes enormous changes in cellular morphology, enzyme activities and gene expression. Genome-wide gene expression profiling provides a comprehensive knowledge of gene expression under nutrient depletion and senescence; however, that of a monocot model plant, rice, under sucrose depletion is still under investigation. Here, the time-course monitoring of gene expression profiles in sucrose-starved rice (Oryza sativa cv Tainung67) suspension cells was investigated by 21495 probes contained in Agilent rice chip. In sucrose-starved cells, the induced vacuolar biogenesis coincided with significantly upregulated transcripts of H+-pyrophosphatase, delta-TIP, one putative alpha-TIP, several vacuolar proteases and proteinase inhibitors, and one OsATG3. To survey the overall metabolic adaptations under sucrose depletion, the genes with significantly altered expression level were incorporated into multiple metabolic pathways. Most genes encoding enzymes involved in biosynthesis and degradation pathways of various macromolecules were comprehensively down-and upregulated, respectively, with sucrose starvation. Transcriptional regulation of gene expression is important for physiological adaptations to environmental stress, and many transcription factors, including bZIPs, NACs, and WRKY, showed significant increase in transcriptional level under sucrose starvation. Concurrently, statistical analysis revealed that their corresponding consensus cis-elements, such as ABA-responsive element, CACG, ACI, ACII and CTTATCC, were frequently found in the promoter regions of many sucrose starvation-upregulated genes. Particle bombardment-mediated and luciferase activity-based transient promoter assays revealed the CTTATCC, derived form TATCCA, and the AC motifs to be promising sucrose-starvation responsive activators in rice suspension cells.

Microarray analysis of transcriptional responses to abscisic acid and osmotic, salt, and drought stress in the moss, Physcomitrella patens

Dehydration tolerance was an adaptive trait necessary for the colonization of land by plants, and remains widespread among bryophytes: the nearest extant relatives of the first land plants. A genome-wide analysis was undertaken of water-stress responses in the model moss Physcomitrella patens to identify stress-responsive genes. An oligonucleotide microarray was used for transcriptomic analysis of Physcomitrella treated with abscisic acid (ABA), or subjected to osmotic, salt and drought stress. Bioinformatic analysis of the Physcomitrella genome identified the responsive genes, and a number of putative stress-related cis-regulatory elements. In protonemal tissue, 130 genes were induced by dehydration, 56 genes by ABA, but only 10 and eight genes, respectively, by osmotic and salt stress. Fifty-one genes were induced by more than one treatment. Seventy-six genes, principally encoding chloroplast proteins, were drought down-regulated. Many ABA- and drought-responsive genes are homologues of angiosperm genes expressed during drought stress and seed development. These ABA- and drought-responsive genes include those encoding a number of late embryogenesis abundant (LEA) proteins, a 'DREB' transcription factor and a Snf-related kinase homologous with the Arabidopsis ABA signal transduction component 'OPEN STOMATA 1'. Evolutionary capture of conserved stress-regulatory transcription factors by the seed developmental pathway probably accounts for the seed-specificity of desiccation tolerance among angiosperms.

Cis-motifs upstream of the transcription and translation initiation sites are effectively revealed by their positional disequilibrium in eukaryote genomes using frequency distribution curves

Background

The discovery of cis-regulatory motifs still remains a challenging task even though the number of sequenced genomes is constantly growing. Computational analyses using pattern search algorithms have been valuable in phylogenetic footprinting approaches as have expression profile experiments to predict co-occurring motifs. Surprisingly little is known about the nature of cis-regulatory element (CRE) distribution in promoters.

Results

In this paper we used the Motif Mapper open-source collection of visual basic scripts for the analysis of motifs in any aligned set of DNA sequences. We focused on promoter motif distribution curves to identify positional over-representation of DNA motifs. Using differentially aligned datasets from the model species Arabidopsis thaliana, Caenorhabditis elegans, Drosophila melanogaster and Saccharomyces cerevisiae, we convincingly demonstrated the importance of the position and orientation for motif discovery. Analysis with known CREs and all possible hexanucleotides showed that some functional elements gather close to the transcription and translation initiation sites and that elements other than the TATA-box motif are conserved between eukaryote promoters. While a high background frequency usually decreases the effectiveness of such an enumerative investigation, we improved our analysis by conducting motif distribution maps using large datasets.

Conclusion

This is the first study to reveal positional over-representation of CREs and promoter motifs in a cross-species approach. CREs and motifs shared between eukaryotic promoters support the observation that an eukaryotic promoter structure has been conserved throughout evolutionary time. Furthermore, with the information on positional enrichment of a motif or a known functional CRE, it is possible to get a more detailed insight into where an element appears to function. This in turn might accelerate the in depth examination of known and yet unknown cis-regulatory sequences in the laboratory.

ATTED-II: a database of co-expressed genes and cis elements for identifying co-regulated gene groups in Arabidopsis

Publicly available database of co-expressed gene sets would be a valuable tool for a wide variety of experimental designs, including targeting of genes for functional identification or for regulatory investigation. Here, we report the construction of an Arabidopsis thaliana trans-factor and cis-element prediction database (ATTED-II) that provides co-regulated gene relationships based on co-expressed genes deduced from microarray data and the predicted cis elements. ATTED-II () includes the following features: (i) lists and networks of co-expressed genes calculated from 58 publicly available experimental series, which are composed of 1388 GeneChip data in A.thaliana; (ii) prediction of cis-regulatory elements in the 200 bp region upstream of the transcription start site to predict co-regulated genes amongst the co-expressed genes; and (iii) visual representation of expression patterns for individual genes. ATTED-II can thus help researchers to clarify the function and regulation of particular genes and gene networks.

PlantQTL-GE: a database system for identifying candidate genes in rice and Arabidopsis by gene expression and QTL information

We have designed and implemented a web-based database system, called PlantQTL-GE, to facilitate quantitatine traits locus (QTL) based candidate gene identification and gene function analysis. We collected a large number of genes, gene expression information in microarray data and expressed sequence tags (ESTs) and genetic markers from multiple sources of Oryza sativa and Arabidopsis thaliana. The system integrates these diverse data sources and has a uniform web interface for easy access. It supports QTL queries specifying QTL marker intervals or genomic loci, and displays, on rice or Arabidopsis genome, known genes, microarray data, ESTs and candidate genes and similar putative genes in the other plant. Candidate genes in QTL intervals are further annotated based on matching ESTs, microarray gene expression data and cis-elements in regulatory sequences. The system is freely available at .

Rapid transcriptome changes induced by cytosolic Ca2+ transients reveal ABRE-related sequences as Ca2+-responsive cis elements in Arabidopsis

The regulation of gene expression by cellular calcium is crucial for plant defense against biotic and abiotic stresses. However, the number of genes known to respond to specific transient calcium signals is limited, and as yet there is no definition of a calcium-responsive cis element in plants. Here, we generated specific cytosolic calcium transients in intact Arabidopsis thaliana seedlings and linked them to early transcriptome changes, followed by bioinformatic analysis of the responsive genes. A cytosolic calcium transient induced by calmodulin antagonists and blocked by lanthanides was characterized using aequorin-based luminometry and photon imaging. Analysis of transcriptome changes revealed 230 calcium-responsive genes, of which 162 were upregulated and 68 were downregulated. These include known early stress-responsive genes as well as genes of unknown function. Analysis of their upstream regions revealed, exclusively in the upregulated genes, a highly significant occurrence of a consensus sequence (P < 10(-13)) comprising two abscisic acid-specific cis elements: the abscisic acid-responsive element (ABRE; CACGTG[T/C/G]) and its coupling element ([C/A]ACGCG[T/C/G]) [corrected] Finally, we show that a tetramer of the ABRE cis element is sufficient to confer transcriptional activation in response to cytosolic Ca(2+) transients. Thus, at least for some specific Ca(2+) transients and motif combinations, ABREs function as Ca(2+)-responsive cis elements.

Computational prediction and experimental verification of HVA1-like abscisic acid responsive promoters in rice (Oryza sativa)

Abscisic acid (ABA) is one of the central plant hormones, responsible for controlling both maturation and germination in seeds, as well as mediating adaptive responses to desiccation, injury, and pathogen infection in vegetative tissues. Thorough analyses of two barley genes, HVA1 and HVA22, indicate that their response to ABA relies on the interaction of two cis-acting elements in their promoters, an ABA response element (ABRE) and a coupling element (CE). Together, they form an ABA response promoter complex (ABRC). Comparison of promoters of barley HVA1 and it rice orthologue indicates that the structures and sequences of their ABRCs are highly similar. Prediction of ABA responsive genes in the rice genome is then tractable to a bioinformatics approach based on the structures of the well-defined barley ABRCs. Here we describe a model developed based on the consensus, inter-element spacing and orientations of experimentally determined ABREs and CEs. Our search of the rice promoter database for promoters that fit the model has generated a partial list of genes in rice that have a high likelihood of being involved in the ABA signaling network. The ABA inducibility of some of the rice genes identified was validated with quantitative reverse transcription PCR (QPCR). By limiting our input data to known enhancer modules and experimentally derived rules, we have generated a high confidence subset of ABA-regulated genes. The results suggest that the pathways by which cereals respond to biotic and abiotic stresses overlap significantly, and that regulation is not confined to the level transcription. The large fraction of putative regulatory genes carrying HVA1-like enhancer modules in their promoters suggests the ABA signal enters at multiple points into a complex regulatory network that remains largely unmapped.

Unraveling abiotic stress tolerance mechanisms--getting genomics going

Homeostasis, a set-value for metabolism under optimal conditions, is rarely achieved by plants because of the cost exerted by external stress factors: climatic, biotic, and nutrient imbalances. Among these, stresses caused by abiotic conditions, such as temperature extremes (freezing, cold and heat), water availability (drought and ion excess) and ion toxicity (salinity and heavy metals), have been difficult to dissect because defense responses to abiotic factors require regulatory changes to the activation of multiple genes and pathways. Genomics technologies that have emerged during the past decade have been useful in addressing, in an integrated fashion, the multigenicity of the plant abiotic stress response through genome sequences; cell-, organ-, tissue- and stress-specific transcript collections; transcript, protein and metabolite profiles and their dynamic changes; protein interactions; and mutant screens.