Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Disclosing proteins in the leaves of cork oak plants associated with the immune response to Phytophthora cinnamomi inoculation in the roots: A long-term proteomics approach

  • Ana Cristina Coelho ,

    Roles Conceptualization, Data curation, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Visualization, Writing – original draft

    acoelho@ualg.pt

    Affiliations Center for Electronic, Optoelectronic and Telecommunications (CEOT), University of Algarve, Faro, Portugal, Escola Superior de Educação e Comunicação (ESEC), University of Algarve, Faro, Portugal

    ⨯
  • Rosa Pires,

    Roles Investigation, Writing – review & editing

    Affiliation Center for Electronic, Optoelectronic and Telecommunications (CEOT), University of Algarve, Faro, Portugal

    ⨯
  • Gabriela Schütz,

    Roles Data curation, Formal analysis, Writing – review & editing

    Affiliations Center for Electronic, Optoelectronic and Telecommunications (CEOT), University of Algarve, Faro, Portugal, Instituto Superior de Engenharia, University of Algarve, Faro, Portugal

    ⨯
  • Cátia Santa,

    Roles Data curation, Formal analysis, Funding acquisition, Writing – review & editing

    Affiliations CNC—Center for Neuroscience and Cell Biology, University of Coimbra, Coimbra, Portugal, Institute for Interdisciplinary Research (IIIUC), University of Coimbra, Coimbra, Portugal

    ⨯
  • Bruno Manadas,

    Roles Data curation, Formal analysis, Funding acquisition, Writing – review & editing

    Affiliation CNC—Center for Neuroscience and Cell Biology, University of Coimbra, Coimbra, Portugal

    ⨯
  • Patrícia Pinto

    Roles Data curation, Formal analysis, Writing – review & editing

    Affiliation Center for Marine Sciences (CCMAR), University of Algarve, Faro, Portugal

    ⨯

Abstract

The pathological interaction between oak trees and Phytophthora cinnamomi has implications in the cork oak decline observed over the last decades in the Iberian Peninsula. During host colonization, the phytopathogen secretes effector molecules like elicitins to increase disease effectiveness. The objective of this study was to unravel the proteome changes associated with the cork oak immune response triggered by P. cinnamomi inoculation in a long-term assay, through SWATH-MS quantitative proteomics performed in the oak leaves. Using the Arabidopis proteome database as a reference, 424 proteins were confidently quantified in cork oak leaves, of which 80 proteins showed a p-value below 0.05 or a fold-change greater than 2 or less than 0.5 in their levels between inoculated and control samples being considered as altered. The inoculation of cork oak roots with P. cinnamomi increased the levels of proteins associated with protein-DNA complex assembly, lipid oxidation, response to endoplasmic reticulum stress, and pyridine-containing compound metabolic process in the leaves. In opposition, several proteins associated with cellular metabolic compound salvage and monosaccharide catabolic process had significantly decreased abundances. The most significant abundance variations were observed for the Ribulose 1,5-Bisphosphate Carboxylase small subunit (RBCS1A), Heat Shock protein 90–1 (Hsp90-1), Lipoxygenase 2 (LOX2) and Histone superfamily protein H3.3 (A8MRLO/At4G40030) revealing a pertinent role for these proteins in the host-pathogen interaction mechanism. This work represents the first SWATH-MS analysis performed in cork oak plants inoculated with P. cinnamomi and highlights host proteins that have a relevant action in the homeostatic states that emerge from the interaction between the oomycete and the host in the long term and in a distal organ.

Citation: Coelho AC, Pires R, Schütz G, Santa C, Manadas B, Pinto P (2021) Disclosing proteins in the leaves of cork oak plants associated with the immune response to Phytophthora cinnamomi inoculation in the roots: A long-term proteomics approach. PLoS ONE 16(1): e0245148. https://doi.org/10.1371/journal.pone.0245148

Editor: Sara Amancio, Universidade de Lisboa Instituto Superior de Agronomia, PORTUGAL

Received: September 18, 2020; Accepted: December 22, 2020; Published: January 22, 2021

Copyright: © 2021 Coelho et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD021455.

Funding: This work was financially supported by FCT, integrated in projects UID/Multi/00631/2013, UID/Multi/00631/2019 and UIDB/00631/2020 CEOT BASE to CEOT and ACC, GS and RP; UIDB/04326/2020 to CCMAR; fellowship SFRH/BPD/84033/2012 and researcher contract with the University of Algarve under Norma Transitória-DL57/2016/CP1361/CT0015 to PP; contract NIBAP (ALG-01-0247-FEDER-037303) to RP; projects POCI-01-0145-FEDER-007440 (Ref. UIDB/04539/2020), POCI-01-0145-FEDER-016428 (Ref. SAICTPAC/0010/2015), POCI-01-0145-FEDER-029311 (Ref. PTDC/BTM-TEC/29311/2017), POCI-01-0145-FEDER-30943 (Ref. PTDC/MECPSQ/30943/2017) and PTDC/MED-NEU/27946/2017 to CNC, BM and CS. The work at CNC was also funded by the National Mass Spectrometry Network (RNEM) under contract POCI-01-0145-FEDER-402-022125 (Ref. ROTEIRO/0028/2013). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

The soil-borne oomycete Phytophthora cinnamomi infects the roots of cork oak (Quercus suber) plants, induces necrotic lesions, and the loss of fine roots [1,2]. This evidence, combined with other factors, are the hallmark for the decline of the cork oak savanna-like ecosystem in Portugal (cork oak montado) and Spain (cork oak dehesa). Climate changes is reducing water availability (drought) [3], and the effectiveness of roots in absorbing water is affected by the health status of the plant [4,5], which can become less effective in accessing groundwater during drought [6]. Insect colonization [7] and fungal infections [8,9] can weaken the tree's defence system and thus contribute to the decline. To help maintain the sustainability of the cork oak agro-forests, the recommended focus is to adopt good management practices [10].

During inter and intracellular cork oak colonization by P. cinnamomi, small 10 kDa proteins (elicitins) are secreted by the oomycete and increases disease effectiveness. This has been demonstrated by studying a β-cinnamomin silenced P. cinnamomi strain, which acted as a weaker pathogen against cork oak when compared to the virulence revealed by the wild type [11,12]. In the roots of the narrow-leafed lupin (Lupinus angustifolius) infected with P. cinnamomi, the expression of β-cinnamomin starts to be detected as early as 24 h post-inoculation and follows the development of the mycelium into the host, anchored to a mycelial cell wall protein, emphasizing the recognition of these proteins as virulence factors [13]. However, effector molecules from the RxLR, CRN (for Crinkling and Necrosis) and Nep1-like (NLPs) protein families are also potentially secreted, encoded by the 171 RxLR, 72 NLPs and 29 CRN putative genes present in the genomes (78 Mb) of three P. cinnamomi isolates, being able to suppress or bypass the plant basic defence responses [14]. The molecular mechanisms by which the effector molecules act are largely unknown, although the entry of some effector proteins into the plant host cells is known to follow a mechanism of endocytosis after binding to receptor molecules of phosphatidylinositol-3-phosphate (PI-3-P) mediated by the effector RxLR domain [15,16]. In the nucleus, the effectors control reactions that trigger host cell death or hypersensitive responses (HR) [17,18], and in the nucleolus, they can act as modulators of histone acetyltransferases (HAT) to reprogram the plant defence gene expression and promote infection [19].

Following compatible or incompatible reactions with plants, oomycete compounds like lipids or carbohydrates referred to as Pathogen-Associated Molecular Patterns (PAMPs) and effector biomolecules elicit local resistance responses or PAMPs/effector-triggered immunity (PTI/ETI) in their hosts [20]. In Q. suber root cells, during the first 24 h of interaction with P. cinnamomi, metabolic patterns undergo a non-linear variation for compounds with carbohydrate, glycoconjugate and lipid groups [21]. At the transcriptomic level, the differential expression of genes encoding pathogenesis-related proteins was observed in avocado roots challenged with P. cinnamomi [22] and in stem tissues of Eucalyptus nitens infected with P. cinnamomi [23]. In a more detailed analysis of the transcriptome of chestnut roots inoculated with P. cinnamomi, the multiplicity of the defence responses becomes evident with the identification of genes related to the HR (hypersensitive response), cell wall strengthening, synthesis of flavonoids and systemic acquired resistance [24]. Further, resistance (R) genes coding to transmembrane proteins such as LRR receptor-like serine/threonine-protein kinase in two Castanea species [24] and CC-NB-LRR (coiled coil-nucleotide binding-leucine rich repeat) in cork oak [25] are also potentially associated to the recognition of effector molecules, eventually interacting, according to the gene-for-gene model [26]. Activation of these resistance proteins can result in the activation of mitogen-activated protein kinase (MAPK) signal transduction cascades, leading to transcription factor activation and transcription of responsive genes, and these cascades can also be activated by proteins sensitive to the production of reactive oxygen species (ROS, O2-, H2O2) [reviewed by 20,27–29].

Salicylic acid (SA)/salicylate is also a signaling molecule that plays a central role in PAMPs/effector-triggered immunity (PTI/ETI) and in the systemic acquired resistance (SAR). SAR is a type of immunity that extends to the entire plant beyond the site of infection, protecting the plant against a broad spectrum of pathogens [30,31]. The expression of a large number of pathogenesis-related genes is activated by nuclear transcription factors interacting with NPR1 monomers (nonexpressor of pathogenesis related 1), known as the main regulatory molecule of the SA-signaling pathway [32,33].

To overcome the harmful implications of P. cinnamomi on susceptible species of thousands of plants worldwide, one of the current challenges is the identification of molecular markers or physiological processes suitable for recognition of resistant or susceptible host plant species or varieties. Information about the constitutive expression level of pathogenesis-related genes in non-infected hosts and the reaction time mediating the recognition of the invader and the activation of local and systemic defence systems can contribute to this global goal, and was critical for the recognition of Castanea crenata as a less susceptible species than C. sativa [28]. In less susceptible avocado rootstocks, the physical and chemical composition of the host's tissues at the site of infection was critical to the effectiveness of P. cinnamomi zoospore germination and penetration, as the early deposition of callose instead of lignin near the site of hyphae penetration along the cell wall hindered the development of the oomycete's hyphae [34].

The hypothesis of the present study is that after inoculation of plant roots with a pathogen, an immune response is initiated that will lead to a new homeostatic state, with protein changes that can be detectable in the long-term, distally from the infection site. The aim was to identify and quantify proteins in the leaves of cork oak plants inoculated with P. cinnamomi in the roots and compare them to those in the leaves of non-inoculated plants, at 248 days post-inoculation, using SWATH-MS proteomics [35]. SWATH-MS (Sequential Window Acquisition of all Theoretical Mass Spectra) is a quantitative, label-free and unbiased proteomics method that is able to acquire information about virtually every ion (in this case peptides), introduced into the mass spectrometer [36]. SWATH is a promising strategy for the quantitative screening of a large number of proteins that has previously been applied in the field of plant biology [37–39] and recognized as a valuable tool for the comprehensive study of proteins in plants [40,41].

The leaves are a distal organ that can be sampled in a minimally invasive way in adult trees, so they can also be a potential organ for practical monitoring of infection or resistance. Four hundred and twenty-four proteins were identified in the cork oak leaves, and a subset of 80 proteins showed differential levels between inoculated and control plants, being considered responsive to P. cinnamomi. These included 18 proteins associated with several gene ontologies (GO) biological processes, and their potential role in the cork oak immune response is discussed. The GO cellular component “stromules” was also significantly enriched among the differential proteins, indicating that communication between cellular organelles may be important in the cork oak immune response to P. cinnamomi.

Materials and methods

The design of the project included several experimental procedures operated at different time points. In the first phase, the biological material was prepared consisting of twelve cork oak seedlings, germinated from seeds, with half of these plants being inoculated with P. cinnamomi. The following phases started 248 days after inoculation and included the harvesting of the leaves from each plant for protein extraction and subsequent SWATH-MS proteomics. The experimental assay ended with the bioinformatic annotation and quantification of proteins present in the extracts of each plant.

Biological material

Cork oak plants used in this experimental project were germinated from acorns taken from six cork oak trees located in Cachopo, Algarve, Portugal (S1 Fig). Parental cork oak trees referenced as S1.1, S2.1, S4.1, S5.1, S7.1, and S8.1 showed signs of decline at distinct stages of progression, based on visual observation of the canopy defoliation level typical of P. cinnamomi infection. The study included two experimental conditions with six biological replicates: 6 cork oak plants inoculated with the PA45 P. cinnamomi isolate and 6 non-inoculated plants. Seeds from six parental cork oak trees were germinated and were distributed between the control and inoculated groups so that each inoculated plant had a paired control from the same progenitor. S1 Table provides the cork oak references used in the study. PA 45 was isolated from the rhizosphere of cork oak trees that showed symptoms of decline in the Algarve region and its high virulence on cork oak seedlings was extensively studied [11,12,22]. To reconfirm the identity of the isolate as P. cinnamomi, DNA was extracted from PA 45 isolate and was used in PCR reactions with primers (95.422/96.007) designed for a colorimetric molecular assay [42] targeting the elicitin genes (GenBank accession number AJ000071).

For the preparation of control and inoculated plants, twelve 77-day-hold cork oak plants were removed from the germination alveoli, freeing most of the organic substrate that accompanied the roots, and were laid down on trays whose surface was protected with moist absorbent paper. Then, a 2 cm2 agar plug of P. cinnamomi mycelium isolate PA45, grown in clarified V8 (Campbell Soup) semi-solid agar, in the dark at 25°C [11] for 9 days, was placed mycelial surface down on the tap root of 6 cork oak plants—inoculated plants. The roots of the control plants were not exposed to non-colonised semi-solid agar plugs to prevent the growth of microorganisms present in unsterilized roots on the nutritious support (agar surface), whose interaction with plant tissues could elicit defence reactions not present in the natural plants. This situation is prevented in the inoculated samples due to the large amount of P. cinnamomi hyphae present on the surface of the agar plugs avoiding bacteria and other microorganisms from having acess to the nutritious support.

The roots of the inoculated and non-inoculated plants were covered with aluminum foil and kept on the moistened trays at 25°C for 48 hours.

Forty-eight hours after P. cinnamomi inoculation the plugs were removed and all plants were potted into a misture of planting soil (PFLANZ-ERDE) and sand (proportion 2/3 for 1/3) in free-drining plastic containers (Top Ø 16 cm; Base Ø 13 cm; H 33 cm), transferred outside and watered regularly to container capacity. S1 Fig outlines the procedure and timing of the experiment. After 248 days, the leaves of cork oak plants, inoculated and non-inoculated, were collected and immediately frozen in liquid nitrogen and stored at -80°C until further use for protein extraction.

SWATH-MS proteomics

Total protein extraction.

The optimized extraction of proteins from cork oak leaves included eight steps. 1) Leaf tissue (200 mg) was ground in a mortar and pestle in the presence of liquid nitrogen to obtain a fine powder. 2) Buffer 1 (1.25 mL/100 mg sample; DTT–Dithiothreitol 60 mM; 10% TCA-Trichloroacetic acid solubilized in acetone) was added to the mortar and samples were macerated in the presence of the buffer with the pestle. 3) The heterogeneous solution was transferred to a 2 mL microcentrifuge tube and incubated for 1 hour at -80°C. 4) The samples were centrifuged at 15,000× g for 15 min at 4°C and the supernatant was discarded. 5) The pellets were dissolved in 1 mL of Buffer 2 (2.5 ml/100 mg sample; DTT 60 mM solubilized in acetone), and the sample volume was divided into two microcentrifuge tubes followed by the addition of 750 μL of Buffer 2 to each tube. 6) These solutions were incubated for 1 hour at -80°C. Procedures 5 and 6 were repeated until the solution was clear green. 7) The samples were centrifuged at 15,000×g for 15 min at 4°C and the supernatant was discarded. 8) Pellets were dried and resuspended in 250 μL of SDS-PAGE buffer [TRIS Glycine buffer solution (25 mM TRIS; 192 mM Glycine; Sigma-Aldrich); 2% SDS-Sodium dodecyl sulfate] followed by incubation at 95°C for 5 min and centrifugation at 20,000x g for 15 min at 4°C. All reagents used were molecular biology grade.

Protein concentration in the samples was estimated using the 2D-Quant kit (GE Healthcare Life Sciences) with serum albumin as standard [43].

SWATH-MS strategy.

For the proteomic screening, the short GeLC-SWATH-MS strategy was used according to [44] with minor modifications. Briefly, 50 μg of each sample and a pooled sample per group (pool of the protein extracts for the six control or six inoculated samples) were subjected to in-gel digestion after a partial SDS-PAGE run. Then, LC-MS information was acquired in two different acquisition modes: information-dependent acquisition (IDA) of the pooled samples, and SWATH-MS (Sequential Windowed data-independent Acquisition of the Total High-resolution Mass Spectra) of each sample. Protein identification and library construction was performed using ProteinPilot™ (v5.0.1, Sciex), and compared with the Arabidopsis thaliana reference proteome (retrieved from https://www.uniprot.org/ in April 2018). In addition, protein identification was tested against the predicted proteins deduced from the recently published draft genome sequence of cork oak [45], available at http://corkoakdb.org/downloads (fileGCF_002906115.1_CorkOak1.0_protein.faa, accessed in November 2020). The relative quantification was performed using the SWATH™ processing plug-in for PeakView™ (v2.2, Sciex). For each experimental group, the average protein levels, standard deviation and percentage coefficient of variation (% CV) were calculated based on the quantification levels obtained for each individual with six biological replicates per group. The fold change (FC) between inoculated and control plants was calculated by dividing the respective median protein levels for all quantified proteins. Statistical comparisons between protein levels were carried out using the software SPSS v23 (IBM) and the non-parametric Mann Whitney U-test (MW). Proteins were considered as differentially modified when FC was greater than 2 or less than 0.5 or MW p-value was below 0.05.

The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE [42] partner repository with the dataset identifier PXD021455. A detailed description of LC-MS materials and methods are provided as supporting information (S1 File).

Enrichment analyses and hierarchical clustering.

Gene ontology (GO) and pathway (KEGG and Reactome) enrichment analyses were carried out as in [46,47], using Cytoscape v3.5.1 and ClueGO plug-in v2.5.2 [48,49], comparing the list of 80 differentially modified proteins against the Arabidopsis thaliana [organism 3702] set of GO biological process and cellular component databases from November 2017. Enrichment analyses were repeated using the databases updated in 2020 and the same general enriched terms were found (data not shown). The following settings were used for the ClueGO enrichment analysis (right-side): GO levels 3 to 8, Benjamini-Hochberg false discovery rate (FDR) correction with a cut-off at FDR<0.05 and minimum of three genes/4% for terms to be considered significant. The initial group size was set as 1, group merging at 50%, and Kappa-statistics score at 0.4.

Enrichment scores of the functionally related network groups were calculated as -Log2 [group FDR]. The leading terms of each enriched group were those with the lowest term FDR (highest enrichment score), which was used to name the respective group.

The hierarchical clustering of the 80 differentially modified proteins was analysed with Cluster 3.0 at http://bonsai.hgc.jp/~mdehoon/software/cluster/software.htm#ctv using normalized protein levels, applying the uncentered correlation and complete linkage options.

Results and discussion

Observation of the plants

Cork oak plants were in contact with the P. cinnamomi mycelium at the beginning of the experiment for 48 hours, with no (re)inoculation over the next eight months until the end of the experimental assay. During the first 24 hours of inoculation with P. cinnamomi strain PA45, the aerial apex of the inoculated plants wilted and after 48 hours, all the inoculated roots appeared necrotic at the inoculation site (S2 Fig). At this time there was no observable changes in the control plants. Seven months after inoculation, 1 month before the end of the experiment, it was not possible to distinguish control plants from plants inoculated with P. cinnamomi by visual observation of the aerial part (S2 Fig). The vegetative development of the plants looked similar in both experimental conditions, inoculated and non-inoculated. Although no foliar symptoms of P. cinnamomi infection were observed, the infection is expected to have spread beyond the inoculation site through zoospores released from sporangia who migrated into the irrigation water or through root to root contact.

The virulence of the PA45 strain had been previously tested in cork oak roots, inoculated under the same conditions as in the present study for 3 days [11]. Histological studies performed on colonized root tissue demonstrated the ability of the oomycete to invade the epidermis, cortical parenchyma and vascular cylinder both inter- and intra-cellularly, and to destroy host cells [11].

In nature, at infested sites, cork oak trees may succumb (sudden death) after the summer, without showing obvious previous symptoms of decline, or they can remain for years with symptoms of defoliation that slowly worsen over time (slow decline). P. cinnamomi has been isolated from roots of declined symptomatic trees and from infested soils throughout Portugal, and it is important to recognize that oomycete infection can be a determining factor for cork oak decline. However, the recovery of P. cinnamomi from declining trees does not provide information about the plants' responsiveness or vitality over time.

At the end of this experiment, the plants inoculated with P. cinnamomi in the form of a single event were visually asymptomatic for leaf fall or yellowing, plant height or number of leaves. Nevertheless, the molecular interaction between P. cinnamomi and the hosts may have occurred differently in each of the six plants due to the high molecular diversity characteristic of Q. suber species [50]. One of the pertinent questions is how to detect that a plant is or has been invaded/infected when it has no symptoms, avoiding the (re)isolation of the pathogen and the use of invasive methods. As the degree of tree defoliation is a symptom of decline and leaf harvesting is a method minimally invasive to adult trees, the search for molecular markers in the leaves of plants challenged with P. cinnamomi can be a valuable option. The defence responses induced in the host by P. cinnamomi in the long term and distant from the inoculation site establish a homeostatic state adapted to living with the invader. This new homeostatic state stands out when comparing (below) the type and amount of proteins present in the leaves of inoculated and non-inoculated plants.

Leaf cork oak proteome changes in response to P. cinnamomi inoculation

SWATH-MS Proteomics analysis performed on leaf samples was used to characterize the proteome of cork oak plants inoculated with P. cinnamomi, and compared with the proteome of non-inoculated plants. With this technique, 12 individual protein profiles were obtained, and protein abundances were quantified in each of the leaf extracts. Thus, the protein profiles obtained for the six biological replicates in the two experimental conditions (control and inoculated) reflect the genetic variability of the Q. suber species, assuming the average of the results a value closer to reality.

Four proteome databases were used for a comprehensive protein sequence catalogue and to compare their differential abundance. Table 1 shows the number of identified or quantified proteins in the samples with reference to the Plant proteome (containing all plant entries in the SwissProt database), or reference proteomes for Populus trichocarpa and Arabidopsis thaliana contained in the Uniprot database, as well as to the proteins deduced from the first draft genome of Quercus suber from CorkOakDB (release 2018) [45].

thumbnail
Download:
Table 1. Number of proteins identified and quantified.

https://doi.org/10.1371/journal.pone.0245148.t001

The identification and quantification proteomics results are presented for the analyses using the Uniprot reference proteome for Arabidopsis (S2 Table) and for the predicted proteins from the cork oak genome (S3 Table). The later provided a probable annotation to 1,388 predicted proteins obtained by information-dependent acquisition (IDA) from pooled samples for each group. However, 58.8% of these proteins matched protein predictions of low confidence (containing the designations -like, -probable, -uncharacterized or -low quality protein) or corresponded to repeated entries among the CorkOakDB predicted proteins, revealing a high redundancy in this database. Consequently, the quantification of 841 predicted oak proteins was of low confidence, as the shared peptides were not able to be quantified under the quality criteria used for SWATH (S3 Table).

Given the robustness of the Arabidopsis protein database (reference proteome available at Uniprot.org, an highly curated protein database with low frequency of proteins of unknown function) and the availability of substantial functional annotation for gene ontologies and pathways, the Arabidopsis thaliana (considered a model organism for plants) was chosen as a reference for the following analyses. The exercise of inferring a biological meaning for proteins that stand out in the context of the interaction between Q. suber and P. cinnamomi it is only achievable taking as a reference a database with evidence-based functional annotation.

Thus, using the Arabidopis proteome database as reference, 424 proteins were confidently quantified in the cork oak leaves, with six biological replicates for each of the conditions control or inoculated (S2 Table). From these, 80 proteins showed a fold-change greater than 2 (or less than 0.5, in the case of proteins with decreased levels) or a p-value below 0.05 in their median levels between inoculated and control samples (Table 2). The Venn diagram in Fig 1 shows the number of proteins that met one or both criteria.

thumbnail
Download:
Fig 1. Protein groups.

The Venn diagram illustrates the number of proteins with a fold-change greater than 2 or less than 0.5 (yellow colour), those with a p-value below 0.05 (blue colour), and those that meet simultaneously both criteria based on protein levels between inoculated and control samples.

https://doi.org/10.1371/journal.pone.0245148.g001

thumbnail
Download:
Table 2. List of differentially accumulated proteins in Q. suber leaf proteome 8 months after P. cinnamomi inoculation, using the Arabidopis proteome database as a reference.

https://doi.org/10.1371/journal.pone.0245148.t002

Among the 80 proteins with differential levels, 60 proteins increased abundance, and 20 proteins decreased abundance in the leaves of inoculated cork oak plants, 8 months after P. cinnamomi inoculation, compared to the control plants (Table 2).

Hierarchical clustering of differentially produced cork oak proteins

The proteins with differential levels between inoculated and control samples were clustered in a heatmap to allow better visualization of the protein variation patterns (Fig 2). Inoculated plants are clearly distinguished from control plants based on the profiles of this protein dataset. In other words, eight months after a single inoculation of the cork oak root with P. cinnamomi, the inoculated but asymptomatic plants revealed a leaf proteome significantly different from the non-inoculated plants.

thumbnail
Download:
Fig 2. Hierarchical clustering of differentially produced cork oak proteins.

The heat map clusters the expression patterns of the 80 proteins with altered abundances between inoculated and control plants. Each column represents one cork oak plant; the control (C) plants are the first 6 columns on the left, and the 6 columns on the right are the inoculated (I) plants. The first 60 lines starting from the top of the heat map are proteins with an increased level in the inoculated samples (red color code) and the 20 lines towards the bottom, are proteins with a decreased level in the inoculated samples (green color code). The color scale of the heat map ranges from -3 to 3 (from light green to red).

https://doi.org/10.1371/journal.pone.0245148.g002

Two scenarios are possible for the inoculated plants: 1) the development of the oomycete was restricted to the inoculation site, with no spread of the infection; or 2) the development of the oomycete took place beyond the inoculation site, invading other tissues, but the infection still did not affect the vegetative state of the host. For the first hypothesis, the protein profiles observed in the leaves may be the result of the activation of the systemic defence system, maintained in memory over time. But, for the second hypothesis, the protein profiles of the inoculated plants may denote a homeostatic state of continuous interaction with P. cinnamomi.

The evaluation of cork oak infection by P. cinnamomi are always assessed at the root level, in a qualitative way, requiring an experienced technician for the identification of necrosis and/or absence of feeder roots. But, in this experiment, attempts were made to mimic field conditions, which are hampered by limitations regarding the detection and quantification of the oomycete in the rhizosphere of the trees. All the inoculated plants were used at the end of the experiment, and there was no selection based on the re-isolation of the oomycete or the existence of infection symptoms like leaf yellowing and wilting. Assessing cork oak decline in the field is based on the degree of canopy defoliation, and even if P. cinnamomi is isolated from the roots of declining trees, it is not possible to know the level or time of infection. Furthermore, the current methods used to isolate and identify P. cinnamomi from the rizosphere of oak roots are based on baiting tecnhiques, pathogen growth in selective media and molecular identification with specific primers. These procedures are time consuming, require expertise and are of relatively low effectiveness. Thus, evaluating cork oak decline through the leaf immune response protein profile induced by P. cinnamomi inoculation establishes a new approach for understanding the importance of this oomycete to cork oak decline.

Association of proteins to GO functional categories and biological pathways

Enrichment analysis was carried out on the selected dataset of 80 differential proteins, which assigned several GO terms to the proteins (based on the Arabidopsis thaliana proteins functional annotations), integrating them into the Gene Ontology functional categories of “Biological Process” (GO_BP) or “Cellular Component” (GO_CC). Within each classification, the significantly enriched terms (FDR<0.05; see lists in S4 and S5 Tables) were assembled into groups of functionally related terms by Cytoscape/ClueGO analysis and the most significantly enriched groups are presented in Fig 3 and summarized in S6 and S7 Tables.

thumbnail
Download:
Fig 3. Enrichment analysis applied to the subset of 80 differential proteins.

The bars represents the groups with higher enrichment score [−log2 (group FDR)] obtained for each group of functional related GO terms; the enrichment of the Gene Ontology category of GO_BP are showed in panel a), and those for GO_CC in panel b). Each group is labelled by the most significant (<FDR) enriched term, used as representative of the total enriched terms in each group that can be consulted in detail in S4 and S5 Tables.

https://doi.org/10.1371/journal.pone.0245148.g003

Thirty six GO Biological Process terms were significantly enriched among the differential proteins (S4 Table), which were grouped into 6 groups of functional related GO terms (summarized in Fig 3A); the group Pyridine-containing compound metabolic process has the highest enrichment score of all groups (12.9, corresponding to an FDR of 1.29 x 10−4), followed by Monosaccharide catabolic process (11.5), Cellular metabolic compound salvage (10.7), Protein-DNA complex assembly (10.1), Lipid oxidation (10.0) and Response to endoplasmic reticulum stress (10.0).

In the Cellular Component category, 15 GO terms were significantly enriched among the differential proteins (S5 Table), which were grouped into 4 groups of functional related GO terms (summarized in Fig 3B); the groups with the highest enrichment scores belonged to Cytosolic part (82.1), Stromule (22.8), Nucleosome (14.8) and Photosystem (14.1).

Furthermore, the enrichment analysis uncovered the most representative KEGG or REACTOME biological pathways in this dataset of the 80 differential proteins, which are listed in Table 3. Comparing the results from KEGG and REACTOME pathways, the protein subset is enriched in Ribosome and SRP-dependent cotranslational protein targeting to membrane, with the highest enrichment (lowest FDR) scores, respectively, and Glycolysis/Gluconeogenesis and Glucose metabolism with the lowest enrichment scores, respectively.

thumbnail
Download:
Table 3. Biological pathways uncovered for the selected dataset.

https://doi.org/10.1371/journal.pone.0245148.t003

The significantly enriched pathway with Reactome code R-ATH: 3371497- HSP90 chaperone cycle for steroid hormone receptors (SHR) is very relevant in the context of this investigation and of the available bibliography. The innate immunity and plant defence in Arabidopsis are biological events that are associated with the biological function of heat shock protein 90–2 as a molecular chaperone, involved in RPM1-mediated resistance and component of the RPM1/RAR1/SGT1 complex [51]. To circumvent the autoimmunity associated with high levels of immunity receptors, HSP90 proteins may assist in the formation of protein complexes that target the immune receptors SNC1, RPS2, and RPS4 for degradation [52].

Protein abundance patterns associated with GO functional categories

GO biological process category.

In the enrichment analysis, 18 proteins contributed significantly to certain biological processes (GO_BP groups of enriched GO_BP terms), which are detailed in Table 4. Of the six highlighted GO-BP groups, three stand out based on the constant patterns of variation in the abundance of the associated proteins. These are: Cellular metabolic compound salvage (GO group 3), with four down-accumulated proteins in the inoculated plants; Protein-DNA complex assembly (GO group 5), with four up-accumulated proteins in the inoculated plants and Response to endoplasmic reticulum stress (GO group 2), with two up-accumulated proteins in the inoculated plants. In the GO_BP groups Pyridine-containing compound metabolic process, Monosaccharide catabolic process and Lipid oxidation, proteins with different variation forms were housed in the same group.

thumbnail
Download:
Table 4. Proteins from the selected dataset assigned to Biological process GO terms with the highest enrichment scores.

https://doi.org/10.1371/journal.pone.0245148.t004

Within the proteins mapping to Cellular metabolic compound salvage (GO group 3), the protein Ribulose 1,5-bisphosphate carboxylase/oxygenase (RubisCO) small subunit 1A (P10795; RBCS1A) was the one showing the most expressive negative variation between inoculated and control samples (Log2FC = -3.1). Further, this protein showed a very high coefficient of variation (% CV), both in control (151%) and inoculated samples (152%). CV values may reflect the natural biological variability observed for RubisCO in Q. suber and the corresponding RubisCO patchiness in the host's response to the oomycete. In Quercus ilex inoculated with P. cinnamomi a decrease in the abundance of RubisCO proteins was also found, which was correlated with the tolerance/susceptibility of the provenances, being more accentuated in susceptible provenances [53].

RubisCO is very abundant in plants and the amount of this protein in the leaves is considered an indicator of the photosynthetic vigour and nitrogen availability. RBCS1A is a member of the multigene family RBCS from Arabidopsis and Isumy and colleagues [54] reported the additive effect of the expression of RBCS1A and RBCS3C genes on RubisCO accumulation in Arabidopsis leaves. Therefore, low levels of RBCS1A in the P. cinnamomi inoculated cork oak plants may indicate decreased levels of total RBCS mRNA and a smaller content of RubisCO accumulated in the leaves. It will be interesting to evaluate if the leaves of cork oak plants inoculated with P. cinnamomi have a reduced photosynthetic activity but, if this was a decrease, it appears not to have significantly affected plant growth that was similar between inoculated and control plants.

Like RubisCO, the chloroplastic phosphoribulokinase (P25697/PRK) is specifically associated with the Calvin-Benson cycle and catalyses D-ribulose 1,5-bisphosphate formation, used by RubisCO with CO2 or O2 to form 3-phosphoglycerate (3-PGA) and 2-phosphoglycolate (2-PG) [55,56]. Similarly, this protein also showed lower levels of accumulation in the P. cinnamomi inoculated plants in this study. The protein D-ribulose-5-phosphate-3-epimerase protein (Q9SAU2/RPE), mapping to GO groups 4 and 6 and participating in the carbon photoassimilation cycle, also showed lower levels of accumulation in the inoculated plants like RubisCO and PRK. It is likely that the decrease in the accumulation of these proteins can compromise the levels of carbon assimilated by the plant and the sequent synthesis of sugars, proteins, lipids and nucleic acids. A similar effect was observed in Q. ilex inoculated with P. cinnamomi where many proteins involved in the Calvin-Benson cycle, such as RubisCO large and small subunits, phosphoglycerate kinase, glyceraldehyde-3-phosphate dehydrogenase B and transketolase 1 were also decreased [53]. The complementarity of these data makes sense if we think about the existence of macromolecular complexes formed by phosphoribulokinase and glyceraldehyde-3-phosphate dehydrogenase interacting with the small peptide CP12, with relevance to the regulation of photosynthesis in the chloroplasts [57,58].

The two other proteins mapping to this GO group, (S)-2-hydroxy-acid oxidase, with the alternative name of glycolate oxidase 1 (Q9LRR9; GOX1/GLO1) and malate dehydrogenase (F4KDZ4; PMDH2), are known to be located in the peroxisomes, and both showed decreased levels in the inoculated plants. GOX1 catalyses the conversion of glycolate into glyoxylate with the production of H2O2 in the photorespiration pathway (EC 1.1.3.15). Modulation of hydrogen peroxide accumulation has been suggested as the mechanism adopted by the GOX protein family in Arabidopsis and Nicotiana benthamiana, associated with PAMP-triggered immunity (PTI), host and nonhost defence responses [59–61]. Furthermore, the defence pathways activated by different GOX genes vary between plant species and depend on the type of interaction that occurs between plants and pathogens or elicitors and GOX1-dependent defence responses may involve salicylic acid (SA) and WRKY62-mediated pathways [59–61].

Concerning malate dehydrogenase (F4KDZ4; PMDH2), in 2007, Pracharoenwattana and colleagues proposed a model in which the action of this enzyme is the production of malate from oxaloacetate with NADH oxidation, recruited for fatty acid ß-oxidation.

In summary, the decrease in the accumulation of these four proteins rebound on photosynthesis and concomitant photorespiration, and may affect sugar metabolism. Triacylglycerides oxidation may be used as an alternative source of energy and supply of gluconeogenic intermediates. Changes in the redox state of the cells resulting from the production of reactive oxygen species (ROS) are perceived and result in the activation of the defence system. Compromising the photosynthetic efficiency in source tissues may result in a reduced supply of sugars to sink tissues and less accumulation of soluble sugars.

Focusing on the Pyridine-containing compound metabolic process (GO group 4) and Monosaccharide catabolic process (GO group 6), hexose sugars like glucose are central molecules in plant metabolism and in sugar signaling. Cytosolic resources of phosphate hexoses originating from starch mobilization and sucrose hydrolysis are channeled for energy-producing and synthesis of biomolecule precursors [62]. Proteins PGMC1/PGM1 (Probable phosphoglucomutase), PMG1/iPGAM1 (2,3-bisphosphoglycerate-independent phosphoglycerate mutase 1), ALFC8/FBA8 (fructose-bisphosphate aldolase 8) and HXK1 (Hexokinase 1) from GO groups 4 and 6 can be grouped into two pairs according to their quantification pattern but also to the role they play in sugar metabolism and as sugar sensors. PGMC1/PGM1 is an enzyme that participates in both the breakdown and synthesis of glucose (EC: 5.4.2.2.) and PMG1/iPGAM1 is involved in the synthesis of pyruvate in glycolysis (EC: 5.4.2.12). These proteins were more abundant in the inoculated samples of this study, revealing a metabolic tendency in favour of energy production and reducing power as opposed to the accumulation of sucrose and carbohydrates as reserve substances. The availability of energy resulting from the functioning of enzymes may hamper energy-depending cell actions, such as the movement of stomata, and this requirement was studied in Arabidopsis through silencing the expression of glycolytic proteins. Silencing iPGAM activity in Arabidopsis is associated with reduced stomatal function and plant phenotypes with delayed development. This, probably results from the decrease in ATP production by the glycolytic pathway and also by tricarboxylic acid (TCA) cycle and oxidative phosphorylation in consequence of the concomitant reduction in the levels of pyruvate provided [63]. HXK1 and FBA8 were less abundant in the inoculated samples and both enzymes are involved in glycolysis (EC: 2.7.1.1; EC: 4.1.2.13) and sucrose metabolism, also being referred to as proteins involved in sugar and stress signaling [64,65]. In Arabidopsis, transcripts levels of AtFBA8 showed increased expression after 24h of glucose, fructose and sucrose treatment and these were responsive to ABA, SA, NaCl and drought stresses [65]. A reduction in the production of these enzymes in the inoculated cork oak plants may be a consequence of the imbalance of the metabolism towards the production of energy associated with the immune response. Knowing the subcellular location of HXK1 is essential to understand the role it plays in response to biotic stresses, because HXK1 located in the nucleus may interact with other proteins regulating the transcription of genes by binding directly to the chromatin and mitochondrial hexokinases can modulate programmed cell death (PCD) [66].

Cytosolic ICDHC/cICDH (Isocitrate dehydrogenase [NADP]) protein, mapping to GO group 4 and increasing its levels in inoculated plants, is potentially responsible for 2-oxoglutarate production for amino acid biosynthesis; however, in Arabidopsis cICDH is not required for plant development and primary metabolism in optimal growth conditions, instead, cICDH contributes to thiol–disulphide homeostasis during oxidative stress [67]. The NADPH produced by cICDH may contribute to activate defence responses to pathogen infection that are triggered by changes in cellular redox state [67].

Reprogramming gene expression in situations of biotic stress requires modulation of transcriptional activity in the nucleus. Focusing now on Protein-DNA complex assembly (GO group 5) and Response to endoplasmic reticulum stress (GO group 2), several proteins were found with increased levels, including two histones. It was suggested that appropriate levels of H3.3 are required to avoid H1 deposition over gene bodies preserving an adequate density of nucleosomes ideal for chromatin unfolding and access to DNA methyltransferases that methylate gene bodies [68]. Besides, the local enrichment of the nuclear histone H3.3 variant was positively correlated with transcription of responsive genes [69,70] and with gene body methylation [68]. Thus, the higher levels of H3.3 protein in the inoculated plants in this study suggest an increased access to DNA, allowing for modulation of transcription of biotic responsive-genes through gene body methylation. Still related to the formation of histone-DNA tetrasome is the Arabidopis chaperone NASP, described to bind to H3-H4 dimers and to stimulate the conversion of dimers to tetramers, in vitro [71]. Furthermore, in tobacco and Arabidopsis cell lines the modifications observed in histone H3 in response to abiotic stresses, with up-regulation of marker genes, happens together with histone H4 acetylation, revealing the parallel intervention of these histones [72]. Also, in the present study, a significant up-accumulation of both H3.3 and H4 proteins were found in the leaves of cork oak plants inoculated with P. cinnamomi.

To reach cellular homeostasis, protein degradation by proteolytic enzymes is a regulated procedure used to adjust protein abundance and efficiency. This activity in the chloroplasts requires Clp protease complexes composed by several protein catalytic (ClpP3 to ClpP6) and non-catalytic (ClpR1 to ClpR4) subunits arranged in ring-like structures (P-ring and R-ring) in Arabidopsis surrounding the proteolytic chamber whose activity is assisted by several chaperone members [for review see 73]. The assembly of the rings of the Clp core complex is compromised if there is an uneven number of subunits available for its formation. It was observed that reducing the abundance of the subunit ClpP6 by 50% caused a reduction in the protein abundance of other P- and R-ring components, interfering with the complex assembly and functionality [74]. Therefore, it is reasonable to expect that the increase in accumulation of the ClpR4 subunit in the inoculated cork oak samples, may point to the importance of protein degradation in the regulation of the photosynthetic process mediated by the Clp complex. Additionally, knockdown of protease subunits in tobacco allowed the identification of putative protease substrates, including proteins involved in photosynthesis like PRK and RPE, which were found to be down-accumulated in the cork oak samples inoculated with P. cinnamomi (Table 4).

Beyond the Clp complex, there are other biological mechanisms that predict protein degradation, through proteasome complexes, to modulate the activity of disease resistance proteins (R) in plant-pathogen interactions and also other processes such as the oxidative burst, hormone signaling, gene induction, and programmed cell death [75]. In tobacco cells challenged by the elicitin cryptogein, the accumulation of 20S proteasome subunits was observed simultaneously with the development of systemic acquired resistance [76]. By analogy, the RPT4A protein highlighted in the leaf proteome of the inoculated samples in the present study can be regarded as a defence-induced subunit of 26S proteasome, with a possible role in plant defence reactions eventually triggered by elicitins produced by P. cinnamomi. It is also possible to assume ATPase activity for Q. suber RTP4A subunit based on the functional characterization of the 26S proteasomal subunit RPT4a from Solanum lycopersicum that has an active ATPase site and can modulate the resistance to the ToLCNDV virus by physically interacting with viral DNA molecules [77].

Modifying the programming of host nuclear gene transcription in response to biotic stress is one of the mechanisms adopted by oomycetes and promoted through effector molecules. The target genes may be those that code for HSP (Heat Shock Proteins) proteins, modulating the role played by these molecular chaperones; these are active partners of numerous enzyme complexes and are responsible for the folding and unfolding of proteins included in protein degradation/renaturation and movement of signaling proteins and transcription factors into cell organelles [78,79]. In Q. suber leaves inoculated with P. cinnamomi, the accumulation of HSP70-5 (Q9S9N1) and HSP90-1 (P27323) proteins was higher than in control samples. Knowing the activities of these proteins in other plant species it can be inferred the possible role they may play in the interaction between cork oak and P. cinnamomi. Song et al (2015) reported the identification of a P. sojae intracellular CRN (Crinkler or crinkling- and necrosis-inducing protein) effector which directly interacts with promoters of the genes encoding HSP proteins, preventing the binding of specific transcription factors [80]. The expression of the defence-related genes in Arabidopsis, N. benthamiana and soybean is then changed, unbalancing the host's resistance level to Phytophthora species [80]. More recently, HSP70s have been noted as proteins that interact with RXLR effectors produced by P. infestans and that get involved in N. benthamiana defence response by stimulating programmed cell death mediated by MAPK signaling and suppressing the growth of the pathogen [81]. It was also reported the importance of protein complexes formed between HSP90s and co-chaperones in the activation of defence mechanisms mediated by resistance (R) proteins after the detection of pathogen effector molecules [82]. Preventing the formation of these protein complexes by inhibiting the binding of HSP90 has implications for the accumulation of R proteins and resistance mediated by these proteins [86]. In Arabidopsis, the cytosolic AtHSP90.1 was the only HSP90 isoform significantly induced after inoculation with Pseudomonas (Pst) strains containing avirulence genes (avrRpm1 and avrRpt2) and was required for the full resistance mediated by one of the corresponding R proteins [83]. It is then expected that the host's resistance proteins will recognize the effector molecules secreted by P. cinnamomi and activate the defence responses with the collaboration of chaperones and co-chaperones.

The disclosed Q. suber HSP70-5 is recognized as a homolog of the Arabidopsis BIP1 for Binding Immunoglobulin Protein or Binding Protein (BiP), a chaperone set in the endoplasmic reticulum (ER) lumen, known to bind a membrane-associated transcription factor (TF) under non-stressed conditions [84]. The interaction between BiP1 and the TF is the requirement to turn on or off a protein secretory signaling pathway via ER-Golgi-Nucleus, ending with the transcription of stress response genes [84]. In soybean, the BiP protein was described as a negative regulator of a stress-induced cell death response and, in Arabidopsis, Wang et al. (2005) reported the implications of BiP 2 silencing on the secretion of pathogenesis related-proteins, compromising the systemic acquired resistance against bacterial pathogens [85].

When looking at the enrichment of Lipid oxidation (GO group 1), two differential proteins were assigned to it, lipoxygenase 2 (LOX2; P38418) and thiolase family protein (AACT1; F4JYM8), although with different patterns of variation: LOX2 was more abundant in the inoculated plants and AACT1 was less abundant. By homology to the Arabidopsis ACCT1 isoform, it is expected for Q. suber ACT1 to be located in the peroxisome, based on the presence of two alternative targeting sequences PTS1 and PTS2 motives found in AtACCT1, excluding a metabolic function related to isoprenoid biosynthesis [86,87]. Jin et al (2012) [88] found a strong expression of AtACCT1 in the vascular system of the aerial organs and roots of Arabidopsis, verifying that gene silencing or induction of abiotic stresses did not result in an evident phenotypic response [88]. In Q. suber inoculated by P. cinnamomi, the leaf proteome reveals a reduction in the production of a thiolase AACT1 protein, apparently included in its defence strategy.

The involvement of LOX2 in the cork oak defence response may be associated with the production of jasmonic acid (JA) via the Vick and Zimmerman pathway [89]. Recently, it was confirmed that LOX2 forms a protein complex with AOS (allene oxide synthase) and AOC2 (allene oxide cyclase), two proteins that also participate in JA precursor biosynthesis, located in the inner envelope of the Arabidopsis chloroplasts [90]. The formation of this molecular complex is evident in the effectiveness of JA production to the disadvantage of other products resulting from parallel reactions during oxylipin biosynthesis, guiding the defence response to the activation of genes responsive to JA [90]. Sometimes, depending on the needs of the pathogen, the signaling reactions are a balance between the activation of the salicylic acid (SA) pathway with suppression of the jasmonic acid (JA) pathway or vice versa [91]. Nevertheless, Starý et al. (2019) conclude that the level of resistance induced in different tomato genotypes after β-cryptogein treatment correlated with the upregulation of defence genes and activated ethylene and JA signaling but not SA signaling [92]. Other authors refer to a biphasic defence response in avocado against the hemibiotroph P. cinnamomi, which initially involves SA-mediated gene expression followed by the enrichment of JA-mediated defence from 18 to 24 hours post-inoculation [93].

Finally, when analysing the GO Cellular Component category, the highest enrichment scores were obtained for Cytosolic part (82.1), Stromule (22.8), Nucleosome (14.8) and Photosystem (14.1). In a simplified view, it is recognized that achieving new cell balances during interaction with pathogenic organisms requires the remodeling of physiological processes by the action of cytoplasmic or organelle-associated enzymes and the modulation of transcription factors for nuclear gene expression, which are energy-dependent processes. Three of the obtained categories fit this profile with the exception of Stromule, which is a novelty for the host-P. cinnamomi interactions. Nevertheless, there were previous descriptions for the importance of the communication between cellular organelles during immune responses carried out through stroma-filled tubular structures (stromules) of the chloroplasts-to-nucleus, which use them as a support for the exchange of molecules integrated into the defence response [94,95]. Caplan et al (2015) observed the induction of stromules in response to viral and bacterial effectors after recognition by host receptors (ETI; effector-triggered immunity) [94]. In plants infected with tobacco mosaic virus (TMV) a hypersensitive response is observed at the infection site and in the border regions with increased production of stromules in both areas, probably stimulated by the production of pro-defence signaling molecules like H2O2, O2- and SA [95]. During innate immunity, the cellular relocation of chloroplasts in the nucleus surroundings is dependent on the organization of microtubules in connection with the anchoring points provided by actin filaments to enhance the effectiveness of the communication between these organelles [96]. In 2013, Sghaier-Hammami et al. reported the up-accumulation of actin in holm oak plants inoculated with P. cinnamomi [53]. The published information reinforces the importance of the Stromules GO category highlighted in the present study for the cork oak-P. cinnamomi interaction. Moreover, the cork oak leaf proteome data suggests a possible function for stromules in long-term defence responses, far from the inoculation point, in close connection with the production and transport of signaling molecules.

Conclusions

In this work, the proteomes of cork oaks plants submitted to biotic stress-induced by P. cinnamomi inoculation are revealed for the first time. Among the 424 proteins confidently quantified in the inoculated and non-inoculated plants, a dataset of 80 proteins was selected based on the abundance variability observed between the experimental conditions. The immune response of the plants was analysed eight months after the inoculation event, and, at that moment, there were no evident phenotypic differences between inoculated and non-inoculated plants. Nevertheless, the hierarchical clustering of differentially produced cork oak proteins shows two different groups of plants, matching to the experimental conditions. By comparing protein profiles, it was observed that the number of proteins in which the abundance increased in the inoculated plants is 3 times greater than the number of proteins in which there was a decrease in abundance. Therefore, the defence responses induced in the host by P. cinnamomi in the long term and distant from the inoculation site are inscribed in the proteome of the leaves, reproducing the in progress homeostatic state of the plants. The results obtained in this study increase the possibilities of screening trees infected with P. cinnamomi using protein markers identified in the leaves without the need to isolate the oomycete from the roots of the host or surrounding soil.

The homeostatic state of the inoculated cork oak plants was characterized by protein patterns associated with differential biological processes occurring potentially in different subcellular organelles. When performing the enrichment analysis, eighteen proteins were highlighted, and their possible functions in an immune response context were discussed. In short, the decrease in the accumulation of photosynthesis enzymes and concomitant photorespiration may compromise the levels of carbon assimilated by the plant and its development, although throughout the experiment, no differences in growth were observed. The dynamics of proteins associated with sugar metabolism and sugar signaling reveals a metabolic tendency in favour of energy production and reducing power as opposed to the accumulation of sucrose and carbohydrates as reserve substances.

The reprogramming of gene expression, eventually in response to the action of effector molecules produced by P. cinnamomi is a major function associated with proteins that are in greater abundance in the inoculated plants. It is also clear the participation of proteolytic complexes and chaperones in the cork oak immune response and of proteins sensitive to changes in the redox state of the cell promoted by ROS species.

In addition, the cork oak leaf proteome data suggests the importance of the communication between cellular organelles mediated by stromules in the long-term defence responses. Immune response amplification and effectiveness may be dependent on the repositioning of the chloroplasts close to the nucleus and the transfer of pro-defence molecules such as SA, JA and H2O2.

Supporting information

S1 Fig. Biological material.

Information on biological material and procedures performed in the experimental assay.

https://doi.org/10.1371/journal.pone.0245148.s001

(TIF)

S2 Fig. Visual observation of the plants over the duration of the experiment.

https://doi.org/10.1371/journal.pone.0245148.s002

(TIF)

S1 Table. Cork oak references.

Plant labels and GPS references for location of cork oak parental trees.

https://doi.org/10.1371/journal.pone.0245148.s003

(PDF)

S2 Table. Proteins identified and quantified in cork oak leaves compared to the Arabidopsis thaliana reference proteome.

List of the 608 proteins identified by IDA analysis (Table S2.1.) and list of the 424 proteins quantified by SWATH-MS in control (C) or in inoculated (I) cork oak leaves (Table S2.2.).

https://doi.org/10.1371/journal.pone.0245148.s004

(XLSX)

S3 Table. Proteins identified and quantified in cork oak leaves compared to the proteins deduced from the draft genome of cork oak.

List of the 1388 proteins identified by IDA analysis (Table S3.1.) and list of the 841 proteins quantified by SWATH-MS in control (C) or in inoculated (I) cork oak leaves (Table S3.2.).

https://doi.org/10.1371/journal.pone.0245148.s005

(XLSX)

S4 Table. Significantly enriched GO biological process terms and groups in the list of 80 differential proteins.

https://doi.org/10.1371/journal.pone.0245148.s006

(PDF)

S5 Table. Significantly enriched GO cellular component terms and groups in the list of 80 differential proteins.

https://doi.org/10.1371/journal.pone.0245148.s007

(PDF)

S6 Table. Significantly enriched GO biological process groups in the list of 80 differential proteins.

https://doi.org/10.1371/journal.pone.0245148.s008

(PDF)

S7 Table. Significantly enriched GO cellular component groups in the list of 80 differential proteins.

https://doi.org/10.1371/journal.pone.0245148.s009

(PDF)

S1 File. Description of the SWATH-MS.

Principles and detailed materials and methods.

https://doi.org/10.1371/journal.pone.0245148.s010

(PDF)

References

  1. 1. Serrano MS, Rios P, Gonzalez M, Sanchez ME. Experimental minimum threshold for Phytophthora cinnamomi root disease expression on Quercus suber. Phytopathol Mediterr. 2015:461–4.
  2. 2. Camilo-Alves CDEP, da Clara MIE, Ribeiro NMCD. Decline of Mediterranean oak trees and its association with Phytophthora cinnamomi: a review. Eur J Forest Res. 2013;132(3):411–32.
  3. 3. Vessella F, López-Tirado J, Simeone MC, Schirone B, Hidalgo PJ. A tree species range in the face of climate change: cork oak as a study case for the Mediterranean biome. Eur J For Res. 2017;136(3):555–69.
  4. 4. Camilo-Alves CS, Vaz M, Da Clara MIE, Ribeiro NMDA. Chronic cork oak decline and water status: new insights. New For. 2017;48(6):753–72.
  5. 5. Avila J, Linares J, García-Nogales A, Sánchez M, Gómez-Aparicio L. Across-scale patterning of plant–soil–pathogen interactions in Quercus suber decline. Eur J For Res. 2017;136(4):677–88.
  6. 6. Mendes MP, Ribeiro L, David TS, Costa A. How dependent are cork oak (Quercus suber L.) woodlands on groundwater? A case study in southwestern Portugal. Forest Ecol Manag. 2016;378:122–30.
  7. 7. Tiberi R, Branco M, Bracalini M, Croci F, Panzavolta T. Cork oak pests: a review of insect damage and management. Ann For Sci. 2016;73(2):219–32.
  8. 8. Henriques J, Barrento MJ, Bonifácio L, Gomes AA, Lima A, Sousa E. Factors affecting the dispersion of Biscogniauxia mediterranea in Portuguese cork oak stands. Silva Lusit. 2014;22(1):83–97.
  9. 9. Moricca S, Linaldeddu BT, Ginetti B, Scanu B, Franceschini A, Ragazzi A. Endemic and emerging pathogens threatening cork oak trees: management options for conserving a unique forest ecosystem. Plant Dis. 2016;100(11):2184–93. pmid:30682920
  10. 10. Arosa ML, Bastos R, Cabral JA, Freitas H, Costa SR, Santos M. Long-term sustainability of cork oak agro-forests in the Iberian Peninsula: model-based approach aimed at supporting the best management options for the montado conservation. Ecol Modell. 2017;343:68–79.
  11. 11. Horta M, Caetano P, Medeira C, Maia I, Cravador A. Involvement of the β-cinnamomin elicitin in infection and colonisation of cork oak roots by Phytophthora cinnamomi. Eur J Plant Pathol. 2010;127(3):427–36.
  12. 12. Horta M, Sousa N, Coelho AC, Neves D, Cravador A. In vitro and in vivo quantification of elicitin expression in Phytophthora cinnamomi. Physiol Mol Plant Pathol. 2008;73(1–3):48–57.
  13. 13. Islam M, Hussain H, Russo R, Chambery A, Amoresano A, Schallmey A, et al. Functional analysis of elicitins and identification of cell wall proteins in Phytophthora cinnamomi. Physiol Mol Plant Pathol. 2019;107:21–32.
  14. 14. Hardham AR, Blackman LM. Phytophthora cinnamomi. Mol Plant Pathol. 2018;19(2):260–85. pmid:28519717
  15. 15. Kale S, Gu B, Capelluto D, Dou D, Feldman E, Rumore A, et al. External phosphatidylinositol-3-phosphate mediates host cell entry by eukaryotic pathogen effectors. Cell. 2010;142(2):284–95. pmid:20655469
  16. 16. Kale SD, Tyler BM. Entry of oomycete and fungal effectors into plant and animal host cells. Cell Microbiol. 2011;13(12):1839–48. pmid:21819515
  17. 17. Wang H, Ren Y, Zhou J, Du J, Hou J, Jiang R, et al. The cell death triggered by the nuclear localized RxLR effector PITG_22798 from Phytophthora infestans is suppressed by the effector AVR3b. Int J Mol Sci. 2017;18(2):409. pmid:28216607
  18. 18. Gascuel Q, Buendia L, Pecrix Y, Blanchet N, Muños S, Vear F, et al. RXLR and CRN effectors from the sunflower downy mildew pathogen Plasmopara halstedii induce hypersensitive-like responses in resistant sunflower lines. Front Plant Sci. 2016;7:1887. pmid:28066456
  19. 19. Kong L, Qiu X, Kang J, Wang Y, Chen H, Huang J, et al. A Phytophthora effector manipulates host histone acetylation and reprograms defense gene expression to promote infection. Current Biol. 2017;27(7):981–91. pmid:28318979
  20. 20. Jones J, Dangl. J. The plant immune system. Nature. 2006;444:323–9. pmid:17108957
  21. 21. Hardoim P, Guerra R, Rosa da Costa A, Serrano M, Sánchez M, Coelho A. Temporal metabolic profiling of the Quercus suber–Phytophthora cinnamomi system by middle‐infrared spectroscopy. For Pathol. 2016;46(2):122–33.
  22. 22. Engelbrecht J, Van den Berg N. Expression of defence-related genes against Phytophthora cinnamomi in five avocado rootstocks. S Afr J Sci. 2013;109(11–12):1–8.
  23. 23. Meyer FE, Shuey LS, Naidoo S, Mamni T, Berger DK, Myburg AA, et al. Dual RNA-sequencing of Eucalyptus nitens during Phytophthora cinnamomi challenge reveals pathogen and host factors influencing compatibility. Front Plant Sci. 2016;7:191. pmid:26973660
  24. 24. Serrazina S, Santos C, Machado H, Pesquita C, Vicentini R, Pais MS, et al. Castanea root transcriptome in response to Phytophthora cinnamomi challenge. Tree Genet Genomes. 2015;11(1):6.
  25. 25. Coelho AC, Horta M, Ebadzad G, Cravador A. Quercus suber–P. cinnamomi interaction: hypothetical molecular mechanism model. NZ J Forestry Sci. 2011;(41S):143–57.
  26. 26. Flor HH. Current status of the gene-for-gene concept. Annu Rev Phytopathol. 1971;9(1):275–96.
  27. 27. van der Hoorn RA, Kamoun S. From guard to decoy: a new model for perception of plant pathogen effectors. Plant Cell. 2008;20(8):2009–17. pmid:18723576
  28. 28. Santos C, Duarte S, Tedesco S, Fevereiro P, Costa RL. Expression profiling of Castanea genes during resistant and susceptible interactions with the oomycete pathogen Phytophthora cinnamomi reveal possible mechanisms of immunity. Front Plant Sci. 2017;8:515. pmid:28443110
  29. 29. Dodds PN, Rathjen JP. Plant immunity: towards an integrated view of plant-pathogen interactions. Nat Rev Genet. 2010;11(8):53948. pmid:20585331
  30. 30. Fu ZQ, Dong X. Systemic acquired resistance: turning local infection into global defense. Annu Rev Plant Biol. 2013;64:839–63. pmid:23373699
  31. 31. Gao Q-M, Zhu S, Kachroo P, Kachroo A. Signal regulators of systemic acquired resistance. Front Plant Sci. 2015;6:228. pmid:25918514
  32. 32. Tada Y, Spoel SH, Pajerowska-Mukhtar K, Mou Z, Song J, Wang C, et al. Plant immunity requires conformational changes of NPR1 via S-nitrosylation and thioredoxins. Science. 2008;321(5891):952–6. pmid:18635760
  33. 33. Zhou J-M, Trifa Y, Silva H, Pontier D, Lam E, Shah J, et al. NPR1 differentially interacts with members of the TGA/OBF family of transcription factors that bind an element of the PR-1 gene required for induction by salicylic acid. Mol Plant Microbe Interact. 2000;13(2):191–202. pmid:10659709
  34. 34. van den Berg N, Christie J, Aveling T, Engelbrecht J. Callose and β‐1, 3‐glucanase inhibit Phytophthora cinnamomi in a resistant avocado rootstock. Plant Pathol. 2018;67(5):1150–60.
  35. 35. Aizat WM, Ibrahim S, Rahnamaie-Tajadod R, Loke K-K, Goh H-H, Noor NM. Proteomics (SWATH-MS) informed by transcriptomics approach of tropical herb Persicaria minor leaves upon methyl jasmonate elicitation. PeerJ. 2018;6:e5525. pmid:30186693
  36. 36. Anjo SI, Santa C, Manadas B. SWATH-MS as a tool for biomarker discovery: From basic research to clinical applications. Proteomics. 2017;17(3–4). pmid:28127880
  37. 37. Zhang H, He D, Yu J, Li M, Damaris RN, Gupta R, et al. Analysis of dynamic protein carbonylation in rice embryo during germination through AP-SWATH. Proteomics. 2016;16(6):989–1000. pmid:26801057
  38. 38. Zhu F-Y, Chan W-L, Chen M-X, Kong RPW, Cai C, Wang Q, et al. SWATH-MS Quantitative Proteomic Investigation Reveals a Role of Jasmonic Acid during Lead Response in Arabidopsis. J Proteome Res. 2016;15(10):3528–39. pmid:27599093
  39. 39. Zhu F-Y, Chen M-X, Chan W-L, Yang F, Tian Y, Song T, et al. SWATH-MS quantitative proteomic investigation of nitrogen starvation in Arabidopsis reveals new aspects of plant nitrogen stress responses. J Proteomics. 2018;187:161–70. pmid:30048775
  40. 40. Chen MX, Zhang Y, Fernie AR, Liu YG, Zhu FY. SWATH-MS-Based Proteomics: Strategies and Applications in Plants. Trends Biotech. 2020. Epub 2020/10/11. pmid:33036785
  41. 41. Zhu FY, Song YC, Zhang KL, Chen X, Chen MX. Quantifying Plant Dynamic Proteomes by SWATH-based Mass Spectrometry. Trends Plant Sci. 2020;25(11):1171–2. pmid:32891562
  42. 42. Coelho AC, Cravador A, Bollen A, Ferraz J, Moreira A, Fauconnier A, et al. Highly specific and sensitive non-radioactive molecular identification of Phytophthora cinnamomi. Mycol Res. 1997;101(12):1499–507.
  43. 43. Bradford MM. Rapid and Sensitive Method for Quantitation of Microgram Quantities of Protein Utilizing Principle of Protein-Dye Binding. Anal Biochem. 1976;72(1–2):248–54. pmid:942051
  44. 44. Anjo SI, Santa C, Manadas B. Short GeLC‐SWATH: A fast and reliable quantitative approach for proteomic screenings. Proteomics. 2015;15(4):757–62. pmid:25418953
  45. 45. Ramos A, Usie A, Barbosa P, Barros P, Capote T, Chaves I. The draft genome sequence of cork oak. Sci Data. 2018;5:180069. pmid:29786699
  46. 46. Pinto PI, Andrade AR, Moreira C, Zapater C, Thorne MA, Santos S, et al. Genistein and estradiol have common and specific impacts on the sea bass (Dicentrarchus labrax) skin-scale barrier. J Steroid Biochem Mol Biol. 2019;195:105448. pmid:31421232
  47. 47. Pinto PI, Andrade AR, Thorne MA, Estêvão MD, Canario AV, Power DM. Estradiol and genistein effects on the sea bass (Dicentrarchus labrax) scales: Transcriptome dataset. Data Brief. 2019;27:104587. pmid:31763380
  48. 48. Isserlin R, Merico D, Voisin V, Bader GD. Enrichment Map–a Cytoscape app to visualize and explore OMICs pathway enrichment results. F1000 Res. 2014;3.
  49. 49. Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, et al. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. 2009;25(8):1091–3. pmid:19237447
  50. 50. Coelho AC, Lima M, Neves D, Cravador A. Genetic diversity of two evergreen oaks [Quercus suber (L.) and Quercus ilex subsp. rotundifolia (Lam.)] in Portugal using AFLP markers. Silvae Genet. 2006;55(1–6):105–18.
  51. 51. Hubert DA, He Y, McNulty BC, Tornero P, Dangl JL. Specific Arabidopsis HSP90. 2 alleles recapitulate RAR1 cochaperone function in plant NB-LRR disease resistance protein regulation. Proc Natl Acad Sci U S A. 2009;106(24):9556–63. pmid:19487680
  52. 52. Huang S, Monaghan J, Zhong X, Lin L, Sun T, Dong OX, et al. HSP 90s are required for NLR immune receptor accumulation in Arabidopsis. Plant J. 2014;79(3):427–39. pmid:24889324
  53. 53. Sghaier-Hammami B, Valero-Galvàn J, Romero-Rodríguez MC, Navarro-Cerrillo RM, Abdelly C, Jorrín-Novo J. Physiological and proteomics analyses of Holm oak (Quercus ilex subsp. ballota [Desf.] Samp.) responses to Phytophthora cinnamomi. Plant Physiol Biochem. 2013;71:191–202. pmid:23962806
  54. 54. Izumi M, Tsunoda H, Suzuki Y, Makino A, Ishida H. RBCS1A and RBCS3B, two major members within the Arabidopsis RBCS multigene family, function to yield sufficient Rubisco content for leaf photosynthetic capacity. J Exp Bot. 2012;63(5):2159–70. pmid:22223809
  55. 55. Kleffmann T, Russenberger D, von Zychlinski A, Christopher W, Sjölander K, Gruissem W, et al. The Arabidopsis thaliana chloroplast proteome reveals pathway abundance and novel protein functions. Curr Biol. 2004;14(5):354–62. pmid:15028209
  56. 56. Maurino VG, Engqvist MK. 2-Hydroxy acids in plant metabolism. The Arabidopsis book: American Society of Plant Biologists; 2015. pmid:26380567
  57. 57. Marri L, Sparla F, Pupillo P, Trost P. Co-ordinated gene expression of photosynthetic glyceraldehyde-3-phosphate dehydrogenase, phosphoribulokinase, and CP12 in Arabidopsis thaliana. J Exp Bot. 2005;56(409):73–80. pmid:15533878
  58. 58. Graciet E, Lebreton S, Gontero B. Emergence of new regulatory mechanisms in the Benson–Calvin pathway via protein–protein interactions: a glyceraldehyde‐3‐phosphate dehydrogenase/CP12/phosphoribulokinase complex. J Exp Bot. 2004;55(400):1245–54. pmid:15047759
  59. 59. Rojas CM, Senthil-Kumar M, Wang K, Ryu CM, Kaundal A, Mysore KS. Glycolate oxidase modulates reactive oxygen species-mediated signal transduction during nonhost resistance in Nicotiana benthamiana and Arabidopsis. Plant Cell. 2012;24(1):336–52. pmid:22286136
  60. 60. Rojas CM, Mysore KS. Glycolate oxidase is an alternative source for H2O2 production during plant defense responses and functions independently from NADPH oxidase. Plant Signal Behav. 2012;7(7):752–5. pmid:22751316
  61. 61. Xu Y-P, Yang J, Cai X-Z. Glycolate oxidase gene family in Nicotiana benthamiana: genome-wide identification and functional analyses in disease resistance. Sci Rep. 2018;8(1):1–11. pmid:29311619
  62. 62. Everard J, Loescher W. Primary Products of Photosynthesis, Sucrose and Other Soluble Carbohydrates. In: Thomas B, Murray BG, Murphy DJ, editors. Encyclopedia of Applied Plant Sciences. 1: Waltmam, MA: Academic Press; 2017. p. 96–104.
    • 63. Zhao Z, Assmann SM. The glycolytic enzyme, phosphoglycerate mutase, has critical roles in stomatal movement, vegetative growth, and pollen production in Arabidopsis thaliana. J Exp Bot. 2011;62(14):5179–89. pmid:21813794
    • 64. Rosa M, Prado C, Podazza G, Interdonato R, González JA, Hilal M, et al. Soluble sugars: Metabolism, sensing and abiotic stress: A complex network in the life of plants. Plant Signal Behav. 2009;4(5):388–93. pmid:19816104
    • 65. Lu W, Tang X, Huo Y, Xu R, Qi S, Huang J, et al. Identification and characterization of fructose 1, 6-bisphosphate aldolase genes in Arabidopsis reveal a gene family with diverse responses to abiotic stresses. Gene. 2012;503(1):65–74. pmid:22561114
    • 66. Kim M, Lim J-H, Ahn CS, Park K, Kim GT, Kim WT, et al. Mitochondria-associated hexokinases play a role in the control of programmed cell death in Nicotiana benthamiana. Plant Cell. 2006;18(9):2341–55. pmid:16920781
    • 67. Mhamdi A, Mauve C, Gouia H, Saindrenan P, Hodges M, Noctor G. Cytosolic NADP‐dependent isocitrate dehydrogenase contributes to redox homeostasis and the regulation of pathogen responses in Arabidopsis leaves. Plant Cell Environ. 2010;33(7):1112–23. pmid:20199623
    • 68. Wollmann H, Stroud H, Yelagandula R, Tarutani Y, Jiang D, Jing L, et al. The histone H3 variant H3. 3 regulates gene body DNA methylation in Arabidopsis thaliana. Genome Biol 2017;18(1):94. pmid:28521766
    • 69. Nie X, Wang H, Li J, Holec S, Berger F. The HIRA complex that deposits the histone H3. 3 is conserved in Arabidopsis and facilitates transcriptional dynamics. Biol Open. 2014;3(9):794–802. pmid:25086063
    • 70. Wollmann H, Holec S, Alden K, Clarke ND, Jacques P-E, Berger F. Dynamic deposition of histone variant H3. 3 accompanies developmental remodeling of the Arabidopsis transcriptome. PLoS Genet. 2012;8(5):e1002658. pmid:22570629
    • 71. Maksimov V, Nakamura M, Wildhaber T, Nanni P, Ramström M, Bergquist J, et al. The H3 chaperone function of NASP is conserved in Arabidopsis. Plant J. 2016;88(3):425–36. pmid:27402088
    • 72. Sokol A, Kwiatkowska A, Jerzmanowski A, Prymakowska-Bosak M. Up-regulation of stress-inducible genes in tobacco and Arabidopsis cells in response to abiotic stresses and ABA treatment correlates with dynamic changes in histone H3 and H4 modifications. Planta. 2007;227(1):245–54. pmid:17721787
    • 73. Clarke AK. The chloroplast ATP‐dependent Clp protease in vascular plants–new dimensions and future challenges. Physiol Plant. 2012;145(1):235–44. pmid:22085372
    • 74. Moreno JC, Martínez-Jaime S, Schwartzmann J, Karcher D, Tillich M, Graf A, et al. Temporal proteomics of inducible RNAi lines of Clp protease subunits identifies putative protease substrates. Plant Physiol. 2018;176(2):1485–508. pmid:29229697
    • 75. Trujillo M, Shirasu K. Ubiquitination in plant immunity. Current opinion in plant biology. 2010;13(4):402–8. pmid:20471305
    • 76. Suty L, Lequeu J, Lançon A, Etienne P, Petitot A-S, Blein J-P. Preferential induction of 20S proteasome subunits during elicitation of plant defense reactions: towards the characterization of “plant defense proteasomes”. Int J Biochem Cell B. 2003;35(5):637–50.
    • 77. Sahu PP, Sharma N, Puranik S, Chakraborty S, Prasad M. Tomato 26S Proteasome subunit RPT4a regulates ToLCNDV transcription and activates hypersensitive response in tomato. Sci Rep. 2016;6(1):1–12. pmid:28442746
    • 78. Liu J, Pang X, Cheng Y, Yin Y, Zhang Q, Su W, et al. The Hsp70 gene family in Solanum tuberosum: genome-wide identification, phylogeny, and expression patterns. Sci Rep. 2018;8(1):1–11. pmid:29311619
    • 79. Cha J-Y, Ahn G, Kim JY, Kang SB, Kim MR, Su'udi M, et al. Structural and functional differences of cytosolic 90-kDa heat-shock proteins (Hsp90s) in Arabidopsis thaliana. Plant Physiol Biochem. 2013;70:368–73. pmid:23827697
    • 80. Song T, Ma Z, Shen D, Li Q, Li W, Su L, et al. An oomycete CRN effector reprograms expression of plant HSP genes by targeting their promoters. PLoS Pathog. 2015;11(12):e1005348. pmid:26714171
    • 81. Lee JH, Lee SE, Oh S, Seo E, Choi D. HSP70s Enhance a Phytophthora infestans Effector-Induced Cell Death via an MAPK Cascade in Nicotiana benthamiana. Mol Plant Microbe Interact. 2018;31(3):356–62. pmid:29140163
    • 82. Botër M, Amigues B, Peart J, Breuer C, Kadota Y, Casais C, et al. Structural and functional analysis of SGT1 reveals that its interaction with HSP90 is required for the accumulation of Rx, an R protein involved in plant immunity. Plant Cell. 2007;19(11):3791–804. pmid:18032631
    • 83. Takahashi A, Casais C, Ichimura K, Shirasu K. HSP90 interacts with RAR1 and SGT1 and is essential for RPS2-mediated disease resistance in Arabidopsis. Proc Natl Acad Sci U S A. 2003;100(20):11777–82. pmid:14504384
    • 84. Srivastava R, Deng Y, Shah S, Rao AG, Howell SH. BINDING PROTEIN is a master regulator of the endoplasmic reticulum stress sensor/transducer bZIP28 in Arabidopsis. Plant Cell. 2013;25(4):1416–29. pmid:23624714
    • 85. Wang D, Weaver ND, Kesarwani M, Dong X. Induction of protein secretory pathway is required for systemic acquired resistance. Science. 2005;308(5724):1036–40. pmid:15890886
    • 86. Ahumada I, Cairó A, Hemmerlin A, González V, Pateraki I, Bach TJ, et al. Characterisation of the gene family encoding acetoacetyl-CoA thiolase in Arabidopsis. Funct Plant Biol. 2008;35(11):1100–11. pmid:32688858
    • 87. Carrie C, Murcha MW, Millar AH, Smith SM, Whelan J. Nine 3-ketoacyl-CoA thiolases (KATs) and acetoacetyl-CoA thiolases (ACATs) encoded by five genes in Arabidopsis thaliana are targeted either to peroxisomes or cytosol but not to mitochondria. Plant Mol Biol. 2007;63(1):97–108. pmid:17120136
    • 88. Jin H, Song Z, Nikolau BJ. Reverse genetic characterization of two paralogous acetoacetyl CoA thiolase genes in Arabidopsis reveals their importance in plant growth and development. Plant J. 2012;70(6):1015–32. pmid:22332816
    • 89. Vick BA, Zimmerman DC. Biosynthesis of jasmonic acid by several plant species. Plant Physiol. 1984;75(2):458–61. pmid:16663643
    • 90. Pollmann S, Springer A, Rustgi S, Von Wettstein D, Kang C, Reinbothe C, et al. Substrate channeling in oxylipin biosynthesis through a protein complex in the plastid envelope of Arabidopsis thaliana. J Exp Bot. 2019;70(5):1483–95. pmid:30690555
    • 91. Tanaka S, Han X, Kahmann R. Microbial effectors target multiple steps in the salicylic acid production and signaling pathway. Front Plant Sci. 2015;6:349. pmid:26042138
    • 92. Starý T, Satková P, Piterková J, Mieslerová B, Luhová L, Mikulík J, et al. The elicitin β-cryptogein’s activity in tomato is mediated by jasmonic acid and ethylene signalling pathways independently of elicitin–sterol interactions. Planta. 2019;249(3):739–49. pmid:30374914
    • 93. Van den Berg N, Mahomed W, Olivier NA, Swart V, Crampton BG. Transcriptome analysis of an incompatible Persea americana-Phytophthora cinnamomi interaction reveals the involvement of SA-and JA-pathways in a successful defense response. PloS one. 2018;13(10):e0205705. pmid:30332458
    • 94. Caplan JL, Kumar AS, Park E, Padmanabhan MS, Hoban K, Modla S, et al. Chloroplast stromules function during innate immunity. Dev Cell. 2015;34(1):45–57. pmid:26120031
    • 95. Nejat N, Rookes J, Mantri NL, Cahill DM. Plant–pathogen interactions: toward development of next-generation disease-resistant plants. Crit Rev Biotechnol. 2017;37(2):229–37. pmid:26796880
    • 96. Kumar AS, Park E, Nedo A, Alqarni A, Ren L, Hoban K, et al. Stromule extension along microtubules coordinated with actin-mediated anchoring guides perinuclear chloroplast movement during innate immunity. eLife. 2018;7:e23625. pmid:29338837