Publisher Correction: Extensive signal integration by the phytohormone protein network.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Extensive signal integration by the phytohormone protein network.

Plant hormones coordinate responses to environmental cues with developmental programs1, and are fundamental for stress resilience and agronomic yield2. The core signalling pathways underlying the effects of phytohormones have been elucidated by genetic screens and hypothesis-driven approaches, and extended by interactome studies of select pathways3. However, fundamental questions remain about how information from different pathways is integrated. Genetically, most phenotypes seem to be regulated by several hormones, but transcriptional profiling suggests that hormones trigger largely exclusive transcriptional programs4. We hypothesized that protein-protein interactions have an important role in phytohormone signal integration. Here, we experimentally generated a systems-level map of the Arabidopsis phytohormone signalling network, consisting of more than 2,000 binary protein-protein interactions. In the highly interconnected network, we identify pathway communities and hundreds of previously unknown pathway contacts that represent potential points of crosstalk. Functional validation of candidates in seven hormone pathways reveals new functions for 74% of tested proteins in 84% of candidate interactions, and indicates that a large majority of signalling proteins function pleiotropically in several pathways. Moreover, we identify several hundred largely small-molecule-dependent interactions of hormone receptors. Comparison with previous reports suggests that noncanonical and nontranscription-mediated receptor signalling is more common than hitherto appreciated.

The Systems Architecture of Molecular Memory in Poplar after Abiotic Stress.

Throughout the temperate zones, plants face combined drought and heat spells in increasing frequency and intensity. Here, we compared periodic (intermittent, i.e., high-frequency) versus chronic (continuous, i.e., high-intensity) drought-heat stress scenarios in gray poplar (Populus× canescens) plants for phenotypic and transcriptomic effects during stress and after recovery. Photosynthetic productivity after stress recovery exceeded the performance of poplar trees without stress experience. We analyzed the molecular basis of this stress-related memory phenotype and investigated gene expression responses across five major tree compartments including organs and wood tissues. For each of these tissue samples, transcriptomic changes induced by the two stress scenarios were highly similar during the stress phase but strikingly divergent after recovery. Characteristic molecular response patterns were found across tissues but involved different genes in each tissue. Only a small fraction of genes showed similar stress and recovery expression profiles across all tissues, including type 2C protein phosphatases, the LATE EMBRYOGENESIS ABUNDANT PROTEIN4-5 genes, and homologs of the Arabidopsis (Arabidopsis thaliana) transcription factor HOMEOBOX7. Analysis of the predicted transcription factor regulatory networks for these genes suggested that a complex interplay of common and tissue-specific components contributes to the coordination of post-recovery responses to stress in woody plants.

Barley susceptibility factor RACB modulates transcript levels of signalling protein genes in compatible interaction with Blumeria graminis f.sp. hordei.

RHO (rat sarcoma homologue) GTPases (guanosine triphosphatases) are regulators of downstream transcriptional responses of eukaryotes to intracellular and extracellular stimuli. For plants, little is known about the function of Rho-like GTPases [called RACs (rat sarcoma-related C botulinum substrate) or ROPs (RHO of plants)] in transcriptional reprogramming of cells. However, in plant hormone response and innate immunity, RAC/ROP proteins influence gene expression patterns. The barley RAC/ROP RACB is required for full susceptibility of barley to the powdery mildew fungus Blumeria graminis f.sp. hordei (Bgh). We compared the transcriptomes of barley plants either silenced for RACB or over-expressing constitutively activated RACB with and without inoculation with Bgh. This revealed a large overlap of the barley transcriptome during the early response to Bgh and during the over-expression of constitutively activated RACB. Global pathway analyses and stringent analyses of differentially expressed genes suggested that RACB influences, amongst others, the expression of signalling receptor kinases. Transient induced gene silencing of RACB-regulated signalling genes (a leucine-rich repeat protein, a leucine-rich repeat receptor-like kinase and an S-domain SD1-receptor-like kinase) suggested that they might be involved in RACB-modulated susceptibility to powdery mildew. We discuss the function of RACB in regulating the transcriptional responses of susceptible barley to Bgh.

Genome sequence of the progenitor of the wheat D genome Aegilops tauschii.

Aegilops tauschii is the diploid progenitor of the D genome of hexaploid wheat (Triticum aestivum, genomes AABBDD) and an important genetic resource for wheat. The large size and highly repetitive nature of the Ae. tauschii genome has until now precluded the development of a reference-quality genome sequence. Here we use an array of advanced technologies, including ordered-clone genome sequencing, whole-genome shotgun sequencing, and BioNano optical genome mapping, to generate a reference-quality genome sequence for Ae. tauschii ssp. strangulata accession AL8/78, which is closely related to the wheat D genome. We show that compared to other sequenced plant genomes, including a much larger conifer genome, the Ae. tauschii genome contains unprecedented amounts of very similar repeated sequences. Our genome comparisons reveal that the Ae. tauschii genome has a greater number of dispersed duplicated genes than other sequenced genomes and its chromosomes have been structurally evolving an order of magnitude faster than those of other grass genomes. The decay of colinearity with other grass genomes correlates with recombination rates along chromosomes. We propose that the vast amounts of very similar repeated sequences cause frequent errors in recombination and lead to gene duplications and structural chromosome changes that drive fast genome evolution.

Examining the Transcriptional Response in Wheat Near-Isogenic Lines to Infection and Deoxynivalenol Treatment.

head blight (FHB) is a disease caused predominantly by the fungal pathogen that affects wheat and other small-grain cereals and can lead to severe yield loss and reduction in grain quality. Trichothecene mycotoxins, such as deoxynivalenol (DON), accumulate during infection and increase pathogen virulence and decrease grain quality. The locus on wheat chromosome 3BS confers Type II resistance to FHB and resistance to the spread of infection on the spike and has been associated with resistance to DON accumulation. To gain a better genetic understanding of the functional role of and resistance or susceptibility to FHB, we examined DON and ergosterol accumulation, FHB resistance, and the whole-genome transcriptomic response using RNA-seq in a near-isogenic line (NIL) pair carrying the resistant and susceptible alleles for during infection and DON treatment. Our results provide a gene expression atlas for the resistant and susceptible wheat- interaction. The DON concentration and transcriptomic results show that the rachis is a key location for conferring Type II resistance. In addition, the wheat transcriptome analysis revealed a set of -responsive genes that may play a role in resistance and a set of DON-responsive genes that may play a role in trichothecene resistance. Transcriptomic results from the pathogen show that the genome responds differently to the host level of resistance. The results of this study extend our understanding of host and pathogen responses in the wheat- interaction.

Brassinosteroids participate in the control of basal and acquired freezing tolerance of plants.

Brassinosteroids (BRs) are growth-promoting plant hormones that play a role in abiotic stress responses, but molecular modes that enable this activity remain largely unknown. Here we show that BRs participate in the regulation of freezing tolerance. BR signaling-defective mutants of Arabidopsis thaliana were hypersensitive to freezing before and after cold acclimation. The constitutive activation of BR signaling, in contrast, enhanced freezing resistance. Evidence is provided that the BR-controlled basic helix-loop-helix transcription factor CESTA (CES) can contribute to the constitutive expression of the C-REPEAT/DEHYDRATION-RESPONSIVE ELEMENT BINDING FACTOR (CBF) transcriptional regulators that control cold responsive (COR) gene expression. In addition, CBF-independent classes of BR-regulated COR genes are identified that are regulated in a BR- and CES-dependent manner during cold acclimation. A model is presented in which BRs govern different cold-responsive transcriptional cascades through the posttranslational modification of CES and redundantly acting factors. This contributes to the basal resistance against freezing stress, but also to the further improvement of this resistance through cold acclimation.

Ribosome quality control is a central protection mechanism for yeast exposed to deoxynivalenol and trichothecin.

BACKGROUND: The trichothecene mycotoxins deoxynivalenol (DON) and trichothecin (TTC) are inhibitors of eukaryotic protein synthesis. Their effect on cellular homeostasis is poorly understood. We report a systematic functional investigation of the effect of DON and TTC on the yeast Saccharomyces cerevisiae using genetic array, network and microarray analysis. To focus the genetic analysis on intracellular consequences of toxin action we eliminated the PDR5 gene coding for a potent pleiotropic drug efflux protein potentially confounding results. We therefore used a knockout library with a pdr5Δ strain background. RESULTS: DON or TTC treatment creates a fitness bottleneck connected to ribosome efficiency. Genes isolated by systematic genetic array analysis as contributing to toxin resistance function in ribosome quality control, translation fidelity, and in transcription. Mutants in the E3 ligase Hel2, involved in ribosome quality control, and several members of the Rpd3 histone deacetylase complex were highly sensitive to DON. DON and TTC have similar genetic profiles despite their different toxic potency. Network analysis shows a coherent and tight network of genetic interactions among the DON and TTC resistance conferring gene products. The networks exhibited topological properties commonly associated with efficient processing of information. Many sensitive mutants have a "slow growth" gene expression signature. DON-exposed yeast cells increase transcripts of ribosomal protein and histone genes indicating an internal signal for growth enhancement. CONCLUSIONS: The combination of gene expression profiling and analysis of mutants reveals cellular pathways which become bottlenecks under DON and TTC stress. These are generally directly or indirectly connected to ribosome biosynthesis such as the general secretory pathway, cytoskeleton, cell cycle delay, ribosome synthesis and translation quality control. Gene expression profiling points to an increased demand of ribosomal components and does not reveal activation of stress pathways. Our analysis highlights ribosome quality control and a contribution of a histone deacetylase complex as main sources of resistance against DON and TTC.

A protocol for the systematic and quantitative measurement of protein-lipid interactions using the liposome-microarray-based assay.

Lipids organize the activity of the cell's proteome through a complex network of interactions. The assembly of comprehensive atlases embracing all protein-lipid interactions is an important challenge that requires innovative methods. We recently developed a liposome-microarray-based assay (LiMA) that integrates liposomes, microfluidics and fluorescence microscopy and which is capable of measuring protein recruitment to membranes in a quantitative and high-throughput manner. Compared with previous assays that are labor-intensive and difficult to scale up, LiMA improves the protein-lipid interaction assay throughput by at least three orders of magnitude. Here we provide a step-by-step LiMA protocol that includes the following: (i) the serial and generic production of the liposome microarray; (ii) its integration into a microfluidic format; (iii) the measurement of fluorescently labeled protein (either purified proteins or from cell lysate) recruitment to liposomal membranes using high-throughput microscopy; (iv) automated image analysis pipelines to quantify protein-lipid interactions; and (v) data quality analysis. In addition, we discuss the experimental design, including the relevant quality controls. Overall, the protocol-including device preparation, assay and data analysis-takes 6-8 d. This protocol paves the way for protein-lipid interaction screens to be performed on the proteome and lipidome scales.

PGSB PlantsDB: updates to the database framework for comparative plant genome research.

PGSB (Plant Genome and Systems Biology: formerly MIPS) PlantsDB (http://pgsb.helmholtz-muenchen.de/plant/index.jsp) is a database framework for the comparative analysis and visualization of plant genome data. The resource has been updated with new data sets and types as well as specialized tools and interfaces to address user demands for intuitive access to complex plant genome data. In its latest incarnation, we have re-worked both the layout and navigation structure and implemented new keyword search options and a new BLAST sequence search functionality. Actively involved in corresponding sequencing consortia, PlantsDB has dedicated special efforts to the integration and visualization of complex triticeae genome data, especially for barley, wheat and rye. We enhanced CrowsNest, a tool to visualize syntenic relationships between genomes, with data from the wheat sub-genome progenitor Aegilops tauschii and added functionality to the PGSB RNASeqExpressionBrowser. GenomeZipper results were integrated for the genomes of barley, rye, wheat and perennial ryegrass and interactive access is granted through PlantsDB interfaces. Data exchange and cross-linking between PlantsDB and other plant genome databases is stimulated by the transPLANT project (http://transplantdb.eu/).

Identification and Characterization of Carboxylesterases from Brachypodium distachyon Deacetylating Trichothecene Mycotoxins.

Increasing frequencies of 3-acetyl-deoxynivalenol (3-ADON)-producing strains of Fusarium graminearum (3-ADON chemotype) have been reported in North America and Asia. 3-ADON is nearly nontoxic at the level of the ribosomal target and has to be deacetylated to cause inhibition of protein biosynthesis. Plant cells can efficiently remove the acetyl groups of 3-ADON, but the underlying genes are yet unknown. We therefore performed a study of the family of candidate carboxylesterases (CXE) genes of the monocot model plant Brachypodium distachyon. We report the identification and characterization of the first plant enzymes responsible for deacetylation of trichothecene toxins. The product of the BdCXE29 gene efficiently deacetylates T-2 toxin to HT-2 toxin, NX-2 to NX-3, both 3-ADON and 15-acetyl-deoxynivalenol (15-ADON) into deoxynivalenol and, to a lesser degree, also fusarenon X into nivalenol. The BdCXE52 esterase showed lower activity than BdCXE29 when expressed in yeast and accepts 3-ADON, NX-2, 15-ADON and, to a limited extent, fusarenon X as substrates. Expression of these Brachypodium genes in yeast increases the toxicity of 3-ADON, suggesting that highly similar genes existing in crop plants may act as susceptibility factors in Fusarium head blight disease.

Joint Transcriptomic and Metabolomic Analyses Reveal Changes in the Primary Metabolism and Imbalances in the Subgenome Orchestration in the Bread Wheat Molecular Response to Fusarium graminearum.

Fusarium head blight is a prevalent disease of bread wheat (Triticum aestivum L.), which leads to considerable losses in yield and quality. Quantitative resistance to the causative fungus Fusarium graminearum is poorly understood. We integrated transcriptomics and metabolomics data to dissect the molecular response to the fungus and its main virulence factor, the toxin deoxynivalenol in near-isogenic lines segregating for two resistance quantitative trait loci, Fhb1 and Qfhs.ifa-5A. The data sets portrait rearrangements in the primary metabolism and the translational machinery to counter the fungus and the effects of the toxin and highlight distinct changes in the metabolism of glutamate in lines carrying Qfhs.ifa-5A. These observations are possibly due to the activity of two amino acid permeases located in the quantitative trait locus confidence interval, which may contribute to increased pathogen endurance. Mapping to the highly resolved region of Fhb1 reduced the list of candidates to few genes that are specifically expressed in presence of the quantitative trait loci and in response to the pathogen, which include a receptor-like protein kinase, a protein kinase, and an E3 ubiquitin-protein ligase. On a genome-scale level, the individual subgenomes of hexaploid wheat contribute differentially to defense. In particular, the D subgenome exhibited a pronounced response to the pathogen and contributed significantly to the overall defense response.

Lipid Cooperativity as a General Membrane-Recruitment Principle for PH Domains.

Many cellular processes involve the recruitment of proteins to specific membranes, which are decorated with distinctive lipids that act as docking sites. The phosphoinositides form signaling hubs, and we examine mechanisms underlying recruitment. We applied a physiological, quantitative, liposome microarray-based assay to measure the membrane-binding properties of 91 pleckstrin homology (PH) domains, the most common phosphoinositide-binding target. 10,514 experiments quantified the role of phosphoinositides in membrane recruitment. For most domains examined, the observed binding specificity implied cooperativity with additional signaling lipids. Analyses of PH domains with similar lipid-binding profiles identified a conserved motif, mutations in which-including some found in human cancers-induced discrete changes in binding affinities in vitro and protein mislocalization in vivo. The data set reveals cooperativity as a key mechanism for membrane recruitment and, by enabling the interpretation of disease-associated mutations, suggests avenues for the design of small molecules targeting PH domains.

Brassinosteroids Are Master Regulators of Gibberellin Biosynthesis in Arabidopsis.

Plant growth and development are highly regulated processes that are coordinated by hormones including the brassinosteroids (BRs), a group of steroids with structural similarity to steroid hormones of mammals. Although it is well understood how BRs are produced and how their signals are transduced, BR targets, which directly confer the hormone's growth-promoting effects, have remained largely elusive. Here, we show that BRs regulate the biosynthesis of gibberellins (GAs), another class of growth-promoting hormones, in Arabidopsis thaliana. We reveal that Arabidopsis mutants deficient in BR signaling are severely impaired in the production of bioactive GA, which is correlated with defective GA biosynthetic gene expression. Expression of the key GA biosynthesis gene GA20ox1 in the BR signaling mutant bri1-301 rescues many of its developmental defects. We provide evidence that supports a model in which the BR-regulated transcription factor BES1 binds to a regulatory element in promoters of GA biosynthesis genes in a BR-induced manner to control their expression. In summary, our study underscores a role of BRs as master regulators of GA biosynthesis and shows that this function is of major relevance for the growth and development of vascular plants.

The deca-GX3 proteins Yae1-Lto1 function as adaptors recruiting the ABC protein Rli1 for iron-sulfur cluster insertion.

Cytosolic and nuclear iron-sulfur (Fe-S) proteins are involved in many essential pathways including translation and DNA maintenance. Their maturation requires the cytosolic Fe-S protein assembly (CIA) machinery. To identify new CIA proteins we employed systematic protein interaction approaches and discovered the essential proteins Yae1 and Lto1 as binding partners of the CIA targeting complex. Depletion of Yae1 or Lto1 results in defective Fe-S maturation of the ribosome-associated ABC protein Rli1, but surprisingly no other tested targets. Yae1 and Lto1 facilitate Fe-S cluster assembly on Rli1 in a chain of binding events. Lto1 uses its conserved C-terminal tryptophan for binding the CIA targeting complex, the deca-GX3 motifs in both Yae1 and Lto1 facilitate their complex formation, and Yae1 recruits Rli1. Human YAE1D1 and the cancer-related ORAOV1 can replace their yeast counterparts demonstrating evolutionary conservation. Collectively, the Yae1-Lto1 complex functions as a target-specific adaptor that recruits apo-Rli1 to the generic CIA machinery.

TopKLists: a comprehensive R package for statistical inference, stochastic aggregation, and visualization of multiple omics ranked lists.

High-throughput sequencing techniques are increasingly affordable and produce massive amounts of data. Together with other high-throughput technologies, such as microarrays, there are an enormous amount of resources in databases. The collection of these valuable data has been routine for more than a decade. Despite different technologies, many experiments share the same goal. For instance, the aims of RNA-seq studies often coincide with those of differential gene expression experiments based on microarrays. As such, it would be logical to utilize all available data. However, there is a lack of biostatistical tools for the integration of results obtained from different technologies. Although diverse technological platforms produce different raw data, one commonality for experiments with the same goal is that all the outcomes can be transformed into a platform-independent data format - rankings - for the same set of items. Here we present the R package TopKLists, which allows for statistical inference on the lengths of informative (top-k) partial lists, for stochastic aggregation of full or partial lists, and for graphical exploration of the input and consolidated output. A graphical user interface has also been implemented for providing access to the underlying algorithms. To illustrate the applicability and usefulness of the package, we integrated microRNA data of non-small cell lung cancer across different measurement techniques and draw conclusions. The package can be obtained from CRAN under a LGPL-3 license.

chromoWIZ: a web tool to query and visualize chromosome-anchored genes from cereal and model genomes.

BACKGROUND: Over the last years reference genome sequences of several economically and scientifically important cereals and model plants became available. Despite the agricultural significance of these crops only a small number of tools exist that allow users to inspect and visualize the genomic position of genes of interest in an interactive manner. DESCRIPTION: We present chromoWIZ, a web tool that allows visualizing the genomic positions of relevant genes and comparing these data between different plant genomes. Genes can be queried using gene identifiers, functional annotations, or sequence homology in four grass species (Triticum aestivum, Hordeum vulgare, Brachypodium distachyon, Oryza sativa). The distribution of the anchored genes is visualized along the chromosomes by using heat maps. Custom gene expression measurements, differential expression information, and gene-to-group mappings can be uploaded and can be used for further filtering. CONCLUSIONS: This tool is mainly designed for breeders and plant researchers, who are interested in the location and the distribution of candidate genes as well as in the syntenic relationships between different grass species. chromoWIZ is freely available and online accessible at http://mips.helmholtz-muenchen.de/plant/chromoWIZ/index.jsp.

Bacteria-triggered systemic immunity in barley is associated with WRKY and ETHYLENE RESPONSIVE FACTORs but not with salicylic acid.

Leaf-to-leaf systemic immune signaling known as systemic acquired resistance is poorly understood in monocotyledonous plants. Here, we characterize systemic immunity in barley (Hordeum vulgare) triggered after primary leaf infection with either Pseudomonas syringae pathovar japonica (Psj) or Xanthomonas translucens pathovar cerealis (Xtc). Both pathogens induced resistance in systemic, uninfected leaves against a subsequent challenge infection with Xtc. In contrast to systemic acquired resistance in Arabidopsis (Arabidopsis thaliana), systemic immunity in barley was not associated with NONEXPRESSOR OF PATHOGENESIS-RELATED GENES1 or the local or systemic accumulation of salicylic acid. Instead, we documented a moderate local but not systemic induction of abscisic acid after infection of leaves with Psj. In contrast to salicylic acid or its functional analog benzothiadiazole, local applications of the jasmonic acid methyl ester or abscisic acid triggered systemic immunity to Xtc. RNA sequencing analysis of local and systemic transcript accumulation revealed unique gene expression changes in response to both Psj and Xtc and a clear separation of local from systemic responses. The systemic response appeared relatively modest, and quantitative reverse transcription-polymerase chain reaction associated systemic immunity with the local and systemic induction of two WRKY and two ETHYLENE RESPONSIVE FACTOR (ERF)-like transcription factors. Systemic immunity against Xtc was further associated with transcriptional changes after a secondary/systemic Xtc challenge infection; these changes were dependent on the primary treatment. Taken together, bacteria-induced systemic immunity in barley may be mediated in part by WRKY and ERF-like transcription factors, possibly facilitating transcriptional reprogramming to potentiate immunity.

Genetic analysis of DEFECTIVE KERNEL1 loop function in three-dimensional body patterning in Physcomitrella patens.

DEFECTIVE KERNEL1 (DEK1) of higher plants plays an essential role in position-dependent signaling and consists of a large transmembrane domain (MEM) linked to a protease catalytic domain and a regulatory domain. Here, we show that the postulated sensory Loop of the MEM domain plays an important role in the developmental regulation of DEK1 activity in the moss Physcomitrella patens. Compared with P. patens lacking DEK1 (∆dek1), the dek1∆loop mutant correctly positions the division plane in the bud apical cell. In contrast with an early developmental arrest of ∆dek1 buds, dek1∆loop develops aberrant gametophores lacking expanded phyllids resulting from misregulation of mitotic activity. In contrast with the highly conserved sequence of the protease catalytic domain, the Loop is highly variable in land plants. Functionally, the sequence from Marchantia polymorpha fully complements the dek1∆loop phenotype, whereas sequences from maize (Zea mays) and Arabidopsis (Arabidopsis thaliana) give phenotypes with retarded growth and affected phyllid development. Bioinformatic analysis identifies MEM as a member of the Major Facilitator Superfamily, membrane transporters reacting to stimuli from the external environment. Transcriptome analysis comparing wild-type and ∆dek1 tissues identifies an effect on two groups of transcripts connected to dek1 mutant phenotypes: transcripts related to cell wall remodeling and regulation of the AINTEGUMENTA, PLETHORA, and BABY BOOM2 (APB2) and APB3 transcription factors known to regulate bud initiation. Finally, sequence data support the hypothesis that the advanced charophyte algae that evolved into ancestral land plants lost cytosolic calpains, retaining DEK1 as the sole calpain in the evolving land plant lineage.

Genome interplay in the grain transcriptome of hexaploid bread wheat.

Allohexaploid bread wheat (Triticum aestivum L.) provides approximately 20% of calories consumed by humans. Lack of genome sequence for the three homeologous and highly similar bread wheat genomes (A, B, and D) has impeded expression analysis of the grain transcriptome. We used previously unknown genome information to analyze the cell type-specific expression of homeologous genes in the developing wheat grain and identified distinct co-expression clusters reflecting the spatiotemporal progression during endosperm development. We observed no global but cell type- and stage-dependent genome dominance, organization of the wheat genome into transcriptionally active chromosomal regions, and asymmetric expression in gene families related to baking quality. Our findings give insight into the transcriptional dynamics and genome interplay among individual grain cell types in a polyploid cereal genome.