Proteome Profiling of Paulownia Seedlings Infected with Phytoplasma

Paulownia Witches' Broom Disease: A Comprehensive Review

Phytoplasmas are insect-transmitted bacterial pathogens associated with diseases in a wide range of host plants, resulting in significant economic and ecological losses. Perennial deciduous trees in the genus Paulownia are widely planted for wood harvesting and ornamental purposes. Paulownia witches' broom (PaWB) disease, associated with a 16SrI-D subgroup phytoplasma, is a destructive disease of paulownia in East Asia. The PaWB phytoplasmas are mainly transmitted by insect vectors in the Pentatomidae (stink bugs), Miridae (mirid bugs) and Cicadellidae (leafhoppers) families. Diseased trees show typical symptoms, such as branch and shoot proliferation, which together are referred to as witches' broom. The phytoplasma presence affects the physiological and anatomical structures of paulownia. Gene expression in paulownia responding to phytoplasma presence have been studied at the transcriptional, post-transcriptional, translational and post-translational levels by high throughput sequencing techniques. A PaWB pathogenic mechanism frame diagram on molecular level is summarized. Studies on the interactions among the phytoplasma, the insect vectors and the plant host, including the mechanisms underlying how paulownia effectors modify processes of gene expression, will lead to a deeper understanding of the pathogenic mechanisms and to the development of efficient control measures.

Comparative proteome analysis reveals the role of negative floral regulators and defense-related genes in phytoplasma infected sesame

"Candidatus Phytoplasma australiense" is associated with floral malformations in sesame but the interaction remains largely unexplored. A label-free quantitative shotgun proteomics approach through liquid chromatography-mass spectrometry quadruple time-of-flight was used to analyze changes in the proteome of asymptomatic (control) and symptomatic (phytoplasma-infected) sesame plants to identify proteins differentially expressed during phytoplasma infection at early stages of flower development. A total of 3457 and 1704 proteins were identified from asymptomatic and symptomatic samples respectively through proteome profiling with three runs per sample. Several differentially abundant proteins (DAPs) were identified which might be involved in sesame-phytoplasma interaction. The DAPs identified were related to transcription, cell division, chromosome partitioning, defense mechanisms, negative regulation of flower development, amino acid transport and metabolism, signal transduction and RNA processing, and its modifications. Of these proteins, 21 were downregulated while 212 were significantly upregulated in symptomatic sesame plants compared to the control plants. The floral development-related proteins like UBP16 and DCAF1 were found to be downregulated while negative regulators/repressors of floral development genes, HUA2, PIE1, and ICU2, were upregulated in symptomatic samples indicating phytoplasma's role in altering the expression of these genes. Validation of these genes through quantitative retro-transcripted PCR suggested that the DAPs observed in symptomatic sesame might be induced by phytoplasma presence to suppress flowering via negative regulation of flower development.

Increased susceptibility to Chrysanthemum Yellows phytoplasma infection in Atcals7ko plants is accompanied by enhanced expression of carbohydrate transporters

Loss of CALS7 appears to confer increased susceptibility to phytoplasma infection in Arabidopsis, altering expression of genes involved in sugar metabolism and membrane transport. Callose deposition around sieve pores, under control of callose synthase 7 (CALS7), has been interpreted as a mechanical response to limit pathogen spread in phytoplasma-infected plants. Wild-type and Atcals7ko mutants were, therefore, employed to unveil the mode of involvement of CALS7 in the plant's response to phytoplasma infection. The fresh weights of healthy and CY-(Chrysanthemum Yellows) phytoplasma-infected Arabidopsis wild type and mutant plants indicated two superimposed effects of the absence of CALS7: a partial impairment of photo-assimilate transport and a stimulated phytoplasma proliferation as illustrated by a significantly increased phytoplasma titre in Atcal7ko mutants. Further studies solely dealt with the effects of CALS7 absence on phytoplasma growth. Phytoplasma infection affected sieve-element substructure to a larger extent in mutants than in wild-type plants, which was also true for the levels of some free carbohydrates. Moreover, infection induced a similar upregulation of gene expression of enzymes involved in sucrose cleavage (AtSUS5, AtSUS6) and transmembrane transport (AtSWEET11) in mutants and wild-type plants, but an increased gene expression of carbohydrate transmembrane transporters (AtSWEET12, AtSTP13, AtSUC3) in infected mutants only. It remains still unclear how the absence of AtCALS7 leads to gene upregulation and how an increased intercellular mobility of carbohydrates and possibly effectors contributes to a higher susceptibility. It is also unclear if modified sieve-pore structures in mutants allow a better spread of phytoplasmas giving rise to higher titre.

Conformational Changes in Three-Dimensional Chromatin Structure in Paulownia fortunei After Phytoplasma Infection

Higher-order chromatin structures play important roles in regulating multiple biological processes such as growth and development as well as biotic and abiotic stress response. However, little is known about three-dimensional chromatin structures in Paulownia or about whole-genome chromatin conformational changes that occur in response to Paulownia witches' broom (PaWB) disease. We used high-throughput chromosome conformation capture (Hi-C) to obtain genome-wide profiles of chromatin conformation in both healthy and phytoplasma-infected Paulownia fortunei genome. The heat map results indicated that the strongest interactions between chromosomes were in the telomeres. We confirmed that the main structural characteristics of A/B compartments, topologically associated domains, and chromatin loops were prominent in the Paulownia genome and were clearly altered in phytoplasma-infected plants. By combining chromatin immunoprecipitation sequencing, Hi-C signals, and RNA sequencing data, we inferred that the chromatin structure changed and the modification levels of three histones (H3K4me3/K9ac/K36me3) increased in phytoplasma-infected P. fortunei, which was associated with changes of transcriptional activity. We concluded that for epigenetic modifications, transcriptional activity might function in combination to shape chromatin packing in healthy and phytoplasm-infected Paulownia. Finally, 11 genes (e.g., RPN6, Sec61 subunit-α) that were commonly located at specific topologically associated domain boundaries, A/B compartment switching and specific loops, and had been associated with histone marks were identified and considered as closely related to PaWB stress. Our results provide new insights into the nexus between gene regulation and chromatin conformational alterations in nonmodel plants upon phytopathogen infection and plant disease resistance.

Discovery of potential protein biomarkers associated with sugarcane white leaf disease susceptibility using a comparative proteomic approach

Sugarcane white leaf disease (SCWLD) is caused by phytoplasma, a serious sugarcane phytoplasma pathogen, which causes significant decreases in crop yield and sugar quality. The identification of proteins involved in the defense mechanism against SCWLD phytoplasma may help towards the development of varieties resistant to SCWLD. We investigated the proteomes of four sugarcane varieties with different levels of susceptibility to SCWLD phytoplasma infection, namely K88-92 and K95-84 (high), KK3 (moderate), and UT1 (low) by quantitative label-free nano-liquid chromatography-tandem mass spectrometry (nano LC-MS/MS). A total of 248 proteins were identified and compared among the four sugarcane varieties. Two potential candidate protein biomarkers for reduced susceptibility to SCWLD phytoplasma were identified as proteins detected only in UT1. The functions of these proteins are associated with protein folding, metal ion binding, and oxidoreductase. The candidate biomarkers could be useful for further study of the sugarcane defense mechanism against SCWLD phytoplasma, and in molecular and conventional breeding strategies for variety improvement.

Integration of metabolomics and existing omics data reveals new insights into phytoplasma-induced metabolic reprogramming in host plants

Phytoplasmas are cell wall-less bacteria that induce abnormal plant growth and various diseases, causing severe economic loss. Phytoplasmas are highly dependent on nutrients imported from host cells because they have lost many genes involved in essential metabolic pathways during reductive evolution. However, metabolic crosstalk between phytoplasmas and host plants and the mechanisms of phytoplasma nutrient acquisition remain poorly understood. In this study, using metabolomics approach, sweet cherry virescence (SCV) phytoplasma-induced metabolite alterations in sweet cherry trees were investigated. A total of 676 metabolites were identified in SCV phytoplasma-infected and mock inoculated leaves, of which 187 metabolites were differentially expressed, with an overwhelming majority belonging to carbohydrates, fatty acids/lipids, amino acids, and flavonoids. Available omics data of interactions between plant and phytoplasma were also deciphered and integrated into the present study. The results demonstrated that phytoplasma infection promoted glycolysis and pentose phosphate pathway activities, which provide energy and nutrients, and facilitate biosynthesis of necessary low-molecular metabolites. Our findings indicated that phytoplasma can induce reprograming of plant metabolism to obtain nutrients for its own replication and infection. The findings from this study provide new insight into interactions of host plants and phytoplasmas from a nutrient acquisition perspective.

Integrated dual RNA-seq and dual iTRAQ of infected tissue reveals the functions of a diguanylate cyclase gene of Pseudomonas plecoglossicida in host-pathogen interactions with Epinephelus coioides

The interactions between host and pathogen is exceedingly complex, which involves alterations at multiple molecular layers. However, research to simultaneously monitor the alterations of transcriptome and proteome between a bacterial pathogen and aquatic animal host through integrated dual RNA-seq and dual iTRAQ of tissue during infection is currently lacking. The important role of a diguanylate cyclase gene (L321_RS15240) in pathogenicity of Pseudomonas plecoglossicida against Epinephelus coioides was suggested by previous dual RNA-seq of our lab. Then L321_RS15240-RNAi strains of P. plecoglossicida were constructed with pCM130/tac, and the mutant with the best silencing effect was selected for follow-up study. The RNAi of L321_RS15240 resulted in a significant decrease in bacterial virulence of P. plecoglossicida. The E. coioides spleens infected by wild type strain or L321_RS15240-RNAi strain of P. plecoglossicida were subjected to dual RNA-seq and dual iTRAQ, respectively. The results showed that: RNAi of L321_RS15240 led to 1)alterations of host transcriptome associated with complement and coagulation cascades, ribosome, arginine and proline metabolism, and oxidative phosphorylation; 2)high expression of host proteins which related to phagosome and metabolism responses (metabolism of glutathione, amino sugar and nucleotide sugar); 3)the highly differentially expression of host lncRNAs and miRNAs. The differentially expressed proteins and mRNAs of pathogen were different after infection, but the functions of these proteins and mRNAs were mainly related to metabolism and virulence. This study provides a new insight to comprehensively understand the gene functions of pathogens and hosts at multiple molecular layers during in vivo infection.

Phytoplasma-induced Changes in the Acetylome and Succinylome of Paulownia tomentosa Provide Evidence for Involvement of Acetylated Proteins in Witches' Broom Disease

Lysine acetylation and succinylation are post-translational modifications of proteins that have been shown to play roles in plants response to pathogen infection. Phytoplasma infection can directly alter multiple metabolic processes in the deciduous plant Paulownia and lead to Paulownia witches' broom (PaWB) disease, the major cause of Paulownia mortality worldwide. However, the extent and function of lysine aceylation and succinylation during phytoplasma infection have yet to be explored. Here, we investigated the changes in the proteome, acetylome, and succinylome of phytoplasma-infected Paulownia tomentosa seedlings using quantitative mass spectrometry. In total, we identified 8963 proteins, 2893 acetylated proteins (5558 acetylation sites), and 1271 succinylated proteins (1970 succinylation sites), with 425 (533 sites) simultaneously acetylated and succinylated. Comparative analysis revealed that 276 proteins, 546 acetylated proteins (741 acetylation sites) and 5 succinylated proteins (5 succinylation sites) were regulated in response to phytoplasma infection, suggesting that acetylation may be more important than succinylation in PaWB. Enzymatic assays showed that acetylation of specific sites in protochlorophyllide reductase and RuBisCO, key enzymes in chlorophyll and starch biosynthesis, respectively, modifies their activity in phytoplasma-infected seedlings. On the basis of these results, we propose a model to elucidate the molecular mechanism of responses to PaWB and offer a resource for functional studies on the effects of acetylation on protein function.

iTRAQ-Based Quantitative Proteomics and Transcriptomics Provide Insights Into the Importance of Expansins During Root Development in Carrot

Carrot is an important root vegetable crop with a variety of nutrients. As the main product of carrots, the growth and development of fleshy roots directly determine the yield and quality of carrots. However, molecular mechanism underlying the carrot root formation and expansion is still limited. In our study, isobaric tags for relative and absolute quantification (iTRAQ) was utilized to explore the differentially expressed proteins (DEPs) during different developmental stages of carrot roots. Overall, 2,845 proteins were detected, of which 118 were significantly expressed in all three stages. DEPs that participated in several growth metabolisms were identified, including energy metabolism, defense metabolism, cell growth and shape regulation. Among them, two expansin proteins were obtained. A total of 30 expansin genes were identified based on the carrot genome database. Structure analysis showed that carrot expansin gene family was relatively conserved. Based on the expression analysis, we found that the expression profile of expansins genes was up-regulated during the vigorous growing period of carrot root. Furthermore, there was a consistent relationship between the expression patterns of mRNA and protein. The results indicated that expansin proteins might play important roles during root development in carrot. Our work provided useful information for understanding molecular mechanism of carrot root development.

Genome-wide analysis of three histone marks and gene expression in Paulownia fortunei with phytoplasma infection

BACKGROUND

Paulownia withes'-broom (PaWB) disease caused by phytoplasma is a serious infectious disease for Paulownia. However, the underlying molecular pathogenesis is not fully understood. Recent studies have demonstrated that histone modifications could play a role in plant defense responses to pathogens. But there is still no available genome-wide histone modification data in non-model ligneous species infected with phytoplasma.

RESULTS

Here, we provided the first genome-wide profiles of three histone marks (H3K4me3, H3K36me3 and H3K9ac) in Paulownia fortunei under phytoplasma stress by using chromatin immunoprecipitation sequencing (ChIP-Seq). We found that H3K4me3, H3K36me3 and H3K9ac were mainly enriched in the genic regions in P. fortunei with (PFI) and without (PF) phytoplasma infection. ChIP-Seq analysis revealed 1738, 986, and 2577 genes were differentially modified by H3K4me3, H3K36me3 and H3K9ac marks in PFI under phytoplasma infection, respectively. The functional analysis of these genes suggested that most of them were mainly involved in metabolic pathways, biosynthesis of secondary metabolites, phenylpropanoid biosynthesis, plant-pathogen interaction and plant hormone signal transduction. In addition, the combinational analysis of ChIP-Seq and RNA-Seq showed that differential histone methylation and acetylation only affected a small subset of phytoplasma-responsive genes.

CONCLUSIONS

Taken together, this is the first report of integrated analysis of histone modifications and gene expression involved in Paulownia-phytoplasma interaction. Our results will provide the valuable resources for the mechanism studies of gene regulation in non-model plants upon pathogens attack.

Proteomics, Holm Oak (Quercus ilex L.) and Other Recalcitrant and Orphan Forest Tree Species: How do They See Each Other?

Proteomics has had a big impact on plant biology, considered as a valuable tool for several forest species, such as Quercus, Pines, Poplars, and Eucalyptus. This review assesses the potential and limitations of the proteomics approaches and is focused on Quercus ilex as a model species and other forest tree species. Proteomics has been used with Q. ilex since 2003 with the main aim of examining natural variability, developmental processes, and responses to biotic and abiotic stresses as in other species of the genus Quercus or Pinus. As with the progress in techniques in proteomics in other plant species, the research in Q. ilex moved from 2-DE based strategy to the latest gel-free shotgun workflows. Experimental design, protein extraction, mass spectrometric analysis, confidence levels of qualitative and quantitative proteomics data, and their interpretation are a true challenge with relation to forest tree species due to their extreme orphan and recalcitrant (non-orthodox) nature. Implementing a systems biology approach, it is time to validate proteomics data using complementary techniques and integrate it with the -omics and classical approaches. The full potential of the protein field in plant research is quite far from being entirely exploited. However, despite the methodological limitations present in proteomics, there is no doubt that this discipline has contributed to deeper knowledge of plant biology and, currently, is increasingly employed for translational purposes.

ceRNA Cross-Talk in Paulownia Witches' Broom Disease

Long noncoding RNA (lncRNA), circular RNA (circRNA), and microRNA (miRNA) are important in the regulation of life activities. However, their function is unclear in Paulownia fortunei. To identify lncRNAs, circRNAs, and miRNA, and investigate their roles in the infection progress of Paulownia witches’ broom (PaWB) disease, we performed RNA sequencing of healthy and infected P. fortunei. A total of 3126 lncRNAs, 1634 circRNAs, and 550 miRNAs were identified. Among them, 229 lncRNAs, 65 circRNAs, and 65 miRNAs were differentially expressed in a significant manner. We constructed a competing endogenous RNA (ceRNA) network, which contains 5 miRNAs, 4 circRNAs, 5 lncRNAs, and 15 mRNAs, all of which were differentially expressed between healthy and infected P. fortunei. This study provides the first catalog of candidate ceRNAs in Paulownia and gives a revealing insight into the molecular mechanism responsible for PaWB.

Regulation of Long Noncoding RNAs Responsive to Phytoplasma Infection in Paulownia tomentosa

Paulownia witches' broom caused by phytoplasma infection affects the production of Paulownia trees worldwide. Emerging evidence showed that long noncoding RNAs (lncRNA) play a protagonist role in regulating the expression of genes in plants. So far, the identification of lncRNAs has been limited to a few model plant species, and their roles in mediating responses to Paulownia tomentosa that free of phytoplasma infection are yet to be characterized. Here, whole-genome identification of lncRNAs, based on strand-specific RNA sequencing, from four Paulownia tomentosa samples, was performed and identified 3689 lncRNAs. These lncRNAs showed low conservation among plant species and some of them were miRNA precursors. Further analysis revealed that the 112 identified lncRNAs were related to phytoplasma infection. We predicted the target genes of these phytoplasma-responsive lncRNAs, and our analysis showed that 51 of the predicted target genes were alternatively spliced. Moreover, we found the expression of the lncRNAs plays vital roles in regulating the genes involved in the reactive oxygen species induced hypersensitive response and effector-triggered immunity in phytoplasma-infected Paulownia. This study indicated that diverse sets of lncRNAs were responsive to Paulownia witches' broom, and the results will provide a starting point to understand the functions and regulatory mechanisms of Paulownia lncRNAs in the future.

Genome-wide expression analysis of salt-stressed diploid and autotetraploid Paulownia tomentosa

Paulownia tomentosa is a fast-growing tree species with multiple uses. It is grown worldwide, but is native to China, where it is widely cultivated in saline regions. We previously confirmed that autotetraploid P. tomentosa plants are more stress-tolerant than the diploid plants. However, the molecular mechanism underlying P. tomentosa salinity tolerance has not been fully characterized. Using the complete Paulownia fortunei genome as a reference, we applied next-generation RNA-sequencing technology to analyze the effects of salt stress on diploid and autotetraploid P. tomentosa plants. We generated 175 million clean reads and identified 15,873 differentially expressed genes (DEGs) from four P. tomentosa libraries (two diploid and two autotetraploid). Functional annotations of the differentially expressed genes using the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases revealed that plant hormone signal transduction and photosynthetic activities are vital for plant responses to high-salt conditions. We also identified several transcription factors, including members of the AP2/EREBP, bHLH, MYB, and NAC families. Quantitative real-time PCR analysis validated the expression patterns of eight differentially expressed genes. Our findings and the generated transcriptome data may help to accelerate the genetic improvement of cultivated P. tomentosa and other plant species for enhanced growth in saline soils.