Discovery of genes related to witches broom disease in Paulownia tomentosa × Paulownia fortunei by a De Novo assembled transcriptome

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.

Combined transcriptome and metabolome analysis of Nerium indicum L. elaborates the key pathways that are activated in response to witches' broom disease

BACKGROUND

Nerium indicum Mill. is an ornamental plant that is found in parks, riversides, lakesides, and scenic areas in China and other parts of the world. Our recent survey indicated the prevalence of witches' broom disease (WBD) in Guangdong, China. To find out the possible defense strategies against WBD, we performed a MiSeq based ITS sequencing to identify the possible casual organism, then did a de novo transcriptome sequencing and metabolome profiling in the phloem and stem tip of N. indicum plants suffering from WBD compared to healthy ones.

RESULTS

The survey showed that Wengyuen county and Zengcheng district had the highest disease incidence rates. The most prevalent microbial species in the diseased tissues was Cophinforma mamane. The transcriptome sequencing resulted in the identification of 191,224 unigenes of which 142,396 could be annotated. There were 19,031 and 13,284 differentially expressed genes (DEGs) between diseased phloem (NOWP) and healthy phloem (NOHP), and diseased stem (NOWS) and healthy stem (NOHS), respectively. The DEGs were enriched in MAPK-signaling (plant), plant-pathogen interaction, plant-hormone signal transduction, phenylpropanoid and flavonoid biosynthesis, linoleic acid and α-linoleic acid metabolism pathways. Particularly, we found that N. indicum plants activated the phytohormone signaling, MAPK-signaling cascade, defense related proteins, and the biosynthesis of phenylpropanoids and flavonoids as defense responses to the pathogenic infection. The metabolome profiling identified 586 metabolites of which 386 and 324 metabolites were differentially accumulated in NOHP vs NOWP and NOHS and NOWS, respectively. The differential accumulation of metabolites related to phytohormone signaling, linoleic acid metabolism, phenylpropanoid and flavonoid biosynthesis, nicotinate and nicotinamide metabolism, and citrate cycle was observed, indicating the role of these pathways in defense responses against the pathogenic infection.

CONCLUSION

Our results showed that Guangdong province has a high incidence of WBD in most of the surveyed areas. C. mamane is suspected to be the causing pathogen of WBD in N. indicum. N. indicum initiated the MAPK-signaling cascade and phytohormone signaling, leading to the activation of pathogen-associated molecular patterns and hypersensitive response. Furthermore, N. indicum accumulated high concentrations of phenolic acids, coumarins and lignans, and flavonoids under WBD. These results provide scientific tools for the formulation of control strategies of WBD in N. indicum.

The potential of Paulownia fortunei seedlings for the phytoremediation of manganese slag amended with spent mushroom compost

The use of phytoremediation was an efficient strategy for the restoration of mine slag and the addition of modifier was favorable for improving the phytoremediation efficiency. Herein, spent mushroom compost (SMC) was added in manganese (Mn) slag to reveal the phytoremediation potential of Paulownia fortunei seedlings. The transportation, subcellular distribution and chemical forms of Mn in P. fortunei, the diurnal variation of photosynthesis and antioxidant enzyme activities in P. fortunei leaves were measured to reveal the effect of SMC (mass ratios of 10%, M+) on the phytoremediation of Mn slag. Results showed that the addition of SMC increased the accumulation content of Mn by 408.54% due to the increased biomass of P. fortunei seedlings. After SMC amendment, the maximum net photosynthetic rate (Pn) increased and the superoxide dismutase (SOD) activities decreased significantly (p < 0.05), which was beneficial to the tolerance of leaves to Mn stress. SMC amendment maintained the cell structural integrity of P. fortunei seedlings observed by transmission electron microscope (TEM). Additionally, SMC amendment decreased the damage level of Mn to the cell of P. fortunei seedlings by using function groups (-CH₃ and -COOH) to bond Mn in the cell walls and vacuoles. SMC amendment reduced the Mn toxicity to P. fortunei seedlings and improved the phytoremediation capacity.

'Candidatus Phytoplasma solani' interferes with the distribution and uptake of iron in tomato

Background

‘Candidatus Phytoplasma solani’ is endemic in Europe and infects a wide range of weeds and cultivated plants. Phytoplasmas are prokaryotic plant pathogens that colonize the sieve elements of their host plant, causing severe alterations in phloem function and impairment of assimilate translocation. Typical symptoms of infected plants include yellowing of leaves or shoots, leaf curling, and general stunting, but the molecular mechanisms underlying most of the reported changes remain largely enigmatic. To infer a possible involvement of Fe in the host-phytoplasma interaction, we investigated the effects of ‘Candidatus Phytoplasma solani’ infection on tomato plants (Solanum lycopersicum cv. Micro-Tom) grown under different Fe regimes.

Results

Both phytoplasma infection and Fe starvation led to the development of chlorotic leaves and altered thylakoid organization. In infected plants, Fe accumulated in phloem tissue, altering the local distribution of Fe. In infected plants, Fe starvation had additive effects on chlorophyll content and leaf chlorosis, suggesting that the two conditions affected the phenotypic readout via separate routes. To gain insights into the transcriptional response to phytoplasma infection, or Fe deficiency, transcriptome profiling was performed on midrib-enriched leaves. RNA-seq analysis revealed that both stress conditions altered the expression of a large (> 800) subset of common genes involved in photosynthetic light reactions, porphyrin / chlorophyll metabolism, and in flowering control. In Fe-deficient plants, phytoplasma infection perturbed the Fe deficiency response in roots, possibly by interference with the synthesis or transport of a promotive signal transmitted from the leaves to the roots.

Conclusions

‘Candidatus Phytoplasma solani’ infection changes the Fe distribution in tomato leaves, affects the photosynthetic machinery and perturbs the orchestration of root-mediated transport processes by compromising shoot-to-root communication.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6062-x) contains supplementary material, which is available to authorized users.

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.

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.

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.

Transcriptome profiling analysis revealed co-regulation of multiple pathways in jujube during infection by 'Candidatus Phytoplasma ziziphi'

BACKGROUND

Jujube witches' broom (JWB), caused by a phytoplasma, devastates jujube tree (Ziziphus jujuba) growth and production in Asia. Although host responses to phytoplasmas are studied at the phenotypic, physiological, biochemical and molecular levels, it remains unclear how a host plant responds at the molecular level during the primary stage of infection.

METHODS

To understand the response of the jujube tree to JWB infection, leaves were sampled at different times during the phytoplasma infection. Transcriptomic analyses at six stages were performed to reveal how phytoplasma infection affects Chinese jujube gene expression through the determination of the key differentially expressed genes (DEGs), and their related pathways. Quantitative real-time PCR was applied to validate 10 differentially expressed genes at different JWB phytoplasma infection stages.

RESULTS

A total of 25,067 unigenes were mapped to jujube genome reference sequences. In the first infection stage (0-2 weeks after grafting (WAG), a total of 582 jujube genes were differentially regulated but no visible symptoms appeared. Quite a few DEGs related to abscisic acid (ABA) and cytokinin (CTK) were down-regulated, while some related to jasmonic acid (JA) and salicylic acid (SA) were up-regulated, Genes related to plant-pathogen interaction were also differentially expressed. In the second infection stage (37-39WAG), witches' broom symptoms were visible, and a total of 4373 DEGs were identified. Genes involved in biosynthesis and signal transduction of ABA, brassinosteroid (BR), CTK, ethylene (ET), and auxin (IAA), GA, JA and SA, plant-pathogen interaction, flavonoid biosynthesis genes were significantly regulated, suggesting that jujube trees activated defense factors related to SA, JA, ABA and secondary metabolites to defend against phytoplasma infection. By the third infection stage (48-52WAG), serious symptoms occurred and 3386 DEGs were identified. Most DEGs involved in biosynthesis and signal transduction of JA, SA and GA were up-regulated, while those relating to ABA were down-regulated. Genes involved in plant-pathogen interaction were up- or down-regulated, while phenylpropanoid and flavonoid biosynthesis genes were significantly up-regulated. Meanwhile, DEGs involved in photosynthesis, chlorophyll and peroxisome biosynthesis, and carbohydrate metabolism were down-regulated. These results suggested that phytoplasma infection had completely destroyed jujube trees' defense system and had disturbed chlorophyll synthesis and photosynthetic activity in the infected leaves at the late stage, resulting in yellow leaves and other JWB symptoms.

DISCUSSION

The results in this report suggested that phytohormone biosynthesis and signal transduction, photosynthesis, and secondary metabolism all played important roles in the battle between colonization and defense in the interaction between Ca. Phytoplasma ziziphi and jujube.

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.

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.

Quantitative Proteomic and Transcriptomic Study on Autotetraploid Paulownia and Its Diploid Parent Reveal Key Metabolic Processes Associated with Paulownia Autotetraploidization

Polyploidy plays a very important role in speciation and plant evolution by way of genomic merging and doubling. In the process of polyploidy, rapid genomic, and transcriptomic changes have been observed and researched. However, proteomic divergence caused by the effects of polyploidization is still poorly understood. In the present study, we used iTRAQ coupled with mass spectrometry to quantitatively analyze proteomic changes in the leaves of autotetraploid Paulownia and its diploid parent. A total of 2963 proteins were identified and quantified. Among them, 463 differentially abundant proteins were detected between autotetraploid Paulownia and its diploid parent, and 198 proteins were found to be non-additively abundant in autotetraploid Paulownia, suggesting the presence of non-additive protein regulation during genomic merger and doubling. We also detected 1808 protein-encoding genes in previously published RNA sequencing data. We found that 59 of the genes that showed remarkable changes at mRNA level encoded proteins with consistant changes in their abundance levels, while a further 48 genes that showed noteworthy changes in their expression levels encoded proteins with opposite changes in their abundance levels. Proteins involved in posttranslational modification, protein turnover, and response to stimulus, were significantly enriched among the non-additive proteins, which may provide some of the driving power for variation and adaptation in autopolyploids. Quantitative real-time PCR analysis verified the expression patterns of related protein-coding genes. In addition, we found that the percentage of differentially abundant proteins that matched previously reported differentially expressed genes was relatively low.

Identification of microRNAs and their targets in Paulownia fortunei plants free from phytoplasma pathogen after methyl methane sulfonate treatment

MicroRNAs (miRNAs) play major roles in plant responses to various biotic and abiotic stresses by regulating gene expression at the transcriptional and post-transcriptional levels. Paulownia witches' broom (PaWB) disease caused by phytoplasmas reduces Paulownia production worldwide. In this study, we investigated the miRNA-mediated plant response to PaWB phytoplasma by Illumina sequencing and degradome analysis of Paulownia fortunei small RNAs (sRNAs). The sRNA and degradome libraries were constructed from healthy and diseased P. fortunei plants and the plants free from phytoplasma pathogen after 60 mg L(-1) methyl methane sulfonate treatment. A total of 96 P. fortunei-conserved miRNAs and 83 putative novel miRNAs were identified. Among them, 37 miRNAs (17 conserved, 20 novel) were found to be differentially expressed in response to PaWB phytoplasma infection. In addition, 114 target genes for 18 of the conserved miRNA families and 33 target genes for 15 of the novel miRNAs in P. fortunei were detected. The expression patterns of 14 of the PaWB phytoplasma-responsive miRNAs and 12 target genes were determined by quantitative real-time polymerase chain reaction (qPCR) experiments. A functional analysis of the miRNA targets indicated that these targeted genes may regulate transcription, stress response, nitrogen metabolism, and various other activities. Our results will help identify the potential roles of miRNAs involved in protecting P. fortunei from diseases.

Transcriptome, microRNA, and degradome analyses of the gene expression of Paulownia with phytoplamsa

BACKGROUND

Paulownia witches' broom (PaWB) is a fatal disease of Paulownia caused by a phytoplasma. In previous studies, we found that plants with PaWB symptoms would revert to a healthy morphology after methyl methane sulfonate (MMS) treatment. To completely understand the gene expression profiles of the Paulownia-phytoplasma interaction, three high-throughput sequencing technologies were used to investigate changes of gene expression and microRNAs (miRNAs) in healthy Paulownia tomentosa plantlets, PaWB-infected plantlets, and PaWB-infected plantlets treated with 60 mg · L(-1) MMS.

METHODS

Transcriptome, miRNAs and degradome sequencing were performed to explore the global gene expression profiles in the process of Paulownia tomentosa with phytoplasma infection.

RESULTS

A total of 98,714 all-unigenes, 62 conserved miRNAs, and 35 novel miRNAs were obtained, among which 902 differentially expressed genes (DEGs) and 24 miRNAs were found to be associated with PaWB disease. Subsequently, the target genes of these miRNAs were predicted by degradome sequencing. Interestingly, we found that 19 target genes of these differentially expressed miRNAs were among the 902 DEGs. The targets of pau-miR156g, pau-miR403, and pau-miR166c were significantly up-regulated in the P. tomentosa plantlets infected with phytoplasma. Interaction of miRNA -target genes mediated gene expression related to PaWB were identified.

CONCLUSIONS

The results elucidated the possible roles of the regulation of genes and miRNAs in the Paulownia-phytoplasma interaction, which will enrich our understanding of the mechanisms of PaWB disease in this plant.

Plant-Pathogen Interaction-Related MicroRNAs and Their Targets Provide Indicators of Phytoplasma Infection in Paulownia tomentosa × Paulownia fortunei

Paulownia witches’ broom (PaWB) caused by a phytoplasma, has caused extensive losses in the yields of paulownia timber and resulted in significant economic losses. However, the molecular mechanisms in Paulownia that underlie the phytoplasma stress are poorly characterized. In this study, we use an Illumina platform to sequence four small RNA libraries and four degradome sequencing libraries derived from healthy, PaWB-infected, and PaWB-infected 15 mg·L⁻¹ and 30 mg·L⁻¹ methyl methane sulfonate (MMS)-treated plants. In total, 125 conserved and 118 novel microRNAs (miRNAs) were identified and 33 miRNAs responsive to PaWB disease were discovered. Furthermore, 166 target genes for 18 PaWB disease-related miRNAs were obtained, and found to be involved in plant-pathogen interaction and plant hormone signal transduction metabolic pathways. Eleven miRNAs and target genes responsive to PaWB disease were examined by a quantitative real-time PCR approach. Our findings will contribute to studies on miRNAs and their targets in Paulownia, and provide new insights to further understand plant-phytoplasma interactions.

Discovery of microRNAs and transcript targets related to witches' broom disease in Paulownia fortunei by high-throughput sequencing and degradome approach

Paulownia witches' broom (PaWB) caused by the phytoplasma is a devastating disease of Paulownia trees. It has caused heavy yield losses to Paulownia production worldwide. However, knowledge of the transcriptional and post-transcriptional regulation of gene expression by microRNAs (miRNAs), especially miRNAs responsive to PaWB disease stress, is still rudimentary. In this study, to identify miRNAs and their transcript targets that are responsive to PaWB disease stress, six sequencing libraries were constructed from healthy (PF), PaWB-infected (PFI), and PaWB-infected, 20 mg L(-1) methyl methane sulfonate-treated (PFI20) P. fortunei seedlings. As a result, 95 conserved miRNAs belonging to 18 miRNA families, as well as 122 potential novel miRNAs, were identified. Most of them were found to be a response to PaWB disease-induced stress, and the expression levels of these miRNAs were validated by quantitative real-time PCR analysis. The study simultaneously identified 109 target genes from the P. fortunei for 14 conserved miRNA families and 24 novel miRNAs by degradome sequencing. Furthermore, the functions of the miRNA targets were annotated based on Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis. The results presented here provide the groundwork for further analysis of miRNAs and target genes responsive to the PaWB disease stress, and could be also useful for addressing new questions to better understand the mechanisms of plant infection by phytoplasma in the future.

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.

Transcriptomics-based analysis using RNA-Seq of the coconut (Cocos nucifera) leaf in response to yellow decline phytoplasma infection

Invasive phytoplasmas wreak havoc on coconut palms worldwide, leading to high loss of income, food insecurity and extreme poverty of farmers in producing countries. Phytoplasmas as strictly biotrophic insect-transmitted bacterial pathogens instigate distinct changes in developmental processes and defence responses of the infected plants and manipulate plants to their own advantage; however, little is known about the cellular and molecular mechanisms underlying host-phytoplasma interactions. Further, phytoplasma-mediated transcriptional alterations in coconut palm genes have not yet been identified. This study evaluated the whole transcriptome profiles of naturally infected leaves of Cocos nucifera ecotype Malayan Red Dwarf in response to yellow decline phytoplasma from group 16SrXIV, using RNA-Seq technique. Transcriptomics-based analysis reported here identified genes involved in coconut innate immunity. The number of down-regulated genes in response to phytoplasma infection exceeded the number of genes up-regulated. Of the 39,873 differentially expressed unigenes, 21,860 unigenes were suppressed and 18,013 were induced following infection. Comparative analysis revealed that genes associated with defence signalling against biotic stimuli were significantly overexpressed in phytoplasma-infected leaves versus healthy coconut leaves. Genes involving cell rescue and defence, cellular transport, oxidative stress, hormone stimulus and metabolism, photosynthesis reduction, transcription and biosynthesis of secondary metabolites were differentially represented. Our transcriptome analysis unveiled a core set of genes associated with defence of coconut in response to phytoplasma attack, although several novel defence response candidate genes with unknown function have also been identified. This study constitutes valuable sequence resource for uncovering the resistance genes and/or susceptibility genes which can be used as genetic tools in disease resistance breeding.

Identification of genes related to the phenotypic variations of a synthesized Paulownia (Paulownia tomentosa×Paulownia fortunei) autotetraploid

Paulownia is a fast-growing deciduous tree native to China. It has great economic importance for the pulp and paper industries, as well as ecological prominence in forest ecosystems. Paulownia is of much interest to plant breeder keen to explore new plant varieties by selecting on the basis of phenotype. A newly synthesized autotetraploid Paulownia exhibited advanced characteristics, such as greater yield, and higher resistance than the diploid tree. However, tissue-specific transcriptome and genomic data in public databases are not sufficient to understand the molecular mechanisms associated with genome duplication. To evaluate the effects of genome duplication on the phenotypic variations in Paulownia tomentosa×Paulownia fortunei, the transcriptomes of the autotetraploid and diploid Paulownia were compared. Using Illumina sequencing technology, a total of 82,934 All-unigenes with a mean length of 1109 bp were assembled. The data revealed numerous differences in gene expression between the two transcriptomes, including 718 up-regulated and 667 down-regulated differentially expressed genes between the two Paulownia trees. An analysis of the pathway and gene annotations revealed that genes involved in nucleotide sugar metabolism in plant cell walls were down-regulated, and genes involved in the light signal pathway and the biosynthesis of structural polymers were up-regulated in autotetraploid Paulownia. The differentially expressed genes may contribute to the observed phenotypic variations between diploid and autotetraploid Paulownia. These results provide a significant resource for understanding the variations in Paulownia polyploidization and will benefit future breeding work.

Plant-pathogen interaction, circadian rhythm, and hormone-related gene expression provide indicators of phytoplasma infection in Paulownia fortunei

Phytoplasmas are mycoplasma-like pathogens of witches’ broom disease, and are responsible for serious yield losses of Paulownia trees worldwide. The molecular mechanisms of disease development in Paulownia are of considerable interest, but still poorly understood. Here, we have applied transcriptome sequencing technology and a de novo assembly approach to analyze gene expression profiles in Paulownia fortunei infected by phytoplasmas. Our previous researches suggested that methyl methane sulfonated (MMS) could reverse the effects of the infection. In this study, leaf samples from healthy, infected, and both infected and methyl methane sulfonate treated plants were analyzed. The results showed that the gene expression profile of P. fortunei underwent dramatic changes after Paulownia witches’ broom (PaWB) phytoplasma infection. Genes that encoded key enzymes in plant-pathogen interaction processes were significantly up-regulated in the PaWB-infected Paulownia. Genes involved in circadian rhythm and hormone-related genes were also altered in Paulownia after PaWB infection. However, after the PaWB-infected plants were treated with MMS, the expression profiles of these genes returned to the levels in the healthy controls. The data will help identify potential PaWB disease-resistance genes that could be targeted to inhibit the growth and reproduction of the pathogen and to increase plant resistance.

Morphological changes of Paulownia seedlings infected phytoplasmas reveal the genes associated with witches' broom through AFLP and MSAP

Paulownia witches' broom (PaWB) caused by phytoplasma might result in devastating damage to the growth and wood production of Paulownia. To study the effect of phytoplasma on DNA sequence and to discover the genes related to PaWB occurrence, DNA polymorphisms and DNA methylation levels and patterns in PaWB seedlings, the ones treated with various concentration of methyl methane sulfonate (MMS) and healthy seedlings were investigated with amplified fragment length polymorphism (AFLP) and methylation-sensitive amplification polymorphism (MSAP). Our results indicated that PaWB seedlings recovered a normal morphology, similar to healthy seedlings, after treatment with more than 20 mg·L⁻¹ MMS; Phytoplasma infection did not change the Paulownia genomic DNA sequence at AFLP level, but changed the global DNA methylation levels and patterns; Genes related to PaWB were discovered through MSAP and validated using quantitative real-time PCR (qRT-PCR). These results implied that changes of DNA methylation levels and patterns were closely related to the morphological changes of seedlings infected with phytoplasmas.

Identification of genes related to Paulownia witches' broom by AFLP and MSAP

DNA methylation is believed to play important roles in regulating gene expression in plant growth and development. Paulownia witches' broom (PaWB) infection has been reported to be related to gene expression changes in paulownia plantlets. To determine whether DNA methylation is associated with gene expression changes in response to phytoplasma, we investigated variations in genomic DNA sequence and methylation in PaWB plantlets treated with methyl methane sulfonate (MMS) using amplified fragment length polymorphism (AFLP) and methylation-sensitive amplification polymorphism (MSAP) techniques, respectively. The results indicated that PaWB seedings recovered a normal morphology after treatment with more than 15 mg·L(-1) MMS. PaWB infection did not cause changes of the paulownia DNA sequence at the AFLP level; However, DNA methylation levels and patterns were altered. Quantitative real-time PCR (qRT-PCR) showed that three of the methylated genes were up-regulated and three were down-regulated in the MMS-treated PaWB plantlets that had regained healthy morphology. These six genes might be involved in transcriptional regulation, plant defense, signal transduction and energy. The possible roles of these genes in PaWB are discussed. The results showed that changes of DNA methylation altered gene expression levels, and that MSAP might help identify genes related to PaWB.

Differential transcriptome analysis between Paulownia fortunei and its synthesized autopolyploid

Paulownia fortunei is an ecologically and economically important tree species that is widely used as timber and chemical pulp. Its autotetraploid, which carries a number of valuable traits, was successfully induced with colchicine. To identify differences in gene expression between P. fortunei and its synthesized autotetraploid, we performed transcriptome sequencing using an Illumina Genome Analyzer IIx (GAIIx). About 94.8 million reads were generated and assembled into 383,056 transcripts, including 18,984 transcripts with a complete open reading frame. A conducted Basic Local Alignment Search Tool (BLAST) search indicated that 16,004 complete transcripts had significant hits in the National Center for Biotechnology Information (NCBI) non-redundant database. The complete transcripts were given functional assignments using three public protein databases. One thousand one hundred fifty eight differentially expressed complete transcripts were screened through a digital abundance analysis, including transcripts involved in energy metabolism and epigenetic regulation. Finally, the expression levels of several transcripts were confirmed by quantitative real-time PCR. Our results suggested that polyploidization caused epigenetic-related changes, which subsequently resulted in gene expression variation between diploid and autotetraploid P. fortunei. This might be the main mechanism affected by the polyploidization. Our results represent an extensive survey of the P. fortunei transcriptome and will facilitate subsequent functional genomics research in P. fortunei. Moreover, the gene expression profiles of P. fortunei and its autopolyploid will provide a valuable resource for the study of polyploidization.