Comparative Proteomic Analysis Reveals Novel Insights into the Interaction between Rice and Xanthomonas oryzae pv. oryzae

doi:10.21203/rs.3.rs-41837/v2

Download PDF

Research article

Comparative Proteomic Analysis Reveals Novel Insights into the Interaction between Rice and Xanthomonas oryzae pv. oryzae

https://doi.org/10.21203/rs.3.rs-41837/v2

This work is licensed under a CC BY 4.0 License

Journal Publication

published 14 Dec, 2020

Read the published version in BMC Plant Biology →

You are reading this older preprint version

Read the latest preprint version →

Background: Bacterial blight, which is caused by Xanthomonas oryzae pv. oryzae (Xoo), is a devastating rice disease worldwide. Rice introgression line H471, derived from the recurrent parent Huang-Hua-Zhan (HHZ) and the donor parent PSBRC28, exhibits broad-spectrum resistance to Xoo, including to the highly virulent Xoo strain PXO99^A, whereas its parents are susceptible to PXO99^A. To characterize the responses to Xoo, we compared the proteome profiles of the host and pathogen in the incompatible interaction (H471 inoculated with PXO99^A) and the compatible interaction (HHZ inoculated with PXO99^A).

Results: In this study, a total of 374 rice differentially abundant proteins (DAPs) and 117 Xoo DAPs were detected in the comparison between H471 + PXO99^A and HHZ + PXO99^A. Most of the Xoo DAPs related to pathogen virulence, including the outer member proteins, type III secretion system proteins, TonB-dependent receptors, and transcription activator-like effectors, were less abundant in the incompatible interaction than in the compatible interaction. The rice DAPs were mainly involved in secondary metabolic processes, including phenylalanine metabolism and the biosynthesis of flavonoids and phenylpropanoids. Additionally, some DAPs involved in the phenolic phytoalexin and salicylic acid (SA) biosynthetic pathways accumulated much more in H471 than in HHZ after the inoculation with PXO99^A, suggesting that phytoalexin and SA production was induced faster in H471 than in HHZ. Further analysis revealed that the SA content in H471 increased much more rapidly after inoculation than in HHZ, suggesting that SA signaling pathway was activated faster in the incompatible interaction than in compatible interaction.

Conclusions: Overall, our results indicate that during an incompatible interaction between H471 and PXO99^A, rice plants prevent pathogen invasion and also initiate multi-component defense responses that inhibit disease development.

Plant Physiology and Morphology

Plant Molecular Biology and Genetics

Rice

Xanthomonas oryzae pv. oryzae

compatible interaction

incompatible interaction

quantitative proteome

Bacterial blight (BB) caused by the Gram-negative bacterium Xanthomonas oryzae pv. oryzae (Xoo) is one of the most widely distributed and devasting diseases of rice (Oryza sativa L.), leading to substantial yield losses [1]. Developing and deploying broad-spectrum resistant varieties is the most economically viable BB management strategy. Therefore, further characterizing rice-Xoo interactions will provide insights useful for developing and breeding resistant cultivars.

Several types of Xoo virulence factors have been identified, including proteins associated with exopolysaccharide (EPS) production and motility, outer membrane (OM) proteins and TonB-dependent receptors, hypersensitive response and pathogenicity (Hrp) proteins, and type III (T3) effectors [2].

Plants have developed rapid and sophisticated mechanisms that confer resistance to invading pathogens. The cell surface immune receptors (CSIRs) and intracellular immune receptors (IIRs) directly or indirectly recognize apoplastic or cytoplasmic “invasion molecules” to induce weak or strong immune responses [3]. Activated defense responses consist of callose deposition in cell walls, oxidative burst, and accumulation of phytohormones and antimicrobial chemicals (e.g., secondary metabolites and phytoalexins) [4]. In response to pathogen infections, rice produces many phytoalexins, including diterpenoid compounds and the flavonoid sakuranetin [5]. A new class of phenolic phytoalexins with antimicrobial activities against bacterial and fungal pathogens in rice comprises stress-induced phenylamides, which rely on the shikimate pathway for aromatic L-amino acids and the phenylpropanoid pathway for the phenolic acid moieties in phenylamides and sakuranetin [6]. Salicylic acid (SA) is an important signaling hormone in plant defense responses, but it is also a phenolic compound synthesized by plants in response to diverse pathogens [7]. The triggered immunity involves considerable increases in the endogenous levels of SA and its conjugates prior to the induction of pathogenesis-related (PR) proteins and the onset of local and systemic acquired resistance [8]. Additionally, SA biosynthesis in higher plants occurs through the shikimate phenylpropanoid pathway, and may involve two distinct routes, the isochorismate (IC) pathway and the phenylalanine ammonia-lyase (PAL) pathway [9]. Many investigations proved that the bulk of SA is produced from IC via two reactions catalyzed by isochorismate synthase (ICS) and isochorismate pyruvate lyase (IPL) [10-12].

To date, 43 major resistance (R) genes conferring resistance to BB have been identified in rice [13-16]. Although there has been considerable research aimed at characterizing the molecular mechanism mediated by R genes, the interactions between rice and Xoo at the genome level remain unclear. The rapid development of proteomics technology has enabled researchers to comprehensively explore the changes to protein abundances during physiological and biochemical processes. Over the last decade, researchers have applied diverse proteomic experimental techniques to investigate rice resistance to Xoo. For example, Mahmood et al. (2006) used two-dimensional polyacrylamide gel electrophoresis (2-DE) to separate rice cytosolic and membrane proteins extracted from plants inoculated with compatible or incompatible Xoo races [17]. Chen et al. (2007) combined 2-DE and tandem mass spectrometry (MS/MS) analyses to identify plasma membrane proteins involved in the Xa21 signaling pathway [18]. Wang et al. (2013) used 2-DE coupled with matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS) to identify proteins involved in defense responses to Xoo [19].

We previously identified the rice introgression line (IL) H471, which was derived from a hybridization between the susceptible recurrent parent Huang-Hua-Zhan (HHZ) and the donor parent PSBRC28. Line H471 reportedly undergoes a broad-spectrum hypersensitive response (HR) to BB mediated by a new rice resistance gene, Xa39 [16]. In the present study, we conducted a differential proteome analysis of rice lines (H471 and its recurrent parent HHZ) and the strongly virulent Xoo strain PXO99^A using tandem mass tag (TMT) technology coupled with liquid chromatography-quadrupole mass spectrometry (LC-MS/MS) [20]. Our results provide novel information that clarifies the incompatible interactions between Xoo and rice that lead to HR.

Reactions of H471 and HHZ in response to PXO99^A

Rice introgression line (IL) H471 and its recurrent parent HHZ were evaluated following an inoculation with Xoo PXO99^A at the tillering stage. Distinct brown edges along the clipped sites of H471 were visible at 3 dpi, and the brown necrosis at the infection site became more apparent at 5 dpi. These observations indicated the IL exhibited a typical HR to Xoo. In contrast, on the infected HHZ leaves, chlorosis was visible at 3 dpi, and water-soaked lesions rapidly spread along the clipped sites at 5 dpi (Fig. 1a). The lesion lengths of H471 and HHZ at 14 dpi were 0.3 ± 0.2 cm and 16.4 ± 2.2 cm, respectively (Fig. 1b). Moreover, PXO99^A growth in H471 was markedly lower than that in HHZ as follows: more than 10-fold lower at 2 dpi, more than 20-fold lower at 3 dpi, about 190-fold lower at 4 dpi, and more than 1,600-fold lower at 5 dpi (Fig. 1c). We also observed the ultrastructural changes in leaf cells of H471 and HHZ after inoculation with PXO99^A at 3 dpi using transmission electron microscope (TEM). Before inoculation, there were no significant differences in the structure of mesophyll cells between H471 and HHZ. While, at 3 dpi, HHZ had a little damage to the cell membranes, as chloroplast and mitochondria didn’t have a clearly outline. In contrast, H471 displayed clear outline of the whole cell and organelles, with many starch grains (Fig. 1d). These results indicated that H471 is highly resistant to PXO99^A (incompatible interaction between H471 and PXO99^A), whereas HHZ is highly susceptible (compatible interaction between HHZ and PXO99^A).

Protein identification and quantification

The molecular weight range was greater for the identified rice proteins (4,000–469,100 Da) than for the identified Xoo proteins (5,300–195,500 Da) (Additional files 1 and 2). The isoelectric points were 3.5–12.5 and 4.2–12.4 for most rice proteins and Xoo proteins, respectively (Additional files 1 and 2). To improve analytical precision, all quantitative data for peptides identified in multiple fractions were used to quantify proteins. A total of 1,289,153 MS/MS spectra were obtained and were approximately matched to 196,745 known peptide sequences (Additional file 3). Overall, we identified and quantified 46,076 peptides associated with 8,120 different proteins (7,784 rice proteins and 336 Xoo proteins) in the mixed proteomics samples (H471 + PXO99^A and HHZ + PXO99^A) at 2 and 3 dpi (Additional files 1-4). A clustering analysis involving all identified proteins revealed that samples from the same genetic background clustered in the same groups (Additional file 5).

Differentially abundant proteins between incompatible and compatible interactions

Only unique peptides were considered for quantifying proteins, resulting in the identification of 374 rice DAPs and 117 Xoo DAPs. Additionally, 264 rice DAPs and 3 Xoo DAPs accumulated more in H471 + PXO99^A than in HHZ + PXO99^A, with a threshold fold-change > 1.5 and p < 0.05 in t-tests, whereas 171 rice DAPs and 116 Xoo DAPs accumulated less in H471 + PXO99^A than in HHZ + PXO99^A , with a threshold fold-change < 0.67 and p < 0.05 in t-tests (Fig. 2a, 2b, Additional files 6 and 7). Furthermore, 249 rice DAPs and 36 Xoo DAPs were detected between H471 + PXO99^A and HHZ + PXO99^A at 2 dpi, whereas 262 rice DAPs and 116 Xoo DAPs were detected between H471 + PXO99^A and HHZ + PXO99^A at 3 dpi. A total of 137 rice DAPs and 35 Xoo DAPs were common between the two time-points (Fig. 2a and 2b). The hierarchical clustering of all rice DAPs indicated that rice DAPs in HHZ at 2 and 3 dpi were clustered in one subgroup, whereas those in H471 at 2 and 3 dpi were in another subgroup (Fig. 2c). The clustering of all Xoo DAPs suggested that Xoo DAPs in H471 at 2 and 3 dpi were clustered in one subgroup, whereas DAPs in HHZ at 2 dpi were grouped with those in H471 (Fig. 2d, Additional file 6). On the basis of protein abundances, the 374 rice DAPs (Additional file 7) were classified into five groups (G-I to V). The proteins in G-III and G-V were more abundant in H471 than in HHZ, whereas the opposite trend was observed for the DAPs in G-II (Fig. 2c). We also extracted the total proteins of HHZ and H471 at 0, 0.5, 1, 1.5, 2, 2.5, 3.5, 4, and 4.5 dpi, and used specific antibodies against four of the obtained DAPs to validate their abundances at different time-points. The western blot results were highly consistent with the proteomic experiments (Fig. 3). However, the data obtained in the present proteomic study was poorly correlated with the results of our previous transcriptomic analyses, implying there are specific mechanisms for maintaining proper levels of transcripts and proteins [21].

The application of the KOBAS platform to identify enriched GO terms and pathways revealed that Xoo DAPs in H471 were significantly enriched with 14 GO terms, and ribosomal proteins (22.5%) formed the largest set of modulated proteins (Additional file 8). The functional classification of these proteins indicated that most were involved in structural molecule activity and RNA binding. Most of the proteins assigned to the biological process category were involved in translation, protein folding, and generation of precursor metabolites and energy (Additional file 8). The significantly enriched KEGG pathways in Xoo were those related to the ribosome (xop03010, p = 1.4E−14), RNA degradation (xop03018, p = 1.3E−3), methane metabolism (xop00680, p = 9.51E−03), gluconeogenesis (xop00010, p = 1.71E−02), carbon metabolism (xop01200, p = 2.75E−02), and the TCA cycle (xop00020, p = 2.89E−02) (Additional file 9). Regarding the rice DAPs, the three enriched GO terms were secondary metabolic processes, endopeptidase inhibitor activity, and peptidase inhibitor activity, and the five enriched KEGG pathways were biosynthesis of secondary metabolites (osa01110, p = 6.74E−03), phenylalanine metabolism (osa00360, p = 1.23E−02), flavonoid biosynthesis (osa00941, p = 3.52E−02), phenylpropanoid biosynthesis (osa00940, p = 4.31E−02), and vitamin B6 metabolism (osa00750, p = 4.55E−02) (Additional files 8 and 9).

Differentially abundant Xoo proteins between the incompatible and compatible interactions

On the basis of protein functions, the detected Xoo DAPs were grouped into the following categories: type III secretion system (T3SS), transcription activator-like (TAL) effector, nutrient uptake, and uncharacterized protein. In this study, the abundances of six outer membrane proteins (Omps) and nine TonB-dependent receptor-related proteins were significantly different between H471 and HHZ mainly at 3 dpi (12/15) (Fig. 4a). Similarly, one T3SS-related protein (HrpE, PXO_03411), one TAL effector (talC3b, PXO_00505), and a VirK protein (PXO_03361) accumulated much more in HHZ-3d than in HHZ-2d and H471-3d. In this study, the VirK protein (PXO_03361) accumulated less in H471 than in HHZ at 2 and 3 dpi. Additionally, elongation factors (PXO_04524 and PXO_01131) were also less abundant in H471 than in HHZ at 3 dpi (Fig. 4b). Increases in the contents of these Xoo DAPs were more pronounced in HHZ at 3 dpi than at 2 dpi, and the fold-changes of these protein abundances were significantly greater in HHZ-3d vs HHZ-2d than in H471-3d vs H471-2d.

Differentially abundant rice proteins between the incompatible and compatible interactions

In addition to the unknown proteins, the rice DAPs identified in the comparison between H471 and HHZ were grouped into the following five categories: signal transduction, transcription, phytohormone, phytoalexin, and defense response. Hierarchical clustering provided an overview of the differential abundance patterns of the protein types (Fig. 4c-g).

Protein kinases and phosphatases are key co-regulators of protein phosphorylation, and are particularly prominent in signal transduction pathways. In this study, 21 differentially abundant protein kinases and 12 differentially abundant protein phosphatases were identified between the incompatible and compatible interactions (Fig. 4c). Calcium-dependent protein kinases (CDPKs) function as Ca²⁺ sensors and effectors that relay specific Ca²⁺ signatures to downstream components via CDPK-dependent protein phosphorylation. In this study, the accumulation of OsCDPK7/OsCDPK13 (LOC_Os04g49510) was more than 200 times greater in H471 than in HHZ at 2 dpi (Fig. 4c, Additional file 7). Likewise, a western blot assay indicated that OsCDPK13 accumulated more in H471 than in HHZ at 2.5 dpi (Fig. 3). The abundances of two CDPK-related kinases (CRKs), LOC_Os04g25060/LOC_Os04g25650 and OsCRK5 (LOC_Os04g56430), were also significantly different between H471 and HHZ. Additionally, the accumulation of a somatic embryogenesis receptor kinase (SERK), OsSERK2 (LOC_Os04g38480), was approximately 10-fold greater in H471 than in HHZ at 2 dpi (Fig. 4c, Additional file 7). The mitogen-activated protein kinase (MAPK) cascades comprise three protein kinase components, MAPK, MAPK kinase (MAPKK), and MAPKK kinase (MAPKKK). The OsMKK4 (LOC_Os02g54600) abundance was lower in H471 than in HHZ at 2 dpi (Fig. 4c, Additional file 7). As indicated in Fig. 3, OsMKK4 accumulated less in H471 than in HHZ. Moreover, it was maintained at a relatively high abundance in HHZ, whereas in H471, it accumulated at 0.5 dpi, peaked at 2 dpi, and then decreased at 4 dpi. In the proteomic assays, there was no obvious difference in OsMPK6 accumulation between HHZ and H471 at 2 and 3 dpi. However, the western blot assay involving the OsMPK6-specific antibody suggested that the accumulation of OsMPK6 significantly increased after HHZ and H471 were inoculated with PXO99^A, and this protein accumulated much more in H471 than in HHZ from 0.5 to 4.0 dpi (Fig. 3). Additionally, 14 protein phosphatases differentially accumulated in H471 and HHZ, including PP2C30 (LOC_Os03g16170) and two Ser/Thr protein phosphatases (LOC_Os03g13540 and LOC_Os12g44020).

Transcription factors (TF) play key roles in the large-scale transcriptional reprogramming of plants in response to pathogens. A comparison with HHZ revealed 14 DAPs encoding TFs in H471 infected by Xoo (Fig. 4d), including bZIP, MYB, zinc finger, AP2, homeobox, and HSF family members. Seven and five TFs were more and less abundant, respectively, in H471 than in HHZ at 2 or 3 dpi. The MYB family protein LOC_Os01g74020 accumulated significantly less in H471 than in HHZ at 2 dpi, but more in H471 than in HHZ at 3 dpi. A zinc finger family protein (LOC_Os12g18120) accumulated more and less in H471 than in HHZ at 2 and 3 dpi, respectively (Fig. 4d, Additional file 7). Additionally, two bZIP TFs, OsbZIP23 (LOC_Os02g52780) and LOC_Os03g13614, accumulated much more in H471 than in HHZ at 2 and 3 dpi, respectively (Fig. 4d, Additional file 7). The western blot assay result further confirmed that OsbZIP23 was more abundant in H471 than in HHZ at 2, 2.5, and 3 dpi (Fig. 3). Heat stress TFs (Hsfs) also regulate gene expression in response to environmental stress. In this study, OsHsfA2e (LOC_Os03g58160) and OsVOZ2 (LOC_Os05g43950) accumulated less in H471 than in HHZ at 2 dpi (Fig. 4d, Additional file 7).

Phytohormones are signaling molecules that circulate throughout plants and stimulate responses to various environmental stresses. The abundances of two gibberellin receptors (LOC_Os07g06860 and LOC_Os03g57640) differed between H471 and HHZ. Specifically, LOC_Os07g06860 was less abundant in H471 than in HHZ at both 2 and 3 dpi, whereas LOC_Os03g57640 was more abundant in H471 than in HHZ at 3 dpi. Cytokinin-O-glucosyltransferase 3 (LOC_Os10g09990) was less abundant in H471 than in HHZ at both time-points (Fig. 4e). Notably, the following three PALs accumulated much more in H471 than in HHZ: LOC_Os02g41630 and LOC_Os04g43760 at 2 and 3 dpi, and LOC_Os02g41650 at 2 dpi (Fig. 4e).

Phytoalexin accumulation at an infection site is frequently associated with relatively broad-spectrum antimicrobial activity. Additionally, PAL is a key enzyme in the phenylpropanoid pathway contributing to phytoalexin synthesis, and it is also important for SA synthesis. We observed that several proteins involved in the phenylpropanoid pathway, including Os4CL3 (LOC_Os02g08100), two O-methyltransferases (OMTs) (LOC_Os10g02880 and LOC_Os08g06100), and two caffeoyl-CoA O-methyltransferases (COMTs) (LOC_Os06g06980 and LOC_Os08g38900), were considerably more abundant in H471 than in HHZ at one or both time-points after inoculation (Fig. 4f). Moreover, two DAPs associated with the shikimate pathway were identified. The accumulation of chorismate mutase (LOC_Os01g55870) was lower in H471 than in HHZ at 2 dpi, but higher at 3 dpi. In contrast, shikimate kinase (LOC_Os02g51410) accumulated much more in H471 than in HHZ at 2 dpi, but less at 3 dpi. Chalcone-flavonone isomerase (LOC_Os11g02440), which is related to sakuranetin biosynthesis, also accumulated more in H471 than in HHZ at 3 dpi (Fig. 4f). Furthermore, three proteins (LOC_Os03g15050, LOC_Os01g07730, and LOC_Os03g27230) involved in phosphoenol pyruvate biosynthesis accumulated differently in H471 and HHZ. The abundance of LOC_Os03g15050 was lower in H471 than in HHZ at 3 dpi, whereas the abundances of LOC_Os01g07730 and LOC_Os03g27230 were higher in H471 than in HHZ at 2 and 3 dpi. (Fig. 4f).

Three and two defense-related proteins were identified as less and more abundant, respectively, in H471 than in HHZ plants inoculated with Xoo. Two NB-ARC domain-containing protein family members (LOC_Os11g44990 and LOC_Os11g45090) and an NBS-LRR disease resistance protein (LOC_Os11g38580) accumulated less in H471 than in HHZ. Both OsPR1b (LOC_Os01g28450) and a harpin-induced protein (LOC_Os12g06220) accumulated much more in H471 than in HHZ. Notably, the accumulation of LOC_Os12g06220 was more than 160-times higher in H471 than in HHZ at 2 dpi (Fig. 4g, Additional file 7). Additionally, a harpin-induced protein (LOC_Os04g33990) accumulated more and less at 2 and 3 dpi, respectively, in H471 than in HHZ (Fig. 4g).

Abundances of proteins associated with SA signaling during interactions between rice and Xoo

To further investigate whether SA contributes to the disease resistance of H471, we analyzed the abundances of the SA-inducible OsPR1a and OsPR1b. Western blot assays indicated that in plants inoculated with PXO99^A, OsPR1a levels increased at 1 and 0.5 dpi in HHZ and H471, respectively, and were much higher in H471 than in HHZ from 2.5 to 4 dpi (Additional file 10a). Because the accumulation of OsPR1b was too low to detect, we were unable to quantify its abundance in HHZ, and only unclear bands were detected at 2.5 and 4.5 dpi for H471 (Additional file 10a). After inoculation with PXO99^A, the SA content in H471 increased rapidly and by twice as much as in HHZ at 1.5 dpi. The SA content was much higher in H471 than in HHZ at 2.5 dpi and 3.5 dpi, and it reduced in HHZ at 4.5 dpi (Fig. 5). These results suggested that the SA signaling pathway was activated in response to the Xoo infection in H471.

Proteomic strategies have recently been successfully applied to analyze plant disease resistance. A number of proteins associated with BB resistance responses in rice have been identified. However, little information is available on the global changes in pathogen proteins during plant-pathogen interactions. In this study, the DAPs of the plant host and pathogen were examined simultaneously. Specifically, quantitative proteomic technology was used to compare incompatible (H471 and PXO99^A) and compatible (HHZ and PXO99^A) interactions. The differences in the proteome profiles of rice and Xoo between the incompatible and compatible interactions revealed interesting molecular mechanisms underlying the broad-spectrum resistance of rice to BB.

The abundances of proteins related to Xoo invasion and growth in rice were lower during the incompatible interaction than during the compatible interaction

Entry into host cells is a crucial step in bacterial infections. Bacterial pathogens use diverse virulence-related factors, including the T3SS, for host cell binding, the triggering of membrane ruffles, and cooperative invasions. The Omps of Gram-negative bacteria are essential for bacterial survival during eukaryotic cell invasions [22, 23]. The OmpA protein is required for bacterial infection/colonization and cellular adhesion/development in host cells [24]. Bacterial TonB-dependent receptors, which are involved in the uptake of nutrients from the surrounding environment, and extracellular enzymes are also important for the interactions between Xoo and rice [19]. An elongation factor is essential for viability and is required for protein synthesis. In this study, DAPs related to six Omps, nine TonB-dependent receptors, and two elongation factors were detected, and all of them accumulated less in H471 than in HHZ. These results suggested that PXO99^A invasion and growth were significantly more inhibited in the incompatible interaction than in the compatible interaction. The dramatically decreased Xoo population in H471 was consistent with the differential protein abundances.

The T3SS, encoded by hrp genes, plays an important role in the interaction between Xoo and rice because it injects T3 effectors into plant cells [25]. The Xoo T3 effectors include two collections of proteins, namely TAL effectors and Xanthomonas outer proteins (Xop). The TAL effectors can recognize and bind to specific DNA sequences within the promoters of corresponding host genes, resulting in the transcriptional activation of the genes mediating disease susceptibility or resistance [26]. In the current study, HrpE (PXO_03411) accumulated less in H471 than in HHZ. Notably, the T3 effectors talC3b (PXO_00505) and OsVOZ2 (LOC_Os05g43950) also accumulated less in H471 than in HHZ. A previous study proved that interactions between XopN_KXO85 and OsVOZ2 increase rice susceptibility to Xoo [27]. However, in the present study, we were unable to identify the target gene of talC3b based on the list of rice DAPs. Additionally, the PXO99^A XopN protein abundance was similar in H471 and HHZ at the two examined time-points. This might have been due to a shift in the timing of gene expression as the host-pathogen interaction developed. We speculated that the effectors talC3b and XopN may serve as signaling molecules that modulate the host response to infection when they are released in plant tissues, wherein they activate their target genes whose protein products contribute to the susceptibility of the host plants in a compatible interaction. In contrast, these proteins do not mediate a susceptible response in incompatible interactions because of their lower abundances. Clarifying the functions of the two effectors and their host targets will be a key step for further characterizing the interaction between H471 and Xoo.

Diverse rice proteins related to signal transduction and transcriptional regulation are involved in the plant-pathogen interaction

During plant-pathogen interactions, plants not only perceive and prevent the pathogen from invading, they also initiate defense responses by activating multi-component systems via modulated expression or abundance and/or post-translational modifications of the associated proteins. After the inoculation with Xoo, the accumulation of diverse categories of rice protein kinases and TFs differed significantly between the incompatible and compatible interactions.

Receptor-like kinases (RLKs) reportedly regulate defense processes [28]. On the basis of amino acid sequence and structural differences, RLKs have been categorized into several subfamilies, including leucine-rich repeat RLKs (LRR-RLKs), cysteine-rich repeat RLKs (CRKs), domain of unknown function 26 RLKs, and S-domain RLKs, among others [29]. As common key signaling molecules, CRKs influence diverse stress responses. For example, OsCRK5 expression is up-regulated in rice in response to the blast fungus or brown planthopper [30, 31]. In the current study, OsCRK5 was more abundant in H471 than in HHZ, in contrast to another CRK, LOC_Os04g25060/LOC_Os04g25650, which was less abundant in H471 than in HHZ. Genes encoding CDPKs are critical for abiotic stress responses [32, 33], but they are also expressed in response to biotic stresses in plants [34]. A previous study indicated that OsCDPK7 (OsCDPK13) expression can be induced by an exposure to chilling stress and high salt concentrations [35]. In the present study, OsCDPK7/OsCDPK13 (LOC_Os04g49510) was significantly more abundant in H471 than in HHZ. The MAPK cascade has been implicated in signaling during plant defense responses to pathogens and various environmental stimuli [36]. In rice, chitin elicitors may activate OsMPK1, OsMPK5, OsMKK4, and the OsMKK4-OsMPK6 cascade to regulate defense activities, including antimicrobial compound biosynthesis, and induce plant cell death [37, 38]. In the current study, OsMKK4 and OsMPK6 accumulated much more in H471 than in HHZ at 1.5–3 dpi. Moreover, 14 phosphatase DAPs, including PP2C30 and two Ser/Thr protein phosphatases, accumulated differently in H471 and HHZ. A previous investigation revealed that PP2C30 interacts with the ABA receptor PYL/RCAR5 and SAPK2 to regulate ABA-dependent gene expression [39]. It also functions as an upstream regulator of HOX12 expression, acting directly through EUI1 to control panicle exsertion in rice [40].

The bZIP TF family participates in plant responses to abiotic stress via the ABA signaling pathway. For example, OsbZIP71 RNAi knockdown transgenic lines are highly sensitive to salt, PEG-induced osmotic stress, and ABA [41]. Additionally, OsbZIP23 is a central regulator of ABA signaling and biosynthesis as well as drought tolerance in rice [42]. Moreover, OsbZIP46 positively regulates ABA signaling and drought tolerance in rice, depending on whether it is activated [43]. Our data indicated that two bZIP TFs, namely OsbZIP23 (LOC_Os02g52780) and LOC_Os03g13614, accumulated much more in H471 than in HHZ at 2 and 3 dpi. Furthermore, western blot assays confirmed that OsbZIP23 accumulated more in H471 than in HHZ. These results imply that bZIP TFs are important for the incompatible interaction between rice and Xoo.

Phytoalexin and SA biosynthetic pathways were more enhanced in rice during the incompatible interaction than during the compatible interaction following the Xoo infection

In response to pathogens, plants may produce many phytoalexins, including the flavonoid sakuranetin [5]. Sakuranetin is the only known phenolic phytoalexin in rice, and phenylamides are involved in its biosynthesis. The accumulating phenylamides or sakuranetin in plants help reinforce the cell wall and are toxic to pathogens at the infection site [6, 44-46]. Additionally, SA is a small phenolic compound that has an important regulatory role in multiple physiological processes. There are two SA biosynthetic pathways in plants: one involves cinnamate and includes a reaction catalyzed by PAL, whereas the other requires chorismate and includes a reaction catalyzed by ICS [47]. The accumulation of SA always coincides with the up-regulated expression of antimicrobial PR genes, which enhances disease resistance responses [7]. Methylated derivatives of SA spread from the infected tissue to distal tissues, thereby inducing systemic acquired resistance [48].

Phenylalanine ammonia-lyases are key enzymes in the phenylpropanoid pathway, which mediates the biosynthesis of the flavonoid-type phytoalexin, sakuranetin, and SA [48]. In the current study, 14 DAPs associated with phenolic phytoalexin biosynthesis were detected. Three PALs (OsPALs) accumulated much more in H471 than in HHZ, whereas no proteins contributing to the chorismite-based SA biosynthetic pathway were identified as a DAP. Furthermore, SA content in H471 increased rapidly in H471 after inoculation with PXO99^A, and the accumulation level of OsPR1b was about two-fold higher in H471 than in HHZ. These results suggest that the phytoalexin and SA biosynthetic pathways were activated faster during the incompatible interaction than during the compatible interaction, implying that phytoalexins and SA are involved in the HR of H471 (Fig. 6).

The results of this study confirm that analyzing proteome-level differences is useful for characterizing plant host-pathogen interactions. The proteomic data presented herein provide a global overview of the protein abundance changes in rice and Xoo during their interactions. Further clarifying the functions of DAPs related to Xoo virulence and host responses will lead to a more comprehensive understanding of the mechanisms underlying plant host-pathogen interactions.

Plant materials and artificial inoculation

A rice introgression line (IL) (H471) with broad-spectrum BB resistance was selected from a BC₁F₆ population derived from a cross between the recurrent parent HHZ and the donor parent PSBRC28. Lines H471 and HHZ were used in this study [16]. Seeds were sown and the germinated seedlings were grown in a seedling nursery with a 13-h light (30°C): 11-h dark (26°C) photoperiod. The 30-day-old seedlings were transplanted to a screenhouse at the Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, China. Lines H471 and HHZ were artificially inoculated at the early tillering stage with a bacterial solution comprising 10⁵ cells L⁻¹ Xoo strain PXO99^A (Philippines race 6) according to the scissors-dip method [16]. Lesion lengths and the growth curves of PXO99^A in H471 and HHZ were evaluated as previously described [16].

Transmission electron microscopy (TEM)

Leaves of H471 and HHZ inoculated by Xoo strain PXO99^A at 0 dpi and 3 dpi were collected for TEM analysis. A 2 mm leaf fragment under 1 mm of the lesions was fixed in 2% glutaraldehyde in 0.1 mol L⁻¹phosphate buffer (pH 7.2) for 4 h at 4°C and then washed 3 times in the same buffer, which was followed by post-fixation in 1% osmium tetroxide (OsO4) for 2 h at room temperature and 3 rinses using the same buffer. After staining with 2.5% uranyl acetate, specimens were further dehydrated in a graded ethanol series and embedded in Spurr’s medium. Polymerization took place for 24 h at 40°C, followed by 24 h at 60°C. Specimens were cut to a thickness of 40-60 nm on RMC PowerTome-XL ultramicrotome (RMC, Arizona, USA), ultrathin sections were double-stained with 2% (w/v) uranylacetate and 2.6% (w/v) lead citrate aqueous solution and examined with a Hitachi H-7500 transmission electron microscope (Hitachi Company, Japan).

Preparation of plant and pathogen samples for the proteomic analysis

In this study, 10 plants per rice line were included in the proteomic analysis. Specifically, 5-10 of the uppermost leaves of each plant were inoculated with PXO99^A. At 2- and 3-days post-inoculation (dpi), leaf tips (3 cm long) were collected from three plants (more than three tips per plant) and combined to form one biological replicate, with three independent biological replicates prepared for each genotype at each time-point. The H471 and HHZ leaf samples collected at 2 dpi were designated as H471-2d and HHZ-2d, respectively, whereas the leaf samples collected at 3 dpi were named H471-3d and HHZ-3d, respectively. The reference sample consisted of a mixture of an equal amount of the 12 samples (HHZ-2d-1/2/3, HHZ-3d-1/2/3, H471-2d-1/2/3, and H471-3d-1/2/3).

Protein extraction and digestion

Tissues were ground in liquid nitrogen and suspended in pre-cooled acetone (–20°C) containing 10% (w/v) trichloroacetic acid. After vigorously mixing the solutions, the proteins were precipitated for 16 h at –20°C and then collected by centrifuging at 16,000 × g (Beckman Allegra 22R) for 30 min at 4°C. The pelleted proteins were washed three times with pre-cooled acetone. After a lyophilization, the pellets were dissolved in extraction buffer [8 M urea, 50 mM triethylammonium bicarbonate (TEAB), pH 8.5, 1% SDS, and protease inhibitor cocktail (Roche, Indianapolis, USA)]. The soluble protein extracts were collected by centrifuging at 16,000 × g for 20 min at 4°C, after which the protein concentrations were determined with the BCA assay kit (Pierce, Rockford, IL, USA).

Sample processing for the LC-MS analysis

Tandem mass tags with different reporter ions (126–131 Da) were conjugated to proteins as isobaric tags with the TMT10plex kit (Pierce) for a subsequent relative quantification. Briefly, 100 mM TEAB was added to the protein solution (100 μg total protein) for a final volume of 100 μL. After adding 5 μL 200 mM Tris(2-carboxyethyl)phosphine (TCEP), the mixture was incubated at 55°C for 1 h. Immediately before use, one tube of iodoacetamide (9 mg) was dissolved in 132 μL 100 mM TEAB to prepare the 375 mM iodoacetamide solution. A 5-μL aliquot of 375 mM iodoacetamide was added to the protein sample and the mixture was incubated for 30 min at room temperature in darkness. Proteins were precipitated with six volumes of pre-chilled (–20°C) acetone. The solution was then centrifuged at 8,000 × g for 10 min at 4°C and the resulting pellets were washed twice, once with pre-chilled (–20°C) 90% (v/v) methanol and then with pre-chilled (–20°C) 100% methanol. Protein pellets were briefly dried and resuspended in an adequate amount of 100 mM TEAB (pH 8.5) containing 5% acetonitrile for a concentration of 1 μg/μL. Proteins were digested overnight at 37°C with 2.5 μg trypsin (Sigma, USA) and then stored at –80°C until analyzed. One tube of TMT Label Reagent (Pierce) was added to each 100 μg sample and the labeling reaction was completed at room temperature for 1 h. The reaction was quenched by adding 8 μL 5% hydroxylamine. The samples were labeled with TMT10plex tags. The reference sample was labeled with TMT¹⁰-131 (Additional file 11). The 12 labeled samples (HHZ-2d-1/2/3, HHZ-3d-1/2/3, H471-2d-1/2/3, and H471-3d-1/2/3) were divided into pools A2 (HHZ-3d-1/2/3 and H471-2d-1/2/3) and B2 (HHZ-2d-1/2/3 and H471-3d-1/2/3). Pools A2 and B2 (30 μg protein per sample) were separately combined with 30 μg labeled reference sample protein in new microcentrifuge tubes.

The two pooled protein samples were lyophilized, reconstituted in solvent A (2% acetonitrile, pH 10), loaded onto an XBridge PST C18 column (130 Å, 5 μm, 250 × 4.6 mm) (Waters, USA), resolved by basic reversed-phase liquid chromatography (off-line) with a gradient of 5%–95% solvent B (90% acetonitrile, pH 10) over 40 min, and divided into 20 fractions per pool. The 40 fractions of both pools were collected according to the chromatogram map, concatenated to 20 fractions, centrifuged under a vacuum, and stored at −80°C until analyzed by LC-MS/MS.

LC-MS/MS analysis

The LC-MS/MS analysis involving a Q Exactive hybrid quadrupole orbitrap mass spectrometer (Thermo Scientific, CA, USA) was performed according to the manufacturer-recommended protocols by CapitalBio Technology (China). Briefly, the peptide mixture was separated by reversed-phase chromatography with a DIONEX nano-UPLC system comprising an Acclaim C18 PepMap100 nano-Trap column (75 μm × 2 cm, 2 μm particle size) (Thermo Scientific) connected to an Acclaim PepMap RSLC C18 analytical column (75 μm × 25 cm, 2 μm particle size) (Thermo Scientific). Before loading, samples were dissolved in sample buffer containing 4% acetonitrile and 0.1% formic acid. A linear gradient of mobile phase B (0.1% formic acid in 99.9% acetonitrile) from 3% to 30% in 43 min followed by an increase to 80% mobile phase B in 1 min was applied at a flow rate of 300 nL min⁻¹. The nano-LC was coupled online with the Q Exactive mass spectrometer using a stainless-steel Emitter linked to a nanospray ion source.

The MS survey scan was performed with full scans [350–1,600 mass-to-charge ratio (m/z)] acquired with an Orbitrap mass analyzer. The mass resolution was 70,000 at 400 m/z. The 20 most intense precursor ions from a survey scan were selected for MS/MS for each duty cycle and detected at a mass resolution of 35,000 at 400 m/z. All tandem mass spectra were produced by higher-energy collision dissociation. The dynamic exclusion was set for 20 s.

Protein identification, quantification, and analysis

Proteome Discoverer (version 1.4) software (Thermo Scientific) was used to search the Oryza sativa Nipponbare reference genome RGAP 7.0 database [49] and the Xanthomonas oryzae pv. oryzae PXO99^A RefSeq database [50] (32,431 total proteins) with the Sequest algorithms. The precursor mass tolerance and fragment mass tolerance were set to 20 ppm and 0.02 Da, respectively. With the exception of TMT modifications (N-terminus and lysine residues), carbamidomethylation of cysteine was set as a fixed modification and oxidation at methionine was set as a dynamic modification. The digesting enzyme was trypsin, with two missed cleavages allowed. Proteins were quantified with the TMT10plex Isobaric Label Reagent Set (Pierce). Specifically, the unique peptides were used to quantify proteins, with 0 or missing values replaced with the minimum intensity. Additionally, values were rejected if any quantification channels were missing. The false discovery rates (FDRs) of the identified proteins were estimated with the target-decoy approach of Proteome Discoverer. A cutoff value of 1% was used to identify proteins. The protein normalization was conducted as previously described [51]. All protein ratios were normalized against the median protein ratio, which was 1 after the data normalization. The differentially abundant proteins (DAPs) between two samples were identified following a comparison of the mean abundances of three replicates with Student’s t-test. The DAPs between H471 and HHZ at each time-point after the PXO99^A inoculation (i.e., H471-2d vs HHZ-2d and H471-3d vs HHZ-3d) were screened based on a 1.5-fold change in mean abundance and p < 0.05 [52, 53]. All raw mass spectrometry data and search results have been deposited in the iProX database (https://www.iprox.org/), with the data set identifier IPX0001741000 [54].

Functional analyses based on gene ontology and pathway enrichment

Proteins were functionally analyzed with Gene Ontology (GO) annotations with KOBAS 2.0 [55]. Specifically, the proteins were classified according to the molecular function, biological process, and cellular component categories. The Kyoto Encyclopedia of Genes and Genomes (KEGG) database was screened to identify the enriched pathways among the proteins [56].

For the hierarchical clustering based on different protein functional classifications (GO terms and KEGG pathways), the cluster membership was visualized as a heat map, which was prepared with the heatmap.2 function of the gplots R package.

Protein immunoblot analysis

Proteins were extracted from HHZ and H471 seedlings with the Plant Protein Extraction Kit (CWBIO, Beijing, China) and quantified with the BSA Protein Assay Kit (CWBIO). Proteins were separated by 10% polyacrylamide gel electrophoresis with protein markers. Polyclonal anti-OsbZIP23, anti-OsCDPK13, anti-OsMKK4, anti-OsMPK6, anti-OsPR1a, and anti-OsPR1b (Beijing Protein Innovation, Beijing, China) were used as primary antibodies (1:2,000 dilution), whereas horseradish peroxidase-conjugated goat anti-rabbit (Abmart, Shanghai, China) served as the secondary antibody (1:10,000 dilution). Coomassie Brilliant Blue staining of the large subunit of Rubisco was used as a loading control. The western blot experiment was independently repeated at least three times, with essentially the same results. Representative results are presented.

Quantification of endogenous salicylic acid

Leaf tips (3 cm long) of HHZ and H471 were harvested at 0 dpi, 1.5 dpi, 2.5 dpi, 3.5 dpi, and 4.5 dpi after inoculation with PXO99^A and stored at −80°C. Quantification of endogenous SA was conducted as previously described [57]. Each time-point was performed with two biological replications.

Xoo: Xanthomonas oryzae pv. oryzae; HR: hypersensitive response; DAPs: Differentially abundant proteins; TMT: Tandem mass tag; LC-MS: Liquid chromatography-quadrupole mass spectrometry; TAL: Transcription activator-like; Omps: Outer membrane proteins; SA: Salicylic acid; NBS-LRR: Nucleotide-binding domain leucine rich repeat; NB-ARC: Nucleotide-binding adaptor shared by APAF-1, R proteins and CED-4

Ethics approval and consent to participate

No applicable.

Consent for publication

No applicable.

Availability of data and materials

All raw mass spectrometry data and search results have been deposited at the iProX (https://www.iprox.org/) with the data set identifiers IPX0001741000.

Competing interests

The authors declare that they have no competing interests.

Funding

This research was supported by grants from the National Natural Science Foundation of China (31661143009, 31571632), the National High-tech Program of China (2014AA10A603), the Bill & Melinda Gates Foundation (OPP1130530), and the CAAS Innovative Team Award. These funding bodies have no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author contributions

Y.Z., Z.L., and C.V. designed the research; F.Z.^1,4, L.H., and D.Z. performed the research; F.Z.¹ and F.Z.^1,4 analyzed the data; F.Z.¹, F.Z.^1,4, and Y.Z. wrote the manuscript. All authors have read and approved the manuscript.

Acknowledgments

We would like to thank Huw Tyson, PhD, from Liwen Bianji, Edanz Editing China (www.liwenbianji.cn/ac), for editing the English text of a draft of this manuscript.

Nino-Liu DO, Ronald PC, Bogdanove AJ: Xanthomonas oryzae pathovars: model pathogens of a model crop. Molecular plant pathology 2006, 7(5):303-324.
Zhang F, Du Z, Huang L, Vera Cruz C, Zhou Y, Li Z: Comparative transcriptome profiling reveals different expression patterns in Xanthomonas oryzae pv. oryzae strains with putative virulence-relevant genes. PloS one 2013, 8(5):e64267.
Kanyuka K, Rudd JJ: Cell surface immune receptors: the guardians of the plant's extracellular spaces. Curr Opin Plant Biol 2019, 50:1-8.
Liu WD, Liu JL, Triplett L, Leach JE, Wang GL: Novel Insights into Rice Innate Immunity Against Bacterial and Fungal Pathogens. Annual Review Of Phytopathology 2014, 52:213-241.
Ahuja I, Kissen R, Bones AM: Phytoalexins in defense against pathogens. Trends In Plant Science 2012, 17(2):73-90.
Park HL, Yoo Y, Hahn TR, Bhoo SH, Lee SW, Cho MH: Antimicrobial activity of UV-induced phenylamides from rice leaves. Molecules 2014, 19(11):18139-18151.
Vlot AC, Dempsey DA, Klessig DF: Salicylic Acid, a multifaceted hormone to combat disease. Annu Rev Phytopathol 2009, 47:177-206.
De Vleesschauwer D, Xu J, Hofte M: Making sense of hormone-mediated defense networking: from rice to Arabidopsis. Front Plant Sci 2014, 5:611.
Dempsey DA, Vlot AC, Wildermuth MC, Klessig DF: Salicylic Acid biosynthesis and metabolism. Arabidopsis Book 2011, 9:e0156.
Wildermuth MC, Dewdney J, Wu G, Ausubel FM: Isochorismate synthase is required to synthesize salicylic acid for plant defence. Nature 2001, 414(6863):562-565.
Dempsey DA, Klessig DF: How does the multifaceted plant hormone salicylic acid combat disease in plants and are similar mechanisms utilized in humans? BMC Biol 2017, 15(1):23.
Klessig DF, Choi HW, Dempsey DA: Systemic Acquired Resistance and Salicylic Acid: Past, Present, and Future. Mol Plant Microbe Interact 2018, 31(9):871-888.
Kim SM, Suh JP, Qin Y, Noh TH, Reinke RF, Jena KK: Identification and fine-mapping of a new resistance gene, Xa40, conferring resistance to bacterial blight races in rice (Oryza sativa L.). Theor Appl Genet 2015, 128(10):1933-1943.
Hutin M, Sabot F, Ghesquiere A, Koebnik R, Szurek B: A knowledge-based molecular screen uncovers a broad-spectrum OsSWEET14 resistance allele to bacterial blight from wild rice. Plant J 2015, 84(4):694-703.
Busungu C, Taura S, Sakagami JI, Ichitani K: Identification and linkage analysis of a new rice bacterial blight resistance gene from XM14, a mutant line from IR24. Breed Sci 2016, 66(4):636-645.
Zhang F, Zhuo DL, Huang LY, Wang WS, Xu JL, Cruz CV, Li ZK, Zhou YL: Xa39, a novel dominant gene conferring broad-spectrum resistance to Xanthomonas oryzae pv. oryzae in rice. Plant Pathology 2015, 64(3):568-575.
Mahmood T, Jan A, Kakishima M, Komatsu S: Proteomic analysis of bacterial-blight defense-responsive proteins in rice leaf blades. Proteomics 2006, 6(22):6053-6065.
Chen F, Yuan YX, Li Q, He ZH: Proteomic analysis of rice plasma membrane reveals proteins involved in early defense response to bacterial blight. Proteomics 2007, 7(9):1529-1539.
Wang YM, Kim SG, Wu JN, Huh HH, Lee SJ, Rakwal R, Agrawal GK, Park ZY, Young Kang K, Kim ST: Secretome analysis of the rice bacterium Xanthomonas oryzae (Xoo) using in vitro and in planta systems. Proteomics 2013, 13(12-13):1901-1912.
Rauniyar N, Yates JR, 3rd: Isobaric labeling-based relative quantification in shotgun proteomics. J Proteome Res 2014, 13(12):5293-5309.
Sharabiani MTA, Siermala M, Lehtinen TO, Vihinen M: Dynamic covariation between gene expression and proteome characteristics. Bmc Bioinformatics 2005, 6.
Miller VL, Bliska JB, Falkow S: Nucleotide sequence of the Yersinia enterocolitica ail gene and characterization of the Ail protein product. J Bacteriol 1990, 172(2):1062-1069.
Tommassen J, Stoorvogel J, van Bussel MJ, van de Klundert JA: Molecular characterization of an Enterobacter cloacae outer membrane protein (OmpX). J Bacteriol 1991, 173(1):156-160.
Rott PC, Fleites L, Marlow G, Royer M, Gabriel DW: An OmpA family outer membrane protein is required for both disease symptom development and sugarcane stalk colonization by Xanthomonas albilineans. Phytopathology 2009, 99(6):S110-S111.
Tang JL, Feng JX, Li QQ, Wen HX, Zhou DL, Wilson TJ, Dow JM, Ma QS, Daniels MJ: Cloning and characterization of the rpfC gene of Xanthomonas oryzae pv. oryzae: involvement in exopolysaccharide production and virulence to rice. Mol Plant Microbe Interact 1996, 9(7):664-666.
Romer P, Recht S, Strauss T, Elsaesser J, Schornack S, Boch J, Wang S, Lahaye T: Promoter elements of rice susceptibility genes are bound and activated by specific TAL effectors from the bacterial blight pathogen, Xanthomonas oryzae pv. oryzae. New Phytol 2010, 187(4):1048-1057.
Cheong H, Kim CY, Jeon JS, Lee BM, Moon JS, Hwang I: Xanthomonas oryzae pv. oryzae Type III Effector XopN Targets OsVOZ2 and a Putative Thiamine Synthase as a Virulence Factor in Rice. Plos One 2013, 8(9).
Osakabe Y, Maruyama K, Seki M, Satou M, Shinozaki K, Yamaguchi-Shinozaki K: Leucine-rich repeat receptor-like kinase1 is a key membrane-bound regulator of abscisic acid early signaling in Arabidopsis. Plant Cell 2005, 17(4):1105-1119.
Shiu SH, Bleecker AB: Receptor-like kinases from Arabidopsis form a monophyletic gene family related to animal receptor kinases. Proc Natl Acad Sci U S A 2001, 98(19):10763-10768.
Kim ST, Kim SG, Hwang DH, Kang SY, Kim HJ, Lee BH, Lee JJ, Kang KY: Proteomic analysis of pathogen-responsive proteins from rice leaves induced by rice blast fungus, Magnaporthe grisea. Proteomics 2004, 4(11):3569-3578.
Wei Z, Hu W, Lin Q, Cheng X, Tong M, Zhu L, Chen R, He G: Understanding rice plant resistance to the Brown Planthopper (Nilaparvata lugens): a proteomic approach. Proteomics 2009, 9(10):2798-2808.
Ludwig AA, Romeis T, Jones JD: CDPK-mediated signalling pathways: specificity and cross-talk. J Exp Bot 2004, 55(395):181-188.
Saijo Y, Hata S, Kyozuka J, Shimamoto K, Izui K: Over-expression of a single Ca2+-dependent protein kinase confers both cold and salt/drought tolerance on rice plants. Plant J 2000, 23(3):319-327.
Romeis T, Herde M: From local to global: CDPKs in systemic defense signaling upon microbial and herbivore attack. Curr Opin Plant Biol 2014, 20:1-10.
Wan B, Lin Y, Mou T: Expression of rice Ca(2+)-dependent protein kinases (CDPKs) genes under different environmental stresses. FEBS Lett 2007, 581(6):1179-1189.
Meng X, Zhang S: MAPK cascades in plant disease resistance signaling. Annu Rev Phytopathol 2013, 51:245-266.
Kishi-Kaboshi M, Okada K, Kurimoto L, Murakami S, Umezawa T, Shibuya N, Yamane H, Miyao A, Takatsuji H, Takahashi A et al: A rice fungal MAMP-responsive MAPK cascade regulates metabolic flow to antimicrobial metabolite synthesis. Plant Journal 2010, 63(4):599-612.
Yoo SJ, Kim SH, Kim MJ, Ryu CM, Kim YC, Cho BH, Yang KY: Involvement of the OsMKK4-OsMPK1 Cascade and its Downstream Transcription Factor OsWRKY53 in the Wounding Response in Rice. Plant Pathology J 2014, 30(2):168-177.
Kim H, Hwang H, Hong JW, Lee YN, Ahn IP, Yoon IS, Yoo SD, Lee S, Lee SC, Kim BG: A rice orthologue of the ABA receptor, OsPYL/RCAR5, is a positive regulator of the ABA signal transduction pathway in seed germination and early seedling growth. J Exp Bot 2012, 63(2):1013-1024.
Gao SP, Fang J, Xu F, Wang W, Chu CC: Rice HOX12 Regulates Panicle Exsertion by Directly Modulating the Expression of ELONGATED UPPERMOST INTERNODE1. Plant Cell 2016, 28(3):680-695.
Zhang CY, Liu J, Zhao T, Gomez A, Li C, Yu CS, Li HY, Lin JZ, Yang YZ, Liu B et al: A Drought-Inducible Transcription Factor Delays Reproductive Timing in Rice. Plant Physiology 2016, 171(1):334-343.
Xiang Y, Tang N, Du H, Ye HY, Xiong LZ: Characterization of OsbZIP23 as a Key Player of the Basic Leucine Zipper Transcription Factor Family for Conferring Abscisic Acid Sensitivity and Salinity and Drought Tolerance in Rice. Plant Physiology 2008, 148(4):1938-1952.
Tang N, Zhang H, Li XH, Xiao JH, Xiong LZ: Constitutive Activation of Transcription Factor OsbZIP46 Improves Drought Tolerance in Rice. Plant Physiology 2012, 158(4):1755-1768.
Edreva AM, Velikova VB, Tsonev TD: Phenylamides in plants. Russ J Plant Physl+ 2007, 54(3):287-301.
von Roepenack-Lahaye E, Newman MA, Schornack S, Hammond-Kosack KE, Lahaye T, Jones JDG, Daniels MJ, Dow JM: p-coumaroylnoradrenaline, a novel plant metabolite implicated in tomato defense against pathogens. Journal Of Biological Chemistry 2003, 278(44):43373-43383.
Yang Q, He Y, Kabahuma M, Chaya T, Kelly A, Borrego E, Bian Y, El Kasmi F, Yang L, Teixeira P et al: A gene encoding maize caffeoyl-CoA O-methyltransferase confers quantitative resistance to multiple pathogens. Nat Genet 2017.
Yang L, Li BS, Zheng XY, Li JG, Yang M, Dong XN, He GM, An CC, Deng XW: Salicylic acid biosynthesis is enhanced and contributes to increased biotrophic pathogen resistance in Arabidopsis hybrids. Nature Communications 2015, 6.
Duan L, Liu HB, Li XH, Xiao JH, Wang SP: Multiple phytohormones and phytoalexins are involved in disease resistance to Magnaporthe oryzae invaded from roots in rice. Physiol Plantarum 2014, 152(3):486-500.
Kawahara Y, de la Bastide M, Hamilton JP, Kanamori H, McCombie WR, Ouyang S, Schwartz DC, Tanaka T, Wu J, Zhou S et al: Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data. Rice (N Y) 2013, 6(1):4.
Salzberg SL, Sommer DD, Schatz MC, Phillippy AM, Rabinowicz PD, Tsuge S, Furutani A, Ochiai H, Delcher AL, Kelley D et al: Genome sequence and rapid evolution of the rice pathogen Xanthomonas oryzae pv. oryzae PXO99A. BMC Genomics 2008, 9:204.
Ross PL, Huang YN, Marchese JN, Williamson B, Parker K, Hattan S, Khainovski N, Pillai S, Dey S, Daniels S et al: Multiplexed protein quantitation in Saccharomyces cerevisiae using amine-reactive isobaric tagging reagents. Mol Cell Proteomics 2004, 3(12):1154-1169.
Hu J, Ren B, Dong S, Liu P, Zhao B, Zhang J: Comparative proteomic analysis reveals that exogenous 6-benzyladenine (6-BA) improves the defense system activity of waterlogged summer maize. BMC Plant Biol 2020, 20(1):44.
Wang X, Shi T, Zhao Z, Hou H, Zhang L: Proteomic analyses of sheep (ovis aries) embryonic skeletal muscle. Sci Rep 2020, 10(1):1750.
Ma J, Chen T, Wu S, Yang C, Bai M, Shu K, Li K, Zhang G, Jin Z, He F et al: iProX: an integrated proteome resource. Nucleic Acids Res 2019, 47(D1):D1211-D1217.
Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, Kong L, Gao G, Li CY, Wei L: KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res 2011, 39(Web Server issue):W316-322.
Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res 2000, 28(1):27-30.
Chen ML, Fu XM, Liu JQ, Ye TT, Hou SY, Huang YQ, Yuan BF, Wu Y, Feng YQ: Highly sensitive and quantitative profiling of acidic phytohormones using derivatization approach coupled with nano-LC-ESI-Q-TOF-MS analysis. J Chromatogr B 2012, 905:67-74.

Additionalfile1.xlsx
Additional file 1. The identified rice proteins of two genotypes at two time points after inoculation of PXO99A.
Additionalfile2.xlsx
Additional file 2. The Xanthomonas oryzae pv. oryzae (Xoo) proteins detected in two genotypes at two time points after inoculation of PXO99A.
Additionalfile3.xlsx
Additional file 3. The information of all identified peptides.
Additionalfile4.xlsx
Additional file 4. The information of all identified proteins.
Additionalfile5.tif
Additional file 5. Hierarchical cluster analysis of all identified proteins. Clustering of expressed Xoo (a) and rice (b) proteins detected in H471 and HHZ at 2 dpi and 3 dpi with PXO99A with three biological replications. Each horizontal line represents a single protein and the colored line indicates the protein abundances relative to the median in a specific sample: high accumulation in red, low accumulation in green. H471, the rice introgression line with a broad-spectrum hypersensitive reaction mediated by a new rice resistance gene Xa39; HHZ, Huang-Hua-Zhan, the recurrent parent of H471.
Additionalfile6.xlsx
Additional file 6. The differentially abundant Xanthomonas oryzae pv. oryzae (Xoo) proteins compared H741 with HHZ at two time points after inoculation of PXO99A.
Additionalfile7.xlsx
Additional file 7. The differentially abundant rice proteins compared H741 with HHZ at two time points after inoculation of PXO99A.
Additionalfile8.xlsx
Additional file 8. GO enrichment analysis of differentially abundant Xanthomonas oryzae pv. oryzae (Xoo) and rice proteins between H471 and HHZ.
Additionalfile9.xlsx
Additional file 9. KEGG pathway enrichment analysis of differentially abundant Xanthomonas oryzae pv. oryzae (Xoo) and rice proteins between H471 and HHZ.
Additionalfile10.tif
Additional file 10. Protein accumulations of OsPR1a and OsPR1b in HHZ and H471 at different time-points after inoculation with PXO99A. a OsPR1a. b OsPR1b.
Additionalfile11.tif
Additional file 11. General workflow of this study using isobaric tandem mass tag labelling. After inoculation with PXO99A, three independent biological replicates of leaf samples were taken for each genotype (H471 and HHZ) at 2 d and 3 d. After protein extraction and quantification of 12 samples (HHZ-2d-1/2/3, HHZ-3d-1/2/3, H471-2d-1/2/3, and H471-3d-1/2/3), an equal amount of 12 samples was mixed, which is referred as the reference sample. The 13 samples were prepared, digested, and labelled with TMT10plex tags in parallel. The 12 labeled samples (HHZ-2d-1/2/3, HHZ-3d-1/2/3, H471-2d-1/2/3, and H471-3d-1/2/3) were divided into pools A2 (HHZ-3d-1/2/3 and H471-2d-1/2/3) and B2 (HHZ-2d-1/2/3 and H471-3d-1/2/3). Pools A2 and B2 (30 μg protein per sample) were separately combined with 30 μg labeled reference sample protein in new microcentrifuge tubes. The two pooled protein samples were first lyophilized, reconstituted, and fractionated. Then the peptide mixtures performed LC-MS/MS analysis.
Additionalfile12.pptx
Additional file 12. All original and uncropped blot images for Figure 3 and Additional file 10. The left images are the protein abundances of OsbZIP23, OsCDPK13, OsMKK4, OsMPK6, PR1a, and PR1b in HHZ after the inoculation with PXO99A at 0d, 0.5d, 1d, 1.5d, 2d, 2.5d, 3.5d, 4d, and 4.5d respectively. The right images are the protein abundances of OsbZIP23, OsCDPK13, OsMKK4, OsMPK6, PR1a, and PR1b in H471 after the inoculation with PXO99A at 0d, 0.5d, 1d, 1.5d, 2d, 2.5d, 3.5d, 4d, and 4.5d respectively. The arrows indicate these proteins.

Download PDF

Journal Publication

published 14 Dec, 2020

Read the published version in BMC Plant Biology →

Review #1 received at journal
30 Oct, 2020
Editorial decision: Minor revision
30 Oct, 2020
Reviewer #2 agreed at journal
12 Oct, 2020
Review #2 received at journal
12 Oct, 2020
Reviewer #1 agreed at journal
11 Oct, 2020
Editor assigned by journal
09 Oct, 2020
Reviewers invited by journal
09 Oct, 2020
Submission checks completed at journal
08 Oct, 2020
Editor invited by journal
08 Oct, 2020

You are reading this older preprint version

Read the latest preprint version →

Comparative Proteomic Analysis Reveals Novel Insights into the Interaction between Rice and Xanthomonas oryzae pv. oryzae

Status:

Journal Publication

Version 2

Abstract

Figures

Background

Results

Discussion

Conclusion

Methods

Abbreviations

Declarations

References

Supplementary Files

Status:

Journal Publication

Version 2