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

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.

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

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

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.

Transcriptomic Analyses of Phytoplasmas

Transcriptomic analyses addressed to study phytoplasma gene expression may present few difficulties due to the uncultivable nature of these intracellular, obligate pathogens. While RNA extraction from insect vectors does not imply any particular adaptation of the protocols used in most commercial kits, RNA isolation from phytoplasma-infected plants can be a challenging task, given the high levels of polyphenol contents and accumulation of sucrose and starch in the different plant tissues. Here, we describe two different transcriptomic approaches, one focused on RNA phytoplasma sequencing and the other on phytoplasma quantitative gene expression in relation to pathogen load.

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.

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.

Genome of Paulownia (Paulownia fortunei) illuminates the related transcripts, miRNA and proteins for salt resistance

Polyploidy in plants can bestow long-term evolutionary flexibility and resistance to biotic and abiotic stresses. The upstream activation mechanisms of salt response remain unknown. Here we integrated transcriptome, miRNA and proteome data to describe the link between abscisic acid (ABA)-effectors and salt resistance against the background of Paulownia genome. Combing GO and KEGG pathway annotation of differentially expressed genes and proteins, as well as differentially expressed miRNA, these results reflect endogenous signal ABA activate the downstream effectors, such as ion channel effectors and oxido-reduction effectors, to maintain the homeostasis of Paulownia’s growth. The cascaded metabolic network involved ABA biosynthesis, signaling transduction and the response of effectors. Our results will contribute to a comprehensive understanding of the genetic basis of salt tolerance, which may help to expand the available arable land for P. fortunei cultivation.

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.

Transcriptome-Wide Analysis of Botrytis elliptica Responsive microRNAs and Their Targets in Lilium Regale Wilson by High-Throughput Sequencing and Degradome Analysis

MicroRNAs, as master regulators of gene expression, have been widely identified and play crucial roles in plant-pathogen interactions. A fatal pathogen, Botrytis elliptica, causes the serious folia disease of lily, which reduces production because of the high susceptibility of most cultivated species. However, the miRNAs related to Botrytis infection of lily, and the miRNA-mediated gene regulatory networks providing resistance to B. elliptica in lily remain largely unexplored. To systematically dissect B. elliptica-responsive miRNAs and their target genes, three small RNA libraries were constructed from the leaves of Lilium regale, a promising Chinese wild Lilium species, which had been subjected to mock B. elliptica treatment or B. elliptica infection for 6 and 24 h. By high-throughput sequencing, 71 known miRNAs belonging to 47 conserved families and 24 novel miRNA were identified, of which 18 miRNAs were downreguleted and 13 were upregulated in response to B. elliptica. Moreover, based on the lily mRNA transcriptome, 22 targets for 9 known and 1 novel miRNAs were identified by the degradome sequencing approach. Most target genes for elliptica-responsive miRNAs were involved in metabolic processes, few encoding different transcription factors, including ELONGATION FACTOR 1 ALPHA (EF1a) and TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR 2 (TCP2). Furthermore, the expression patterns of a set of elliptica-responsive miRNAs and their targets were validated by quantitative real-time PCR. This study represents the first transcriptome-based analysis of miRNAs responsive to B. elliptica and their targets in lily. The results reveal the possible regulatory roles of miRNAs and their targets in B. elliptica interaction, which will extend our understanding of the mechanisms of this disease in lily.

Transcriptome and Degradome of microRNAs and Their Targets in Response to Drought Stress in the Plants of a Diploid and Its Autotetraploid Paulownia australis

MicroRNAs (miRNAs) are small, non-coding RNAs that play vital roles in plant growth, development, and stress response. Increasing numbers of studies aimed at discovering miRNAs and analyzing their functions in plants are being reported. In this study, we investigated the effect of drought stress on the expression of miRNAs and their targets in plants of a diploid and derived autotetraploid Paulownia australis. Four small RNA (sRNA) libraries and four degradome libraries were constructed from diploid and autotetraploid P. australis plants treated with either 75% or 25% relative soil water content. A total of 33 conserved and 104 novel miRNAs (processing precision value > 0.1) were identified, and 125 target genes were identified for 36 of the miRNAs by using the degradome sequencing. Among the identified miRNAs, 54 and 68 were differentially expressed in diploid and autotetraploid plants under drought stress (25% relative soil water content), respectively. The expressions of miRNAs and target genes were also validated by quantitative real-time PCR. The results showed that the relative expression trends of the randomly selected miRNAs were similar to the trends predicted by Illumina sequencing. And the correlations between miRNAs and their target genes were also analyzed. Furthermore, the functional analysis showed that most of these miRNAs and target genes were associated with plant development and environmental stress response. This study provided molecular evidence for the possible involvement of certain miRNAs in the drought response and/or tolerance in P. australis, and certain level of differential expression between diploid and autotetraploid plants.

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.

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.

Small RNA Sequencing Based Identification of MiRNAs in Daphnia magna

Small RNA molecules are short, non-coding RNAs identified for their crucial role in post-transcriptional regulation. A well-studied example includes miRNAs (microRNAs) which have been identified in several model organisms including the freshwater flea and planktonic crustacean Daphnia. A model for epigenetic-based studies with an available genome database, the identification of miRNAs and their potential role in regulating Daphnia gene expression has only recently garnered interest. Computational-based work using Daphnia pulex, has indicated the existence of 45 miRNAs, 14 of which have been experimentally verified. To extend this study, we took a sequencing approach towards identifying miRNAs present in a small RNA library isolated from Daphnia magna. Using Perl codes designed for comparative genomic analysis, 815,699 reads were obtained from 4 million raw reads and run against a database file of known miRNA sequences. Using this approach, we have identified 205 putative mature miRNA sequences belonging to 188 distinct miRNA families. Data from this study provides critical information necessary to begin an investigation into a role for these transcripts in the epigenetic regulation of Daphnia magna.