Integrative Proteomic and Phosphoproteomic Analyses of Pattern- and Effector-Triggered Immunity in Tomato

Tomato roots exhibit distinct, development-specific responses to bacterial-derived peptides

Plants possess cell-surface recognition receptors that detect molecular patterns from microbial invaders and initiate an immune response. Understanding the conservation of pattern-triggered immunity within different plant organs and across species is crucial to its sustainable and effective use in plant disease management but is currently unclear. We examined the activation and immune response patterns of three pattern recognition receptors (PRRs: Sl FLS2, Sl FLS3, and Sl CORE) in different developmental regions of roots and in leaves of multiple accessions of domesticated and wild tomato ( Solanum lycopersicum and S. pimpinellifolium ) using biochemical and genetic assays. Roots from different tomato accessions differed in the amplitude and dynamics of their immune response, but all exhibited developmental-specific PTI responses in which the root early differentiation zone was the most sensitive to molecular patterns. PRR signaling pathways also showed distinct but occasionally overlapping responses downstream of each immune receptor in tomato roots.These results reveal that each PRR initiates a unique PTI pathway and suggest that the specificity and complexity of tomato root immunity are tightly linked to the developmental stage, emphasizing the importance of spatial and temporal regulation in PTI.

Quantitative phosphoproteomics reveals molecular pathway network in wheat resistance to stripe rust

Protein phosphorylation plays an important role in immune signaling transduction in plant resistance to pathogens. Wheat stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), severely devastates wheat production. Nonetheless, the molecular mechanism of wheat resistance to stripe rust remains limited. In this study, quantitative phosphoproteomics was employed to investigate the protein phosphorylation changes in wheat challenged by Pst. A total of 1537 and 2470 differentially accumulated phosphoproteins (DAPs) were identified from four early infection stage (6, 12, 18 and 24 h post-inoculation) in incompatible and compatible wheat-Pst interactions respectively. KEGG analysis revealed that Oxidative Phosphorylation, Phosphatidylinositol Signaling, and MAPK signaling processes are distinctively enriched in incompatible interaction, while Biosynthesis of secondary metabolites and RNA degradation process were significantly enriched in compatible interactions. In particular, abundant changes in phosphorylation levels of chloroplast proteins were identified, suggesting the regulatory role of photosynthesis in wheat-Pst interaction, which is further emphasized by protein-protein interaction (PPI) network analysis. Motif-x analysis identified [xxxxSPxxxx] motif, likely phosphorylation sites for defensive response-related kinases, and a new [xxxxSSxxxx] motif significantly enriched in incompatible interaction. The results shed light on the early phosphorylation events contributing to wheat resistance against Pst. Moreover, our study demonstrated that the phosphorylation levels of Nucleoside diphosphate kinase TaNAPK1 are upregulated at 12 hpi with CYR23 and at 24 hpi with CYR31. Transient silencing of TaNAPK1 was able to attenuate wheat resistance to CYR23 and CYR31. Our study provides new insights into the mechanisms underlying Pst-wheat interactions and may provide database to find potential targets for the development of new resistant varieties.

Phosphorylation of 399S at CsHsp70 of Cymbidium sinense is essential to maintain chlorophyll stability

The Chinese orchids symbolise nobility and gentility in China, and the variation of leaf color makes Cymbidium sinense more diversified and valuable. However, its color variations especially at the protein level still remain largely unexplored. In this study, the proteomics and phosphoproteomics of Cymbidium sinense leaf color variation mutants were studied. A total of 1059 differentially abundant proteins (DAPs) and 1127 differentially abundant phosphorylation sites belonging to 644 phosphoproteins (DAPPs) were identified in the yellow section of leaf variegation mutant of Cymbidium sinense (MY) compared with the green section (MG). Moreover, 349 co-expressing proteins were found in both omics' datasets, while only 26 proteins showed the same expression patterns in the two omics. The interaction network analysis of kinases and phosphatases showed that DAPs and DAPPs in photosynthesis, response to hormones, pigment metabolic process, phosphorylation, glucose metabolic process, and dephosphorylation might contribute to leaf color variation. The abundance of 28 Hsps and 28 phosphorylation sites belonging to 10 Hsps showed significant differences between MG and MY. CsHsp70 was selected to explore the function in Cymbidium sinense leaf variegation. The results showed CsHsp70 is essential for maintaining photosynthetic pigment content and the 399S phosphorylation site is crucial to the function of CsHsp70. Collectively, our findings construct a comprehensive coverage of protein and protein phosphorylation in leaf variegation of C. sinense, providing valuable insights into its formation mechanisms.

Effects of Enhanced Resistance and Transcriptome Analysis of Twig Blight Disease by Exogenous Brassinolide in Myrica rubra

Twig blight disease is the primary disease that affects the production of Myrica rubra in China. It was reported that exogenous brassinolide (BL) can improve disease resistance in plants. Here, we examined the effects of exogenous BL on disease resistance, chlorophyll contents, antioxidant enzyme activity, ROS accumulation, and key gene expression of M. rubra to analyze the mechanism of BR-induced resistance of twig blight disease in M. rubra. The results demonstrated that 2.0 mg·L-1 of BL could significantly lessen the severity of twig blight disease in M. rubra. Exogenous BL increased the contents of chlorophyll a, chlorophyll b, carotenoids, and total chlorophyll. Moreover, exogenous BL also significantly enhanced the activity of antioxidant enzymes such as superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT), and decreased malondialdehyde (MDA) content and reactive oxygen species (ROS) accumulation in leaves, such as H2O2 and O2·-. Additionally, exogenous BL dramatically up-regulated the expression of pathogenesis-related (PR) genes such as MrPR1, MrPR2, and MrPR10, as well as important genes such as MrBAK1, MrBRI1, and MrBZR1 involved in brassinosteroid (BR) signaling pathway. The transcriptome analysis revealed that a total of 730 common differentially expressed genes (DEGs) under BL treatment were found, and these DEGs were primarily enriched in four Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Based on these findings, nine important candidate genes related to the resistance of twig blight disease under BL treatment were further identified. In this study, we elucidated the effects of exogenous BL on enhancing the resistance of M. rubra to twig blight disease and preliminary analyzed the potential mechanism of resistance induction, which will provide a crucial foundation for the management and prevention of twig blight disease in M. rubra.

Large-scale phosphoproteome analysis in wheat seedling leaves provides evidence for extensive phosphorylation of regulatory proteins during CWMV infection

BACKGROUND

Chinese wheat mosaic virus (CWMV) often causes severe damage to wheat (Triticum aestivum L.) growth and yield. It is well known that a successful infection in plants depends on a complex interaction between the host plant and the pathogen. Post-translational modification (PTM) of proteins is considered to be one of the main processes that decides the outcome of the plant-pathogen arms race during this interaction. Although numerous studies have investigated PTM in various organisms, there has been no large-scale phosphoproteomic analysis of virus-infected wheat plants. We therefore aimed to investigate the CWMV infection-induced phosphoproteomics changes in wheat by high-resolution liquid chromatography-tandem mass spectroscopy (LC-MS/MS) using affinity-enriched peptides followed by comprehensive bioinformatics analysis.

RESULTS

Through this study, a total of 4095 phosphorylation sites have been identified in 1968 proteins, and 11.6% of the phosphorylated proteins exhibited significant changes (PSPCs) in their phosphorylation levels upon CWMV infection. The result of Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that most of the PSPCs were associated with photosynthesis, plant-pathogen interactions, and MAPK signaling pathways. The protein-protein interaction (PPI) network analysis result showed that these PSPCs were mainly participated in the regulation of biosynthesis and metabolism, protein kinase activities, and transcription factors. Furthermore, the phosphorylation levels of TaChi1 and TaP5CS, two plant immunity-related enzymes, were significantly changed upon CWMV infection, resulting in a significant decrease in CWMV accumulation in the infected plants.

CONCLUSIONS

Our results indicate that phosphorylation modification of protein plays a critical role in wheat resistance to CWMV infection. Upon CWMV infection, wheat plants will regulate the levels of extra- and intra-cellular signals and modifications of enzyme activities via protein phosphorylation. This novel information about the strategies used by wheat to resist CWMV infection will help researchers to breed new CWMV-resistant cultivars and to better understand the arms race between wheat and CWMV.

Integrative transcriptomic, proteomic, and phosphoproteomic analysis on the defense response to Magnaporthe oryzae reveals different expression patterns at the molecular level of durably resistant rice cultivar Mowanggu

Rice blast, caused by Magnaporthe oryzae is one of the most destructive diseases of rice (Oryza sativa L.) in most rice-cultivated areas worldwide. Mowanggu (MWG) is a traditional landrace rice variety in Yunnan with broad-spectrum and durable blast resistance against rice blast fungus. However, the underlying disease-resistance mechanisms remain unknown. An integrative transcriptomic, proteomic, and phosphoproteomic analysis of MWG was performed after inoculation with M. oryzae in this study. The transcriptomic and proteomic results revealed that MWG was moderately correlated at the transcriptional and protein levels. Differentially expressed genes and proteins were up-regulated and significantly enriched in protein phosphorylation, peroxisome, plant-pathogen interactions, phenylpropanoid metabolism and phenylalanine biosynthesis pathways. The phosphoproteomic profile and phosphorylated-protein-interaction network revealed that the altered phosphoproteins were primarily associated with reactive oxygen species (ROS), glycolysis, MAPK signaling pathways, and amino acid biosynthesis. In addition, a series of physiological and biochemical parameters, including ROS, soluble sugars, soluble protein and callus accumulation and defense-related enzyme activities, were used to validate the possible blast resistance mechanisms of MWG. The integrative transcriptomic, proteomic, and phosphoproteomic analysis revealed the different expression patterns at the molecular level of the durably resistant rice cultivar MWG after inoculation with M. oryzae, which provides insight into the molecular mechanisms of rice blast resistance.

Roles of long non-coding RNAs in plant immunity

Robust plant immune systems are fine-tuned by both protein-coding genes and non-coding RNAs. Long non-coding RNAs (lncRNAs) refer to RNAs with a length of more than 200 nt and usually do not have protein-coding function and do not belong to any other well-known non-coding RNA types. The non-protein-coding, low expression, and non-conservative characteristics of lncRNAs restrict their recognition. Although studies of lncRNAs in plants are in the early stage, emerging studies have shown that plants employ lncRNAs to regulate plant immunity. Moreover, in response to stresses, numerous lncRNAs are differentially expressed, which manifests the actions of low-expressed lncRNAs and makes plant-microbe/insect interactions a convenient system to study the functions of lncRNAs. Here, we summarize the current advances in plant lncRNAs, discuss their regulatory effects in different stages of plant immunity, and highlight their roles in diverse plant-microbe/insect interactions. These insights will not only strengthen our understanding of the roles and actions of lncRNAs in plant-microbe/insect interactions but also provide novel insight into plant immune responses and a basis for further research in this field.

Quantitative phosphoproteomic analysis of chitin-triggered immune responses in the plasma membrane of Arabidopsis

Plant diseases seriously damage crop production, and most plant diseases are caused by fungi. Fungal cell walls contain chitin, a highly conserved component that is widely recognised by plants as a PAMP (pathogen-associated molecular pattern) to induce defence responses. The molecular mechanisms that function downstream of chitin-triggered intracellular phosphorylation remain largely unknown. In this study, we performed quantitative phosphoproteomics analysis to study protein phosphorylation changes in the plasma membrane after chitin treatment in Arabidopsis thaliana L. seedlings. Proteins with altered phosphorylation status after chitin treatment participated in biological processes ranging from signalling, localisation, and transport, to biogenesis, processing, and metabolism, suggesting that PAMP signalling targets multiple processes to coordinate the immune response. These results provide important insights into the molecular mechanism of chitin-induced plant immunity.

Pseudomonas syringae DC3000 infection increases glucosylated N-glycans in Arabidopsis thaliana

Studying the interaction between the hemibiotrophic bacterium Pseudomonas syringae pv. tomato DC3000 and Arabidopsis thaliana has shed light onto the various forms of mechanisms plants use to defend themselves against pathogen attack. While a lot of emphasis has been put on investigating changes in protein expression in infected plants, only little information is available on the effect infection plays on the plants N-glycan composition. To close this gap in knowledge, total N-glycans were enriched from P. syringae DC3000-infected and mock treated Arabidopsis seedlings and analyzed via MALDI-TOF-MS. Additionally, fluorescently labelled N-glycans were quantified via HPLC-FLD. N-glycans from infected plants were overall less processed and displayed increased amounts of oligomannosidic N-glycans. As multiple peaks for certain oligomannosidic glycoforms were detected upon separation via liquid chromatography, a porous graphitic carbon (PGC)-analysis was conducted to separate individual N-glycan isomers. Indeed, multiple different N-glycan isomers with masses of two N-acetylhexosamine residues plus 8, 9 or 10 hexoses were detected in the infected plants which were absent in the mock controls. Treatment with jack bean α-mannosidase resulted in incomplete removal of hexoses from these N-glycans, indicating the presence of glucose residues. This hints at the accumulation of misfolded glycoproteins in the infected plants, likely because of endoplasmic reticulum (ER) stress. In addition, poly-hexose structures susceptible to α-amylase treatment were found in the DC3000-infected plants, indicating alterations in starch metabolism due to the infection process.

Proteomics of Reproductive Development, Fruit Ripening, and Stress Responses in Tomato

The fruits of the tomato crop (Solanum lycopersicum L.) are increasingly consumed by humans worldwide. Due to their rich nutritional quality, pharmaceutical properties, and flavor, tomato crops have gained a salient role as standout crops among other plants. Traditional breeding and applied functional research have made progress in varying tomato germplasms to subdue biotic and abiotic stresses. Proteomic investigations within a span of few decades have assisted in consolidating the functional genomics and transcriptomic research. However, due to the volatility and dynamicity of proteins in the regulation of various biosynthetic pathways, there is a need for continuing research in the field of proteomics to establish a network that could enable a more comprehensive understanding of tomato growth and development. With this view, we provide a comprehensive review of proteomic studies conducted on the tomato plant in past years, which will be useful for future breeders and researchers working to improve the tomato crop.

Comprehensive Analysis of NAC Genes Reveals Differential Expression Patterns in Response to Pst DC3000 and Their Overlapping Expression Pattern during PTI and ETI in Tomato

NAC (NAM/ATAF/CUC) transcription factors belong to a unique gene family in plants, which play vital roles in regulating diverse biological processes, including growth, development, senescence, and in response to biotic and abiotic stresses. Tomato (Solanum lycopersicum), as the most highly valued vegetable and fruit crop worldwide, is constantly attacked by Pseudomonas syringae pv. tomato DC3000 (Pst DC3000), causing huge losses in production. Thus, it is essential to conduct a comprehensive identification of the SlNAC genes involved in response to Pst DC3000 in tomato. In this study, a complete overview of this gene family in tomato is presented, including genome localization, protein domain architectures, physical and chemical features, and nuclear location score. Phylogenetic analysis identified 20 SlNAC genes as putative stress-responsive genes, named SSlNAC 1–20. Expression profiles analysis revealed that 18 of these 20 SSlNAC genes were significantly induced in defense response to Pst DC3000 stress. Furthermore, the RNA-seq data were mined and analyzed, and the results revealed the expression pattern of the 20 SSlNAC genes in response to Pst DC3000 during the PTI and ETI. Among them, SSlNAC3, SSlNAC4, SSlNAC7, SSlNAC8, SSlNAC12, SSlNAC17, and SSlNAC19 were up-regulated against Pst DC3000 during PTI and ETI, which suggested that these genes may participate in both the PTI and ETI pathway during the interaction between tomato and Pst DC3000. In addition, SSlNAC genes induced by exogenous hormones, including indole-3-acetic acid (IAA), abscisic acid (ABA), salicylic acid (SA), and methyl jasmonic acid (MeJA), were also recovered. These results implied that SSlNAC genes may participate in the Pst DC3000 stress response by multiple regulatory pathways of the phytohormones. In all, this study provides important clues for further functional analysis and of the regulatory mechanism of SSlNAC genes under Pst DC3000 stress.

Integrated volatile metabolomic and transcriptomic analysis provides insights into the regulation of floral scents between two contrasting varieties of Lonicera japonica

Lonicera japonica Thunb., belonging to the Caprifoliaceae family, is an important traditional Chinese medicinal plant. The L. japonica flower (LJF) is widely used in medicine, cosmetics, drinks, and food due to its medicinal and sweet-smelling properties. Considerable efforts have been devoted to investigating the pharmacological activities of LJF; however, the regulatory mechanism of the floral scents remains unknown. We previously selected and bred an elite variety of L. japonica var. chinensis Thunb. called 'Yujin2', which has a strong aroma and is used in functional drinks and cosmetics. In order to reveal the regulatory mechanism of the floral scents of LJF, volatile metabolomic and transcriptomic analyses of the LJF at the silver flowering stage of 'Yujin2' (strong aroma) and 'Fengjin1' (bland odor) were performed. Our results revealed that a total of 153 metabolites and 9,523 genes were differentially regulated in LJF between 'Yujin2' and 'Fengjin1'. The integrated analysis of omics data indicated that the biosynthetic pathways of terpenoids (i.e., monoterpenoids, including geraniol and alpha-terpineol; sesquiterpenoids, including farnesol, farnesal, and alpha-farnesene; triterpenoid squalene), tryptophan and its derivatives (methyl anthranilate), and fatty acid derivatives, were major contributors to the stronger aroma of 'Yujin2' compared to 'Fengjin1'. Moreover, several genes involved in the terpenoid biosynthetic pathway were characterized using quantitative real-time PCR. These results provide insights into the metabolic mechanisms and molecular basis of floral scents in LJF, enabling future screening of genes related to the floral scent regulation, such as alpha-terpineol synthase, geranylgeranyl diphosphate synthase, farnesyl pyrophosphate synthase, anthranilate synthase, as well as transcription factors such as MYB, WRKY, and LFY. The knowledge from this study will facilitate the breeding of quality-improved and more fragrant variety of L. japonica for ornamental purpose and functional beverages and cosmetics.

Protein glycosylation changes during systemic acquired resistance in Arabidopsis thaliana

N-glycosylation, an important post-translational modification of proteins in all eukaryotes, has been clearly shown to be involved in numerous diseases in mammalian systems. In contrast, little is known regarding the role of protein N-glycosylation in plant defensive responses to pathogen infection. We identified, for the first time, glycoproteins related to systemic acquired resistance (SAR) in an Arabidopsis thaliana model, using a glycoproteomics platform based on high-resolution mass spectrometry. 407 glycosylation sites corresponding to 378 glycopeptides and 273 unique glycoproteins were identified. 65 significantly changed glycoproteins with 80 N-glycosylation sites were detected in systemic leaves of SAR-induced plants, including numerous GDSL-like lipases, thioglucoside glucohydrolases, kinases, and glycosidases. Functional enrichment analysis revealed that significantly changed glycoproteins were involved mainly in N-glycan biosynthesis and degradation, phenylpropanoid biosynthesis, cutin and wax biosynthesis, and plant-pathogen interactions. Comparative analysis of glycoproteomics and proteomics data indicated that glycoproteomics analysis is an efficient method for screening proteins associated with SAR. The present findings clarify glycosylation status and sites of A. thaliana proteins, and will facilitate further research on roles of glycoproteins in SAR induction.