N6-methyladenosine modification changes during the recovery processes for Paulownia witches' broom disease under the methyl methanesulfonate treatment.

Phytoplasmas induce diseases in more than 1000 plant species and cause substantial ecological damage and economic losses, but the specific pathogenesis of phytoplasma has not yet been clarified. N 6-methyladenosine (m6A) is the most common internal modification of the eukaryotic Messenger RNA (mRNA). As one of the species susceptible to phytoplasma infection, the pathogenesis and mechanism of Paulownia has been extensively studied by scholars, but the m6A transcriptome map of Paulownia fortunei (P. fortunei) has not been reported. Therefore, this study aimed to explore the effect of phytoplasma infection on m6A modification of P. fortunei and obtained the whole transcriptome m6A map in P. fortunei by m6A-seq. The m6A-seq results of Paulownia witches' broom (PaWB) disease and healthy samples indicate that PaWB infection increased the degree of m6A modification of P. fortunei. The correlation analysis between the RNA-seq and m6A-seq data detected that a total of 315 differentially methylated genes were predicted to be significantly differentially expressed at the transcriptome level. Moreover, the functions of PaWB-related genes were predicted by functional enrichment analysis, and two genes related to maintenance of the basic mechanism of stem cells in shoot apical meristem were discovered. One of the genes encodes the receptor protein kinase CLV2 (Paulownia_LG2G000076), and the other gene encodes the homeobox transcription factor STM (Paulownia_LG15G000976). In addition, genes F-box (Paulownia_LG17G000760) and MSH5 (Paulownia_LG8G001160) had exon skipping and mutually exclusive exon types of alternative splicing in PaWB-infected seedling treated with methyl methanesulfonate, and m6A modification was found in m6A-seq results. Moreover, Reverse Transcription-Polymerase Chain Reaction (RT-PCR) verified that the alternative splicing of these two genes was associated with m6A modification. This comprehensive map provides a solid foundation for revealing the potential function of the mRNA m6A modification in the process of PaWB. In future studies, we plan to verify genes directly related to PaWB and methylation-related enzymes in Paulownia to elucidate the pathogenic mechanism of PaWB caused by phytoplasma invasion.