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

Four Arabidopsis AREB/ABF transcription factors function predominantly in gene expression downstream of SnRK2 kinases in abscisic acid signalling in response to osmotic stress

Under osmotic stress conditions such as drought and high salinity, the plant hormone abscisic acid (ABA) plays important roles in stress-responsive gene expression mainly through three bZIP transcription factors, AREB1/ABF2, AREB2/ABF4 and ABF3, which are activated by SNF1-related kinase 2s (SnRK2s) such as SRK2D/SnRK2.2, SRK2E/SnRK2.6 and SRK2I/SnRK2.3 (SRK2D/E/I). However, since the three AREB/ABFs are crucial, but not exclusive, for the SnRK2-mediated gene expression, transcriptional pathways governed by SRK2D/E/I are not fully understood. Here, we show that a bZIP transcription factor, ABF1, is a functional homolog of AREB1, AREB2 and ABF3 in ABA-dependent gene expression in Arabidopsis. Despite lower expression levels of ABF1 than those of the three AREB/ABFs, the areb1 areb2 abf3 abf1 mutant plants displayed increased sensitivity to drought and decreased sensitivity to ABA in primary root growth compared with the areb1 areb2 abf3 mutant. Genome-wide transcriptome analyses revealed that expression of downstream genes of SRK2D/E/I, which include many genes functioning in osmotic stress responses and tolerance such as transcription factors and LEA proteins, was mostly impaired in the quadruple mutant. Thus, these results indicate that the four AREB/ABFs are the predominant transcription factors downstream of SRK2D/E/I in ABA signalling in response to osmotic stress during vegetative growth.

Identification of microRNAs in six solanaceous plants and their potential link with phosphate and mycorrhizal signaling

To date, only a limited number of solanaceous miRNAs have been deposited in the miRNA database. Here, genome-wide bioinformatic identification of miRNAs was performed in six solanaceous plants (potato, tomato, tobacco, eggplant, pepper, and petunia). A total of 2,239 miRNAs were identified following a range of criteria, of which 982 were from potato, 496 from tomato, 655 from tobacco, 46 from eggplant, 45 were from pepper, and 15 from petunia. The sizes of miRNA families and miRNA precursor length differ in all the species. Accordingly, 620 targets were predicted, which could be functionally classified as transcription factors, metabolic enzymes, RNA and protein processing proteins, and other proteins for plant growth and development. We also showed evidence for miRNA clusters and sense and antisense miRNAs. Additionally, five Pi starvation- and one arbuscular mycorrhiza (AM)-related cis-elements were found widely distributed in the putative promoter regions of the miRNA genes. Selected miRNAs were classified into three groups based on the presence or absence of P1BS and MYCS cis-elements, and their expression in response to Pi starvation and AM symbiosis was validated by quantitative reverse transcription-polymerase chain reaction (qRT-PCR). These results show that conserved miRNAs exist in solanaceous species and they might play pivotal roles in plant growth, development, and stress responses.

Natural variation in abiotic stress responsive gene expression and local adaptation to climate in Arabidopsis thaliana

Gene expression varies widely in natural populations, yet the proximate and ultimate causes of this variation are poorly known. Understanding how variation in gene expression affects abiotic stress tolerance, fitness, and adaptation is central to the field of evolutionary genetics. We tested the hypothesis that genes with natural genetic variation in their expression responses to abiotic stress are likely to be involved in local adaptation to climate in Arabidopsis thaliana. Specifically, we compared genes with consistent expression responses to environmental stress (expression stress responsive, "eSR") to genes with genetically variable responses to abiotic stress (expression genotype-by-environment interaction, "eGEI"). We found that on average genes that exhibited eGEI in response to drought or cold had greater polymorphism in promoter regions and stronger associations with climate than those of eSR genes or genomic controls. We also found that transcription factor binding sites known to respond to environmental stressors, especially abscisic acid responsive elements, showed significantly higher polymorphism in drought eGEI genes in comparison to eSR genes. By contrast, eSR genes tended to exhibit relatively greater pairwise haplotype sharing, lower promoter diversity, and fewer nonsynonymous polymorphisms, suggesting purifying selection or selective sweeps. Our results indicate that cis-regulatory evolution and genetic variation in stress responsive gene expression may be important mechanisms of local adaptation to climatic selective gradients.

Sequence and functional analyses of the aldehyde dehydrogenase 7B4 gene promoter in Arabidopsis thaliana and selected Brassicaceae: regulation patterns in response to wounding and osmotic stress

Aldehyde dehydrogenases metabolise a wide range of aliphatic and aromatic aldehydes, which become cytotoxic at high levels. Family 7 aldehyde dehydrogenase genes, often described as antiquitins or turgor-responsive genes in plants, are broadly conserved across all domains. Despite the high conservation of the plant ALDH7 proteins and their importance in stress responses, their regulation has not been investigated. Here, we compared ALDH7 genes of different Brassicaceae and found that, in contrast to the gene organisation and protein coding sequences, similarities in the promoter sequences were limited to the first few hundred nucleotides upstream of the translation start codon. The function of this region was studied by isolating the core promoter of the Arabidopsis thaliana ALDH7B4 gene, taken as model. The promoter was found to be responsive to wounding in addition to salt and dehydration stress. Cis-acting elements involved in stress responsiveness were analysed and two conserved ACGT-containing motifs proximal to the translation start codon were found to be essential for the responsiveness to osmotic stress in leaves and in seeds. The integrity of an upstream ACGT motif and a dehydration-responsive element/C-repeat-low temperature-responsive element was found to be necessary for ALDH7B4 expression in seeds and induction by salt, dehydration and ABA in leaves. The comparison of the gene expression in selected Arabidopsis mutants demonstrated that osmotic stress-induced ALDH7B4 expression in leaves and seeds involves both ABA- and lipid-signalling components.

Genome-wide analysis and identification of stress-responsive genes of the CCCH zinc finger family in Solanum lycopersicum

Zinc finger genes comprise a large and diverse gene family. Based on their individual finger structures and spacing, zinc finger proteins are further divided into different families according to their specific molecular functions. Genes in the CCCH family encode zinc finger proteins containing a motif with three cysteines and one histidine. They play important roles in plant growth and development, and in response to biotic and abiotic stresses. However, the limited analysis of the genome sequence has meant that there is no detailed information concerning the CCCH zinc finger family in tomato (Solanum lycopersicum). Here, we identified 80 CCCH zinc finger protein genes in the tomato genome. A complete overview of this gene family in tomato was presented, including the chromosome locations, gene duplications, phylogeny, gene structures and protein motifs. Promoter sequences and expression profiles of putative stress-responsive members were also investigated. These results revealed that, with the exception of four genes, the 80 CCCH genes are distributed over all 12 chromosomes with different densities, and include six segmental duplication events. The CCCH family in tomato could be divided into 12 groups based on their different CCCH motifs and into eight subfamilies by phylogenetic analysis. Analysis showed that almost all CCCH genes contain putative stress-responsive cis-elements in their promoter regions. Nine CCCH genes chosen for further quantitative real-time PCR analysis showed differential expression patterns in three representative tomato tissues. In addition, their expression levels indicated that these genes are mostly involved in the response to mannitol, heat, salicylic acid, ethylene or methyl jasmonate treatments. To the best of our knowledge, this is the first report of a genome-wide analysis of the tomato CCCH zinc finger family. Our data provided valuable information on tomato CCCH proteins and form a foundation for future studies of these proteins, especially for those members that may play important roles in stress responses.

Abscisic acid perception and signaling: structural mechanisms and applications

Adverse environmental conditions are a threat to agricultural yield and therefore exert a global effect on livelihood, health and the economy. Abscisic acid (ABA) is a vital plant hormone that regulates abiotic stress tolerance, thereby allowing plants to cope with environmental stresses. Previously, attempts to develop a complete understanding of the mechanisms underlying ABA signaling have been hindered by difficulties in the identification of bona fide ABA receptors. The discovery of the PYR/PYL/RCAR family of ABA receptors therefore represented a major milestone in the effort to overcome these roadblocks; since then, many structural and functional studies have provided detailed insights into processes ranging from ABA perception to the activation of ABA-responsive gene transcription. This understanding of the mechanisms of ABA perception and signaling has served as the basis for recent, preliminary developments in the genetic engineering of stress-resistant crops as well as in the design of new synthetic ABA agonists, which hold great promise for the agricultural enhancement of stress tolerance.

'In silico expression analysis', a novel PathoPlant web tool to identify abiotic and biotic stress conditions associated with specific cis-regulatory sequences

Using bioinformatics, putative cis-regulatory sequences can be easily identified using pattern recognition programs on promoters of specific gene sets. The abundance of predicted cis-sequences is a major challenge to associate these sequences with a possible function in gene expression regulation. To identify a possible function of the predicted cis-sequences, a novel web tool designated ‘in silico expression analysis’ was developed that correlates submitted cis-sequences with gene expression data from Arabidopsis thaliana. The web tool identifies the A. thaliana genes harbouring the sequence in a defined promoter region and compares the expression of these genes with microarray data. The result is a hierarchy of abiotic and biotic stress conditions to which these genes are most likely responsive. When testing the performance of the web tool, known cis-regulatory sequences were submitted to the ‘in silico expression analysis’ resulting in the correct identification of the associated stress conditions. When using a recently identified novel elicitor-responsive sequence, a WT-box (CGACTTTT), the ‘in silico expression analysis’ predicts that genes harbouring this sequence in their promoter are most likely Botrytis cinerea induced. Consistent with this prediction, the strongest induction of a reporter gene harbouring this sequence in the promoter is observed with B. cinerea in transgenic A. thaliana.

Database URL:http://www.pathoplant.de/expression_analysis.php.

Determining common insertion sites based on retroviral insertion distribution across tumors

A CIS (common insertion site) indicates a genome region that is hit more frequently by retroviral insertions than expected by chance. Such a region is strongly related to cancer gene loci, which leads to the detection of cancer genes. An algorithm for detecting CISs should satisfy the following: (1) it does not require any prior knowledge of underlying insertion distribution; (2) it can resolve the insertion biases caused by hotspots; (3) it can detect CISs of any biological width; (4) it can identify noises resulting from statistic mistakes and non-CIS insertions; and (5) it can identify the widths of CISs as accurately as possible. We develop a method to resolve these difficulties. We verify a region's significance from two perspectives: distribution width and distribution depth. The former indicates how many insertions in a region while the latter evaluates the insertion distribution across the tumors in a region. We compare our method with kernel density estimation and sliding window on the simulated data, showing that our method not only identifies cancer-related insertions effectively, but also filters noises correctly. The experiments on the real data show that taking insertion distribution into account can highlight significant CISs. We detect 53 novel CISs, some of which have been proven correct by the biological literature.

Cross-family translational genomics of abiotic stress-responsive genes between Arabidopsis and Medicago truncatula

Cross-species translation of genomic information may play a pivotal role in applying biological knowledge gained from relatively simple model system to other less studied, but related, genomes. The information of abiotic stress (ABS)-responsive genes in Arabidopsis was identified and translated into the legume model system, Medicago truncatula. Various data resources, such as TAIR/AtGI DB, expression profiles and literatures, were used to build a genome-wide list of ABS genes. tBlastX/BlastP similarity search tools and manual inspection of alignments were used to identify orthologous genes between the two genomes. A total of 1,377 genes were finally collected and classified into 18 functional criteria of gene ontology (GO). The data analysis according to the expression cues showed that there was substantial level of interaction among three major types (i.e., drought, salinity and cold stress) of abiotic stresses. In an attempt to translate the ABS genes between these two species, genomic locations for each gene were mapped using an in-house-developed comparative analysis platform. The comparative analysis revealed that fragmental colinearity, represented by only 37 synteny blocks, existed between Arabidopsis and M. truncatula. Based on the combination of E-value and alignment remarks, estimated translation rate was 60.2% for this cross-family translation. As a prelude of the functional comparative genomic approaches, in-silico gene network/interactome analyses were conducted to predict key components in the ABS responses, and one of the sub-networks was integrated with corresponding comparative map. The results demonstrated that core members of the sub-network were well aligned with previously reported ABS regulatory networks. Taken together, the results indicate that network-based integrative approaches of comparative and functional genomics are important to interpret and translate genomic information for complex traits such as abiotic stresses.

Responses of Populus trichocarpa galactinol synthase genes to abiotic stresses

Galactinol synthase (GolS; EC 2.4.1.123) is a member of the glycosyltransferase eight family that catalyzes the first step in the biosynthesis pathway of the raffinose family of oligosaccharides (RFOs). The accumulation of RFOs in response to abiotic stress indicates a role for RFOs in stress adaptation. To obtain information on the roles of RFOs in abiotic stress adaptation in trees, we investigated the expression patterns of nine Populus trichocarpa GolS (PtrGolS) genes with special reference to stress responses. PtrGolS genes were differentially expressed in different organs, and the expressions of PtrGolS4 and PtrGolS6 were relatively high in all tested organs. The expression levels of all PtrGolS genes, except PtrGolS9, changed in response to abiotic stress in gene- and stress-type-specific manners. Moreover, short- and long-term stress treatments revealed that induction of PtrGolS by salt stress is obvious only in the early period of treatment (within 24 h), whereas water-deficit stress treatments continued to upregulate PtrGolS gene expression after two days of treatment, in addition to induction within 24 h of treatment. Consistent with these expression patterns, the galactinol content in leaves increased after four days of drought stress, but not under salt stress. Our findings suggest divergent roles for PtrGolS genes in abiotic stress responses in poplars.

Identifying Cis-Regulatory Elements and Modules Using Conditional Random Fields

Accurate identification of cis-regulatory elements and their correlated modules is essential for analysis of transcriptional regulation, which is a challenging problem in computational biology. Unsupervised learning has the advantage of compensating for missing annotated data, and is thus promising to be effective to identify cis-regulatory elements and modules. We introduced a Conditional Random Fields model, referred to as CRFEM, to integrate sequence features and long-range dependency of genomic sequences such as epigenetic features to identify cis-regulatory elements and modules at the same time. The proposed method is able to automatically learn model parameters with no labeled data and explicitly optimize the predictive probability of cis-regulatory elements and modules. In comparison with existing methods, our method is more accurate and can be used for genome-wide studies of gene regulation.

Bidirectional promoters in seed development and related hormone/stress responses

BACKGROUND

Bidirectional promoters are common in genomes but under-studied experimentally, particularly in plants. We describe a targeted identification and selection of a subset of putative bidirectional promoters to identify genes involved in seed development and to investigate possible coordinated responses of gene pairs to conditions important in seed maturation such as desiccation and ABA-regulation.

RESULTS

We combined a search for 100-600 bp intergenic regions in the Arabidopsis genome with a cis-element based selection for those containing multiple copies of the G-box motif, CACGTG. One of the putative bidirectional promoters identified also contained a CE3 coupling element 5 bp downstream of one G-box and is identical to that characterized previously in the HVA1 promoter of barley. CE3 elements are significantly under-represented and under-studied in Arabidopsis. We further characterized the pair of genes associated with this promoter and uncovered roles for two small, previously uncharacterized, plant-specific proteins in Arabidopsis seed development and stress responses.

CONCLUSIONS

Using bioinformatics we identified putative bidirectional promoters involved in seed development and analysed expression patterns for a pair of plant-specific genes in various tissues and in response to hormones/stress. We also present preliminary functional analysis of these genes that is suggestive of roles in seed development.

The role of ABA and MAPK signaling pathways in plant abiotic stress responses

As sessile organisms, plants have developed specific mechanisms that allow them to rapidly perceive and respond to stresses in the environment. Among the evolutionarily conserved pathways, the ABA (abscisic acid) signaling pathway has been identified as a central regulator of abiotic stress response in plants, triggering major changes in gene expression and adaptive physiological responses. ABA induces protein kinases of the SnRK family to mediate a number of its responses. Recently, MAPK (mitogen activated protein kinase) cascades have also been shown to be implicated in ABA signaling. Therefore, besides discussing the role of ABA in abiotic stress signaling, we will also summarize the evidence for a role of MAPKs in the context of abiotic stress and ABA signaling.

TreeTFDB: an integrative database of the transcription factors from six economically important tree crops for functional predictions and comparative and functional genomics

Crop plants, whose productivity is affected by a wide range of growing and environmental conditions, are grown for economic purposes. Transcription factors (TFs) play central role in regulation of many biological processes, including plant development and responses to environmental stimuli, by activating or repressing spatiotemporal gene expression. Here, we describe the TreeTFDB (http://treetfdb.bmep.riken.jp/index.pl) that houses the TF repertoires of six economically important tree crop species: Jatropha curcas, papaya, cassava, poplar, castor bean and grapevine. Among these, the TF repertoire of J. curcas has not been reported by any other TF databases. In addition to their basic information, such as sequence and domain features, domain alignments, gene ontology assignment and sequence comparison, information on available full-length cDNAs, identity and positions of all types of known cis-motifs found in the promoter regions, gene expression data are provided. With its newly designed and friendly interface and its unique features, TreeTFDB will enable research community to predict the functions and provide access to available genetic resources for performing comparative and functional genomics of the crop TFs, either individually or at whole family level, in a comprehensive and convenient manner.

Evolution of variety-specific regulatory schema for expression of osa-miR408 in indica rice varieties under drought stress

Evolution of differential regulatory mechanisms can lead to quite distinct physiological attributes. In the present study, we have identified one such regulatory schema that regulates osa-miR408 and responds differentially in drought-sensitive and -tolerant indica rice varieties. A comparison of the drought stress response in drought-sensitive (Pusa Basmati 1 and IR64) and drought-tolerant (Nagina 22 and Vandana) indica rice varieties revealed that, during drought stress, levels of miR408 transcript decrease significantly in sensitive cultivars, whereas they remain elevated in the tolerant cultivars. The trend is reflected in young seedlings, as well as in flag leaf and spikelets of adult plants (heading stage). Members of the plastocyanin-like protein family targeted by miR408 also show the inverse expression profile and thus accumulate at a lower level in tolerant cultivars during drought. Interestingly, some members of this family are implicated in maintaining the cellular redox state and spikelet fertility in Arabidopsis. An investigation of miR408 loci (including promoter) in all four cultivars did not reveal any significant sequence variation indicating an involvement of the upstream regulatory schema. Indeed, a similar variety-specific stress response was found in the Oryza sativa squamosa promoter-binding-like 9 transcription factor that regulates miR408 expression. We further demonstrate that drought-mediated induction of miR408 in Nagina 22 is regulated by [Ca(2+)]cyt levels. However, [Ca(2+)]cyt does not appear to regulate miR408 levels in Pusa Basmati 1, suggesting a variety-specific evolution of regulatory schema in rice.

Insight into differential responses of upland and paddy rice to drought stress by comparative expression profiling analysis

In this study, the drought responses of two genotypes, IRAT109 and Zhenshan 97 (ZS97), representing upland and paddy rice, respectively, were systematically compared at the morphological, physiological and transcriptional levels. IRAT109 has better performance in traits related to drought avoidance, such as leaf rolling, root volumes, the ratio of leaf water loss and relative conductivity. At the transcriptional level, more genes were induced by drought in IRAT109 at the early drought stage, but more genes had dynamic expression patterns in ZS97 at different drought degrees. Under drought conditions, more genes related to reproductive development and establishment of localization were repressed in IRAT109, but more genes involved in degradation of cellular components were induced in ZS97. By checking the expression patterns of 36 drought-responsive genes (located in 14 quantitative trail loci [QTL] intervals) in ZS97, IRAT109 and near isogenic lines (NILs) of the QTL intervals, we found that more than half of these genes had their expression patterns or expression levels changed in the NILs when compared to that in ZS97 or IRAT109. Our results may provide valuable information for dissecting the genetic bases of traits related to drought resistance, as well as for narrowing the candidate genes for the traits.

Pivotal role of the AREB/ABF-SnRK2 pathway in ABRE-mediated transcription in response to osmotic stress in plants

Water availability is one of the main limiting factors for plant growth and development. The phytohormone abscisic acid (ABA) fulfills a critical role in coordinating the responses to reduced water availability as well as in multiple developmental processes. Endogenous ABA levels increase in response to osmotic stresses such as drought and high salinity, and ABA activates the expression of many genes via ABA-responsive elements (ABREs) in their promoter regions. ABRE-binding protein/ABRE-binding factor (AREB/ABF) transcription factors (TFs) regulate the ABRE-mediated transcription of downstream target genes. Three subclass III sucrose non-fermenting-1 related protein kinase 2 (SnRK2) protein kinases (SRK2D/SnRK2.2, SRK2E/SnRK2.6/OST1 and SRK2I/SnRK2.3) phosphorylate and positively control the AREB/ABF TFs. Substantial progress has been made in our understanding of the ABA-sensing system mediated by Pyrabactin resistance1/PYR1-like/regulatory components of ABA receptor (PYR/PYL/RCAR)-protein phosphatase 2C complexes. In addition to PP2C-PYR/PYL/RCAR ABAreceptor complex, the AREB/ABF-SnRK2 pathway, which is well conserved in land plants, was recently shown to play a major role as a positive regulator of ABA/stress signaling through ABRE-mediated transcription of target genes implicated in the osmotic stress response. This review focuses on current progress in the study of the AREB/ABF-SnRK2 positive regulatory pathway in plants and describes additional signaling factors implicated in the AREB/ABF-SnRK2 pathway. Moreover, to help promote the link between basic and applied studies, the nomenclature and phylogenetic relationships between the AREB/ABFs and SnRK2s are summarized and discussed.

Identification and expression analysis of cytokinin metabolic genes in soybean under normal and drought conditions in relation to cytokinin levels

Cytokinins (CKs) mediate cellular responses to drought stress and targeted control of CK metabolism can be used to develop drought-tolerant plants. Aiming to manipulate CK levels to improve drought tolerance of soybean cultivars through genetic engineering of CK metabolic genes, we surveyed the soybean genome and identified 14 CK biosynthetic (isopentenyltransferase, GmIPT) and 17 CK degradative (CK dehydrogenase, GmCKX) genes. Comparative analyses of GmIPTs and GmCKXs with Arabidopsis counterparts revealed their similar architecture. The average numbers of abiotic stress-inducible cis-elements per promoter were 0.4 and 1.2 for GmIPT and GmCKX genes, respectively, suggesting that upregulation of GmCKXs, thereby reduction of CK levels, maybe the major events under abiotic stresses. Indeed, the expression of 12 GmCKX genes was upregulated by dehydration in R2 roots. Overall, the expressions of soybean CK metabolic genes in various tissues at various stages were highly responsive to drought. CK contents in various organs at the reproductive (R2) stage were also determined under well-watered and drought stress conditions. Although tRNA-type GmIPT genes were highly expressed in soybean, cis-zeatin and its derivatives were found at low concentrations. Moreover, reduction of total CK content in R2 leaves under drought was attributable to the decrease in dihydrozeatin levels, suggesting a role of this molecule in regulating soybean's responses to drought stress. Our systematic analysis of the GmIPT and GmCKX families has provided an insight into CK metabolism in soybean under drought stress and a solid foundation for in-depth characterization and future development of improved drought-tolerant soybean cultivars by manipulation of CK levels via biotechnological approach.

CCCH-type zinc finger family in maize: genome-wide identification, classification and expression profiling under abscisic acid and drought treatments

Background

CCCH-type zinc finger proteins comprise a large protein family. Increasing evidence suggests that members of this family are RNA-binding proteins with regulatory functions in mRNA processing. Compared with those in animals, functions of CCCH-type zinc finger proteins involved in plant growth and development are poorly understood.

Methodology/Principal Findings

Here, we performed a genome-wide survey of CCCH-type zinc finger genes in maize (Zea mays L.) by describing the gene structure, phylogenetic relationships and chromosomal location of each family member. Promoter sequences and expression profiles of putative stress-responsive members were also investigated. A total of 68 CCCH genes (ZmC3H1-68) were identified in maize and divided into seven groups by phylogenetic analysis. These 68 genes were found to be unevenly distributed on 10 chromosomes with 15 segmental duplication events, suggesting that segmental duplication played a major role in expansion of the maize CCCH family. The Ka/Ks ratios suggested that the duplicated genes of the CCCH family mainly experienced purifying selection with limited functional divergence after duplication events. Twelve maize CCCH genes grouped with other known stress-responsive genes from Arabidopsis were found to contain putative stress-responsive cis-elements in their promoter regions. Seven of these genes chosen for further quantitative real-time PCR analysis showed differential expression patterns among five representative maize tissues and over time in response to abscisic acid and drought treatments.

Conclusions

The results presented in this study provide basic information on maize CCCH proteins and form the foundation for future functional studies of these proteins, especially for those members of which may play important roles in response to abiotic stresses.

An omics approach to understand the plant abiotic stress

Abiotic stress can lead to changes in development, productivity, and severe stress and may even threaten survival of plants. Several environmental stresses cause drastic changes in the growth, physiology, and metabolism of plants leading to the increased accumulation of secondary metabolites. As medicinal plants are important sources of drugs, steps are taken to understand the effect of stress on the physiology, biochemistry, genomic, proteomic, and metabolic levels. The molecular responses of plants to abiotic stress are often considered as a complex process. They are mainly based on the modulation of transcriptional activity of stress-related genes. Many genes have been induced under stress conditions. The products of stress-inducible genes protecting against these stresses includes the enzymes responsible for the synthesis of various osmoprotectants. Genetic engineering of tolerance to abiotic stresses help in molecular understanding of pathways induced in response to one or more of the abiotic stresses. Systems biology and virtual experiments allow visualizing and understanding how plants work to overcome abiotic stress. This review discusses the omic approach to understand the plant response to abiotic stress with special emphasis on medicinal plant.

An ABRE promoter sequence is involved in osmotic stress-responsive expression of the DREB2A gene, which encodes a transcription factor regulating drought-inducible genes in Arabidopsis

In plants, osmotic stress-responsive transcriptional regulation depends mainly on two major classes of cis-acting elements found in the promoter regions of stress-inducible genes: ABA-responsive elements (ABREs) and dehydration-responsive elements (DREs). ABRE has been shown to perceive ABA-mediated osmotic stress signals, whereas DRE is known to be involved in an ABA-independent pathway. Previously, we reported that the transcription factor DRE-BINDING PROTEIN 2A (DREB2A) regulates DRE-mediated transcription of target genes under osmotic stress conditions in Arabidopsis (Arabidopsis thaliana). However, the transcriptional regulation of DREB2A itself remains largely uncharacterized. To elucidate the transcriptional mechanism associated with the DREB2A gene under osmotic stress conditions, we generated a series of truncated and base-substituted variants of the DREB2A promoter and evaluated their transcriptional activities individually. We found that both ABRE and coupling element 3 (CE3)-like sequences located approximately -100 bp from the transcriptional initiation site are necessary for the dehydration-responsive expression of DREB2A. Coupling our transient expression analyses with yeast one-hybrid and chromatin immunoprecipitation (ChIP) assays indicated that the ABRE-BINDING PROTEIN 1 (AREB1), AREB2 and ABRE-BINDING FACTOR 3 (ABF3) bZIP transcription factors can bind to and activate the DREB2A promoter in an ABRE-dependent manner. Exogenous ABA application induced only a modest accumulation of the DREB2A transcript when compared with the osmotic stress treatment. However, the osmotic stress-induced DREB2A expression was found to be markedly impaired in several ABA-deficient and ABA-insensitive mutants. These results suggest that in addition to an ABA-independent pathway, the ABA-dependent pathway plays a positive role in the osmotic stress-responsive expression of DREB2A.

Specific expression of LATERAL SUPPRESSOR is controlled by an evolutionarily conserved 3' enhancer

Aerial plant architecture is largely based on the activity of axillary meristems (AMs), initiated in the axils of leaves. The Arabidopsis gene LATERAL SUPPRESSOR (LAS), which is expressed in well-defined domains at the adaxial boundary of leaf primordia, is a key regulator of AM formation. The precise definition of organ boundaries is an essential step for the formation of new organs in general and for meristem initiation; however, mechanisms leading to these specific patterns are not well understood. To increase understanding of how the highly specific transcript accumulation in organ boundary regions is established, we investigated the LAS promoter. Analysis of deletion constructs revealed that an essential enhancer necessary for complementation is situated about 3.2 kb downstream of the LAS open reading frame. This enhancer is sufficient to confer promoter specificity as upstream sequences in LAS could be replaced by non-specific promoters, such as the 35S minimal promoter. Further promoter swapping experiments using the PISTILLATA or the full 35S promoter demonstrated that the LAS 3' enhancer also has suppressor functions, largely overwriting the activity of different 5' promoters. Phylogenetic analyses suggest that LAS function and regulation are evolutionarily highly conserved. Homologous elements in downstream regulatory sequences were found in all LAS orthologs, including grasses. Transcomplementation experiments demonstrated the functional conservation of non-coding sequences between Solanum lycopersicum (tomato) and Arabidopsis. In summary, our results show that a highly conserved enhancer/suppressor element is the main regulatory module conferring the boundary-specific expression of LAS.

Prediction of drought-resistant genes in Arabidopsis thaliana using SVM-RFE

BACKGROUND

Identifying genes with essential roles in resisting environmental stress rates high in agronomic importance. Although massive DNA microarray gene expression data have been generated for plants, current computational approaches underutilize these data for studying genotype-trait relationships. Some advanced gene identification methods have been explored for human diseases, but typically these methods have not been converted into publicly available software tools and cannot be applied to plants for identifying genes with agronomic traits.

METHODOLOGY

In this study, we used 22 sets of Arabidopsis thaliana gene expression data from GEO to predict the key genes involved in water tolerance. We applied an SVM-RFE (Support Vector Machine-Recursive Feature Elimination) feature selection method for the prediction. To address small sample sizes, we developed a modified approach for SVM-RFE by using bootstrapping and leave-one-out cross-validation. We also expanded our study to predict genes involved in water susceptibility.

CONCLUSIONS

We analyzed the top 10 genes predicted to be involved in water tolerance. Seven of them are connected to known biological processes in drought resistance. We also analyzed the top 100 genes in terms of their biological functions. Our study shows that the SVM-RFE method is a highly promising method in analyzing plant microarray data for studying genotype-phenotype relationships. The software is freely available with source code at http://ccst.jlu.edu.cn/JCSB/RFET/.

In silico analysis of transcription factor repertoires and prediction of stress-responsive transcription factors from six major gramineae plants

The interactions between transcription factors (TFs) and cis-regulatory DNA sequences control gene expression, constituting the essential functional linkages of gene regulatory networks. The aim of this study is to identify and integrate all putative TFs from six grass species: Brachypodium distachyon, maize, rice, sorghum, barley, and wheat with significant information into an integrative database (GramineaeTFDB) for comparative genomics and functional genomics. For each TF, sequence features, promoter regions, domain alignments, GO assignment, FL-cDNA information, if available, and cross-references to various public databases and genetic resources are provided. Additionally, GramineaeTFDB possesses a tool which aids the users to search for putative cis-elements located in the promoter regions of TFs and predict the functions of the TFs using cis-element-based functional prediction approach. We also supplied hyperlinks to expression profiles of those TF genes of maize, rice, and barley, for which data are available. Furthermore, information about the availability of FOX and Ds mutant lines for rice and maize TFs, respectively, are also accessible through hyperlinks. Our study provides an important user-friendly public resource for functional analyses and comparative genomics of grass TFs, and understanding of the architecture of transcriptional regulatory networks and evolution of the TFs in agriculturally important cereal crops.

ABA-mediated transcriptional regulation in response to osmotic stress in plants

The plant hormone abscisic acid (ABA) plays a pivotal role in a variety of developmental processes and adaptive stress responses to environmental stimuli in plants. Cellular dehydration during the seed maturation and vegetative growth stages induces an increase in endogenous ABA levels, which control many dehydration-responsive genes. In Arabidopsis plants, ABA regulates nearly 10% of the protein-coding genes, a much higher percentage than other plant hormones. Expression of the genes is mainly regulated by two different families of bZIP transcription factors (TFs), ABI5 in the seeds and AREB/ABFs in the vegetative stage, in an ABA-responsive-element (ABRE) dependent manner. The SnRK2-AREB/ABF pathway governs the majority of ABA-mediated ABRE-dependent gene expression in response to osmotic stress during the vegetative stage. In addition to osmotic stress, the circadian clock and light conditions also appear to participate in the regulation of ABA-mediated gene expression, likely conferring versatile tolerance and repressing growth under stress conditions. Moreover, various other TFs belonging to several classes, including AP2/ERF, MYB, NAC, and HD-ZF, have been reported to engage in ABA-mediated gene expression. This review mainly focuses on the transcriptional regulation of ABA-mediated gene expression in response to osmotic stress during the vegetative growth stage in Arabidopsis.

Microarray-based analysis of cadmium-responsive microRNAs in rice (Oryza sativa)

MicroRNAs (miRNAs) are a class of small non-coding RNAs that negatively regulate specific target mRNAs at the post-transcriptional level. Plant miRNAs have been implicated in developmental processes and adaptations to environmental stresses. Cadmium (Cd) is a non-essential heavy metal that is highly toxic to plants. To investigate the responsive functions of miRNAs under Cd stress, miRNA expression in Cd-stressed rice (Oryza sativa) was profiled using a microarray assay. A total of 19 Cd-responsive miRNAs were identified, of which six were further validated experimentally. Target genes were also predicted for these Cd-responsive miRNAs, which encoded transcription factors, and proteins associated with metabolic processes or stress responses. In addition, the mRNA levels of several targets were negatively correlated with the corresponding miRNAs under Cd stress. Promoter analysis showed that metal stress-responsive cis-elements tended to occur more frequently in the promoter regions of Cd-responsive miRNAs. These findings suggested that miRNAs played an important role in Cd tolerance in rice, and highlighted a novel molecular mechanism of heavy metal tolerance in plants.

Prediction of new abiotic stress genes in Arabidopsis thaliana and Oryza sativa according to enumeration-based statistical analysis

Plants undergo an extensive change in gene regulation during abiotic stress. It is of great agricultural importance to know which genes are affected during stress response. The genome sequence of a number of plant species has been determined, among them Arabidopsis and Oryza sativa, whose genome has been annotated most completely as of yet, and are well-known organisms widely used as experimental systems. This paper applies a statistical algorithm for predicting new stress-induced motifs and genes by analyzing promoter sets co-regulated by abiotic stress in the previously mentioned two species. After identifying characteristic putative regulatory motif sequence pairs (dyads) in the promoters of 125 stress-regulated Arabidopsis genes and 87 O. sativa genes, these dyads were used to screen the entire Arabidopsis and O. sativa promoteromes to find related stress-induced genes whose promoters contained a large number of these dyads found by our algorithm. We were able to predict a number of putative dyads, characteristic of a large number of stress-regulated genes, some of them newly discovered by our algorithm and serve as putative transcription factor binding sites. Our new motif prediction algorithm comes complete with a stand-alone program. This algorithm may be used in motif discovery in the future in other species. The more than 1,200 Arabidopsis and 1,700 Orzya sativa genes found by our algorithm are good candidates for further experimental studies in abiotic stress.

Genome-wide expression profiling of soybean two-component system genes in soybean root and shoot tissues under dehydration stress

Two-component systems (TCSs) play vital functions in the adaptation of plants to environmental stresses. To identify soybean TCS genes involved in the regulation of drought stress response, we performed tissue-specific expression profiling of all 83 putative TCS genes in plants subjected to dehydration. Under well-watered conditions, the majority of soybean TCS genes were expressed higher in the root tissues. Additionally, a high variability in transcript abundance was observed for the TCS genes in both roots and shoots. Under dehydration, TCS genes were more responsive in shoots than in roots. Further analysis indicated that 50% more TCS genes were repressed by dehydration than induced. Specifically, 18 genes were induced by 2-fold or more, whereas 33 genes were down-regulated at least 2-fold by dehydration. TCS genes putatively involved in cytokinin and ethylene signallings strongly responded to dehydration, suggesting that crosstalk exists between different hormonal and stress pathways. Our study provides the first glance into the complex regulatory roles of soybean TCSs underlying their functions in response to dehydration. Additionally, these systematic expression analyses identified excellent dehydration-responsive candidate genes to further clarify soybean TCS functions in drought response and to enable the development of improved drought tolerance in transgenic soybeans.

Molecular basis of the core regulatory network in ABA responses: sensing, signaling and transport

ABA is a major phytohormone that regulates a broad range of plant traits and is especially important for adaptation to environmental conditions. Our understanding of the molecular basis of ABA responses in plants improved dramatically in 2009 and 2010, banner years for ABA research. There are three major components; PYR/PYL/ RCAR (an ABA receptor), type 2C protein phosphatase (PP2C; a negative regulator) and SNF1-related protein kinase 2 (SnRK2; a positive regulator), and they offer a double negative regulatory system, [PYR/PYL/RCAR-| PP2C-| SnRK2]. In the absence of ABA, PP2C inactivates SnRK2 by direct dephosphorylation. In response to environmental or developmental cues, ABA promotes the interaction of PYR/PYL/RCAR and PP2C, resulting in PP2C inhibition and SnRK2 activation. This signaling complex can work in both the nucleus and cytosol, as it has been shown that SnRK2 phosphorylates basic-domain leucine zipper (bZIP) transcription factors or membrane proteins. Several structural analyses of PYR/PYL/RCAR have provided the mechanistic basis for this 'core signaling' model, by elucidating the mechanism of ABA binding of receptors, or the 'gate-latch-lock' mechanism of interaction with PP2C in inhibiting activity. On the other hand, intercellular ABA transport had remained a major issue, as had intracellular ABA signaling. Recently, two plasma membrane-type ABC transporters were identified and shed light on the influx/efflux system of ABA, resolving how ABA is transported from cell to cell in plants. Our knowledge of ABA responses in plants has been greatly expanded from intracellular signaling to intercellular transport of ABA.

Characterization of the beta-carotene hydroxylase gene DSM2 conferring drought and oxidative stress resistance by increasing xanthophylls and abscisic acid synthesis in rice

Drought is a major limiting factor for crop production. To identify critical genes for drought resistance in rice (Oryza sativa), we screened T-DNA mutants and identified a drought-hypersensitive mutant, dsm2. The mutant phenotype was caused by a T-DNA insertion in a gene encoding a putative β-carotene hydroxylase (BCH). BCH is predicted for the biosynthesis of zeaxanthin, a carotenoid precursor of abscisic acid (ABA). The amounts of zeaxanthin and ABA were significantly reduced in two allelic dsm2 mutants after drought stress compared with the wild type. Under drought stress conditions, the mutant leaves lost water faster than the wild type and the photosynthesis rate, biomass, and grain yield were significantly reduced, whereas malondialdehyde level and stomata aperture were increased in the mutant. The mutant is also hypersensitive to oxidative stresses. The mutant had significantly lower maximal efficiency of photosystem II photochemistry and nonphotochemical quenching capacity than the wild type, indicating photoinhibition in photosystem II and decreased capacity for eliminating excess energy by thermal dissipation. Overexpression of DSM2 in rice resulted in significantly increased resistance to drought and oxidative stresses and increases of the xanthophylls and nonphotochemical quenching. Some stress-related ABA-responsive genes were up-regulated in the overexpression line. DSM2 is a chloroplast protein, and the response of DSM2 to environmental stimuli is distinctive from the other two BCH members in rice. We conclude that the DSM2 gene significantly contributes to control of the xanthophyll cycle and ABA synthesis, both of which play critical roles in the establishment of drought resistance in rice.

Genome-wide identification and analysis of drought-responsive microRNAs in Oryza sativa

In addition to regulating growth and development, the most important function of microRNAs (miRNAs) in plants is the regulation of a variety of cellular processes underlying plant adaptation to environmental stresses. To gain a deep understanding of the mechanism of drought tolerance in rice, genome-wide profiling and analysis of miRNAs was carried out in drought-challenged rice across a wide range of developmental stages, from tillering to inflorescence formation, using a microarray platform. Among the 30 miRNAs identified as significantly down- or up-regulated under the drought stress, 11 down-regulated miRNAs (miR170, miR172, miR397, miR408, miR529, miR896, miR1030, miR1035, miR1050, miR1088, and miR1126) and eight up-regulated miRNAs (miR395, miR474, miR845, miR851, miR854, miR901, miR903, and miR1125) were revealed for the first time to be induced by drought stress in plants, and nine (miR156, miR168, miR170, miR171, miR172, miR319, miR396, miR397, and miR408) showed opposite expression to that observed in drought-stressed Arabidopsis. The most conserved down-regulated miRNAs were ath-miR170, the miR171 family, and ath-miR396, and the most conserved up-regulated miRNAs were ptc-miR474 and ath-miR854a. The identification of differentially expressed novel plant miRNAs and their target genes, and the analysis of cis-elements provides molecular evidence for the possible involvement of miRNAs in the process of drought response and/or tolerance in rice.

Genome-wide analysis of two-component systems and prediction of stress-responsive two-component system members in soybean

In plants, the two-component systems (TCSs) play important roles in regulating diverse biological processes, including responses to environmental stress stimuli. Within the soybean genome, the TCSs consist of at least 21 histidine kinases, 13 authentic and pseudo-phosphotransfers and 18 type-A, 15 type-B, 3 type-C and 11 pseudo-response regulator proteins. Structural and phylogenetic analyses of soybean TCS members with their Arabidopsis and rice counterparts revealed similar architecture of their TCSs. We identified a large number of closely homologous soybean TCS genes, which likely resulted from genome duplication. Additionally, we analysed tissue-specific expression profiles of those TCS genes, whose data are available from public resources. To predict the putative regulatory functions of soybean TCS members, with special emphasis on stress-responsive functions, we performed comparative analyses from all the TCS members of soybean, Arabidopsis and rice and coupled these data with annotations of known abiotic stress-responsive cis-elements in the promoter region of each soybean TCS gene. Our study provides insights into the architecture and a solid foundation for further functional characterization of soybean TCS elements. In addition, we provide a new resource for studying the conservation and divergence among the TCSs within plant species and/or between plants and other organisms.

Early abscisic acid signal transduction mechanisms: newly discovered components and newly emerging questions

The plant hormone abscisic acid (ABA) regulates many key processes in plants, including seed germination and development and abiotic stress tolerance, particularly drought resistance. Understanding early events in ABA signal transduction has been a major goal of plant research. The recent identification of the PYRABACTIN (4-bromo-N-[pyridin-2-yl methyl]naphthalene-1-sulfonamide) RESISTANCE (PYR)/REGULATORY COMPONENT OF ABA RECEPTOR (RCAR) family of ABA receptors and their biochemical mode of action represents a major breakthrough in the field. The solving of PYR/RCAR structures provides a context for resolving mechanisms mediating ABA control of protein-protein interactions for downstream signaling. Recent studies show that a pathway based on PYR/RCAR ABA receptors, PROTEIN PHOSPHATASE 2Cs (PP2Cs), and SNF1-RELATED PROTEIN KINASE 2s (SnRK2s) forms the primary basis of an early ABA signaling module. This pathway interfaces with ion channels, transcription factors, and other targets, thus providing a mechanistic connection between the phytohormone and ABA-induced responses. This emerging PYR/RCAR-PP2C-SnRK2 model of ABA signal transduction is reviewed here, and provides an opportunity for testing novel hypotheses concerning ABA signaling. We address newly emerging questions, including the potential roles of different PYR/RCAR isoforms, and the significance of ABA-induced versus constitutive PYR/RCAR-PP2C interactions. We also consider how the PYR/RCAR-PP2C-SnRK2 pathway interfaces with ABA-dependent gene expression, ion channel regulation, and control of small molecule signaling. These exciting developments provide researchers with a framework through which early ABA signaling can be understood, and allow novel questions about the hormone response pathway and possible applications in stress resistance engineering of plants to be addressed.

Arabidopsis thaliana plant defensin AtPDF1.1 is involved in the plant response to biotic stress

*Previously, it was shown that the Arabidopsis thaliana plant defensins AtPDF1.1 (At1g75830) and AtPDF1.2a (At5g44420) exert in vitro antimicrobial properties and that their corresponding genes are expressed in seeds and induced in leaves upon pathogen attack, respectively. *In this study, the expression profile of both AtPDF1.1 and AtPDF1.2a is analysed in wild-type plants upon different stress-related treatments and the effect of modulation of their expression in transgenic plants is examined in both host and nonhost resistance. *AtPDF1.1, which was originally considered to be seed-specific, is demonstrated to be locally induced in leaves upon fungal attack and exhibits an expression profile distinct from that of AtPDF1.2a, a gene frequently used as marker for the ethylene/jasmonate-mediated signaling pathway. Transgenic plants with modulated AtPDF1.1 or AtPDF1.2a gene expression show no altered phenotype upon Botrytis cinerea inoculation. However, constitutive overexpression of AtPDF1.1 in A. thaliana leads to a reduction in symptoms caused by the nonhost Cercospora beticola causing non-spreading spots on A. thaliana leaves. *These results indicate that AtPDF1.1 and AtPDF1.2a clearly differ regarding their expression profile and functionality in planta. It emphasizes the additional level of complexity and fine-tuning within the highly redundant plant defensin genes in A. thaliana.

Identification and prediction of abiotic stress responsive transcription factors involved in abiotic stress signaling in soybean

Abiotic stresses such as extreme temperature, drought, high salinity, cold and waterlogging often result in significant losses to the yields of economically important crops such as soybean (Glycine max L.). Transcription factors (TFs) which bind to DNA through specific cis-regulatory sequences either activate or repress gene transcription have been reported to act as control switches in stress signaling. Recent completion of the soybean genomic sequence has open wide opportunities for large-scale identification and annotations of regulatory TFs in soybean for functional studies. Within the soybean genome, we identified 5,035 TF models which grouped into 61 families. Detailed annotations of soybean TF genes can be accessed at SoybeanTFDB (soybeantfdb.psc.riken.jp). Moreover, we have reported a new idea of high throughput prediction and selection of abiotic stress responsive TFs based on the existence of known stress responsive cis-element(s) located in the promoter regions of respective TFs and GO annotations. We, therefore, have provided a basic platform for the genome-wide analysis of regulatory mechanisms underlying abiotic stress responses and a reliable tool for prediction and selection of stress responsive TFs for further functional studies and genetic engineering.

AREB1, AREB2, and ABF3 are master transcription factors that cooperatively regulate ABRE-dependent ABA signaling involved in drought stress tolerance and require ABA for full activation

A myriad of drought stress-inducible genes have been reported, and many of these are activated by abscisic acid (ABA). In the promoter regions of such ABA-regulated genes, conserved cis-elements, designated ABA-responsive elements (ABREs), control gene expression via bZIP-type AREB/ABF transcription factors. Although all three members of the AREB/ABF subfamily, AREB1, AREB2, and ABF3, are upregulated by ABA and water stress, it remains unclear whether these are functional homologs. Here, we report that all three AREB/ABF transcription factors require ABA for full activation, can form hetero- or homodimers to function in nuclei, and can interact with SRK2D/SnRK2.2, an SnRK2 protein kinase that was identified as a regulator of AREB1. Along with the tissue-specific expression patterns of these genes and the subcellular localization of their encoded proteins, these findings clearly indicate that AREB1, AREB2, and ABF3 have largely overlapping functions. To elucidate the role of these AREB/ABF transcription factors, we generated an areb1 areb2 abf3 triple mutant. Large-scale transcriptome analysis, which showed that stress-responsive gene expression is remarkably impaired in the triple mutant, revealed novel AREB/ABF downstream genes in response to water stress, including many LEA class and group-Ab PP2C genes and transcription factors. The areb1 areb2 abf3 triple mutant is more resistant to ABA than are the other single and double mutants with respect to primary root growth, and it displays reduced drought tolerance. Thus, these results indicate that AREB1, AREB2, and ABF3 are master transcription factors that cooperatively regulate ABRE-dependent gene expression for ABA signaling under conditions of water stress.

Identification of iron-deficiency responsive microRNA genes and cis-elements in Arabidopsis

MicroRNAs (miRNAs) are a class of endogenous non-coding small RNAs that bind to their target mRNAs to repress their translation or induce their degradation. Recent studies have shown that several miRNAs regulate plant adaptation to sulfate and phosphate deficiency. However, whether miRNAs are involved in regulation of stress response to iron (Fe) deficiency is unknown. In this study, we carried out a survey of Arabidopsis miRNA genes in response to Fe deficiency and identified IDE1/IDE2 (Iron-deficiency responsive cis-Element 1 and 2) in their promoter regions. We constructed a small RNA library from Arabidopsis seedlings under Fe deficiency. Sequence analysis revealed 8 conserved miRNA genes in 5 families, all of which were up-regulated during Fe deficiency. Further, we analyzed cis-regulatory elements upstream of all miRNA genes in Arabidopsis and found 24 miRNA genes containing IDE1/IDE2 motifs in their promoter regions. Transcriptional analysis using RT-PCR showed that 70.8% (17/24) of the IDE-containing miRNA genes were expressed in response to Fe deficiency. We presented a putative interaction model between protein-coding genes and miRNA genes under Fe deficiency. Our analytic approach is useful and efficient because it is applicable to cis-element finding for miRNAs responding to other abiotic stresses. Also, the data obtained in this study may aid our understanding of the role of Fe deficiency responsive specific sequences upstream of miRNA genes and the functional implications of miRNA genes in response to Fe stress in plants.

LegumeTFDB: an integrative database of Glycine max, Lotus japonicus and Medicago truncatula transcription factors

We have established a database named LegumeTFDB to provide access to transcription factor (TF) repertoires of three major legume species: soybean (Glycine max), Lotus japonicus and Medicago truncatula. LegumeTFDB integrates unique information for each TF gene and family, including sequence features, gene promoters, domain alignments, gene ontology (GO) assignment and sequence comparison data derived from comparative analysis with TFs found within legumes, in Arabidopsis, rice and poplar as well as with proteins in NCBI nr and UniProt. We also analyzed the promoter regions for all of the TFs to identify all types of cis-motifs provided by the PLACE database. Additionally, we supply hyperlinks to make available expression data of 2411 soybean TF genes. LegumeTFDB provides an important user-friendly public resource for comparative genomics and understanding of transcriptional regulation in agriculturally important legumes.

AVAILABILITY

http://legumetfdb.psc.riken.jp/.

SUPPLEMENTARY INFORMATION

Supplementary data available at Bioinformatics online.

In silico analysis of transcription factor repertoire and prediction of stress responsive transcription factors in soybean

Sequence-specific DNA-binding transcription factors (TFs) are often termed as 'master regulators' which bind to DNA and either activate or repress gene transcription. We have computationally analysed the soybean genome sequence data and constructed a proper set of TFs based on the Hidden Markov Model profiles of DNA-binding domain families. Within the soybean genome, we identified 4342 loci encoding 5035 TF models which grouped into 61 families. We constructed a database named SoybeanTFDB (http://soybeantfdb.psc.riken.jp) containing the full compilation of soybean TFs and significant information such as: functional motifs, full-length cDNAs, domain alignments, promoter regions, genomic organization and putative regulatory functions based on annotations of gene ontology (GO) inferred by comparative analysis with Arabidopsis. With particular interest in abiotic stress signalling, we analysed the promoter regions for all of the TF encoding genes as a means to identify abiotic stress responsive cis-elements as well as all types of cis-motifs provided by the PLACE database. SoybeanTFDB enables scientists to easily access cis-element and GO annotations to aid in the prediction of TF function and selection of TFs with functions of interest. This study provides a basic framework and an important user-friendly public information resource which enables analyses of transcriptional regulation in soybean.

The word landscape of the non-coding segments of the Arabidopsis thaliana genome

Background

Genome sequences can be conceptualized as arrangements of motifs or words. The frequencies and positional distributions of these words within particular non-coding genomic segments provide important insights into how the words function in processes such as mRNA stability and regulation of gene expression.

Results

Using an enumerative word discovery approach, we investigated the frequencies and positional distributions of all 65,536 different 8-letter words in the genome of Arabidopsis thaliana. Focusing on promoter regions, introns, and 3' and 5' untranslated regions (3'UTRs and 5'UTRs), we compared word frequencies in these segments to genome-wide frequencies. The statistically interesting words in each segment were clustered with similar words to generate motif logos. We investigated whether words were clustered at particular locations or were distributed randomly within each genomic segment, and we classified the words using gene expression information from public repositories. Finally, we investigated whether particular sets of words appeared together more frequently than others.

Conclusion

Our studies provide a detailed view of the word composition of several segments of the non-coding portion of the Arabidopsis genome. Each segment contains a unique word-based signature. The respective signatures consist of the sets of enriched words, 'unwords', and word pairs within a segment, as well as the preferential locations and functional classifications for the signature words. Additionally, the positional distributions of enriched words within the segments highlight possible functional elements, and the co-associations of words in promoter regions likely represent the formation of higher order regulatory modules. This work is an important step toward fully cataloguing the functional elements of the Arabidopsis genome.

Cis-regulatory elements in plant cell signaling

Plant cell signaling pathways are in part dependent on transcriptional regulatory networks comprising circuits of transcription factors (TFs) and regulatory DNA elements that control the expression of target genes. Here, we describe experimental and bioinformatic approaches for identifying potential cis-regulatory elements. We also discuss recent integrative genomics studies aimed at elucidating the functions of cis-regulatory elements in aspects of plant biology, including the circadian clock, interactions with the environment, stress responses, and regulation of growth and development by phytohormones. Finally, we discuss emerging technologies and approaches that offer great potential for accelerating the discovery and functional characterization of cis-elements and interacting TFs--which will help realize the promise of systems biology.

CYP86B1 is required for very long chain omega-hydroxyacid and alpha, omega -dicarboxylic acid synthesis in root and seed suberin polyester

Suberin composition of various plants including Arabidopsis (Arabidopsis thaliana) has shown the presence of very long chain fatty acid derivatives C20 in addition to the C16 and C18 series. Phylogenetic studies and plant genome mining have led to the identification of putative aliphatic hydroxylases belonging to the CYP86B subfamily of cytochrome P450 monooxygenases. In Arabidopsis, this subfamily is represented by CYP86B1 and CYP86B2, which share about 45% identity with CYP86A1, a fatty acid omega-hydroxylase implicated in root suberin monomer synthesis. Here, we show that CYP86B1 is located to the endoplasmic reticulum and is highly expressed in roots. Indeed, CYP86B1 promoter-driven beta-glucuronidase expression indicated strong reporter activities at known sites of suberin production such as the endodermis. These observations, together with the fact that proteins of the CYP86B type are widespread among plant species, suggested a role of CYP86B1 in suberin biogenesis. To investigate the involvement of CYP86B1 in suberin biogenesis, we characterized an allelic series of cyp86B1 mutants of which two strong alleles were knockouts and two weak ones were RNA interference-silenced lines. These root aliphatic plant hydroxylase lines had a root and a seed coat aliphatic polyester composition in which C22- and C24-hydroxyacids and alpha,omega-dicarboxylic acids were strongly reduced. However, these changes did not affect seed coat permeability and ion content in leaves. The presumed precursors, C22 and C24 fatty acids, accumulated in the suberin polyester. These results demonstrate that CYP86B1 is a very long chain fatty acid hydroxylase specifically involved in polyester monomer biosynthesis during the course of plant development.

Cis-element- and transcriptome-based screening of root hair-specific genes and their functional characterization in Arabidopsis

Understanding the cellular differentiation of multicellular organisms requires the characterization of genes whose expression is modulated in a cell type-specific manner. The Arabidopsis (Arabidopsis thaliana) root hair cell is one model for studying cellular differentiation. In this study, root hair cell-specific genes were screened by a series of in silico and experimental filtration procedures. This process included genome-wide screening for genes with a root hair-specific cis-element in their promoters, filtering root-specific genes from the root hair-specific cis-element-containing genes, further filtering of genes that were suppressed in root hair-defective plant lines, and experimental confirmation by promoter assay. These procedures revealed 19 root hair-specific genes, including many protein kinases and cell wall-related genes, most of which have not been characterized thus far. Functional analyses of these root hair-specific genes with loss-of-function mutants and overexpressing transformants revealed that they play roles in hair growth and morphogenesis. This study demonstrates that a defined cis-element can serve as a filter to screen certain cell type-specific genes and implicates many new root hair-specific genes in root hair development.

Comparative expression and transcript initiation of three peach dehydrin genes

Dehydrin genes encode proteins with demonstrated cryoprotective and antifreeze activity, and they respond to a variety of abiotic stress conditions that have dehydration as a common component. Two dehydrins from peach (Prunus persica L. [Batsch.]) have been previously characterized; here, we describe the characterization of a third dehydrin from peach bark, PpDhn3, isolated by its response to low temperature. The expression of all three dehydrin genes was profiled by semi-quantitative reverse transcription PCR, and transcript initiation was mapped for all three genes using the RNA ligase-mediated 5' rapid amplification of cDNA ends technique. PpDhn3 transcripts from bark collected in December or July, as well as transcripts from developing fruit, initiated at a single site. Although most of the PpDhn1 transcripts initiated at a similar position, those from young fruit initiated much further upstream of the consensus TATA box. Bark and fruit transcripts encoding PpDhn2 initiated ca. 30 bases downstream of a consensus TATA box; however, transcripts from ripe fruit initiated further upstream. Ripe fruit transcripts of PpDhn2 contain a 5' leader intron which is predicted to add some 34 amino acids to the N-terminal methionine of the cognate protein when properly processed. Secondary structure prediction of sequences surrounding the TATA box suggests that conformational transitions associated with decreasing temperature contribute to the regulation of expression of the cold-responsive dehydrin genes. Taken together these results reveal new, unexpected levels of gene regulation contributing to the overall expression pattern of peach dehydrins.

Unraveling transcriptional control in Arabidopsis using cis-regulatory elements and coexpression networks

Analysis of gene expression data generated by high-throughput microarray transcript profiling experiments has demonstrated that genes with an overall similar expression pattern are often enriched for similar functions. This guilt-by-association principle can be applied to define modular gene programs, identify cis-regulatory elements, or predict gene functions for unknown genes based on their coexpression neighborhood. We evaluated the potential to use Gene Ontology (GO) enrichment of a gene's coexpression neighborhood as a tool to predict its function but found overall low sensitivity scores (13%-34%). This indicates that for many functional categories, coexpression alone performs poorly to infer known biological gene functions. However, integration of cis-regulatory elements shows that 46% of the gene coexpression neighborhoods are enriched for one or more motifs, providing a valuable complementary source to functionally annotate genes. Through the integration of coexpression data, GO annotations, and a set of known cis-regulatory elements combined with a novel set of evolutionarily conserved plant motifs, we could link many genes and motifs to specific biological functions. Application of our coexpression framework extended with cis-regulatory element analysis on transcriptome data from the cell cycle-related transcription factor OBP1 yielded several coexpressed modules associated with specific cis-regulatory elements. Moreover, our analysis strongly suggests a feed-forward regulatory interaction between OBP1 and the E2F pathway. The ATCOECIS resource (http://bioinformatics.psb.ugent.be/ATCOECIS/) makes it possible to query coexpression data and GO and cis-regulatory element annotations and to submit user-defined gene sets for motif analysis, providing an access point to unravel the regulatory code underlying transcriptional control in Arabidopsis (Arabidopsis thaliana).

Between a rock and a dry place: the water-stressed moss

The earliest land plants faced a suite of abiotic stresses largely unknown to their aquatic algal ancestors. The descendants of these plants evolved two general mechanisms for survival in the relatively arid aerial environment. While the vascular plants or 'tracheophytes' developed tissue specializations to transport and retain water, the other main lineages of land plants, the bryophytes, retained a simple, nonvascular morphology. The bryophytes--mosses, hornworts, and liverworts--continually undergo a co-equilibration of their water content with the surrounding environment and rely to a great extent on intrinsic cellular mechanisms to mitigate damage due to water stress. This short review will focus on the cellular and molecular responses to dehydration and rehydration in mosses, and offer insights into general plant responses to water stress.

An algorithm for identifying novel targets of transcription factor families: application to hypoxia-inducible factor 1 targets

Efficient and effective analysis of the growing genomic databases requires the development of adequate computational tools. We introduce a fast method based on the suffix tree data structure for predicting novel targets of hypoxia-inducible factor 1 (HIF-1) from huge genome databases. The suffix tree data structure has two powerful applications here: one is to extract unknown patterns from multiple strings/sequences in linear time; the other is to search multiple strings/sequences using multiple patterns in linear time. Using 15 known HIF-1 target gene sequences as a training set, we extracted 105 common patterns that all occur in the 15 training genes using suffix trees. Using these 105 common patterns along with known subsequences surrounding HIF-1 binding sites from the literature, the algorithm searches a genome database that contains 2,078,786 DNA sequences. It reported 258 potentially novel HIF-1 targets including 25 known HIF-1 targets. Based on microarray studies from the literature, 17 putative genes were confirmed to be upregulated by HIF-1 or hypoxia inside these 258 genes. We further studied one of the potential targets, COX-2, in the biological lab; and showed that it was a biologically relevant HIF-1 target. These results demonstrate that our methodology is an effective computational approach for identifying novel HIF-1 targets.

Genome-wide targeted prediction of ABA responsive genes in rice based on over-represented cis-motif in co-expressed genes

Abscisic acid (ABA), the popular plant stress hormone, plays a key role in regulation of sub-set of stress responsive genes. These genes respond to ABA through specific transcription factors which bind to cis-regulatory elements present in their promoters. We discovered the ABA Responsive Element (ABRE) core (ACGT) containing CGMCACGTGB motif as over-represented motif among the promoters of ABA responsive co-expressed genes in rice. Targeted gene prediction strategy using this motif led to the identification of 402 protein coding genes potentially regulated by ABA-dependent molecular genetic network. RT-PCR analysis of arbitrarily chosen 45 genes from the predicted 402 genes confirmed 80% accuracy of our prediction. Plant Gene Ontology (GO) analysis of ABA responsive genes showed enrichment of signal transduction and stress related genes among diverse functional categories.

An in silico strategy identified the target gene candidates regulated by dehydration responsive element binding proteins (DREBs) in Arabidopsis genome

Identification of downstream target genes of stress-relating transcription factors (TFs) is desirable in understanding cellular responses to various environmental stimuli. However, this has long been a difficult work for both experimental and computational practices. In this research, we presented a novel computational strategy which combined the analysis of the transcription factor binding site (TFBS) contexts and machine learning approach. Using this strategy, we conducted a genome-wide investigation into novel direct target genes of dehydration responsive element binding proteins (DREBs), the members of AP2-EREBPs transcription factor super family which is reported to be responsive to various abiotic stresses in Arabidopsis. The genome-wide searching yielded in total 474 target gene candidates. With reference to the microarray data for abiotic stresses-inducible gene expression profile, 268 target gene candidates out of the total 474 genes predicted, were induced during the 24-h exposure to abiotic stresses. This takes about 57% of total predicted targets. Furthermore, GO annotations revealed that these target genes are likely involved in protein amino acid phosphorylation, protein binding and Endomembrane sorting system. The results suggested that the predicted target gene candidates were adequate to meet the essential biological principle of stress-resistance in plants.

Arabidopsis transcriptome reveals control circuits regulating redox homeostasis and the role of an AP2 transcription factor

Sensors and regulatory circuits that maintain redox homeostasis play a central role in adjusting plant metabolism and development to changing environmental conditions. We report here control networks in Arabidopsis (Arabidopsis thaliana) that respond to photosynthetic stress. We independently subjected Arabidopsis leaves to two commonly used photosystem II inhibitors: high light (HL) and 3-(3,4-dichlorophenyl)-1,1-dimethylurea (DCMU). Microarray analysis of expression patterns during the period of redox adjustment to these inhibitors reveals that 20% and 8% of the transcriptome are under HL and DCMU regulation, respectively. Approximately 6% comprise a subset of genes common to both perturbations, the redox responsive genes (RRGs). A redox network was generated in an attempt to identify genes whose expression is tightly coordinated during adjustment to homeostasis, using expression of these RRGs under HL conditions. Ten subnetworks were identified from the network. Hierarchal subclustering of subnetworks responding to the DCMU stress identified novel groups of genes that were tightly controlled while adjusting to homeostasis. Upstream analysis of the promoters of the genes in these clusters revealed different motifs for each subnetwork, including motifs that were previously identified with responses to other stresses, such as light, dehydration, or abscisic acid. Functional categorization of RRGs demonstrated involvement of genes in many metabolic pathways, including several families of transcription factors, especially those in the AP2 family. Using a T-DNA insertion in one AP2 transcription factor (redox-responsive transcription factor 1 [RRTF1]) from the RRGs, we showed that the genes predicted to be within the subnetwork containing RRTF1 were changed in this insertion line (Deltarrtf1). Furthermore, Deltarrtf1 showed greater sensitivity to photosynthetic stress compared to the wild type.