Original Research ARTICLE
Partial Activation of SA- and JA-Defensive Pathways in Strawberry upon Colletotrichum acutatum Interaction
- 1Departamento de Bioquímica y Biología Molecular e Instituto Andaluz de Biotecnología, Edificio Severo Ochoa (C6), Universidad de Córdoba, Córdoba, Spain
- 2Instituto de Biología Molecular y Celular de Plantas, Universidad Politécnica de Valencia–Consejo Superior de Investigaciones Científicas, Valencia, Spain
- 3Departamento de Arquitectura de Computación, Universidad de Málaga, Campus de Teatinos, Málaga, Spain
- 4Centro Andalusian Institute of Agricultural and Fishering Research and Training (IFAPA) Las Torres-Tomejil, CAPMA–Junta de Andalucía, Sevilla, Spain
- 5Departamento de Biología Vegetal, Facultad de Ciencias, Universidad de Málaga, Campus de Teatinos, Málaga, Spain
Understanding the nature of pathogen host interaction may help improve strawberry (Fragaria × ananassa) cultivars. Plant resistance to pathogenic agents usually operates through a complex network of defense mechanisms mediated by a diverse array of signaling molecules. In strawberry, resistance to a variety of pathogens has been reported to be mostly polygenic and quantitatively inherited, making it difficult to associate molecular markers with disease resistance genes. Colletotrichum acutatum spp. is a major strawberry pathogen, and completely resistant cultivars have not been reported. Moreover, strawberry defense network components and mechanisms remain largely unknown and poorly understood. Assessment of the strawberry response to C. acutatum included a global transcript analysis, and acidic hormones SA and JA measurements were analyzed after challenge with the pathogen. Induction of transcripts corresponding to the SA and JA signaling pathways and key genes controlling major steps within these defense pathways was detected. Accordingly, SA and JA accumulated in strawberry after infection. Contrastingly, induction of several important SA, JA, and oxidative stress-responsive defense genes, including FaPR1-1, FaLOX2, FaJAR1, FaPDF1, and FaGST1, was not detected, which suggests that specific branches in these defense pathways (those leading to FaPR1-2, FaPR2-1, FaPR2-2, FaAOS, FaPR5, and FaPR10) were activated. Our results reveal that specific aspects in SA and JA dependent signaling pathways are activated in strawberry upon interaction with C. acutatum. Certain described defense-associated transcripts related to these two known signaling pathways do not increase in abundance following infection. This finding suggests new insight into a specific putative molecular strategy for defense against this pathogen.
Strawberry fruit (Fragaria × ananassa Duch.) is of great importance throughout the world due to its flavor and nutritious qualities [FAOSTAT (http://faostat.fao.org/)] (Freeman et al., 2001). However, strawberry exhibits wide diversity in its susceptibility to a large variety of phytopathogenic organisms, including Colletotrichum spp., which are major pathogens of this crop (Simpson, 1991; Maas, 1998), requiring excessive use of chemical agents for disease control. Breeding for resistance by crossing in natural resistance mechanisms found in related genotypes lead to more sustainable farming with fewer chemical inputs. However, resistance to a variety of pathogens is generally polygenic and quantitatively inherited (Amil-Ruiz et al., 2011), making it difficult to associate single molecular markers with disease resistance genes.
Understanding the molecular interplay between plant and microbes has led to identification of candidate genes that have been used in developing transgenic strategies and breeding efforts to increase resistance against specific pathogens in many plants (Tohidfar and Khosravi, 2015). Plant resistance to pathogenic agents usually operates through a complex defense mechanism network. Compounds such as salicylic acid (SA), jasmonic acid (JA), and ethylene (ET) regulate plant defense pathways to trigger appropriate responses to different pathogens (Robert-Seilaniantz et al., 2011). Whereas the SA signaling pathway is mainly activated against biotrophic pathogens, the JA/ET signaling pathway is activated against necrotrophic pathogens. Interplay and antagonism between these signaling pathways has also been described (Robert-Seilaniantz et al., 2011). Pathogens have evolved to overcome plant immunity mechanisms by disrupting the fine crosstalk between these defense pathways (El Oirdi et al., 2011; Rahman et al., 2012; Chung et al., 2013).
In Arabidopsis, genetic markers have been identified to denote activation of canonical SA- and JA defense pathways. Both PAD4 and EDS1 activate SA biosynthesis, and both are also SA induced (Rustérucci et al., 2001; Venugopal et al., 2009). PR1 and PR2 are known pathogen- and SA-responsive genes, which are well-established markers for the Arabidopsis defense responses against P. syringae (Uknes et al., 1992). WRKY70 encodes a transcription factor and is an important regulator in the interplay of SA- and JA-related plant defense responses (Li et al., 2004). AOS encodes an allene oxide synthase, which is responsible for JA synthesis, and its expression is also regulated by this hormone (Turner et al., 2002). PDF1.2 and JAR1 are genes that are commonly used to monitor JA responses (Staswick and Tiryaki, 2004; Pieterse et al., 2009), whereas LOX2, which is involved in JA biosynthesis, is activated by a positive feedback loop (Sasaki et al., 2001). GST1 is part of the array of defense-related genes induced in response to oxidative burst produced after pathogen infection (Bhattacharjee, 2012). This gene encodes a glutathione S-transferase that is known to play a key role in reactive oxygen species (ROS) detoxification and reduction.
In strawberry, isolation of individual genes related to plant defense has been previously reported (Amil-Ruiz et al., 2011). Furthermore, Casado-Díaz et al. (2006) first reported isolation of a large set of genes with altered expression during strawberry and C. acutatum interaction. Over the last decade, microarrays have proved to be a valuable tool to analyze the expression of thousands of genes on a simultaneous basis, contributing to elucidate the underlying networks of gene regulation that lead to a wide variety of defense responses. Thus, Maleck et al. (2000) have provided a comprehensive description of the SAR genes from Arabidopsis thaliana. Wang et al. (2005) reported that SAR requires protein secretory pathway induction. Further, microarrays have led to numerous findings of key regulatory genes for defense signaling as well as valuable end-point genes whose products display direct action against pest and diseases (Wang et al., 2006; Sarowar et al., 2011). In addition, microarray analysis has demonstrated a substantial crosstalk among different defense signaling pathways (Schenk et al., 2003) and ultimately among genes and their products, and the entire pathways are not always tuned by signaling (Lodha and Basak, 2012). In strawberry, microarrays have also been used to analyze gene expression differences between white and red fruit after 24 h interaction with C. acutatum and provision of certain data from factors associated with pathogen quiescence during fruit immature stages (Guidarelli et al., 2011). However, strawberry defense network components remain largely unknown or poorly understood, and accurate mechanisms remain elusive.
In this report, a transcriptomic approach has been used to identify pathogen-responsive genes in F × ananassa strawberry crown and petiole, key sites of C. acutatum infection. Transcript accumulation was also monitored in response to SA and JA treatments. Synthesis of SA and JA signaling molecules was also analyzed. These trials reveal changes in the transcriptomic profile of plants challenged with C. acutatum and present a hypothesis about a strategy potentially used by this pathogen to overcome strawberry defense. These data deepen our understanding of the complex genetic and molecular mechanisms of strawberry defense.
Materials and Methods
Plant Materials, Pathogen Inoculation, and Hormonal Treatments
Plant culture (Fragaria × ananassa cultivar Camarosa) and growth conditions, C. acutatum (isolate CECT 20240) inoculation, and treatments with chemicals have been previously described (Encinas-Villarejo et al., 2009). In breaf, 8-week-old strawberry plantlets were placed in 20 cm diameter plastic pots containing sterilized peat and grown for a minimum of 6 additional weeks prior to mock or pathogen inoculation by spraying a spore suspension of 104 conidia·ml−1. Crown was collected 5 days after treatment (spray-infected and mock-treated) for microarray studies. At this time, under our experimental conditions, plants still looked healthy, and no visible disease symptoms were easily detected, even in petioles, crowns or leaves. For RTqPCR analysis, crowns and petioles were collected 1, 3, 5, 7, and 9 days after treatment, as previously described in Casado-Díaz et al. (2006). For treatment and hormonal content analysis, axenic in-vitro plantlets were used and aseptically sprayed with either MeJa (2 mM) or SA (5 mM) solutions, or inoculated with C. acutatum conidia suspension (104 conidia·ml−1), respectively. All collected plants were immediately frozen in liquid nitrogen and stored at −80°C until use. For transmission electron microscopy, pathogen was spot inoculated in crown by applying 50-μL droplets of conidia suspension (106 conidia·ml−1) (Arroyo et al., 2005). Control plants were similarly inoculated with 50 μL of sterile distilled water. The position of each infection site was marked for reference. Inoculated plants were enclosed in plastic bags for 48 h to maintain high relative humidity and incubated in similar conditions as described above. Hence, disease is forced to progress quicker than by spraying.
Light Microscopic Fungal Development Observation
Light microscopy analyses of C. acutatum development were performed on strawberry leaf discs (10 mm diameter) randomly excised from infected leaflet at 1, 3, 5, 7, 9 dpi, using a method modified from Debode et al. (2009). In breaf, leaf discs were cleared in 0.15% trichloroacetic acid (TCA) in a 3:1 (v/v) mixture of ethanol and chloroform for 48 h, with at least three changes of the bleaching solution, rinsed briefly in lactoglycerol and incubated at room temperature for 1 h in lactophenol blue (Sigma), and washed 3 times in lactoglycerol. For ROS detection, leaf discs were infiltrated in a 10 mg/ml DAB (3,3′ -diaminobenzidine, Sigma) solution for 10 min and subsequently incubated overnight at room temperature in dark, and cleared for 24 h in 3:1 (v/v) ethanol: glacial acetic acid with three changes. Treated leaf disc was mounted in 50% fresh glycerol on glass slides and examined using a Leica DM5000B microscope. Images were captured with a Leica DC500 digital camera. Three leaftlets were sampled from each plant, and 3 plants were observed at each time point. The overall number of conidia (germinated and non-germinated) and appresoria were counted per leaf disc taken up to 9 dpi.
Transmission Electron Microscopy (TEM)
Strips of tissue ~1 mm thick and 1–2 mm long were removed from beneath the inoculation droplets and fixed in 4% (v/v) glutaraldehyde in 0.1 mol·L –1cacodylate buffer (pH 7.2) for 3 h at 4°C. Upon rinsing in the same buffer, tissues were post-fixed in 1% (w/v) osmium tetroxide for 2 h at 4°C and subsequently dehydrated in a graded acetone series and embedded in EMBED-812 (Polysciences, Warrington, Penn.) according to manufacturer's instructions. Slides with semithin sections (0.5 μm) were placed on a hotplate at 50°C, stained for 1 min with 0.1% aqueous toluidine blue O and examined using a light microscope (Leitz Aristoplan). Ultrathin 60–80 nm sections were made with a Reichert-Jung Ultracut E ultramicrotome and a diamond knife, and collected on 300-mesh copper grids (Dashek and Mayfield 2000). Grids were stained with 7% aqueous uranyl acetate and lead citrate. Sections were observed and images were collected using a Philips CM-10 transmission electron microscope (TEM).
Total RNA Extraction and Real-Time qPCR
Total RNA from strawberry tissues was isolated as described previously (Casado-Díaz et al., 2006), treated with DnaseI (Invitrogen) for residual DNA removal, and further purified with the RNeasy MinElute Cleanup Kit (QIAGEN). Purified RNA was quantified by NanoDrop 1000 Spectrophotometer (Thermo scientific). RNA integrity was checked using the Agilent 2100 Bioanalyzer (Agilent Technologies, Deutschland). First-strand cDNA synthesis was carried out using 1 μg of purified total RNA as template for a 20 μL reaction [iScript cDNA Synthesis kit (Bio-Rad)]. RT reactions were diluted 5-fold with nuclease-free water prior to qPCR.
Specific primer pairs set were designed using Oligo Primer Analysis software version 6.65, tested by dissociation curve analysis, and verified for absence of non-specific amplification (Table S1). FaGAPDH2 gene was used for normalization (Khan and Shih, 2004; Amil-Ruiz et al., 2013). RTqPCR runs were performed using two technical replicates in the same run and three biological replicates in different runs, as described previously (Encinas-Villarejo et al., 2009), using SsoAdvanced™ SYBR® Green supermix, and MyIQ v1.004 and iCycler v3.1 real-time PCR systems (Bio-Rad).
Microarray Analysis and Strawberry Gene Annotation
For microarray analysis, strawberry samples were collected 5 days after treatment (spray-infected and mock-treated). Crowns were used to make biological replicates, and overall RNA was isolated from three independent biological replicates for hybridization against a proprietary microarray representing approximately 2529 predicted unigenes from F. vesca (Shulaev et al., 2011), previously identified from strawberry libraries (Casado-Díaz et al., 2006; and JL Caballero unpublished). Microarray data with accession GSE56296 were deposited in the NCBI Gene Expression Omnibus. Quality control, labeling, hybridization, and scanning were carried out by the SCAI, University of Córdoba (http://www.uco.es/servicios/scai/index.html), following the Genomic Unit guidelines. Microarray images were analyzed using GenePix 6.0 software (Molecular Devices). Data were transformed using an intensity-based Lowess function (Yang et al., 2002) with Acuity 4.0 software (Axon Instruments). Genes were considered as differentially expressed if they fulfilled both a FDR < 0.05 after a SAM test analysis (Tusher et al., 2001), and the fold-change (up or down) was above 1.75-fold.
To assign a putative biological function to every detected differentially expressed gene, their respective orthologous genes from the wild species F. vesca, which genome has been recently released (Shulaev et al., 2011), were identified by blasting the EST sequence associated with each singular spot within the array to the overall collection of F. vesca predicted genes (Altschul et al., 1990; Shulaev et al., 2011; http://www.rosaceae.org/). In order to enrich this process, A. thaliana putative orthologs were also identified for every F × ananassa gene as vast functional information is available for the former species (TAIR10: http://www.arabidopsis.org/). FunCat and GO terms assignments were first used to perform an automated functional categorization of differentially expressed genes (Ashburner et al., 2000; Ruepp et al., 2004). FatiGO tool (a web tool to find significant associations of Gene Ontology terms with gene groups; Ashburner et al., 2000; Ruepp et al., 2004) was used to perform a comprehensive Singular Enrichment Analysis (SEA) to extract relevant GO terms associated with up-regulated genes. Briefly stated, it takes two lists of genes (ideally a group of interest and the remaining genes in the experiment, while any two groups, formed in any way, can be tested against each other) and converts them into two lists of GO annotations using the corresponding gene or protein—term annotation table. A Fisher's exact test for 2 × 2 contingency tables is subsequently used to check for significant over-representation of GO annotations in one of the sets with respect to the other. Multiple test correction to account for the multiple hypotheses tested (one for each functional term) is applied. The terms are considered to be relevant by the application of statistical tests, as described in Al-Shahrour et al. (2004).
Hormone Determination in Strawberry Tissues
Extraction and purification procedures and chromatographic analysis have been previously described (Durgbanshi et al., 2005). In short, 3 grams of frozen green tissue were lyophilized and immediately homogenized in 5 mL of ultrapure water. After centrifugation (5000 g, 10 min), supernatant pH was adjusted to 2.8 with 15% (v/v) CH3COOH, and the supernatant was partitioned twice against an equal volume of diethyl ether. The aqueous phase was discarded and the organic fraction was vacuum evaporated at room temperature. The solid residue was resuspended in 1 mL of a 90:10 (v/v) water/methanol solution and subsequently filtered through a cellulose acetate filter (0.22 μm). A 20 μL aliquot of this solution was then injected into the high performance liquid chromatography (HPLC) system from Waters, Milford MA (Alliance 2690 system). Aliquots were injected on a Nucleosil ODS reversed-phase column. Phytohormones were eluted with a gradient of methanol and 0.01% CH3COOH in water that started from 10:90 (v/v) and linearly reached 60:40 (v/v) in 10 min. In the following 4 min, the gradient was increased to 80:20 (v/v). Isocratic conditions of 80:20 (v/v) were then retained during the last 2 min of the run. Initial conditions were restored and allowed to equilibrate for 5 min, for a total time of 21 min per sample. The solvent flow rate was 0.3 mL/min, with working pressures at around 70–100 bar.
Quoted plant hormone endogenous contents are mean values from 2 measurements of each of all 3 biological replicates. The One-way Analysis of Variance (ANOVA) with a Bonferroni Multiple Comparisons Test was performed using GraphPad InStat3 for Windows (GraphPad Software, La Jolla California USA, www.graphpad.com) to calculate the significant differences between control and inoculated plants.
Monitoring C. acutatum Development in Plant Tissue
The histopathology of the interaction strawberry-C. acutatum has been previously reported (Curry et al., 2002), and the transition to the necrotrophic phase was established after 4 dpi (Horowitz et al., 2002). However, fungal colonization in plant may vary with infection conditions. Thus, we monitored the infection progress of C. acutatum on strawberry plant after artificial inoculation (spraying) using microscopic analysis.
To identify the initial infection stages, conidial development was monitored. The number of ungerminated conidia, germinated conidia without forming appresoria, and germinated conidia with fully developed appressoria, was monitored 1, 3, 5, 7 and 9 dpi (days post-inoculation; Figure 1). No visual disease symptoms were detected in the plant during this period of time. At 5 dpi 44.76% of conidia registered fully developed appresoria (Figure 1A). Other conidia did not germinate (34.98%); others, while they did, they did not produce appresoria (20.24%). Early mycelium formation was observed microscopically at 5 dpi (Figure 1D), whereas more extended and abundant mycelium was detected at 7 and 9 dpi, respectively (Figures 1E,F).
Figure 1. Colletotrichum acutatum development in strawberry tissue. (A) Progression of C. acutatum isolate CECT 20240 conidial germination. (B–F) Photomicrographs of C. acutatum structures formed during infection of strawberry leaves over a 9-days period. (B) Non germinated conidia at 1 dpi; (C) Germinated conidia after 3 dpi; (D,E) Appressorium and mycelium development at 5 and 7 dpi, respectively; (F) Abundant mycelium at 9 dpi. All bars 25 μm.
As pathogen development in crown it is difficult to monitor by optical microscopy in the same way as aerial tissues, the progress of pathogen development in crown was monitored using TEM with the inoculation methods presented in Arroyo et al. (2005) (Figure 2). Such higher concentration of inoculum expedited infection development. Figure 2A shows penetration peg formation through the host cuticle and a small infection vesicle, 36 h post-inoculation (hpi), which reflects the establishment of the biotrophic stage of C. acutatum. No morphological signs of cuticular component degradation were observed during this penetration phase. Necrotic signals appeared as early as 4 dpi (Figure 2B), and pathogen development clearly expanded throughout vascular tissue at 7 dpi (Figure 2C).
Figure 2. Transmission Electronic Microscopy (TEM) in strawberry crown. (A) Penetration of cuticle and biotrophic stage of C. acutatum by subcuticular hyphae on wall epidermal cell 36 h post-inoculation. (B) Necrotrophic signals of development of C. acutatum by intercellular hypha 4 days post-inoculation. (C) Necrotrophic development of C. acutatum by intracellular hyphae in xylem tissues at 7 dpi. A, appresorium; C, plant cuticle; IV, infection vesicle; W, wall plant cell; IM, intramural o intercellular hypha; FW, fungal cell wall; CD, cell debris; IC, intracellular hyphae; X, xylem. Bar = 1 μm.
Oxidative Stress Analyses
ROS production was monitored by diaminobenzidine (DAB) staining throughout the plant infection time-course. No clear accumulation of DAB (reddish-brown color visible to the naked-eye) was evident at any stage of the infection analyzed here (data not shown). Furthermore, under light microscopy, DAB staining beneath fungal appressoria and surrounding the penetration pegs was very rarely observed and only after 7 dpi (Figures 3A–C). FaGST1 transcripts (the strawberry AtGST1 orthologous gene) did not increase across the time points examined, neither in crown (Figure 3D) nor in leaves (Figure 3E). A slight increase in the expression of this gene was detected in petiole tissue 3–5 dpi (Figure 3F).
Figure 3. Light microscopic analysis of DAB stained preparations of C. acutatum infected strawberry tissue at 7 dpi, and relative expression values by RTqPCR analysis of FaGST1 (At1g02930 ortholog) gene at different time points of infection. (A) Tissue sample showing no DAB reaction, (B,C) very rarely detected reduction of DAB beneath fungal appressoria, and surrounding the penetration pegs, at upper and lower focal planes, respectively. All bars 25 μm. (D–F) Relative FaGST1 expression values. At each time point, every inoculated sample was compared with its corresponding mock treated sample. In the graphics, standard value 1 at T0 was added to better illustrate changes.
Based on these results, and the observation that spray inoculation with 104 conidia·ml−1 more closely mimics natural infection conditions, the 5 dpi time point was chosen for transcriptome analysis. The assumption was that most of the changes in transcript accumulation in response to early C. acutatum infection would be detected around this time point, including those in response to early stages of switching to necrotrophic growth. In addition, analysis of known defense-related transcripts was performed at 1, 3, 5, 7, and 9 dpi to assess plant response to a wider range of fungal developmental stages and disease progression.
Expression Profiling of C. acutatum-Infected Strawberry
Transcript abundance was measured in crown tissue 5 days after C. acutatum infection and in mock-treated plants. Tables 1A,B and Tables S1, S2 show a summary of identified transcripts with the highest induction or repression levels after crown infection. A total of 147 transcripts varied in accumulation more than 1.75-fold following criteria described under Material and Methods. Of these, 118 genes were induced, and 29 genes were repressed. Verification of gene expression changes by real-time RTqPCR was performed on eleven up-regulated genes and two down-regulated genes, representing different categories shown in Tables 1A,B (marked with #). The expression pattern of analyzed genes after C. acutatum inoculation was consistent with that obtained by microarray analysis.
Table 1A. Up-regulated genes by Colletotrichum acutatum in crown tissue of Fragaria x ananassa, cultivar Camarosa.
Table 1B. Down-regulated genes by Colletotrichum acutatum in crown tissue of Fragaria x ananassa, cultivar Camarosa.
Automated functional analysis by FunCat and GO terms assignments showed that many of these up- and down-regulated genes described in Tables 1A,B belong to plant defense and stress response-related categories (Figure S1 and Tables S3–S5). However, when no obvious functional role was annotated within the corresponding F. vesca ortholog genes, a thorough search through the references available in the database from many plant species was performed (see Table S3). A wider range of the strawberry transcripts matched defense and biotic stress annotations with the number of up- and down-regulated genes changing 89.93 and 65.51%, respectively.
Transcripts representing five subsets of putative molecular function were identified. These include plant receptors, signal transduction mechanisms under hormonal control (protein modification and degradation), transcriptional changes (transcription factors), new protein synthesis, and secretion of active defense components (PR proteins, degradative enzymes or chemical defenses).
Singular Enrichment Analysis (SEA) Shows Key Components of SA-Mediated Signaling Pathway are Up-regulated
A comprehensive Singular Enrichment Analysis (SEA) using FatiGO (Al-Shahrour et al., 2004) identified key processes altered in strawberry after C. acutatum infection (Figure 4, and Table S6). Transcripts increasing in abundance include (p < 0.005) those associated with Systemic Acquired Resistance and SA-mediated signaling pathways, responding to bacterium and fungus, and activating the immune response. Strawberry genes within these enriched categories are genes FaEDS1 (EDS1-936EST, AtEDS1-like) and FaPAD4 (M4E10EST, AtPAD4-like), which are known to be involved in PRR- and R-mediated pathogen-induced SA accumulation in other plants; genes FaWRKY70-1 and FaWRKY70-2 (M17H1EST, and M12E12EST, respectively, two AtWRKY70-like genes); gene FaMeSA1 (M9D5EST, a methyl salicylate esterase); gene FaPBS1 (M13C5EST, a SA-dependent Ser/Thr kinase); and gene FaGRX1 (M30F8EST, similar to a member of the glutaredoxin family that regulates the protein redox state), which are major downstream components of the SA signal transduction pathway, and known to be activators of SA-dependent defense in many plants.
Figure 4. Singular Enrichment Analysis (SEA) on strawberry up-regulated genes. FatiGO was used to extract relevant GO terms for biological processes (BP) at p < 0.005. The terms are considered to be relevant by the application of statistical tests as described in Al-Shahrour et al. (2004). Data are presented as a heat map, as prompts color intensity correlates with adj. p-value, the highest intensity, the lowest adj. p-value. See Table S1 for a detailed list of further genes belonging to overrepresented functions.
RTqPCR Analysis Indicates that Components of Jasmonic Acid Defense Signaling Pathway are Also Induced in Strawberry after C. acutatum Infection
Analysis of a representative set of up-regulated genes, including those used for microarray validation (see Tables 1A,B), was carried out by RTqPCR after crown infection (Figure 5). All genes tested showed clear upregulation after infection (Figure 5A). Similar results were also found on petiole tissue analysis after infection (Figure S2A). In addition, two main expression patterns were detected. Transcripts corresponding to FaWRKY1 (J_4-9EST), FaLIP-1 (M12C12EST), FaCHI4-2 (M16D12EST), FaPR5-2 (EPR5-77EST), and FaPR10-4 (M22A10EST) reached a maximum level at 5 dpi, while transcripts corresponding to FaWRKY2 (M21B3EST), FaLRR1 (ELRR-39EST), FaGLN-2 (M24B7EST), and FaPR5-1 (EPR5-284EST) reached their maximum expression level at later times. Transcripts representing FaWRKY1 (J_4-9EST) and FaWRKY2 (M21B3EST), which belong to the same family of transcription factors (orthologs to Arabidopsis AtWRKY75), registered a different timing in their response. Other variations in timing were also detected for the three members of the PR5 family (FaPR5-1, FaPR5-2, FaPR5-3).
Figure 5. Relative expression values by RTqPCR analysis of upregulated strawberry genes during C. acutatum infection. (A) Relevant strawberry genes in this study, (B) JA-responsive marker genes, and (C) SA-responsive marker genes. Strawberry crowns were harvested at different days post-treatment (dpi) either with mock or C. acutatum spore suspension. At each time point, every inoculated sample was compared with its corresponding mock treated sample. In the graphics, standard value 1 at T0 was added to better illustrate changes. Asterisk marks genes not present in the Array dataset. Arabidopsis orthologous genes are AT5G13080 (AtWRKY75), AT5G42650 (AtAOS), AT2G38470 (AtWRKY33), AT3G45140 (AtLOX2), AT2G46370 (AtJAR1), AT5G44420 (AtPDF1.2), AT3G48090 (AtEDS1), AT3G52430 (AtPAD4), AT1G28480 (AtGRX480), AT3G56400 (AtWRKY70), AT2G14610 (AtPR1), AT3G57260 (AtPR2).
An additional experiment was performed using in-vitro plantlets to describe networks associated with SA or JA treatment. Plantlets were used because they show enhanced sensitivity and a faster response to hormone treatment against mature plants (Figure 6). Almost all C. acutatum induced transcripts also showed significant induction after MeJA, but not so after SA treatment (Figure 6A). Transcripts corresponding to PR genes such as FaGLN-2 (M24B7EST), FaCHI4-2 (M16D12EST), FaPR5-2 (EPR5-77EST), FaPR5-1 (EPR5-284EST), FaPR5-3 (M1F10EST), and FaPR10-4 (M22A10EST), as well as the WRKY75-like transcription factors [genes FaWRKY1 (J_4-9EST) and FaWRKY2 (M21B3EST)] increased in abundance in response to JA.
Figure 6. Relative expression values by RTqPCR analysis of upregulated strawberry genes by hormone treatment. (A) Relevant strawberry genes in this study, (B) JA-responsive marker genes, and (C) SA-responsive marker genes. Strawberry plantlets were treated with mock, SA (5 mM) and MeJA (2 mM) elicitors, and harvested at different hours post-treatment (htp). At each time point, every elicited sample was compared against its corresponding mock treated sample. Left and right legends represent relative expression values for SA and JA treatments, respectively. In the graphics, standard value 1 at T0 was added to better illustrate changes. Asterisk marks genes not present in the Array dataset. Arabidopsis orthologous are AT5G13080 (AtWRKY75), AT5G42650 (AtAOS), AT2G38470 (AtWRKY33), AT3G45140 (AtLOX2), AT2G46370 (AtJAR1), AT5G44420 (AtPDF1.2), AT3G48090 (AtEDS1), AT3G52430 (AtPAD4), AT1G28480 (AtGRX480), AT3G56400 (AtWRKY70), AT2G14610 (AtPR1), AT3G57260 (AtPR2).
Incomplete Activation of SA and JA Pathways Occurs during C. acutatum Infection
To further investigate how SA- and JA- hormone-dependent pathways are responding to C. acutatum infection, transcript levels of their described signaling pathways were measured by RTqPCR. Transcripts corresponding to JA-associated genes FaWRKY33-1 (M8H3EST) and FaWRKY33-2 (M1C12EST) (two orthologs to AtWRKY33), FaAOS-1 (M28A2EST, AtAOS ortholog), FaLOX2-1 (AtLOX2 ortholog), FaJAR1 (AtJAR1 ortholog), and FaPDF1 (AtPDF1.2 ortholog), and SA-associated transcripts such as FaEDS1 (EDS1-936EST, AtEDS1 ortholog), FaPAD4 (M4E10EST, AtPAD4 ortholog), FaGRX1 (M30F8EST, AtGRX480 ortholog), FaWRKY70-1 and FaWRKY70-2 (M17H1EST, and M12E12EST, respectively), FaPR1-1 and FaPR1-2 (AtPR1 ortholog), FaPR2-1 and FaPR2-2 (AtPR2 ortholog), were analyzed in crown tissue in response to C. acutatum inoculation (Figures 5B,C), and after MeJA or SA exogenous applications (Figures 6B,C).
None of known JA inducible pathways, except for gene regulators FaWRKY33-1, FaWRKY33-2, and FaAOS-1 (whose induction was also detected by microarray, see Tables 1A,B) were infection induced (Figure 5B). Similar results were also obtained on petiole tissue analysis after infection (Figure S2B). On the contrary, all transcripts were induced (FaAOS-1, FaWRKY33-1, FaWRKY33-2, FaLOX2-1, FaJAR1, and FaPDF1) in response to MeJA treatment (Figure 6B).
On the other hand all SA-pathway associated orthologs were induced in crown after C. acutatum infection (Figure 5C) except for FaPR1-1, a well-known SA-pathway-associated gene. The same result was found in petiole (Figure S2C). Nearly all strawberry SA-associated transcripts including FaPR1-1, increased in abundance in response to SA treatment. Two classical SA-associated PR orthologous genes, FaPR1-2 and FaPR2-2, are shown to be mainly JA-dependent in strawberry (Figure 6C).
Level of SA and JA during the Strawberry/C. acutatum Interaction
Salicylic acid (SA) and jasmonic acid (JA) content in strawberry plantlets cv. Camarosa was measured at 3 and 5 days after inoculation with C. acutatum. A significant 272% increase in free SA concentration was detected in aerial tissues only at 3 dpi (202.21 ng/g dw, in infected plantlets vs. 74.42 ng/g dw, in mock-treated plantlets; Figure 7A). In addition, free SA in infected plantlets increased up to 678% at 5 dpi compared with that of mock treatment (354.77 ng/g dw vs. 52.33 ng/g dw, respectively). A significant 241% increase in free JA concentration was detected in infected plantlets at 3 dpi compared with that of mock treatment (771.39 ng/g dw, vs. 320.42 ng/g dw, respectively; Figure 7B). This increase was even higher after 5 dpi (425%), compared with that of mock treatment (1707.03 ng/g dw vs. 401.93 ng/g dw, respectively).
Figure 7. Quantification of (A) SA and (B) JA [ng g-1 (dry weigh)] in aerial tissues of in-vitro strawberry plantlets at 3 and 5 days after mock (white bars) and pathogen inoculation (black bars). Data are the mean of three biological replicates, and error bars represent the SD. Numbers inside arrows indicate percentage of increase of infected samples vs. mock. Asterisks indicate confidence of ANOVA-Bonferroni Multiple Comparisons Test (***p < 0.001; **p < 0.01).
This report characterizes the molecular events related to strawberry response to C. acutatum infection largely prior to plant disease symptom occurrence. In other pathosystems, main changes in gene expression are known to occur before significant lesion development (Windram et al., 2012). Previous research indicates that C. acutatum may persist and grow extensively on strawberry tissues without causing visible symptoms (Freeman et al., 1998; Leandro et al., 2001; Debode et al., 2009). Under inoculation conditions (104 conidia·ml−1) used herein, no symptoms were visible during the time course of the experiment. Conditions mirror those used in previously published work (Casado-Díaz et al., 2006; Encinas-Villarejo et al., 2009; Amil-Ruiz et al., 2012, 2013). Transmission electron microscopy results (Figure 2) are roughly consistent with the timing of pathogen development observed after spray inoculation, taking into account that spot inoculation with a higher spore concentration (106 conidia·ml−1) is forcing a much quicker infection process. Indeed, this inoculum concentration is also used by Horowitz et al. (2002) to study the developmental stages of this pathogen in strawberry, and necrotrophic transition was first identified after 4 dpi, which is consistent with our TEM observations. In addition, the penetration phase of C. acutatum on strawberry tissues is not a synchronous process, as previously described for this pathogen and other Colletotrichum spp. on other hosts (Makowski and Mortensen, 1998; Arroyo et al., 2005). Furthermore, at 5 days after inoculation, less than fifty percent of the conidia showed fully formed appresoria, and limited primary hyphae formation was observed (Figure 1). To strengthen the transcriptomic analysis, a more detailed gene expression study (RTqPCR) was performed throughout the C. acutatum infection stages from 1 to 9 days after inoculation, which certainly should reflect changes in host in response to most pathogen development stages. A similar time-window for gene expression has been monitored in other pathosystems (Doehlemann et al., 2008; Marcel et al., 2010; Vargas et al., 2012). Previous studies on C. acutatum-strawberry interaction have also shown that main transcript changes occur in host along the period of time herein analyzed (Casado-Díaz et al., 2006; Encinas-Villarejo et al., 2009). All in all, these observations informed about the appropriate time points for transcriptomic analyses.
Molecular Components of the Strawberry Response to C. acutatum Identified in this Study
Following C. acutatum infection, the transcriptomic analysis shows a wide range of responses to this pathogen. Many of the identified transcripts encode proteins with demonstrated roles in resistance and defense functions. Other transcript classes identified include members of plant pathogen perception and sensing apparatus, signal transduction machinery, transcriptional factors and regulatory genes, and protein synthesis and secretion mechanisms. For the purposes of this report, only components belonging to known SA and JA signaling pathways are discussed.
Signaling Transduction Pathways: Downstream Responses against C. acutatum
In model plants, one of the big gaps in plant immunity understanding lies in the signaling pathways operating immediately downstream of PRR and R protein activation. However, partially understood pathways have been established (Dodds and Rathjen, 2010). We have found that members of kinase, phosphatase, ubiquitin and calcium-associated gene families related to signal transduction pathways were induced upon interaction with C. acutatum (Tables 1A,B and Table S3). Known components of both SA- and JA-dependent defense signaling pathways were also up-regulated. No induction was observed of important components of MAP Kinases pathway, which seems to play a role in the interplay between SA and JA-defense signaling pathways (Rodriguez et al., 2010; Caarls et al., 2015). The reasons behind this are unclear, but it is possible that the pathogen could be manipulating part of this pathway.
Enrichment of transcripts corresponding to specific SA-pathway members was detected. Thus, the expression of genes FaEDS1 (EDS1-936EST) and FaPAD4 (M4E10EST) is induced by C. acutatum. The lipase-like protein EDS1 is an important node acting upstream of SA in PAMP-triggered immunity (PTI) after PRRs stimulation and is also required for signaling of all TIR-NB-LRRs tested to date (Wiermer et al., 2005; Heidrich et al., 2011). These findings suggest that specific effector-triggered immunity (ETI) through TIR domain signaling might also be acting in strawberry against this pathogen. EDS1 is known to physically interact with two other positive regulators, PAD4 and SAG101, both of which are putative lipases, although hydrolase activity has not been demonstrated for either protein (Wiermer et al., 2005). Interestingly, the expression of Arabidopsis EDS1 is positively regulated by WRKY70 transcription factors (Li et al., 2004), and enrichment in WRKY70 orthologs has also been detected in strawberry (see further below). Moreover, a strawberry PAD4 ortholog (FaPAD4) was indeed upregulated. PAD4 affect SA accumulation (Wang et al., 2011). Thus, the dissociated forms of EDS1 and PAD4 are fully competent in signaling receptor triggered localized cell death at infection locations (Rustérucci et al., 2001; Aviv et al., 2002) but, by contrast, an EDS1–PAD4 complex is necessary for basal resistance involving transcriptional up-regulation of PAD4 itself and mobilization of salicylic acid defenses (Rietz et al., 2011).
Local SA production has been shown to trigger defenses in the surrounding cells downstream of EDS1 and PAD4 activity. This SAR is activated thorough a systemic signal that primes distal tissues against similar invaders. The SA derivative methyl salicylate (MeSA) is believed to serve as a long-distance phloem-mobile SAR signal in plants (Liu et al., 2011; Dempsey and Klessig, 2012). Once in the distal, uninfected tissue, MeSA must be converted into biologically active SA by esterase activity (Dempsey and Klessig, 2012). Interestingly, induction of the strawberry gene M9D5EST encoding a methyl salicylate (MeSA) esterase similar to the Arabidopsis AtMES9 was detected (Tables 1A,B), which suggests that this signaling mechanism might also be activated in strawberry during C. acutatum interaction. Curiously, the Arabidopsis AtMES9 presents in-vitro activity with MeSA, MeJA and MeIAA (Yang et al., 2008) but it showed preference for MeSA as a substrate (Vlot et al., 2008; Dempsey and Klessig, 2012).
Induction of other transcripts acting downstream of SA was detected. Two WRKY70-like genes, FaWRKY70-1 (M17H1EST) and FaWRKY70-2 (M12E12EST), and a glutaredoxin GRX480-like gene, FaGRX1 (M30F8EST), which have been described as essential components for SA-dependent defense activation, were detected in strawberry after infection. In addition, transcript accumulation of strawberry orthologs to classical SA associated genes such as PRs FaPR1-2, FaPR2-1, and FaPR2-2 was highly induced in crown or petiole (Figure 5C and Figure S2C). Recently the FaPR2-2 transcript has been reported as a reliable indicator of SA-dependent defenses in strawberry, as it was induced by C. acutatum, C. fragariae, and SA (Zamora et al., 2012). However, although FaPR1-2, FaPR2-1 and FaPR2-2 were detected in this study after SA treatment, both FaPR1-2 and FaPR2-2 turned out to be also induced after JA treatment (Figure 6C). Therefore, we propose that these two later genes should not be considered as SA-specific transcripts in strawberry.
C. acutatum infection also increased M8G7EST transcript, encoding a protein that resembles the HSF-like transcription factor TBF1, a member of an extensive family of heat responsive proteins (Baniwal et al., 2004; Ikeda and Ohme-Takagi, 2009). These proteins are associated with diverse functions, including heat stress response (Charng et al., 2007; Ikeda et al., 2011), and plant development (Pernas et al., 2010; ten Hove et al., 2010). Interestingly, the TBF1 protein has recently been shown to be a major molecular switch for plant growth-to-defense transition in Arabidopsis (Pajerowska-Mukhtar et al., 2012). This transcription factor is a positive regulator of immune responses induced by salicylic acid and PAMPs, and it binds to the TL1 (GAAGAAGAA) cis element of NPR1-dependent ER-resident genes required for antimicrobial protein secretion.
Transcripts related to the JA-mediated signaling pathway were also induced in strawberry after C. acutatum infection, including FaAOS-1 (M28A2EST), FaWRKY33-1 (M8H3EST), and FaWRKY33-2 (M1C12EST). FaAOS-1 encodes an allene oxide synthase, a member of the cytochrome p450 CYP74 gene family (Song et al., 1993) that functions as a key enzyme in the initial steps of the JA biosynthetic pathway (Peña-Cortés et al., 2004; Leon-Reyes et al., 2010), thus generating signaling molecules that are essential for host immunity and plant development (Acosta and Farmer, 2010; Gfeller et al., 2010; Bak et al., 2011). Interestingly, while only a single copy of AOS gene exists in Arabidopsis (Kubigsteltig et al., 1999), a small AOS gene family with five members can be detected in F. vesca genome (unpublished), and three AOS members have been detected in tomato (López-Ráez et al., 2010), suggesting a more complex regulation of this pathway in fruiting plants. In addition, FaWRKY33-1 and FaWRKY33-2 are duplicated in Fragaria and can be paralogous, and both are similar to the well-known WRKY33 transcription factor from Arabidopsis. This important transcription factor acts downstream of JA and regulates expression of classical JA-dependent defense genes, such as those encoding glucanases, chitinases, and thaumatin-like proteins, which have been extensively described as JA-associated genes in other plants. Accordingly, many strawberry orthologs to these JA-induced proteins, such as FaGLN-2, FaCHI4-2, FaPR10-4, FaPR5-1, FaPR5-2, and FaPR5-3, were strongly induced by C. acutatum (Figure 5A).
Upregulation of genetic components needed for SA and JA synthesis in strawberry is accompanied by a concomitant increase in concentration of such phytohormones in response to inoculation by C. acutatum (Figure 7). Therefore, taken together, these results clearly demonstrate that both SA and JA defense signaling pathways are activated in strawberry during C. acutatum infection.
In addition, main transcript accumulation was observed between 3 and 5 dpi for genes FaWRKY1, FaCHI4-2, FaPR10-4, FaPR5-2, FaPR5-3, FaLIP1, FaPR2-1, and FaPR2-2, coding for known PR defense proteins (Figure 5). This increase roughly coincides with the break time observed during conidial germination and appressoria formation of C. acutatum (Figure 1), and it may reflect an early attempt of strawberry plant to halt pathogen development. After 5 dpi, expression of these genes decreases gradually while a progressive increase in the percentage of germinated conidia with fully developed appressoria occurs (Figure 1). The latter may reflect the process whereby pathogen hijacks plant defenses to its own benefit.
Evidence of C. acutatum Influence on SA- and JA-Dependent Defense Pathways to Promote Pathogen Development in Strawberry
The SA and JA signaling pathways are activated in strawberry challenged with C. acutatum. However, results indicate that expression of well-known components of both SA- and JA-dependent defense pathways are not activated during this interaction. FaPR1-1 has been described as SA inducible in strawberry in cv. Pájaro, challenged with the avirulent strain M23 of C. fragariae but not after infection with virulent strain M11 of C. acutatum (Grellet-Bournonville et al., 2012). Infection with the avirulent strain M23 induced temporary SA accumulation (nearly 2-fold) in strawberry plants that was accompanied by induction of FaPR1-1 transcripts and protection against a subsequent infection with C. acutatum.
Here the infection of cv. Camarosa with C. acutatum induced the SA accumulation to a higher level (nearly 3-fold), but no significant accumulation of FaPR1-1 transcript above basal levels in crown and petiole was observed (Figure 5C and Figure S2C). In contrast, many other SA-responsive PRs (FaPR1-2, FaPR2-1, and FaPR2-1) were upregulated. The FaPR1-1 transcript did increase in abundance after SA treatment in “Camarosa,” which means that strawberry cv. Camarosa has the ability to induce this gene (Figure 6C).
No significant increase in FaGST1 expression (the strawberry AtGST1 heterolog, a well-known oxidative stress inducible transcript) was observed in crown and leaf tissue (Figures 3D,E) at any time after pathogen inoculation. Only a slight increase was observed in petioles after 5 dpi. These results suggest that little, if any, ROS production occurs during pathogen development. No significant DAB staining was evident in strawberry tissue during this time course, which is consistent with the results. It is known that ROS is generated in response to abiotic stress, particularly H2O2, which is an active signaling molecule triggering a variety of defense responses, with GST induction being the most significant (Bhattacharjee, 2012).
Production of H2O2 and activation of gene PR-1 gene occur simultaneously, and it has been suggested that H2O2 acts downstream from SA in the pathogenesis-related (PR-1) gene induction . Intriguingly, no ROS production and no significant FaGST and FaPR1.1 transcript induction were detected in strawberry after pathogen inoculation, even though SA increased. Our results strongly correlate with those of Grellet-Bournonville et al. (2012), who also reported that virulent C. acutatum strain M11 showed no ROS accumulation and no FaPR1.1 gene induction after strawberry infection.
On the other hand, expression of strawberry orthologs to JA-associated defense related genes, such as FaPDF1, FaLOX2-1, and FaJAR1, also remained unchanged after infection with C. acutatum, despite the fact that many other components of the JA-mediated signaling pathway were induced (FaAOS-1, FaWRKY33-1, FaWRKY33-2; Figure 5B and Figure S2B).
Our results indicate that a number of known components of both SA- and JA-dependent defense pathways are not activated in strawberry during interaction with C. acutatum. Furthermore, our results support a hypothesis that C. acutatum might be interfering with certain branches of both SA- and JA-dependent defense pathways in strawberry. Absence of a significant plant defense response would allow successful colonization of host tissue by this pathogen. In this sense, recent results reported on the tomato-Botrytis system have shown that the exopolysaccharide production by this pathogen [EPS, known as b-(1,3)(1,6)-D-glucan] acted as elicitor of the tomato SA biosynthesis pathway and that inappropriate induction of SA by this pathogen impaired tomato JA-dependent defenses by interrupting the JA signaling pathway downstream of JA production (El Oirdi et al., 2011; Rahman et al., 2012). Consequently, the fungus could gradually spread through tomato plant tissues.
It is tempting to speculate that C. acutatum may be able to interact with strawberry defense response. Based on our results, and considering the canonical SA JA crosstalk model proposed in Arabidopsis (Spoel et al., 2007; Spoel and Dong, 2012), an integrated model has been devised (Figure 8). In this model, in which activation of both SA and JA pathways and increased amount of SA and JA signal molecules occurs after C. acutatum infection, negative crosstalk between these two pathways should somehow be expected. Spoel et al. (2007) showed that simultaneous A. thaliana inoculation with a biotrophic and a necrotrophic pathogen resulted in impaired resistance to the necrotrophic pathogen, and demonstrated that the SA pathway that was activated by the biotrophus suppressed the level of JA-dependent resistance against the necrotroph. Indeed, SA-mediated suppression of JA-responsive gene expression has been reported to be targeted downstream of JA biosynthesis (Leon-Reyes et al., 2010). Thus, GRX480 is an NPR1 dependent-SA-inducible class III glutaredoxin (Rouhier et al., 2006; Krinke et al., 2007) specific to land plants (Ziemann et al., 2009), which interacts with TGA factors and suppresses JA-responsive PDF1.2 transcription (Ndamukong et al., 2007; Zander et al., 2011). In addition, NPR1 and WRKY70 also acts downstream of the SA molecule as a node of convergence for JA-mediated and SA-mediated signals (Dong, 2004; Li et al., 2004; Wang et al., 2006; Ren et al., 2008), balancing the JA- and SA-dependent responses (Li et al., 2006). Interestingly, strawberry orthologs FaGRX1, FaWRKY70-1 and FaWRKY70-2 were specifically induced during interaction with C. acutatum, but no induction on FaPDF1 and other major components of JA-dependent signaling pathway, such as genes FaLOX2-1 and FaJAR1, was detected. Moreover, an increase in JA synthesis and upregulation of FaAOS-1, the strawberry ortholog to Arabidopsis AtAOS, a well-known JA-associated marker gene encoding a key enzyme for JA synthesis, was also reported after C. acutatum infection, supporting the idea that in strawberry repression of some JA-responsive genes is targeted downstream of the JA biosynthesis. Indeed, AtAOS has been described as a MeJA-inducible gene but not suppressed by WRKY70 (Li et al., 2006). Very interestingly, the expression of a second group of known JA-responsive genes, such as FaGLN-2, FaCHI4-2, FaPR10-4, FaPR5-1, FaPR5-2, FaPR5-3, increased after being challenged with this pathogen, suggesting the presence in strawberry of a second GRX480/WRKY70-independent JA-dependent defense branch.
Figure 8. Hypothetical model of SA- and JA-dependent defense pathways activated in strawberry in response to C. acutatum. This model is based on the canonical pathways described in model plant. Upon interaction with C. acutatum, the strawberry plant activates upstream components of SA and JA defense pathways. Thus, synthesis of these signal molecules increases and main downstream key components for SA (FaNPR1, FaWRKY70, FaGRX) and JA (FaWRKY33) are activated. Unlike Arabidopsis WRKY33, FaWRKY33 does not act as a negative regulator of the entire SA-dependent defense signaling pathway, either by a direct or an indirect effect of fungal activity, but only for some components (FaGST1, FaPR1.1). That allows FaGRX, together with FaNPR1 and FaWRKY70, to act as negative regulators of JA responsive genes similarly to their Arabidopsis orthologs. As a result, important JA-responsive defense marker genes, such as FaLOX2, FaJAR1, and FaPDF1, are not induced. These impaired mechanisms might provide some advantage for fungal spreading.
Importantly, two other JA-associated AtWRKY33-like genes, FaWRKY33-1 and FaWRKY33-2, were also upregulated in strawberry by C. acutatum. The JA-associated component AtWRKY33 has recently been reported as a key transcriptional regulator of defense responses to necrotrophus (Birkenbihl et al., 2012). Indeed, AtWRKY33 acts as a negative regulator of the SA-defense pathway upon pathogen infection and negatively controls expression of many important genes, including those responsible for SA biosynthesis and accumulation, positive regulatory proteins EDS1 and PAD4, and SA responsive genes PR1, PR2, and PR3. Interestingly, in strawberry, expression of the SA-dependent orthologous genes FaGST and FaPR1-1 remained unaltered but very intriguingly, the synthesis of SA and the expression of orthologs to components of SA-mediated signaling pathway acting upstream (FaEDS1 and FaPAD4), and downstream of SA (FaGRX1, FaWRKY70-1, FaWRKY70-2, FaPR1-2, FaPR2-1, and FaPR2-2), was remarkably induced during the infection with C. acutatum, despite the fact that FaWRKY33-1 and FaWRKY33-2 were clearly upregulated. Unlike what has been previously described for AtWRKY33 (Birkenbihl et al., 2012), these results suggest that repressive control of many known components of the SA-pathway through these FaWRKY33 transcription factors does not work in strawberry during interaction with C. acutatum (Figure 8).
In conclusion, our results demonstrate both that SA and JA increase in strawberry upon C. acutatum infection, and that known plant defenses through SA and JA dependent signaling pathways are partially promoted during such interaction. Indeed, major plant defense marker genes are not up-regulated in strawberry after infection with this pathogen, which might evidence a putative molecular strategy used by this pathogen to overcome strawberry plant defense. This is in line with the new emerging paradigm that a key pathogen virulence strategy involves modulation of plant hormone signaling (El Oirdi et al., 2011; Chung et al., 2013).
Results from our research are of great use to further our understanding of the strawberry immune system to enable future disease control through biotechnological and breeding strategies.
FAR and supervisors JC and JMB conceived and designed the experiments. FAR, JC, and JMB contributed to reagents/materials and analysis tools. FAR, JGG, JG, RB, AM, OT, BD, FR, JAM, FP carried out experiments, analyzed and interpreted the data. FAR, JC contributed to drafting the manuscript. AA, FA, JGG, FAR, JC performed the light and TEM microscopic observations. JC, JMB, JG revising the article critically. All authors approved the final version of the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Authors are grateful to Dr. JM Lopez-Aranda (IFAPA-Centro de Churriana) for providing micropropagated strawberry plants and to Nicolás García-Caparrós for technical assistance. Authors also want to thank Kevin M. Folta for his insightful comments on the paper. This work was supported by Junta de Andalucía, Spain [Proyectos de Excelencia P07-AGR-02482/P12-AGR-2174, and grants to Grupo-BIO278].
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.01036
EPS, exopolysaccharide; ER-resident, endoplasmic reticulumresident; ET, ethylene; ETI, effector-triggered immunity; FDR, false discovery rate; JA, jasmonic acid; PAMP, pathogen-associated molecular pattern; PR proteins, pathogenesis related proteins; PRR, pattern recognition receptor; PTI, PAMPtriggered immunity; qPCR, quantitative PCR; R protein, resistance protein; SA, salicylic acid; SEA, singular enrichment analysis
Al-Shahrour, F., Díaz-Uriarte, R., and Dopazo, J. (2004). FatiGO: a web tool for finding significant associations of Gene Ontology terms with groups of genes. Bioinformatics 20, 578–580. doi: 10.1093/bioinformatics/btg455
Amil-Ruiz, F., Blanco-Portales, R., Muñoz-Blanco, J., and Caballero, J. L. (2011). The strawberry plant defence mechanism: a molecular review. Plant Cell Physiol. 52, 1873–1903. doi: 10.1093/pcp/pcr136
Amil-Ruiz, F., Encinas-Villarejo, S., de los Santos, B., Muñoz-Mérida, A., Mercado, J. A., Trelles, O., et al. (2012). “Distinctive transcriptome response of two strawberry (Fragaria x ananassa) cultivars to Colletotrichum acutatum infection,” in XXVIII International Horticultural Congress on Science and Horticulture for People (IHC2010): International Symposium on Genomics and Genetic Transformation of Horticultural Crops, eds R. E. Litz, K. Folta, M. Talon, and F. Pliego Alfaro (Lisbon: ISHS), 47–50.
Amil-Ruiz, F., Garrido-Gala, J., Blanco-Portales, R., Folta, K. M., Muñoz-Blanco, J., and Caballero, J. L. (2013). Identification and validation of reference genes for transcript normalization in strawberry (Fragaria × ananassa) defense responses. PLoS ONE 8:e70603. doi: 10.1371/journal.pone.0070603
Arroyo, F. T., Moreno, J., García-Herdugo, G., de los Santos, B., Barrau, C., Porras, M., et al. (2005). Ultrastructure of the early stages of Colletotrichum acutatum infection of strawberry tissues. Can. J. Bot. 83, 491–500. doi: 10.1139/b05-022
Aviv, D. H., Rustérucci, C., Iii, B. F. H., Dietrich, R. A., Parker, J. E., and Dangl, J. L. (2002). Runaway cell death, but not basal disease resistance, in lsd1 is SA- and NIM1/NPR1-dependent. Plant J. 29, 381–391. doi: 10.1046/j.0960-7412.2001.01225.x
Baniwal, S. K., Bharti, K., Chan, K. Y., Fauth, M., Ganguli, A., Kotak, S., et al. (2004). Heat stress response in plants: a complex game with chaperones and more than twenty heat stress transcription factors. J. Biosci. 29, 471–487. doi: 10.1007/BF02712120
Birkenbihl, R. P., Diezel, C., and Somssich, I. E. (2012). Arabidopsis WRKY33 Is a key transcriptional regulator of hormonal and metabolic responses toward botrytis cinerea infection. Plant Physiol. 159, 266–285. doi: 10.1104/pp.111.192641
Casado-Díaz, A., Encinas-Villarejo, S., de los Santos, B., Schilirò, E., Yubero-Serrano, E.-M., Amil-Ruiz, F., et al. (2006). Analysis of strawberry genes differentially expressed in response to Colletotrichum infection. Physiol. Plant. 128, 633–650. doi: 10.1111/j.1399-3054.2006.00798.x
Charng, Y. Y., Liu, H. C., Liu, N. Y., Chi, W. T., Wang, C. N., Chang, S. H., et al. (2007). A heat-inducible transcription factor, HsfA2, is required for extension of acquired thermotolerance in Arabidopsis. Plant Physiol. 143, 251–262. doi: 10.1104/pp.106.091322
Chung, S. H., Rosa, C., Scully, E. D., Peiffer, M., Tooker, J. F., Hoover, K., et al. (2013). Herbivore exploits orally secreted bacteria to suppress plant defenses. Proc. Natl. Acad. Sci. U.S.A. 110, 15728–15733. doi: 10.1073/pnas.1308867110
Curry, K. J., Abril, M., Avant, J. B., and Smith, B. J. (2002). Strawberry anthracnose: histopathology of Colletotrichum acutatum and C. fragariae. Phytopathology 92, 1055–1063. doi: 10.1094/PHYTO.2002.92.10.1055
Debode, J., Van Hemelrijck, W., Baeyen, S., Creemers, P., Heungens, K., and Maes, M. (2009). Quantitative detection and monitoring of Colletotrichum acutatum in strawberry leaves using real-time PCR. Plant Pathol. 58, 504–514. doi: 10.1111/j.1365-3059.2008.01987.x
Doehlemann, G., Wahl, R., Horst, R. J., Voll, L. M., Usadel, B., Poree, F., et al. (2008). Reprogramming a maize plant: transcriptional and metabolic changes induced by the fungal biotroph Ustilago maydis. Plant J. 56, 181–195. doi: 10.1111/j.1365-313X.2008.03590.x
Durgbanshi, A., Arbona, V., Pozo, O., Miersch, O., Sancho, J. V., and Gomez-Cadenas, A. (2005). Simultaneous determination of multiple phytohormones in plant extracts by liquid chromatography−electrospray tandem mass spectrometry. J. Agric. Food Chem. 53, 8437–8442. doi: 10.1021/jf050884b
El Oirdi, M., El Rahman, T. A., Rigano, L., El Hadrami, A., Rodriguez, M. C., Daayf, F., et al. (2011). Botrytis cinerea manipulates the antagonistic effects between immune pathways to promote disease development in tomato. Plant Cell Online 23, 2405–2421. doi: 10.1105/tpc.111.083394
Encinas-Villarejo, S., Maldonado, A. M., Amil-Ruiz, F., de los Santos, B., Romero, F., Pliego-Alfaro, F., et al. (2009). Evidence for a positive regulatory role of strawberry (Fragaria x ananassa) Fa WRKY1 and Arabidopsis At WRKY75 proteins in resistance. J. Exp. Bot. 60, 3043–3065. doi: 10.1093/jxb/erp152
Freeman, S., Horowitz, S., and Sharon, A. (2001). Pathogenic and nonpathogenic lifestyles in Colletotrichum acutatum from strawberry and other plants. Phytopathology 91, 986–992. doi: 10.1094/PHYTO.2001.91.10.986
Freeman, S., Katan, T., and Shabi, E. (1998). Characterization of Colletotrichum species responsible for anthracnose diseases of various fruits. Plant Dis. 82, 596–605. doi: 10.1094/PDIS.19220.127.116.116
Grellet-Bournonville, C. F., Martinez-Zamora, M. G., Castagnaro, A. P., and Díaz-Ricci, J. C. (2012). Temporal accumulation of salicylic acid activates the defense response against Colletotrichum in strawberry. Plant Physiol. Biochem. 54, 10–16. doi: 10.1016/j.plaphy.2012.01.019
Guidarelli, M., Carbone, F., Mourgues, F., Perrotta, G., Rosati, C., Bertolini, P., et al. (2011). Colletotrichum acutatum interactions with unripe and ripe strawberry fruits and differential responses at histological and transcriptional levels. Plant Pathol. 60, 685–697. doi: 10.1111/j.1365-3059.2010.02423.x
Heidrich, K., Wirthmueller, L., Tasset, C., Pouzet, C., Deslandes, L., and Parker, J. E. (2011). Arabidopsis EDS1 connects pathogen effector recognition to cell compartment–specific immune responses. Science 334, 1401–1404. doi: 10.1126/science.1211641
Horowitz, S., Freeman, S., and Sharon, A. (2002). Use of green fluorescent protein-transgenic strains to study pathogenic and nonpathogenic lifestyles in Colletotrichum acutatum. Phytopathology 92, 743–749. doi: 10.1094/PHYTO.2002.92.7.743
Ikeda, M., Mitsuda, N., and Ohme-Takagi, M. (2011). Arabidopsis HsfB1 and HsfB2b act as repressors of the expression of heat-inducible hsfs but positively regulate the acquired thermotolerance. Plant Physiol. 157, 1243–1254. doi: 10.1104/pp.111.179036
Khan, A. A., and Shih, D. S. (2004). Molecular cloning, characterization, and expression analysis of two class II chitinase genes from the strawberry plant. Plant Sci. 166, 753–762. doi: 10.1016/j.plantsci.2003.11.015
Krinke, O., Ruelland, E., Valentová, O., Vergnolle, C., Renou, J.-P., Taconnat, L., et al. (2007). Phosphatidylinositol 4-Kinase activation is an early response to salicylic acid in Arabidopsis suspension cells. Plant Physiol. 144, 1347–1359. doi: 10.1104/pp.107.100842
Leandro, L. F. S., Gleason, M. L., Nutter, F. W., Wegulo, S. N., and Dixon, P. M. (2001). Germination and sporulation of Colletotrichum acutatum on symptomless strawberry leaves. Phytopathology 91, 659–664. doi: 10.1094/PHYTO.2001.91.7.659
Leon-Reyes, A., Van der Does, D., De Lange, E., Delker, C., Wasternack, C., Van Wees, S., et al. (2010). Salicylate-mediated suppression of jasmonate-responsive gene expression in Arabidopsis is targeted downstream of the jasmonate biosynthesis pathway. Planta 232, 1423–1432. doi: 10.1007/s00425-010-1265-z
Li, J., Brader, G., and Palva, E. T. (2004). The WRKY70 transcription factor: a node of convergence for Jasmonate-mediated and Salicylate-mediated signals in plant defense. Plant Cell 16, 319–331. doi: 10.1105/tpc.016980
Liu, P.-P., von Dahl, C. C., Park, S.-W., and Klessig, D. F. (2011). Interconnection between methyl salicylate and lipid-based long-distance signaling during the development of systemic acquired resistance in Arabidopsis and tobacco. Plant Physiol. 155, 1762–1768. doi: 10.1104/pp.110.171694
López-Ráez, J. A., Verhage, A., Fernández, I., García, J. M., Azcón-Aguilar, C., Flors, V., et al. (2010). Hormonal and transcriptional profiles highlight common and differential host responses to arbuscular mycorrhizal fungi and the regulation of the oxylipin pathway. J. Exp. Bot. 61, 2589–2601. doi: 10.1093/jxb/erq089
Makowski, R. M. D., and Mortensen, K. (1998). Latent infections and penetration of the bioherbicide agent Colletotrichum gloeosporioides f. sp. malvae in non-target field crops under controlled environmental conditions. Mycological Res. 102, 1545–1552. doi: 10.1017/S0953756298006960
Maleck, K., Levine, A., Eulgem, T., Morgan, A., Schmid, J., Lawton, K. A., et al. (2000). The transcriptome of Arabidopsis thaliana during systemic acquired resistance. Nat. Genet. 26, 403–410. doi: 10.1038/82521
Marcel, S., Sawers, R., Oakeley, E., Angliker, H., and Paszkowski, U. (2010). Tissue-adapted invasion strategies of the rice blast fungus Magnaporthe oryzae. Plant Cell 22, 3177–3187. doi: 10.1105/tpc.110.078048
Ndamukong, I., Abdallat, A. A., Thurow, C., Fode, B., Zander, M., Weigel, R., et al. (2007). SA-inducible Arabidopsis glutaredoxin interacts with TGA factors and suppresses JA-responsive PDF1.2 transcription. Plant J. 50, 128–139. doi: 10.1111/j.1365-313X.2007.03039.x
Pajerowska-Mukhtar, K. M., Wang, W., Tada, Y., Oka, N., Tucker, C. L., Fonseca, J. P., et al. (2012). The HSF-like transcription factor TBF1 is a major molecular switch for plant growth-to-defense transition. Curr. Biol. 22, 103–112. doi: 10.1016/j.cub.2011.12.015
Peña-Cortés, H., Barrios, P., Dorta, F., Polanco, V., Sánchez, C., Sánchez, E., et al. (2004). Involvement of jasmonic acid and derivatives in plant response to pathogen and insects and in fruit ripening. J. Plant Growth Regul. 23, 246–260. doi: 10.1007/s00344-004-0035-1
Rahman, T. A. E., Oirdi, M. E., Gonzalez-Lamothe, R., and Bouarab, K. (2012). Necrotrophic pathogens use the salicylic acid signaling pathway to promote disease development in tomato. Mol. Plant Microbe Interact. 25, 1584–1593. doi: 10.1094/MPMI-07-12-0187-R
Ren, C.-M., Zhu, Q., Gao, B.-D., Ke, S.-Y., Yu, W.-C., Xie, D.-X., et al. (2008). Transcription factor WRKY70 displays important but no indispensable roles in jasmonate and salicylic acid signaling. J. Integr. Plant Biol. 50, 630–637. doi: 10.1111/j.1744-7909.2008.00653.x
Rietz, S., Stamm, A., Malonek, S., Wagner, S., Becker, D., Medina-Escobar, N., et al. (2011). Different roles of Enhanced Disease Susceptibility1 (EDS1) bound to and dissociated from Phytoalexin Deficient4 (PAD4) in Arabidopsis immunity. New Phytol. 191, 107–119. doi: 10.1111/j.1469-8137.2011.03675.x
Robert-Seilaniantz, A., Grant, M., and Jones, J. D. G. (2011). Hormone crosstalk in plant disease and defense: more than just JASMONATE-SALICYLATE antagonism. Annu. Rev. Phytopathol. 49, 317–343. doi: 10.1146/annurev-phyto-073009-114447
Ruepp, A., Zollner, A., Maier, D., Albermann, K., Hani, J., Mokrejs, M., et al. (2004). The FunCat, a functional annotation scheme for systematic classification of proteins from whole genomes. Nucleic Acids Res. 32, 5539–5545. doi: 10.1093/nar/gkh894
Rustérucci, C., Aviv, D. H., Holt, B. F., Dangl, J. L., and Parker, J. E. (2001). The disease resistance signaling components eds1 and pad4 are essential regulators of the cell death pathway controlled by LSD1 in Arabidopsis. Plant Cell Online 13, 2211–2224. doi: 10.1105/tpc.13.10.2211
Sarowar, S., Zhao, Y., Soria-Guerra, R. E., Ali, S., Zheng, D., Wang, D., et al. (2011). Expression profiles of differentially regulated genes during the early stages of apple flower infection with Erwinia amylovora. J. Exp. Bot. 62, 4851–4861. doi: 10.1093/jxb/err147
Sasaki, Y., Asamizu, E., Shibata, D., Nakamura, Y., Kaneko, T., Awai, K., et al. (2001). Monitoring of methyl jasmonate-responsive Genes in Arabidopsis by cDNA macroarray: Self-activation of jasmonic acid biosynthesis and crosstalk with other phytohormone signaling pathways. DNA Res. 8, 153–161. doi: 10.1093/dnares/8.4.153
Schenk, P. M., Kazan, K., Manners, J. M., Anderson, J. P., Simpson, R. S., Wilson, I. W., et al. (2003). Systemic gene expression in Arabidopsis during an Incompatible Interaction with Alternaria brassicicola. Plant Physiol. 132, 999–1010. doi: 10.1104/pp.103.021683
Shulaev, V., Sargent, D. J., Crowhurst, R. N., Mockler, T. C., Folkerts, O., Delcher, A. L., et al. (2011). The genome of woodland strawberry (Fragaria vesca). Nat. Genet. 43, 109–116. doi: 10.1038/ng.740
Song, W. C., Funk, C. D., and Brash, A. R. (1993). Molecular cloning of an allene oxide synthase: a cytochrome P450 specialized for the metabolism of fatty acid hydroperoxides. Proc. Natl. Acad. Sci. U.S.A. 90, 8519–8523. doi: 10.1073/pnas.90.18.8519
Spoel, S. H., Johnson, J. S., and Dong, X. (2007). Regulation of tradeoffs between plant defenses against pathogens with different lifestyles. Proc. Natl. Acad. Sci. U.S.A. 104, 18842–18847. doi: 10.1073/pnas.0708139104
Staswick, P. E., and Tiryaki, I. (2004). The oxylipin signal jasmonic acid is activated by an enzyme that conjugates it to isoleucine in Arabidopsis. Plant Cell Online 16, 2117–2127. doi: 10.1105/tpc.104.023549
ten Hove, C. A., Willemsen, V., de Vries, W. J., van Dijken, A., Scheres, B., and Heidstra, R. (2010). SCHIZORIZA encodes a nuclear factor regulating asymmetry of stem cell divisions in the Arabidopsis root. Curr. Biol. 20, 452–457. doi: 10.1016/j.cub.2010.01.018
Tusher, V. G., Tibshirani, R., and Chu, G. (2001). Significance analysis of microarrays applied to the ionizing radiation response. Proc. Natl. Acad. Sci. U.S.A. 98, 5116–5121. doi: 10.1073/pnas.091062498
Vargas, W. A., Martín, J. M. S., Rech, G. E., Rivera, L. P., Benito, E. P., Díaz-Mínguez, J. M., et al. (2012). Plant defense mechanisms are activated during biotrophic and necrotrophic development of colletotricum graminicola in maize. Plant Physiol. 158, 1342–1358. doi: 10.1104/pp.111.190397
Venugopal, S. C., Jeong, R.-D., Mandal, M. K., Zhu, S., Chandra-Shekara, A. C., Xia, Y., et al. (2009). Enhanced disease susceptibility 1 and salicylic acid act redundantly to regulate resistance gene-mediated signaling. PLoS Genet. 5:e1000545. doi: 10.1371/journal.pgen.1000545
Vlot, A. C., Liu, P.-P., Cameron, R. K., Park, S.-W., Yang, Y., Kumar, D., et al. (2008). Identification of likely orthologs of tobacco salicylic acid-binding protein 2 and their role in systemic acquired resistance in Arabidopsis thaliana. Plant J. 56, 445–456. doi: 10.1111/j.1365-313X.2008.03618.x
Wang, D., Amornsiripanitch, N., and Dong, X. (2006). A genomic approach to identify regulatory nodes in the transcriptional network of systemic acquired resistance in plants. PLoS Pathog. 2:e123. doi: 10.1371/journal.ppat.0020123
Wang, G.-F., Seabolt, S., Hamdoun, S., Ng, G., Park, J., and Lu, H. (2011). Multiple Roles of WIN3 in regulating disease resistance, cell death, and flowering time in Arabidopsis. Plant Physiol. 156, 1508–1519. doi: 10.1104/pp.111.176776
Windram, O., Madhou, P., McHattie, S., Hill, C., Hickman, R., Cooke, E., et al. (2012). Arabidopsis defense against botrytis cinerea: chronology and regulation deciphered by high-resolution temporal transcriptomic analysis. Plant Cell Online 24, 3530–3557. doi: 10.1105/tpc.112.102046
Yang, Y. H., Dudoit, S., Luu, P., Lin, D. M., Peng, V., Ngai, J., et al. (2002). Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res. 30:e15. doi: 10.1093/nar/30.4.e15
Yang, Y., Xu, R., Ma, C. J., Vlot, A. C., Klessig, D. F., and Pichersky, E. (2008). Inactive methyl Indole-3-Acetic acid ester can be hydrolyzed and activated by several esterases belonging to the AtMES esterase family of Arabidopsis. Plant Physiol. 147, 1034–1045. doi: 10.1104/pp.108.118224
Zamora, M. G. M., Bournonville, C. G., Castagnaro, A. P., and Ricci, J. C. D. (2012). Identification and characterisation of a novel class I endo-β-1,3-glucanase regulated by salicylic acid, ethylene and fungal pathogens in strawberry. Funct. Plant Biol. 39, 412–420. doi: 10.1071/FP11275
Zander, M., Chen, S., Imkampe, J., Thurow, C., and Gatz, C. (2011). Repression of the Arabidopsis thaliana jasmonic acid/ethylene-induced defense pathway by TGA-Interacting Glutaredoxins depends on their C-Terminal ALWL Motif. Mol. Plant. 5, 831–840. doi: 10.1093/mp/ssr113
Keywords: Colletotrichum acutatum, Fragaria × ananassa, quantification of gene expression, salicylic and jasmonic acid, strawberry defense response
Citation: Amil-Ruiz F, Garrido-Gala J, Gadea J, Blanco-Portales R, Muñoz-Mérida A, Trelles O, de los Santos B, Arroyo FT, Aguado-Puig A, Romero F, Mercado J-Á, Pliego-Alfaro F, Muñoz-Blanco J and Caballero JL (2016) Partial Activation of SA- and JA-Defensive Pathways in Strawberry upon Colletotrichum acutatum Interaction. Front. Plant Sci. 7:1036. doi: 10.3389/fpls.2016.01036
Received: 26 April 2016; Accepted: 01 July 2016;
Published: 15 July 2016.
Edited by:Abdul Latif Khan, University of Nizwa, Oman
Reviewed by:Lei Zhang, North Carolina State University, USA
Marcos Egea-Cortines, Universidad Politécnica Cartagena, Spain
Copyright © 2016 Amil-Ruiz, Garrido-Gala, Gadea, Blanco-Portales, Muñoz-Mérida, Trelles, de los Santos, Arroyo, Aguado-Puig, Romero, Mercado, Pliego-Alfaro, Muñoz-Blanco and Caballero. 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) or licensor 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: José L. Caballero, firstname.lastname@example.org