Arabidopsis defense against Botrytis cinerea: chronology and regulation deciphered by high-resolution temporal transcriptomic analysis

Negative regulation of ABA signaling by WRKY33 is critical for Arabidopsis immunity towards Botrytis cinerea 2100

The Arabidopsis mutant wrky33 is highly susceptible to Botrytis cinerea. We identified >1680 Botrytis-induced WRKY33 binding sites associated with 1576 Arabidopsis genes. Transcriptional profiling defined 318 functional direct target genes at 14 hr post inoculation. Comparative analyses revealed that WRKY33 possesses dual functionality acting either as a repressor or as an activator in a promoter-context dependent manner. We confirmed known WRKY33 targets involved in hormone signaling and phytoalexin biosynthesis, but also uncovered a novel negative role of abscisic acid (ABA) in resistance towards B. cinerea 2100. The ABA biosynthesis genes NCED3 and NCED5 were identified as direct targets required for WRKY33-mediated resistance. Loss-of-WRKY33 function resulted in elevated ABA levels and genetic studies confirmed that WRKY33 acts upstream of NCED3/NCED5 to negatively regulate ABA biosynthesis. This study provides the first detailed view of the genome-wide contribution of a specific plant transcription factor in modulating the transcriptional network associated with plant immunity.

Negative regulation of ABA signaling by WRKY33 is critical for Arabidopsis immunity towards Botrytis cinerea 2100

The Arabidopsis mutant wrky33 is highly susceptible to Botrytis cinerea. We identified >1680 Botrytis-induced WRKY33 binding sites associated with 1576 Arabidopsis genes. Transcriptional profiling defined 318 functional direct target genes at 14 hr post inoculation. Comparative analyses revealed that WRKY33 possesses dual functionality acting either as a repressor or as an activator in a promoter-context dependent manner. We confirmed known WRKY33 targets involved in hormone signaling and phytoalexin biosynthesis, but also uncovered a novel negative role of abscisic acid (ABA) in resistance towards B. cinerea 2100. The ABA biosynthesis genes NCED3 and NCED5 were identified as direct targets required for WRKY33-mediated resistance. Loss-of-WRKY33 function resulted in elevated ABA levels and genetic studies confirmed that WRKY33 acts upstream of NCED3/NCED5 to negatively regulate ABA biosynthesis. This study provides the first detailed view of the genome-wide contribution of a specific plant transcription factor in modulating the transcriptional network associated with plant immunity.

Chloroplasts play a central role in plant defence and are targeted by pathogen effectors

Microbe associated molecular pattern (MAMP) receptors in plants recognize MAMPs and activate basal defences; however a complete understanding of the molecular and physiological mechanisms conferring immunity remains elusive. Pathogens suppress active defence in plants through the combined action of effector proteins. Here we show that the chloroplast is a key component of early immune responses. MAMP perception triggers the rapid, large-scale suppression of nuclear encoded chloroplast-targeted genes (NECGs). Virulent Pseudomonas syringae effectors reprogramme NECG expression in Arabidopsis, target the chloroplast and inhibit photosynthetic CO2 assimilation through disruption of photosystem II. This activity prevents a chloroplastic reactive oxygen burst. These physiological changes precede bacterial multiplication and coincide with pathogen-induced abscisic acid (ABA) accumulation. MAMP pretreatment protects chloroplasts from effector manipulation, whereas application of ABA or the inhibitor of photosynthetic electron transport, DCMU, abolishes the MAMP-induced chloroplastic reactive oxygen burst, and enhances growth of a P. syringae hrpA mutant that fails to secrete effectors.

The toolbox of Trichoderma spp. in the biocontrol of Botrytis cinerea disease

Botrytis cinerea is a necrotrophic fungal pathogen causing disease in many plant species, leading to economically important crop losses. So far, fungicides have been widely used to control this pathogen. However, in addition to their detrimental effects on the environment and potential risks for human health, increasing fungicide resistance has been observed in the B. cinerea population. Biological control, that is the application of microbial organisms to reduce disease, has gained importance as an alternative or complementary approach to fungicides. In this respect, the genus Trichoderma constitutes a promising pool of organisms with potential for B. cinerea control. In the first part of this article, we review the specific mechanisms involved in the direct interaction between the two fungi, including mycoparasitism, the production of antimicrobial compounds and enzymes (collectively called antagonism), and competition for nutrients and space. In addition, biocontrol has also been observed when Trichoderma is physically separated from the pathogen, thus implying an indirect systemic plant defence response. Therefore, in the second part, we describe the consecutive steps leading to induced systemic resistance (ISR), starting with the initial Trichoderma-plant interaction and followed by the activation of downstream signal transduction pathways and, ultimately, the defence response resulting in ISR (ISR-prime phase). Finally, we discuss the ISR-boost phase, representing the effect of ISR priming by Trichoderma spp. on plant responses after additional challenge with B. cinerea.

Functional studies of transcription factors involved in plant defenses in the genomics era

Plant transcription factors (TFs) play roles in diverse biological processes including defense responses to pathogens. Here, we provide an overview of recent studies of plant TFs with regard to defense responses. TFs play roles in plant innate immunity by regulating genes related to pathogen-associated molecular pattern-triggered immunity, effector-triggered immunity, hormone signaling pathways and phytoalexin synthesis. Currently, genome-wide phylogenetic and transcriptomic analyses are as important as functional analyses in the study of plant TFs. The integration of genomics information with the knowledge obtained from functional studies provides new insights into the regulation of plant defense mechanisms as well as engineering crops with improved resistance to invading pathogens.

Interplays between the cell wall and phytohormones in interaction between plants and necrotrophic pathogens

The plant cell wall surrounds every cell in plants. During microbial infection, the cell wall provides a dynamic interface for interaction with necrotrophic phytopathogens as a rich source of carbohydrates for the growth of pathogens, as a physical barrier restricting the progression of the pathogens, and as an integrity sensory system that can activate intracellular signaling cascades and ultimately lead to a multitude of inducible host defense responses. Studies over the last decade have provided evidence of interplays between the cell wall and phytohormone signaling. This review summarizes the current state of knowledge about the cell wall-phytohormone interplays, with the focus on auxin, cytokinin, brassinosteroids, and abscisic acid, and discuss how they impact the outcome of plant-necrotrophic pathogen interaction.

Transcriptional networks-crops, clocks, and abiotic stress

Several factors affect the yield potential and geographical range of crops including the circadian clock, water availability, and seasonal temperature changes. In order to sustain and increase plant productivity on marginal land in the face of both biotic and abiotic stresses, we need to more efficiently generate stress-resistant crops through marker-assisted breeding, genetic modification, and new genome-editing technologies. To leverage these strategies for producing the next generation of crops, future transcriptomic data acquisition should be pursued with an appropriate temporal design and analyzed with a network-centric approach. The following review focuses on recent developments in abiotic stress transcriptional networks in economically important crops and will highlight the utility of correlation-based network analysis and applications.

Reexamination of chlorophyllase function implies its involvement in defense against chewing herbivores

Chlorophyllase (CLH) is a common plant enzyme that catalyzes the hydrolysis of chlorophyll to form chlorophyllide, a more hydrophilic derivative. For more than a century, the biological role of CLH has been controversial, although this enzyme has been often considered to catalyze chlorophyll catabolism during stress-induced chlorophyll breakdown. In this study, we found that the absence of CLH does not affect chlorophyll breakdown in intact leaf tissue in the absence or the presence of methyl-jasmonate, which is known to enhance stress-induced chlorophyll breakdown. Fractionation of cellular membranes shows that Arabidopsis (Arabidopsis thaliana) CLH is located in the endoplasmic reticulum and the tonoplast of intact plant cells. These results indicate that CLH is not involved in endogenous chlorophyll catabolism. Instead, we found that CLH promotes chlorophyllide formation upon disruption of leaf cells, or when it is artificially mistargeted to the chloroplast. These results indicate that CLH is responsible for chlorophyllide formation after the collapse of cells, which led us to hypothesize that chlorophyllide formation might be a process of defense against chewing herbivores. We found that Arabidopsis leaves with genetically enhanced CLH activity exhibit toxicity when fed to Spodoptera litura larvae, an insect herbivore. In addition, purified chlorophyllide partially suppresses the growth of the larvae. Taken together, these results support the presence of a unique binary defense system against insect herbivores involving chlorophyll and CLH. Potential mechanisms of chlorophyllide action for defense are discussed.

Multidimensional approaches for studying plant defence against insects: from ecology to omics and synthetic biology

The biggest challenge for modern biology is to integrate multidisciplinary approaches towards understanding the organizational and functional complexity of biological systems at different hierarchies, starting from the subcellular molecular mechanisms (microscopic) to the functional interactions of ecological communities (macroscopic). The plant-insect interaction is a good model for this purpose with the availability of an enormous amount of information at the molecular and the ecosystem levels. Changing global climatic conditions are abruptly resetting plant-insect interactions. Integration of discretely located heterogeneous information from the ecosystem to genes and pathways will be an advantage to understand the complexity of plant-insect interactions. This review will present the recent developments in omics-based high-throughput experimental approaches, with particular emphasis on studying plant defence responses against insect attack. The review highlights the importance of using integrative systems approaches to study plant-insect interactions from the macroscopic to the microscopic level. We analyse the current efforts in generating, integrating and modelling multiomics data to understand plant-insect interaction at a systems level. As a future prospect, we highlight the growing interest in utilizing the synthetic biology platform for engineering insect-resistant plants.

Apoplastic Nucleoside Accumulation in Arabidopsis Leads to Reduced Photosynthetic Performance and Increased Susceptibility Against Botrytis cinerea

Interactions between plant and pathogen often occur in the extracellular space and especially nucleotides like ATP and NAD have been identified as key players in this scenario. Arabidopsis mutants accumulating nucleosides in the extracellular space were generated and studied with respect to susceptibility against Botrytis cinerea infection and general plant fitness determined as photosynthetic performance. The mutants used are deficient in the main nucleoside uptake system ENT3 and the extracellular nucleoside hydrolase NSH3. When grown on soil but not in hydroponic culture, these plants markedly accumulate adenosine and uridine in leaves. This nucleoside accumulation was accompanied by reduced photosystem II efficiency and altered expression of photosynthesis related genes. Moreover, a higher susceptibility toward Botrytis cinerea infection and a reduced induction of pathogen related genes PR1 and WRKY33 was observed. All these effects did not occur in hydroponically grown plants substantiating a contribution of extracellular nucleosides to these effects. Whether reduced general plant fitness, altered pathogen response capability or more direct interactions with the pathogen are responsible for these observations is discussed.

Improving crop disease resistance: lessons from research on Arabidopsis and tomato

One of the great challenges for food security in the 21st century is to improve yield stability through the development of disease-resistant crops. Crop research is often hindered by the lack of molecular tools, growth logistics, generation time and detailed genetic annotations, hence the power of model plant species. Our knowledge of plant immunity today has been largely shaped by the use of models, specifically through the use of mutants. We examine the importance of Arabidopsis and tomato as models in the study of plant immunity and how they help us in revealing a detailed and deep understanding of the various layers contributing to the immune system. Here we describe examples of how knowledge from models can be transferred to economically important crops resulting in new tools to enable and accelerate classical plant breeding. We will also discuss how models, and specifically transcriptomics and effectoromics approaches, have contributed to the identification of core components of the defense response which will be key to future engineering of durable and sustainable disease resistance in plants.

Deciphering the molecular bases for drought tolerance in Arabidopsis autotetraploids

Whole genome duplication (autopolyploidy) is common in many plant species and often leads to better adaptation to adverse environmental conditions. However, little is known about the physiological and molecular mechanisms underlying these adaptations. Drought is one of the major environmental conditions limiting plant growth and development. Here, we report that, in Arabidopsis thaliana, tetraploidy promotes alterations in cell proliferation and organ size in a tissue-dependent manner. Furthermore, it potentiates plant tolerance to salt and drought stresses and decreases transpiration rate, likely through controlling stomata density and closure, abscisic acid (ABA) signalling and reactive oxygen species (ROS) homeostasis. Our transcriptomic analyses revealed that tetraploidy mainly regulates the expression of genes involved in redox homeostasis and ABA and stress response. Taken together, our data have shed light on the molecular basis associated with stress tolerance in autopolyploid plants.

Transcriptome analysis reveals genes commonly induced by Botrytis cinerea infection, cold, drought and oxidative stresses in Arabidopsis

Signaling pathways controlling biotic and abiotic stress responses may interact synergistically or antagonistically. To identify the similarities and differences among responses to diverse stresses, we analyzed previously published microarray data on the transcriptomic responses of Arabidopsis to infection with Botrytis cinerea (a biotic stress), and to cold, drought, and oxidative stresses (abiotic stresses). Our analyses showed that at early stages after B. cinerea inoculation, 1498 genes were up-regulated (B. cinerea up-regulated genes; BUGs) and 1138 genes were down-regulated (B. cinerea down-regulated genes; BDGs). We showed a unique program of gene expression was activated in response each biotic and abiotic stress, but that some genes were similarly induced or repressed by all of the tested stresses. Of the identified BUGs, 25%, 6% and 12% were also induced by cold, drought and oxidative stress, respectively; whereas 33%, 7% and 5.5% of the BDGs were also down-regulated by the same abiotic stresses. Coexpression and protein-protein interaction network analyses revealed a dynamic range in the expression levels of genes encoding regulatory proteins. Analysis of gene expression in response to electrophilic oxylipins suggested that these compounds are involved in mediating responses to B. cinerea infection and abiotic stress through TGA transcription factors. Our results suggest an overlap among genes involved in the responses to biotic and abiotic stresses in Arabidopsis. Changes in the transcript levels of genes encoding components of the cyclopentenone signaling pathway in response to biotic and abiotic stresses suggest that the oxylipin signal transduction pathway plays a role in plant defense. Identifying genes that are commonly expressed in response to environmental stresses, and further analyzing the functions of their encoded products, will increase our understanding of the plant stress response. This information could identify targets for genetic modification to improve plant resistance to multiple stresses.

Increased resistance to biotrophic pathogens in the Arabidopsis constitutive induced resistance 1 mutant is EDS1 and PAD4-dependent and modulated by environmental temperature

The Arabidopsis constitutive induced resistance 1 (cir1) mutant displays salicylic acid (SA)-dependent constitutive expression of defence genes and enhanced resistance to biotrophic pathogens. To further characterise the role of CIR1 in plant immunity we conducted epistasis analyses with two key components of the SA-signalling branch of the defence network, ENHANCED DISEASE SUSCEPTIBILITY1 (EDS1) and PHYTOALEXIN DEFICIENT4 (PAD4). We demonstrate that the constitutive defence phenotypes of cir1 require both EDS1 and PAD4, indicating that CIR1 lies upstream of the EDS1-PAD4 regulatory node in the immune signalling network. In light of this finding we examined EDS1 expression in cir1 and observed increased protein, but not mRNA levels in this mutant, suggesting that CIR1 might act as a negative regulator of EDS1 via a post-transcriptional mechanism. Finally, as environmental temperature is known to influence the outcome of plant-pathogen interactions, we analysed cir1 plants grown at 18, 22 or 25°C. We found that susceptibility to Pseudomonas syringae pv. tomato (Pst) DC3000 is modulated by temperature in cir1. Greatest resistance to this pathogen (relative to PR-1:LUC control plants) was observed at 18°C, while at 25°C no difference in susceptibility between cir1 and control plants was apparent. The increase in resistance to Pst DC3000 at 18°C correlated with a stunted growth phenotype, suggesting that activation of defence responses may be enhanced at lower temperatures in the cir1 mutant.

Modelling transcriptional networks in leaf senescence

The process of leaf senescence is induced by an extensive range of developmental and environmental signals and controlled by multiple, cross-linking pathways, many of which overlap with plant stress-response signals. Elucidation of this complex regulation requires a step beyond a traditional one-gene-at-a-time analysis. Application of a more global analysis using statistical and mathematical tools of systems biology is an approach that is being applied to address this problem. A variety of modelling methods applicable to the analysis of current and future senescence data are reviewed and discussed using some senescence-specific examples. Network modelling with a senescence transcriptome time course followed by testing predictions with gene-expression data illustrates the application of systems biology tools.

Reversal of senescence by N resupply to N-starved Arabidopsis thaliana: transcriptomic and metabolomic consequences

Leaf senescence is a developmentally controlled process, which is additionally modulated by a number of adverse environmental conditions. Nitrogen shortage is a well-known trigger of precocious senescence in many plant species including crops, generally limiting biomass and seed yield. However, leaf senescence induced by nitrogen starvation may be reversed when nitrogen is resupplied at the onset of senescence. Here, the transcriptomic, hormonal, and global metabolic rearrangements occurring during nitrogen resupply-induced reversal of senescence in Arabidopsis thaliana were analysed. The changes induced by senescence were essentially in keeping with those previously described; however, these could, by and large, be reversed. The data thus indicate that plants undergoing senescence retain the capacity to sense and respond to the availability of nitrogen nutrition. The combined data are discussed in the context of the reversibility of the senescence programme and the evolutionary benefit afforded thereby. Future prospects for understanding and manipulating this process in both Arabidopsis and crop plants are postulated.

Correcting for link loss in causal network inference caused by regulator interference

MOTIVATION

There are a number of algorithms to infer causal regulatory networks from time series (gene expression) data. Here we analyse the phenomena of regulator interference, where regulators with similar dynamics mutually suppress both the probability of regulating a target and the associated link strength; for instance, interference between two identical strong regulators reduces link probabilities by ∼50%.

RESULTS

We construct a robust method to define an interference-corrected causal network based on an analysis of the conditional link probabilities that recovers links lost through interference. On a large real network (Streptomyces coelicolor, phosphate depletion), we demonstrate that significant interference can occur between regulators with a correlation as low as 0.865, losing an estimated 34% of links by interference. However, levels of interference cannot be predicted from the correlation between regulators alone and are data specific. Validating against known networks, we show that high numbers of functional links are lost by regulator interference. Performance against other methods on DREAM4 data is excellent.

AVAILABILITY AND IMPLEMENTATION

The method is implemented in R and is publicly available as the NIACS package at http://www2.warwick.ac.uk/fac/sci/systemsbiology/research/software.

Environmental History Modulates Arabidopsis Pattern-Triggered Immunity in a HISTONE ACETYLTRANSFERASE1-Dependent Manner

In nature, plants are exposed to a fluctuating environment, and individuals exposed to contrasting environmental factors develop different environmental histories. Whether different environmental histories alter plant responses to a current stress remains elusive. Here, we show that environmental history modulates the plant response to microbial pathogens. Arabidopsis thaliana plants exposed to repetitive heat, cold, or salt stress were more resistant to virulent bacteria than Arabidopsis grown in a more stable environment. By contrast, long-term exposure to heat, cold, or exposure to high concentrations of NaCl did not provide enhanced protection against bacteria. Enhanced resistance occurred with priming of Arabidopsis pattern-triggered immunity (PTI)-responsive genes and the potentiation of PTI-mediated callose deposition. In repetitively stress-challenged Arabidopsis, PTI-responsive genes showed enrichment for epigenetic marks associated with transcriptional activation. Upon bacterial infection, enrichment of RNA polymerase II at primed PTI marker genes was observed in environmentally challenged Arabidopsis. Finally, repetitively stress-challenged histone acetyltransferase1-1 (hac1-1) mutants failed to demonstrate enhanced resistance to bacteria, priming of PTI, and increased open chromatin states. These findings reveal that environmental history shapes the plant response to bacteria through the development of a HAC1-dependent epigenetic mark characteristic of a primed PTI response, demonstrating a mechanistic link between the primed state in plants and epigenetics.

Resistance to Botrytis cinerea in Solanum lycopersicoides involves widespread transcriptional reprogramming

BACKGROUND

Tomato (Solanum lycopersicum), one of the world's most important vegetable crops, is highly susceptible to necrotrophic fungal pathogens such as Botrytis cinerea and Alternaria solani. Improving resistance through conventional breeding has been hampered by a shortage of resistant germplasm and difficulties in introgressing resistance into elite germplasm without linkage drag. The goal of this study was to explore natural variation among wild Solanum species to identify new sources of resistance to necrotrophic fungi and dissect mechanisms underlying resistance against B. cinerea.

RESULTS

Among eight wild species evaluated for resistance against B. cinerea and A. solani, S. lycopersicoides expressed the highest levels of resistance against both pathogens. Resistance against B. cinerea manifested as containment of pathogen growth. Through next-generation RNA sequencing and de novo assembly of the S. lycopersicoides transcriptome, changes in gene expression were analyzed during pathogen infection. In response to B. cinerea, differentially expressed transcripts grouped into four categories: genes whose expression rapidly increased then rapidly decreased, genes whose expression rapidly increased and plateaued, genes whose expression continually increased, and genes with decreased expression. Homology-based searches also identified a limited number of highly expressed B. cinerea genes. Almost immediately after infection by B. cinerea, S. lycopersicoides suppressed photosynthesis and metabolic processes involved in growth, energy generation, and response to stimuli, and simultaneously induced various defense-related genes, including pathogenesis-related protein 1 (PR1), a beta-1,3-glucanase (glucanase), and a subtilisin-like protease, indicating a shift in priority towards defense. Moreover, cluster analysis revealed novel, uncharacterized genes that may play roles in defense against necrotrophic fungal pathogens in S. lycopersicoides. The expression of orthologous defense-related genes in S. lycopersicum after infection with B. cinerea revealed differences in the onset and intensity of induction, thus illuminating a potential mechanism explaining the increased susceptibility. Additionally, metabolic pathway analyses identified putative defense-related categories of secondary metabolites.

CONCLUSIONS

In sum, this study provided insight into resistance against necrotrophic fungal pathogens in the Solanaceae, as well as novel sequence resources for S. lycopersicoides.

Orchestration of plant defense systems: genes to populations

Research over the past decades has made immense progress in identifying some genes and mechanisms underlying plant defense against biotic organisms. The recent movement towards systems biology approaches has increased mechanistic knowledge, revealing a need for understanding how all the genes and mechanisms integrate to create a response to any given biotic interaction. This begins with evidence that diverse molecular patterns converge, suggesting that the plant perceives signals not the interacting species. These signals then coordinate across regulatory networks via molecular interactions and cause non-cell autonomous responses in neighboring and systemic cells. Finally, the identification of transporters is showing that plant defenses are harmonized across tissues and even show the potential for coordination across individuals within a population.

Machine learning-based differential network analysis: a study of stress-responsive transcriptomes in Arabidopsis

Machine learning (ML) is an intelligent data mining technique that builds a prediction model based on the learning of prior knowledge to recognize patterns in large-scale data sets. We present an ML-based methodology for transcriptome analysis via comparison of gene coexpression networks, implemented as an R package called machine learning-based differential network analysis (mlDNA) and apply this method to reanalyze a set of abiotic stress expression data in Arabidopsis thaliana. The mlDNA first used a ML-based filtering process to remove nonexpressed, constitutively expressed, or non-stress-responsive "noninformative" genes prior to network construction, through learning the patterns of 32 expression characteristics of known stress-related genes. The retained "informative" genes were subsequently analyzed by ML-based network comparison to predict candidate stress-related genes showing expression and network differences between control and stress networks, based on 33 network topological characteristics. Comparative evaluation of the network-centric and gene-centric analytic methods showed that mlDNA substantially outperformed traditional statistical testing-based differential expression analysis at identifying stress-related genes, with markedly improved prediction accuracy. To experimentally validate the mlDNA predictions, we selected 89 candidates out of the 1784 predicted salt stress-related genes with available SALK T-DNA mutagenesis lines for phenotypic screening and identified two previously unreported genes, mutants of which showed salt-sensitive phenotypes.

RNA-Seq derived identification of differential transcription in the chrysanthemum leaf following inoculation with Alternaria tenuissima

BACKGROUND

A major production constraint on the important ornamental species chrysanthemum is black spot which is caused by the necrotrophic fungus Alternaria tenuissima. The molecular basis of host resistance to A. tenuissima has not been studied as yet in any detail. Here, high throughput sequencing was taken to characterize the transcriptomic response of the chrysanthemum leaf to A. tenuissima inoculation.

RESULTS

The transcriptomic data was acquired using RNA-Seq technology, based on the Illumina HiSeq™ 2000 platform. Four different libraries derived from two sets of leaves harvested from either inoculated or mock-inoculated plants were characterized. Over seven million clean reads were generated from each library, each corresponding to a coverage of >350,000 nt. About 70% of the reads could be mapped to a set of chrysanthemum unigenes. Read frequency was used as a measure of transcript abundance and therefore as an identifier of differential transcription in the four libraries. The differentially transcribed genes identified were involved in photosynthesis, pathogen recognition, reactive oxygen species generation, cell wall modification and phytohormone signalling; in addition, a number of varied transcription factors were identified. A selection of 23 of the genes was transcription-profiled using quantitative RT-PCR to validate the RNA-Seq output.

CONCLUSIONS

A substantial body of chrysanthemum transcriptomic sequence was generated, which led to a number of insights into the molecular basis of the host response to A. tenuissima infection. Although most of the differentially transcribed genes were up-regulated by the presence of the pathogen, those involved in photosynthesis were down-regulated.

Network modeling to understand plant immunity

Deciphering the networks that underpin complex biological processes using experimental data remains a significant, but promising, challenge, a task made all the harder by the added complexity of host-pathogen interactions. The aim of this article is to review the progress in understanding plant immunity made so far by applying network modeling algorithms and to show how this computational/mathematical strategy is facilitating a systems view of plant defense. We review the different types of network modeling that have been used, the data required, and the type of insight that such modeling can provide. We discuss the current challenges in modeling the regulatory networks that underlie plant defense and the future developments that may help address these challenges.

Wigwams: identifying gene modules co-regulated across multiple biological conditions

Motivation: Identification of modules of co-regulated genes is a crucial first step towards dissecting the regulatory circuitry underlying biological processes. Co-regulated genes are likely to reveal themselves by showing tight co-expression, e.g. high correlation of expression profiles across multiple time series datasets. However, numbers of up- or downregulated genes are often large, making it difficult to discriminate between dependent co-expression resulting from co-regulation and independent co-expression. Furthermore, modules of co-regulated genes may only show tight co-expression across a subset of the time series, i.e. show condition-dependent regulation.

Results: Wigwams is a simple and efficient method to identify gene modules showing evidence for co-regulation in multiple time series of gene expression data. Wigwams analyzes similarities of gene expression patterns within each time series (condition) and directly tests the dependence or independence of these across different conditions. The expression pattern of each gene in each subset of conditions is tested statistically as a potential signature of a condition-dependent regulatory mechanism regulating multiple genes. Wigwams does not require particular time points and can process datasets that are on different time scales. Differential expression relative to control conditions can be taken into account. The output is succinct and non-redundant, enabling gene network reconstruction to be focused on those gene modules and combinations of conditions that show evidence for shared regulatory mechanisms. Wigwams was run using six Arabidopsis time series expression datasets, producing a set of biologically significant modules spanning different combinations of conditions.

Availability and implementation: A Matlab implementation of Wigwams, complete with graphical user interfaces and documentation, is available at: warwick.ac.uk/wigwams.

Contact:k.j.denby@warwick.ac.uk

Supplementary Data:Supplementary data are available at Bioinformatics online.

Arabidopsis phosphoglycerate dehydrogenase1 of the phosphoserine pathway is essential for development and required for ammonium assimilation and tryptophan biosynthesis

In plants, two independent serine biosynthetic pathways, the photorespiratory and glycolytic phosphoserine (PS) pathways, have been postulated. Although the photorespiratory pathway is well characterized, little information is available on the function of the PS pathway in plants. Here, we present a detailed characterization of phosphoglycerate dehydrogenases (PGDHs) as components of the PS pathway in Arabidopsis thaliana. All PGDHs localize to plastids and possess similar kinetic properties, but they differ with respect to their sensitivity to serine feedback inhibition. Furthermore, analysis of pgdh1 and phosphoserine phosphatase mutants revealed an embryo-lethal phenotype and PGDH1-silenced lines were inhibited in growth. Metabolic analyses of PGDH1-silenced lines grown under ambient and high CO2 conditions indicate a direct link between PS biosynthesis and ammonium assimilation. In addition, we obtained several lines of evidence for an interconnection between PS and tryptophan biosynthesis, because the expression of PGDH1 and phosphoserine aminotransferase1 is regulated by MYB51 and MYB34, two activators of tryptophan biosynthesis. Moreover, the concentration of tryptophan-derived glucosinolates and auxin were reduced in PGDH1-silenced plants. In essence, our results provide evidence for a vital function of PS biosynthesis for plant development and metabolism.

RNAseq-based transcriptome analysis of Lactuca sativa infected by the fungal necrotroph Botrytis cinerea

The fungal pathogen Botrytis cinerea establishes a necrotrophic interaction with its host plants, including lettuce (Lactuca sativa), causing it to wilt, collapse and eventually dry up and die, which results in serious economic losses. Global expression profiling using RNAseq and the newly sequenced lettuce genome identified a complex network of genes involved in the lettuce-B. cinerea interaction. The observed high number of differentially expressed genes allowed us to classify them according to the biological pathways in which they are implicated, generating a holistic picture. Most pronounced were the induction of the phenylpropanoid pathway and terpenoid biosynthesis, whereas photosynthesis was globally down-regulated at 48 h post-inoculation. Large-scale comparison with data available on the interaction of B. cinerea with the model plant Arabidopsis thaliana revealed both general and species-specific responses to infection with this pathogen. Surprisingly, expression analysis of selected genes could not detect significant systemic transcriptional alterations in lettuce leaves distant from the inoculation site. Additionally, we assessed the response of these lettuce genes to a biotrophic pathogen, Bremia lactucae, revealing that similar pathways are induced during compatible interactions of lettuce with necrotrophic and biotrophic pathogens.

Making open data work for plant scientists

Despite the clear demand for open data sharing, its implementation within plant science is still limited. This is, at least in part, because open data-sharing raises several unanswered questions and challenges to current research practices. In this commentary, some of the challenges encountered by plant researchers at the bench when generating, interpreting, and attempting to disseminate their data have been highlighted. The difficulties involved in sharing sequencing, transcriptomics, proteomics, and metabolomics data are reviewed. The benefits and drawbacks of three data-sharing venues currently available to plant scientists are identified and assessed: (i) journal publication; (ii) university repositories; and (iii) community and project-specific databases. It is concluded that community and project-specific databases are the most useful to researchers interested in effective data sharing, since these databases are explicitly created to meet the researchers' needs, support extensive curation, and embody a heightened awareness of what it takes to make data reuseable by others. Such bottom-up and community-driven approaches need to be valued by the research community, supported by publishers, and provided with long-term sustainable support by funding bodies and government. At the same time, these databases need to be linked to generic databases where possible, in order to be discoverable to the majority of researchers and thus promote effective and efficient data sharing. As we look forward to a future that embraces open access to data and publications, it is essential that data policies, data curation, data integration, data infrastructure, and data funding are linked together so as to foster data access and research productivity.

A mutation in the expansin-like A2 gene enhances resistance to necrotrophic fungi and hypersensitivity to abiotic stress in Arabidopsis thaliana

Expansins are cell wall loosening agents, known for their endogenous function in cell wall extensibility. The Arabidopsis expansin-like A2 (EXLA2) gene was identified by its down-regulation in response to infection by the necrotrophic pathogen Botrytis cinerea, and by the reduced susceptibility of an exla2 mutant to the same pathogen. The exla2 mutant was equally susceptible to Pseudomonas syringae pv. tomato, but was more resistant to the necrotrophic fungus Alternaria brassicicola, when compared with the wild-type or with transgenic, ectopic EXLA2-overexpressing lines. The exla2 mutants also enhanced tolerance to the phytoprostane-A1 . This suggests that the absence or down-regulation of EXLA2 leads to increased resistance to B. cinerea in a CORONATINE INSENSITIVE 1 (COI1)-dependent manner, and this down-regulation can be achieved by phytoprostane-A1 treatment. EXLA2 is induced significantly by salinity and cold, and by the exogenous application of abscisic acid. The exla2 mutant also showed hypersensitivity towards increased salt and cold, and this hypersensitivity required a functional abscisic acid pathway. The differential temporal expression of EXLA2 and the phenotypes in transgenic plants with altered expression of EXLA2 indicate that plant cell wall structure is an important player during Arabidopsis developmental stages. Our results indicate that EXLA2 appears to be important in response to various biotic and abiotic stresses, particularly in the pathogenesis of necrotrophic pathogens and in the tolerance to abiotic stress.

Alternative splicing at the intersection of biological timing, development, and stress responses

High-throughput sequencing for transcript profiling in plants has revealed that alternative splicing (AS) affects a much higher proportion of the transcriptome than was previously assumed. AS is involved in most plant processes and is particularly prevalent in plants exposed to environmental stress. The identification of mutations in predicted splicing factors and spliceosomal proteins that affect cell fate, the circadian clock, plant defense, and tolerance/sensitivity to abiotic stress all point to a fundamental role of splicing/AS in plant growth, development, and responses to external cues. Splicing factors affect the AS of multiple downstream target genes, thereby transferring signals to alter gene expression via splicing factor/AS networks. The last two to three years have seen an ever-increasing number of examples of functional AS. At a time when the identification of AS in individual genes and at a global level is exploding, this review aims to bring together such examples to illustrate the extent and importance of AS, which are not always obvious from individual publications. It also aims to ensure that plant scientists are aware that AS is likely to occur in the genes that they study and that dynamic changes in AS and its consequences need to be considered routinely.

Alternative splicing at the intersection of biological timing, development, and stress responses

High-throughput sequencing for transcript profiling in plants has revealed that alternative splicing (AS) affects a much higher proportion of the transcriptome than was previously assumed. AS is involved in most plant processes and is particularly prevalent in plants exposed to environmental stress. The identification of mutations in predicted splicing factors and spliceosomal proteins that affect cell fate, the circadian clock, plant defense, and tolerance/sensitivity to abiotic stress all point to a fundamental role of splicing/AS in plant growth, development, and responses to external cues. Splicing factors affect the AS of multiple downstream target genes, thereby transferring signals to alter gene expression via splicing factor/AS networks. The last two to three years have seen an ever-increasing number of examples of functional AS. At a time when the identification of AS in individual genes and at a global level is exploding, this review aims to bring together such examples to illustrate the extent and importance of AS, which are not always obvious from individual publications. It also aims to ensure that plant scientists are aware that AS is likely to occur in the genes that they study and that dynamic changes in AS and its consequences need to be considered routinely.

A local regulatory network around three NAC transcription factors in stress responses and senescence in Arabidopsis leaves

A model is presented describing the gene regulatory network surrounding three similar NAC transcription factors that have roles in Arabidopsis leaf senescence and stress responses. ANAC019, ANAC055 and ANAC072 belong to the same clade of NAC domain genes and have overlapping expression patterns. A combination of promoter DNA/protein interactions identified using yeast 1-hybrid analysis and modelling using gene expression time course data has been applied to predict the regulatory network upstream of these genes. Similarities and divergence in regulation during a variety of stress responses are predicted by different combinations of upstream transcription factors binding and also by the modelling. Mutant analysis with potential upstream genes was used to test and confirm some of the predicted interactions. Gene expression analysis in mutants of ANAC019 and ANAC055 at different times during leaf senescence has revealed a distinctly different role for each of these genes. Yeast 1-hybrid analysis is shown to be a valuable tool that can distinguish clades of binding proteins and be used to test and quantify protein binding to predicted promoter motifs.

High resolution imaging of temporal and spatial changes of subcellular ascorbate, glutathione and H₂O₂ distribution during Botrytis cinerea infection in Arabidopsis

In order to study the mechanisms behind the infection process of the necrotrophic fungus Botrytis cinerea, the subcellular distribution of hydrogen peroxide (H₂O₂) was monitored over a time frame of 96 h post inoculation (hpi) in Arabidopsis thaliana Col-0 leaves at the inoculation site (IS) and the area around the IS which was defined as area adjacent to the inoculation site (AIS). H₂O₂ accumulation was correlated with changes in the compartment-specific distribution of ascorbate and glutathione and chloroplast fine structure. This study revealed that the severe breakdown of the antioxidative system, indicated by a drop in ascorbate and glutathione contents at the IS at later stages of infection correlated with an accumulation of H₂O₂ in chloroplasts, mitochondria, cell walls, nuclei and the cytosol which resulted in the development of chlorosis and cell death, eventually visible as tissue necrosis. A steady increase of glutathione contents in most cell compartments within infected tissues (up to 600% in chloroplasts at 96 hpi) correlated with an accumulation of H₂O₂ in chloroplasts, mitochondria and cell walls at the AIS indicating that high glutathione levels could not prevent the accumulation of reactive oxygen species (ROS) which resulted in chlorosis. Summing up, this study reveals the intracellular sequence of events during Botrytis cinerea infection and shows that the breakdown of the antioxidative system correlated with the accumulation of H₂O₂ in the host cells. This resulted in the degeneration of the leaf indicated by severe changes in the number and ultrastructure of chloroplasts (e.g. decrease of chloroplast number, decrease of starch and thylakoid contents, increase of plastoglobuli size), chlorosis and necrosis of the leaves.

MEME-LaB: motif analysis in clusters

Summary: Genome-wide expression analysis can result in large numbers of clusters of co-expressed genes. Although there are tools for ab initio discovery of transcription factor-binding sites, most do not provide a quick and easy way to study large numbers of clusters. To address this, we introduce a web tool called MEME-LaB. The tool wraps MEME (an ab initio motif finder), providing an interface for users to input multiple gene clusters, retrieve promoter sequences, run motif finding and then easily browse and condense the results, facilitating better interpretation of the results from large-scale datasets.

Availability: MEME-LaB is freely accessible at: http://wsbc.warwick.ac.uk/wsbcToolsWebpage/.

Contact:p.e.brown@warwick.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.

Integration of genome-scale modeling and transcript profiling reveals metabolic pathways underlying light and temperature acclimation in Arabidopsis

Understanding metabolic acclimation of plants to challenging environmental conditions is essential for dissecting the role of metabolic pathways in growth and survival. As stresses involve simultaneous physiological alterations across all levels of cellular organization, a comprehensive characterization of the role of metabolic pathways in acclimation necessitates integration of genome-scale models with high-throughput data. Here, we present an integrative optimization-based approach, which, by coupling a plant metabolic network model and transcriptomics data, can predict the metabolic pathways affected in a single, carefully controlled experiment. Moreover, we propose three optimization-based indices that characterize different aspects of metabolic pathway behavior in the context of the entire metabolic network. We demonstrate that the proposed approach and indices facilitate quantitative comparisons and characterization of the plant metabolic response under eight different light and/or temperature conditions. The predictions of the metabolic functions involved in metabolic acclimation of Arabidopsis thaliana to the changing conditions are in line with experimental evidence and result in a hypothesis about the role of homocysteine-to-Cys interconversion and Asn biosynthesis. The approach can also be used to reveal the role of particular metabolic pathways in other scenarios, while taking into consideration the entirety of characterized plant metabolism.

A temporal switch model for estimating transcriptional activity in gene expression

Motivation: The analysis and mechanistic modelling of time series gene expression data provided by techniques such as microarrays, NanoString, reverse transcription–polymerase chain reaction and advanced sequencing are invaluable for developing an understanding of the variation in key biological processes. We address this by proposing the estimation of a flexible dynamic model, which decouples temporal synthesis and degradation of mRNA and, hence, allows for transcriptional activity to switch between different states.

Results: The model is flexible enough to capture a variety of observed transcriptional dynamics, including oscillatory behaviour, in a way that is compatible with the demands imposed by the quality, time-resolution and quantity of the data. We show that the timing and number of switch events in transcriptional activity can be estimated alongside individual gene mRNA stability with the help of a Bayesian reversible jump Markov chain Monte Carlo algorithm. To demonstrate the methodology, we focus on modelling the wild-type behaviour of a selection of 200 circadian genes of the model plant Arabidopsis thaliana. The results support the idea that using a mechanistic model to identify transcriptional switch points is likely to strongly contribute to efforts in elucidating and understanding key biological processes, such as transcription and degradation.

Contact:B.F.Finkenstadt@Warwick.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.

Activation of Defense Mechanisms against Pathogens in Mosses and Flowering Plants

During evolution, plants have developed mechanisms to cope with and adapt to different types of stress, including microbial infection. Once the stress is sensed, signaling pathways are activated, leading to the induced expression of genes with different roles in defense. Mosses (Bryophytes) are non-vascular plants that diverged from flowering plants more than 450 million years ago, allowing comparative studies of the evolution of defense-related genes and defensive metabolites produced after microbial infection. The ancestral position among land plants, the sequenced genome and the feasibility of generating targeted knock-out mutants by homologous recombination has made the moss Physcomitrella patens an attractive model to perform functional studies of plant genes involved in stress responses. This paper reviews the current knowledge of inducible defense mechanisms in P. patens and compares them to those activated in flowering plants after pathogen assault, including the reinforcement of the cell wall, ROS production, programmed cell death, activation of defense genes and synthesis of secondary metabolites and defense hormones. The knowledge generated in P. patens together with comparative studies in flowering plants will help to identify key components in plant defense responses and to design novel strategies to enhance resistance to biotic stress.

Quiescent and necrotrophic lifestyle choice during postharvest disease development

Insidious fungal infections by postharvest pathogens remain quiescent during fruit growth until, at a particular phase during fruit ripening and senescence, the pathogens switch to the necrotrophic lifestyle and cause decay. During ripening, fruits undergo physiological processes, such as activation of ethylene biosynthesis, cuticular changes, and cell-wall loosening-changes that are accompanied by a decline of antifungal compounds, both those that are preformed and those that are inducible secondary metabolites. Pathogen infection of the unripe host fruit initiates defensive signal-transduction cascades, culminating in accumulation of antifungal proteins that limit fungal growth and development. In contrast, development of the same pathogens during fruit ripening and storage activates a substantially different signaling network, one that facilitates aggressive fungal colonization. This review focuses on responses induced by the quiescent pathogens of postharvest diseases in unripe host fruits. New genome-scale experimental approaches have begun to delineate the complex and multiple networks of host and pathogen responses activated to maintain or to facilitate the transition from the quiescent to the necrotrophic lifestyle.

Integrating nitric oxide into salicylic acid and jasmonic acid/ ethylene plant defense pathways

Plant defense against pests and pathogens is known to be conferred by either salicylic acid (SA) or jasmonic acid (JA)/ethylene (ET) pathways, depending on infection or herbivore-grazing strategy. It is well attested that SA and JA/ET pathways are mutually antagonistic allowing defense responses to be tailored to particular biotic stresses. Nitric oxide (NO) has emerged as a major signal influencing resistance mediated by both signaling pathways but no attempt has been made to integrate NO into established SA/JA/ET interactions. NO has been shown to act as an inducer or suppressor of signaling along each pathway. NO will initiate SA biosynthesis and nitrosylate key cysteines on TGA-class transcription factors to aid in the initiation of SA-dependent gene expression. Against this, S-nitrosylation of NONEXPRESSOR OF PATHOGENESIS-RELATED PROTEINS1 (NPR1) will promote the NPR1 oligomerization within the cytoplasm to reduce TGA activation. In JA biosynthesis, NO will initiate the expression of JA biosynthetic enzymes, presumably to over-come any antagonistic effects of SA on JA-mediated transcription. NO will also initiate the expression of ET biosynthetic genes but a suppressive role is also observed in the S-nitrosylation and inhibition of S-adenosylmethionine transferases which provides methyl groups for ET production. Based on these data a model for NO action is proposed but we have also highlighted the need to understand when and how inductive and suppressive steps are used.

Impacts of resistance gene genetics, function, and evolution on a durable future

Studies on resistance gene function and evolution lie at the confluence of structural and molecular biology, genetics, and plant breeding. However, knowledge from these disparate fields has yet to be extensively integrated. This review draws on ideas and information from these different fields to elucidate the influences driving the evolution of different types of resistance genes in plants and the concurrent evolution of virulence in pathogens. It provides an overview of the factors shaping the evolution of recognition, signaling, and response genes in the context of emerging functional information along with a consideration of the new opportunities for durable resistance enabled by high-throughput DNA sequencing technologies.