ORIGINAL RESEARCH article
Phaeophyceaean (Brown Algal) Extracts Activate Plant Defense Systems in Arabidopsis thaliana Challenged With Phytophthora cinnamomi
- 1School of Life and Environmental Sciences, Deakin University, Waurn Ponds Campus, Geelong, VIC, Australia
- 2Department of Plant Pathology, Faculty of Agriculture, Sher-e-Bangla Agricultural University, Dhaka, Bangladesh
- 3Seasol International R&D Department, Bayswater, VIC, Australia
Seaweed extracts are important sources of plant biostimulants that boost agricultural productivity to meet current world demand. The ability of seaweed extracts based on either of the Phaeophyceaean species Ascophyllum nodosum or Durvillaea potatorum to enhance plant growth or suppress plant disease have recently been shown. However, very limited information is available on the mechanisms of suppression of plant disease by such extracts. In addition, there is no information on the ability of a combination of extracts from A. nodosum and D. potatorum to suppress a plant pathogen or to induce plant defense. The present study has explored the transcriptome, using RNA-seq, of Arabidopsis thaliana following treatment with extracts from the two species, or a mixture of both, prior to inoculation with the root pathogen Phytophthora cinnamomi. Following inoculation, five time points (0−24 h post-inoculation) that represented early stages in the interaction of the pathogen with its host were assessed for each treatment and compared with their respective water controls. Wide scale transcriptome reprogramming occurred predominantly related to phytohormone biosynthesis and signaling, changes in metabolic processes and cell wall biosynthesis, there was a broad induction of proteolysis pathways, a respiratory burst and numerous defense-related responses were induced. The induction by each seaweed extract of defense-related genes coincident with the time of inoculation showed that the plants were primed for defense prior to infection. Each seaweed extract acted differently in inducing plant defense-related genes. However, major systemic acquired resistance (SAR)-related genes as well as salicylic acid-regulated marker genes (PR1, PR5, and NPR1) and auxin associated genes were found to be commonly up-regulated compared with the controls following treatment with each seaweed extract. Moreover, each seaweed extract suppressed P. cinnamomi growth within the roots of inoculated A. thaliana by the early induction of defense pathways and likely through ROS-based signaling pathways that were linked to production of ROS. Collectively, the RNA-seq transcriptome analysis revealed the induction by seaweed extracts of suites of genes that are associated with direct or indirect plant defense in addition to responses that require cellular energy to maintain plant growth during biotic stress.
Plants have evolved marvelous interactive and adaptive systems to grow in challenging and changing environments, including the activation of plant defense response systems. Plants are continually exposed to adverse conditions in their environment whether they be under cultivation or as part of a natural system. Adverse growing conditions can lead to compromised plant growth, reproduction and productivity, and can be abiotic and biotic factors that may occur simultaneously. Abiotic factors are those such as drought and soil salinity while biotic factors include insect herbivory and disease caused by various pathogens. Employing new ways to activate plant-defense-response systems to counteract adverse factors could be transformative for agriculture and for enhancing biodiverse landscapes.
In this regard plant biostimulants, such as those made from seaweed extracts, are unique. Plant biostimulants are defined by a biological mode of action that utilizes plant mechanisms to provide their benefits such as enhanced tolerance to stresses, enhanced nutrient use and productivity (Brown and Saa, 2015). In Europe, the general principles used to justify plant biostimulant claims highlight that their effect is independent of nutrient content (Ricci et al., 2019). Plant biostimulants are used at low rates of application which differentiates their mode of action from synthetic nutritional fertilizers. The low dosage range is consistent with plant biostimulants having properties that accentuate plant response systems for better plant growth and improved tolerances. Many published studies have shown that biostimulants provide a multitude of plant growth benefits such as improved tolerances to abiotic and biotic stresses (Khan et al., 2009; Calvo et al., 2014). To achieve such a wide range of plant benefits across diverse plant families, implies that the molecular mechanisms underlying the plant responses are conserved, complex and pleotropic in character. Despite these insights their modes of action remain elusive.
Seaweed extracts are used successfully to improve agricultural productivity (Calvo et al., 2014; Arioli et al., 2015). A greater understanding of their biological modes of action will further enhance productivity gains in the future. There are a range of commercial seaweed-based products which are available off-the-shelf for commercial and home garden care and the majority of these claim that their use promotes plant growth, improves soil quality and/or enhances resistance against biotic and abiotic stress (Arioli et al., 2015). The effect of several of these products on abiotic and biotic stress mitigation and their mechanism of action has been explored in various studies (Khan et al., 2009; Shukla et al., 2019).
Our hypothesis was that plants treated with seaweed extract would increase their tolerances to subsequent stresses through the activation of a combination of plant defense responses. We envisaged this type of mode of action could be extended by combining different types of seaweed extracts.
The conditioning of plants to stress is an important feature for enhancing crop resilience and reducing productivity losses due to abiotic and biotic stresses (Kerchev et al., 2019). Plant conditioning is based on the molecular activation and priming of plant molecular defense systems so enhanced plant tolerance is exhibited upon subsequent stress occurrences (Martinez-Medina et al., 2016). Importantly, pre-treatment of plants using seaweed extracts is a practical approach to proactively initiate the conditioning phenomenon and was incorporated in our experimental design.
This study used a unique combination of approaches for new insights into the effect of different types of seaweed extracts on the activation of plant defense systems. Here, we used three different seaweed extracts, different plant response time points for assessments, a plant pre-treatment approach (to apply the different seaweed extracts) and the well-studied model system of Arabidopsis thaliana with the root pathogen Phytophthora cinnamomi. Three seaweed-derived extracts, namely extracts from the brown algae Durvillaea potatorum (native to the southern hemisphere) and Ascophyllum nodosum (native to the northern hemisphere), were used either separately or in a mixture, to treat Arabidopsis plants and then to compare activation of their plant defense systems. The model system of Arabidopsis thaliana with the root pathogen Phytophthora cinnamomi was utilized to trigger and synchronize the abiotic attack.
Phytophthora cinnamomi is an oomycete pathogen with an extremely wide host range. It is a notoriously aggressive forest pathogen and is considered a major threat to natural ecosystems in Australia (Cahill et al., 2008; Hardham and Blackman, 2018; Costa et al., 2020). P. cinnamomi is also a serious threat to horticultural, ornamental and nursery industries and, for example, causes one of the most damaging diseases of avocado. The pathogen infects the feeder roots and often the trunk of larger species causing disease that leads to branch die-back, loss of production and eventual death (Reeksting et al., 2016).
To characterize the plant defense response systems upon pathogen inoculation we used molecular and cellular techniques. High throughput RNA-sequencing was used at key plant response time points. This approach complimented the excellent transcriptomics reports on the action of biostimulant extracts (Nair et al., 2012; GonñI et al., 2016; Jithesh et al., 2019; Omidbakhshfard et al., 2020), particularly with respect to the ability of biostimulants to alter the outcome of root pathogen infection. To further confirm the defense transcriptome induction, Reactive Oxygen Species (ROS) were investigated by staining for hydrogen peroxide in the seaweed extract treated plants at 12 h post-inoculation (hpi) with P. cinnamomi. Microscopic analysis was performed to confirm that the pathogen had infected the inner root cell layers, and the extent of infection quantified using quantitative PCR. A plate assay was used to confirm that the seaweed extracts had no direct effect on pathogen growth.
We report that the extracts from two different brown seaweeds and their combination, activated plant defense responses upon pathogen-induced stress. Plants treated with each extract had different but overlapping transcriptomic gene expression profiles, and showed higher ROS levels that coincided with the activation of plant defense.
Materials and Methods
Arabidopsis thaliana Growth Conditions
Seeds of Arabidopsis thaliana ecotype Ler (LEHLE, TX, United States1) were surface-sterilized within a 1.5 mL microcentrifuge tube that contained 50% v/v ethanol (Chem-supply, Australia), 5% of H2O2 30% solution (Sigma-Aldrich, Australia) for 5 min and subsequently rinsed three times in sterile distilled water (sdH2O) and suspended in 0.2% (w/v) water agar. The seed suspension was stored in the dark at 4°C for 2−3 days. The stratified seeds were then seeded into Petri dishes (9-cm-diameter) containing Murashige and Skoog basal medium 0.44% (w/v) (Sigma-Aldrich, Australia) supplemented with 3% sucrose (w/v) (Chem-supply, Australia) and 0.8% (w/v) bacteriological agar and adjusted to pH 5.7 with 1 M potassium phosphate dibasic or potassium phosphate monobasic (Sigma-Aldrich, Australia). Seeds were evenly distributed in Petri dishes by placing individual seeds on the agar surface with a 1000 μL pipette tip, 120 seeds per plate. Petri dishes were transferred to a plant growth chamber (Thermoline Scientific, Australia) under cool white fluorescent lights (100 μmol photons m–1 s –1) with a 16:8 h (day: night cycle) at 21 ± 2°C for 14 days. Plants of uniform size were then selected for further use in experiments.
Plant Growth and Treatment With Seaweed Extracts
Plants were grown in a sand-based tube system that used commercial propagation sand (Bunnings, Australia) that was autoclaved and sterilized before adding to the tubes. The tubes used were 5 mL plastic disposable pipette tubes (Axygen, Australia) with a piece of cotton wool (Woolworths, Australia) inserted into the narrow end to form a plug that held the sand in place. Each tube was filled with sand to within 5 mm of the top and then 1 mL of diluted (1:400) seaweed extract or water as the control, was added at the top of the tube to just moisten the sand. The seaweeds used in this study are different, so the extracts are not identical, therefore we standardized the testing approach. The 1:400 dilution of the extracts were chosen because of (i) greenhouse and field studies demonstrating the efficacy for this dilution (Mattner et al., 2013, 2018), and (ii) by testing for root growth efficacy using the dilution as described previously (Arioli et al., 2015). For alignment with our previous greenhouse and field studies and the root growth testing, each of the seaweed extracts where standardized to 16% (w/w) soluble solids before preparing the 1 in 400 dilution for each seaweed extract. Three seaweed extracts designated as “AN” (an alkaline hydrolysis product from Ascophyllum nodosum), “DP” (an alkaline hydrolysis product from Durvillaea potatorum) and “AN/DP” (an alkaline hydrolysis product from both A. nodosum and D. potatorum, SeasolTM) were used in this study. Single plants of A. thaliana seedlings were gently removed from the MS plate and the roots carefully placed within a 10 mm deep hole made by pushing the narrow end of another 5 mL tube into the sand. Following placement of the seedling roots within the hole a further 1 mL of diluted seaweed extract (1:400) or water was added to enclose the root system by the sand. Tubes were then placed in a plastic rack and transferred to the plant growth chamber under the conditions described in section “Arabidopsis thaliana Growth Conditions.” Each day, and up until 6 days after placing the seedlings in the growth chamber, 700 μL of seaweed extract (1:400) was added to each tube or distilled water for the control.
Infection With Phytophthora cinnamomi Zoospores
Zoospores of P. cinnamomi were produced according to Islam et al. (2017) and the zoospore density adjusted to 1 × 105 zoospores/mL. Inoculation of the roots of plants grown in tubes took place on day seven whereby 700 μL of the zoospore suspension was carefully dispensed by pipette against the side wall of the plant growth tube just above the sand surface. The inoculated plants (8 plants/replicate/treatment) were then harvested at 0 h (i.e., immediately) and then at 3, 6, 12 and 24 h post-inoculation (hpi). To remove individual plants from a growth tube whilst avoiding damage to the root system a tube was briefly submerged in distilled water held within a container and the tube gently tapped to remove sand and the whole plant. The intact plant was then immediately placed with its roots submerged in water within a square plastic culture dish (10 × 10 cm) and the roots agitated gently to remove residual sand particles. Whole plants were gently and briefly dried on absorbent paper and frozen in liquid nitrogen followed by storage at −80°C. To confirm that roots had been inoculated, roots of eight plants from each treatment were sampled at 24 hpi and placed on PARPH medium (Islam et al., 2017) within 9cm-in-diameter Petri plates and examined for typical P. cinnamomi hyphal growth after 72 h incubation at 24°C in the dark. Images of whole root systems were also captured using a digital camera at 7 days after transferring the plants into the sand system and root length measured on individual plants with the aid of imageJ software. Final root growth data represent the mean of three biological replicates (each replicate contained 8 plants) from two independent repeats.
Gene Expression Analysis by Semi-Quantitative PCR
The plants were grown and treated with seaweed extracts as described in section “Arabidopsis thaliana Growth Conditions” and “Plant Growth and Treatment With Seaweed Extracts,” respectively. Then the plants were inoculated and harvested (8 plants for each time point for each treatment) at 0 h and then every 3 h until 9 hpi, as described in section “Infection With Phytophthora cinnamomi Zoospores.” Total RNA was isolated from plant tissues using a TRIzol®-based RNA extraction system (Life Technologies, United States) according to the manufacturer’s instructions. The quantity of RNA was measured using a Nanodrop® spectrophotometer (Thermo Fisher Scientific, United States) and the ratio (>1.8) of sample absorbance at 260/280 was used to determine the purity of samples. All samples extracted were of high yield and purity. The isolated RNA samples were then treated with DNAse-1 (Invitrogen, Carlsbad, CA, United States) according to the manufacturer’s instructions to remove any residual gDNA. Then the Tetro cDNA synthesis kit (Bioline, Australia) was used to synthesize cDNA from isolated RNA according to the manufacturer’s instructions.
Semi-Quantitative RT-PCR Conditions
The expression of genes involved in SA and JA mediated pathogen resistance pathways (PR1, PR5, NPR1, PDF1.2, and THI2.1) described by Lemarié et al. (2015) were examined using semi-quantitative PCR. The actin and EF-1 alpha genes were used as internal controls. The primer sequences of NPR1, PDF1.2, THI2.1, and Actin were as described by Eshraghi et al. (2011) and the primers for PR1, PR5, and EF-1 alpha were designed using primer3plus (Supplementary Table 1). PCR reactions were performed with GoTaq green master mix (Promega, United States) and each reaction contained 2 μL of cDNA and 0.5 μM of the respective primers. PCR cycles consisted of an initial denaturation step of 3 min at 95°C, followed by repetitions (28−36 cycles, depending on the primer set) of the following three steps: a 30 s denaturation step at 95°C, 30 s annealing step ranging between 54°C and 60°C, 1 min elongation step at 72°C and a final extension step at 72° for 5 min. Initial reactions were performed to determine the annealing temperature of each primer set and the appropriate cycle number of the PCR reaction. PCR products were analyzed on a 1% agarose gel with 0.5 × TBE buffer and visualized using gel red staining with a gel doc system. The final gel images are representative of two biological replicates from two experimental repeats.
Examination of the Host Transcriptome Using RNA-Seq
Plant Growth, Seaweed Extract Treatments, and Plant Inoculation
Arabidopsis thaliana plants were grown as described in section “Arabidopsis thaliana Growth Conditions” and seaweed extracts were applied as described in section “Plant Growth and Treatment With Seaweed Extracts.” The plants were inoculated as described in section “Infection With Phytophthora cinnamomi Zoospores” and harvested at 0, 3, 6, 12, and 24 hpi (8 plants at each time point for each treatment). Control plants, treated with water alone, were also harvested at these times as described in section “Plant Growth and Treatment With Seaweed Extracts.” The harvested plants were immediately frozen in liquid nitrogen and kept in a −80°C freezer until RNA was extracted. This experiment contained three biological replicates at each time point for each treatment (5 time points of four treatments including the control resulted in 60 samples for analysis).
RNA Extraction and cDNA Library Preparation
The total RNA was extracted from whole plants using a commercial kit (RNeasy Plant Mini Kit, Qiagen, Valencia, CA, United States) following the manufacturer’s instructions. RNA concentration and integrity was then determined first using spectrophotometry (NanoDrop ND-1000 spectrophotometer, Thermo Fisher Scientific, United States) with absorbance ratio of A260/280 nm and A260/230 nm. RNA integrity was also confirmed with the 4200 TapeStation system (Agilent Technologies, Santa Clara, CA, United States). Only RNA samples with an A260/280 nm ratio between 2.0 and 2.1 and RNA integrity number (RIN) greater than 7 (Sigurgeirsson et al., 2014) were used for further analysis. DNA libraries were constructed from total RNA of control and inoculated samples using the NEBNext® UltraTM II Directional RNA Library Prep Kit for Illumina® according to the manufacturer’s instructions (New England BioLabs, United States). The quality of the libraries was assessed by using a 4200 TapeStation 6000 system and their quantities were measured by Qubit dsDNA BR assay kit (Thermo Fisher Scientific, United States). The library was sequenced using the NovaSeq 6000 platform (Illumina) in a paired-end 150 bp run.
Processing of Sequenced Reads
The sequencer-generated raw reads were pre-processed and mapped to the reference genome using CLC Genomics Workbench (version 8.5.1, CLC Bio, Arhus, Denmark). During pre-processing of RNA-Seq data, adapter sequences, reads with >10% of unknown bases, low quality reads (sequences with more than 50% bases with quality value ≤ 5) and ambiguous bases were removed to obtain high quality reads for further analysis. High quality reads were mapped to the reference A. thaliana genome using the default parameters of the Workbench software to generate normalized gene expression values in the form of reads per kilobase of transcript per million mapped reads (RPKM). Proportion-based statistical analysis of differentially expressed reads was performed using default parameters to identify differentially expressed genes (DEGs) between water treated and each seaweed treated sample harvested at different hours post-infection by P. cinnamomi. DEGs were then filtered based on a FDR-corrected P-value of <0.05 and a fold change of ≥1.5 for up-regulated DEGs and ≤−1.5 for downregulated genes. The Illumina RNA-Seq datasets analyzed for this study have been deposited in the SRA database with the accession number of PRJNA6095902.
Functional Classification of Up-Regulated DEGs
Gene ontology and KEGG
The Blast2go 5 PRO (B2G) program was used to perform GO functional classification into biological process, molecular functions and cellular components to analyze the up-regulated gene function distribution at a macro level. Further, the B2G program was also used for Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations for up-regulated DEGs by searching against the KEGG database3.
Hierarchical clustering and heatmap visualization
For each treatment type, 30 genes with the largest sum of absolute t-test Z-scores were selected. With R version 3.6.2 (R Core Team, 2019), the heatmap.2 function was used to display Z-scores across the five time points. For multidimensional scaling analysis, the uniquely mapped read counts for each gene in each sample underwent library size normalization and distance estimation using the cmd scale function. Hierarchical clustering was performed using the hclust function with Spearman correlation values and a tree cut parameter of 0.67.
The Gaussian based t-test Z statistic (CLC Bio) was used to rank genes from most up-regulated to most down-regulated prior to Gene Set Enrichment Analysis (GSEA) with gene sets obtained from the Reactome and MapMan databases (Thimm et al., 2004; Naithani et al., 2016). Enrichment analysis was performed using the FGSEA R package version 1.12.0 (Korotkevich et al., 2019) with a significance threshold of FDR adjusted p-Values less than 0.05.
Quantitative PCR to Validate RNA-Seq Expression
The Tetro cDNA synthesis kit (Bioline) was used to synthesize cDNA from previously isolated RNA that was used for RNA-seq analysis. Briefly, 1 μg RNA was mixed with 1 μL random hexamer, 1 μL 10 mM dNTP mix (final concentration 0.5 mM), 1 μL ribosafe RNase inhibitor (final concentration 0.5 u/μL), 4 μL 5 × RT buffer, 1 μL tetro reverse transcriptase (final concentration 10 u/μL), and DEPC-treated water up to a total of 20 μL. Then the mix was incubated in the PCR machine according to the following order: initial incubation 25°C for 10 min followed by 45°C for 30 min and then the reaction terminated at 85°C for 5 min.
The primers of all tested genes were designed using primer3plus software (Supplementary Table 2) and annealing temperature of each primer pair was selected using gradient qPCR. The resulting qPCR product was analyzed via gel-electrophoresis to check that the correct gene product was obtained based on the primer design. Moreover, PCR efficiency of all genes was determined by a standard curve analysis of cDNA samples using a series of 10-fold dilutions of cDNA to determine the gene-specific PCR amplification efficiency for each primer pair used in RT-qPCR experiments. The real time PCR amplifications were carried out using SYBR Green detection chemistry. cDNAs were run in triplicate on 96 well reaction plates with the CFX Connect real time PCR (Bio-Rad, United States). 10 μL of reaction mixture containing 5 μL of iTaqTM universal SYBR Green Mix (Bio-Rad, United States), 0.4 μL of each 10 μM of primer and 2 μL of diluted cDNA as template and 2.2 μL RNase/DNase free sterile water (Sigma-Aldrich, Australia). The following amplification program was used in all RT-qPCR reactions: 95°C for 3 min, 40 cycles of 95°C for 10 s and annealing temperature (54−60°C) for 30 s at optimized temperatures for specific candidate genes. The specificity of each amplification reaction was verified by a melting curve analysis after 40 cycles. No template controls (NTC) were included for each primer pair to avoid possible contamination of assay reagents. Three biological replicates were used for each time point and each reaction was run in triplicate for each target and reference gene. All samples were run in parallel with actin reference genes (ACT2 and ACT8) to normalize cDNA loading. The relative expression values for each target gene were calculated against reference genes using the following equation according to Livak and Schmittgen (2001): ΔΔCT = (CT, Target – CT, reference gene) Time X − (CT, Target – CT, reference gene) Time 0.
Microscopic Examination of P. cinnamomi Infection of Roots of A. thaliana
To monitor the root infection process following inoculation with motile zoospores, the plants were grown with seaweed extracts or water as a control and inoculated with P. cinnamomi as described in section “Arabidopsis thaliana Growth Conditions,” “Plant Growth and Treatment With Seaweed Extracts,” and “Infection With Phytophthora cinnamomi Zoospores.” Then the plants were harvested at 12, 24, 48, and 72 hpi for each extract treatment and the roots removed and stained using a trypan blue staining protocol (Wang et al., 2011). Briefly, harvested roots were transferred into 2 mL Eppendorf tubes containing diluted trypan blue solution (10 g phenol, 10 mL glycerol, 10 mL lactic acid, 10 mL water and 10 mg of trypan blue). The tubes were incubated in a heated water bath and boiled for 4 min. After cooling to room temperature, the samples were de-stained by replacing the staining solution with chloral hydrate solution (5 g chloral hydrate/2 mL water) for 24 h. The samples were finally mounted in 50% glycerol and viewed with a light microscope (Zeiss, Göttingen, Germany) using bright field illumination. Images were captured with a digital camera (Zeiss, Göttingen, Germany) mounted on the microscope. The final images are representative of three biological replicates (each with at least 5 plants) at each time point for each treatment.
P. cinnamomi Quantification Using qPCR
To quantify the amount of P. cinnamomi in A. thaliana roots following inoculation plants were first grown with seaweed extracts (AN, DP, and AN/DP) for 6 days and inoculated with P. cinnamomi as described in section “Arabidopsis thaliana Growth Conditions,” “Plant Growth and Treatment With Seaweed Extracts,” and “Infection With Phytophthora cinnamomi Zoospores.” Then, the plant roots were harvested at 12, 24, 48, 72, and 96 h post-inoculation (hpi). The pathogen was quantified from harvested roots according to Engelbrecht et al. (2013). Briefly, DNA from a P. cinnamomi culture was extracted using PrepMan Ultra Reagent (Thermo Fisher Scientific) according to the manufacturer’s protocol and kept at −20°C until further use. DNA from root samples was extracted using the CTAB based method (Supplementary Method 1). The amount of plant genomic DNA present within the sample was quantified first by real-time PCR using primers (Supplementary Method 1) amplifying the A. thaliana actin gene. A normal one-step real-time PCR was conducted for the plant actin gene. The amount of plant DNA was calculated using a standard curve developed from a series of known concentrations of A. thaliana genomic DNA. The amount of P. cinnamomi DNA present within A. thaliana root samples was quantified using a nested real-time PCR method. LPV3 primers (Engelbrecht et al., 2013) were used in the outer first round of PCR, then LPV3N primers were used for the second round nested PCR to bind within the outer PCR product (Supplementary Method 1). The amount of pathogen DNA was calculated based on a standard curve developed from a series of known concentrations of P. cinnamomi DNA. Finally the quantity of P. cinnamomi was determined as ng of P. cinnamomi DNA/100 ng of A. thaliana DNA. The final amount determined represents the mean of three biological replicates from two independent experimental repeats.
Histochemical Localization of Hydrogen Peroxide (H2O2)
For hydrogen peroxide (H2O2) detection, the plants were grown and inoculated as described in section “Arabidopsis thaliana Growth Conditions,” “Plant Growth and Treatment With Seaweed Extracts,” and “Infection With Phytophthora cinnamomi Zoospores”. At 12 hpi, the whole plants were harvested and were placed in a 2 mL micro centrifuge tube with 1 mL diaminobenzidine (DAB) (1 mg/mL). The seedlings were incubated in the dark at room temperature for 3 h for H2O2 detection (Thordal-Christensen et al., 1997;Venus and Oelmüller, 2013). Samples were then transferred to and incubated in a decoloring solution (EtOH: lactic acid: glycerol = 1:1:1) at 80°C for 20 min (Venus and Oelmüller, 2013). Seedling roots were then visualized using a light microscope (Zeiss, Göttingen, Germany) with bright field illumination and images captured with a digital camera mounted on the microscope. The optical density of the colored precipitate was measured using imageJ software. The final images were representative of three biological replicates each with 5 plants.
The Influence of Seaweed Extracts on A. thaliana Root Growth and Infection by P. cinnamomi
Arabidopsis thaliana Ler seedlings were grown in a sand growth system (distilled water as a control and 1:400 seaweed extract as a treatment) and root lengths were measured 7 days after the commencement of treatment. All three seaweed extracts (AN, DP, and AN/DP) significantly enhanced root growth compared with the controls (Supplementary Figure 1). Moreover, root growth was observed to be significantly higher in extract-treated plants at the time of inoculation and the difference in root growth rate was maintained up to 96 hpi (Data not shown).
The quantitative measurement of the amount of P. cinnamomi in A. thaliana roots grown with seaweed extracts showed that overall there was less pathogen growth compared with the water control at 24, 48, and 72 hpi (Figure 1A). However, at 12 hpi, the amount of P. cinnamomi was higher in all three seaweed extract-treated roots compared to the water control. For the AN and DP seaweed extract- treated roots about the same amount of pathogen was found at 96 hpi, as was found in the water treated controls. Notably the combined extract (AN/DP) showed a plateauing of the amount of P. cinnamomi from 24 h onward to a level that was sustained well below that of the controls (Figure 1A).
Figure 1. Effect of seaweed extracts on P. cinnamomi infection in roots of A. thaliana. Plants were grown in a sand culture system with either seaweed extracts (AN, DP, and AN/DP) or water as a control for 6 days and then inoculated with P. cinnamomi on day 7. (A) Nested real time PCR quantification of P. cinnamomi DNA. Plant roots were harvested from 12 to 96 hpi. Data presented are from two experimental repeats each with three biological replicates. Error bars represent the standard error of means. *significant difference for amount of pathogen in different treatment compared to water control at P = 0.05 according to Duncan’s multiple range test. (B) Whole roots of A. thaliana infected with P. cinnamomi zoospores following treatment of roots with seaweed extracts (AN, DP, or AN/DP). Images were captured at 12 and 24 hpi. Scale bar = 20 μM. Each image is representative of three biological replicates.
Analysis of SA and JA Related Gene Expression
The expression of SA and JA-related genes, PR1, PR5, NPR1, PDF1.2 and THI2.1 was analyzed using semi-quantitative PCR and the Actin gene was used as an internal control to confirm even loading of DNA and reaction efficiencies of all cDNA samples prepared. A differential expression pattern was found for each resistance-related gene in all three seaweed extract treatments. Results showed that Actin expressed equally in all tested cDNA samples indicating the quality of cDNA and equal loading on the gel (Supplementary Figure 2). A higher expression of PR1 was found in AN-treated plants at 3 hpi and the expression increased at 6 and 9 hpi. A similar trend was found for the AN/DP treatment. However, consistently higher expression from 3 hpi was observed for the DP treatment. A similar expression pattern (up-regulated at 3 hpi) was recorded for PR5. Moreover, the expression of NPR1 was found to be consistently induced in all three treatments. The expression of PDF1.2 was found to be higher only in those plants treated with DP and AN/DP. In addition, the expression of THI2.1 was up-regulated at 3 hpi in plants treated with AN, and at 6 and 9 hpi in those plants treated with DP and AN/DP (Supplementary Figure 2).
Transcriptome Analysis of Plants Treated With Seaweed Extracts and Then Infected With P. cinnamomi
Confirmation of Pathogen Infection in Roots
Microscopic examination of A. thaliana roots grown with seaweed extracts and inoculated with P. cinnamomi revealed the different patterns of penetration and establishment of infection (Figure 1B). This microscopic analysis confirmed that the system established and optimized for this study was one in which the plants were successfully inoculated and that the pathogen grew both on the root surface and within the root.
Overview of RNA-Seq Data and Mapping to the A. thaliana Reference Genome
The Nova-Seq platform generated 43−155 million reads (average length 151 bp, paired end reads) and these reads were processed through CLC genomics workbench to remove adapters and ambiguous reads from the samples. A read refers to the sequence of a cluster that is obtained after the end of the sequencing process which is ultimately the sequence of a section of a unique fragment. After trimming, more than 98% of reads were recovered as high quality reads to proceed for mapping to the reference genome (Table 1).
On average, more than 87% of the total reads were mapped to the reference A. thaliana genome. The reads not mapped to the reference genome were expected to be pathogen reads as the samples were inoculated with P. cinnamomi. In addition, there were few reads that were mapped as broken pairs (Supplementary Table 3).
Overview of Differentially Expressed Genes (DEGs)
Differentially expressed genes (DEGs) were identified based on their expression value normalized through Reads Per Kilobase per Million mapped reads (RPKM) according to a previous study (Kasirajan et al., 2018). RPKM estimates the gene expression level of a gene normalized for both transcript length and library sequencing depth, allowing a direct comparison of expression levels within and between samples. The following parameters were considered to filter the DEGs of each treatment: FDR corrected P-value ≤ 0.05, fold change ≥ 1.5 for up-regulated genes, fold change ≤ 1.5 for down-regulated genes, each sample was compared with the respective water control. The highest number of up-regulated genes was found in the AN-treated samples harvested at 12 hpi followed by the AN/DP-treated samples harvested at 24 hpi. The highest number of down-regulated genes was found in the AN/DP-treated samples harvested at 3 hpi followed by those harvested at 6 hpi (Supplementary Table 4).
Functional Analysis of Up-Regulated DEGs
Comparison of number of DEGs between treatments and time points
Venn diagram analysis revealed that most of the DEGs from each of the three treatments were uniquely expressed at a specific time point (Figures 2A–C). For example, a total of 882 DEGs were found to be expressed in the AN treatment harvested at 12 hpi and among them 711 DEGs were uniquely expressed at this time point. Whereas there were 56 DEGs that were common to the 6 hpi time point, 37 DEGs were commonly expressed at the 24 hpi time point, 19 DEGs were commonly expressed at the 3 hpi time point and 13 DEGs at the 0 hpi time point.
Figure 2. Venn diagrams that show commonalities and differences among up-regulated DEGs at five time points following treatment of A. thaliana with (A) AN, (B) DP, (C) AN/DP, and (D) up-regulated DEGs (at at least one time point) for each of the three treatments.
Figure 2D shows that a total of 3,704 unique DEGs (expressed at least at one time point, filtered according to the criteria mentioned above) were found in the three treatments. Out of them, 926, 704 and 1027 were expressed in the AN, DP and AN/DP treatments, respectively. Moreover, 349 DEGs were commonly expressed in the AN & DP treatments whereas 277 DEGs were common to the AN & AN/DP treatments and 158 DEGs were common to the DP & AN/DP treatments. In addition, 263 DEGs were commonly up-regulated across all three treatments.
Gene ontology (GO)
The functional classification of up-regulated DEGs was analyzed using Gene Ontology (GO) and classified into three broad categories: molecular function, biological process and cellular component. These broad categories are very useful for identifying the key changes brought about by the treatments.
In the molecular function GO category, protein-binding and metal ion-binding were highly represented for all three treatments (Figures 3–5). Most importantly, the categories associated with plant defense pathways such as hydrolase activity, kinase activity, protein serine/threonine kinase activity, transcription factor binding, receptor serine/threonine kinase binding and terpene synthase activity were enriched in the analysis. For example, protein kinases play a central role in signaling in pathogen recognition and the subsequent activation of plant defense mechanisms (Romeis, 2001). Moreover, the genes identified as having hydrolase activity are likely to be involved in hydrolyzing the pathogen cell wall (Serrazina et al., 2015). The highest number of genes that were represented in the different molecular function categories was found at the 12 hpi time point for the AN treatment whereas for the other two treatments it was at the 24 hpi time point (Figures 3–5).
Figure 3. Gene ontology (GO) enrichment analysis of differentially expressed genes (DEGs) from the AN treatment. The DEGs were categorized into panels (A) Molecular function, (B) Biological process, and (C) Cellular component.
Figure 4. Gene ontology (GO) enrichment analysis of differentially expressed genes (DEGs) from the DP treatment. The DEGs were categorized into panels (A) Molecular function, (B) Biological process, and (C) Cellular component.
Figure 5. Gene ontology (GO) enrichment analysis of differentially expressed genes (DEGs) from the AN/DP treatment. The DEGs were categorized into panel (A) Molecular function, (B) Biological process, and (C) Cellular component.
Most of the up-regulated transcripts in all three seaweed extract-treated A. thaliana plants fell into the biological process categories of cell wall organization, oxidation-reduction process and phosphorylation (Figures 3–5). In terms of the most important categories related to plant defense pathway processes the following were identified: defense response, hydrogen peroxide catabolic process, response to salicylic acid, response to auxin, innate immune response, response to abscisic acid, auxin activated signaling pathway, ethylene-activated signaling pathway and response to jasmonic acid, all were enriched in the biological process category. Classical defense phytohormones such as salicylic acid (SA), jasmonic acid (JA), ethylene (ET) and more recently, growth-related phytohormones, such as auxins, cytokinins (CKs), brassinosteroids (BRs), abscisic acid (ABA), and gibberellins (GAs) have all been shown to modulate plant immune defenses (Han and Kahmann, 2019).
Interestingly, in a comparison between the three treatments, an up-regulation of the expression of genes related to systemic acquired resistance was found only in the DP treatment at different time points after inoculation (Figures 3–5). However, the induction of expression of three major SAR genes (PR1, PR5 and NPR1) was found in the qPCR validation of the RNA-Seq results (Figures 6–8).
Figure 6. Validation of the differentially expressed genes by RT-qPCR for A. thaliana plants treated with AN extract. Samples were collected from the plants grown with the seaweed extract and harvested at 0, 3, 6, 12, and 24 h after P. cinnamomi inoculation. All data were normalized to the expression level of actin 2 (ACT2) and actin 8 (ACT8). The data represent the fold change at each time point in the infected samples vs. the control sample. Bars show the standard error of the mean from three biological replicates.
Figure 7. Validation of the differentially expressed genes by RT-qPCR for A. thaliana plants treated with the DP extract. Samples were collected from the plants grown with the seaweed extract and harvested at 0, 3, 6, 12, and 24 h after P. cinnamomi inoculation. All data were normalized to the expression level of actin 2 (ACT2) and actin 8 (ACT8). The data represent the fold change at each time point in the infected samples vs. the control sample. Bars show the standard error of the mean from three biological replicates.
Figure 8. Validations of the differentially expressed genes by RT-qPCR for A. thaliana plants treated with the AN/DP extract. Samples were collected from the plants grown with the seaweed extract and harvested at 0, 3, 6, 12, and 24 h after P. cinnamomi inoculation. All data were normalized to the expression level of actin 2 (ACT2) and actin 8 (ACT8). The data represent the fold change at each time point in the infected samples vs. the control sample. Bars show the standard error of the mean from three biological replicates.
In the cellular component GO category, up-regulated DEGs in all three of the seaweed extract-treated and inoculated plants were principally assigned to the categories membrane, nucleus, integral component of membrane, plasma membrane and extracellular region (Figures 3–5). A similar predominance of these GO categories was also found in the resistance of plants to incompatible pathogens (Song et al., 2019). For example, cell surface receptors are trans-membrane proteins that bind signal molecules in the extracellular space and generate different intracellular signals on the opposite side of the plasma membrane (Alberts et al., 2014).
KEGG pathway analysis
Kyoto Encyclopedia of Genes and Genomes analysis showed that different biosynthetic and metabolic pathways were up-regulated in response to P. cinnamomi infection of A. thaliana grown with the seaweed extracts. The most highly represented top five pathways that contained the largest numbers of up-regulated genes were purine metabolism, biosynthesis of antibiotics, thiamine metabolism, starch and sucrose metabolism and phenylpropanoid biosynthesis. The up-regulation of gene expression related to antibiotic biosynthesis in all three treatments revealed that the plant may use “antibiotic compounds” to combat the pathogen. Moreover, several genes were up-regulated that are associated with phenylalanine metabolism and terpenoid backbone biosynthesis which lead to the synthesis of antimicrobial phytoalexins, phytoanticipins and phenolic compounds that are known to be involved in plant defense against pathogens (Cho and Lee, 2015; Supplementary Figure 3).
Clustering and heatmap visualization of DEGs
To visualize the expression pattern of DEGs we performed a hierarchical clustering of the DEGs that were extracted at each time point for each treatment and the respective control. The clustering heatmap (Supplementary Figure 4) showed a complex pattern of gene expression at each time point for each treatment compared to the water control. The heatmap showed that the expression pattern of the following groups were similar: DP-12 and DP-24, AN/DP-3 and AN/DP -6, AN/DP -0 and DP-6, AN/DP -12 and AN-12, H-6 and AN-6 (Supplementary Figure 4). Further, we visualized the top 30 DEGs across the time points for each treatment. The result showed that most of the genes were up-regulated at 12 hpi followed by 24 hpi for the AN treatment. Moreover, WRKY42 and CML8 showed the most consistent up-regulation across the time series (Supplementary Figure 5). However, for those genes that were up-regulated most were only up-regulated at 12 and 24 hpi for the DP treatment (Supplementary Figure 6). In addition, the expression patterns of selected genes were slightly different for the AN/DP treatment where the genes were found to be more highly up-regulated at each time point except at 3 hpi (Supplementary Figure 7).
Gene set enrichment analysis (GSEA) using Reactome and MapMan databases
The list of DEGs were mapped using GSEA to Reactome and MapMan databases to reveal any pathways that contained a large proportion of genes. The GSEA categorized the DEGs at each time point for each treatment into a number of functional groups. Among those groups obtained from the Reactome database for the AN treatment batch, auxin signaling was dominant at all time points. Also, ABA- and ET-associated genes were enriched at various time points. In addition, a number of important categories related to plant defense reactions such as SA signaling, recognition of fungal and bacterial pathogen and immunity responses were dominant at 12 hpi for the AN treatment (Supplementary Figure 8). Similarly, auxin signaling was dominant at 0 hpi as well as at the early infection stages (6 and 24 hpi) for the DP treatment. The steroid phytohormone, i.e., brassinosteroids-group was enriched at 3, 6, and 24 hpi for DP-treated plants. Most importantly, the SA signaling was dominant at 12 hpi for the DP-treated plants. In addition, the ethylene signaling group was found to be represented at only 6 and 24 hpi (Supplementary Figure 9). The ABA and ET signaling and biosynthesis groups were dominant at 0 hpi and at the early infection time of 3 hpi for the AN/DP treatment whereas brassinosteroids, SA and secondary metabolism groups were found at 3 and up to 12 hpi. The auxin signaling group was found to be dominant at 12 hpi in the AN/DP-treated plants (Supplementary Figure 10).
The different GSEA categories defined using the MapMan databases of DEGs from different treatments are presented in Supplementary Figures 11–13. The biotic stress_PR proteins_Plant defensins group was dominant at 0 and 3 hpi in the AN treatment whereas the biotic stress- associated group was found at later time points (6−24 hpi). The calcium signaling, protein degrading serine proteases, signaling G-proteins and signaling MAP kinases dominated at different time points. Most importantly, secondary metabolism of phenylpropanoids and peroxidases were highly represented at 12 hpi for the AN treatment (Supplementary Figure 11). Similar categories were also found for the DP treatment. However, protein degrading aspartate proteases, BHLH transcription factor and transcription regulator categories were found at different hpi for the DP treatments. A BZIP transcription factor category was highly dominant at 12 hpi for the DP treatment (Supplementary Figure 12). Many of these categories were also similarly found in AN/DP treatments at different hpi, except for 24 hpi where biotic stress was dominated by metabolite transporter, leucine-rich repeat signaling receptor kinases and MYB transcription factor family proteins (Supplementary Figure 13).
Understanding the expression pattern of important stress-related genes at different time points following infection with the pathogen is necessary for pinpointing their specific contribution to plant defense. A closer look at the MapMan profile in regards to biotic stress pathways affected by each seaweed extract treatment clearly showed that R genes, proteolysis, cell wall, beta glucanase, phytohormones, respiratory burst, heat shock proteins, secondary metabolites and transcription factor- associated genes were up-regulated at all time points for each treatment (Figure 9 and Supplementary Figure 14). In terms of phytohormones, Auxin, BRs, SA, ABA and ET-associated genes were represented and up-regulated at most of the time points for each treatment. However, JA associated genes were only up-regulated in the AN treatments at 6 hpi. In terms of respiratory burst, redox state- and peroxidases-associated genes were induced in all treatments. However, glutathione-S-transferase (GST) was only up-regulated at 12 hpi for both the AN and DP treatments and at 24 hpi for treatment with AN/DP. Importantly, the greatest number of genes in each category were found to be up-regulated at 12 and 24 hpi for each extract treatment (Figure 9 and Supplementary Figure 14).
Figure 9. Mapman overview of DEGs related to hormone, stress and metabolic responses in plants of A. thaliana following seaweed treatment and after inoculation with P. cinnamomi (at 3 and 12 hpi). The average fold change of genes are indicated by the color scale (red represents up-regulated genes and blue represents down-regulated genes).
Validation of RNA-Seq Expression
Based on their known involvement in plant defense pathways five genes from the DEG list for each treatment and three SAR-related genes were selected based on their initial gene expression analysis (see section “Analysis of SA and JA Related Gene Expression”) for RT-qPCR using their specific primers to confirm the reliability of expression of DEGs obtained from RNA sequencing. The relative expression levels of the selected genes were determined at 0, 3, 6, 12, and 24 h after P. cinnamomi inoculation. All of the selected genes from the three seaweed extract treatments showed trends closely aligned to those found for the RNA-seq data (Figures 6–8). In most cases, the RT-qPCR relative expression was higher than that found for the RNA-seq data for both up-regulated and down-regulated genes across the various time points examined. These results confirmed and further indicated that genes in A. thaliana related to a plant defense response, phytohormone signaling and transduction and systemic acquired resistance (Tables 2–4) were induced by the seaweed extracts and that they may function together against infection by P. cinnamomi.
Table 2. Possible function in defense of individual DEGs significantly induced at least at one time point for each seaweed extract treatment.
Key Plant Defense-Related Genes Significantly Up-Regulated in at Least One Time Point
The most important genes that are related to plant defense that were found to be significantly up-regulated in expression are shown in Tables 2–4. For the AN seaweed extract treatment the analysis of key defense-related genes revealed the presence of SA biosynthetic process and signaling associated genes, JA biosynthetic process and signaling-associated genes, pattern recognition receptors, plant defensin family gene and resistance-related gene active against oomycetes (Tables 2, 3). From an analysis of the DP seaweed extract treatment key defense-related genes revealed were those for SA and JA signaling, R protein encoded, auxin biosynthetic process-associated, phytoalexin production regulating, SAR regulating, basal resistance-related and hydrogen peroxide production associated genes (Tables 2,3). For the AN/DP seaweed extract treatment key defense-related genes found were for SA and JA biosynthetic process-associated genes, ABA-signaling genes, receptor-like protein kinases, an innate immune response inducer gene, SAR-inducible gene, a transcriptional regulator gene and ethylene signaling gene (Tables 2,3). In addition to those sets of genes specific to individual extract treatments there were those that were commonly found across the treatments. Noticeably the three key SAR-associated genes (PR1, PR5, and NPR1) were all up-regulated at 12 hpi following the AN and DP treatments showing a clear involvement of a SA-stimulated pathway. However, for the combined extract treatment at 12 hpi only NPR1 was up-regulated. In addition, the expression of auxin transporter, hydrogen peroxide responsive, receptor-like protein kinase and secondary metabolite biosynthetic genes were commonly found to be up-regulated across all three seaweed extract treatments (Table 4).
Table 3. Possible function in defense of individual DEGs significantly induced at least at one time point of at least at two treatments.
Table 4. Commonly up-regulated (at least at one time point) candidate genes following treatment with the three seaweed extracts.
Detection of Hydrogen Peroxide in A. thaliana Roots
The production of H2O2 was identified as a reddish-brown precipitate that resulted from DAB staining of the A. thaliana roots. At 12 hpi, H2O2 was detected in the roots grown with all three seaweed extracts and inoculated with the pathogen (Figures 10D,F,H). No H2O2 was detected in control roots grown with water or mock-inoculated with water (Figure 10A). Moreover, a minimal level of H2O2 was found in all three extract-treated and mock inoculated roots (Figures 10C,E,G). In addition, image analysis of DAB stained roots showed significantly higher stain in each seaweed extract-treated and inoculated root compared to either non-inoculated of each seaweed extract-treated root or water control (Supplementary Figure 15).
Figure 10. Hydrogen peroxide detection in A. thaliana roots grown with seaweed extracts or water as the control and inoculated with P. cinnamomi or mock-inoculated with water. Hydrogen peroxide was detected using the 3,3 O-diaminobenzidine tetrachloride (DAB) stain, which resulted in a reddish-brown precipitate in the root tissue. (A) Control root grown with water and mock inoculated with water showing no H2O2 production. (B) Control root grown with water and inoculated with the pathogen showing no H2O2 production. (C) Infected root grown with seaweed extract AN and mock inoculated with water showing a low level of H2O2 production. (D) Infected root grown with seaweed extract AN and then inoculated with the pathogen showing H2O2 production. (E) Control root grown with seaweed extract DP and mock inoculated with water showing low level H2O2 production. (F) Infected root grown with seaweed extract DP and inoculated with the pathogen showing H2O2 production. (G) Control root grown with seaweed extract AN/DP and mock inoculated with water showing low level of H2O2 production. (H) Infected root grown with seaweed extract AN/DP and inoculated with the pathogen showing H2O2 production. Scale bar = 20 μM. Each image is representative of three biological replicates.
Verification of the Arabidopsis thaliana-Phytophthora cinnamomi Plant-Pathogen System
This study used the model plant A. thaliana and the generalist, globally devastating pathogen P. cinnamomi, to examine the impacts of treatment of plants with two selected brown algal extract-based biostimulants or their combination, on pathogen growth and development in roots. Several previous studies have shown that applications of various brown algal extracts, either to soil or to foliage, enhanced root growth and plant development (Arioli et al., 2015; Mattner et al., 2018). Here we have used a sand culture system to grow A. thaliana with extracts from A. nodosum (“AN”), or D. potatorum (“DP”), or their combination (“AN/DP”). Studies on biostimulants and their impacts on plant disease establishment and progress have been reported (see for example, Gunupuru et al., 2019) although a comprehensive time course study of a root pathogen, and in this case an oomycete root pathogen, in the model plant A. thaliana has not been undertaken. The advantage of using A. thaliana, apart from its incredibly well detailed and characterized genome, is that there is a growing body of information around the interaction of this host with a range of oomycete pathogens including P. cinnamomi (Robinson and Cahill, 2003; Rookes et al., 2008), Phytophthora porri (Roetschi et al., 2001) and P. parasitica (Le Berre et al., 2017), Hyaloperonospora arabidopsidis (Kunz et al., 2008; Ried et al., 2019), and Albugo candida (Cooper et al., 2008), but none on their interactions with biostimulants.
In other host−pathogen systems treatment with seaweed extract-based biostimulants have indicated that disease incidence and severity may be reduced following infection. For example, commercial seaweed extracts from A. nodosum and D. potatorum were found to suppress disease caused by Plasmodiophora brassicae in broccoli (Wite et al., 2015) and an extract from A. nodosum reduced the severity of Fusarium head blight caused by F. graminearum in wheat (Gunupuru et al., 2019). It is worth noting the diversity in seaweed extracts. Liquid seaweed extracts are processed from seaweed biomass using different chemical approaches (such as acid and alkaline extraction) and cellular disruption under pressure (Arioli et al., 2015). The extracts comprise diverse molecules that are heterogenous in nature and representative of the extraction process, which emphasizes the need to characterize their properties. The results of our study show that the extent of colonization by P. cinnamomi of roots of A. thaliana was suppressed by pre-treatment of roots with the alkaline-based extracts from both A. nodosum and D. potatorum and a mixture of both. The generalist pathogen, P. cinnamomi, is an aggressive pathogen that is able to infect close to 5000 plant species (Hardham and Blackman, 2018). Therefore, the suppression of this pathogen by seaweed extracts is a significant finding that demanded further investigation of the details of the potential resistance mechanisms stimulated by different seaweed extracts against P. cinnamomi infection.
The availability of genetic and genomic tools for the model plant A. thaliana makes it a very good system in which to investigate the in planta action of seaweed extracts. Infections by Phytophthora spp. in A. thaliana have not been found under natural conditions, but have been achieved for several Phytophthora species under laboratory conditions (Herlihy et al., 2019). Ecotypic variation to infection by P. cinnamomi was described in an earlier study where ecotype Ler was found to be moderately susceptible (Robinson and Cahill, 2003). The microscopic analysis performed in the present study has confirmed that the system that was established and optimized was one in which the plants were successfully inoculated and that the pathogen grew both on the root surface and within roots.
Effect of Seaweed Extracts on Key Regulatory Resistance-Related Genes
Various seaweeds are a rich source of unique bioactive compounds, for example fucans, carrageenans, ulvans, and laminarins that have been shown to induce plant defense against a variety of pathogens (Cluzet et al., 2004; Shukla et al., 2016, 2019). These elicitor-like molecules may act as priming molecules or pathogen-associated molecular patterns (PAMPs) and thereby activate induced systemic resistance (ISR) and SAR responses. To first test this hypothesis in our system, we examined the expression of three SA- and two JA/ET-responsive marker genes (PR1, NPR1, PR5 and PDF1.2, THI2.1, respectively) that are related to SAR (Eshraghi et al., 2011). In our study, each seaweed extract was found to enhance the expression of the key SA-marker genes from the earliest time point tested after inoculation following treatment with seaweed extracts. The JA-marker genes showed variation in expression depending on extract type and time after inoculation. Moreover, the genes were not induced when the plants were treated with the seaweed extracts alone. Other studies have shown similar upregulation of these genes at a single time point. For example, it had been shown earlier that the expression of the PR1 gene in A. thaliana was up-regulated at 24 h post-treatment with an A. nodosum-based extract (Cook et al., 2018; Shukla et al., 2019). Furthermore, carrot plants primed by A. nodosum-derived extracts induced the accumulation of transcripts of the same or similar genes (Jayaraj et al., 2008). In contrast to these limited studies, the current study has identified the induction of key regulatory genes across a range of time points after pathogen infection therefore providing a post-infection, spatio-temporal analysis of induction following various seaweed extract treatments.
Transcriptome Analysis Revealed the Complexity of Resistance Induced by Seaweed Extracts
Summary of Transcriptional Changes Induced by Each Seaweed Extract With or Without Root Infection by P. cinnamomi
Transcriptome analysis using RNA-seq was performed to explore the whole plant transcriptome to reveal correlations between seaweed extract treatment and pathogen suppression that was found following quantitative analysis of the amount of pathogen within roots. Overall the results show that there was a large number of genes that were either up-regulated or down-regulated following exposure of plant roots to each of the extracts. In this study we have specifically concentrated on genes that were up-regulated in these interactions. Three major SAR-related genes were found to be up-regulated in common between extracts and were confirmed by both RNA-Seq and qRT-PCR validation. Equally importantly, each seaweed extract was found to exert its effect through different subsets of genes. A number of studies explored plant transcriptomes following abiotic and biotic stress (see for example: Tommasini et al., 2008; Allardyce et al., 2013; GonñI et al., 2016; Islam et al., 2018; Jithesh et al., 2019), however, the present study is the first report of a comprehensive transcriptome analysis, following treatment of plants with seaweed extracts, root pathogen infection and analysis over multiple time points.
The up-regulated genes were broadly identified as being involved in phytohormone signaling, defense responses, hydrolase activity and the biosynthesis of antibiotics, and also transcription factors and transcription regulators that were involved in metabolite biosynthesis. The GSEA using both Reactome and MapMan databases indicated that the DEGs were involved in a diverse range of activities during seaweed extract-induced plant defenses. For example, brassinosteroid (BR) signaling was commonly found across all three treatments. BRs are plant steroidal hormones that play vital roles in not only plant growth and development but also in plant defense through coordination with other phytohormones (Saini et al., 2015). Another example, proteases (serine or aspartic proteases), were commonly enriched in both AN and DP treatments. Plant genomes encode a large number of proteases which play a regulatory role in a number of processes that are essential for immune responses, more specifically, programmed cell death (PCD) (Balakireva and Zamyatnin, 2018). Most importantly, an array of proteolysis-related genes and their increased expression was commonly found at all time points for each seaweed extract treatment. Proteolysis machinery acts mainly in a housekeeping role to remove non-functional proteins, however, proteolysis has also been shown to play a key role in the recognition of pathogens and the subsequently induced effective defense responses (van der Hoorn and Jones, 2004). Therefore, the results of our study indicated the deployment of multiple phytohormones and proteolytic machinery in seaweed extract-induced defense against P. cinnamomi. However, there is considerable scope to further investigate the role of individual proteases in seaweed extract-induced defenses.
The plant cell wall is a dynamic and highly controlled structure that is essential for growth and development. It is considered to be a passive defense barrier against a variety of attackers. Plants have mechanisms that maintain cell wall integrity which comprise a set of so-called “plasma membrane-resident sensors” and “pattern recognition receptors” (Bacete et al., 2018; De Lorenzo et al., 2019). When a pathogen alters the cell wall integrity during epidermal penetration or through deeper colonization of sub-epidermal cells, plants activate suites of genes for cell wall biosynthesis and remodeling as repair and defense responses. This activation of genes and the production of their downstream products is very effective at stopping, or slowing down, pathogen ingress. Several studies, including those that have used overexpressor mutants, have demonstrated the central importance of cell wall-related genes in enabling increased disease resistance (Miedes et al., 2014; Bacete et al., 2018). In our study, cell-wall associated genes were dominant at all-time points following infection for each seaweed extract treatment.
Cell walls are the first line of defense and their modification a very early response to pathogen attack. Each seaweed extract stimulated cell wall-related gene activity following pathogen attack that was well above that for the water-treated control. The induction of these genes at early stages of infection for each seaweed extract treatment was strongly indicative of their contribution toward strengthening the cell wall against pathogen penetration. The induction of similar genes at later stages suggested their contribution to cell wall repair and the fortification of new cell walls. For example, the MYB46 transcription factor that was up-regulated at 12 hpi in the AN treatment, is directly involved in regulation of the expression of genes responsible for secondary cell wall formation including lignin and cellulose biosynthesis (Miedes et al., 2014). Another example, CALS5 (Callose synthase 5) that was up-regulated at 6 hpi in the AN/DP treatment, is involved in callose synthesis and was also a pathogen-induced gene in A. thaliana infected with the downy mildew pathogen Hyaloperonospora arabidopsis (Dong et al., 2008). Callose is a well-known plant defense component and is considered an effective barrier against pathogen invasion including in various A. thaliana ecotypes infected by P. cinnamomi (Robinson and Cahill, 2003).
The baseline of plant defense is the activation of PRRs localized in the plasma membrane upon recognition of PAMPs/MAMPs (Bigeard et al., 2015). Indeed the induction of RLK1 in AD-treated plants indicated the activation of PAMP-triggered immunity in these plants against the pathogen. The plant hormones SA, JA and ET have a significant role in plant defense against pathogens. The SA signaling pathway that activates programmed cell death is effective against biotrophic pathogens whereas JA and ET signaling pathways are effective against necrotrophic pathogens (Glazebrook, 2005). The upregulation of both JA and SA biosynthetic or signaling genes suggested the activation by seaweed extracts of both pathways in response to P. cinnamomi. The phytohormone, auxin, is well known to be a regulator of plant growth and development. However, auxin is also being recognized as a key regulator of plant defense (Wang and Fu, 2011; Wang et al., 2011). In our study, for example, the upregulation of TAA1, a gene involved in auxin biosynthesis, in DP-treated plants indicated the involvement of auxin signaling pathways in response to the pathogen. Similarly, ABA is mainly involved in abiotic stress tolerance as well as in biotic stress but it also may promote plant defense in a complicated network of synergistic and antagonistic interactions (Ton et al., 2009). The induction of an ABA biosynthesis-related transcription factor gene (WRKY40) and ABA responsive gene (ABR1) in plants treated with AN/DP, along with other key phytohormone-related genes suggested the activation of multiple phytohormone signaling pathways following seaweed extract treatment.
The second layer of plant defense is based on plant disease resistance, (R) gene, mediated resistance by recognition of the products of pathogen avirulence genes and subsequent effector-triggered immunity (ETI) (Andersen et al., 2018). For example, in the current study the induction of CRK5, which likely functions as a receptor-like kinase (Chen et al., 2003), in DP-treated plants indicated that ETI may have been triggered. WRKY transcription factors are encoded by a large gene superfamily with a broad range of roles in plants and several groups have reported that proteins containing a short VQ motif interact with WRKY motifs. One of these, VQ25, was reported by Cheng et al. (2012), to be involved in resistance to the necrotrophic pathogen Botrytis cinerea. The induction of the candidate resistance-related gene VQ25 in plants treated with AN thus indicated a contribution of this gene to AN-induced plant defense.
The recent review published by Shukla et al. (2019) presented additional information on some of the plant defense components activated by different extracts from A. nodosum. The bioactive compounds present in the prepared A. nodosum alkali extract (ANE) were proposed to elicit defense responses to pathogens. The application of ANE enhanced the activation of various enzymes including peroxidases and phenylalanine ammonia-lyase. In addition, ANE also induced ISR against P. capsici, another oomycete pathogen, that caused disease in tomato. Further, ANE induced SA-related genes and several JA-related genes such as PDF1.2 and plant immune response genes such as WRKY30 and CYP71A12, that we also highlight in our study. The review emphasized the information gap around the role of phytohormones in activating defense-related genes that we have now gone some way to fill. For example, the up-regulation of candidate genes, such as PCC1, ACD6, GR1, ERF014, AOC3, ACS9, and ACS11 all hormone-related by the different extracts derived from both A. nodosum and D. potatorum and their combination used in our study.
The array of plant defense responses that are activated during pathogen invasion requires an abundant supply of energy which is predominantly derived from primary metabolic processes. These primary metabolic pathways are used by plants not only as a source of energy to drive diverse defense responses, but also as a source of signaling molecules to directly or indirectly, trigger defense responses (Rojas et al., 2014). In the current study primary metabolic pathway activation following pathogen infection was a key outcome of seaweed extract treatment and presumably acted as an energy provider and regulator of Arabidopsis defense responses. For example, purine metabolites provide an ongoing source of nitrogen for A. thaliana growth. One of the purine metabolites, allantoin, plays a role in a JA-signaling pathway, suggesting that the role of purine metabolism not only underpins normal plant growth but, as others have found, is also a player in stress hormone homeostasis and signaling (Takagi et al., 2016). In our study, a large number of DEGs from each treatment were classified into purine metabolism through KEGG analysis. Therefore, this result indicated that purine metabolites acted to maintain plant growth during pathogen infection as well as contributing to defense-related hormone signaling pathways. In addition, at 12hpi the highest number of up-regulated purine metabolism genes was found for both AN and DP treatments whereas it was only at 24 hpi for the AN/DP treatment. This difference may have indicated a more sustained defense activation and supply of energy in the AN/DP-treated plants during pathogen infection.
Thiamine metabolism has an important function in many metabolic reactions including in glycolysis, the pentose phosphate pathway and the tricarboxylic acid cycle. In addition, thiamine is also related to the induction of SAR and is involved in plant adaptation toward biotic and abiotic stresses (Kamarudin et al., 2017). For example, several studies reported that thiamine treatment of plants, including A. thaliana, activated plant defense and enhanced resistance to disease (Ahn et al., 2007; Boubakri et al., 2012; Kamarudin et al., 2017). Therefore, up-regulation of thiamine metabolism which was demonstrated in our study has strong implications for its involvement in the induction of defense, as well as adaptation, during infection by P. cinnamomi.
Common Plant Defense-Related DEGs That Were Up-Regulated Following Inoculation With P. cinnamomi in Extract-Treated Plants
Two hundred and sixty three genes (1.3% of the genome) were commonly found to be up-regulated for at least one time point following inoculation with P. cinnamomi of extract-treated plants. A number of candidate resistance-related genes were found to be up-regulated across all treatments including PIN2, GRI, UGT73B4, CRK15, and MLO-8 which have been implicated in diverse resistance-related roles in different host and pathogen combinations. In addition to these genes, even though not above the cutoff by our RNASeq analysis, the pathogenesis-related genes PR1, NPR1, and PR5 were confirmed to be commonly up-regulated following treatment with the extracts through our preliminary semi-quantitative PCR analysis as well as in the RNA-seq validation that used quantitative PCR. The PR proteins are a group of proteins that are induced by phytopathogens through activation of specific defense-signaling pathways and are fundamental components of resistance regulation (Backer et al., 2019). After pathogen infection, activation of defense-signaling pathways, such as those regulated by SA and JA take place which further leads to the accumulation of PR proteins that stops pathogen growth and development within host tissues. The SA pathway is especially active following infection by biotrophic pathogens and which stimulates the transcription of NPR1 which in turn leads to activation, as well as accumulation, of SA-induced PR signature gene (PR1, PR2, and PR5) products locally and systemically that leads to SAR (Ali et al., 2018; Backer et al., 2019). PIN proteins are responsible for polar localization in the plasma membrane that determines the direction and rate of intercellular auxin flow (Sun et al., 2011). Moreover, GR1 plays a crucial role in coordinating gene expression through both SA- and JA-signaling pathways (Mhamdi et al., 2010). The induction of these genes and other phytohormone-related genes in our study suggested that all three extracts induced defense against P. cinnamomi that was dependent on the activation of multiple phytohormone signaling pathways. In addition to all the above defense interactors, receptor-like kinases such as CRK15 found to be up-regulated across treatments in our study, are fundamental signaling components that regulate a variety of cellular processes (Lee et al., 2017).
Plant secondary metabolites have numerous functions in plant−pathogen interactions and experimental evidence has demonstrated their important contributions in plant innate immunity (Piasecka et al., 2015). Plant-produced antibiotics are antimicrobial secondary metabolites and can be broadly classified as phytoalexins and phytoanticipins (Morrissey and Osbourn, 1999). UGT for example, plays an essential role in the biosynthesis of secondary metabolites in plants (Guo et al., 2016) and the induction of UGT73B4 and other genes associated with the biosynthesis of antibiotics in all seaweed extract treatments indicated the synthesis of potentially novel antimicrobial compounds as a reaction to infection by P. cinnamomi.
The diverse patterns of differential gene expression found in our study were consistent with seaweed extracts having complex and pleiotropic modes of action that involved a cascade of gene activation for different plant responses. The commonality in the transcriptome profiles suggested that, at least for the seaweed extracts derived from the brown seaweeds used in the current study, behaved in a similar, but not identical, way.
Novel Genes That Were Up-Regulated That Provide Insight Into the Mechanisms of Action of Seaweed Extracts Against P. cinnamomi
WRKY transcription factors play important roles in plant responses to various biotic and abiotic stresses. WRKYs act as substrates of calcium-dependent protein kinases and calmodulin (CaM) is a Ca2+ -binding protein that is involved in various cellular functions (Bakshi and Oelmüller, 2014). The function of calmodulin-like (CML) proteins is largely unknown. However, one of these, CML8 was found to be up-regulated in our study and has been shown to be involved in Arabidopsis immunity against Pseudomonas syringae (Zhu et al., 2017). The strong and consistent upregulation of WRKY42 and CML8 in the AN treatment indicated a correlation with, and the increased involvement of, calcium signaling in defense activation.
The A. thaliana genome has four jasmonate-induced oxygenases (JOXs) and one of them hydroxylates jasmonic acid to 12-OH-JA (Caarls et al., 2017). In our study the higher expression of JOX1 in the DP treatment at 6 hpi indicated the involvement of other phytohormone signaling pathways at this stage of the interaction with the pathogen. The expression of the POLARIS (PLS) gene that encodes a 36-amino acid peptide that regulates plant root growth and vascular development, is related to auxin transport and coordinates the ethylene signaling pathway (Chilley et al., 2006). The strong up-regulation of PLS expression in AN/DP treatments suggested that PLS also contributed to A. thaliana root growth, as well as functioning in regulation of phytohormone-induced signaling pathways, that resulted in suppression of P. cinnamomi. In addition, the strong up-regulation of other uncharacterized genes, in all three treatments, suggested the contribution of unknown novel mechanisms in AN and DP extract-induced defense.
Other publications have compared the transcriptional profiles of plants treated with seaweed extracts derived from the brown seaweed A. nodosum (Nair et al., 2012; GonñI et al., 2016; Santaniello et al., 2017; Jithesh et al., 2019; Omidbakhshfard et al., 2020). The extracts used in these studies varied in their chemical nature (including alkaline, neutral and acid extracts) and extraction approaches. Despite the differences observed among the transcriptional profiles following extract-treatment of plants, the overall results demonstrate the highly dynamic and responsive nature of plants to different types of seaweed extracts, and the inherent capacity for the seaweed extracts to simultaneously enhance plant growth and tolerance to biotic and abiotic stresses.
Role of ROS in Seaweed Induced Plant Defense
To examine the production of defense-related components prior to and following the upregulation of defense-related transcripts, the production of hydrogen peroxide (H2O2), a reactive oxygen species (ROS), was examined in A. thaliana roots grown with each seaweed extract and infected, or not, with P. cinnamomi. The induction of hydrogen peroxide was found only in those roots treated with the three seaweed extracts individually and infected with P. cinnamomi. This result is somewhat different to that of a previous study (Cook et al., 2018) that showed the induction of reactive oxygen species in seedlings treated only with an A. nodosum extract. This variation in results between the two studies may reflect differences in preparation of the extract and the treatment and analysis methods. The method used in our study gave a direct visualization of the location and intensity of ROS in the roots, something that is not possible using alternative assays.
The PEP1 gene, which was found to be up-regulated in our study in DP-treated plants, is involved in activation of the synthesis of enzymes associated with hydrogen peroxide formation (Huffaker et al., 2006). Also, plant peroxidases participate in various physiological processes, such as lignification, suberisation, auxin catabolism and defense mechanisms that are activated during pathogen infection. They are considered to catalyze the generation of aromatic oxyl radicals from several aromatic compounds and the peroxidase-dependent production of such organic radicals often results in the generation of reactive oxygen species (Kagan et al., 1990; Kawano, 2003). The induction of peroxidases in all extract treatments suggested their involvement in ROS generation and potentially other aspects of extract-induced defense mechanisms.
The Trade-Off Between Growth and Defense and Priming for Defense
New insights into how plants balance growth while responding to stress has implications for advanced agriculture. The compromise between growth and stress response is based on plants having limited resources which need to be prioritized for growth, or toward responses to the abiotic and/or biotic stresses they encounter (Huot et al., 2014; Karasov et al., 2017). The trade-off concept is supported by research that has demonstrated that plant-fitness costs are associated with the induction of defense genes (Huot et al., 2014; Karasov et al., 2017), and that maintaining activated plant response systems is metabolically costly (Karasov et al., 2017; Buswell et al., 2018). Nevertheless, emerging research has uncovered chemical priming of immunity that provides defense without costs to plant growth (Buswell et al., 2018).
In our research we observed a balanced trade-off between root growth and the activation of specific defense pathways. Our transcriptomic analysis identified the up-regulation of specific defense-associated pathways (such as those regulated by SA, JA, ET, and PPP) and genes associated with plant resistance (for example PR1, MLO, and others), while root length growth continued, despite the interior of the roots being actively infected by P. cinnamomi. In our experimental design pretreating the plants with seaweed extract was an important pre-requisite. This approach may have contributed to a favorable trade-off that utilized a natural plant priming system. More generally, extracts of different seaweeds have been shown to activate broad spectrum defense systems in plants (Kerchev et al., 2019; Shukla et al., 2019). Collectively the research supports the notion that seaweed extracts may act as a plant priming stimulant, particularly if pre-applied.
Plant-priming is an adaptive and low-cost defensive mechanism that, upon activation by a priming stimulus, results in a faster and/or stronger induction of inducible defenses. Plant-priming occurs in a wide range of plant species and is often associated with enhanced abiotic and biotic stress tolerance (Martinez-Medina et al., 2016). The idea that priming is the result of treatment with a specific seaweed extract is supported by our transcriptomics analysis particularly based on the molecular and cellular GO categories: for example, genes up-regulated for (i) DNA Binding Transcription Factor Activity (ii) Transcription Regulatory Region DNA Binding, (iii) the increased number of transcripts found in the nucleus, (iv) and up regulation of genes associated with redox signaling and sensing. Also ROS are key molecules involved in the priming process (Borges et al., 2014). Our data showing the accumulation of hydrogen peroxide in root cells at 12 hpi, in plants pretreated with seaweed extract, was consistent with ROS acting as a latent signal involved in priming plant resistance (Gonzalez et al., 2009). The priming of plant resistance can also be achieved by exogenous application of synthetic and natural compounds (Aranega-Bou et al., 2014). Hexanoic acid, for example, is a natural primer (Aranega-Bou et al., 2014) which is a component of one of the seaweed extracts used in this study (AN/DP, data not published).
This study showed that A. thaliana was a very useful model plant for studies on the impact that a seaweed extract-based biostimulant had on interactions at the molecular level with a root pathogen. We also demonstrated the up-regulation of key SAR-related genes and phytohormone- associated genes at various critical time points post-inoculation following treatment with extracts of the selected brown seaweeds. Importantly, each seaweed extract induced multiple defense-related pathways prior to penetration and infection by the pathogen. These observations were characteristic of a primed response, and closely associated with ROS production. Transcriptomic analysis has proven to be a powerful approach to elucidate the timing of activation of defense-related mechanisms and the subsequent suppression of pathogen growth. Our results can now be used in future studies that use specific plant mutants that are impaired in various resistance-related pathways or, for example, gene edited hosts to investigate the role of individual defense components in seaweed extract-induced defense. Further, we propose that the approach used in the current study could be applied to agriculturally important crop species to investigate the impact of seaweed extract-treatment on their reaction to a pathogen.
Data Availability Statement
The datasets presented in this study have been deposited in online repositories. The names of the repositories and their accession number(s) can be found either in this article or in the/Supplementary Material.
MI, DC, and TA conceptualized and designed the project. MI performed the laboratory work, analyzed the RNA-seq transcriptome data, and wrote and revised the manuscript. HG contributed to RNA-seq data processing. MZ performed the functional analysis of DEGs. HH contributed to the critical analysis of the RNA-seq data and manuscript draft. DC and TA contributed to the manuscript draft and final version. All authors revised and approved the final version to be published and agreed with all aspects of the work.
This project was funded through the Faculty of Science, Engineering and Built Environment (SEBE), Deakin University, Australia post-doctoral support for industry projects scheme and jointly with Seasol International Pty Ltd., Australia.
Conflict of Interest
Seasol® was manufactured by the company Seasol International (SI). TA and HH are researchers in the SI R&D Department and employed by SI. SI provided in-kind support as salary for TA and HH.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationship that could be construed as a potential conflict of interest.
We would like to thank Drs. Pavani Nadiminti and Xinhua Lu, Deakin University, for assistance with growing plants and other aspects of the laboratory experiments for the project.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2020.00852/full#supplementary-material
Ahn, I. P., Kim, S., Lee, Y. H., and Suh, S. C. (2007). Vitamin B1-induced priming is dependent on hydrogen peroxide and the NPR1 gene in Arabidopsis. Plant Physiol. 143, 838–848. doi: 10.1104/pp.106.092627
Alberts, B., Johnson, A., Lewis, J., Raff, M., Roberts, K., and Walter, P. (eds). (2014). “General principles of cell communication,” in Molecular Biology of the Cell, 4th Edn, (New York, NY: Garland Science).
Ali, S., Ganai, B. A., Kamili, A. N., Bhat, A. A., Mir, Z. A., Bhat, J. A., et al. (2018). Pathogenesis-related proteins and peptides as promising tools for engineering plants with multiple stress tolerance. Microbiol. Res. 212, 29–37. doi: 10.1016/j.micres.2018.04.008
Allardyce, J. A., Rookes, J. E., Hussain, H. I., and Cahill, D. M. (2013). Transcriptional profiling of Zea mays roots reveals roles for jasmonic acid and terpenoids in resistance against Phytophthora cinnamomi. Funct. Integr. Genomic. 13, 217–228. doi: 10.1007/s10142-013-0314-7
Aranega-Bou, P., De La, O., Leyva, M., Finiti, I., García-Agustín, P., and González-Bosch, C. (2014). Priming of plant resistance by natural compounds. Hexanoic acid as a model. Front. Plant Sci. 5:488. doi: 10.3389/fpls.2014.00488
Arioli, T., Mattner, S. W., and Winberg, P. C. (2015). Applications of seaweed extracts in Australian agriculture: past, present and future. J. Appl. Phycol. 27, 2007–2015. doi: 10.1007/s10811-015-0574-9
Backer, R., Naidoo, S., and Van Den Berg, N. (2019). The nonexpressor of pathogenesis-related genes 1 (NPR1) and related Family: mechanistic insights in plant disease resistance. Front. Plant Sci. 10:102. doi: 10.3389/fpls.2019.00102
Berkey, R., Zhang, Y., Ma, X., King, H., Zhang, Q., Wang, W., et al. (2017). Homologues of the RPW8 resistance protein are localized to the extrahaustorial membrane that is likely synthesized de novo. Plant Physiol. 173, 600–613. doi: 10.1104/pp.16.01539
Bi, D., Cheng, Y. T., Li, X., and Zhang, Y. (2010). Activation of plant immune responses by a gain-of-function mutation in an atypical receptor-like kinase. Plant Physiol. 153, 1771–1779. doi: 10.1104/pp.110.158501
Borevitz, J. O., Xia, Y., Blount, J., Dixon, R. A., and Lamb, C. (2000). Activation tagging identifies a conserved MYB regulator of phenylpropanoid biosynthesis. Plant Cell 12, 2383–2393. doi: 10.1105/tpc.12.12.2383
Borges, A. A., Jiménez-Arias, D., Expósito-Rodríguez, M., Sandalio, L. M., and Pérez, J. A. (2014). Priming crops against biotic and abiotic stresses: MSB as a tool for studying mechanisms. Front. Plant Sci. 5:642. doi: 10.3389/fpls.2014.00642
Boubakri, H., Wahab, M. A., Chong, J., Bertsch, C., Mliki, A., and Soustre-Gacougnolle, I. (2012). Thiamine induced resistance to Plasmopara viticola in grapevine and elicited host–defense responses, including HR like-cell death. Plant Physiol. Biochem. 57, 120–133. doi: 10.1016/j.plaphy.2012.05.016
Breitenbach, H. H., Wenig, M., Wittek, F., Jordá, L., Maldonado-Alconada, A. M., Sarioglu, H., et al. (2014). Contrasting roles of the apoplastic aspartyl protease apoplastic, enhanced disease susceptibility1-dependent1 and legume lectin-like protein1 in Arabidopsis systemic acquired resistance. Plant Physiol. 165, 791–809. doi: 10.1104/pp.114.239665
Buswell, W., Schwarzenbacher, R. E., Luna, E., Sellwood, M., Chen, B., Flors, V., et al. (2018). Chemical priming of immunity without costs to plant growth. New Phytol. 218, 1205–1216. doi: 10.1111/nph.15062
Caarls, L., Elberse, J., Awwanah, M., Ludwig, N. R., De Vries, M., Zeilmaker, T., et al. (2017). Arabidopsis JASMONATE-INDUCED OXYGENASES down-regulate plant immunity by hydroxylation and inactivation of the hormone jasmonic acid. Proc. Natl. Acad Sci. U.S.A. 114, 6388–6393. doi: 10.1073/pnas.1701101114
Cahill, D. M., Rookes, J. E., Wilson, B. A., Gibson, L., and Mcdougall, K. L. (2008). Phytophthora cinnamomi and Australia’s biodiversity: impacts, predictions and progress towards control. Aust. J. Bot. 56, 279–310. doi: 10.1071/BT07159
Chen, H., Lai, Z., Shi, J., Xiao, Y., Chen, Z., and Xu, X. (2010). Roles of Arabidopsis WRKY18, WRKY40 and WRKY60 transcription factors in plant responses to abscisic acid and abiotic stress. BMC Plant Biol. 10:281. doi: 10.1186/1471-2229-10-281
Chen, K., Du, L., and Chen, Z. (2003). Sensitization of defense responses and activation of programmed cell death by a pathogen-induced receptor-like protein kinase in Arabidopsis. Plant Mol. Biol. 53, 61–74. doi: 10.1023/b:plan.0000009265.72567.58
Cheng, Y., Zhou, Y., Yang, Y., Chi, Y. J., Zhou, J., Chen, J. Y., et al. (2012). Structural and functional analysis of VQ motif-containing proteins in Arabidopsis as interacting proteins of WRKY transcription factors. Plant Physiol. 159, 810–825. doi: 10.1104/pp.112.196816
Chilley, P. M., Casson, S. A., Tarkowski, P., Hawkins, N., Wang, K. L. C., Hussey, P. J., et al. (2006). The POLARIS peptide of Arabidopsis regulates auxin transport and root growth via effects on ethylene signaling. Plant Cell 18, 3058–3072. doi: 10.1105/tpc.106.040790
Choi, D. S., and Hwang, B. K. (2011). Proteomics and functional analyses of pepper abscisic acid–responsive 1 (ABR1), which is involved in cell death and defense signaling. Plant Cell 23, 823–842. doi: 10.1105/tpc.110.082081
Cluzet, S., Torregrosa, C., Jacquet, C., Lafitte, C., Fournier, J., Mercier, L., et al. (2004). Gene expression profiling and protection of Medicago truncatula against a fungal infection in response to an elicitor from green algae Ulva spp. Plant Cell Environ 27, 917–928. doi: 10.1111/j.1365-3040.2004.01197.x
Cook, J., Zhang, J., Norrie, J., Blal, B., and Cheng, Z. (2018). Seaweed extract (Stella Maris®) activates innate immune responses in Arabidopsis thaliana and protects host against bacterial pathogens. Mar. Drugs 16:221. doi: 10.3390/md16070221
Cooper, A., Latunde-Dada, A., Woods-Tör, A., Lynn, J., Lucas, J., Crute, I., et al. (2008). Basic compatibility of Albugo candida in Arabidopsis thaliana and Brassica juncea causes broad-spectrum suppression of innate immunity. Mol. Plant Microbe. Interact. 21, 745–756. doi: 10.1094/MPMI-21-6-0745
Costa, R., Domínguez, A., and Choupina, A. (2020). Cloning and expression analysis of an endo-1, 3-β-d-glucosidase from Phytophthora cinnamomi. Mol. Biol. Rep. 47, 935–942. doi: 10.1007/s11033-019-05185-9
De Lorenzo, G., Ferrari, S., Giovannoni, M., Mattei, B., and Cervone, F. (2019). Cell wall traits that influence plant development, immunity, and bioconversion. Plant J. 97, 134–147. doi: 10.1111/tpj.14196
de Oliveira, M. V., Xu, G., Li, B., De Souza Vespoli, L., Meng, X., Chen, X., et al. (2016). Specific control of Arabidopsis BAK1/SERK4-regulated cell death by protein glycosylation. Nat. Plants 2:15218. doi: 10.1038/nplants.2015.218
Devoto, A., Hartmann, H. A., Piffanelli, P., Elliott, C., Simmons, C., Taramino, G., et al. (2003). Molecular phylogeny and evolution of the plant-specific seven-transmembrane MLO family. J. Mol. Evol. 56, 77–88. doi: 10.1007/s00239-002-2382-5
Dong, X., Hong, Z., Chatterjee, J., Kim, S., and Verma, D. (2008). Expression of callose synthase genes and its connection with Npr1 signaling pathway during pathogen infection. Planta 229, 87–98. doi: 10.1007/s00425-008-0812-3
Engelbrecht, J., Duong, T. A., and Van Den Berg, N. (2013). Development of a nested quantitative real-time PCR for detecting Phytophthora cinnamomi in Persea americana rootstocks. Plant Dis. 97, 1012–1017. doi: 10.1094/PDIS-11-12-1007-RE
Eshraghi, L., Anderson, J., Aryamanesh, N., Shearer, B., Mccomb, J., Hardy, G. S., et al. (2011). Phosphite primed defence responses and enhanced expression of defence genes in Arabidopsis thaliana infected with Phytophthora cinnamomi. Plant Pathol. 60, 1086–1095. doi: 10.1111/j.1365-3059.2011.02471.x
GonñI, O., Fort, A., Quille, P., Mckeown, P. C., Spillane, C., and O’connell, S. (2016). Comparative transcriptome analysis of two Ascophyllum nodosum extract biostimulants: same seaweed but different. J. Agric. Food Chem. 64, 2980–2989. doi: 10.1021/acs.jafc.6b00621
Gunupuru, L., Patel, J., Sumarah, M., Renaud, J., Mantin, E., and Prithiviraj, B. (2019). A plant biostimulant made from the marine brown algae Ascophyllum nodosum and chitosan reduce Fusarium head blight and mycotoxin contamination in wheat. PLoS One 14:220562. doi: 10.1371/journal.pone.0220562
Guo, D. D., Liu, F., Tu, Y. H., He, B. X., Gao, Y., and Guo, M. L. (2016). Expression patterns of three UGT genes in different chemotype safflower lines and under MeJA stimulus revealed their potential role in flavonoid biosynthesis. PLoS One 11:e0158159. doi: 10.1371/journal.pone.0158159
Herrera-Vásquez, A., Carvallo, L., Blanco, F., Tobar, M., Villarroel-Candia, E., Vicente-Carbajosa, J., et al. (2015). Transcriptional control of glutaredoxin GRXC9 expression by a salicylic acid-dependent and NPR1-independent pathway in Arabidopsis. Plant Mol. Biol. Rep. 33, 624–637. doi: 10.1007/s11105-014-0782-5
Huffaker, A., Pearce, G., and Ryan, C. A. (2006). An endogenous peptide signal in Arabidopsis activates components of the innate immune response. Proc. Natl. Acad Sci. U.S.A. 103, 10098–10103. doi: 10.1073/pnas.0603727103
Islam, M., Hussain, H., Rookes, J., and Cahill, D. (2018). Transcriptome analysis, using RNA-Seq of Lomandra longifolia roots infected with Phytophthora cinnamomi reveals the complexity of the resistance response. Plant Biol. 20, 130–142. doi: 10.1111/plb.12624
Islam, M. T., Rookes, J. E., and Cahill, D. M. (2017). Active defence by an Australian native host, Lomandra longifolia, provides resistance against Phytophthora cinnamomi. Funct. Plant Biol. 44, 386–399. doi: 10.1071/FP16266
Jithesh, M., Shukla, P. S., Kant, P., Joshi, J., Critchley, A. T., and Prithiviraj, B. (2019). Physiological and transcriptomics analyses reveal that Ascophyllum nodosum extracts induce salinity tolerance in Arabidopsis by regulating the expression of stress responsive genes. J. Plant Growth Regul. 38, 463–478. doi: 10.1007/s00344-018-9861-4
Kamarudin, A. N., Idris, A., and Balia Yusof, Z. (2017). Thiamine biosynthesis gene expression analysis in Elaeis guineensis during interactions with Hendersonia toruloidea. J. Oil Palm Res. 29, 218–226. doi: 10.21894/jopr.2017.2902.06
Kasirajan, L., Hoang, N. V., Furtado, A., Botha, F. C., and Henry, R. J. (2018). Transcriptome analysis highlights key differentially expressed genes involved in cellulose and lignin biosynthesis of sugarcane genotypes varying in fiber content. Sci. Rep. 8, 1–16. doi: 10.1038/s41598-018-30033-4
Kerchev, P., Van Der Meer, T., Sujeeth, N., Verlee, A., Stevens, C. V., Van Breusegem, F., et al. (2019). Molecular priming as an approach to induce tolerance against abiotic and oxidative stresses in crop plants. Biotechnol. Adv. 40:107503. doi: 10.1016/j.biotechadv.2019.107503
Khan, W., Rayirath, U. P., Subramanian, S., Jithesh, M. N., Rayorath, P., Hodges, D. M., et al. (2009). Seaweed extracts as biostimulants of plant growth and development. J. Plant Growth Regul. 28, 386–399. doi: 10.1007/s00344-009-9103-x
Knoth, C., and Eulgem, T. (2008). The oomycete response gene LURP1 is required for defense against Hyaloperonospora parasitica in Arabidopsis thaliana. Plant J. 55, 53–64. doi: 10.1111/j.1365-313X.2008.03486.x
Kunz, B. A., Dando, P. K., Grice, D. M., Mohr, P. G., Schenk, P. M., and Cahill, D. M. (2008). UV-induced DNA damage promotes resistance to the biotrophic pathogen Hyaloperonospora parasitica in Arabidopsis. Plant Physiol. 148, 1021–1031. doi: 10.1104/pp.108.125435
Le Berre, J. Y., Gourgues, M., Samans, B., Keller, H., Panabières, F., and Attard, A. (2017). Transcriptome dynamic of Arabidopsis roots infected with Phytophthora parasitica identifies VQ29, a gene induced during the penetration and involved in the restriction of infection. PLoS One 12:190341. doi: 10.1371/journal.pone.0190341
Lee, D. S., Kim, Y. C., Kwon, S. J., Ryu, C. M., and Park, O. K. (2017). The Arabidopsis cysteine-rich receptor-like kinase CRK36 regulates immunity through interaction with the cytoplasmic kinase BIK1. Front. Plant Sci. 8:1856. doi: 10.3389/fpls.2017.01856
Lemarié, S., Robert-Seilaniantz, A., Lariagon, C., Lemoine, J., Marnet, N., Jubault, M., et al. (2015). Both the jasmonic acid and the salicylic acid pathways contribute to resistance to the biotrophic clubroot agent Plasmodiophora brassicae in Arabidopsis. Plant Cell Physiol. 56, 2158–2168. doi: 10.1093/pcp/pcv127
Lu, H., Rate, D. N., Song, J. T., and Greenberg, J. T. (2003). ACD6, a novel ankyrin protein, is a regulator and an effector of salicylic acid signaling in the Arabidopsis defense response. Plant Cell 15, 2408–2420. doi: 10.1105/tpc.015412
Ma, W., Smigel, A., Tsai, Y. C., Braam, J., and Berkowitz, G. A. (2008). Innate immunity signaling: cytosolic Ca2+ elevation is linked to downstream nitric oxide generation through the action of calmodulin or a calmodulin-like protein. Plant Physiol. 148, 818–828. doi: 10.1104/pp.108.125104
Martinez-Medina, A., Flors, V., Heil, M., Mauch-Mani, B., Pieterse, C. M., Pozo, M. J., et al. (2016). Recognizing plant defense priming. Trends Plant Sci. 21, 818–822. doi: 10.1016/j.tplants.2016.07.009
Mattner, S. W., Milinkovic, M., and Arioli, T. (2018). Increased growth response of strawberry roots to a commercial extract from Durvillaea potatorum and Ascophyllum nodosum. J. Appl. Phycol. 30, 2943–2951. doi: 10.1007/s10811-017-1387-9
Mattner, S. W., Wite, D., Riches, D. A., Porter, I. J., and Arioli, T. (2013). The effect of kelp extract on seedling establishment of broccoli on contrasting soil types in southern Victoria, Australia. Biol. Agric. Hortic. 29, 258–270. doi: 10.1080/01448765.2013.830276
Meier, S., Ruzvidzo, O., Morse, M., Donaldson, L., Kwezi, L., and Gehring, C. (2010). The Arabidopsis wall associated kinase-like 10 gene encodes a functional guanylyl cyclase and is co-expressed with pathogen defense related genes. PLoS One 5:8904. doi: 10.1371/journal.pone.0008904
Mhamdi, A., Hager, J., Chaouch, S., Queval, G., Han, Y., Taconnat, L., et al. (2010). Arabidopsis GLUTATHIONE REDUCTASE1 plays a crucial role in leaf responses to intracellular hydrogen peroxide and in ensuring appropriate gene expression through both salicylic acid and jasmonic acid signaling pathways. Plant Physiol. 153, 1144–1160. doi: 10.1104/pp.110.153767
Mir, R., Hernández, M. L., Abou-Mansour, E., Martínez-Rivas, J. M., Mauch, F., Métraux, J.-P., et al. (2013). Pathogen and Circadian Controlled 1 (PCC1) regulates polar lipid content, ABA-related responses, and pathogen defence in Arabidopsis thaliana. J. Exp. Bot. 64, 3385–3395. doi: 10.1093/jxb/ert177
Nair, P., Kandasamy, S., Zhang, J., Ji, X., Kirby, C., Benkel, B., et al. (2012). Transcriptional and metabolomic analysis of Ascophyllum nodosum mediated freezing tolerance in Arabidopsis thaliana. BMC Genomics 13:643. doi: 10.1186/1471-2164-13-643
Naithani, S., Preece, J., D’eustachio, P., Gupta, P., Amarasinghe, V., Dharmawardhana, P. D., et al. (2016). Plant Reactome: a resource for plant pathways and comparative analysis. Nucleic Acids Res. 45, 1029–1039. doi: 10.1093/nar/gkw932
Najafi, J., Brembu, T., Vie, A. K., Viste, R., Winge, P., Somssich, I. E., et al. (2020). PAMP-INDUCED SECRETED PEPTIDE 3 modulates immunity in Arabidopsis. J. Exp. Bot. 71, 850–864. doi: 10.1093/jxb/erz482
Návarová, H., Bernsdorff, F., Döring, A.-C., and Zeier, J. (2012). Pipecolic acid, an endogenous mediator of defense amplification and priming, is a critical regulator of inducible plant immunity. Plant Cell 24, 5123–5141. doi: 10.1105/tpc.112.103564
Omidbakhshfard, M. A., Sujeeth, N., Gupta, S., Omranian, N., Guinan, K. J., Brotman, Y., et al. (2020). A Biostimulant Obtained from the Seaweed Ascophyllum nodosum protects Arabidopsis thaliana from severe oxidative stress. Int. J. Mol. Sci. 21:474. doi: 10.3390/ijms21020474
Pye, V. E., Christensen, C. E., Dyer, J. H., Arent, S., and Henriksen, A. (2010). Peroxisomal plant 3-ketoacyl-CoA thiolase structure and activity are regulated by a sensitive redox switch. J. Biol. Chem. 285, 24078–24088. doi: 10.1074/jbc.M110.106013
Reeksting, B. J., Olivier, N. A., and Van Den Berg, N. (2016). Transcriptome responses of an ungrafted Phytophthora root rot tolerant avocado (Persea americana) rootstock to flooding and Phytophthora cinnamomi. BMC Plant Biol. 16:205. doi: 10.1186/s12870-016-0893-2
Ried, M. K., Banhara, A., Hwu, F. Y., Binder, A., Gust, A. A., Höfle, C., et al. (2019). A set of Arabidopsis genes involved in the accommodation of the downy mildew pathogen Hyaloperonospora arabidopsidis. PLoS Pathog. 15:7747. doi: 10.1371/journal.ppat.1007747
Roetschi, A., Si-Ammour, A., Belbahri, L., Mauch, F., and Mauch-Mani, B. (2001). Characterization of an Arabidopsis–Phytophthora pathosystem: resistance requires a functional PAD2 gene and is independent of salicylic acid, ethylene and jasmonic acid signalling. Plant J. 28, 293–305. doi: 10.1046/j.1365-313x.2001.01148.x
Rojas, C. M., Senthil-Kumar, M., Tzin, V., and Mysore, K. (2014). Regulation of primary plant metabolism during plant-pathogen interactions and its contribution to plant defense. Front. Plant Sci. 5:17. doi: 10.3389/fpls.2014.00017
Rookes, J. E., Wright, M. L., and Cahill, D. M. (2008). Elucidation of defence responses and signalling pathways induced in Arabidopsis thaliana following challenge with Phytophthora cinnamomi. Physiol. Mol. Plant 72, 151–161. doi: 10.1016/j.pmpp.2008.08.005
Saini, S., Sharma, I., and Pati, P. K. (2015). Versatile roles of brassinosteroid in plants in the context of its homoeostasis, signaling and crosstalks. Front. Plant Sci. 6:950. doi: 10.3389/fpls.2015.00950
Santaniello, A., Scartazza, A., Gresta, F., Loreti, E., Biasone, A., Di Tommaso, D., et al. (2017). Ascophyllum nodosum seaweed extract alleviates drought stress in Arabidopsis by affecting photosynthetic performance and related gene expression. Front. Plant Sci. 8:1362. doi: 10.3389/fpls.2017.01362
Schaller, F., Biesgen, C., Müssig, C., Altmann, T., and Weiler, E. W. (2000). 12-Oxophytodienoate reductase 3 (OPR3) is the isoenzyme involved in jasmonate biosynthesis. Planta 210, 979–984. doi: 10.1007/s004250050706
Sels, J., Mathys, J., De Coninck, B. M., Cammue, B. P., and De Bolle, M. F. (2008). Plant pathogenesis-related (PR) proteins: a focus on PR peptides. Plant Physiol. Biochem. 46, 941–950. doi: 10.1016/j.plaphy.2008.06.011
Seo, H. S., Song, J. T., Cheong, J. J., Lee, Y. H., Lee, Y. W., Hwang, I., et al. (2001). Jasmonic acid carboxyl methyltransferase: a key enzyme for jasmonate-regulated plant responses. Proc. Natl. Acad Sci. U.S.A. 98, 4788–4793. doi: 10.1073/pnas.081557298
Serrazina, S., Santos, C., Machado, H., Pesquita, C., Vicentini, R., Pais, M. S., et al. (2015). Castanea root transcriptome in response to Phytophthora cinnamomi challenge. Tree Genet. Genomes 11:6. doi: 10.1007/s11295-014-0829-7
Shine, M., Yang, J. W., El-Habbak, M., Nagyabhyru, P., Fu, D. Q., Navarre, D., et al. (2016). Cooperative functioning between phenylalanine ammonia lyase and isochorismate synthase activities contributes to salicylic acid biosynthesis in soybean. New Phytol. 212, 627–636. doi: 10.1111/nph.14078
Shukla, P. S., Borza, T., Critchley, A. T., and Prithiviraj, B. (2016). Carrageenans from red seaweeds as promoters of growth and elicitors of defense response in plants. Front. Mar. Sci. 3:81. doi: 10.3389/fmars.2016.00081
Shukla, P. S., Mantin, E. G., Adil, M., Bajpai, S., Critchley, A. T., and Prithiviraj, B. (2019). Ascophyllum nodosum-based biostimulants: sustainable applications in agriculture for the stimulation of plant growth, stress tolerance, and disease management. Front. Plant Sci. 10:655. doi: 10.3389/fpls.2019.00655
Simon, C., Langlois-Meurinne, M., Didierlaurent, L., Chaouch, S., Bellvert, F., Massoud, K., et al. (2014). The secondary metabolism glycosyltransferases UGT73B3 and UGT73B5 are components of redox status in resistance of Arabidopsis to Pseudomonas syringae pv. tomato. Plant Cell Environ. 37, 1114–1129. doi: 10.1111/pce.12221
Song, J. T., Lu, H., Mcdowell, J. M., and Greenberg, J. T. (2004). A key role for ALD1 in activation of local and systemic defenses in Arabidopsis. Plant J. 40, 200–212. doi: 10.1111/j.1365-313X.2004.02200.x
Song, W., Qi, N., Liang, C., Duan, F., and Zhao, H. (2019). Soybean root transcriptome profiling reveals a nonhost resistant response during Heterodera glycines infection. PLoS One 14:e0217130. doi: 10.1371/journal.pone.0217130
Stepanova, A. N., Robertson-Hoyt, J., Yun, J., Benavente, L. M., Xie, D.-Y., Doležal, K., et al. (2008). TAA1-mediated auxin biosynthesis is essential for hormone crosstalk and plant development. Cell 133, 177–191. doi: 10.1016/j.cell.2008.01.047
Sun, J., Chen, Q., Qi, L., Jiang, H., Li, S., Xu, Y., et al. (2011). Jasmonate modulates endocytosis and plasma membrane accumulation of the Arabidopsis PIN2 protein. New Phytol. 191, 360–375. doi: 10.1111/j.1469-8137.2011.03713.x
Takagi, H., Ishiga, Y., Watanabe, S., Konishi, T., Egusa, M., Akiyoshi, N., et al. (2016). Allantoin, a stress-related purine metabolite, can activate jasmonate signaling in a MYC2-regulated and abscisic acid-dependent manner. J. Exp. Bot. 67, 2519–2532. doi: 10.1093/jxb/erw071
Thimm, O., Bläsing, O., Gibon, Y., Nagel, A., Meyer, S., Krüger, P., et al. (2004). MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 37, 914–939. doi: 10.1111/j.1365-313x.2004.02016.x
Thordal-Christensen, H., Zhang, Z., Wei, Y., and Collinge, D. B. (1997). Subcellular localization of H2O2 in plants. H2O2 accumulation in papillae and hypersensitive response during the barley-powdery mildew interaction. Plant J. 11, 1187–1194. doi: 10.1046/j.1365-313X.1997.11061187.x
Tommasini, L., Svensson, J. T., Rodriguez, E. M., Wahid, A., Malatrasi, M., Kato, K., et al. (2008). Dehydrin gene expression provides an indicator of low temperature and drought stress: transcriptome-based analysis of barley (Hordeum vulgare L.). Funct. Integr. Genom. 8, 387–405. doi: 10.1007/s10142-008-0081-z
Venus, Y., and Oelmüller, R. (2013). Arabidopsis ROP1 and ROP6 influence germination time, root morphology, the formation of F-actin bundles, and symbiotic fungal interactions. Mol. Plant. 6, 872–886. doi: 10.1093/mp/sss101
Wang, Y., Meng, Y., Zhang, M., Tong, X., Wang, Q., Sun, Y., et al. (2011). Infection of Arabidopsis thaliana by Phytophthora parasitica and identification of variation in host specificity. Mol. Plant Pathol. 12, 187–201. doi: 10.1111/j.1364-3703.2010.00659.x
Wite, D., Mattner, S., Porter, I., and Arioli, T. (2015). The suppressive effect of a commercial extract from Durvillaea potatorum and Ascophyllum nodosum on infection of broccoli by Plasmodiophora brassicae. J. Appl. Phycol. 27, 2157–2161. doi: 10.1007/s10811-015-0564-y
Xie, Y. D., Li, W., Guo, D., Dong, J., Zhang, Q., Fu, Y., et al. (2010). The Arabidopsis gene SIGMA FACTOR-BINDING PROTEIN 1 plays a role in the salicylate-and jasmonate-mediated defence responses. Plant Cell Environ. 33, 828–839. doi: 10.1111/j.1365-3040.2009.02109.x
Xin, X.-F., Nomura, K., Ding, X., Chen, X., Wang, K., Aung, K., et al. (2015). Pseudomonas syringae effector avirulence protein E localizes to the host plasma membrane and down-regulates the expression of the nonrace-specific disease resistance1/harpin-induced1-like13 gene required for antibacterial immunity in Arabidopsis. Plant Physiol. 169, 793–802. doi: 10.1104/pp.15.00547
Xu, J., Li, Y., Wang, Y., Liu, H., Lei, L., Yang, H., et al. (2008). Activation of MAPK kinase 9 induces ethylene and camalexin biosynthesis and enhances sensitivity to salt stress in Arabidopsis. J. Biol. Chem. 283, 26996–27006. doi: 10.1074/jbc.M801392200
Yang, Y., Zhang, Y., Ding, P., Johnson, K., Li, X., and Zhang, Y. (2012). The ankyrin-repeat transmembrane protein BDA1 functions downstream of the receptor-like protein SNC2 to regulate plant immunity. Plant Physiol. 159, 1857–1865. doi: 10.1104/pp.112.197152
Yeh, Y. H., Panzeri, D., Kadota, Y., Huang, Y. C., Huang, P. Y., Tao, C. N., et al. (2016). The Arabidopsis malectin-like/LRR-RLK IOS1 is critical for BAK1-dependent and BAK1-independent pattern-triggered immunity. Plant Cell 28, 1701–1721. doi: 10.1105/tpc.16.00313
Zhang, H., Huang, L., Dai, Y., Liu, S., Hong, Y., Tian, L., et al. (2015). Arabidopsis AtERF15 positively regulates immunity against Pseudomonas syringae pv. tomato DC3000 and Botrytis cinerea. Front. Plant Sci. 6:686. doi: 10.3389/fpls.2015.00686
Zhang, J., Li, W., Xiang, T., Liu, Z., Laluk, K., Ding, X., et al. (2010). Receptor-like cytoplasmic kinases integrate signaling from multiple plant immune receptors and are targeted by a Pseudomonas syringae effector. Cell Host Microb. 7, 290–301. doi: 10.1016/j.chom.2010.03.007
Keywords: seaweed, Ascophyllum nodosum, Durvillaea potatorum, RNA-Seq, Phytophthora cinnamomi, Arabidopsis thaliana
Citation: Islam MT, Gan HM, Ziemann M, Hussain HI, Arioli T and Cahill D (2020) Phaeophyceaean (Brown Algal) Extracts Activate Plant Defense Systems in Arabidopsis thaliana Challenged With Phytophthora cinnamomi. Front. Plant Sci. 11:852. doi: 10.3389/fpls.2020.00852
Received: 31 March 2020; Accepted: 27 May 2020;
Published: 07 July 2020.
Edited by:Paolo Carletti, University of Padua, Italy
Reviewed by:Dongdong Niu, Nanjing Agricultural University, China
Alan T. Critchley, Cape Breton University, Canada
Copyright © 2020 Islam, Gan, Ziemann, Hussain, Arioli and Cahill. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: David Cahill, firstname.lastname@example.org