Transcription dynamics in plant immunity

GmBTB/POZ promotes the ubiquitination and degradation of LHP1 to regulate the response of soybean to Phytophthora sojae

Phytophthora sojae is a pathogen that causes stem and root rot in soybean (Glycine max [L.] Merr.). We previously demonstrated that GmBTB/POZ, a BTB/POZ domain-containing nuclear protein, enhances resistance to P. sojae in soybean, via a process that depends on salicylic acid (SA). Here, we demonstrate that GmBTB/POZ associates directly with soybean LIKE HETEROCHROMATIN PROTEIN1 (GmLHP1) in vitro and in vivo and promotes its ubiquitination and degradation. Both overexpression and RNA interference analysis of transgenic lines demonstrate that GmLHP1 negatively regulates the response of soybean to P. sojae by reducing SA levels and repressing GmPR1 expression. The WRKY transcription factor gene, GmWRKY40, a SA-induced gene in the SA signaling pathway, is targeted by GmLHP1, which represses its expression via at least two mechanisms (directly binding to its promoter and impairing SA accumulation). Furthermore, the nuclear localization of GmLHP1 is required for the GmLHP1-mediated negative regulation of immunity, SA levels and the suppression of GmWRKY40 expression. Finally, GmBTB/POZ releases GmLHP1-regulated GmWRKY40 suppression and increases resistance to P. sojae in GmLHP1-OE hairy roots. These findings uncover a regulatory mechanism by which GmBTB/POZ-GmLHP1 modulates resistance to P. sojae in soybean, likely by regulating the expression of downstream target gene GmWRKY40.

Transcriptomic profiling of susceptible and resistant flax seedlings after Fusarium oxysporum lini infection

In this study transcriptome was analyzed on two fibrous varieties of flax: the susceptible Regina and the resistant Nike. The experiment was carried out on 2-week-old seedlings, because in this phase of development flax is the most susceptible to infection. We analyzed the whole seedlings, which allowed us to recognize the systemic response of the plants to the infection. We decided to analyze two time points: 24h and 48h, because our goal was to learn the mechanisms activated in the initial stages of infection, these points were selected based on the previous analysis of chitinase gene expression, whose increase in time of Fusarium oxysporum lini infection has been repeatedly confirmed both in the case of flax and other plant species. The results show that although qualitatively the responses of the two varieties are similar, it is the degree of the response that plays the role in the differences of their resistance to F. oxysporum.

Ca14-3-3 Interacts With CaWRKY58 to Positively Modulate Pepper Response to Low-Phosphorus Starvation

Low-phosphorus stress (LPS) and pathogen attack are two important stresses frequently experienced by plants in their natural habitats, but how plant respond to them coordinately remains under-investigated. Here, we demonstrate that CaWRKY58, a known negative regulator of the pepper (Capsicum annuum) response to attack by Ralstonia solanacearum, is upregulated by LPS. Virus-induced gene silencing (VIGS) and overexpression of CaWRKY58 in Nicotiana benthamiana plants in combination with chromatin immunoprecipitation (ChIP) and electrophoretic mobility shift assays (EMSA) demonstrated that CaWRKY58 positively regulates the response of pepper to LPS by directly targeting and regulating genes related to phosphorus-deficiency tolerance, including PHOSPHATE STARVATION RESPONSE1 (PHR1). Yeast two-hybrid assays revealed that CaWRKY58 interacts with a 14-3-3 protein (Ca14-3-3); this interaction was confirmed by pull-down, bimolecular fluorescence complementation (BiFC), and microscale thermophoresis (MST) assays. The interaction between Ca14-3-3 and CaWRKY58 enhanced the activation of PHR1 expression by CaWRKY58, but did not affect the expression of the immunity-related genes CaNPR1 and CaDEF1, which are negatively regulated by CaWRKY58 in pepper upon Ralstonia solanacearum inoculation. Collectively, our data indicate that CaWRKY58 negatively regulates immunity against Ralstonia solanacearum, but positively regulates tolerance to LPS and that Ca14-3-3 transcriptionally activates CaWRKY58 in response to LPS.

R2R3-MYB Transcription Factor MdMYB73 Confers Increased Resistance to the Fungal Pathogen Botryosphaeria dothidea in Apples via the Salicylic Acid Pathway

MYB transcription factors (TFs) participate in many biological processes. However, the molecular mechanisms by which MYB TFs affect plant resistance to apple ring rot remain poorly understood. Here, the R2R3-MYB gene MdMYB73 was cloned from "Royal Gala" apples and functionally characterized as a positive regulator of the defense response to Botryosphaeria dothidea. qRT-PCR and GUS staining demonstrated that MdMYB73 was strongly induced in apple fruits and transgenic calli after inoculation with B. dothidea. MdMYB73 overexpression improved resistance to B. dothidea in apple calli and fruits, while MdMYB73 suppression weakened. Increased resistance to B. dothidea was also observed in MdMYB73-expressing Arabidopsis thaliana. Interestingly, salicylic acid (SA) contents and the expression levels of genes related with SA synthesis and signaling were greater in MdMYB73-overexpressing plant materials compared to wild-type controls after inoculation, suggesting that MdMYB73 might enhance resistance to B. dothidea via the SA pathway. Finally, we discovered that MdMYB73 interacts with MdWRKY31, a positive regulator of B. dothidea. Together, MdWRKY31 and MdMYB73 enhanced B. dothidea resistance in apples. Our results clarify the mechanisms by which MdMYB73 improves resistance to B. dothidea and suggest that resistance may be affected by regulating the SA pathway.

A Bifunctional Synthetic Peptide With Antimicrobial and Plant Elicitation Properties That Protect Tomato Plants From Bacterial and Fungal Infections

The hybrid peptide BP178 (KKLFKKILKYLAGPAGIGKFLHSAKKDEL-OH), derived from BP100 (KKLFKKILKYL) and magainin (1-10), and engineered for plant expression, had a strong bactericidal activity but not fungicidal. Moreover, the preventive spray of tomato plants with BP178 controlled infections by the plant pathogenic bacteria Pseudomonas syringae pv. tomato and Xanthomonas campestris pv. vesicatoria, as well as the fungus Botrytis cinerea. The treatment of tomato plants with BP178 induced the expression of several genes according to microarray and RT-qPCR analysis. Upregulated genes coded for several pathogenesis-related proteins, including PR1, PR2, PR3, PR4, PR5, PR6, PR7, PR9, PR10, and PR14, as well as transcription factors like ethylene transcription factors, WRKY, NAC and MYB, involved in the salicylic acid, jasmonic acid, and ethylene-signaling pathways. BP178 induced a similar gene expression pattern to flg15 according to RT-qPCR analysis, whereas the parent peptide BP100 did not trigger such as a strong plant defense response. It was concluded that BP178 was a bifunctional peptide protecting the plant against pathogen infection through a dual mechanism of action consisting of antimicrobial activity against bacterial pathogens and plant defense elicitation on plant host.

A Meta-Analysis of Quantitative Trait Loci Associated with Multiple Disease Resistance in Rice (Oryza sativa L.)

Rice blast, sheath blight and bacterial leaf blight are major rice diseases found worldwide. The development of resistant cultivars is generally perceived as the most effective way to combat these diseases. Plant disease resistance is a polygenic trait where a combinatorial effect of major and minor genes affects this trait. To locate the source of this trait, various quantitative trait loci (QTL) mapping studies have been performed in the past two decades. However, investigating the congruency between the reported QTL is a daunting task due to the heterogeneity amongst the QTLs studied. Hence, the aim of our study is to integrate the reported QTLs for resistance against rice blast, sheath blight and bacterial leaf blight and objectively analyze and consolidate the location of QTL clusters in the chromosomes, reducing the QTL intervals and thus identifying candidate genes within the selected meta-QTL. A total of twenty-seven studies for resistance QTLs to rice blast (8), sheath blight (15) and bacterial leaf blight (4) was compiled for QTL projection and analyses. Cumulatively, 333 QTLs associated with rice blast (114), sheath blight (151) and bacterial leaf blight (68) resistance were compiled, where 303 QTLs could be projected onto a consensus map saturated with 7633 loci. Meta-QTL analysis on 294 QTLs yielded 48 meta-QTLs, where QTLs with membership probability lower than 60% were excluded, reducing the number of QTLs within the meta-QTL to 274. Further, three meta-QTL regions (MQTL2.5, MQTL8.1 and MQTL9.1) were selected for functional analysis on the basis that MQTL2.5 harbors the highest number of QTLs; meanwhile, MQTL8.1 and MQTL9.1 have QTLs associated with all three diseases mentioned above. The functional analysis allows for determination of enriched gene ontology and resistance gene analogs (RGAs) and other defense-related genes. To summarize, MQTL2.5, MQTL8.1 and MQTL9.1 have a considerable number of R-genes that account for 10.21%, 4.08% and 6.42% of the total genes found in these meta-QTLs, respectively. Defense genes constitute around 3.70%, 8.16% and 6.42% of the total number of genes in MQTL2.5, MQTL8.1 and MQTL9.1, respectively. This frequency is higher than the total frequency of defense genes in the rice genome, which is 0.0096% (167 defense genes/17,272 total genes). The integration of the QTLs facilitates the identification of QTL hotspots for rice blast, sheath blight and bacterial blight resistance with reduced intervals, which helps to reduce linkage drag in breeding. The candidate genes within the promising regions could be utilized for improvement through genetical engineering.

CaASR1 promotes salicylic acid- but represses jasmonic acid-dependent signaling to enhance the resistance of Capsicum annuum to bacterial wilt by modulating CabZIP63

CabZIP63 acts positively in the resistance of pepper (Capsicum annuum) to bacterial wilt caused by Ralstonia solanacearum or tolerance to high-temperature/high-humidity stress, but it is unclear how CabZIP63 achieves its functional specificity against R. solanacearum. Here, CaASR1, an abscisic acid-, stress-, and ripening-inducible protein of C. annuum, was functionally characterized in modulating the functional specificity of CabZIP63 during the defense response of pepper to R. solanacearum. In pepper plants inoculated with R. solanacearum, CaASR1 was up-regulated before 24 h post-inoculation but down-regulated thereafter, and was down-regulated by high-temperature/high-humidity stress. Data from gene silencing and transient overexpression experiments indicated that CaASR1 acts as a positive regulator in the immunity of pepper against R. solanacearum and a negative regulator of thermotolerance. Pull-down combined with mass spectrometry revealed that CaASR1 interacted with CabZIP63 upon R. solanacearum infection; the interaction was confirmed by microscale thermophoresis and bimolecular fluorescence complementation assays.CaASR1 silencing upon R. solanacearum inoculation repressed CabZIP63-mediated transcription from the promoters of the salicylic acid (SA)-dependent CaPR1 and CaNPR1, but derepressed transcription of CaHSP24 and the jasmonic acid (JA)-dependent CaDEF1. Our findings suggest that CaASR1 acts as a positive regulator of the defense response of pepper to R. solanacearum by interacting with CabZIP63, enabling it to promote SA-dependent but repress JA-dependent immunity and thermotolerance during the early stages of infection.

Transcriptomic analysis reveals insights into the response to Hop stunt viroid (HSVd) in sweet cherry (Prunus avium L.) fruits

Hop stunt viroid (HSVd) is a member of the genus Hostuviroid of the family Pospiviroidae and has been found in a wide range of herbaceous and woody hosts. It causes serious dapple fruit symptoms on infected sweet cherry, notably inducing cherry tree decay. In order to better understand the molecular mechanisms of HSVd infection in sweet cherry fruit, transcriptome analysis of HSVd-infected and healthy sweet cherry fruits was carried out. A total of 1,572 differentially expressed genes (DEGs) were identified, involving 961 upregulated DEGs and 611 downregulated DEGs. Functional analysis indicated that the DEGs were mainly involved in plant hormone signal transduction, plant-pathogen interactions, secondary metabolism, and the MAPK signaling pathway. In addition, C2H2 zinc finger, MYB, bHLH, AP2/ERF, C2C2-dof, NAC and WRKY transcription factors can respond to HSVd infection. In order to confirm the high-throughput sequencing results, 16 DEGs were verified by RT-qPCR analysis. The results provided insight into the pathways and genes of sweet cherry fruit in response to HSVd infection.

The overexpression of OsACBP5 protects transgenic rice against necrotrophic, hemibiotrophic and biotrophic pathogens

The most devastating diseases in rice (Oryza sativa) are sheath blight caused by the fungal necrotroph Rhizoctonia solani, rice blast by hemibiotrophic fungus Magnaporthe oryzae, and leaf blight by bacterial biotroph Xanthomonas oryzae (Xoo). It has been reported that the Class III acyl-CoA-binding proteins (ACBPs) such as those from dicots (Arabidopsis and grapevine) play a role in defence against biotrophic pathogens. Of the six Arabidopsis (Arabidopsis thaliana) ACBPs, AtACBP3 conferred protection in transgenic Arabidopsis against Pseudomonas syringae, but not the necrotrophic fungus, Botrytis cinerea. Similar to Arabidopsis, rice possesses six ACBPs, designated OsACBPs. The aims of this study were to test whether OsACBP5, the homologue of AtACBP3, can confer resistance against representative necrotrophic, hemibiotrophic and biotrophic phytopathogens and to understand the mechanisms in protection. Herein, when OsACBP5 was overexpressed in rice, the OsACBP5-overexpressing (OsACBP5-OE) lines exhibited enhanced disease resistance against representative necrotrophic (R. solani & Cercospora oryzae), hemibiotrophic (M. oryzae & Fusarium graminearum) and biotrophic (Xoo) phytopathogens. Progeny from a cross between OsACBP5-OE9 and the jasmonate (JA)-signalling deficient mutant were more susceptible than the wild type to infection by the necrotroph R. solani. In contrast, progeny from a cross between OsACBP5-OE9 and the salicylic acid (SA)-signalling deficient mutant was more susceptible to infection by the hemibiotroph M. oryzae and biotroph Xoo. Hence, enhanced resistance of OsACBP5-OEs against representative necrotrophs appears to be JA-dependent whilst that to (hemi)biotrophs is SA-mediated.

The fungal NADPH oxidase is an essential element for the molecular dialog between Trichoderma and Arabidopsis

Members of the fungal genus Trichoderma stimulate growth and reinforce plant immunity. Nevertheless, how fungal signaling elements mediate the establishment of a successful Trichoderma-plant interaction is largely unknown. In this work, we analyzed growth, root architecture and defense in an Arabidopsis-Trichoderma co-cultivation system, including the wild-type (WT) strain of the fungus and mutants affected in NADPH oxidase. Global gene expression profiles were assessed in both the plant and the fungus during the establishment of the interaction. Trichoderma atroviride WT improved root branching and growth of seedling as previously reported. This effect diminished in co-cultivation with the ∆nox1, ∆nox2 and ∆noxR null mutants. The data gathered of the Arabidopsis interaction with the ∆noxR strain showed that the seedlings had a heightened immune response linked to jasmonic acid in roots and shoots. In the fungus, we observed repression of genes involved in complex carbohydrate degradation in the presence of the plant before contact. However, in the absence of NoxR, such repression was lost, apparently due to a poor ability to adequately utilize simple carbon sources such as sucrose, a typical plant exudate. Our results unveiled the critical role played by the Trichoderma NoxR in the establishment of a fine-tuned communication between the plant and the fungus even before physical contact. In this dialog, the fungus appears to respond to the plant by adjusting its metabolism, while in the plant, fungal perception determines a delicate growth-defense balance.

Chromatin dynamics during interphase and cell division: similarities and differences between model and crop plants

Genetic information in the cell nucleus controls organismal development and responses to the environment, and finally ensures its own transmission to the next generations. To achieve so many different tasks, the genetic information is associated with structural and regulatory proteins, which orchestrate nuclear functions in time and space. Furthermore, plant life strategies require chromatin plasticity to allow a rapid adaptation to abiotic and biotic stresses. Here, we summarize current knowledge on the organization of plant chromatin and dynamics of chromosomes during interphase and mitotic and meiotic cell divisions for model and crop plants differing as to genome size, ploidy, and amount of genomic resources available. The existing data indicate that chromatin changes accompany most (if not all) cellular processes and that there are both shared and unique themes in the chromatin structure and global chromosome dynamics among species. Ongoing efforts to understand the molecular mechanisms involved in chromatin organization and remodeling have, together with the latest genome editing tools, potential to unlock crop genomes for innovative breeding strategies and improvements of various traits.

CaCBL1 Acts as a Positive Regulator in Pepper Response to Ralstonia solanacearum

Bacterial wilt caused by Ralstonia solanacearum is an important disease of pepper (Capsicum annuum), an economically important solanaceous vegetable worldwide, in particular, under high temperature (HT) conditions. However, the molecular mechanism underlying pepper immunity against bacterial wilt remains poorly understood. Herein, CaCBL1, a putative calcineurin B-like protein, was functionally characterized in the pepper response to R. solanacearum inoculation (RSI) under HT (RSI/HT). CaCBL1 was significantly upregulated by RSI at room temperature (RSI/RT), HT, or RSI/HT. CaCBL1-GFP fused protein targeted to whole epidermal cells of Nicotiana benthamiana when transiently overexpressed. CaCBL1 silencing by virus-induced gene silencing significantly enhanced pepper susceptibility to RSI under RT or HT, while its transient overexpression triggered hypersensitive response mimic cell death and upregulation of immunity-associated marker genes, including CabZIP63, CaWRKY40, and CaCDPK15, the positive regulators in the pepper response to RSI or HT found in our previous studies. In addition, by chromatin immunoprecipitation PCR and electrophoretic mobility shift assay, CaCBL1 was found to be directly targeted by CaWRKY40, although not by CaWRKY27 or CaWRKY58, via the W-box-2 within its promoter, and its transcription was found to be downregulated by silencing of CaWRKY40 while it was enhanced by its transient overexpression. These results suggest that CaCBL1 acts as a positive regulator in pepper immunity against R. solanacearum infection, constituting a positive feedback loop with CaWRKY40.

A systematic analysis of apple root resistance traits to Pythium ultimum infection and the underpinned molecular regulations of defense activation

Apple replant disease (ARD), caused by a pathogen complex, significantly impacts apple orchard establishment. The molecular regulation on ARD resistance has not been investigated until recently. A systematic phenotyping effort and a series of transcriptomic analyses were performed to uncover the underpinned molecular mechanism of apple root resistance to P. ultimum, a representative member in ARD pathogen complex. Genotype-specific plant survival rates and biomass reduction corresponded with microscopic features of necrosis progression patterns along the infected root. The presence of defined boundaries separating healthy and necrotic sections likely caused delayed necrosis expansion in roots of resistant genotypes compared with swift necrosis progression and profuse hyphae growth along infected roots of susceptible genotypes. Comprehensive datasets from a series of transcriptome analyses generated the first panoramic view of genome-wide transcriptional networks of defense activation between resistant and susceptible apple roots. Earlier and stronger molecular defense activation, such as pathogen perception and hormone signaling, may differentiate resistance from susceptibility in apple root. Delayed and interrupted activation of multiple defense pathways could have led to an inadequate resistance response. Using the panel of apple rootstock germplasm with defined resistant and susceptible phenotypes, selected candidate genes are being investigated by transgenic manipulation including CRISPR/Cas9 tools for their specific roles during apple root defense toward P. ultimum infection. Individual apple genes with validated functions regulating root resistance responses can be exploited for developing molecular tools for accurate and efficient incorporation of resistance traits into new apple rootstocks.

HOS15 and HDA9 negatively regulate immunity through histone deacetylation of intracellular immune receptor NLR genes in Arabidopsis

Plant immune responses need to be tightly controlled for growth-defense balance. The mechanism underlying this tight control is not fully understood. Here we identify epigenetic regulation of nucleotide-binding leucine rich repeat or Nod-Like Receptor (NLR) genes as an important mechanism for immune responses. Through a sensitized genetic screen and molecular studies, we identified and characterized HOS15 and its associated protein HDA9 as negative regulators of immunity and NLR gene expression. The loss-of-function of HOS15 or HDA9 confers enhanced resistance to pathogen infection accompanied with increased expression of one-third of the 207 NLR genes in Arabidopsis thaliana. HOS15 and HDA9 are physically associated with some of these NLR genes and repress their expression likely through reducing the acetylation of H3K9 at these loci. In addition, these NLR genes are repressed by HOS15 under both pathogenic and nonpathogenic conditions but by HDA9 only under infection condition. Together, this study uncovers a previously uncharacterized histone deacetylase complex in plant immunity and highlights the importance of epigenetic regulation of NLR genes in modulating growth-defense balance.

The Lifecycle of the Plant Immune System

Throughout their life span, plants confront an endless barrage of pathogens and pests. To successfully defend against biotic threats, plants have evolved a complex immune system responsible for surveillance, perception, and the activation of defense. Plant immunity requires multiple signaling processes, the outcome of which vary according to the lifestyle of the invading pathogen(s). In short, these processes require the activation of host perception, the regulation of numerous signaling cascades, and transcriptome reprograming, all of which are highly dynamic in terms of temporal and spatial scales. At the same time, the development of a single immune event is subjective to the development of plant immune system, which is co-regulated by numerous processes, including plant ontogenesis and the host microbiome. In total, insight into each of these processes provides a fuller understanding of the mechanisms that govern plant-pathogen interactions. In this review, we will discuss the "lifecycle" of plant immunity: the development of individual events of defense, including both local and distal processes, as well as the development and regulation of the overall immune system by ontogenesis regulatory genes and environmental microbiota. In total, we will integrate the output of recent discoveries and theories, together with several hypothetical models, to present a dynamic portrait of plant immunity.

Multiomics analysis of tolerant interaction of potato with potato virus Y

Potato virus Y (PVY) is the most economically important viral pathogen of potato worldwide. Different potato cultivars react to the pathogen differently, resulting in resistant, tolerant or disease outcome of the interaction. Here we focus on tolerant interaction between potato cv. Désirée and PVYNTN. To capture the response in its full complexity, we analyzed the dynamic changes on multiple molecular levels, including transcriptomics, sRNAomics, degradomics, proteomics and hormonomics. The analysis was complemented by the measurements of viral accumulation, photosynthetic activity and phenotypisation of the symptoms. Besides cv. Désirée we also studied its transgenic counterpart depleted for the accumulation of salicylic acid (NahG-Désirée). This multiomics analysis provides better insights into the mechanisms leading to tolerant response of potato to viral infection and can be used as a base in further studies of plant immunity regulation.

Effect of temperature on Pi54-mediated leaf blast resistance in rice

Assessment of temperature effect on plant resistance against diseases has become essential under climate change scenario as temperature rise is anticipated to modify host resistance. To determine temperature influence on resistance gene, a pair of near-isogenic rice lines differing for the Pi54 resistance gene was assessed against leaf blast. Blast resistance was determined as the extent of infection efficiency (IE) and sporulation (SP) at suboptimal (22 °C and 32 °C) and optimal temperature (27 °C) of pathogen aggressiveness. Relative resistance for IE and SP was higher at suboptimal temperature as compared to that of optimal temperature. Maximum level of resistance was at 22 °C where higher levels of expression of Pi54 and defence-regulatory transcription factor WRKY45 were also noted. At 32 °C, although some level of resistance noted, but level of Pi54 and WRKY45 expression was too low, suggesting that resistance recorded at higher temperature was due to reduced pathogen aggressiveness. At the optimal temperature for pathogen aggressiveness, comparatively lower levels of Pi54 and WRKY45 expression suggest possible temperature-induced interruption of the defence processes. The variation in resistance patterns modulated by temperature is appeared to be due to pathogen's sensitivity to temperature that leads to varying levels of Pi54 gene activation. Quick and violent activity of the pathogen at optimal temperature came into sight for the interruption of defence process activated by Pi54 gene. Evaluation of blast resistance genes under variable temperature conditions together with weather data could be applied in screening rice genotypes for selection of resistance having resilience to temperature rise.

RNA-Seq analysis of gene expression changes triggered by Xanthomonas oryzae pv. oryzae in a susceptible rice genotype

BACKGROUND

Xanthomonas oryzae pv. oryzae (Xoo) is a destructive disease in most of the rice growing regions worldwide. Xoo injects the transcriptional activator-like (TAL) effector protein into the host cell to induce the susceptibility (S) gene(s) for spreading the disease. In the current study, a susceptible rice genotype, JG30, was inoculated with wild Xoo strain PXO99^A and its mutant PH without any TAL effector, to retrieve the differentially expressed genes (DEGs) having a role in susceptibility.

RESULTS

RNA-Seq data analysis showed that 1143 genes were significantly differentially expressed (p-value ≤0.05) at 12, 24, 36 and 48 h post inoculation (hpi). Expression patterns, evaluated by quantitative real-time PCR (qRT-PCR), of randomly selected eight genes were similar to the RNA-Seq data. KEGG pathway classified the DEGs into photosynthesis and biosynthesis of phenylpropanoid pathway. Gene ontology (GO) analysis categorized the DEGs into the biological pathway, cellular component, and molecular function. We identified 43 differentially expressed transcription factors (TFs) belonging to different families. Also, clusters of the DEGs representing kinase and peroxidase responsive genes were retrieved. MapMan pathway analysis representing the expression pattern of genes expressed highly in biotic stress and metabolic pathways after PXO99^A infection relative to PH.

CONCLUSIONS

DEGs were identified in susceptible rice genotype inoculated with PXO99^A relative to mutant strain PH. The identified 1143 DEGs were predicted to be included in the different biological processes, signaling mechanism and metabolic pathways. The Jasmonic acid (JA) responsive genes were identified to be downregulated in PXO99^A infected leaves. This study would be useful for the researchers to reveal the potential functions of genes involved in the rice susceptibility to PXO99^A infection.

Soybean root transcriptome profiling reveals a nonhost resistant response during Heterodera glycines infection

Heterodera glycines (soybean cyst nematode, SCN) is one of the most devastating pathogens of soybean worldwide. The compatible and in compatible interactions between soybean and SCN have well documented. Nevertheless, the molecular mechanism of a nonhost resistant response in soybean against SCN infection remains obscure. Toward this end, a global transcriptional comparison was conducted between susceptible and resistant reactions of soybean roots infected by taking advantage of finding a new pathotype of SCN (SCNT). The soybean cultivar Lee, which exhibits resistant to SCNT and susceptible to HG 1.2.3.4.7 (SCNs) was utilized in the expriments. The results highlighted a nonhost resistant response of soybean. Transcriptome analysis indicated that the number of differentially expressed genes (DEGs) in the resistant interaction (3746) was much larger than that in the susceptible interaction (602). A great number of genes acting as intrinsic component of membrane, integral component of membrane, cell periphery and plasma membrance were remarkably enriched only in the resistant interaction, while the taurine and hypotaurine, phenylpropanoid pathway, plant-pathogen interaction and transcript factors were modulated in both interactions. This is the first study to examine genes expression patterns in a soybean genotype in response to invasion by a virulent and avirulent SCN population at the transcriptional level, which will provide insights into the complicate molecular mechanism of the nonhost resistant interaction.

Revisiting the Role of Transcription Factors in Coordinating the Defense Response Against Citrus Bark Cracking Viroid Infection in Commercial Hop (Humulus Lupulus L.)

Transcription factors (TFs) play a major role in controlling gene expression by intricately regulating diverse biological processes such as growth and development, the response to external stimuli and the activation of defense responses. The systematic identification and classification of TF genes are essential to gain insight into their evolutionary history, biological roles, and regulatory networks. In this study, we performed a global mining and characterization of hop TFs and their involvement in Citrus bark cracking viroid CBCVd infection by employing a digital gene expression analysis. Our systematic analysis resulted in the identification of a total of 3,818 putative hop TFs that were classified into 99 families based on their conserved domains. A phylogenetic analysis classified the hop TFs into several subgroups based on a phylogenetic comparison with reference TF proteins from Arabidopsis thaliana providing glimpses of their evolutionary history. Members of the same subfamily and subgroup shared conserved motif compositions. The putative functions of the CBCVd-responsive hop TFs were predicted using their orthologous counterparts in A. thaliana. The analysis of the expression profiling of the CBCVd-responsive hop TFs revealed a massive differential modulation, and the expression of the selected TFs was validated using qRT-PCR. Together, the comprehensive integrated analysis in this study provides better insights into the TF regulatory networks associated with CBCVd infections in the hop, and also offers candidate TF genes for improving the resistance in hop against viroids.

Genome-wide analysis of three histone marks and gene expression in Paulownia fortunei with phytoplasma infection

BACKGROUND

Paulownia withes'-broom (PaWB) disease caused by phytoplasma is a serious infectious disease for Paulownia. However, the underlying molecular pathogenesis is not fully understood. Recent studies have demonstrated that histone modifications could play a role in plant defense responses to pathogens. But there is still no available genome-wide histone modification data in non-model ligneous species infected with phytoplasma.

RESULTS

Here, we provided the first genome-wide profiles of three histone marks (H3K4me3, H3K36me3 and H3K9ac) in Paulownia fortunei under phytoplasma stress by using chromatin immunoprecipitation sequencing (ChIP-Seq). We found that H3K4me3, H3K36me3 and H3K9ac were mainly enriched in the genic regions in P. fortunei with (PFI) and without (PF) phytoplasma infection. ChIP-Seq analysis revealed 1738, 986, and 2577 genes were differentially modified by H3K4me3, H3K36me3 and H3K9ac marks in PFI under phytoplasma infection, respectively. The functional analysis of these genes suggested that most of them were mainly involved in metabolic pathways, biosynthesis of secondary metabolites, phenylpropanoid biosynthesis, plant-pathogen interaction and plant hormone signal transduction. In addition, the combinational analysis of ChIP-Seq and RNA-Seq showed that differential histone methylation and acetylation only affected a small subset of phytoplasma-responsive genes.

CONCLUSIONS

Taken together, this is the first report of integrated analysis of histone modifications and gene expression involved in Paulownia-phytoplasma interaction. Our results will provide the valuable resources for the mechanism studies of gene regulation in non-model plants upon pathogens attack.

CaLRR-RLK1, a novel RD receptor-like kinase from Capsicum annuum and transcriptionally activated by CaHDZ27, act as positive regulator in Ralstonia solanacearum resistance

BACKGROUND

Bacterial wilt caused by Ralstonia solanacearum is one of the most important diseases in pepper worldwide, however, the molecular mechanism underlying pepper resistance to bacterial wilt remains poorly understood.

RESULTS

Herein, a novel RD leucine-rich repeat receptor-like kinase, CaLRR-RLK1, was functionally characterized in immunity against R. solanacearum. CaLRR-RLK1 was targeted exclusively to plasma membrane and was up-regulated by R. solanacearum inoculation (RSI) as well as by the exogenous application of salicylic acid (SA), methyl jasmonate (MeJA) or ethephon (ETH). The silencing of CaLRR-RLK1 led to enhanced susceptibility of pepper plants to RSI, accompanied by down-regulation of immunity-related genes including CaACO1, CaHIR1, CaPR4 and CaPO2. In contrast, transient overexpression of CaLRR-RLK1 triggered hypersensitive response (HR)-like cell death and H2O2 accumulation in pepper leaves, manifested by darker trypan blue and DAB staining respectively. In addition, the ectopic overexpression of CaLRR-RLK1 in tobacco plants enhanced resistance R. solanacearum, accompanied with the immunity associated marker genes including NtPR2, NtPR2, NtHSR203 and NtHSR515. Furthermore, it was found that CaHDZ27, a positive regulator in pepper response to RSI in our previous study, transcriptionally activated CaLRR-RLK1 by direct targeting its promoter probably in a CAATTATTG dependent manner.

CONCLUSION

The study revealed that CaLRR-RLK1 confers pepper resistance to R. solanacearum as the direct targeting of CaHDZ27.

Genotype-specific suppression of multiple defense pathways in apple root during infection by Pythium ultimum

The genotype-specific defense activation in the roots of perennial tree crops to soilborne necrotrophic pathogens remains largely unknown. A recent phenotyping study indicated that the apple rootstock genotypes B.9 and G.935 have contrasting resistance responses to infection by Pythium ultimum. In the current study, a comparative transcriptome analysis by Illumina Solexa HiSeq 3000 platform was carried out to identify the global transcriptional regulation networks between the susceptible B.9 and the resistant G.935 to P. ultimum infection. Thirty-six libraries were sequenced to cover three timepoints after pathogen inoculation, with three biological replicates for each sample. The transcriptomes in the roots of the susceptible genotype B.9 were reflected by overrepresented differentially expressed genes (DEGs) with downregulated patterns and systematic suppression of cellular processes at 48 h post inoculation (hpi). In contrast, DEGs with annotated functions, such as kinase receptors, MAPK signaling, JA biosynthesis enzymes, transcription factors, and transporters, were readily induced at 24 hpi and continued up-regulation at 48 hpi in G.935 roots. The earlier and stronger defense activation is likely associated with an effective inhibition of necrosis progression in G.935 roots. Lack of effector-triggered immunity or existence of a susceptibility gene could contribute to the severely disturbed transcriptome and susceptibility in B.9 roots. The identified DEGs constitute a valuable resource for hypothesis-driven studies to elucidate the resistance/tolerance mechanisms in apple roots and validating their potential association with resistance traits.

Salicylic acid activates poplar defense against the biotrophic rust fungus Melampsora larici-populina via increased biosynthesis of catechin and proanthocyanidins

Poplar trees synthesize flavan-3-ols (catechin and proanthocyanidins) as a defense against foliar rust fungi, but the regulation of this defense response is poorly understood. Here, we investigated the role of hormones in regulating flavan-3-ol accumulation in poplar during rust infection. We profiled levels of defense hormones, signaling genes, and flavan-3-ol metabolites in black poplar leaves at different stages of rust infection. Hormone levels were manipulated by external sprays, genetic engineering, and drought to reveal their role in rust fungal defenses. Levels of salicylic acid (SA), jasmonic acid, and abscisic acid increased in rust-infected leaves and activated downstream signaling, with SA levels correlating closely with those of flavan-3-ols. Pretreatment with the SA analog benzothiadiazole increased flavan-3-ol accumulation by activating the MYB-bHLH-WD40 complex and reduced rust proliferation. Furthermore, transgenic poplar lines overproducing SA exhibited higher amounts of flavan-3-ols constitutively via the same transcriptional activation mechanism. These findings suggest a strong association among SA, flavan-3-ol biosynthesis, and rust resistance in poplars. Abscisic acid also promoted poplar defense against rust infection, but likely through stomatal immunity independent of flavan-3-ols. Jasmonic acid did not confer any apparent defense responses to the fungal pathogen. We conclude that SA activates flavan-3-ol biosynthesis in poplar against rust infection.

Differential and tissue-specific activation pattern of the AtPROPEP and AtPEPR genes in response to biotic and abiotic stress in Arabidopsis thaliana

In Arabidopsis thaliana AtPEPR1 and AtPEPR2 act as the receptors for the endogenous AtPROPEP-derived Pep peptides and subsequently initiate defense-signaling cascades. In the previous work,⁹ the expression pattern of the genes encoding the PEPR receptors and the AtPROPEP peptide precursor proteins was studied using promoter-GUS reporter constructs. Here, using the same constructs to study their expression pattern under biotic and abiotic stress, including AtPep1, flg22, methyl jasmonate (MeJA), and NaCl treatments, we observed that in response to AtPep1 and flg22, the activation of AtPEPR1 promoter was different from AtPEPR2. We also found that these promoters were differentially activated in response to NaCl. Remarkably, we showed that it is possible to classify the genes of the AtPROPEP family, based on the response of their promoters to the various stimuli employed: thus, we classify AtPROPEP1 in one group; AtPROPEP2 and AtPROPEP3 in a second group; AtPROPEP4, AtPROPEP7 and AtPROPEP8 in a third group and AtPROPEP5 in a fourth group. Our finding, confirm non-redundant roles among the members of the AtPROPEP family and their corresponding receptors.

Capsicum annuum HsfB2a Positively Regulates the Response to Ralstonia solanacearum Infection or High Temperature and High Humidity Forming Transcriptional Cascade with CaWRKY6 and CaWRKY40

The responses of pepper (Capsicum annuum) plants to inoculation with the pathogenic bacterium Ralstonia solanacearum and to high-temperature-high-humidity (HTHH) conditions were previously found to be coordinated by the transcription factors CaWRKY6 and CaWRKY40; however, the underlying molecular mechanism was unclear. Herein, we identified and functionally characterized CaHsfB2a, a nuclear-localized heat shock factor involved in pepper immunity to R. solanacearum inoculation (RSI) and tolerance to HTHH. CaHsfB2a is transcriptionally induced in pepper plants by RSI or HTHH and by exogenous application of salicylic acid (SA), methyl jasmonate (MeJA), ethylene (ETH), or abscisic acid (ABA). Virus-induced gene silencing (VIGS) of CaHsfB2a significantly impaired pepper immunity to RSI, hampered HTHH tolerance, and curtailed expression of immunity- and thermotolerance-associated marker genes such as CaHIR1, CaNPR1, CaABR1, and CaHSP24. Likewise, transient overexpression of CaHsfB2a in pepper leaves induced hypersensitive response (HR)-like cell death and H2O2 accumulation and upregulated the above-mentioned marker genes as well as CaWRKY6 and CaWRKY40. Chromatin immunoprecipitation (ChIP) and microscale thermophoresis (MST) analysis revealed that CaHsfB2a bound the promoters of both CaWRKY6 and CaWRKY40. In a parallel experiment, we determined by ChIP-PCR and MST that CaHsfB2a was regulated directly by CaWRKY40 but indirectly by CaWRKY6. Cumulatively, our results suggest that CaHsfB2a positively regulates plant immunity against RSI and tolerance to HTHH, via transcriptional cascades and positive feedback loops involving CaWRKY6 and CaWRKY40.

Transcriptomic basis for reinforcement of elm antiherbivore defence mediated by insect egg deposition

Plant responses to insect egg depositions are known to shape subsequent defensive responses to larvae hatching from the eggs. Elm (Ulmus minor) leaves, on which elm leaf beetles laid their eggs, mount a more efficient defence against larvae hatching from the eggs. However, the molecular mechanisms of this egg-mediated, improved defence are insufficiently understood and have so far only been studied in annual plants. We analysed the dynamics of transcriptomic changes in larval feeding-damaged elm leaves with and without prior egg deposition using de novo assembled RNA-seq data. Compared to egg-free leaves, egg deposition-treated leaves showed earlier and/or faster transcriptional regulations, as well as slightly enhanced differential transcriptional regulation after the onset of larval feeding. These early responding transcripts were overrepresented in gene ontology terms associated with post-translational protein modification, signalling and stress (defence) responses. We found evidence of transcriptional memory in initially egg deposition-induced transcripts whose differential expression was reset prior to larval hatching, but was more rapidly induced again by subsequent larval feeding. This potential memory effect of prior egg deposition, as well as the earlier/faster and enhanced feeding-induced differential regulation of transcripts in egg deposition-treated leaves, may contribute to the egg-mediated reinforcing effect on the elm's defence against larvae. Hence, our study shows that a plant's experience of a stress-indicating environmental cue (here: insect eggs) can push the dynamics of the plant's transcriptomic response to subsequent stress (here: larval feeding). Such experience-mediated acceleration of a stress-induced plant response may result in improved stress resistance.

Identification of a strawberry NPR-like gene involved in negative regulation of the salicylic acid-mediated defense pathway

Hormonal modulation plays a central role in triggering various resistant responses to biotic and abiotic stresses in plants. In cultivated strawberry (Fragaria x ananassa), the salicylic acid (SA)-dependent defense pathway has been associated with resistance to Colletotrichum spp. and the other pathogens. To better understand the SA-mediated defense mechanisms in strawberry, we analyzed two strawberry cultivars treated with SA for their resistance to anthracnose and gene expression profiles at 6, 12, 24, and 48 hr post-treatment. Strawberry genes related to SA biosynthesis, perception, and signaling were identified from SA-responsive transcriptomes of the two cultivars, and the induction of 17 candidate genes upon SA treatment was confirmed by qRT-PCR. Given the pivotal role of the non-expressor of pathogenesis-related (NPR) family in controlling the SA-mediated defense signaling pathway, we then analyzed NPR orthologous genes in strawberry. From the expression profile, FaNPRL-1 [ortholog of FvNPRL-1 (gene20070 in F. vesca)] was identified as an NPR-like gene significantly induced after SA treatment in both cultivars. With a conserved BTB/POZ domain, ankyrin repeat domain, and nuclear localization signal, FvNPRL-1 was found phylogenetically closer to NPR3/NPR4 than NPR1 in Arabidopsis. Ectopic expression of FvNPRL-1 in the Arabidopsis thaliana wild type suppressed the SA-mediated PR1 expression and the resistance to Pseudomonas syringae pv. tomato DC3000. Transient expression of FvNPRL-1 fused with green fluorescent protein in Arabidopsis protoplasts showed that SA affected nuclear translocation of FvNPRL-1. FvNPRL-1 likely functions similar to Arabidopsis NPR3/NPR4 as a negative regulator of the SA-mediated defense.

l-lysine metabolism to N-hydroxypipecolic acid: an integral immune-activating pathway in plants

l-lysine catabolic routes in plants include the saccharopine pathway to α-aminoadipate and decarboxylation of lysine to cadaverine. The current review will cover a third l-lysine metabolic pathway having a major role in plant systemic acquired resistance (SAR) to pathogen infection that was recently discovered in Arabidopsis thaliana. In this pathway, the aminotransferase AGD2-like defense response protein (ALD1) α-transaminates l-lysine and generates cyclic dehydropipecolic (DP) intermediates that are subsequently reduced to pipecolic acid (Pip) by the reductase SAR-deficient 4 (SARD4). l-pipecolic acid, which occurs ubiquitously in the plant kingdom, is further N-hydroxylated to the systemic acquired resistance (SAR)-activating metabolite N-hydroxypipecolic acid (NHP) by flavin-dependent monooxygenase1 (FMO1). N-hydroxypipecolic acid induces the expression of a set of major plant immune genes to enhance defense readiness, amplifies resistance responses, acts synergistically with the defense hormone salicylic acid, promotes the hypersensitive cell death response and primes plants for effective immune mobilization in cases of future pathogen challenge. This pathogen-inducible NHP biosynthetic pathway is activated at the transcriptional level and involves feedback amplification. Apart from FMO1, some cytochrome P450 monooxygenases involved in secondary metabolism catalyze N-hydroxylation reactions in plants. In specific taxa, pipecolic acid might also serve as a precursor in the biosynthesis of specialized natural products, leading to C-hydroxylated and otherwise modified piperidine derivatives, including indolizidine alkaloids. Finally, we show that NHP is glycosylated in Arabidopsis to form a hexose-conjugate, and then discuss open questions in Pip/NHP-related research.

CaC3H14 encoding a tandem CCCH zinc finger protein is directly targeted by CaWRKY40 and positively regulates the response of pepper to inoculation by Ralstonia solanacearum

Tandem CCCH zinc finger (TZnF) proteins have been implicated in plant defence, but their role in pepper (Capsicum annuum) is unclear. In the present study, the role of CaC3H14, a pepper TZnF protein, in the immune response of pepper plants to Ralstonia solanacearum infection was characterized. When fused to the green fluorescent protein, CaC3H14 was localized exclusively to the nuclei in leaf cells of Nicotiana benthamiana plants transiently overexpressing CaC3H14. Transcript abundance of CaC3H14 was up-regulated by inoculation with R. solanacearum. Virus-induced silencing of CaC3H14 increased the susceptibility of the plants to R. solanacearum and down-regulated the genes associated with the hypersensitive response (HR), specifically HIR1 and salicylic acid (SA)-dependent PR1a. By contrast, silencing resulted in the up-regulation of jasmonic acid (JA)-dependent DEF1 and ethylene (ET) biosynthesis-associated ACO1. Transient overexpression of CaC3H14 in pepper triggered an intensive HR, indicated by cell death and hydrogen peroxide (H₂ O₂ ) accumulation, up-regulated PR1a and down-regulated DEF1 and ACO1. Ectopic overexpression of CaC3H14 in tobacco plants significantly decreased the susceptibility of tobacco plants to R. solanacearum. It also up-regulated HR-associated HSR515, immunity-associated GST1 and the SA-dependent marker genes NPR1 and PR2, but down-regulated JA-dependent PR1b and ET-dependent EFE26. The CaC3H14 promoter and was bound and its transcription was up-regulated by CaWRKY40. Collectively, these results indicate that CaC3H14 is transcriptionally targeted by CaWRKY40, is a modulator of the antagonistic interaction between SA and JA/ET signalling, and enhances the defence response of pepper plants to infection by R. solanacearum.

Transcriptome dynamics associated with resistance and susceptibility against fusarium head blight in four wheat genotypes

Background

Fusarium head blight (FHB) of wheat in North America is caused mostly by the fungal pathogen Fusarium graminearum (Fg). Upon exposure to Fg, wheat initiates a series of cellular responses involving massive transcriptional reprogramming. In this study, we analyzed transcriptomics data of four wheat genotypes (Nyubai, Wuhan 1, HC374, and Shaw), at 2 and 4 days post inoculation (dpi) with Fg, using RNA-seq technology.

Results

A total of 37,772 differentially expressed genes (DEGs) were identified, 28,961 from wheat and 8811 from the pathogen. The susceptible genotype Shaw exhibited the highest number of host and pathogen DEGs, including 2270 DEGs associating with FHB susceptibility. Protein serine/threonine kinases and LRR-RK were associated with susceptibility at 2 dpi, while several ethylene-responsive, WRKY, Myb, bZIP and NAC-domain containing transcription factors were associated with susceptibility at 4 dpi. In the three resistant genotypes, 220 DEGs were associated with resistance. Glutathione S-transferase (GST), membrane proteins and distinct LRR-RKs were associated with FHB resistance across the three genotypes. Genes with unique, high up-regulation by Fg in Wuhan 1 were mostly transiently expressed at 2 dpi, while many defense-associated genes were up-regulated at both 2 and 4 dpi in Nyubai; the majority of unique genes up-regulated in HC374 were detected at 4 dpi only. In the pathogen, most genes showed increased expression between 2 and 4 dpi in all genotypes, with stronger levels in the susceptible host; however two pectate lyases and a hydrolase were expressed higher at 2 dpi, and acetyltransferase activity was highly enriched at 4 dpi.

Conclusions

There was an early up-regulation of LRR-RKs, different between susceptible and resistant genotypes; subsequently, distinct sets of genes associated with defense response were up-regulated. Differences in expression profiles among the resistant genotypes indicate genotype-specific defense mechanisms. This study also shows a greater resemblance in transcriptomics of HC374 to Nyubai, consistent with their sharing of two FHB resistance QTLs on 3BS and 5AS, compared to Wuhan 1 which carries one QTL on 2DL in common with HC374.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5012-3) contains supplementary material, which is available to authorized users.

Transcriptome dynamics associated with resistance and susceptibility against fusarium head blight in four wheat genotypes

BACKGROUND

Fusarium head blight (FHB) of wheat in North America is caused mostly by the fungal pathogen Fusarium graminearum (Fg). Upon exposure to Fg, wheat initiates a series of cellular responses involving massive transcriptional reprogramming. In this study, we analyzed transcriptomics data of four wheat genotypes (Nyubai, Wuhan 1, HC374, and Shaw), at 2 and 4 days post inoculation (dpi) with Fg, using RNA-seq technology.

RESULTS

A total of 37,772 differentially expressed genes (DEGs) were identified, 28,961 from wheat and 8811 from the pathogen. The susceptible genotype Shaw exhibited the highest number of host and pathogen DEGs, including 2270 DEGs associating with FHB susceptibility. Protein serine/threonine kinases and LRR-RK were associated with susceptibility at 2 dpi, while several ethylene-responsive, WRKY, Myb, bZIP and NAC-domain containing transcription factors were associated with susceptibility at 4 dpi. In the three resistant genotypes, 220 DEGs were associated with resistance. Glutathione S-transferase (GST), membrane proteins and distinct LRR-RKs were associated with FHB resistance across the three genotypes. Genes with unique, high up-regulation by Fg in Wuhan 1 were mostly transiently expressed at 2 dpi, while many defense-associated genes were up-regulated at both 2 and 4 dpi in Nyubai; the majority of unique genes up-regulated in HC374 were detected at 4 dpi only. In the pathogen, most genes showed increased expression between 2 and 4 dpi in all genotypes, with stronger levels in the susceptible host; however two pectate lyases and a hydrolase were expressed higher at 2 dpi, and acetyltransferase activity was highly enriched at 4 dpi.

CONCLUSIONS

There was an early up-regulation of LRR-RKs, different between susceptible and resistant genotypes; subsequently, distinct sets of genes associated with defense response were up-regulated. Differences in expression profiles among the resistant genotypes indicate genotype-specific defense mechanisms. This study also shows a greater resemblance in transcriptomics of HC374 to Nyubai, consistent with their sharing of two FHB resistance QTLs on 3BS and 5AS, compared to Wuhan 1 which carries one QTL on 2DL in common with HC374.

Comparative transcriptome analysis of the interaction between Actinidia chinensis var. chinensis and Pseudomonas syringae pv. actinidiae in absence and presence of acibenzolar-S-methyl

BACKGROUND

Since 2007, bacterial canker caused by Pseudomonas syringae pv. actinidiae (Psa) has become a pandemic disease leading to important economic losses in every country where kiwifruit is widely cultivated. Options for controlling this disease are very limited and rely primarily on the use of bactericidal compounds, such as copper, and resistance inducers. Among the latter, the most widely studied is acibenzolar-S-methyl. To elucidate the early molecular reaction of kiwifruit plants (Actinidia chinensis var. chinensis) to Psa infection and acibenzolar-S-methyl treatment, a RNA seq analysis was performed at different phases of the infection process, from the epiphytic phase to the endophytic invasion on acibenzolar-S-methyl treated and on non-treated plants. The infection process was monitored in vivo by confocal laser scanning microscopy.

RESULTS

De novo assembly of kiwifruit transcriptome revealed a total of 39,607 transcripts, of which 3360 were differentially expressed during the infection process, primarily 3 h post inoculation. The study revealed the coordinated changes of important gene functional categories such as signaling, hormonal balance and transcriptional regulation. Among the transcription factor families, AP2/ERF, MYB, Myc, bHLH, GATA, NAC, WRKY and GRAS were found differentially expressed in response to Psa infection and acibenzolar-S-methyl treatment. Finally, in plants treated with acibenzolar-S-methyl, a number of gene functions related to plant resistance, such as PR proteins, were modulated, suggesting the set-up of a more effective defense response against the pathogen. Weighted-gene coexpression network analysis confirmed these results.

CONCLUSIONS

Our work provides an in-depth description of the plant molecular reactions to Psa, it highlights the metabolic pathway related to acibenzolar-S-methyl-induced resistance and it contributes to the development of effective control strategies in open field.

Transcriptional and temporal response of Populus stems to gravi-stimulation

Plants modify development in response to external stimuli, to produce new growth that is appropriate for environmental conditions. For example, gravi-stimulation of leaning branches in angiosperm trees results in modifications of wood development, to produce tension wood that pulls leaning stems upright. Here, we use gravi-stimulation and tension wood response to dissect the temporal changes in gene expression underlying wood formation in Populus stems. Using time-series analysis of seven time points over a 14-d experiment, we identified 8,919 genes that were differentially expressed between tension wood (upper) and opposite wood (lower) sides of leaning stems. Clustering of differentially expressed genes showed four major transcriptional responses, including gene clusters whose transcript levels were associated with two types of tissue-specific impulse responses that peaked at about 24-48 h, and gene clusters with sustained changes in transcript levels that persisted until the end of the 14-d experiment. Functional enrichment analysis of those clusters suggests they reflect temporal changes in pathways associated with hormone regulation, protein localization, cell wall biosynthesis and epigenetic processes. Time-series analysis of gene expression is an underutilized approach for dissecting complex developmental responses in plants, and can reveal gene clusters and mechanisms influencing development.

The Arabidopsis homolog of human G3BP1 is a key regulator of stomatal and apoplastic immunity

Arabidopsis homolog of human G3BP1 negatively regulates plant immunity and defense responses. Atg3bp1 mutant lines show constitutive stomata closure, expression of a number of key defense marker genes, and accumulate salicylic acid.

Mammalian Ras-GTPase–activating protein SH3-domain–binding proteins (G3BPs) are a highly conserved family of RNA-binding proteins that link kinase receptor-mediated signaling to RNA metabolism. Mammalian G3BP1 is a multifunctional protein that functions in viral immunity. Here, we show that the Arabidopsis thaliana homolog of human G3BP1 negatively regulates plant immunity. Arabidopsis g3bp1 mutants showed enhanced resistance to the virulent bacterial pathogen Pseudomonas syringae pv. tomato. Pathogen resistance was mediated in Atg3bp1 mutants by altered stomatal and apoplastic immunity. Atg3bp1 mutants restricted pathogen entry into stomates showing insensitivity to bacterial coronatine–mediated stomatal reopening. AtG3BP1 was identified as a negative regulator of defense responses, which correlated with moderate up-regulation of salicylic acid biosynthesis and signaling without growth penalty.

Differential Suppression of Nicotiana benthamiana Innate Immune Responses by Transiently Expressed Pseudomonas syringae Type III Effectors

The plant pathogen Pseudomonas syringae injects about 30 different virulence proteins, so-called effectors, via a type III secretion system into plant cells to promote disease. Although some of these effectors are known to suppress either pattern-triggered immunity (PTI) or effector-triggered immunity (ETI), the mode of action of most of them remains unknown. Here, we used transient expression in Nicotiana benthamiana, to test the abilities of type III effectors of Pseudomonas syringae pv. tomato (Pto) DC3000 and Pseudomonas syringae pv. tabaci (Pta) 11528 to interfere with plant immunity. We monitored the sequential and rapid bursts of cytoplasmic Ca2+ and reactive oxygen species (ROS), the subsequent induction of defense gene expression, and promotion of cell death. We found that several effector proteins caused cell death, but independently of the known plant immune regulator NbSGT1, a gene essential for ETI. Furthermore, many effectors delayed or blocked the cell death-promoting activity of other effectors, thereby potentially contributing to pathogenesis. Secondly, a large number of effectors were able to suppress PAMP-induced defense responses. In the majority of cases, this resulted in suppression of all studied PAMP responses, suggesting that these effectors target common elements of PTI. However, effectors also targeted different steps within defense pathways and could be divided into three major groups based on their suppressive activities. Finally, the abilities of effectors of both Pto DC3000 and Pta 11528 to suppress plant immunity was conserved in most but not all cases. Overall, our data present a comprehensive picture of the mode of action of these effectors and indicate that most of them suppress plant defenses in various ways.

CaWRKY22 Acts as a Positive Regulator in Pepper Response to RalstoniaSolanacearum by Constituting Networks with CaWRKY6, CaWRKY27, CaWRKY40, and CaWRKY58

The WRKY web, which is comprised of a subset of WRKY transcription factors (TFs), plays a crucial role in the regulation of plant immunity, however, the mode of organization and operation of this network remains obscure, especially in non-model plants such as pepper (Capsicum annuum). Herein, CaWRKY22, a member of a subgroup of IIe WRKY proteins from pepper, was functionally characterized in pepper immunity against Ralstonia Solanacearum. CaWRKY22 was found to target the nuclei, and its transcript level was significantly upregulated by Ralstonia Solanacearum inoculation (RSI) and exogenously applied salicylic acid (SA), Methyl jasmonate (MeJA), or ethephon (ETH). Loss-of-function CaWRKY22, caused by virus-induced gene silencing (VIGS), enhanced pepper’s susceptibility to RSI. In addition, the silencing of CaWRKY22 perturbed the hypersensitive response (HR)-like cell death elicited by RSI and downregulated defense-related genes including CaPO2, CaPR4, CaACC, CaBPR1, CaDEF1, CaHIR1, and CaWRKY40. CaWRKY22 was found to directly bind to the promoters of CaPR1, CaDEF1, and CaWRKY40 by chromatin immuno-precipitation (ChIP) analysis. Contrastingly, transient overexpression of CaWRKY22 in pepper leaves triggered significant HR-like cell death and upregulated the tested immunity associated maker genes. Moreover, the transient overexpression of CaWRKY22 upregulated the expression of CaWRKY6 and CaWRKY27 while it downregulated of the expression of CaWRKY58. Conversely, the transient overexpression of CaWRKY6, CaWRKY27, and CaWRKY40 upregulated the expression of CaWRKY22, while transient overexpression of CaWRKY58 downregulated the transcript levels of CaWRKY22. These data collectively recommend the role of CaWRKY22 as a positive regulator of pepper immunity against R. Solanacearum, which is regulated by signaling synergistically mediated by SA, jasmonic acid (JA), and ethylene (ET), integrating into WRKY networks with WRKY TFs including CaWRKY6, CaWRKY27, CaWRKY40, and CaWRKY58.

CaWRKY40b in Pepper Acts as a Negative Regulator in Response to Ralstonia solanacearum by Directly Modulating Defense Genes Including CaWRKY40

WRKY transcription factors (TFs) have been implicated in plant growth, development, and in response to environmental cues; however, the function of the majority of pepper WRKY TFs remains unclear. In the present study, we functionally characterized CaWRKY40b, a homolog of AtWRKY40, in pepper immunity. Ralstonia solanacearum inoculation (RSI) in pepper plants resulted in downregulation of CaWRKY40b transcript, and green fluorescent protein (GFP)-tagged CaWRKY40b was localized to the nuclei when transiently overexpressed in the leaves of Nicotiana benthamiana. Virus-induced gene silencing (VIGS) of CaWRKY40b significantly decreased pepper’ susceptibility to RSI. Consistently, the transient over-expression of CaWRKY40b-SRDX (chimeric repressor version of CaWRKY40b) triggered cell death, as indicated by darker trypan blue and DAB staining. CaWRKY40b targets a number of immunity-associated genes, including CaWRKY40 JAR, RLK1, EIN3, FLS2, CNGIC8, CDPK13, and heat shock cognate protein 70 (HSC70), which were identified by ChIP-seq and confirmed using ChIP-real time PCR. Among these target genes, the negative regulator HSC70 was upregulated by transient overexpression of CaWRKY40b and downregulated by silencing of CaWRKY40b, whereas other positive regulators as well as two non-target genes, CaNPR1 and CaDEF1, were downregulated by the transient overexpression of CaWRKY40b and upregulated by CaWRKY40b silencing or transient overexpression of CaWRKY40b-SRDX. In addition, CaWRKY40b exhibited a positive feedback regulation at transcriptional level by directly targeting the promoter of itself. In conclusion, the findings of the present study suggest that CaWRKY40b acts as a negative regulator in pepper immunity against R. solanacearum by transcriptional modulation of a subset of immunity-associated genes; it also represses immunity in the absence of a pathogen, and derepresses immunity upon pathogen challenge.

Novel Fungal Pathogenicity and Leaf Defense Strategies Are Revealed by Simultaneous Transcriptome Analysis of Colletotrichum fructicola and Strawberry Infected by This Fungus

Colletotrichum fructicola, which is part of the C. gloeosporioides species complex, can cause anthracnose diseases in strawberries worldwide. However, the molecular interactions between C. fructicola and strawberry are largely unknown. A deep RNA-sequencing approach was applied to gain insights into the pathogenicity mechanisms of C. fructicola and the defense response of strawberry plants at different stages of infection. The transcriptome data showed stage-specific transcription accompanied by a step-by-step strawberry defense response and the evasion of this defense system by fungus. Fungal genes involved in plant cell wall degradation, secondary metabolism, and detoxification were up-regulated at different stage of infection. Most importantly, C. fructicola infection was accompanied by a large number of highly expressed effectors. Four new identified effectors function in the suppression of Bax-mediated programmed cell death. Strawberry utilizes pathogen-associated molecular patterns (PAMP)-triggered immunity and effector-triggered immunity to prevent C. fructicola invasion, followed by the initiation of downstream innate immunity. The up-regulation of genes related to salicylic acid provided evidence that salicylic acid signaling may serve as the core defense signaling mechanism, while jasmonic acid and ethylene pathways were largely inhibited by C. fructicola. The necrotrophic stage displayed a significant up-regulation of genes involved in reactive oxygen species activation. Collectively, the transcriptomic data of both C. fructicola and strawberry shows that even though plants build a multilayered defense against infection, C. fructicola employs a series of escape or antagonizing mechanisms to successfully infect host cells.

Calcium Signalling in Plant Biotic Interactions

Calcium (Ca2+) is a universal second messenger involved in various cellular processes, leading to plant development and to biotic and abiotic stress responses. Intracellular variation in free Ca2+ concentration is among the earliest events following the plant perception of environmental change. These Ca2+ variations differ in their spatio-temporal properties according to the nature, strength and duration of the stimulus. However, their conversion into biological responses requires Ca2+ sensors for decoding and relaying. The occurrence in plants of calmodulin (CaM) but also of other sets of plant-specific Ca2+ sensors such as calmodulin-like proteins (CMLs), Ca2+-dependent protein kinases (CDPKs) and calcineurin B-like proteins (CBLs) indicate that plants possess specific tools and machineries to convert Ca2+ signals into appropriate responses. Here, we focus on recent progress made in monitoring the generation of Ca2+ signals at the whole plant or cell level and their long distance propagation during biotic interactions. The contribution of CaM/CMLs and CDPKs in plant immune responses mounted against bacteria, fungi, viruses and insects are also presented.

Dickeya dadantii pectic enzymes necessary for virulence are also responsible for activation of the Arabidopsis thaliana innate immune system

Soft-rot diseases of plants attributed to Dickeya dadantii result from lysis of the plant cell wall caused by pectic enzymes released by the bacterial cell by a type II secretion system (T2SS). Arabidopsis thaliana can express several lines of defence against this bacterium. We employed bacterial mutants with defective envelope structures or secreted proteins to examine early plant defence reactions. We focused on the production of AtrbohD-dependent reactive oxygen species (ROS), callose deposition and cell death as indicators of these reactions. We observed a significant reduction in ROS and callose formation with a bacterial mutant in which genes encoding five pectate lyases (Pels) were disrupted. Treatment of plant leaves with bacterial culture filtrates containing Pels resulted in ROS and callose production, and both reactions were dependent on a functional AtrbohD gene. ROS and callose were produced in response to treatment with a cellular fraction of a T2SS-negative mutant grown in a Pels-inducing medium. Finally, ROS and callose were produced in leaves treated with purified Pels that had also been shown to induce the expression of jasmonic acid-dependent defence genes. Pel catalytic activity is required for the induction of ROS accumulation. In contrast, cell death observed in leaves infected with the wild-type strain appeared to be independent of a functional AtrbohD gene. It was also independent of the bacterial production of pectic enzymes and the type III secretion system (T3SS). In conclusion, the work presented here shows that D. dadantii is recognized by the A. thaliana innate immune system through the action of pectic enzymes secreted by bacteria at the site of infection. This recognition leads to AtrbohD-dependent ROS and callose accumulation, but not cell death.

GSHR, a Web-Based Platform Provides Gene Set-Level Analyses of Hormone Responses in Arabidopsis

Phytohormones regulate diverse aspects of plant growth and environmental responses. Recent high-throughput technologies have promoted a more comprehensive profiling of genes regulated by different hormones. However, these omics data generally result in large gene lists that make it challenging to interpret the data and extract insights into biological significance. With the rapid accumulation of theses large-scale experiments, especially the transcriptomic data available in public databases, a means of using this information to explore the transcriptional networks is needed. Different platforms have different architectures and designs, and even similar studies using the same platform may obtain data with large variances because of the highly dynamic and flexible effects of plant hormones; this makes it difficult to make comparisons across different studies and platforms. Here, we present a web server providing gene set-level analyses of Arabidopsis thaliana hormone responses. GSHR collected 333 RNA-seq and 1,205 microarray datasets from the Gene Expression Omnibus, characterizing transcriptomic changes in Arabidopsis in response to phytohormones including abscisic acid, auxin, brassinosteroids, cytokinins, ethylene, gibberellins, jasmonic acid, salicylic acid, and strigolactones. These data were further processed and organized into 1,368 gene sets regulated by different hormones or hormone-related factors. By comparing input gene lists to these gene sets, GSHR helped to identify gene sets from the input gene list regulated by different phytohormones or related factors. Together, GSHR links prior information regarding transcriptomic changes induced by hormones and related factors to newly generated data and facilities cross-study and cross-platform comparisons; this helps facilitate the mining of biologically significant information from large-scale datasets. The GSHR is freely available at http://bioinfo.sibs.ac.cn/GSHR/.

Transcriptome analysis, using RNA-Seq of Lomandra longifolia roots infected with Phytophthora cinnamomi reveals the complexity of the resistance response

The plant pathogen Phytophthora cinnamon the causal agent of disease in numerous species, is a major threat to natural vegetation and has economic impacts in agriculture. The pathogen principally invades the root system, which, in susceptible species, is rapidly colonised and functionally destroyed. Few species are resistant, however, where resistance is expressed the pathogen is restricted to small, localised lesions. The molecular mechanisms that underpin this response in resistant species are not well understood. Lomandra longifolia, an Australian native species, is highly resistant to P. cinnamomi. In an earlier study, we showed induction of resistance-related components such as callose, lignin and hydrogen peroxide (H₂ O₂ ) in L. longifolia roots that had been inoculated with P. cinnamomi. Here, in order to further identify, during the very early stages of infection, the molecular components and regulatory networks that may trigger resistance, a comprehensive root transcriptome analysis was performed using next generation sequencing. Overall, 18 cDNA libraries were produced generating 52.8 GB 126 base pair reads, which were de novo assembled into contigs. Differentially expressed genes (DEGs) were identified allowing the identification of infection-responsive candidate genes that were putatively related to resistance, and from this set ten were selected for qRT-PCR to validate the RNA-Seq expression value. Further analysis of individual candidates revealed that many were involved in PAMP-triggered immunity (PTI; pattern recognition receptors, glutathione S-transferase, callose synthases, pathogenesis-related protein-1, mitogen activated protein kinases) and effector-triggered immunity (ETI) (NBS-LRR, signalling genes, transcription factors and anti-pathogenic compound synthase genes). As these candidate genes or mediated components activate different defence signalling systems, they may have potential for investigation of novel approaches to disease control and in transgenic approaches for improvement, in susceptible species, of resistance to P. cinnamomi.

Gene expression trade-offs between defence and growth in English elm induced by Ophiostoma novo-ulmi

Wilt diseases caused by vascular pathogens include some of the most damaging stresses affecting trees. Dutch elm disease (DED), caused by the fungus Ophiostoma novo-ulmi, destroyed most of North American and European elm populations in the 20th century. The highly susceptible English elm, also known as Atinian clone, suffered the highest mortality rates during the last pandemic event, probably due to its lack of genetic diversity. To study the DED pathosystem, we inoculated English elm ramets with O. novo-ulmi and evaluated xylem anatomy, molecular response, and disease symptoms. The high DED susceptibility of the clone was linked to xylem structure. The transcript levels changed significantly for 1,696 genes during O. novo-ulmi invasion. Genes covering different steps of the plant immune system were identified, many of which showed homology with Arabidopsis thaliana genes involved in systemic acquired resistance. Induction of several pathogenesis-related proteins and repression of fasciclin-like arabinogalactan proteins and other cell wall biosynthesis pathways evidence unbalanced costs between growth and defence mechanisms far from the inoculation point. This study sheds light on elm molecular defence mechanisms against DED.

Histochemical Analyses Reveal That Stronger Intrinsic Defenses in Gossypium barbadense Than in G. hirsutum Are Associated With Resistance to Verticillium dahliae

Verticillium wilt, caused by Verticillium dahliae Kleb., is a serious threat to cotton (Gossypium spp.) crop production. To enhance our understanding of the plant's complex defensive mechanism, we examined colonization patterns and interactions between V. dahliae and two cotton species, the resistant G. barbadense and the susceptible G. hirsutum. Microscopic examinations and grafting experiments showed that the progression of infection was restricted within G. barbadense. At all pre- and postinoculation sampling times, levels of salicylic acid (SA) were also higher in that species than in G. hirsutum. Comparative RNA-Seq analyses indicated that infection induced dramatic changes in the expression of thousands of genes in G. hirsutum, whereas those changes were fewer and weaker in G. barbadense. Investigations of the morphological and biochemical nature of cell-wall barriers demonstrated that depositions of lignin, phenolic compounds, and callose were significantly higher in G. barbadense. To determine the contribution of a known resistance gene to these processes, we silenced GbEDS1 and found that the transformed plants had decreased SA production, which led to the upregulation of PLASMODESMATA-LOCATED PROTEIN (PDLP) 1 and PDLP6. This was followed by a decline in callose deposition in the plasmodesmata, which then led to increased pathogen susceptibility. This comparison between resistant and susceptible species indicated that both physical and chemical mechanisms play important roles in the defenses of cotton against V. dahliae.

Reactions of tobacco genotypes with different antioxidant capacities to powdery mildew and Tobacco mosaic virus infections

The interactions of powdery mildew (Golovinomyces orontii) and Tobacco mosaic virus (TMV) with tobacco lines having down or upregulated antioxidants were investigated. Xanthi-nc, its salicylic acid-deficient NahG mutant, a paraquat-sensitive Samsun (PS) and its paraquat tolerant (PT) mutant were used. Cell membrane damage caused by H2O2 was significantly higher in NahG than Xanthi, whereas it was lower in PT than in PS. Leakage of ions from PT was reduced by the powdery mildew infection. On the other hand TMV inoculation led to a 6-fold and 2-fold elevation of ion leakage from hypersensitive resistant NahG and Xanthi leaves, respectively, whereas ion leakage increased slightly from susceptible PS leaves. G. orontii infection induced ribonuclease (RNase) enzyme activity in extracts from Xanthi and NahG (about 200-250% increase) and weakly (about 20-30% increase) from PS and PT lines. Pre-treatment with protein kinase inhibitor staurosporine or protein phosphatase inhibitor okadaic acid very strongly inhibited mildew development on tobacco lines. Our experiments suggest that protein kinases inhibited by staurosporine seem to be important factors, while protein phosphatases inhibited by okadaic acid play less significant role in TMV-induced lesion development. Both powdery mildew and TMV infections up-regulated the expression of PR-1b, PR-1c and WRKY12 genes in all tobacco lines to various extents.

Functional and Promoter Analysis of ChiIV3, a Chitinase of Pepper Plant, in Response to Phytophthora capsici Infection

Despite the involvement of many members of the chitinase family in plant immunity, the precise functions of the majority of the members remain poorly understood. Herein, the gene ChiIV3 in Capsicum annuum encoding a chitinase protein containing a chitin binding domain and targeting to the plasma membrane was found to be induced by Phytophthora capsici inoculation (PCI) and applied chitin treatment. Besides its direct inhibitory effect on growth of Phytophthora capsici (P. capsici), ChiIV3 was also found by virus-induced gene silencing (VIGS) and transient overexpression (TOE) in pepper plants to act as a positive regulator of plant cell death and in triggering defense signaling and upregulation of PR (pathogenesis related) genes against PCI. A 5' deletion assay revealed that pChiIV3-712 to -459 bp was found to be sufficient for ChiIV3' response to PCI. Furthermore, a mutation assay indicated that W-box-466 to -461 bp in pChiIV3-712 to -459 bp was noted to be the PCI-responsible element. These results collectively suggest that ChiIV3 acts as a likely antifungal protein and as a receptor for unidentified chitin in planta to trigger cell death and defense signaling against PCI.

Transcriptome Analysis of the Sm-Mediated Hypersensitive Response to Stemphylium lycopersici in Tomato

Gray leaf spot disease caused by Stemphylium lycopersici is a major disease in cultivated tomato plants and threatens tomato-growing areas worldwide. Sm is a single dominant gene that confers resistance to tomato gray leaf spot disease agent. However, the underlying molecular mechanism remains unclear. Here, resistant (cv. Motelle, containing the Sm gene) and susceptible (cv. Moneymaker) plants were inoculated with virulent Stemphylium lycopersici isolate at a time point at which both cultivars showed a strong response to S. lycopersici infection. Transcriptome analyses were performed in both cultivars using RNA-seq. The number of differentially expressed genes (DEGs) was higher in Motelle than Moneymaker. Functional classification revealed that most DEGs were involved in plant-pathogen interactions, plant hormone signal transduction, regulation of autophagy, glycerophospholipid metabolism, and α-linolenic acid metabolism. Moreover, the genes that were significantly up-regulated in Sm tomatoes were involved in plant-pathogen interaction pathways. A total of 26 genes were selected for confirmation of differentially expressed levels by quantitative real-time PCR. This knowledge will yield new insights into the molecular mechanism of Sm responses to S. lycopersici infection.

Comparison between Proteome and Transcriptome Response in Potato (Solanum tuberosum L.) Leaves Following Potato Virus Y (PVY) Infection

Plant diseases caused by viral infection are affecting all major crops. Being an obligate intracellular organisms, chemical control of these pathogens is so far not applied in the field except to control the insect vectors of the viruses. Understanding of molecular responses of plant immunity is therefore economically important, guiding the enforcement of crop resistance. To disentangle complex regulatory mechanisms of the plant immune responses, understanding system as a whole is a must. However, integrating data from different molecular analysis (transcriptomics, proteomics, metabolomics, smallRNA regulation etc.) is not straightforward. We evaluated the response of potato (Solanum tuberosum L.) following the infection with potato virus Y (PVY). The response has been analyzed on two molecular levels, with microarray transcriptome analysis and mass spectroscopy-based proteomics. Within this report, we performed detailed analysis of the results on both levels and compared two different approaches for analysis of proteomic data (spectral count versus MaxQuant). To link the data on different molecular levels, each protein was mapped to the corresponding potato transcript according to StNIB paralogue grouping. Only 33% of the proteins mapped to microarray probes in a one-to-one relation and additionally many showed discordance in detected levels of proteins with corresponding transcripts. We discussed functional importance of true biological differences between both levels and showed that the reason for the discordance between transcript and protein abundance lies partly in complexity and structure of biological regulation of proteome and transcriptome and partly in technical issues contributing to it.

Large-scale transcriptome analysis reveals arabidopsis metabolic pathways are frequently influenced by different pathogens

Through large-scale transcriptional data analyses, we highlighted the importance of plant metabolism in plant immunity and identified 26 metabolic pathways that were frequently influenced by the infection of 14 different pathogens. Reprogramming of plant metabolism is a common phenomenon in plant defense responses. Currently, a large number of transcriptional profiles of infected tissues in Arabidopsis (Arabidopsis thaliana) have been deposited in public databases, which provides a great opportunity to understand the expression patterns of metabolic pathways during plant defense responses at the systems level. Here, we performed a large-scale transcriptome analysis based on 135 previously published expression samples, including 14 different pathogens, to explore the expression pattern of Arabidopsis metabolic pathways. Overall, metabolic genes are significantly changed in expression during plant defense responses. Upregulated metabolic genes are enriched on defense responses, and downregulated genes are enriched on photosynthesis, fatty acid and lipid metabolic processes. Gene set enrichment analysis (GSEA) identifies 26 frequently differentially expressed metabolic pathways (FreDE_Paths) that are differentially expressed in more than 60% of infected samples. These pathways are involved in the generation of energy, fatty acid and lipid metabolism as well as secondary metabolite biosynthesis. Clustering analysis based on the expression levels of these 26 metabolic pathways clearly distinguishes infected and control samples, further suggesting the importance of these metabolic pathways in plant defense responses. By comparing with FreDE_Paths from abiotic stresses, we find that the expression patterns of 26 FreDE_Paths from biotic stresses are more consistent across different infected samples. By investigating the expression correlation between transcriptional factors (TFs) and FreDE_Paths, we identify several notable relationships. Collectively, the current study will deepen our understanding of plant metabolism in plant immunity and provide new insights into disease-resistant crop improvement.