AthaMap web tools for the analysis of transcriptional and posttranscriptional regulation of gene expression in Arabidopsis thaliana

Identification of gene function based on models capturing natural variability of Arabidopsis thaliana lipid metabolism

Lipids play fundamental roles in regulating agronomically important traits. Advances in plant lipid metabolism have until recently largely been based on reductionist approaches, although modulation of its components can have system-wide effects. However, existing models of plant lipid metabolism provide lumped representations, hindering detailed study of component modulation. Here, we present the Plant Lipid Module (PLM) which provides a mechanistic description of lipid metabolism in the Arabidopsis thaliana rosette. We demonstrate that the PLM can be readily integrated in models of A. thaliana Col-0 metabolism, yielding accurate predictions (83%) of single lethal knock-outs and 75% concordance between measured transcript and predicted flux changes under extended darkness. Genome-wide associations with fluxes obtained by integrating the PLM in diel condition- and accession-specific models identify up to 65 candidate genes modulating A. thaliana lipid metabolism. Using mutant lines, we validate up to 40% of the candidates, paving the way for identification of metabolic gene function based on models capturing natural variability in metabolism.

XAANTAL1 Reveals an Additional Level of Flowering Regulation in the Shoot Apical Meristem in Response to Light and Increased Temperature in Arabidopsis

Light and photoperiod are environmental signals that regulate flowering transition. In plants like Arabidopsis thaliana, this regulation relies on CONSTANS, a transcription factor that is negatively posttranslational regulated by phytochrome B during the morning, while it is stabilized by PHYA and cryptochromes 1/2 at the end of daylight hours. CO induces the expression of FT, whose protein travels from the leaves to the apical meristem, where it binds to FD to regulate some flowering genes. Although PHYB delays flowering, we show that light and PHYB positively regulate XAANTAL1 and other flowering genes in the shoot apices. Also, the genetic data indicate that XAL1 and FD participate in the same signaling pathway in flowering promotion when plants are grown under a long-day photoperiod at 22 °C. By contrast, XAL1 functions independently of FD or PIF4 to induce flowering at higher temperatures (27 °C), even under long days. Furthermore, XAL1 directly binds to FD, SOC1, LFY, and AP1 promoters. Our findings lead us to propose that light and temperature influence the floral network at the meristem level in a partially independent way of the signaling generated from the leaves.

Arabidopsis TIE1 and TIE2 transcriptional repressors dampen cytokinin response during root development

Cytokinin plays critical roles in root development. Cytokinin signaling depends on activation of key transcription factors known as type B Arabidopsis response regulators (ARRs). However, the mechanisms underlying the finely tuned regulation of type B ARR activity remain unclear. In this study, we demonstrate that the ERF-associated amphiphilic repression (EAR) motif-containing protein TCP interactor containing ear motif protein2 (TIE2) forms a negative feedback loop to finely tune the activity of type B ARRs during root development. Disruption of TIE2 and its close homolog TIE1 causes severely shortened roots. TIE2 interacts with type B ARR1 and represses transcription of ARR1 targets. The cytokinin response is correspondingly enhanced in tie1-1 tie2-1. We further show that ARR1 positively regulates TIE1 and TIE2 by directly binding to their promoters. Our findings demonstrate that TIEs play key roles in controlling plant development and reveal an important negative feedback regulation mechanism for cytokinin signaling.

Sugar export from Arabidopsis leaves: actors and regulatory strategies

Plant acclimation and stress responses depend on the dynamic optimization of carbon balance between source and sink organs. This optimization also applies to the leaf export rate of photosynthetically produced sugars. So far, investigations into the molecular mechanisms of how the rate is controlled have focused on sugar transporters responsible for loading sucrose into the phloem sieve element-companion cell complex of leaf veins. Here, we take a broader view of the various proteins with potential direct influence on the leaf sugar export rate in the model plant Arabidopsis thaliana, helped by the cell type-specific transcriptome data that have recently become available. Furthermore, we integrate current information on the regulation of these potential target proteins. Our analysis identifies putative control points and units of transcriptionally and post-transcriptionally co-regulated genes. Most notable is the potential regulatory unit of sucrose transporters (SUC2, SWEET11, SWEET12, and SUC4) and proton pumps (AHA3 and AVP1). Our analysis can guide future research aimed at understanding the regulatory network controlling leaf sugar export by providing starting points for characterizing regulatory strategies and identifying regulatory factors that link sugar export rate to the major signaling pathways.

Evolution and identification of DREB transcription factors in the wheat genome: modeling, docking and simulation of DREB proteins associated with salt stress

Soil salinity and the resulting salt stress it imposes on crop plants is a major problem for modern agriculture. Understanding how salt tolerance mechanisms in plants are regulated is therefore important. One regulatory mechanism is the APETALA2/Ethylene Responsive Factor (AP2/ERF) transcription factor family, including dehydration responsive element binding (DREB) transcription factors. By binding to DNA, specifically upstream of genes that play roles in salt tolerance pathways, DREB proteins upregulate expression of these genes. DREB in Triticum aestivum (wheat) cluster in sub-groups and in this study by scanning the recently extended predicted proteome of wheat for DREB, we increased the number of members of this sub-family. Using the wheat genome, we identified 576 genes coding for the AP2 domain of which 508 were identified to have one AP2 domain, a characteristic of the DREB/ERF subfamily. We confirmed the existing four sub-groups by sequence-based phylogenetic analyses but also identified 32 new DREB subfamily members, not belonging to any known sub-group. Transcription factor profile inference analysis identified two genes, TraesCS2B02G002700 and TraesCS2D02G015200, being homologous to DREB1A of Arabidopsis thaliana. Based on molecular simulation (25 ns) analysis, TraesCS2B02G002700 with a CCGAC motif was observed to interact very stably with DNA. In silico mutational analysis at the 19^th position in the DREB domain of TraesCS2B02G002700-DNA complex indicated this as a stable part for recognizing and forming interaction with DNA. Moreover, six target genes were predicted having an upstream CCGAC motif regulated by TraesCS2B02G002700. Our study provides an overall framework for exploring the transcription factors in plants and identifying e.g. potential salt stress target genes.Communicated by Ramaswamy H. Sarma.

Network Analysis of Gene Transcriptions of Arabidopsis thaliana in Spaceflight Microgravity

The transcriptomic datasets of the plant model organism Arabidopsis thaliana grown in the International Space Station provided by GeneLab have been mined to isolate the impact of spaceflight microgravity on gene expressions related to root growth. A set of computational tools is used to identify the hub genes that respond differently in spaceflight with controlled lighting compared to on the ground. These computational tools based on graph-theoretic approaches are used to infer gene regulatory networks from the transcriptomic datasets. The three main algorithms used for network analyses are Least Absolute Shrinkage and Selection Operator (LASSO), Pearson correlation, and the Hyperlink-Induced Topic Search (HITS) algorithm. Graph-based spectral analyses reveal distinct properties of the spaceflight microgravity networks for the Wassilewskija (WS), Columbia (Col)-0, and mutant phytochromeD (phyD) ecotypes. The set of hub genes that are significantly altered in spaceflight microgravity are mainly involved in cell wall synthesis, protein transport, response to auxin, stress responses, and catabolic processes. Network analysis highlights five important root growth-regulating hub genes that have the highest outdegree distribution in spaceflight microgravity networks. These concerned genes coding for proteins are identified from the Gene Regulatory Networks (GRNs) corresponding to spaceflight total light environment. Furthermore, network analysis uncovers genes that encode nucleotide-diphospho-sugar interconversion enzymes that have higher transcriptional regulation in spaceflight microgravity and are involved in cell wall biosynthesis.

Unusual DNA-binding properties of the Arabidopsis thaliana WRKY50 transcription factor at target gene promoters

WRKY50 from A. thaliana requires WT-boxes at target gene promoters for activation and binding. Based on the genome-wide prediction of WRKY50 target genes and the similarity of a WRKY50 binding site to WT-boxes in microbe-associated molecular pattern (MAMP)-responsive cis-regulatory modules (CRM), four WT-box containing CRMs from the promoter region of three WRKY50 target genes were investigated for their interaction with WRKY50. These target genes are DJ1E, WRKY30 and ATBBE4. Two of the four CRMs, one from DJ1E and one from WRKY30, were able to activate reporter gene expression in the presence of WRKY50. Activation requires the WT-boxes GGACTTTT, GGACTTTG from DJ1E and GGACTTTC from WRKY30. WRKY50 does not activate a second CRM from WRKY30 and the CRM from ATBBE4, both containing the WT-box TGACTTTT. In vitro gel-shift assays demonstrate WT-box-specific binding of the WRKY50 DNA-binding domain to all four CRMs. This work shows a high flexibility of WRKY50 binding site recognition beyond the classic W-box TTGACC/T.

Genome-wide binding analysis reveals that ANAC060 directly represses sugar-induced transcription of ABI5 in Arabidopsis

The sugar status of a plant acts as a signal affecting growth and development. The phenomenon by which high levels of sugars inhibit seedling establishment has been widely used to gain insight into sugar-signaling pathways. Natural allelic variation has been identified at the ANAC060 locus. The Arabidopsis Columbia ecotype produces a short ANAC060 protein without a transmembrane domain that is constitutively located to the nucleus, causing sugar insensitivity when overexpressed. In this study, we generated a genome-wide DNA-binding map of ANAC060 via chromatin immunoprecipitation sequencing using transgenic lines that express a functional ANAC060-GFP fusion protein in an anac060 background. A total of 3282 genes associated with ANAC060-binding sites were identified. These genes were enriched in biotic and abiotic stress responses, and the G-box binding motif was highly enriched in ANAC060-bound genomic regions. Expression microarray analysis resulted in the identification of 8350 genes whose activities were altered in the anac060 mutant and upon sugar treatment. Cluster analysis revealed that ANAC060 attenuates sugar-regulated gene expression. Direct target genes of ANAC060 included equivalent numbers of genes that were upregulated or downregulated by ANAC060. The various functions of these target genes indicate that ANAC060 has several functions. Our results demonstrate that ANAC060 directly binds to the promoter of ABI5 and represses the sugar-induced transcription of ABI5. Genetic data indicate that ABI5 is epistatic to ANAC060 in both sugar and abscisic acid responses.

PCSD: a plant chromatin state database

Genome-wide maps of chromatin states have become a powerful representation of genome annotation and regulatory activity. We collected public and in-house plant epigenomic data sets and applied a Hidden Markov Model to define chromatin states, which included 290 553 (36 chromatin states), 831 235 (38 chromatin states) and 3 936 844 (26 chromatin states) segments across the whole genome of Arabidopsis thaliana, Oryza sativa and Zea mays, respectively. We constructed a Plant Chromatin State Database (PCSD, http://systemsbiology.cau.edu.cn/chromstates) to integrate detailed information about chromatin states, including the features and distribution of states, segments in states and related genes with segments. The self-organization mapping (SOM) results for these different chromatin signatures and UCSC Genome Browser for visualization were also integrated into the PCSD database. We further provided differential SOM maps between two epigenetic marks for chromatin state comparison and custom tools for new data analysis. The segments and related genes in SOM maps can be searched and used for motif and GO analysis, respectively. In addition, multi-species integration can be used to discover conserved features at the epigenomic level. In summary, our PCSD database integrated the identified chromatin states with epigenetic features and may be beneficial for communities to discover causal functions hidden in plant chromatin.

AtSPX1-mediated transcriptional regulation during leaf senescence in Arabidopsis thaliana

Leaf senescence is the final stage of leaf growth, a highly coordinated and complicated process. Phosphorus as an essential macronutrient for plant growth is remobilized from senescing leaves to other vigorous parts of the plant. In this study, through data mining, we found some phosphate starvation induced genes such as AtSPX1, were significantly induced in aging leaves in Arabidopsis. We applied a reverse genetics approach to investigate the phenotypes of transgenic plants and mutant plants, and the results showed that the overexpression of AtSPX1 accelerated leaf senescence, suppressed Pi accumulation, promoted SA production and H₂O₂ levels in leaves, while the mutant lines of AtSPX1 showed slightly delayed leaf senescence. We conducted RNA-seq-based transcriptome analysis together with GO and GSEA enrichment analyses for transgenic vs. wild-type plants to elucidate the possible underlying regulatory mechanism. The 558 genes that were up-regulated in the overexpression plants 35S::AtSPX1/WT, were significantly enriched in the process of leaf senescence, Pi starvation responses and SA signaling pathways, as were the target genes of some transcription factors such as WRKYs and NACs. In a word, we characterized AtSPX1 as a key regulator, which mediated the crosstalks among leaf senescence, Pi starvation and SA signaling pathways in Arabidopsis thaliana.

AppleMDO: A Multi-Dimensional Omics Database for Apple Co-Expression Networks and Chromatin States

As an economically important crop, apple is one of the most cultivated fruit trees in temperate regions worldwide. Recently, a large number of high-quality transcriptomic and epigenomic datasets for apple were made available to the public, which could be helpful in inferring gene regulatory relationships and thus predicting gene function at the genome level. Through integration of the available apple genomic, transcriptomic, and epigenomic datasets, we constructed co-expression networks, identified functional modules, and predicted chromatin states. A total of 112 RNA-seq datasets were integrated to construct a global network and a conditional network (tissue-preferential network). Furthermore, a total of 1,076 functional modules with closely related gene sets were identified to assess the modularity of biological networks and further subjected to functional enrichment analysis. The results showed that the function of many modules was related to development, secondary metabolism, hormone response, and transcriptional regulation. Transcriptional regulation is closely related to epigenetic marks on chromatin. A total of 20 epigenomic datasets, which included ChIP-seq, DNase-seq, and DNA methylation analysis datasets, were integrated and used to classify chromatin states. Based on the ChromHMM algorithm, the genome was divided into 620,122 fragments, which were classified into 24 states according to the combination of epigenetic marks and enriched-feature regions. Finally, through the collaborative analysis of different omics datasets, the online database AppleMDO (http://bioinformatics.cau.edu.cn/AppleMDO/) was established for cross-referencing and the exploration of possible novel functions of apple genes. In addition, gene annotation information and functional support toolkits were also provided. Our database might be convenient for researchers to develop insights into the function of genes related to important agronomic traits and might serve as a reference for other fruit trees.

PlantEAR: Functional Analysis Platform for Plant EAR Motif-Containing Proteins

The Ethylene-responsive element binding factor-associated Amphiphilic Repression (EAR) motifs, which were initially identified in members of the Arabidopsis ethylene response factor (ERF) family, are transcriptional repression motifs in plants and are defined by the consensus sequence patterns of either LxLxL or DLNxxP. EAR motif-containing proteins can function as transcription repressors, thus interacting with co-repressors, such as TOPLESS and AtSAP18, affecting the structure of chromatin by histone modifications and thereby repressing gene transcription. EAR motif-containing proteins are highly conserved across diverse plant species and play important roles in hormone signal transduction, stress responses and development, but they have not been identified in most plants. In this study, we identified 20,542 EAR motif-containing proteins from 71 plant species based on a Hidden Markov Model and orthologous gene search, and then we constructed a functional analysis platform for plant EAR motif-containing proteins (PlantEAR, http://structuralbiology.cau.edu.cn/plantEAR) by integrating a variety of functional annotations and processed data. Several tools were provided as functional support for EAR motif-containing proteins, such as browse, search, co-expression and protein-protein interaction (PPI) network analysis as well as cis-element analysis and gene set enrichment analysis (GSEA). In addition, basing on the identified EAR motif-containing proteins, we also explored their distribution in various species and found that the numbers of EAR motif-containing proteins showed an increasing trend in evolution from algae to angiosperms.

Circadian clock components control daily growth activities by modulating cytokinin levels and cell division-associated gene expression in Populus trees

Trees are carbon dioxide sinks and major producers of terrestrial biomass with distinct seasonal growth patterns. Circadian clocks enable the coordination of physiological and biochemical temporal activities, optimally regulating multiple traits including growth. To dissect the clock's role in growth, we analysed Populus tremula × P. tremuloides trees with impaired clock function due to down-regulation of central clock components. late elongated hypocotyl (lhy-10) trees, in which expression of LHY1 and LHY2 is reduced by RNAi, have a short free-running period and show disrupted temporal regulation of gene expression and reduced growth, producing 30-40% less biomass than wild-type trees. Genes important in growth regulation were expressed with an earlier phase in lhy-10, and CYCLIN D3 expression was misaligned and arrhythmic. Levels of cytokinins were lower in lhy-10 trees, which also showed a change in the time of peak expression of genes associated with cell division and growth. However, auxin levels were not altered in lhy-10 trees, and the size of the lignification zone in the stem showed a relative increase. The reduced growth rate and anatomical features of lhy-10 trees were mainly caused by misregulation of cell division, which may have resulted from impaired clock function.

High-quality genome assembly of Capsella bursa-pastoris reveals asymmetry of regulatory elements at early stages of polyploid genome evolution

Polyploidization and subsequent sub- and neofunctionalization of duplicated genes represent a major mechanism of plant genome evolution. Capsella bursa-pastoris, a widespread ruderal plant, is a recent allotetraploid and, thus, is an ideal model organism for studying early changes following polyploidization. We constructed a high-quality assembly of C. bursa-pastoris genome and a transcriptome atlas covering a broad sample of organs and developmental stages (available online at http://travadb.org/browse/Species=Cbp). We demonstrate that expression of homeologs is mostly symmetric between subgenomes, and identify a set of homeolog pairs with discordant expression. Comparison of promoters within such pairs revealed emerging asymmetry of regulatory elements. Among them there are multiple binding sites for transcription factors controlling the regulation of photosynthesis and plant development by light (PIF3, HY5) and cold stress response (CBF). These results suggest that polyploidization in C. bursa-pastoris enhanced its plasticity of response to light and temperature, and allowed substantial expansion of its distribution range.

Genome-wide mapping of DNase I hypersensitive sites reveals chromatin accessibility changes in Arabidopsis euchromatin and heterochromatin regions under extended darkness

Light, as the energy source in photosynthesis, is essential for plant growth and development. Extended darkness causes dramatic gene expression changes. In this study, we applied DNase-seq (DNase I hypersensitive site sequencing) to study changes of chromatin accessibility in euchromatic and heterochromatic regions under extended darkness in Arabidopsis. We generated 27 Gb DNase-seq and 67.6 Gb RNA-seq data to investigate chromatin accessibility changes and global gene expression under extended darkness and control condition in Arabidopsis. We found that ~40% DHSs (DNaseI hypersensitive sites) were diminished under darkness. In non-TE regions, the majority of DHS-changed genes were DHS-diminished under darkness. A total of 519 down-regulated genes were associated with diminished DHSs under darkness, mainly involved in photosynthesis process and retrograde signaling, and were regulated by chloroplast maintenance master regulators such as GLK1. In TE regions, approximately half of the DHS-changed TEs were DHS-increased under darkness and were primarily associated with the LTR/Gypsy retrotransposons in the heterochromatin flanking the centromeres. In contrast, DHS-diminished TEs under darkness were enriched in Copia, LINE, and MuDR dispersed across chromosomes. Together, our results indicated that extended darkness resulted in more increased chromatin compaction in euchromatin and decompaction in heterochromatin, thus further leading to gene expression changes in Arabidopsis.

Co-expression network analyses identify functional modules associated with development and stress response in Gossypium arboreum

Cotton is an economically important crop, essential for the agriculture and textile industries. Through integrating transcriptomic data, we discovered that multi-dimensional co-expression network analysis was powerful for predicting cotton gene functions and functional modules. Here, the recently available transcriptomic data on Gossypium arboreum, including data on multiple growth stages of tissues and stress treatment samples were applied to construct a co-expression network exploring multi-dimensional expression (development and stress) through multi-layered approaches. Based on differential gene expression and network analysis, a fibre development regulatory module of the gene GaKNL1 was found to regulate the second cell wall through repressing the activity of REVOLUTA, and a tissue-selective module of GaJAZ1a was examined in response to water stress. Moreover, comparative genomics analysis of the JAZ1-related regulatory module revealed high conservation across plant species. In addition, 1155 functional modules were identified through integrating the co-expression network, module classification and function enrichment tools, which cover functions such as metabolism, stress responses, and transcriptional regulation. In the end, an online platform was built for network analysis (http://structuralbiology.cau.edu.cn/arboreum), which could help to refine the annotation of cotton gene function and establish a data mining system to identify functional genes or modules with important agronomic traits.

ccNET: Database of co-expression networks with functional modules for diploid and polyploid Gossypium

Plant genera with both diploid and polyploid species are a common evolutionary occurrence. Polyploids, especially allopolyploids such as cotton and wheat, are a great model system for heterosis research. Here, we have integrated genome sequences and transcriptome data of Gossypium species to construct co-expression networks and identified functional modules from different cotton species, including 1155 and 1884 modules in G. arboreum and G. hirsutum, respectively. We overlayed the gene expression results onto the co-expression network. We further provided network comparison analysis for orthologous genes across the diploid and allotetraploid Gossypium. We also constructed miRNA-target networks and predicted PPI networks for both cotton species. Furthermore, we integrated in-house ChIP-seq data of histone modification (H3K4me3) together with cis-element analysis and gene sets enrichment analysis tools for studying possible gene regulatory mechanism in Gossypium species. Finally, we have constructed an online ccNET database (http://structuralbiology.cau.edu.cn/gossypium) for comparative gene functional analyses at a multi-dimensional network and epigenomic level across diploid and polyploid Gossypium species. The ccNET database will be beneficial for community to yield novel insights into gene/module functions during cotton development and stress response, and might be useful for studying conservation and diversity in other polyploid plants, such as T. aestivum and Brassica napus.

Histone modifications facilitate the coexpression of bidirectional promoters in rice

Background

Bidirectional gene pairs are highly abundant and mostly co-regulated in eukaryotic genomes. The structural features of bidirectional promoters (BDPs) have been well studied in yeast, humans and plants. However, the underlying mechanisms responsible for the coexpression of BDPs remain understudied, especially in plants.

Results

Here, we characterized chromatin features associated with rice BDPs. Several unique chromatin features were present in rice BDPs but were missing from unidirectional promoters (UDPs), including overrepresented active histone marks, canonical nucleosomes and underrepresented H3K27me3. In particular, overrepresented active marks (H3K4ac, H4K12ac, H4K16ac, H3K4me2 and H3K36me3) were truly overrepresented in type I BDPs but not in the other two BDPs, based on a Kolmogorov-Smirnov test.

Conclusions

Our analyses indicate that active marks (H3K4ac, H4K12ac, H4K16ac, H3K4me3, H3K9ac and H3K27ac) may coordinate with repressive marks (H3K27me3 and H3K9me1/3) to build a unique chromatin structure that favors the coregulation of bidirectional gene pairs. Thus, our findings help to enhance the understanding of unique epigenetic mechanisms that regulate bidirectional gene pairs and may improve the manipulation of gene pairs for crop bioengineering.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3125-0) contains supplementary material, which is available to authorized users.

SorghumFDB: sorghum functional genomics database with multidimensional network analysis

Sorghum (Sorghum bicolor [L.] Moench) has excellent agronomic traits and biological properties, such as heat and drought-tolerance. It is a C₄ grass and potential bioenergy-producing plant, which makes it an important crop worldwide. With the sorghum genome sequence released, it is essential to establish a sorghum functional genomics data mining platform. We collected genomic data and some functional annotations to construct a sorghum functional genomics database (SorghumFDB). SorghumFDB integrated knowledge of sorghum gene family classifications (transcription regulators/factors, carbohydrate-active enzymes, protein kinases, ubiquitins, cytochrome P450, monolignol biosynthesis related enzymes, R-genes and organelle-genes), detailed gene annotations, miRNA and target gene information, orthologous pairs in the model plants Arabidopsis, rice and maize, gene loci conversions and a genome browser. We further constructed a dynamic network of multidimensional biological relationships, comprised of the co-expression data, protein–protein interactions and miRNA-target pairs. We took effective measures to combine the network, gene set enrichment and motif analyses to determine the key regulators that participate in related metabolic pathways, such as the lignin pathway, which is a major biological process in bioenergy-producing plants.

Database URL:http://structuralbiology.cau.edu.cn/sorghum/index.html.

AtWRKY22 promotes susceptibility to aphids and modulates salicylic acid and jasmonic acid signalling

Aphids induce many transcriptional perturbations in their host plants, but the signalling cascades responsible and the effects on plant resistance are largely unknown. Through a genome-wide association (GWA) mapping study in Arabidopsis thaliana, we identified WRKY22 as a candidate gene associated with feeding behaviour of the green peach aphid, Myzus persicae The transcription factor WRKY22 is known to be involved in pathogen-triggered immunity, and WRKY22 gene expression has been shown to be induced by aphids. Assessment of aphid population development and feeding behaviour on knockout mutants and overexpression lines showed that WRKY22 increases susceptibility to M. persicae via a mesophyll-located mechanism. mRNA sequencing analysis of aphid-infested wrky22 knockout plants revealed the up-regulation of genes involved in salicylic acid (SA) signalling and down-regulation of genes involved in plant growth and cell-wall loosening. In addition, mechanostimulation of knockout plants by clip cages up-regulated jasmonic acid (JA)-responsive genes, resulting in substantial negative JA-SA crosstalk. Based on this and previous studies, WRKY22 is considered to modulate the interplay between the SA and JA pathways in response to a wide range of biotic and abiotic stimuli. Its induction by aphids and its role in suppressing SA and JA signalling make WRKY22 a potential target for aphids to manipulate host plant defences.

Boosting AthaMap Database Content with Data from Protein Binding Microarrays

The AthaMap database generates a map of predicted transcription factor binding sites (TFBS) and small RNA target sites for the whole Arabidopsis thaliana genome. With the advent of protein binding microarrays (PBM), the number of known TFBS for A. thaliana transcription factors (TFs) has increased dramatically. Using 113 positional weight matrices (PWMs) generated from a single PBM study and representing a total number of 68 different TFs, the number of predicted TFBS in AthaMap was now increased by about 3.8 × 10(7) to 4.9 × 10(7). The number of TFs with PWM-predicted TFBS annotated in AthaMap has increased to 126, representing a total of 29 TF families and newly including ARF, AT-Hook, YABBY, LOB/AS2 and SRS. Furthermore, links from all Arabidopsis TFs and genes to the newly established Arabidopsis Information Portal (AIP) have been implemented. With this qualitative and quantitative update, the improved AthaMap increases its value as a resource for the analysis of A. thaliana gene expression regulation at www.athamap.de.

PlantPAN 2.0: an update of plant promoter analysis navigator for reconstructing transcriptional regulatory networks in plants

Transcription factors (TFs) are sequence-specific DNA-binding proteins acting as critical regulators of gene expression. The Plant Promoter Analysis Navigator (PlantPAN; http://PlantPAN2.itps.ncku.edu.tw) provides an informative resource for detecting transcription factor binding sites (TFBSs), corresponding TFs, and other important regulatory elements (CpG islands and tandem repeats) in a promoter or a set of plant promoters. Additionally, TFBSs, CpG islands, and tandem repeats in the conserve regions between similar gene promoters are also identified. The current PlantPAN release (version 2.0) contains 16 960 TFs and 1143 TF binding site matrices among 76 plant species. In addition to updating of the annotation information, adding experimentally verified TF matrices, and making improvements in the visualization of transcriptional regulatory networks, several new features and functions are incorporated. These features include: (i) comprehensive curation of TF information (response conditions, target genes, and sequence logos of binding motifs, etc.), (ii) co-expression profiles of TFs and their target genes under various conditions, (iii) protein-protein interactions among TFs and their co-factors, (iv) TF-target networks, and (v) downstream promoter elements. Furthermore, a dynamic transcriptional regulatory network under various conditions is provided in PlantPAN 2.0. The PlantPAN 2.0 is a systematic platform for plant promoter analysis and reconstructing transcriptional regulatory networks.

How to Isolate a Plant's Hypomethylome in One Shot

Genome assembly remains a challenge for large and/or complex plant genomes due to their abundant repetitive regions resulting in studies focusing on gene space instead of the whole genome. Thus, DNA enrichment strategies facilitate the assembly by increasing the coverage and simultaneously reducing the complexity of the whole genome. In this paper we provide an easy, fast, and cost-effective variant of MRE-seq to obtain a plant's hypomethylome by an optimized methyl filtration protocol followed by next generation sequencing. The method is demonstrated on three plant species with knowingly large and/or complex (polyploid) genomes: Oryza sativa, Picea abies, and Crocus sativus. The identified hypomethylomes show clear enrichment for genes and their flanking regions and clear reduction of transposable elements. Additionally, genomic sequences around genes are captured including regulatory elements in introns and up- and downstream flanks. High similarity of the results obtained by a de novo assembly approach with a reference based mapping in rice supports the applicability for studying and understanding the genomes of nonmodel organisms. Hence we show the high potential of MRE-seq in a wide range of scenarios for the direct analysis of methylation differences, for example, between ecotypes, individuals, within or across species harbouring large, and complex genomes.

Expression of the floral repressor miRNA156 is positively regulated by the AGAMOUS-like proteins AGL15 and AGL18

The regulation of flowering time has crucial implications for plant fitness. MicroRNA156 (miR156) represses the floral transition in Arabidopsis thaliana, but the mechanisms regulating its transcription remain unclear. Here, we show that two AGAMOUS-like proteins, AGL15 and AGL18, act as positive regulators of the expression of MIR156. Small RNA northern blot analysis revealed a significant decrease in the levels of mature miR156 in agl15 agl18 double mutants, but not in the single mutants, suggesting that AGL15 and AGL18 co-regulate miR156 expression. Histochemical analysis further indicated that the double mutants showed a reduction in MIR156 promoter strength. The double mutants also showed reduced abundance of pri-miR156a and pri-miR156c, two of the primary transcripts from MIR156 genes. Electrophoretic mobility shift assays demonstrated that AGL15 directly associated with the CArG motifs in the MIR156a/c promoters. AGL18 did not show binding affinity to the CArG motifs, but pull-down and yeast two-hybrid assays showed that AGL18 forms a heterodimer with AGL15. GFP reporter assays and bimolecular fluorescence complementation (BiFC) showed that AGL15 and AGL18 co-localize in the nucleus and confirmed their in vivo interaction. Overexpression of miR156 did not affect the levels of AGL15 and AGL18 transcripts. Taking these data together, we present a model for the transcriptional regulation of MIR156. In this model, AGL15 and AGL18 may form a complex along with other proteins, and bind to the CArG motifs of the promoters of MIR156 to activate the MIR156 expression.

Evolution of proline biosynthesis: enzymology, bioinformatics, genetics, and transcriptional regulation

Proline is not only an essential component of proteins but it also has important roles in adaptation to osmotic and dehydration stresses, redox control, and apoptosis. Here, we review pathways of proline biosynthesis in the three domains of life. Pathway reconstruction from genome data for hundreds of eubacterial and dozens of archaeal and eukaryotic organisms revealed evolutionary conservation and variations of this pathway across different taxa. In the most prevalent pathway of proline synthesis, glutamate is phosphorylated to γ-glutamyl phosphate by γ-glutamyl kinase, reduced to γ-glutamyl semialdehyde by γ-glutamyl phosphate reductase, cyclized spontaneously to Δ(1)-pyrroline-5-carboxylate and reduced to proline by Δ(1)-pyrroline-5-carboxylate reductase. In higher plants and animals the first two steps are catalysed by a bi-functional Δ(1) -pyrroline-5-carboxylate synthase. Alternative pathways of proline formation use the initial steps of the arginine biosynthetic pathway to ornithine, which can be converted to Δ(1)-pyrroline-5-carboxylate by ornithine aminotransferase and then reduced to proline or converted directly to proline by ornithine cyclodeaminase. In some organisms, the latter pathways contribute to or could be fully responsible for the synthesis of proline. The conservation of proline biosynthetic enzymes and significance of specific residues for catalytic activity and allosteric regulation are analysed on the basis of protein structural data, multiple sequence alignments, and mutant studies, providing novel insights into proline biosynthesis in organisms. We also discuss the transcriptional control of the proline biosynthetic genes in bacteria and plants.