The biocontrol agent Pseudomonas chlororaphis PA23 primes Brassica napus defenses through distinct gene networks
BMC Genomics volume 18, Article number: 467 (2017)
The biological control agent Pseudomonas chlororaphis PA23 is capable of protecting Brassica napus (canola) from the necrotrophic fungus Sclerotinia sclerotiorum via direct antagonism. While we have elucidated bacterial genes and gene products responsible biocontrol, little is known about how the host plant responds to bacterial priming on the leaf surface, including global changes in gene activity in the presence and absence of S. sclerotiorum.
Application of PA23 to the aerial surfaces of canola plants reduced the number of S. sclerotiorum lesion-forming petals by 91.1%. RNA sequencing of the host pathogen interface showed that pretreatment with PA23 reduced the number of genes upregulated in response to S. sclerotiorum by 16-fold. By itself, PA23 activated unique defense networks indicative of defense priming. Genes encoding MAMP-triggered immunity receptors detecting flagellin and peptidoglycan were downregulated in PA23 only-treated plants, consistent with post-stimulus desensitization. Downstream, we observed reactive oxygen species (ROS) production involving low levels of H2O2 and overexpression of genes associated with glycerol-3-phosphate (G3P)-mediated systemic acquired resistance (SAR). Leaf chloroplasts exhibited increased thylakoid membrane structures and chlorophyll content, while lipid metabolic processes were upregulated.
In addition to directly antagonizing S. sclerotiorum, PA23 primes the plant defense response through induction of unique local and systemic defense networks. This study provides novel insight into the effects of biocontrol agents applied to the plant phyllosphere. Understanding these interactions will aid in the development of biocontrol systems as an alternative to chemical pesticides for protection of important crop systems.
Plants have evolved intricate defense mechanisms to thwart attack from devastating fungal pathogens. Morphological and structural barriers such as a waxy cuticle and tough cell wall are part of an innate defense mechanism against both living organisms and abiotic forces [1, 2]. Successfully bypassing these barriers causes activation of a defense response via detection of microbial- or pathogen-associated molecular patterns (MAMPs/PAMPs), or effector molecules [2, 3]. Activation of MAMP and PAMP pathways are often species-specific and occur through carefully orchestrated signal transduction networks. Following a targeted immune response to PAMPs/MAMPs, systemic resistance may be conferred through two distinct pathways: systemic acquired resistance (SAR) or induced systemic resistance (ISR). While ISR is activated upon colonization of plant roots by nonpathogenic rhizobacteria and fungi, SAR is induced by pathogen attack as well as other elicitors of targeted immune responses [4,5,6]. The specificity of the response is mediated by changes in gene expression and crosstalk between pathways [7,8,9]. Nonpathogenic organisms eliciting systemic resistance in plants boost host defense strategies leading to stronger and quicker reactions to future threats [5, 6, 10,11,12,13]. Such “priming” of plant defenses involves short- and long-term cellular changes .
Plants are exposed to high numbers of nonpathogenic organisms when biocontrol systems are employed. In addition to direct pathogen antibiosis, biocontrol agents (BCAs) can have the added benefit of priming host defenses [12, 15,16,17]. The majority of bacterial species that exhibit these dual roles are ISR-inducing Pseudomonas spp. or Bacillus spp. . For example, Pseudomonas putida WCS358 suppresses soil-borne pathogens through siderophore-mediated competition for iron, but can also induce ISR in Arabidopsis thaliana via host detection of flagellin, pseudobactin and lipopolysaccharides . Such microorganisms are good candidates to replace chemical pesticides, and the number of commercially available BCAs is steadily increasing [18, 20, 21]. In order to successfully implement BCAs in the field, a complete understanding of biocontrol system interactions, including their impact on the host plant, is required.
Brassica napus (canola) is an economically important crop of global significance. Despite attempts to breed cultivars with broad resistance traits, canola remains susceptible to a variety of pathogens. The necrotrophic fungus Sclerotinia sclerotiorum represents a particularly devastating pathogen to which no immune or highly resistant germplasm has been identified [22, 23]. As the causal agent of canola stem rot, this fungus can infect over 400 plant species worldwide . Such a wide host range has resulted in a heavy reliance on chemical pesticides for managing disease. Because of unwanted effects on the environment along with consumer safety concerns, biological control has emerged as an attractive alternative for crop protection.
Pseudomonas chlororaphis PA23 is a BCA capable of preventing S. sclerotiorum growth in vitro and in planta in susceptible B. napus cultivars [25, 26]. We have previously shown that PA23 directly antagonizes S. sclerotiorum through the excretion of antifungal metabolites including phenazines, pyrrolnitrin, proteases and lipases [27, 28], with pyrrolnitrin being the primary compound responsible for antagonism . Despite our growing understanding of the molecular mechanisms underlying PA23 antifungal activity, we have yet to understand how the presence of PA23 affects the host organism.
In the current study, we examine gene activity genome-wide in leaf tissues of B. napus in response to PA23 in the presence and absence of S. sclerotiorum. Global RNA sequencing (RNA-seq) was employed to identify differentially expressed genes indicative of B. napus defense responses. By itself, PA23 activated gene networks associated with defense priming. Moreover, changes in leaf architecture and increased chlorophyll content were observed in plants treated with PA23 alone. The presence of PA23 diminished S. sclerotiorum-induced defense pathways, including production of reactive oxygen species (ROS) and SAR induction. Collectively, these findings show that in addition to direct S. sclerotiorum antagonism, PA23 exerts a protective effect through host priming of defense networks.
P. chlororaphis PA23 reduces S. sclerotiorum infection rates in B. napus
To understand how B. napus responds to PA23 and how PA23 protects the host plant from S. sclerotiorum infection, we compared infection rates at the 30–50% flowering stage in the presence or absence of PA23. When comparing the rate of infection as the proportion of lesion-forming petals to total petals fallen onto the plant canopy, application of PA23 reduced the number of lesions by 91.1% (Fig. 1a) and sustained pathogen suppression for at least 7 days post treatment. In this infection model, leaf necrosis was visible under lesion-forming petals at 24 h post S. sclerotiorum infection in plants receiving the pathogen only treatment (Fig 1b).
Global patterns of gene expression in B. napus treated with combinations of PA23 and S. sclerotiorum
Next, we studied global patterns of gene activity using RNA-seq to better understand how B. napus responds to PA23 in the presence or absence of S. sclerotiorum. Principal component analysis (PCA) identified relationships between transcriptomes of plants receiving combinations of PA23 and S. sclerotiorum. Biological replicates for each treatment grouped together, with PA23-only treatments grouping closely to the water controls (Fig. 1c). Replicates representing treatment with S. sclerotiorum only (Ss) or PA23 + Ss clustered into distinct groups, with Ss clustering farthest from the water control (Fig. 1c).
An average of 81.42% of reads mapped to the B. napus cv. Darmor genome (Additional file 1: Table S1). Unmapped reads from all treatments except Ss were composed primarily (58.2%) of Brassicaceae, chloroplast, and mitochondrial sequences (Additional file 2: Figure S1). Unmapped reads from the Ss treatment are derived primarily (95.4%) from Sclerotiniaceae and other fungi. To identify transcripts that may be unique to B. napus cv. Westar we performed a transcriptome assembly of unmapped reads (Additional file 3: Dataset S1). The majority of assembled transcript fragments were identified as B. oleracea and suggest cv. Westar may have been recently outcrossed with this species (Additional file 4: Dataset S2).
Figure 1d summarizes mRNA detection and distribution of transcript abundance in treatment groups. A total of 48,454 genes with Fragments Per Kilobase of transcript per Million mapped reads (FPKM) ≥ 1 were detected across all samples, representing 48% of the predicted B. napus gene models . The number of B. napus transcripts detected was similar across treatments at an average of 43,007 expressed genes. We divided expression levels into low (FPKM ≥1, < 5), moderate (FPKM 5–25), and high (FPKM >25). The proportion of expressed genes falling into each category was similar across all treatment groups (39.0–43.5%, low; 35.3–44.3% moderate; 11.5–13.2%, high), where plants treated with S. sclerotiorum alone resulted in the greatest number of highly accumulating transcripts (Fig. 1d). In contrast, the PA23 + Ss treatment group had the greatest number of moderately and lowly accumulating transcripts across treatments.
Because we observed marked differences on the leaf surface when S. sclerotiorum was present with or without PA23, we compared differentially expressed genes (DEGs) of the treatment versus water control to identify similarities and differences at the RNA level. Figure 2a shows shared and specifically-upregulated DEGs in treatment groups compared to the water control. Plants treated with S. sclerotiorum alone had the greatest number of specifically-upregulated DEGs at 8237 genes. This trend held for both up- and downregulated genes (Additional file 5: Figure S2A). Plants treated with a combination of PA23 and S. sclerotiorum had the fewest specifically-upregulated DEGs at 515 genes, a 16-fold reduction compared to the S. sclerotiorum treatment group. The majority of upregulated DEGs observed in PA23 + Ss were shared with the pathogen only (Ss) group (3159 genes). While the number of upregulated DEGs in plants treated with PA23 alone was comparatively small (1361 genes), 556 genes were specific to PA23 treatment alone (Fig. 2a). DEGs of significance from the PA23 treated plants are listed in Table 1. Several markers of systemic acquired response (SAR) including pathogenesis-related proteins PR-1 (BnaC03G45470D) and PR-2 (BnaC08G28150D), lipid transporter protein DIR1 (BnaA03G11410D, BnaC03G14230D, BnaA10G09640D) and EARLI1 (BnaC03g29580D, BnaA09g20900D) are upregulated in the presence of PA23. All DEGs were specific to B. napus. When the sequence reads were mapped to the S. sclerotiorum genome, B. napus tissues treated with S. sclerotiorum alone showed an appreciable mapping rate compared to water control, PA23 treatment or PA23 + Ss (Additional file 1: Table S1).
To identify biological processes activated by the different treatments, we used the custom gene ontology (GO) term enrichment function of ChipEnrich with gene sets identified in Fig. 2a . Figure 2b summarizes GO terms of interest from this analysis. A heatmap of relevant enriched GO terms for genes downregulated in these groups is available in Additional file 5: Figure S2B and the complete list of GO terms with associated p-values and genes can be found in Additional file 6: Dataset S3. Response to ROS (log10 p-value < −4) was significantly enriched for all of the treatments. Subsets of genes belonging to this category accumulated in PA23 + Ss and Ss but not in plants treated with PA23 alone, while others were upregulated in all three treatment groups. In the latter, response to ROS involved upregulation of FERRETIN 1 (FER1, BnaC03G00160D), FERRETIN 3 (FER3, BnaC06G15730D), and two homologs of HEAT SHOCK TRANSCRIPTION FACTOR A4A (HSFA4A, BnaC01G11370D, BnaC03G62890D). In PA23 + Ss and Ss treatment groups, four homologs of ASCORBATE PEROXIDASE 1 (APX1, BnaA06G04380D, BnaA09G49190D, BnaC05G05550D, BnaC08G43490D), three homologs of ETHYLENE RESPONSIVE ELEMENT BINDING FACTOR 6 (ERF6, BnaA01G34910D, BnaA08G08300D, BnaC01G10080D), two homologs of XANTHINE DEHYDROGENASE 1 (XDH1, BnaA09G00610D, BnaCNNG46690D) and two additional homologs of HEAT SHOCK TRANSCRIPTION FACTOR A4A (HSFA4A, BnaANNG31620D, BnaC07G35520D) were identified. Similarly, all treatment groups were enriched for the SAR GO term, while the PA23 + Ss and Ss treatment groups were enriched for respiratory burst involved in defense response.
Because reinforcement of the plant cell wall is a marker for priming, it is noteworthy that the PA23 + Ss group was uniquely enriched for GO terms involving cell wall remodeling (phospholipid transfer to membrane, log10 p-value < −6; anchored component of membrane, log10 p-value < −5; cellular glucan metabolic process, log10 p-value < −4) (Fig. 2b). Overrepresented transcripts include lipid transfer proteins and xyloglucan endotransglycosylases. Furthermore, PA23 only and PA23 + Ss groups were both enriched for the proteinaceous extracellular matrix GO term (log10 p-value < −3).
Several GO terms associated with the chloroplast were overrepresented in PA23-treated plants (Fig. 2b). Specifically, thylakoid (log10 p-value < −7), plastid translation (log10 p-value < −4) and plastoglobule (log10 p-value < −3) were enriched in the PA23 only treatment group, and plastid thylakoid membrane (log10 p-value < −3) as well as negative regulation of chlorophyll catabolic process (log10 p-value < −6) were enriched in both PA23 only and PA23 + Ss treatment groups. Transcripts among these groups included subunits of photosystems I and II, plastid-specific ribosomal subunits, and two homologs of the negative regulator of chlorophyll degradation STAY-GREEN 2 (SGR2, BnaA03G24900D and BnaC03G72930D).
PA23 prevents the accumulation of ROS in the leaf
We stained leaf tissues for detection of hydrogen peroxide (H2O2) and superoxide (O2 −) radicals to help validate the RNA seq dataset. Whereas Ss-treated leaves retained both stains in the regions directly surrounding lesions, indicative of H2O2 and O2 − presence (Fig. 3m-p), ROS production was greatly reduced when plants were treated with PA23 prior to the fungal pathogen (Fig. 3i-l). PA23-treated leaves had no regions containing large aggregations of H2O2 and O2 −, although production of the former appeared similar to the PA23 + Ss treatment group suggesting PA23 activates a mild H2O2 reaction (Fig. 3g, h).
Dominant patterns of gene expression reveal treatment-specific responses in B. napus
To identify more complex patterns of expression, we clustered gene activity from all treatment groups into dominant patterns (DPs) of co-expression using a modified fuzzy K means clustering algorithm. We identified six DPs from this analysis (Additional file 7: Figure S3). We focused our attention on three DPs where genes accumulated specifically in response to one of the treatments (Fig. 4a). The number of genes clustering in these DPs ranged from 282 in DP5 to 11,340 in DP1. We identified genes associated with translation, response to fungus, plant-type hypersensitive response, and stomatal closure (log10 p-value < −8) when plants were infected with S. sclerotiorum without protection by PA23 (DP1). The chloroplast GO term is also significantly represented in DP1 (Fig. 4b; log10 p-value < −28), although no specific chloroplast-related processes were identified. In contrast, pre-treatment with PA23 (DP5) induced chloroplast-related components and processes, such as the chloroplast envelope, vitamin E biosynthesis, and starch metabolic processes (log10 p-value < −3). When plants were treated with PA23 alone (DP3), genes associated with GO terms for cytoskeletal changes and auxin-activated signaling were upregulated (log10 p-value < −6), as well as several chloroplast-related terms including plastid translation, thylakoid and chloroplast envelope. Some GO terms were common to more than one DP; in particular chloroplast (DP1 and DP3) and chloroplast envelope (DP3 and DP5), indicating differential expression of genes within these GO terms among treatments (Fig. 4b). A complete list of genes for each DP is available in Additional file 8: Dataset S4.
PA23 treatment results in structural and metabolic changes in the B. napus chloroplast
Since each DP identified a number of chloroplast component- and process-related GO terms, we decided to explore gene expression patterns within this organelle in more detail. The chloroplast GO term was significantly represented in more than one DP indicating differential expression of subsets of associated genes. Accordingly, we compared relative expression levels of genes belonging to chloroplast-related sub-regional and functional GO terms to reveal contrasting expression patterns (Fig.5a). For example, a comparison of relative abundance of genes for chlorophyll catabolite transmembrane transport reveals upregulation in the S. sclerotiorum-treated group. Although PA23 + Ss treatment caused many of the same genes to be expressed, transcript abundance was higher in the absence of PA23. It is possible that the reduced infection load associated with the PA23 + Ss group results in decreased transcript abundance. This trend was also observed for the chloroplast inner membrane and chloroplast stroma GO terms (Fig. 5a, II and IV). Genes downregulated in S. sclerotiorum -treated plants were enriched for chloroplast photosystem I & II GO terms as well as thylakoid-related GO terms, including genes encoding photosystem I subunits (PSAN, BnaA06g23190D and BnaC03g50210D; PSAG, BnaC06g07480D; PSAP, BnaC04g51600D), photosystem II subunits and regulatory proteins (PSBY, BnaA07g38700D and BnaC06g26560D; PSB27, BnaC08g46250D; PSBP-1, BnaC08g44890D) and other photosynthesis-related proteins (CRR23, BnaA07g28860D and BnaC06g31900D; PLASTOCYANIN 2, BnaA06g38550D; PNSL2, BnaA09g45770D). Overall upregulation of the chloroplast thylakoid and plastoglobule GO terms (Fig. 5a, V and VI) in the PA23 only group was also confirmed. A complete list of gene names and fold change values for genes in Fig. 5a is available in Additional file 9: Dataset S5.
Through transmission electron microscopy (TEM), PA23-mediated changes in chloroplast membrane structure were identified validating our transcriptional findings (Fig. 5b-f). Chloroplasts within the uppermost layer of the palisade mesophyll from the water-treated control group contained large starch granules and reduced thylakoid membrane structure. When plants were treated with PA23, the area dedicated to thylakoid structures increased with a concomitant increase in granal stack organization and the accumulation of plastoglobuli 24 and 48 h post inoculation (Fig. 5c, d). While gene expression in plants treated with PA23 + Ss indicated significant upregulation of starch metabolic processes (Fig. 4b), these chloroplasts were similar in appearance to those of the biocontrol only-treated group containing many grana stacks and plastoglobules (Fig. 5e). When plants were inoculated with S. sclerotiorum, the area dedicated to thylakoid structures within the chloroplast appeared to be reduced compared to other treatments (Fig. 5f).
Since structural changes were observed in PA23-treated chloroplasts, we sought to determine whether plant chlorophyll was also affected. Chlorophyll content between treatment groups showed concentrations of total chlorophyll (chlorophyll a and b) significantly increased by 22.8% when plants were treated with PA23. Conversely, plants exposed to PA23 and S. sclerotiorum or S. sclerotiorum alone showed no difference in chlorophyll content when compared to the water control (Fig. 5g). To clarify whether these changes were due to increased chlorophyll production or decreased chlorophyll degradation, we examined expression of genes encoding chlorophyllase. Global transcriptional analysis of genes involved in chlorophyll degradation pathways revealed significant downregulation of CHLOROPHYLLASE 1 (CLH1) in PA23-treated plants, which was confirmed via qRT-PCR (Fig. 5h). In contrast, S. sclerotiorum-treated plants showed significantly increased levels of CLH1 and decreased expression of genes associated with multiple chlorophyll-related pathways (Fig. 5h, Additional file 2: Figure S1B).
PA23 activates unique innate immunity and SAR networks to prime plant defenses
To understand how PA23 triggers plant priming mechanisms, we compared the differential accumulation of transcripts known to be involved in innate immunity among treatment groups. The interactions of significantly up- or down-regulated genes encoding products known to function in triggered immunity are presented in Fig. 6a. Treatment-specific expression levels are represented as a heat map for each gene and homolog. Many of these genes were upregulated in response to the pathogen, S. sclerotiorum. In response to PA23, two prokaryote-specific surface receptor genes were downregulated, namely FLS2 (BnaA09g17950D), encoding a receptor kinase which binds bacterial flagellin, and LYM3 (BnaCnng11350D), encoding a receptor required for detection of peptidoglycan.
Systemic resistance and defense priming are often the result of activation of these innate immunity networks. As such, we investigated downstream responses to discover how PA23 primes the host plant on a systemic level. The SAR GO term was enriched in all treatment groups, including PA23 alone (Fig. 2b). Analysis of RNA-seq data for genes associated with the SAR GO term revealed genes upregulated in the PA23 only group but downregulated in PA23 + Ss and Ss groups (Fig. 6b subgroup 1). Subgroup 1 contained DIR1, CPR5 and ACP4, suggesting that these genes are involved in PA23-induced SAR. We also observed genes induced by S. sclerotiorum only (Fig. 6b subgroup 2), as well as those showing altered expression in all treatment groups compared to the water control (Fig. 6b, subgroup 3). Treatment with S. sclerotiorum alone induced the expression of ALD1, GLIP1, PR4 and DOX1. The genes MES9, PR1, PR2 and CHI were upregulated in all treatment groups, although expression was higher in plants where S. sclerotiorum was present. The same trend could be observed for the downregulation of CPN60B1. Expression levels for several of the genes mentioned above (ALD1, CHI, DOX1, FMO1, PR1, and PR4) were confirmed using qRT-PCR (Additional file 10: Fig. S4).
Many nonpathogenic rhizobacteria are capable of priming plants for an enhanced defense response via induced systemic resistance (ISR) or systemic acquired resistance (SAR) (see 5,6 and  for reviews). However, few of these nonpathogenic microorganisms have been shown to directly antagonize pathogens, especially when applied to the plant phyllosphere . In the current study, we provide compelling evidence that in addition to direct antagonism of S. sclerotiorum mitigating disease progression, PA23 protects B. napus through priming of host defense systems at the RNA level.
The B. napus-S. sclerotiorum pathosystem has recently been transcriptionally profiled under a variety of conditions, providing new insight into genetic processes behind hormone and defense pathway signalling changes that support plant innate resistance [33, 34]. Patterns of gene expression revealed that SAR is weakly activated in plants when treated with the biocontrol agent PA23. The activation of an immune response in PA23-treated plants was reflected in both clustering analysis and in the number of DEGs among groups (Figs. 1c, 2a, Additional file 2: Figure S1A). Given the majority of DEGs in the PA23 + Ss group were shared with the Ss group, we hypothesize that S. sclerotiorum triggers many of the same defense mechanisms in B. napus regardless of PA23. It is important to note, however, that the magnitude of the response is lessened by the presence of the bacteria. This reduction likely results from PA23 antibiosis controlling the S. sclerotiorum infection. Additionally, PA23 may attenuate the host defense response through detection of bacterial effectors leading to priming.
Plant cellular signaling cascades which activate defense priming mechanisms are modulated by PAMP/MAMP-detecting pattern recognition receptors (PRRs). Wu al. (2016) observed differential activation of PRRs between resistant and susceptible lines of canola, highlighting the role of PRR activation in response to S. sclerotiorum. In plants treated with PA23 alone, surface detection of this nonpathogenic bacterium culminates in systemic defense stimulation via SAR. Genes associated with PAMP/MAMP- and DAMP-triggered immunity networks were turned on in response to S. sclerotiorum, but not PA23 (Fig. 6a). The downregulation of FLS2 in plants treated with PA23 alone is consistent with studies showing that FLS2 is degraded after initial ligand binding to prevent continuous stimulation . Thus, it is logical to surmise that detection of PA23 flagellin by the FLS2 receptor complex contributes to downstream induced defense processes. A similar mechanism may be responsible for repression of LYM3. Willmann et al.  found that in Arabidopsis, infection with virulent P. syringae pv. tomato DC3000 lead to strong downregulation of LYM3, FLS2 and to a lesser degree LYM1. LYM1 activity is unchanged in response to PA23. Differences in LYM1 and LYM3 expression may indicate receptor specificity, as these two proteins are not functionally redundant . The role these receptors play in the activation of systemic resistance appears to be variable. For example, flagellin plays an important role in the induction of systemic defense in response to enteric bacteria applied to leaves of Arabidopsis . In the rhizosphere, purified flagellin from nonpathogenic P. putida WCS358 induced ISR in Arabidopsis but not in bean or tomato . Moreover, Pseudomonas fluorescens SS101 induced systemic resistance in Arabidopsis fls2 mutants suggesting that FLS2 is not required for priming during this interaction . In the future, determining the roles of these receptor complexes in PA23-mediated defense priming should provide additional insight into the molecular underpinnings of this tri-partite system at the plant cell surface.
SAR is a classic example of defense priming which is usually associated with localized pathogen attack . While genes associated with SAR were induced in all treatments, we identified genes specifically upregulated in plants treated with PA23 alone. Elevated DIR1, CPR5, and ACP4 activity are indicative of SAR induction via glycerol-3-phosphate (G3P). G3P is one of several inducers of SAR, and is also a precursor for the synthesis of membrane and storage lipids [39, 40]. The lipid transfer protein, DIR1, which is upregulated in PA23-treated plants, leads to G3P accumulation and G3P-induced SAR . DIR1 binds to EARLI1, a paralog of AZI1 which is required for SAR and upregulated in PA23-only treated plants (Table 1) . Expression patterns of ACP4 further support this hypothesis, as ACP4 knockout plants are unable to perceive mobile SAR signals . ACP4 also plays a major role in the biosynthesis of fatty acids within the chloroplast, which indirectly increases G3P levels [41, 44]. While CPR5 is thought to have a role in inducing SAR, its exact function has not been defined. Notwithstanding, CPR5 has been shown to negatively regulate programmed cell death caused by effector-triggered immunity, which supports the hypothesis that PA23 promotes plant growth, as it is known that S. sclerotiorum infection results in cell death [45,46,47]. Taken together, our data support a role for G3P in PA23-mediated defense priming which serves to protect plants from fungal infection.
As a downstream marker of immune responses, ROS function as part of a localized hypersensitive-type reaction to invading pathogens, as well as signals for systemic response initiation . Sclerotinia acutely modifies the redox state of the host at infection sites to promote pathogenesis through production of oxalic acid . Several ROS response genes were identified that are upregulated in either all treatment groups (PA23, PA23 + Ss, Ss), or in those exclusively treated with S. sclerotiorum (Ss, PA23 + Ss), similar to the activation of redox homeostasis genes observed by Yang et al. . HSFA4A is thought to both regulate plant responses to oxidative stress and function as an antiapoptotic factor . Upregulation in PA23 groups suggests that treatment with bacteria alone can induce the plant oxidative stress response resulting in ROS production. This reaction appears to be mild in nature, as only two of the four upregulated HSFA4A homologs accumulated in the PA23 treatment group. In addition, HSFA4A-induced ASCORBATE PEROXIDASE 1 (APX1) was upregulated in S. sclerotiorum-treated plants, and is required for H2O2 scavenging and the prevention of protein oxidation during oxidative stress . Regulation of plant defense through ROS in the presence of PA23 may operate via the expression of ETHYLENE RESPONSE FACTOR 6 (ERF6) and XANTHINE DEHYDROGENASE 1 (XDH1). ERF6 modulates the expression of antioxidant enzymes to control ROS production and XDH1 is thought to be a source of O2 − production [48, 53]. An important difference between the PA23 + Ss and Ss groups is upregulation of superoxide dismutase activity in DP5, which suggests enhanced management of O2 − molecules in PA23 + Ss (Fig. 4a, b). Thus, gene expression changes indicate a mild and controlled production of ROS in response to PA23 alone, while pretreatment with PA23 prior to pathogen exposure allows for enhanced management of the oxidative stress response.
To substantiate the changes observed at the RNA level, we were interested to see if PA23 stimulates ROS production and/or modulates ROS levels on the leaf surface. In plants treated with PA23 alone, ROS production involving low levels of H2O2 was observed (Fig. 3). This is significant as non-toxic levels of H2O2 are key to the activation of priming pathways which reinforce resistance to abiotic and biotic stressors . Large regions staining for both H2O2 and O2 − molecules in Ss-treated leaves indicates a widespread oxidative stress response which appears to be unregulated. Leaves treated with PA23 + Ss had not only H2O2 but O2 − as well, confirming that S. sclerotiorum induces an oxidative stress response involving both molecules. Conversely, O2 − staining was not observed in the PA23 only treatment group (Fig. 3f, j). Collectively, our findings establish that ROS production in PA23-only treated leaves is limited to mild oxidative stress and is likely an outcome of weakly induced upstream innate immunity. Such findings support a role for ROS in the priming of plant basal defenses and initiation of long-distance signaling consistent with downstream PA23-mediated induction of SAR genes. In response to S. sclerotiorum, ROS production is limited in plants pretreated with PA23. We believe this results from a combination of direct fungal antagonism by PA23 and increased capacity to maintain ROS homeostasis during oxidative stress.
Thylakoid membranes in chloroplasts are a major source of ROS because they house the photosynthetic electron transport system . Oxidative stress observed during S. sclerotiorum infection is likely responsible for damage to the chloroplast and a reduction in photosynthetic activity [56, 57]. In Ss-treated plants, the observed upregulation of CLH1 is likely due to chloroplast damage leading to the release of chlorophylls from thylakoid membranes. CLH1 expression is induced for quick degradation of these photoactive molecules, a process known to be elicited by necrotrophic pathogens [58, 59]. TEM imaging substantiated these findings, as the region containing thylakoid tissue was decreased. In contrast, we observed an increase in relative area dedicated to thylakoid structures in chloroplasts within PA23-treated plants coupled with an increase in total chlorophyll content and repression of CLH1. In other studies, BCA application has been associated with increased chlorophyll in the plant host [60,61,62]. The lack of overexpression of CLH1 in the PA23-only treatment group is consistent with evidence showing that CLH1 is only induced after cellular collapse, a defense response not elicited by nonpathogenic organisms . Moreover, active repression of CLH1 may be evidence of a yet-to-be-identified role for chlorophyll in defense priming. Detection of more plastoglobules in the chloroplasts of PA23-treated plants suggests increased plastid lipid metabolism, potentially facilitating synthesis of lipid signaling molecules and remodeling of thylakoid membranes .
Such gene expression patterns together with the observed morphological changes in the chloroplasts of PA23-treated leaves correlate to the activation of SAR-modulated expression and mitigation of S. sclerotiorum disease progression.
Finally, the plant cell wall is an important component of innate immunity and defense priming . Changes in gene expression related to cell wall architecture were exclusively observed for plants in the PA23 + Ss treatment group. GO terms associated with cell wall development including thirteen xyloglucan endotransglycosylases which function to build up the cell wall during cell growth  and three homologs of LIPID TRANSFER PROTEIN 2 (LTP2) were enriched. These proteins may be involved in the deposition of cutin or wax in the extracellular matrix, and may increase tolerance to pathogens as has been shown in tobacco . As these genes were only upregulated in PA23 + Ss, we hypothesize this to be the result of a heightened defense response to S. sclerotiorum primed by PA23 pretreatment.
Findings from the current study broaden our understanding of PA23-mediated control of sclerotina stem rot. Beyond antibiosis, application of PA23 to the plant phyllosphere protects canola through induction of plant innate-immune pathways involving G3P-mediated SAR, ROS signaling molecules and protection of chloroplast integrity. Together, these processes serve to prime the plant and enhance defenses in response to fungal attack. As we move towards more sustainable approaches for crop disease management, it is essential that we fully appreciate the impact that BCAs pose on the plants they protect as well as the surrounding environment. Findings from the current study are an important step in this direction.
Plant and bacterial growth conditions
Brassica napus cv. Westar plants were grown in Sunshine Mix #1 soil in growth chambers at 21 °C with a light/dark photoperiod of 16 h/8 h and 0% humidity. P. chlororaphis PA23 was grown overnight in Luria-Bertani broth at 28 °C in a shaking incubator.
Greenhouse infection assays
One day prior to fungal pathogen exposure, B. napus plants at the 30% flowering stage were sprayed until dripping with a 2 × 108 cfu/mL solution of PA23 resuspended in sterile water supplemented with 0.02% Tween 20 as a surfactant. Plants not receiving biocontrol treatment (water control and Ss only groups) were sprayed with sterile water (0.02% Tween 20). Plants were sealed in clear bags to maintain relative humidity and returned to the growth chamber for 24 h. The following day, bags were removed and plants receiving the pathogen treatment were sprayed with an 8 × 104 spores/mL solution of S. sclerotiorum ascospores resuspended in sterile water (0.02% Tween 20). Control plants and plants to be exposed only to PA23 were sprayed with sterile water (0.02% Tween 20). Plants were transferred to a humidity chamber with humidity levels of 70–90% for 72 h. During this time, plants were gently shaken at the base twice to encourage petals to detach and fall into the plant canopy. Infection rates were quantified by calculating the ratio of petals causing lesions to total petals in the plant canopy. Counts from three plants were pooled for each treatment. This experiment was performed three times.
RNA extraction and sequencing
Infection assays were carried out as described above, and tissue from three biological replicates was collected for RNA extraction. Three leaves per plant and three plants per treatment group were used for each biological replicate. Leaves upon which petals had landed were used for collection, as these are sites of potential infection. The petal was removed from the leaf and approximately 1cm2 area of leaf tissue surrounding the site was collected with a scalpel. For S. sclerotiorum-infected leaves, green tissue immediately surrounding the lesion was collected. Cuttings were flash frozen in liquid nitrogen, and stored at −80 °C for no more than 2 days before processing. Total RNA was extracted using PureLink® Plant RNA Reagent (Invitrogen). DNA contamination was removed with the Turbo DNA-free™ kit (Ambion), following the manufacturer’s instructions. RNA concentration was verified using a NanoVue spectrophotometer (GE Healthcare), and quality was measured with an Agilent 2100 Bioanalyzer with Agilent RNA 6000 Pico and Nano Chips (Agilent Technologies; Santa Clara, CA, USA). RNA-seq libraries were prepared according to the alternative HTR protocol (C2) described by Kumar  with the exception of PCR enrichment of the libraries, where the number of cycles was adjusted to 11. Libraries were validated using the Agilent Bioanalyzer High Sensitivity DNA Assay with DNA chips (Agilent Technologies). The desired fragment sizes of sheared cDNA with ligated adapters were isolated employing the E-Gel® electrophoresis system (Invitrogen). 100 bp single-end RNA sequencing was carried out at Génome Québec (Montreal, Canada) on the Illumina HiSeq 2000 platform.
Sequenced reads were analyzed to remove barcode adapters and low quality reads using the Trimmomatic tool . The parameters for Trimmomatic which maximized mapping efficiency to the B. napus and S. sclerotiorum genomes (B. napus: v.4.1, Chalhoub et al. ; S. sclerotiorum: v1, Amselem et al., ) were determined using FastQC reports for quality control (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) followed by alignment with Tophat2 v.2.1.0 . Reads mapped to the B. napus and S. sclerotiorum genomes as expected, with 81.42% of reads mapping to B. napus across samples (Additional file 1: Table S1). Multiple mapping reads were retained to detect expression of plant defense and resistance genes of interest that may be duplicated within the genome. The value of multiple-mapping reads is determined by dividing the count by the total number of hit locations . For example, a read mapping to four locations will be recorded as 0.25 at each position. This correction statistically favors uniquely mapping reads in downstream analyses. Alignment of reads to these genomes was performed in high-sensitivity mode using B. napus reference annotation v5 from Chalhoub et al.  and the S. sclerotiorum reference annotation from Amselem et al.  as guides. The cufflinks and cuffmerge tools within the Cufflinks v.2.2.1 suite  were employed to construct a transcriptome from the reads and identify novel transcripts. Transdecoder (https://transdecoder.github.io) was used to identify open reading frames (ORFs) within transcript sequences. Genes were identified by aligning translated ORF sequences with proteins in the Arabidopsis TAIR10, NCBI and Uniprot databases using BLAST . Read counts were normalized to FPKM values using the cuffquant and cuffdiff tools in the Cufflinks package with default settings . Significantly differentially expressed genes were identified as those with a corrected p-value <0.05 (false discovery rate = 0.05). This output was used for hierarchical clustering via the pvclust package (https://cran.r-project.org/web/packages/pvclust/pvclust.pdf) and Venn diagram generation via Venny v2.1 (http://bioinfogp.cnb.csic.es/tools/venny/index.html). Dominant patterns (DPs) of expression were identified using the cuffdiff output data via the Fuzzy K-means (FKM) implementation FANNY (https://cran.r-project.org/web/packages/cluster/cluster.pdf) with a K value of 10. Transcripts with a Pearson’s correlation of 0.85 or above were assigned to DPs. Principal component analysis was performed on raw counts using the DESeq2 package . Unmapped reads were converted from BAM format to FASTA with the SAMtools software package (http://samtools.sourceforge.net/). Fasta files were then aligned to the SILVA phylogenetic database (https://www.arb-silva.de/) with a local instance of ncbi-BLAST-2.6.0+ (https://blast.ncbi.nlm.nih.gov/Blast.cgi) to identify the composition of unmapped reads (Additional file 2: Figure S1). To identify genes that may be unique to cv. Westar, unmapped reads from water and PA23 control samples were assembled to transcript fragments (Additional file 3: Dataset S1) with the Trinity package (https://github.com/trinityrnaseq/trinityrnaseq/wiki). Assembled transcripts were then aligned with a local instance of ncbi-BLAST-2.6.0+ to a database of Brassicaceae transcript sequences (Additional file 4: Dataset S2).
Staining for reactive oxygen species
B. napus leaves were stained for hydrogen peroxide (H2O2) and superoxide radical (O2 −) accumulation as per the methods of Kumar et al. . Briefly, leaves were severed with petioles intact and immersed in staining solution overnight in the dark, while avoiding contact between the solution and the severed petiole to reduce staining of the vein conduits. For H2O2 detection, a 1 mg/mL solution of 3,3′-diaminobenzidine was used. A 0.2% solution of nitrotetrazolium blue chloride in 50 mM sodium phosphate buffer at pH 7.5 enabled identification of O2 −. The following day, chlorophyll was removed from the leaves by submersion in 95% ethanol and heating in a boiling water bath.
Concentrations of chlorophyll a and b were determined based on methods described by Arnon  and Porra . Briefly, tissue from two leaves per plant was combined as one biological replicate and ground in liquid nitrogen. Approximately 200 mg of powder was extracted in 80% acetone buffered to pH 7.8. All steps of the extraction were performed in the dark. Concentrations were calculated as described by Porra . Three biological replicates were collected and absorbance values were measured twice for each replicate. The experiment was repeated twice.
Transmission electron microscopy
Leaf tissue was collected as above and processed following the methods of Chan and Belmonte . Tissue was fixed overnight in 3% glutaraldehyde in 0.025 M cacodylate buffer supplemented with 5 mM calcium chloride (pH 7.0). Plant material was rinsed with cacodylate buffer and post-fixed with 2% osmium tetroxide in 0.8% KFe(CN)6 in cacodylate buffer. After post-fixation, the sample was rinsed with distilled water and stained overnight with a 0.5% aqueous uranyl acetate solution. Plant material was rinsed in distilled water and dehydrated in a graded ethanol series. The leaf tissue was further dehydrated in 1:1 absolute ethanol to propylene oxide (v:v) followed by 100% propylene oxide. Finally, tissue was gradually infiltrated and embedded in Spurr’s epoxy resin at 70 °C. All methods prior to embedding were performed at 4 °C. Using a Reichert–Jung Ultracut ultramicrotome, sections were cut (90 nm thickness) with a Diatome diamond knife and mounted on copper grids. The sections were visualized with a Hitachi H-7000 transmission electron microscope at 75 kV and pictures were taken using AMT Image Capture Engine version 601.384.
Differential gene expression verification using qRT-PCR
RNA collected for RNA sequencing was also used to synthesize cDNA for qRT-PCR. qRT-PCR primers are listed in Additional file 11: Table S2. Following integrity checks via the Agilent 2100 Bioanalyzer with Agilent RNA 6000 Pico and Nano Chips (Agilent Technologies; Santa Clara, CA, USA), cDNA was synthesized with the Maxima First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, Inc.). qRT-PCR was performed on a Bio-Rad CFX Connect™ Real-Time System with SsoFast™ EvaGreen® Supermix (Bio-Rad, USA). The following cycling conditions were employed: 95 °C for 30 s, followed by 45 cycles of 95 °C for 2 s and 60 °C for 5 s. Melt curves (0.5 °C increments in a 68–90 °C range) for each gene were performed to assess the sample for non-specific targets, splice variants and primer dimers. CFX Manager v3.1 (Bio-Rad, USA) was used to calculate relative mRNA abundance from three technical replicates per primer set per sample utilizing the ΔΔCt method. ATGP4 was incorporated as an endogenous control and water-treated tissue was used as a reference sample. Three biological replicates for each treatment were measured.
Differentially expressed genes
Dominant pattern; transmission electron microscopy (TEM)
Fragments Per Kilobase of transcript per Million mapped reads
- H2O2 :
Induced systemic resistance
Microbial-associated molecular patterns
- O2 − :
Pathogen-associated molecular patterns
Principal component analysis
Plant growth-promoting rhizobacteria
PAMP/MAMP-detecting pattern recognition receptors
Reactive oxygen species
Systemic acquired resistance
The Arabidopsis Information Resource
Atkinson NJ, Urwin PE. The interaction of plant biotic and abiotic stresses: from genes to the field. J Exp Bot. 2012;63:3523–43.
Muthamilarasan M, Prasad M. Plant innate immunity: an updated insight into defense mechanism. J Biosci. 2013;38:433–49.
Jones JDG, Dangl JL. The plant immune system. Nature. 2006;444:323–9.
Ramamoorthy V, Viswanathan R, Raguchander T, Prakasam V, Samiyappan R. Induction of systemic resistance by plant growth promoting rhizobacteria in crop plants against pests and diseases. Crop Prot. 2001;20:1–11.
Pieterse CMJ, Zamioudis C, Berendsen RL, Weller DM, Van Wees SCM, Bakker PAHM. Induced systemic resistance by beneficial microbes. Annu Rev Phytopathol. 2014;52:347–75.
van Loon LC, Bakker PAHM, Pieterse CM. Systemic resistance induced by rhizosphere bacteria. Annu Rev Phytopathol. 1998;36:453–83.
Derksen H, Rampitsch C, Daayf F. Signaling cross-talk in plant disease resistance. Plant Sci Elsevier Ireland Ltd. 2013;207:79–87.
De Vos M, Van Oosten VR, Van Poecke RMP, Van Pelt JA, Pozo MJ, Mueller MJ, et al. Signal signature and transcriptome changes of Arabidopsis during pathogen and insect attack. Mol Plant-Microbe Interact. 2005;18:923–37.
Glazebrook J, Chen W, Estes B, Chang H-S, Nawrath C, Métraux J-P, et al. Topology of the network integrating salicylate and jasmonate signal transduction derived from global expression phenotyping. Plant J. 2003;34:217–28.
Wei G, Kloepper JW, Tuzun S. Induction of systemic resistance of cucumber to Colletotrichum orbiculare by select strains of plant growth-promoting rhizobacteria. Phytopathology. 1991;81:1508–12.
Van Peer R, Niemann GJ, Schippers S. Induced resistance and phytoalexin accumulation in biological control of Fusarium wilt of carnation by Pseudomonas sp. strain WCS417r. Phytopathology. 1991;81:728–34.
Shoresh M, Harman GE, Mastouri F. Induced systemic resistance and plant responses to fungal biocontrol agents. Annu Rev Phytopathol. 2010;48:21–43.
De Vleesschauwer D, Höfte M. Chapter 6 Rhizobacteria-induced systemic resistance. Adv Bot Res. 2009;51:223–81.
Pastor V, Luna E, Mauch-Mani B, Ton J, Flors V. Primed plants do not forget. Environ. Exp. Bot. 2013;94:46–56.
Hase S, Takahashi S, Takenaka S, Nakaho K, Arie T, Seo S, et al. Involvement of jasmonic acid signalling in bacterial wilt disease resistance induced by biocontrol agent Pythium oligandrum in tomato. Plant Pathol. 2008;57:870–6.
Pieterse CM, van Wees SC, Hoffland E, van Pelt JA, van Loon LC. Systemic resistance in Arabidopsis induced by biocontrol bacteria is independent of salicylic acid accumulation and pathogenesis-related gene expression. Plant Cell. 1996;8:1225–37.
Kurkcuoglu S, Degenhardt J, Lensing J, Al-Masri AN, Gau AE. Identification of differentially expressed genes in Malus domestica after application of the non-pathogenic bacterium Pseudomonas fluorescens Bk3 to the phyllosphere. J Exp Bot. 2006;58:733–41.
Pal KK, McSpadden GB. Biological control of plant pathogens. Instr: Plant Heal; 2006.
Meziane H, Van Der Sluis I, Van Loon LC, Höfte M, Bakker PAHM. Determinants of Pseudomonas putida WCS358 involved in inducing systemic resistance in plants. Mol Plant Pathol. 2005;6:177–85.
Bhattacharyya PN, Jha DK. Plant growth-promoting rhizobacteria (PGPR): emergence in agriculture. World J Microbiol Biotechnol. 2012;28:1327–50.
Fravel DR. Commercialization and implementation of biocontrol. Annu Rev Phytopathol. 2005;43:337–59.
Wu J, Cai G, Tu J, Li L, Liu S, Luo X, et al. Identification of QTLs for resistance to sclerotinia stem rot and BnaC.IGMT5.A as a candidate gene of the major resistant QTL SRC6 in Brassica napus. PLoS One. 2013;8:e67740.
Li CX, Li H, Sivasithamparam K, Fu TD, Li YC, Liu SY, et al. Expression of field resistance under western Australian conditions to Sclerotinia sclerotiorum in Chinese and Australian Brassica napus and Brassica juncea germplasm and its relation with stem diameter. Aust J Agric Res. 2006;57:1131.
Boland GJ, Hall R. Index of plant hosts of Sclerotinia sclerotiorum. Can J Plant Pathol. 1994;16:93–108.
Savchuk S, Fernando WGD. Effect of timing of application and population dynamics on the degree of biological control of Sclerotinia sclerotiorum by bacterial antagonists. FEMS Microbiol Ecol. 2004;49:379–88.
Fernando WGD, Nakkeeran S, Zhang Y, Savchuk S. Biological control of Sclerotinia sclerotiorum (lib.) de Bary by Pseudomonas and Bacillus species on canola petals. Crop Prot. 2007;26:100–7.
Poritsanos N, Selin C, Fernando WGD, Nakkeeran S, de Kievit TR. A GacS deficiency does not affect Pseudomonas chlororaphis PA23 fitness when growing on canola, in aged batch culture or as a biofilm. Can J Microbiol. 2006;52:1177–88.
Zhang Y, Fernando W, de Kievit TR, Berry C, Daayf F, Paulitz T. Detection of antibiotic-related genes from bacterial biocontrol agents with polymerase chain reaction. Can J Microbiol. 2006;52:476–81.
Selin C, Habibian R, Poritsanos N, Athukorala SNP, Fernando D, de Kievit TR. Phenazines are not essential for Pseudomonas chlororaphis PA23 biocontrol of Sclerotinia sclerotiorum, but do play a role in biofilm formation. FEMS Microbiol Ecol. 2010;71:73–83.
Chalhoub B, Denoeud F, Liu S, Parkin IAP, Tang H, Wang X, et al. Plant genetics. Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome. Science. 2014;345:950–3.
Belmonte MF, Kirkbride RC, Stone SL, Pelletier JM, Bui AQ, Yeung EC, et al. Comprehensive developmental profiles of gene activity in regions and subregions of the Arabidopsis seed. Proc Natl Acad Sci U S A. 2013;110:E435–44.
Van Wees SCM, Van der Ent S, Pieterse CMJ. Plant immune responses triggered by beneficial microbes. Curr Opin Plant Biol. 2008;11:443–8.
Wu J, Zhao Q, Yang Q, Liu H, Li Q, Yi X, et al. Comparative transcriptomic analysis uncovers the complex genetic network for resistance to Sclerotinia sclerotiorum in Brassica napus. Sci Rep. 2016;6:19007.
Zhao J, Buchwaldt L, Rimmer SR, Sharpe A, McGregor L, Bekkaoui D, et al. Patterns of differential gene expression in Brassica napus cultivars infected with Sclerotinia sclerotiorum. Mol Plant Pathol. 2009;10:635–49.
Smith JM, Salamango DJ, Leslie ME, Collins CA, Heese A. Sensitivity to Flg22 is modulated by ligand-induced degradation and de novo synthesis of the endogenous flagellin-receptor FLAGELLIN-SENSING2. Plant Physiol. 2014;164:440–54.
Willmann R, Lajunen HM, Erbs G, Newman M-A, Kolb D, Tsuda K, et al. Arabidopsis lysin-motif proteins LYM1 LYM3 CERK1 mediate bacterial peptidoglycan sensing and immunity to bacterial infection. Proc Natl Acad Sci U S A. 2011;108:19824–9.
Seo S, Matthews KR. Influence of the plant defense response to Escherichia coli O157:H7 cell surface structures on survival of that enteric pathogen on plant surfaces. Appl Environ Microbiol. 2012;78:5882–9.
van de Mortel JE, de Vos RCH, Dekkers E, Pineda A, Guillod L, Bouwmeester K, et al. Metabolic and transcriptomic changes induced in Arabidopsis by the rhizobacterium Pseudomonas fluorescens SS101. Plant Physiol. 2012;160:2173–88.
Venugopal SC, Chanda B, Vaillancourt L, Kachroo A, Kachroo P. The common metabolite glycerol-3-phosphate is a novel regulator of plant defense signaling. Plant Signal Behav. 2009;4:746–9.
Shah J, Zeier J. Long-distance communication and signal amplification in systemic acquired resistance. Front Plant Sci. 2013;4:30.
Yu K, Soares JM, Mandal MK, Wang C, Chanda B, Gifford AN, et al. A feedback regulatory loop between G3P and lipid transfer proteins DIR1 and AZI1 mediates azelaic-acid-induced systemic immunity. Cell Rep Elsevier. 2013;3:1266–78.
Cecchini NM, Steffes K, Schläppi MR, Gifford AN, Greenberg JT. Arabidopsis AZI1 family proteins mediate signal mobilization for systemic defence priming. Nat Commun. 2015;6:7658.
Xia Y, Gao QM, Yu K, Lapchyk L, Navarre D, Hildebrand D, et al. An intact cuticle in distal tissues is essential for the induction of systemic acquired resistance in plants. Cell Host Microbe. 2009;5:151–65.
Branen JK, Shintani DK, Engeseth NJ. Expression of antisense acyl carrier protein-4 reduces lipid content in Arabidopsis leaf tissue. Plant Physiol. 2003;132:748–56.
Bowling SA, Clarke JD, Liu Y, Klessig DF, Dong X. The cpr5 mutant of Arabidopsis expresses both NPR1-dependent and NPR1-independent resistance. Plant Cell. 1997;9:1573–84.
Kirik V, Bouyer D, Schöbinger U, Bechtold N, Herzog M, Bonneville J-M, et al. CPR5 is involved in cell proliferation and cell death control and encodes a novel transmembrane protein. Curr Biol. 2001;11:1891–5.
Wang S, Gu Y, Zebell SG, Anderson LK, Wang W, Mohan R, et al. A noncanonical role for the CKI-RB-E2F cell-cycle signaling pathway in plant effector-triggered immunity. Cell Host Microbe. 2014;16:787–94.
Sewelam N, Kazan K, Thomas-Hall SR, Kidd BN, Manners JM, Schenk PM. Ethylene response factor 6 is a regulator of reactive oxygen species signaling in Arabidopsis. PLoS One. 2013;8:e70289.
Kim KS, Min J-Y, Dickman MB. Oxalic acid is an elicitor of plant programmed cell death during Sclerotinia sclerotiorum disease development. Mol Plant-Microbe Interact. 2008;21:605–12.
Yang B, Srivastava S, Deyholos MK, Kav NNV. Transcriptional profiling of canola (Brassica napus L.) responses to the fungal pathogen Sclerotinia sclerotiorum. Plant Sci. 2007;173:156–71.
Pérez-Salamó I, Papdi C, Rigó G, Zsigmond L, Vilela B, Lumbreras V, et al. The heat shock factor A4A confers salt tolerance and is regulated by oxidative stress and the mitogen-activated protein kinases MPK3 and MPK6. Plant Physiol. 2014;165:319–34.
Davletova S, Rizhsky L, Liang H, Shengqiang Z, Oliver DJ, Coutu J, et al. Cytosolic ascorbate peroxidase 1 is a central component of the reactive oxygen gene network of Arabidopsis. Plant Cell. 2005;17:268–81.
Zarepour M, Kaspari K, Stagge S, Rethmeier R, Mendel RR, Bittner F. Xanthine dehydrogenase AtXDH1 from Arabidopsis thaliana is a potent producer of superoxide anions via its NADH oxidase activity. Plant Mol Biol. 2010;72:301–10.
Hossain MA, Bhattacharjee S, Armin S-M, Qian P, Xin W, Li H-Y, et al. Hydrogen peroxide priming modulates abiotic oxidative stress tolerance: insights from ROS detoxification and scavenging. Front Plant Sci. 2015;6:420.
Foyer CH, Lelandais M, Kunert KJ. Photooxidative stress in plants. Physiol Plant. 1994;92:696–717.
Seay M, Hayward AP, Tsao J, Dinesh-Kumar SP. Something old, something new: plant innate immunity and autophagy. Curr Top Microbiol Immunol. 2009;335:287–306.
Kabbage M, Williams B, Dickman MB. PLoS Pathog. 2013;9:e1003287.
Kariola T, Brader G, Li J, Palva ET. Chlorophyllase 1, a damage control enzyme, affects the balance between defense pathways in plants. Plant Cell. 2005;17:282–94.
Takamiya K, Tsuchiya T, Ohta H. Degradation pathway(s) of chlorophyll: what has gene cloning revealed? Trends Plant Sci. 2000;5:426–31.
Zhang F, Ge H, Zhang F, Guo N, Wang Y, Chen L, et al. Biocontrol potential of Trichoderma harzianum isolate T-aloe against Sclerotinia sclerotiorum in soybean. Plant Physiol Biochem. 2016;100:64–74.
Kumar S, Chauhan PS, Agrawal L, Raj R, Srivastava A, Gupta S, et al. Paenibacillus lentimorbus inoculation enhances tobacco growth and extenuates the virulence of Cucumber mosaic virus. PLoS One. 2016;11:e0149980.
Solanki MK, Singh RK, Srivastava S, Kumar S, Kashyap PL, Srivastava AK. Characterization of antagonistic-potential of two Bacillus strains and their biocontrol activity against Rhizoctonia solani in tomato. J Basic Microbiol. 2015;55:82–90.
Hu X, Makita S, Schelbert S, Sano S, Ochiai M, Tsuchiya T, et al. Reexamination of chlorophyllase function implies its involvement in defense against chewing herbivores. Plant Physiol. 2015;167:660–70.
Rottet S, Besagni C, Kessler F. The role of plastoglobules in thylakoid lipid remodeling during plant development. Biochim Biophys Acta. 1847;2015:889–99.
Luna E, Pastor V, Robert J, Flors V, Mauch-Mani B, Ton J. Callose deposition: a multifaceted plant defense response. Mol Plant-Microbe Interact. 2011;24:183–93.
Xu W, Campbell P, Vargheese AK, Braam J. The Arabidopsis XET-related gene family: environmental and hormonal regulation of expression. Plant J. 1996;9:879–89.
Molina A, Garcia-Olmedo F. Enhanced tolerance to bacterial pathogens caused by the transgenic expression of barley lipid transfer protein LTP2. Plant J. 1997;12:669–75.
Kumar R, Ichihashi Y, Kimura S, Chitwood DH, Headland LR, Peng J, et al. A high-throughput method for Illumina RNA-Seq library preparation. Front Plant Sci. 2012;3:202.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.
Amselem J, Cuomo CA, van Kan JAL, Viaud M, Benito EP, Couloux A, et al. Genomic analysis of the necrotrophic fungal pathogens Sclerotinia sclerotiorum and Botrytis Cinerea. PLoS Genet. 2011;7:e1002230.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and cufflinks. Nat Protoc. 2012;7:562–78.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
Arnon DI. Copper enzymes in isolated chloroplasts. Polyphenoloxidase in Beta vulgaris. Plant Physiol. 1949;24:1–15.
Porra RJ. The chequered history of the development and use of simultaneous equations for the accurate determination of chlorophylls a and b. Photosynth Res. 2002;73:149–56.
Chan A, Belmonte MF. Histological and ultrastructural changes in canola (Brassica napus) funicular anatomy during the seed lifecycle. Bot NRC Res Press. 2013;91:671–9.
Financial support for this work was provided by the Natural Sciences and Engineering Research Council (NSERC; TdK) as well as a Manitoba Agriculture Rural Development Growing Forward 2 grant and a Canola Agronomic Research Program (CARP) grant through Canola Council of Canada (TdK, MFB and WGDF). KAD and MGB were supported by NSERC Scholarships and IJG through a Manitoba Graduate Scholarship.
Availability of data and materials
The datasets generated during the current study are available in the GEO repository, accession record GSE84798 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE84798).
KAD, TdK, MFB and WGDF conceived and designed the study. MGB synthesized cDNA libraries and prepared samples for sequencing. MGB and IJG provided scripts and advice for RNA-seq data analysis. JLM prepared tissue samples for transmission electron microscopy and captured images. KAD completed all other experiments and all of the data analysis, and was the major contributor in writing this manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Brassica napus cv. Westar seeds were obtained from the laboratory of WGDF.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
RNA-seq library reads mapped to the Brassica napus and Sclerotinia sclerotiorum genomes. (PDF 14.7 kb)
Identification of contamination within unmapped reads with ncbi-BLAST. For all treatments, unmapped reads were converted to fasta format with SAMtools and aligned to the SILVA rRNA sequence database. Average distribution of alignments are represented in individual pie charts for each treatment group. (PDF 10 kb)
Dataset S1. Trinity transcriptome assembly of unmapped reads from control treatment groups. (FASTA 51567 kb)
Dataset S2. ncbi-BLAST results of unmapped transcripts assembled with Trinity to Brassicaceae nucleotide database. (XLSX 3895 kb)
Downregulated genes in canola treatment groups. A. Venn diagram of B. napus gene counts for uniquely and significantly downregulated genes in treatment groups compared to the water control. B. Heatmap of enriched GO terms selected from genes identified in A. (PDF 347 kb).
Dataset S3. GO terms enriched in shared and uniquely-expressed genes between treatment groups, with associated log10 p-values for each treatment. List of GO terms was generated via analysis of lists of genes grouped from Fig. 2a in CanEnrich software. Red-highlighted cells indicate significant p-values of less than 0.001. (XLSX 968 kb)
Dominant patterns of expression generated from RNA-seq FPKM values using a Fuzzy K-means clustering algorithm. (PDF 293 kb)
Dataset S4. List of genes comprising each dominant pattern (DP) of expression seen in Additional file 2: Figure S2. (XLSX 291 kb)
Dataset S5. Gene names, log2 fold change values and corresponding gene ontology (GO) term for select chloroplast-related GO terms. (XLSX 29 kb)
Relative abundance of select SAR-related gene transcripts as determined by RNA-seq (grey bars) and qRT-PCR (black line). (PDF 100 kb)
qRT-PCR primers used in this study. (PDF 14.7 kb)
About this article
Cite this article
Duke, K.A., Becker, M.G., Girard, I.J. et al. The biocontrol agent Pseudomonas chlororaphis PA23 primes Brassica napus defenses through distinct gene networks. BMC Genomics 18, 467 (2017). https://doi.org/10.1186/s12864-017-3848-6