Shotgun proteomic analysis of the Mexican lime tree infected with "CandidatusPhytoplasma aurantifolia"

A newly identified glycosyltransferase AsRCOM provides resistance to purple curl leaf disease in agave

BACKGROUND

Purple curl leaf disease brings a significant threat to the development of agave industry, the underlying mechanism of disease-resistant Agave sisalana. hybrid 11648 (A. H11648R) is still unknown.

RESULTS

To excavate the crucial disease-resistant genes against purple curl leaf disease, we performed an RNA-seq analysis for A.H11648R and A.H11648 during different stages of purple curl leaf disease. The DEGs (differentially expressed genes) were mainly enriched in linolenic acid metabolism, starch and sucrose mechanism, phenylpropanoid biosynthesis, hypersensitive response (HR) and systemic acquired resistance. Further analysis suggested that eight candidate genes (4'OMT2, ACLY, NCS1, GTE10, SMO2, FLS2, SQE1 and RCOM) identified by WGCNA (weighted gene co-expression network analysis) may mediate the resistance to agave purple curl disease by participating the biosynthesis of benzylisoquinoline alkaloids, steroid, sterols and flavonoids, and the regulation of plant innate immunity and systemic acquired resistance. After qPCR verification, we found that AsRCOM, coding a glycosyltransferase and relevant to the regulation of plant innate immunity and systemic acquired resistance, may be the most critical disease-resistant gene. Finally, the overexpression of AsRCOM gene in agave could significantly enhance the resistance to purple curl disease with abundant reactive oxygen species (ROS) accumulations.

CONCLUSIONS

Integrative RNA-seq analysis found that HR may be an important pathway affecting the resistance to purple curl leaf disease in agave, and identified glycosyltransferase AsRCOM as the crucial gene that could significantly enhance the resistance to purple curl leaf disease in agave, with obvious ROS accumulations.

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.

New Cross-Talks between Pathways Involved in Grapevine Infection with 'Candidatus Phytoplasma solani' Revealed by Temporal Network Modelling

Understanding temporal biological phenomena is a challenging task that can be approached using network analysis. Here, we explored whether network reconstruction can be used to better understand the temporal dynamics of bois noir, which is associated with 'Candidatus Phytoplasma solani', and is one of the most widespread phytoplasma diseases of grapevine in Europe. We proposed a methodology that explores the temporal network dynamics at the community level, i.e., densely connected subnetworks. The methodology offers both insights into the functional dynamics via enrichment analysis at the community level, and analyses of the community dissipation, as a measure that accounts for community degradation. We validated this methodology with cases on experimental temporal expression data of uninfected grapevines and grapevines infected with 'Ca. P. solani'. These data confirm some known gene communities involved in this infection. They also reveal several new gene communities and their potential regulatory networks that have not been linked to 'Ca. P. solani' to date. To confirm the capabilities of the proposed method, selected predictions were empirically evaluated.

Differential Response of Grapevine to Infection with 'Candidatus Phytoplasma solani' in Early and Late Growing Season through Complex Regulation of mRNA and Small RNA Transcriptomes

Bois noir is the most widespread phytoplasma grapevine disease in Europe. It is associated with ‘Candidatus Phytoplasma solani’, but molecular interactions between the causal pathogen and its host plant are not well understood. In this work, we combined the analysis of high-throughput RNA-Seq and sRNA-Seq data with interaction network analysis for finding new cross-talks among pathways involved in infection of grapevine cv. Zweigelt with ‘Ca. P. solani’ in early and late growing seasons. While the early growing season was very dynamic at the transcriptional level in asymptomatic grapevines, the regulation at the level of small RNAs was more pronounced later in the season when symptoms developed in infected grapevines. Most differentially expressed small RNAs were associated with biotic stress. Our study also exposes the less-studied role of hormones in disease development and shows that hormonal balance was already perturbed before symptoms development in infected grapevines. Analysis at the level of communities of genes and mRNA-microRNA interaction networks revealed several new genes (e.g., expansins and cryptdin) that have not been associated with phytoplasma pathogenicity previously. These novel actors may present a new reference framework for research and diagnostics of phytoplasma diseases of grapevine.

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.

Longitudinal Transcriptomic, Proteomic, and Metabolomic Analyses of Citrus sinensis (L.) Osbeck Graft-Inoculated with "Candidatus Liberibacter asiaticus"

"Candidatus Liberibacter asiaticus" (CLas) is the bacterium associated with the citrus disease Huanglongbing (HLB). Current CLas detection methods are unreliable during presymptomatic infection, and understanding CLas pathogenicity to help develop new detection techniques is challenging because CLas has yet to be isolated in pure culture. To understand how CLas affects citrus metabolism and whether infected plants produce systemic signals that can be used to develop improved detection techniques, leaves from Washington Navel orange (Citrus sinensis (L.) Osbeck) plants were graft-inoculated with CLas and longitudinally studied using transcriptomics (RNA sequencing), proteomics (liquid chromatography-tandem mass spectrometry), and metabolomics (proton nuclear magnetic resonance). Photosynthesis gene expression and protein levels were lower in infected plants compared to controls during late infection, and lower levels of photosynthesis proteins were identified as early as 8 weeks post-grafting. These changes coordinated with higher sugar concentrations, which have been shown to accumulate during HLB. Cell wall modification and degradation gene expression and proteins were higher in infected plants during late infection. Changes in gene expression and proteins related to plant defense were observed in infected plants as early as 8 weeks post-grafting. These results reveal coordinated changes in greenhouse navel leaves during CLas infection at the transcript, protein, and metabolite levels, which can inform of biomarkers of early infection.

Proteomic and metabolomic analysis of desiccation tolerance in wheat young seedlings

Young wheat seedlings are desiccation tolerant and have the capacity to withstand long dehydration period. In this study, we characterized the proteome and metabolome of wheat seedlings during desiccation and after recovery. Functional classification of differentially identified proteins revealed dynamic changes in the number and abundance of proteins observed during stress and recovery. Desiccation resulted in a decline in the abundance of proteins associated with photosynthesis and carbohydrate reserves, along with an increase in the presence of proteins associated with stress and defense response, such as peroxiredoxins and antioxidant enzymes. Following recovery, the abundance of stress-responsive proteins returned either partially or completely to their baseline level, confirming their importance to the seedling's desiccation response. Furthermore, proteins involved in carbohydrate metabolism, as well as fructose-bisphosphate aldolase and fructokinase-2 and phosphorylated metabolites as the substrate or the end-product, showed the inverse pattern during desiccation and after re-watering. This may reflect the fact that plants maintained energy supply during stress to protect seedlings from further damage, and for use in subsequent recovery after rewatering period. This study provides novel insights into the molecular mechanisms underlying the desiccation tolerance of wheat seedlings, and paves the way for more detailed molecular analysis of this remarkable phenomenon.

Comparative proteomics of two citrus varieties in response to infection by the fungus Alternaria alternata

Alternaria brown spot (ABS) is a disease caused by the necrotrophic fungus Alternaria alternata, which induces necrotic lesions on fruits and young leaves due to the production of the host-specific ACT toxin by the fungus. To better understand the citrus-A. alternata interaction and to identify putative resistance proteins, as well as the receptor of the ACT toxin, citrus plants susceptible ('Minneola' mandarin) and resistant ('Clemenules' tangor) to A. alternata, infected or not (control) with the pathogen were analyzed by proteomics. Protein changes were observed between citrus genotypes after infection, and 150 candidate proteins were obtained. A general scheme of the metabolic processes involved in susceptible and resistant citrus-A. alternata interactions was designed. Susceptible plants presented a high level of proteins involved in stress response at the final stages of the infection, whereas resistant plants presented high level of ROS proteins, metabolic proteins, and proteins involved in the immune system process. Proteins like ferredoxin and cyclophilin are specific to the susceptible variety and may be good candidates as fungal effector-interacting proteins. This is the first citrus-A. alternata proteomics analysis, which has allowed a better understanding of the molecular bases of the citrus response to ABS disease.

Plant Hormones in Phytoplasma Infected Plants

Phytoplasmas are bacterial plant pathogens that need a plant host and an insect vector for their spread and survival. In plants, the physiological responses that phytoplasmas trigger result in symptom development through effects on hormonal, nutritional, and stress signaling pathways, and the interactions between these. In this review, recent advances on the involvement of plant hormones together with their known and deduced roles in plants infected with phytoplasmas are discussed. Several studies have directly, or in many cases indirectly, addressed plant hormone systems in phytoplasma-infected plants. These have provided accumulating evidence that phytoplasmas extensively affect plant hormone pathways. Phytoplasmas thus, with disturbing complex plant hormone networks, suppress plant immunity and modify plant structure, while optimizing their nutrient acquisition and facilitating their colonization of the plants, and their dissemination among plants by their insect vectors.

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.

Label-free and isobaric tandem mass tag (TMT) multiplexed quantitative proteomic data of two contrasting rice cultivars exposed to drought stress and recovery

Two rice cultivars, IAC1131 (drought tolerant) and Nipponbare (drought sensitive), with contrasting genetic backgrounds and levels of tolerance to drought, were analysed using both label-free and tandem mass tags (TMTs) quantitative proteomics approaches, aiming to elucidate the mechanisms of drought tolerance. Four-week-old seedlings of both cultivars were grown in large soil volumes in the glasshouse under controlled conditions and then exposed to moderate and extreme drought for 7 days, followed by 3 days of re-watering period. Mature leaves were harvested from plants of each treatment for protein extraction and subsequent complementary shotgun proteomic analyses. The data from this study are related to the research article “Quantitative proteomic analysis of two different rice varieties reveals that drought tolerance is correlated with reduced abundance of photosynthetic machinery and increased abundance of ClpD1 protease” (Wu et al., 2016) [1].

Comparative genome analysis of jujube witches'-broom Phytoplasma, an obligate pathogen that causes jujube witches'-broom disease

BACKGROUND

JWB phytoplasma is a kind of insect-transmitted and uncultivable bacterial plant pathogen causeing a destructive Jujube disease. To date, no genome information about JWB phytoplasma has been published, which hindered its characterization at genomic level. To understand its pathogenicity and ecology, the genome of a JWB phytoplasma isolate jwb-nky was sequenced and compared with other phytoplasmas enabled us to explore the mechanisms of genomic rearrangement.

RESULTS

The complete genome sequence of JWB phytoplasma (jwb-nky) was determined, which consisting of one circular chromosome of 750,803 bp with a GC content of 23.3%. 694 protein-encoding genes, 2 operons for rRNA genes and 31 tRNA genes as well as 4 potential mobile units (PMUs) containing clusters of DNA repeats were identified. Based on PHIbaes analysis, a large number of genes were genome-specific and approximately 13% of JWB phytoplasma genes were predicted to be associated with virulence. Although transporters for maltose, dipeptides/oligopeptides, spermidine/putrescine, cobalt, Mn/Zn and methionine were identified, KEGG pathway analysis revealed the reduced metabolic capabilities of JWB phytoplasma. Comparative genome analyses between JWB phytoplasma and other phytoplasmas shows the occurrence of large-scale gene rearrangements. The low synteny with other phytoplasmas indicated that the expansion of multiple gene families/duplication probably occurred separately after differentiation.

CONCLUSIONS

In this study, the complete genome sequence of a JWB phytoplasma isolate jwb-nky that causing JWB disease was reported for the first time and a number of species-specific genes were identified in the genome. The study enhanced our understandings about genomic basis and the pathogenicity mechanism of this pathogen, which will aid in the development of improved strategies for efficient management of JWB diseases.

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.

Growth and protease secretion of Scedosporium aurantiacum under conditions of hypoxia

One of the micro-environmental stresses that fungal pathogens, such as Scedosporium aurantiacum, colonising human lungs encounter in vivo is hypoxia, or deficiency of oxygen. In this work, we studied the impacts of a hypoxic micro-environment (oxygen levels ≤1%) on the growth of a clinical S. aurantiacum isolate (WM 06.482; CBS 136046) and an environmental strain (S. aurantiacum WM 10.136; CBS 136049) on mucin-containing synthetic cystic fibrosis sputum medium. Additionally, profiles of secreted proteases were compared between the two isolates and protease activity was assessed using class-specific substrates and inhibitors. Overall, both isolates grew slower and produced less biomass under hypoxia compared to normoxic conditions. The pH of the medium decreased to 4.0 over the cultivation time, indicating that S. aurantiacum released acidic compounds into the medium. Accordingly, secreted proteases of the two isolates were dominated by acidic proteases, including aspartic and cysteine proteases, with optimal protease activity at pH 4.0 and 6.0 respectively. The clinical isolate produced higher aspartic and cysteine protease activities. Conversely, all serine proteases, including elastase-like, trypsin-like, chymotrypsin-like and subtilisin-like proteases had higher activities in the environmental isolate. Sequence similarities to 13 secreted proteases were identified by mass spectrometry (MS) by searching against other fungal proteases in the NCBI database. Results from MS analysis were consistent with those from activity assays. The clinical highly-virulent, and environmental low-virulence S. aurantiacum isolates responded differently to hypoxia in terms of the type of proteases secreted, which may reflect their different virulence properties.

Comparative assessment of the Euglena gracilis var. saccharophila variant strain as a producer of the β-1,3-glucan paramylon under varying light conditions

Euglena gracilis Z and a "sugar loving" variant strain E. gracilis var. saccharophila were investigated as producers of paramylon, a β-1,3-glucan polysaccharide with potential medicinal and industrial applications. The strains were grown under diurnal or dark growth conditions on a glucose-yeast extract medium supporting high-level paramylon production. Both strains produced the highest paramylon yields (7.4-8 g · L^-1 , respectively) while grown in the dark, but the maximum yield was achieved faster by E. gracilis var. saccharophila (48 h vs. 72 h). The glucose-to-paramylon yield coefficient Y_par/glu  = 0.46 ± 0.03 in the E. gracilis var. saccharophila cultivation, obtained in this study, is the highest reported to date. Proteomic analysis of the metabolic pathways provided molecular clues for the strain behavior observed during cultivation. For example, overexpression of enzymes in the gluconeogenesis/glycolysis pathways including fructokinase-1 and chloroplastic fructose-1,6-bisphosphatase (FBP) may have contributed to the faster rate of paramylon accumulation in E. gracilis var. saccharophila. Differentially expressed proteins in the early steps of chloroplastogenesis pathway including plastid uroporphyrinogen decarboxylases, photoreceptors, and a highly abundant (68-fold increase) plastid transketolase may have provided the E. gracilis var. saccharophila strain an advantage in paramylon production during diurnal cultivations. In conclusion, the variant strain E. gracilis var. saccharophila seems to be well suited for producing large amounts of paramylon. This work has also resulted in the identification of molecular targets for future improvement of paramylon production in E. gracilis, including the FBP and phosophofructokinase 1, the latter being a key regulator of glycolysis.

MiRNA-seq-based profiles of miRNAs in mulberry phloem sap provide insight into the pathogenic mechanisms of mulberry yellow dwarf disease

A wide range of miRNAs have been identified as phloem-mobile molecules that play important roles in coordinating plant development and physiology. Phytoplasmas are associated with hundreds of plant diseases, and the pathogenesis involved in the interactions between phytoplasmas and plants is still poorly understood. To analyse the molecular mechanisms of phytoplasma pathogenicity, the miRNAs profiles in mulberry phloem saps were examined in response to phytoplasma infection. A total of 86 conserved miRNAs and 19 novel miRNAs were identified, and 30 conserved miRNAs and 13 novel miRNAs were differentially expressed upon infection with phytoplasmas. The target genes of the differentially expressed miRNAs are involved in diverse signalling pathways showing the complex interactions between mulberry and phytoplasma. Interestingly, we found that mul-miR482a-5p was up-regulated in the infected phloem saps, and grafting experiments showed that it can be transported from scions to rootstock. Based on the results, the complexity and roles of the miRNAs in phloem sap and the potential molecular mechanisms of their changes were discussed. It is likely that the phytoplasma-responsive miRNAs in the phloem sap modulate multiple pathways and work cooperatively in response to phytoplasma infection, and their expression changes may be responsible for some symptoms in the infected plants.

Comparative Proteomic Analysis of Paulownia fortunei Response to Phytoplasma Infection with Dimethyl Sulfate Treatment

Paulownia fortunei is a widely cultivated economic forest tree species that is susceptible to infection with phytoplasma, resulting in Paulownia witches' broom (PaWB) disease. Diseased P. fortunei is characterized by stunted growth, witches' broom, shortened internodes, and etiolated and smaller leaves. To understand the molecular mechanism of its pathogenesis, we applied isobaric tags for relative and absolute quantitation (iTRAQ) and liquid chromatography coupled with tandem mass spectrometry approaches to study changes in the proteomes of healthy P. fortunei, PaWB-infected P. fortunei, and PaWB-infected P. fortunei treated with 15 mg·L⁻¹ or 75 mg·L⁻¹ dimethyl sulfate. We identified 2969 proteins and 104 and 32 differentially abundant proteins that were phytoplasma infection responsive and dimethyl sulfate responsive, respectively. Based on our analysis of the different proteomes, 27 PaWB-related proteins were identified. The protein-protein interactions of these 27 proteins were analyzed and classified into four groups (photosynthesis-related, energy-related, ribosome-related, and individual proteins). These PaWB-related proteins may help in developing a deeper understanding of how PaWB affects the morphological characteristics of P. fortunei and further establish the mechanisms involved in the response of P. fortunei to phytoplasma.

The use of high-throughput small RNA sequencing reveals differentially expressed microRNAs in response to aster yellows phytoplasma-infection in Vitis vinifera cv. 'Chardonnay'

Phytoplasmas are cell wall-less plant pathogenic bacteria responsible for major crop losses throughout the world. In grapevine they cause grapevine yellows, a detrimental disease associated with a variety of symptoms. The high economic impact of this disease has sparked considerable interest among researchers to understand molecular mechanisms related to pathogenesis. Increasing evidence exist that a class of small non-coding endogenous RNAs, known as microRNAs (miRNAs), play an important role in post-transcriptional gene regulation during plant development and responses to biotic and abiotic stresses. Thus, we aimed to dissect complex high-throughput small RNA sequencing data for the genome-wide identification of known and novel differentially expressed miRNAs, using read libraries constructed from healthy and phytoplasma-infected Chardonnay leaf material. Furthermore, we utilised computational resources to predict putative miRNA targets to explore the involvement of possible pathogen response pathways. We identified multiple known miRNA sequence variants (isomiRs), likely generated through post-transcriptional modifications. Sequences of 13 known, canonical miRNAs were shown to be differentially expressed. A total of 175 novel miRNA precursor sequences, each derived from a unique genomic location, were predicted, of which 23 were differentially expressed. A homology search revealed that some of these novel miRNAs shared high sequence similarity with conserved miRNAs from other plant species, as well as known grapevine miRNAs. The relative expression of randomly selected known and novel miRNAs was determined with real-time RT-qPCR analysis, thereby validating the trend of expression seen in the normalised small RNA sequencing read count data. Among the putative miRNA targets, we identified genes involved in plant morphology, hormone signalling, nutrient homeostasis, as well as plant stress. Our results may assist in understanding the role that miRNA pathways play during plant pathogenesis, and may be crucial in understanding disease symptom development in aster yellows phytoplasma-infected grapevines.

Quantitative proteome-level analysis of paulownia witches' broom disease with methyl methane sulfonate assistance reveals diverse metabolic changes during the infection and recovery processes

Paulownia witches' broom (PaWB) disease caused by phytoplasma is a fatal disease that leads to considerable economic losses. Although there are a few reports describing studies of PaWB pathogenesis, the molecular mechanisms underlying phytoplasma pathogenicity in Paulownia trees remain uncharacterized. In this study, after building a transcriptome database containing 67,177 sequences, we used isobaric tags for relative and absolute quantification (iTRAQ) to quantify and analyze the proteome-level changes among healthy P. fortunei (PF), PaWB-infected P. fortunei (PFI), and PaWB-infected P. fortunei treated with 20 mg L^-1 or 60 mg L^-1 methyl methane sulfonate (MMS) (PFI-20 and PFI-60, respectively). A total of 2,358 proteins were identified. We investigated the proteins profiles in PF vs. PFI (infected process) and PFI-20 vs. PFI-60 (recovered process), and further found that many of the MMS-response proteins mapped to "photosynthesis" and "ribosome" pathways. Based on our comparison scheme, 36 PaWB-related proteins were revealed. Among them, 32 proteins were classified into three functional groups: (1) carbohydrate and energy metabolism, (2) protein synthesis and degradation, and (3) stress resistance. We then investigated the PaWB-related proteins involved in the infected and recovered processes, and discovered that carbohydrate and energy metabolism was inhibited, and protein synthesis and degradation decreased, as the plant responded to PaWB. Our observations may be useful for characterizing the proteome-level changes that occur at different stages of PaWB disease. The data generated in this study may serve as a valuable resource for elucidating the pathogenesis of PaWB disease during phytoplasma infection and recovery stages.

Proteome Profiling of Paulownia Seedlings Infected with Phytoplasma

Phytoplasma is an insect-transmitted pathogen that causes witches' broom disease in many plants. Paulownia witches' broom is one of the most destructive diseases threatening Paulownia production. The molecular mechanisms associated with this disease have been investigated by transcriptome sequencing, but changes in protein abundance have not been investigated with isobaric tags for relative and absolute quantitation. Previous results have shown that methyl methane sulfonate (MMS) can help Paulownia seedlings recover from the symptoms of witches' broom and reinstate a healthy morphology. In this study, a transcriptomic-assisted proteomic technique was used to analyze the protein changes in phytoplasma-infected Paulownia tomentosa seedlings, phytoplasma-infected seedlings treated with 20 and 60 mg·L^-1 MMS, and healthy seedlings. A total of 2,051 proteins were obtained, 879 of which were found to be differentially abundant in pairwise comparisons between the sample groups. Among the differentially abundant proteins, 43 were related to Paulownia witches' broom disease and many of them were annotated to be involved in photosynthesis, expression of dwarf symptom, energy production, and cell signal pathways.

Secretion of Proteases by an Opportunistic Fungal Pathogen Scedosporium aurantiacum

Scedosporium aurantiacum is an opportunistic filamentous fungus increasingly isolated from the sputum of cystic fibrosis patients, and is especially prevalent in Australia. At the moment, very little is known about the infection mechanism of this fungus. Secreted proteases have been shown to contribute to fungal virulence in several studies with other fungi. Here we have compared the profiles of proteases secreted by a clinical isolate Scedosporium aurantiacum (WM 06.482) and an environmental strain (WM 10.136) grown on a synthetic cystic fibrosis sputum medium supplemented with casein or mucin. Protease activity was assessed using class-specific substrates and inhibitors. Subtilisin-like and trypsin-like serine protease activity was detected in all cultures. The greatest difference in the secretion of proteases between the two strains occurred in mucin-supplemented medium, where the activities of the elastase-like, trypsin-like and aspartic proteases were, overall, 2.5–75 fold higher in the clinical strain compared to the environmental strain. Proteases secreted by the two strains in the mucin-supplemented medium were further analyzed by mass spectrometry. Six homologs of fungal proteases were identified from the clinical strain and five from the environmental strain. Of these, three were common for both strains including a subtilisin peptidase, a putative leucine aminopeptidase and a PA-SaNapH-like protease. Trypsin-like protease was identified by mass spectrometry only in the clinical isolate even though trypsin-like activity was present in all cultures. In contrast, high elastase-like activity was measured in the culture supernatant of the clinical strain but could not be identified by mass spectrometry searching against other fungi in the NCBI database. Future availability of an annotated genome will help finalise identification of the S. aurantiacum proteases.

Proteomics analysis reveals novel host molecular mechanisms associated with thermotherapy of 'Ca. Liberibacter asiaticus'-infected citrus plants

BACKGROUND

Citrus Huanglongbing (HLB), which is linked to the bacterial pathogen 'Ca. Liberibacter asiaticus' (Las), is the most devastating disease of citrus plants, and longer-term control measures via breeding or genetic engineering have been unwieldy because all cultivated citrus species are susceptible to the disease. However, the degree of susceptibility varies among citrus species, which has prompted efforts to identify potential Las resistance/tolerance-related genes in citrus plants for application in breeding or genetic engineering programs. Plant exposure to one form of stress has been shown to serendipitously induce innate resistance to other forms of stress and a recent study showed that continuous heat treatment (40 to 42 °C) reduced Las titer and HLB-associated symptoms in citrus seedlings. The goal of the present study was to apply comparative proteomics analysis via 2-DE and mass spectrometry to elucidate the molecular processes associated with heat-induced mitigation of HLB in citrus plants. Healthy or Las-infected citrus grapefruit plants were exposed to room temperature or to continuous heat treatment of 40 °C for 6 days.

RESULTS

An exhaustive total protein extraction process facilitated the identification of 107 differentially-expressed proteins in response to Las and/or heat treatment, which included a strong up-regulation of chaperones including small (23.6, 18.5 and 17.9 kDa) heat shock proteins, a HSP70-like protein and a ribulose-1,5-bisphosphate carboxylase oxygenase (RuBisCO)-binding 60 kDa chaperonin, particularly in response to heat treatment. Other proteins that were generally down-regulated due to Las infection but up-regulated in response to heat treatment include RuBisCO activase, chlorophyll a/b binding protein, glucosidase II beta subunit-like protein, a putative lipoxygenase protein, a ferritin-like protein, and a glutathione S-transferase.

CONCLUSIONS

The differentially-expressed proteins identified in this study highlights a premier characterization of the molecular mechanisms potentially involved in the reversal of Las-induced pathogenicity processes in citrus plants and are hence proposed targets for application towards the development of cisgenic Las-resistant/tolerant citrus plants.

Genome-Wide Identification and Analysis of MicroRNAs Involved in Witches'-Broom Phytoplasma Response in Ziziphus jujuba

MicroRNAs (miRNAs) play an important role in responding to biotic and abiotic stresses in plants. Jujube witches'-broom a phytoplasma disease of Ziziphus jujuba is prevalent in China and is a serious problem to the industry. However, the molecular mechanism of the disease is poorly understood. In this study, genome-wide identification and analysis of microRNAs in response to witches'-broom was performed. A total of 85 conserved miRNA unique sequences belonging to 32 miRNA families and 24 novel miRNA unique sequences, including their complementary miRNA* strands were identified from small RNA libraries derived from a uninfected and witches'-broom infected Z. jujuba plant. Differentially expressed miRNAs associated with Jujube witches'-broom disease were investigated between the two libraries, and 12 up-regulated miRNAs and 10 down- regulated miRNAs identified with more than 2 fold changes. Additionally, 40 target genes of 85 conserved miRNAs and 49 target genes of 24 novel miRNAs were predicted and their putative functions assigned. Using the modified 5'-RACE method, we confirmed that SPL and MYB were cleaved by miR156 and miR159, respectively. Our results provide insight into the molecular mechanisms of witches'-broom disease in Z. jujuba.

Shotgun Label-free Proteomic Analysis of Clubroot (Plasmodiophora brassicae) Resistance Conferred by the Gene Rcr1 in Brassica rapa

Clubroot, caused by the plasmodiophorid pathogen Plasmodiophora brassicae, is one of the most serious diseases on Brassica crops worldwide and a major threat to canola production in western Canada. Host resistance is the key strategy for clubroot management on canola. Several clubroot resistance (CR) genes have been identified, but the mechanisms associated with these CR genes are poorly understood. In the current study, a label-free shotgun proteomic approach was used to profile and compare the proteomes of Brassica rapa carrying and not carrying the CR gene Rcr1 in response to P. brassicae infection. A total of 527 differentially accumulated proteins (DAPs) were identified between the resistant (with Rcr1) and susceptible (without Rcr1) samples, and functional annotation of these DAPs indicates that the perception of P. brassicae and activation of defense responses are triggered via an unique signaling pathway distinct from common modes of recognition receptors reported with many other plant-pathogen interactions; this pathway appears to act in a calcium-independent manner through a not-well-defined cascade of mitogen-activated protein kinases and may require the ubiquitin-26S proteasome found to be related to abiotic stresses, especially the cold-stress tolerance in other studies. Both up-regulation of defense-related and down-regulation of pathogenicity-related metabolism was observed in plants carrying Rcr1, and these functions may all contribute to the CR mediated by Rcr1. These results, combined with those of transcriptomic analysis reported earlier, improved our understanding of molecular mechanisms associated with Rcr1 and CR at large, and identified candidate metabolites or pathways related to specific resistance mechanisms. Deploying CR genes with different modes of action may help improve the durability of CR.

Comparative shotgun proteomic analysis of wild and domesticated Opuntia spp. species shows a metabolic adaptation through domestication

The Opuntia genus is widely distributed in America, but the highest richness of wild species are found in Mexico, as well as the most domesticated Opuntia ficus-indica, which is the most domesticated species and an important crop in agricultural economies of arid and semiarid areas worldwide. During domestication process, the Opuntia morphological characteristics were favored, such as less and smaller spines in cladodes and less seeds in fruits, but changes at molecular level are almost unknown. To obtain more insights about the Opuntia molecular changes through domestication, a shotgun proteomic analysis and database-dependent searches by homology was carried out. >1000 protein species were identified and by using a label-free quantitation method, the Opuntia proteomes were compared in order to identify differentially accumulated proteins among wild and domesticated species. Most of the changes were observed in glucose, secondary, and 1C metabolism, which correlate with the observed protein, fiber and phenolic compounds accumulation in Opuntia cladodes. Regulatory proteins, ribosomal proteins, and proteins related with response to stress were also observed in differential accumulation. These results provide new valuable data that will help to the understanding of the molecular changes of Opuntia species through domestication.

BIOLOGICAL SIGNIFICANCE

Opuntia species are well adapted to dry and warm conditions in arid and semiarid regions worldwide, and they are highly productive plants showing considerable promises as an alternative food source. However, there is a gap regarding Opuntia molecular mechanisms that enable them to grow in extreme environmental conditions and how the domestication processes has changed them. In the present study, a shotgun analysis was carried out to characterize the proteomes of five Opuntia species selected by its domestication degree. Our results will help to a better understanding of proteomic features underlying the selection and specialization under evolution and domestication of Opuntia and will provide a platform for basic biology research and gene discovery.

Omics studies of citrus, grape and rosaceae fruit trees

Recent advance of bioinformatics and analytical apparatuses such as next generation DNA sequencer (NGS) and mass spectrometer (MS) has brought a big wave of comprehensive study to biology. Comprehensive study targeting all genes, transcripts (RNAs), proteins, metabolites, hormones, ions or phenotypes is called genomics, transcriptomics, proteomics, metabolomics, hormonomics, ionomics or phenomics, respectively. These omics are powerful approaches to identify key genes for important traits, to clarify events of physiological mechanisms and to reveal unknown metabolic pathways in crops. Recently, the use of omics approach has increased dramatically in fruit tree research. Although the most reported omics studies on fruit trees are transcriptomics, proteomics and metabolomics, and a few is reported on hormonomics and ionomics. In this article, we reviewed recent omics studies of major fruit trees, i.e. citrus, grapevine and rosaceae fruit trees. The effectiveness and prospects of omics in fruit tree research will as well be highlighted.

Citrus tristeza virus (CTV) Causing Proteomic and Enzymatic Changes in Sweet Orange Variety "Westin"

Citrus Tristeza disease, caused by CTV (Citrus tristeza virus), committs citrus plantations around the world and specifically attacks phloem tissues of the plant. The virus exists as a mixture of more or less severe variants, which may or may not cause symptoms of Tristeza. The objective of this study was to analyze the changes caused by CTV in the proteome of stems of sweet orange, as well as in the activity and gene expression of antioxidant enzymes. The CTV-infected sweet orange displayed mild symptoms, which were characterized by the presence of sparse stem pitting throughout their stems. The presence of virus was confirmed by RT-PCR. Proteomic analysis by 2DE-PAGE-MS / MS revealed the identity of 40 proteins differentially expressed between CTV- infected and -non-infected samples. Of these, 33 were up-regulated and 7 were down-regulated in CTV-infected samples. Among the proteins identified stands out a specific from the virus, the coat protein. Other proteins identified are involved with oxidative stress and for this their enzymatic activity was measured. The activity of superoxide dismutase (SOD) was higher in CTV-infected samples, as catalase (CAT) showed higher activity in uninfected samples. The activity of guaiacol peroxidase (GPX) did not vary significantly between samples. However, ascorbate peroxidase (APX) was more active in the infected samples. The relative expression of the genes encoding CAT, SOD, APX and GPX was analyzed by quantitative real time PCR (RT-qPCR). The CTV-infected samples showed greater accumulation of transcripts, except for the CAT gene. This gene showed higher expression in the uninfected samples. Taken together, it can be concluded that the CTV affects the protein profile and activity and gene expression of antioxidant enzymes in plants infected by this virus.

In-Depth Transcriptome Sequencing of Mexican Lime Trees Infected with Candidatus Phytoplasma aurantifolia

Witches' broom disease of acid lime greatly affects the production of Mexican lime in Iran. It is caused by a phytoplasma (Candidatus Phytoplasma aurantifolia). However, the molecular mechanisms that underlie phytoplasma pathogenicity and the mode of interactions with host plants are largely unknown. Here, high-throughput transcriptome sequencing was conducted to explore gene expression signatures associated with phytoplasma infection in Mexican lime trees. We assembled 78,185 unique transcript sequences (unigenes) with an average length of 530 nt. Of these, 41,805 (53.4%) were annotated against the NCBI non-redundant (nr) protein database using a BLASTx search (e-value ≤ 1e-5). When the abundances of unigenes in healthy and infected plants were compared, 2,805 transcripts showed significant differences (false discovery rate ≤ 0.001 and log2 ratio ≥ 1.5). These differentially expressed genes (DEGs) were significantly enriched in 43 KEGG metabolic and regulatory pathways. The up-regulated DEGs were mainly categorized into pathways with possible implication in plant-pathogen interaction, including cell wall biogenesis and degradation, sucrose metabolism, secondary metabolism, hormone biosynthesis and signalling, amino acid and lipid metabolism, while down-regulated DEGs were predominantly enriched in ubiquitin proteolysis and oxidative phosphorylation pathways. Our analysis provides novel insight into the molecular pathways that are deregulated during the host-pathogen interaction in Mexican lime trees infected by phytoplasma. The findings can be valuable for unravelling the molecular mechanisms of plant-phytoplasma interactions and can pave the way for engineering lime trees with resistance to witches' broom disease.

RNA-Seq profile of flavescence dorée phytoplasma in grapevine

BACKGROUND

The phytoplasma-borne disease flavescence dorée is still a threat to European viticulture, despite mandatory control measures and prophylaxis against the leafhopper vector. Given the economic importance of grapevine, it is essential to find alternative strategies to contain the spread, in order to possibly reduce the current use of harmful insecticides. Further studies of the pathogen, the vector and the mechanisms of phytoplasma-host interactions could improve our understanding of the disease. In this work, RNA-Seq technology followed by three de novo assembly strategies was used to provide the first comprehensive transcriptomics landscape of flavescence dorée phytoplasma (FD) infecting field-grown Vitis vinifera leaves.

RESULTS

With an average of 8300 FD-mapped reads per library, we assembled 347 sequences, corresponding to 215 annotated genes, and identified 10 previously unannotated genes, 15 polycistronic transcripts and three genes supposedly localized in the gaps of the FD92 draft genome. Furthermore, we improved the annotation of 44 genes with the addition of 5'/3' untranslated regions. Functional classification revealed that the most expressed genes were either related to translation and protein biosynthesis or hypothetical proteins with unknown function. Some of these hypothetical proteins were predicted to be secreted, so they could be bacterial effectors with a potential role in modulating the interaction with the host plant. Interestingly, qRT-PCR validation of the RNA-Seq expression values confirmed that a group II intron represented the FD genomic region with the highest expression during grapevine infection. This mobile element may contribute to the genomic plasticity that is necessary for the phytoplasma to increase its fitness and endorse host-adaptive strategies.

CONCLUSIONS

The RNA-Seq technology was successfully applied for the first time to analyse the FD global transcriptome profile during grapevine infection. Our results provided new insights into the transcriptional organization and gene structure of FD. This may represent the starting point for the application of high-throughput sequencing technologies to study differential expression in FD and in other phytoplasmas with an unprecedented resolution.

Complete chloroplast genome sequence of Omani lime (Citrus aurantiifolia) and comparative analysis within the rosids

The genus Citrus contains many economically important fruits that are grown worldwide for their high nutritional and medicinal value. Due to frequent hybridizations among species and cultivars, the exact number of natural species and the taxonomic relationships within this genus are unclear. To compare the differences between the Citrus chloroplast genomes and to develop useful genetic markers, we used a reference-assisted approach to assemble the complete chloroplast genome of Omani lime (C. aurantiifolia). The complete C. aurantiifolia chloroplast genome is 159,893 bp in length; the organization and gene content are similar to most of the rosids lineages characterized to date. Through comparison with the sweet orange (C. sinensis) chloroplast genome, we identified three intergenic regions and 94 simple sequence repeats (SSRs) that are potentially informative markers with resolution for interspecific relationships. These markers can be utilized to better understand the origin of cultivated Citrus. A comparison among 72 species belonging to 10 families of representative rosids lineages also provides new insights into their chloroplast genome evolution.

Transcriptomics assisted proteomic analysis of Nicotiana occidentalis infected by Candidatus Phytoplasma mali strain AT

Phytoplasmas are pathogenic bacteria within the class of Mollicutes, which are associated with more than 1000 plant diseases. In this study, we applied quantitative mass spectrometry to analyse affected pathways of the model plant tobacco (Nicotiana occidentalis) upon Candidatus Phytoplasma mali strain AT infection. Using tissue obtained from leaf midribs, 1466 plant-assigned proteins were identified. For 1019 of these proteins, we could reproducibly quantify the expression changes of infected versus noninfected plants, of which 157 proteins were up- and 173 proteins were downregulated. Differential expression took place in a number of pathways, among others strong downregulation of porphyrin and chlorophyll metabolism and upregulation of alpha-linolenic acid metabolism, which was consistent with observed increased levels of jasmonic acid, a key signal molecule of plant defence. Our data shed light on the molecular networks that are involved in defence of plants against phytoplasma infection and provide a resource for further studies.

Phytoplasma-Responsive microRNAs Modulate Hormonal, Nutritional, and Stress Signalling Pathways in Mexican Lime Trees

BACKGROUND

Witches' broom disease of Mexican lime (Citrus aurantifolia L.), which is associated to the phytoplasma 'Candidatus Phytoplasma aurantifolia', is a devastating disease that results in significant economic losses. Plants adapt to biotic stresses by regulating gene expression at the transcriptional and post-transcriptional levels. MicroRNAs (miRNAs) are a recently identified family of molecules that regulate plant responses to environmental stresses through post-transcriptional gene silencing.

METHODS

Using a high-throughput approach to sequence small RNAs, we compared the expression profiles of miRNAs in healthy Mexican lime trees and in plants infected with 'Ca. P. aurantifolia'.

RESULTS

Our results demonstrated the involvement of different miRNAs in the response of Mexican lime trees to infection by 'Ca. P. aurantifolia'. We identified miRNA families that are expressed differentially upon infection with phytoplasmas. Most of the miRNAs had variants with small sequence variations (isomiRs), which are expressed differentially in response to pathogen infection.

CONCLUSIONS

It is likely that the miRNAs that are expressed differentially in healthy and phytoplasma-infected Mexican lime trees are involved in coordinating the regulation of hormonal, nutritional, and stress signalling pathways, and the complex interactions between them. Future research to elucidate the roles of these miRNAs should improve our understanding of the level of diversity of specific plant responses to phytoplasmas.