ORIGINAL RESEARCH article

Front. Plant Sci., 15 July 2016

Sec. Plant Pathogen Interactions

Volume 7 - 2016 | https://doi.org/10.3389/fpls.2016.01036

Partial Activation of SA- and JA-Defensive Pathways in Strawberry upon Colletotrichum acutatum Interaction

  • 1. Departamento de Bioquímica y Biología Molecular e Instituto Andaluz de Biotecnología, Edificio Severo Ochoa (C6), Universidad de Córdoba Córdoba, Spain

  • 2. Instituto de Biología Molecular y Celular de Plantas, Universidad Politécnica de Valencia–Consejo Superior de Investigaciones Científicas Valencia, Spain

  • 3. Departamento de Arquitectura de Computación, Universidad de Málaga, Campus de Teatinos Málaga, Spain

  • 4. Centro Andalusian Institute of Agricultural and Fishering Research and Training (IFAPA) Las Torres-Tomejil, CAPMA–Junta de Andalucía Sevilla, Spain

  • 5. Departamento de Biología Vegetal, Facultad de Ciencias, Universidad de Málaga, Campus de Teatinos Málaga, Spain

Abstract

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.

Introduction

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.

Results

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

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

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

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

Fxananassa gene IDF. vesca orthologA. thaliana orthologGene descriptionRelation with defense/biological functionCC vs. CI
Fold change (RTqPCR)FDR q-value
INVASION SENSING
M13C5*gene07245AT5G13160Serine/threonine-protein kinase PBS1Receptor kinase, R protein-guard model4.170
M19F7*gene15497AT4G33210SLOMO (SLOw MOtion) F-box/LRR-repeat proteinFbox/LRR protein, plant receptor, Proteasome complex2.657.93E-03
M2F10*gene19270AT4G00340Receptor-like protein kinase 4Receptor kinase, Signal transduction regulation2.437.93E-03
M14D5gene13911AT1G30240Proline-, glutamic acid- and leucine-rich protein 1LRR protein, plant receptor2.357.93E-03
M6C2gene01890AT5G42090Lung seven transmembrane receptor family proteinPlant receptor2.357.93E-03
ELRR-39#gene25524AT5G21090CPR30 Leucine-rich repeat (LRR) family proteinLRR protein, plant receptor2.14 (2.82 ± 0.54)7.93E-03
M29F3gene16731AT3G14460LRR and NB-ARC domains-containing disease resistance proteinCC-NBS-LRR class of R proteins, plant receptor2.007.93E-03
M18E3gene20858AT3G14460Proline-, glutamic acid- and leucine-rich protein 1CC-NBS-LRR class of R proteins, plant receptor1.847.93E-03
SIGNAL TRANSDUCTION
M23A9gene14522AT4G35790Phospholipase D deltaPhospholipase D, Transduction of stress responses8.260
M27D3gene18784AT5G01160RING/U-box E3 ubiquitin-protein ligaseE3 ligase, Proteasome complex7.440
M16B7gene00744AT1G69960Serine/threonine-protein phosphatase PP2A catalytic subunitSer/Thr protein phosphatase 2A catalytic subunit5.200
M19D11gene10418AT3G03940Casein Serine/threonine-protein kinaseCalcium binding kinase4.590
M13C5*gene07245AT5G13160Serine/threonine-protein kinase PBS1Receptor kinase, R protein-guard model4.170
EDS1-936*#gene09503AT3G48090EDS1-specific diacylglycerol lipase alphaLipase, SA pathway regulator3.82 (4.4 ± 1.13)0
M23A6gene32391AT4G11740Ubiquitin-like superfamily proteinUbiquitin, Proteasome complex3.497.60E-03
M27C10gene30942.3utrAT5G25510Serine/threonine protein phosphatase 2A regulatory B subunitSer/Thr protein phosphatase 2A catalytic subunit3.440
M8G2gene10067AT4G30960CIPK-Serine/threonine-protein kinase 6Calcium binding kinase SOS23.180
M4F10gene21532AT1G65430E3 ubiquitin-protein ligase ARI8E3 ligase, Proteasome complex3.150
M8G7gene24036AT4G36990TBF1 Heat shock factor proteinMajor molecular switch for plant growth-to-defense transition3.067.93E-03
M24D7*gene28350AT5G40150Peroxidase superfamily proteinClass III peroxidase2.767.93E-03
M10E2gene02575AT1G27460NPGR1-No pollen germination related 1Calmoduling binding protein2.737.93E-03
M3D5gene23778AT1G05180NEDD8-activating enzyme E1 regulatory subunit / AXR1 (Auxin resistant 1) NAD(P)-binding Rossmann-fold superfamily proteinFbox, JA signaling, Proteasome complex2.737.93E-03
M19F7*gene15497AT4G33210SLOMO (SLOw MOtion) F-box/LRR-repeat proteinFbox/LRR protein, plant receptor, Proteasome complex2.657.93E-03
M25E7gene01516AT1G15780Bromodomain-containing proteinInteract with calcium binding protein kinase2.477.93E-03
M8D11*gene06214AT1G60490Phosphatidylinositol 3-kinaseProtein kinase, Protein trafficking, Secretory Pathway2.467.93E-03
M13H9gene12681AT5G57020Myristoyl-CoA:protein N-myristoyltransferaseCo-traslational addition of myristic acid2.457.93E-03
M2F10*gene19270AT4G00340Receptor-like protein kinase 4Receptor kinase, Signal transduction regulation2.437.93E-03
M7G11gene04753AT1G69640Sphingoid base hydroxylase 1 (SBH1)Synthesis of membrane components2.427.93E-03
M4E10*#gene16110AT3G52430Phytoalexin deficient 4, LipaseLipase, Chemical defenses, SA pathway regulator2.33 (2.67 ± 0.64)7.93E-03
M4C3gene15015AT5G10930CIPK-Serine/threonine-protein kinase 5Calcium binding kinase2.257.93E-03
M14H1gene07894AT3G51860Vacuolar cation/proton exchanger 3Proton/Calcium antiporter2.207.93E-03
M7B6gene05859AT1G80210BRCC36A - homologous recombinationHomologous recombination, deubiquitinating activity, Proteasome complex2.207.93E-03
M21H5gene01441AT5G56180Actin-related protein 8Fbox/Actin/helicase domain, proteasome complex, XXXRNAmetabolism2.057.93E-03
M4E6gene12959AT4G332401-phosphatidylinositol-4-phosphate 5-kinaseProtein kinase, protein trafficking, endomembrane homeostasis2.047.93E-03
M28C8*gene12445AT1G05260Peroxidase superfamily proteinClass III peroxidase1.987.93E-03
M17E3gene06367AT4G24830Argininosuccinate synthaseNO synthesis, signal transduction1.987.93E-03
M3E6*#gene27591AT1G71695Peroxidase superfamily protein (Prx12)Class III peroxidase1.92 (5.5 ± 0.7)7.93E-03
M10B6gene01594AT3G13460YTH domain family protein 2Calcium transport to nucleus, regulate gene expression1.867.93E-03
M13F3gene28416AT3G27925Protease DegP1Protease1.797.93E-03
M1H8gene12874AT5G53360E3 Ubiquitin protein ligase SINAT3E3 ligase, proteasome complex1.759.42E-03
NEW PROTEIN SYNTHESIS AND SECRETION
M21B3#gene01340AT5G13080WRKY DNA-binding protein 75Transcription factor5.79 (19.47 ± 5.2)0
M8H8gene10702AT4G17960ATP-dependent RNA helicase DBP10RNA metabolism5.610
M26G7gene31909AT2G25970RNA binding KH domain-containing proteinRNA metabolism5.350
M22D9gene22758AT3G51980Armadillo repeat superfamily protein-Hsp70 nucleotide exchange factor fes1Protein folding4.800
J_4-9#gene07210AT5G13080WRKY DNA-binding protein 75Transcription factor3.89 (56.68 ± 7.93)0
M11C6gene03828AT1G6962060S Ribosomal protein L34Protein synthesis3.790
M6G7gene32154AT3G48030Hypoxia-responsive Zinc finger (C3HC4-type RING finger) family proteinTranscription factor3.750
M10C12gene08531AT1G75780Tubulin beta-1 chainCitosqueleton3.620
M9F6gene29752AT1G28420Homeobox protein orthopediaTranscription factor3.440
M1A2gene24354AT1G62020Coatomer subunit alphaProtein transport3.310
M23C4gene02623AT4G37750AINTEGUMENTA gene - AP2 like transcription factorTranscription factor3.207.60E-03
M18A9gene30367AT5G46190RNA-binding KH domain-containing proteinRNA metabolism2.857.93E-03
M7G4gene23202AT3G52250Duplicated homeodomain-like superfamily proteinRNA metabolism2.757.93E-03
M23C7gene25539.3utrAT4G3386540S ribosomal protein S29Protein synthesis2.577.93E-03
M17H1*#gene13547AT3G56400WRKY DNA-binding protein 70Transcription factor, SA-JA crosstalk2.53 (4.99 ± 0.54)7.93E-03
M18F1gene09051AT1G47490RNA-binding protein 47CRNA metabolism2.497.93E-03
M8D11*gene06214AT1G60490Phosphatidylinositol 3-kinaseProtein kinase, Protein trafficking, Secretory Pathway2.467.93E-03
M11H4gene22626AT3G12110Actin 11Citosqueleton2.427.93E-03
M8H3*#gene13803AT2G38470WRKY DNA-binding protein 33Transcription factor, JA pathway2.41 (3.58 ± 1.52)7.93E-03
M14B5gene29081AT1G59740Peptide transporter PTRProtein secretion2.397.93E-03
M5B8gene24582AT5G22950Vacuolar protein sorting-associated protein 24Protein secretion2.227.93E-03
M12E12*#gene21365AT3G56400WRKY DNA-binding protein 70Transcription factor, SA-JA crosstalk2.19 (2.96 ± 0.54)7.93E-03
M3A1gene30880AT3G16060Kinesin-related proteinCitosqueleton2.127.93E-03
M19E4gene05323AT2G44710RNA-binding (RRM/RBD/RNP motifs) family proteinRNA metabolism2.127.93E-03
M18C5gene04135AT1G66140Zinc finger protein 4Transcription factor2.087.93E-03
M3E11gene25805AT1G18650Plasmodesmata callose-binding endo-1,3-beta-glucosidase protein 3 (PdCB3)Cell-to-cell trafficking2.027.93E-03
M12B6no hit foundAT3G25940DNA-directed RNA polymerase TFIIB zinc-binding proteinRNA metabolism2.017.93E-03
M7D1gene10625AT3G0559060S ribosomal protein L18-2Protein synthesis2.007.93E-03
M20A3gene21473AT5G16715Valyl-tRNA synthetaseProtein synthesis1.987.93E-03
M8A6gene00998AT1G77030DEAD-box ATP-dependent RNA helicase 29RNA metabolism1.937.93E-03
M9E2gene15731AT1G80070Pre-mRNA-processing-splicing factor SUS2RNA metabolism1.927.93E-03
M28B7gene16235.5utrAT2G22430Homeobox-leucine zipper protein ATHB-6Transcription factor1.897.93E-03
M1C12*#gene28174AT2G38470WRKY DNA-binding protein 33Transcription factor, JA pathway1.86 (9.10 ± 1)7.93E-03
M6A9gene00185AT5G67300Transcription factor MYB44Transcription factor1.839.42E-03
M4C6gene20572AT3G62310RNA helicase family proteinRNA metabolism1.799.42E-03
DIRECT DEFENSES
M24B7#gene14817AT4G16260Glycosyl hydrolase superfamily proteinCell wall degradation, PR protein family47.54 (30.54 ± 16.25)0
M16D12#gene02717AT3G54420Chitinase class IVPR protein family7.93 (116.84 ± 22.54)0
EPR5-77#gene32423AT4G11650Pathogenesis-related 5 family proteinPR protein family7.52 (59.11 ± 10.05)0
M5B6gene24296.3utrAT5G09360LaccaseLignin biosynthesis7.480
M23A10gene07086AT1G24020Pathogenesis-related 10 family proteinPR protein family7.080
M12C12#gene31975AT5G14180Triacylglycerol lipase 2Lipase, Chemical defenses6.60 (19.59 ± 2.7)0
M6G11gene26351AT4G34135Flavonol 7-O-glucosyltransferaseSecondary metabolism4.340
M6B9gene05185AT1G24020Pathogenesis-related 10 family proteinPR protein family3.890
EPR5-284#gene32422AT4G11650Pathogenesis-related 5 family proteinPR protein family3.88 (5.63 ± 1.93)0
M1F10#gene09812AT1G20030Pathogenesis-related 5 family proteinPR protein family3.69 (17.58 ± 5.7)0
M22A10#gene07085AT1G24020Pathogenesis-related 10 family proteinPR protein family3.20 (8.87 ± 2.38)0
M24D7*gene28350AT5G40150Peroxidase superfamily proteinClass III peroxidase2.767.93E-03
M5G8gene07082AT1G24020Pathogenesis-related 10 family proteinPR protein family2.677.93E-03
M10C5gene00687AT1G24020Pathogenesis-related 10 family proteinPR protein family2.667.93E-03
M26E5gene32023AT5G17000Zinc-binding dehydrogenase family protein / oxidoreductaseRedox protection2.657.93E-03
M4F3gene27555AT1G22750D-serine/D-alanine/glycine transporterSecondary metabolism2.657.93E-03
M25D10gene07087AT1G24020Pathogenesis-related 10 family proteinPR protein family2.447.93E-03
M5C8gene11632AT4G32320L-ascorbate peroxidase 6Antioxidant defenses2.367.93E-03
M4E10*gene16110AT3G52430Phytoalexin deficient 4, LipaseLipase, Chemical defenses, SA pathway regulator2.33 (2.67 ± 0.64)7.93E-03
M23D11gene20700AT4G37990Cinnamyl alcohol dehydrogenaseLignin biosynthesis2.147.93E-03
M29A9gene21697AT3G54420Endochitinase PR4PR protein family2.027.93E-03
M25D11gene17437AT3G07320O-Glycosyl hydrolases family 17, (1->3)-beta-glucanaseCell wall degradation, PR protein family1.987.93E-03
M28C8*gene12445AT1G05260Peroxidase superfamily proteinClass III peroxidase1.987.93E-03
M3E6*#gene27591AT1G71695Peroxidase superfamily protein (Prx12)Class III peroxidase1.92 (5.5 ± 0.7)7.93E-03
M10D7gene07065AT1G24020Fra a 2 allergenPR protein family1.867.93E-03
M26G2gene31048AT2G30370CHAL secreted proteinInhibite stomatal production1.797.93E-03
M21G5gene04724AT1G69530Expansin-A1Stomatal movement1.769.42E-03
HORMONE-DEPENDENT PATHWAYS
EDS1-936*#gene09503AT3G48090EDS1-specific diacylglycerol lipase alphaLipase, SA pathway regulator3.82 (4.4 ± 1.13)0
M12E4gene32179AT1G27500Tetratricopeptide repeat (TPR)-like superfamily proteinTetratricopeptide repeat3.320
M22A6gene05545AT1G80360Pyridoxal phosphate (PLP)-dependent transferases superfamily proteinPyridoxal-phosphate, oxidative stress response2.847.93E-03
M14G2gene31738AT4G39820Tetratricopeptide repeat (TPR)-like superfamily proteinTetratricopeptide repeat2.697.93E-03
M8H2gene09899AT5G642502-nitropropane dioxygenaseJA pathway2.677.93E-03
M26D3gene18908AT4G01100Adenine nucleotide transporter 1 (ADNT1)Purine transporter, Signaling2.567.93E-03
M17H1*#gene13547AT3G56400WRKY DNA-binding protein 70Transcription factor, SA-JA crosstalk2.53 (4.99 ± 0.54)7.93E-03
M25B1gene23034AT3G13790Cell wall Invertase 1 (AtcwINV1): Glycosyl hydrolases family 32 proteinCell wall invertase, signaling2.487.93E-03
M9E10gene03078AT1G44750Purine permease 11Purine transporter, Signaling2.447.93E-03
M8H3*#gene13803AT2G38470WRKY DNA-binding protein 33Transcription factor, JA pathway2.41 (3.58 ± 1.52)7.93E-03
M4E10*#gene16110AT3G52430Phytoalexin deficient 4, LipaseLipase, Chemical defenses, SA pathway regulator2.33 (2.67 ± 0.64)7.93E-03
M12E12*#gene21365AT3G56400WRKY DNA-binding protein 70Transcription factor, SA-JA crosstalk2.19 (2.96 ± 0.54)7.93E-03
M23C11gene08617AT1G76180Dehydrin cold-regulated 47ABA responsive2.157.93E-03
M16H1gene14094.3utrno hit foundAuxin response factorAuxin responsive2.147.93E-03
M9D5gene29393AT4G37150Methyl salicylate (MeSA) esterase 9SA release from MeSA2.037.93E-03
M30F8#gene29769.3utrAT1G28480Glutaredoxin GRX480SA pathway, REDOX signaling1.92 (4.52 ± 0.82)7.93E-03
M1C12*#gene28174AT2G38470WRKY DNA-binding protein 33Transcription factor, JA pathway1.86 (9.10 ± 1)7.93E-03
M28A2#gene15063AT5G42650Allene oxide synthaseJA synthesis1.75 (1.5 ± 0.64)9.42E-03
NO OBVIOUSLY RELATED TO DEFENSE RESPONSE
M22B1gene01044AT2G25660Embryo defective 24108.250
M18E11gene27435AT1G34550Embryo defective 27566.550
M21E9gene24023AT2G24960MRG family protein, chromatin binding3.310
M7B12gene07388AT2G21170Triosephosphate isomerase2.967.60E-03
M24C11gene32086AT1G64385Unknown protein, endomembrane system2.887.60E-03
M13A4gene23331AT5G13520Aminopeptidase M1 family protein / Leukotriene A-4 hydrolase2.397.93E-03
M27A2gene13677AT1G32060Phosphoribulokinase2.087.93E-03
M4E4gene05017AT5G49930Embryo defective 14412.077.93E-03
M25G5gene06563.3utrAT4G13930Serine hydroxymethyltransferase 42.047.93E-03
M3F5gene13777AT3G08890Protein of unknown function1.977.93E-03
M22G7gene09933.3utrAT5G41835non-LTR retrotransposon family1.933.64E-02
M4F8gene15022.3utrAT2G25140Casein lytic proteinase B4/heat shock protein1.917.93E-03

Up-regulated genes by Colletotrichum acutatum in crown tissue of Fragaria x ananassa, cultivar Camarosa.

Genes were considered as differentially expressed if they fulfilled a FDR < 0.05 after a SAM test analysis and the fold-change was higher that 1.75-fold between the compared conditions. Fold change values represent the ratio of cv. Camarosa mock (CC) vs. infected (CI). (#) indicates genes further analyzed by real time RTqPCR to validate microarray result. Their relative expression value at 5dpi is shown as (media ± SD) in the “Fold Change” column. (

) marks no obvious detection of F. vesca ortolog gene due to putative fails by automated gene prediction (see also Table S3). 3utr, and 5utr, indicate F x ananassa sequences matching untranslated regions of the corresponding F. vesca gene (see also Table S3). Regulated genes are grouped in sections accordingly to their role in different steps of the defense response against C. acutatum (see Table S4 for associated references). Asterisk marks genes which take part in more than one unique functional group.

Table 1B

Fxananassa gene IDF. vesca orthologA. thaliana orthologGene DescriptionRelation with defense/biological functionCC vs. CI
Fold ChangeFDR q-value
INVASION SENSING
M6F8gene29223AT1G57680G-Protein coupled receptor 1G-protein coupled receptor−1.993.95E-02
M20C3gene24345AT2G32240Leucine-rich repeat-containing proteinLRR protein, plant receptor−1.933.95E-02
SIGNAL TRANSDUCTION
M18F3gene21849AT5G43010Regulatory particle AAA-ATPase 4A/Proteasome complexRegulatory ATPase, Proteasome complex−2.023.95E-02
M29G3#gene25430AT2G22990Serine carboxypeptidasePeptidase, Glucosinolate and phenylpropanoid pathway−1.88 (−5.1 ± 2.2)3.95E-02
M5E3gene12921AT1G74960Beta-ketoacyl-ACP synthaseFatty acid biosynthesis−1.803.95E-02
M26F4gene09121AT5G67090Subtilisin-like serine endopeptidasePeptidase−1.783.95E-02
M22F5gene18417AT5G02310Protein ubiquitination component of the N-end ruleUbiquitin ligase, Proteasome complex−1.763.95E-02
NEW PROTEIN SYNTHESIS AND SECRETION
M10H10gene17514AT2G32700LEUNIG_homolog transcriptional correpresorTranscription represor−2.393.95E-02
M28F7gene25662AT5G0296040S Ribosomal protein S12/S23Protein synthesis−2.153.95E-02
M22E3gene12861AT5G53430Histone methyltransferaseIndirect transcription regulation−1.863.95E-02
M22E11gene15974.3utrAT1G15750TOPLESS transcriptional correpresorTranscription represor−1.853.95E-02
M22D5gene31183.3utrAT1G22910RNA-binding (RRM/RBD/RNP motifs) family proteinRNA metabolism−1.783.95E-02
M21G2gene29663AT1G29170SCAR family memberCitoesqueleton−1.753.95E-02
DIRECT DEFENSES
M29H6#gene32347AT4G22880Leucoanthocyanidin dioxygenase (LDOX)Secondary metabolism−1.91 (−4.3 ± 2.11)3.95E-02
M21F3gene11045AT1G36370Serine hydroxymetyltransferaseREDOX production−1.903.95E-02
M29C12gene21346AT5G05270Chalcone-flavanone isomeraseSecondary metabolism−1.893.95E-02
M19C6gene26641AT5G15870Glycosyl hydrolase family 81 proteinCell wall degradation, PR protein family−1.763.95E-02
HORMONE-DEPENDENT PATHWAYS
M18H1gene14092AT1G07590Tetratricopeptide repeat (TPR)-like superfamily proteinTetratricopeptide repeat−1.823.95E-02
M15G5gene02397AT4G03550Glucan / Callose synthaseNegative regulator SA dependent defenses−1.803.95E-02
NO OBVIOUSLY RELATED TO DEFENSE RESPONSE
M8D2gene14995AT5G17920Methionine synthase−2.203.95E-02
M9F8gene16275AT4G39970Haloacid dehalogenase-like hydrolase−2.023.95E-02
M7B2gene10408AT3G03890Flavin mononucleotide binding−1.943.95E-02
M14A10gene29476AT5G52820WD-40 repeat CUL4 RING ubiquitin ligase complex−1.943.95E-02
M5B7gene09169AT1G48380DNA binding protein ROOT HAIRLESS 1, component of the topoisomerase VI complex−1.923.95E-02
M18D12gene20804AT2G22530Alkaline-phosphatase-like family protein−1.833.95E-02
M18A11gene08921AT5G47470Nodulin transporter family protein−1.833.95E-02
M28A7gene15006AT5G10840Endomembrane protein 70 protein family−1.813.95E-02
M26H5gene18624AT1G01090Pyruvate dehydrogenase alpha−1.783.95E-02
M11B2gene07537AT3G13990Kinase-related protein−1.763.95E-02

Down-regulated genes by Colletotrichum acutatum in crown tissue of Fragaria x ananassa, cultivar Camarosa.

Genes were considered as differentially expressed if they fulfilled a FDR < 0.05 after a SAM test analysis and the fold-change was higher that 1.75-fold between the compared conditions. Fold change values represent the ratio of cv. Camarosa mock (CC) vs. infected (CI), transformed by: -1/fold-change for better understanding of values. (#) indicates genes further analyzed by real time RTqPCR to validate microarray result. Their relative expression value at 5dpi is shown as (media ± SD) in the “Fold Change” column. (

) marks no obvious detection of F. vesca ortolog gene due to putative fails by automated gene prediction (see also Table S3). 3utr, and 5utr, indicate F × ananassa sequences matching untranslated regions of the corresponding F. vesca gene (see also Table S3). Regulated genes are grouped in sections accordingly to their role in different steps of the defense response against C. acutatum (see Table S4 for associated references).

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 S3S5). 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

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

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

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

Discussion

Disease progress

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.

SA-signaling 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.

JA-signaling pathway

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 [73]. 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

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.

Statements

Author contributions

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.

Acknowledgments

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].

Conflict of interest

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.

Supplementary material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2016.01036

    Abbreviations

  • EPS

    exopolysaccharide

  • ER-resident

    endoplasmic reticulum-resident

  • 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

    PAMP-triggered immunity

  • qPCR

    quantitative PCR

  • R protein

    resistance protein

  • SA

    salicylic acid

  • SEA

    singular enrichment analysis

References

  • 1

    AcostaI. F.FarmerE. E. (2010). Jasmonates. Arabidopsis Book8:e0129. 10.1199/tab.0129

  • 2

    Al-ShahrourF.Díaz-UriarteR.DopazoJ. (2004). FatiGO: a web tool for finding significant associations of Gene Ontology terms with groups of genes. Bioinformatics20, 578580. 10.1093/bioinformatics/btg455

  • 3

    AltschulS. F.GishW.MillerW.MyersE. W.LipmanD. J. (1990). Basic local alignment search tool. J. Mol. Biol.215, 403410. 10.1016/S0022-2836(05)80360-2

  • 4

    Amil-RuizF.Blanco-PortalesR.Muñoz-BlancoJ.CaballeroJ. L. (2011). The strawberry plant defence mechanism: a molecular review. Plant Cell Physiol.52, 18731903. 10.1093/pcp/pcr136

  • 5

    Amil-RuizF.Encinas-VillarejoS.de los SantosB.Muñoz-MéridaA.MercadoJ. A.TrellesO.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 LitzR. E.FoltaK.TalonM.Pliego AlfaroF. (Lisbon: ISHS), 4750.

  • 6

    Amil-RuizF.Garrido-GalaJ.Blanco-PortalesR.FoltaK. M.Muñoz-BlancoJ.CaballeroJ. L. (2013). Identification and validation of reference genes for transcript normalization in strawberry (Fragaria × ananassa) defense responses. PLoS ONE8:e70603. 10.1371/journal.pone.0070603

  • 7

    ArroyoF. T.MorenoJ.García-HerdugoG.de los SantosB.BarrauC.PorrasM.et al. (2005). Ultrastructure of the early stages of Colletotrichum acutatum infection of strawberry tissues. Can. J. Bot.83, 491500. 10.1139/b05-022

  • 8

    AshburnerM.BallC. A.BlakeJ. A.BotsteinD.ButlerH.CherryJ. M.et al. (2000). Gene ontology: tool for the unification of biology. Nat. Genet.25, 2529. 10.1038/75556

  • 9

    AvivD. H.RustérucciC.IiiB. F. H.DietrichR. A.ParkerJ. E.DanglJ. L. (2002). Runaway cell death, but not basal disease resistance, in lsd1 is SA- and NIM1/NPR1-dependent. Plant J.29, 381391. 10.1046/j.0960-7412.2001.01225.x

  • 10

    BakS.BeissonF.BishopG.HambergerB.HöferR.PaquetteS.et al. (2011). Cytochromes P450. Arabidopsis Book9:e0144. 10.1199/tab.0144

  • 11

    BaniwalS. K.BhartiK.ChanK. Y.FauthM.GanguliA.KotakS.et al. (2004). Heat stress response in plants: a complex game with chaperones and more than twenty heat stress transcription factors. J. Biosci.29, 471487. 10.1007/BF02712120

  • 12

    BhattacharjeeS. (2012). The language of reactive oxygen species signaling in plants. J. Bot.2012, 22. 10.1155/2012/985298

  • 13

    BirkenbihlR. P.DiezelC.SomssichI. E. (2012). Arabidopsis WRKY33 Is a key transcriptional regulator of hormonal and metabolic responses toward botrytis cinerea infection. Plant Physiol.159, 266285. 10.1104/pp.111.192641

  • 14

    CaarlsL.PieterseC. M. J.Van WeesS. C. M. (2015). How salicylic acid takes transcriptional control over jasmonic acid signaling. Front. Plant Sci.6:170. 10.3389/fpls.2015.00170

  • 15

    Casado-DíazA.Encinas-VillarejoS.de los SantosB.SchiliròE.Yubero-SerranoE.-M.Amil-RuizF.et al. (2006). Analysis of strawberry genes differentially expressed in response to Colletotrichum infection. Physiol. Plant.128, 633650. 10.1111/j.1399-3054.2006.00798.x

  • 16

    CharngY. Y.LiuH. C.LiuN. Y.ChiW. T.WangC. N.ChangS. H.et al. (2007). A heat-inducible transcription factor, HsfA2, is required for extension of acquired thermotolerance in Arabidopsis. Plant Physiol.143, 251262. 10.1104/pp.106.091322

  • 17

    ChungS. H.RosaC.ScullyE. D.PeifferM.TookerJ. F.HooverK.et al. (2013). Herbivore exploits orally secreted bacteria to suppress plant defenses. Proc. Natl. Acad. Sci. U.S.A.110, 1572815733. 10.1073/pnas.1308867110

  • 18

    CurryK. J.AbrilM.AvantJ. B.SmithB. J. (2002). Strawberry anthracnose: histopathology of Colletotrichum acutatum and C. fragariae. Phytopathology92, 10551063. 10.1094/PHYTO.2002.92.10.1055

  • 19

    DebodeJ.Van HemelrijckW.BaeyenS.CreemersP.HeungensK.MaesM. (2009). Quantitative detection and monitoring of Colletotrichum acutatum in strawberry leaves using real-time PCR. Plant Pathol.58, 504514. 10.1111/j.1365-3059.2008.01987.x

  • 20

    DempseyD. M. A.KlessigD. F. (2012). SOS – too many signals for systemic acquired resistance?Trends Plant Sci.17, 538545. 10.1016/j.tplants.2012.05.011

  • 21

    DoddsP. N.RathjenJ. P. (2010). Plant immunity: towards an integrated view of plant–pathogen interactions. Nat. Rev. Genet.11, 539548. 10.1038/nrg2812

  • 22

    DoehlemannG.WahlR.HorstR. J.VollL. M.UsadelB.PoreeF.et al. (2008). Reprogramming a maize plant: transcriptional and metabolic changes induced by the fungal biotroph Ustilago maydis. Plant J.56, 181195. 10.1111/j.1365-313X.2008.03590.x

  • 23

    DongX. (2004). NPR1, all things considered. Curr. Opin. Plant Biol.7, 547552. 10.1016/j.pbi.2004.07.005

  • 24

    DurgbanshiA.ArbonaV.PozoO.MierschO.SanchoJ. V.Gomez-CadenasA. (2005). Simultaneous determination of multiple phytohormones in plant extracts by liquid chromatography−electrospray tandem mass spectrometry. J. Agric. Food Chem.53, 84378442. 10.1021/jf050884b

  • 25

    El OirdiM.El RahmanT. A.RiganoL.El HadramiA.RodriguezM. C.DaayfF.et al. (2011). Botrytis cinerea manipulates the antagonistic effects between immune pathways to promote disease development in tomato. Plant Cell Online23, 24052421. 10.1105/tpc.111.083394

  • 26

    Encinas-VillarejoS.MaldonadoA. M.Amil-RuizF.de los SantosB.RomeroF.Pliego-AlfaroF.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, 30433065. 10.1093/jxb/erp152

  • 27

    FreemanS.HorowitzS.SharonA. (2001). Pathogenic and nonpathogenic lifestyles in Colletotrichum acutatum from strawberry and other plants. Phytopathology91, 986992. 10.1094/PHYTO.2001.91.10.986

  • 28

    FreemanS.KatanT.ShabiE. (1998). Characterization of Colletotrichum species responsible for anthracnose diseases of various fruits. Plant Dis.82, 596605. 10.1094/PDIS.1998.82.6.596

  • 29

    GfellerA.DubugnonL.LiechtiR.FarmerE. E. (2010). Jasmonate biochemical pathway. Sci. Signal.3:cm3. 10.1126/scisignal.3109cm3

  • 30

    Grellet-BournonvilleC. F.Martinez-ZamoraM. G.CastagnaroA. P.Díaz-RicciJ. C. (2012). Temporal accumulation of salicylic acid activates the defense response against Colletotrichum in strawberry. Plant Physiol. Biochem.54, 1016. 10.1016/j.plaphy.2012.01.019

  • 31

    GuidarelliM.CarboneF.MourguesF.PerrottaG.RosatiC.BertoliniP.et al. (2011). Colletotrichum acutatum interactions with unripe and ripe strawberry fruits and differential responses at histological and transcriptional levels. Plant Pathol.60, 685697. 10.1111/j.1365-3059.2010.02423.x

  • 32

    HeidrichK.WirthmuellerL.TassetC.PouzetC.DeslandesL.ParkerJ. E. (2011). Arabidopsis EDS1 connects pathogen effector recognition to cell compartment–specific immune responses. Science334, 14011404. 10.1126/science.1211641

  • 33

    HorowitzS.FreemanS.SharonA. (2002). Use of green fluorescent protein-transgenic strains to study pathogenic and nonpathogenic lifestyles in Colletotrichum acutatum. Phytopathology92, 743749. 10.1094/PHYTO.2002.92.7.743

  • 34

    IkedaM.MitsudaN.Ohme-TakagiM. (2011). Arabidopsis HsfB1 and HsfB2b act as repressors of the expression of heat-inducible hsfs but positively regulate the acquired thermotolerance. Plant Physiol.157, 12431254. 10.1104/pp.111.179036

  • 35

    IkedaM.Ohme-TakagiM. (2009). A novel group of transcriptional repressors in Arabidopsis. Plant Cell Physiol.50, 970975. 10.1093/pcp/pcp048

  • 36

    KhanA. A.ShihD. S. (2004). Molecular cloning, characterization, and expression analysis of two class II chitinase genes from the strawberry plant. Plant Sci.166, 753762. 10.1016/j.plantsci.2003.11.015

  • 37

    KrinkeO.RuellandE.ValentováO.VergnolleC.RenouJ.-P.TaconnatL.et al. (2007). Phosphatidylinositol 4-Kinase activation is an early response to salicylic acid in Arabidopsis suspension cells. Plant Physiol.144, 13471359. 10.1104/pp.107.100842

  • 38

    KubigsteltigI.LaudertD.WeilerE. W. (1999). Structure and regulation of the Arabidopsis thaliana allene oxide synthase gene. Planta208, 463471. 10.1007/s004250050583

  • 39

    LeandroL. F. S.GleasonM. L.NutterF. W.WeguloS. N.DixonP. M. (2001). Germination and sporulation of Colletotrichum acutatum on symptomless strawberry leaves. Phytopathology91, 659664. 10.1094/PHYTO.2001.91.7.659

  • 40

    Leon-ReyesA.Van der DoesD.De LangeE.DelkerC.WasternackC.Van WeesS.et al. (2010). Salicylate-mediated suppression of jasmonate-responsive gene expression in Arabidopsis is targeted downstream of the jasmonate biosynthesis pathway. Planta232, 14231432. 10.1007/s00425-010-1265-z

  • 41

    LiJ.BraderG.KariolaT.Tapio PalvaE. (2006). WRKY70 modulates the selection of signaling pathways in plant defense. Plant J.46, 477491. 10.1111/j.1365-313X.2006.02712.x

  • 42

    LiJ.BraderG.PalvaE. T. (2004). The WRKY70 transcription factor: a node of convergence for Jasmonate-mediated and Salicylate-mediated signals in plant defense. Plant Cell16, 319331. 10.1105/tpc.016980

  • 43

    LiuP.-P.von DahlC. C.ParkS.-W.KlessigD. 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, 17621768. 10.1104/pp.110.171694

  • 44

    LodhaT. D.BasakJ. (2012). Plant–pathogen interactions: what microarray tells about it?Mol. Biotechnol.50, 8797. 10.1007/s12033-011-9418-2

  • 45

    López-RáezJ. A.VerhageA.FernándezI.GarcíaJ. M.Azcón-AguilarC.FlorsV.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, 25892601. 10.1093/jxb/erq089

  • 46

    MaasJ. (1998). Compendium of Strawberry Diseases. St Paul, MN: American Phytopathological Society.

  • 47

    MakowskiR. M. D.MortensenK. (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, 15451552. 10.1017/S0953756298006960

  • 48

    MaleckK.LevineA.EulgemT.MorganA.SchmidJ.LawtonK. A.et al. (2000). The transcriptome of Arabidopsis thaliana during systemic acquired resistance. Nat. Genet.26, 403410. 10.1038/82521

  • 49

    MarcelS.SawersR.OakeleyE.AnglikerH.PaszkowskiU. (2010). Tissue-adapted invasion strategies of the rice blast fungus Magnaporthe oryzae. Plant Cell22, 31773187. 10.1105/tpc.110.078048

  • 50

    NdamukongI.AbdallatA. A.ThurowC.FodeB.ZanderM.WeigelR.et al. (2007). SA-inducible Arabidopsis glutaredoxin interacts with TGA factors and suppresses JA-responsive PDF1.2 transcription. Plant J.50, 128139. 10.1111/j.1365-313X.2007.03039.x

  • 51

    Pajerowska-MukhtarK. M.WangW.TadaY.OkaN.TuckerC. L.FonsecaJ. P.et al. (2012). The HSF-like transcription factor TBF1 is a major molecular switch for plant growth-to-defense transition. Curr. Biol.22, 103112. 10.1016/j.cub.2011.12.015

  • 52

    Peña-CortésH.BarriosP.DortaF.PolancoV.SánchezC.SánchezE.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, 246260. 10.1007/s00344-004-0035-1

  • 53

    PernasM.RyanE.DolanL. (2010). SCHIZORIZA controls tissue system complexity in plants. Curr. Biol.20, 818823. 10.1016/j.cub.2010.02.062

  • 54

    PieterseC. M. J.Leon-ReyesA.Van der EntS.Van WeesS. C. M. (2009). Networking by small-molecule hormones in plant immunity. Nat. Chem. Biol.5, 308316. 10.1038/nchembio.164

  • 55

    RahmanT. A. E.OirdiM. E.Gonzalez-LamotheR.BouarabK. (2012). Necrotrophic pathogens use the salicylic acid signaling pathway to promote disease development in tomato. Mol. Plant Microbe Interact.25, 15841593. 10.1094/MPMI-07-12-0187-R

  • 56

    RenC.-M.ZhuQ.GaoB.-D.KeS.-Y.YuW.-C.XieD.-X.et al. (2008). Transcription factor WRKY70 displays important but no indispensable roles in jasmonate and salicylic acid signaling. J. Integr. Plant Biol.50, 630637. 10.1111/j.1744-7909.2008.00653.x

  • 57

    RietzS.StammA.MalonekS.WagnerS.BeckerD.Medina-EscobarN.et al. (2011). Different roles of Enhanced Disease Susceptibility1 (EDS1) bound to and dissociated from Phytoalexin Deficient4 (PAD4) in Arabidopsis immunity. New Phytol.191, 107119. 10.1111/j.1469-8137.2011.03675.x

  • 58

    Robert-SeilaniantzA.GrantM.JonesJ. D. G. (2011). Hormone crosstalk in plant disease and defense: more than just JASMONATE-SALICYLATE antagonism. Annu. Rev. Phytopathol.49, 317343. 10.1146/annurev-phyto-073009-114447

  • 59

    RodriguezM. C. S.PetersenM.MundyJ. (2010). Mitogen-activated protein kinase signaling in plants. Annu. Rev. Plant Biol.61, 621649. 10.1146/annurev-arplant-042809-112252

  • 60

    RouhierN.CouturierJ.JacquotJ.-P. (2006). Genome-wide analysis of plant glutaredoxin systems. J. Exp. Bot.57, 16851696. 10.1093/jxb/erl001

  • 61

    RueppA.ZollnerA.MaierD.AlbermannK.HaniJ.MokrejsM.et al. (2004). The FunCat, a functional annotation scheme for systematic classification of proteins from whole genomes. Nucleic Acids Res.32, 55395545. 10.1093/nar/gkh894

  • 62

    RustérucciC.AvivD. H.HoltB. F.DanglJ. L.ParkerJ. 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 Online13, 22112224. 10.1105/tpc.13.10.2211

  • 63

    SarowarS.ZhaoY.Soria-GuerraR. E.AliS.ZhengD.WangD.et al. (2011). Expression profiles of differentially regulated genes during the early stages of apple flower infection with Erwinia amylovora. J. Exp. Bot.62, 48514861. 10.1093/jxb/err147

  • 64

    SasakiY.AsamizuE.ShibataD.NakamuraY.KanekoT.AwaiK.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, 153161. 10.1093/dnares/8.4.153

  • 65

    SchenkP. M.KazanK.MannersJ. M.AndersonJ. P.SimpsonR. S.WilsonI. W.et al. (2003). Systemic gene expression in Arabidopsis during an Incompatible Interaction with Alternaria brassicicola. Plant Physiol.132, 9991010. 10.1104/pp.103.021683

  • 66

    SimpsonD. W. (1991). Resistance to Botrytis cinerea in pistillate genotypes of the cultivated strawberry Fragaria ananassa. J. Hortic. Sci. Biotechnol.66, 719724

  • 67

    ShulaevV.SargentD. J.CrowhurstR. N.MocklerT. C.FolkertsO.DelcherA. L.et al. (2011). The genome of woodland strawberry (Fragaria vesca). Nat. Genet.43, 109116. 10.1038/ng.740

  • 68

    SongW. C.FunkC. D.BrashA. 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, 85198523. 10.1073/pnas.90.18.8519

  • 69

    SpoelS. H.DongX. (2012). How do plants achieve immunity? Defence without specialized immune cells. Nat. Rev. Immunol.12, 89100. 10.1038/nri3141

  • 70

    SpoelS. H.JohnsonJ. S.DongX. (2007). Regulation of tradeoffs between plant defenses against pathogens with different lifestyles. Proc. Natl. Acad. Sci. U.S.A.104, 1884218847. 10.1073/pnas.0708139104

  • 71

    StaswickP. E.TiryakiI. (2004). The oxylipin signal jasmonic acid is activated by an enzyme that conjugates it to isoleucine in Arabidopsis. Plant Cell Online16, 21172127. 10.1105/tpc.104.023549

  • 72

    ten HoveC. A.WillemsenV.de VriesW. J.van DijkenA.ScheresB.HeidstraR. (2010). SCHIZORIZA encodes a nuclear factor regulating asymmetry of stem cell divisions in the Arabidopsis root. Curr. Biol.20, 452457. 10.1016/j.cub.2010.01.018

  • 73

    TohidfarM.KhosraviS. (2015). Transgenic crops with an improved resistance to biotic stresses. A review. Biotechnol. Agron. Soc. Environ.19, 6270.

  • 74

    TurnerJ. G.EllisC.DevotoA. (2002). The jasmonate signal pathway. Plant Cell Online14(suppl. 1), S153S164. 10.1105/tpc.000679

  • 75

    TusherV. G.TibshiraniR.ChuG. (2001). Significance analysis of microarrays applied to the ionizing radiation response. Proc. Natl. Acad. Sci. U.S.A.98, 51165121. 10.1073/pnas.091062498

  • 76

    UknesS.Mauch-ManiB.MoyerM.PotterS.WilliamsS.DincherS.et al. (1992). Acquired resistance in Arabidopsis. Plant Cell Online4, 645656. 10.1105/tpc.4.6.645

  • 77

    VargasW. A.MartínJ. M. S.RechG. E.RiveraL. P.BenitoE. P.Díaz-MínguezJ. M.et al. (2012). Plant defense mechanisms are activated during biotrophic and necrotrophic development of colletotricum graminicola in maize. Plant Physiol.158, 13421358. 10.1104/pp.111.190397

  • 78

    VenugopalS. C.JeongR.-D.MandalM. K.ZhuS.Chandra-ShekaraA. C.XiaY.et al. (2009). Enhanced disease susceptibility 1 and salicylic acid act redundantly to regulate resistance gene-mediated signaling. PLoS Genet.5:e1000545. 10.1371/journal.pgen.1000545

  • 79

    VlotA. C.LiuP.-P.CameronR. K.ParkS.-W.YangY.KumarD.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, 445456. 10.1111/j.1365-313X.2008.03618.x

  • 80

    WangD.AmornsiripanitchN.DongX. (2006). A genomic approach to identify regulatory nodes in the transcriptional network of systemic acquired resistance in plants. PLoS Pathog.2:e123. 10.1371/journal.ppat.0020123

  • 81

    WangD.WeaverN. D.KesarwaniM.DongX. (2005). Induction of protein secretory pathway is required for systemic acquired resistance. Science308, 10361040. 10.1126/science.1108791

  • 82

    WangG.-F.SeaboltS.HamdounS.NgG.ParkJ.LuH. (2011). Multiple Roles of WIN3 in regulating disease resistance, cell death, and flowering time in Arabidopsis. Plant Physiol.156, 15081519. 10.1104/pp.111.176776

  • 83

    WiermerM.FeysB. J.ParkerJ. E. (2005). Plant immunity: the EDS1 regulatory node. Curr. Opin. Plant Biol.8, 383389. 10.1016/j.pbi.2005.05.010

  • 84

    WindramO.MadhouP.McHattieS.HillC.HickmanR.CookeE.et al. (2012). Arabidopsis defense against botrytis cinerea: chronology and regulation deciphered by high-resolution temporal transcriptomic analysis. Plant Cell Online24, 35303557. 10.1105/tpc.112.102046

  • 85

    YangY. H.DudoitS.LuuP.LinD. M.PengV.NgaiJ.et al. (2002). Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res.30:e15. 10.1093/nar/30.4.e15

  • 86

    YangY.XuR.MaC. J.VlotA. C.KlessigD. F.PicherskyE. (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, 10341045. 10.1104/pp.108.118224

  • 87

    ZamoraM. G. M.BournonvilleC. G.CastagnaroA. P.RicciJ. 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, 412420. 10.1071/FP11275

  • 88

    ZanderM.ChenS.ImkampeJ.ThurowC.GatzC. (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, 831840. 10.1093/mp/ssr113

  • 89

    ZiemannM.BhaveM.ZachgoS. (2009). Origin and Diversification of Land Plant CC-Type Glutaredoxins. Genome Biol. Evol.1, 265277. 10.1093/gbe/evp025

Summary

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

Volume

7 - 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

Updates

Copyright

*Correspondence: José L. Caballero

This article was submitted to Plant Biotic Interactions, a section of the journal Frontiers in Plant Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics