Supporting Information
Figure S1. Experimental flow chart of nanopore direct RNA sequencing.
Figure S2. Overall transcriptomic changes in rice in response to BPH infestation.
Figure S3. The overall m6A sites distribution along the mRNA in each group.
Figure S4. Histograms show the number of m6A methylation positions in un-infested and BPH-infested plants.
Figure S5. Differentiated m6A modification in BPH-infested rice plants.
Figure S6. Consensus motifs containing m6A methylation sites in BPH-infested and un-infested plants.
Figure S7. Euclidian distance coefficients among gene transcript profiles based on DRS analysis of two treatments.
Figure S8. Correlation analysis between m6A methylation and transcription abundance of target gene bodies in rice under BPH infestation.
Figure S9. Correlation analysis of m6A modification and differentially expressed transcripts in BPH-infested rice plants.
Figure S10. The number of gene bodies decreased along with the m6A modification sites increased.
Figure S11. GO and KEGG analyses of differentially m6A methylated genes in BPH-infested rice plants.
Figure S12. Expression analysis of SA pathway-related genes and SA quantification in BPH-infested rice.
Figure S13. RNA loading control that corresponding to the dot-blot analyses.
Table S1. Distribution of m6A modification sites on rice chromosomes in each treatment group.
Table S2. The number of m6A modification sites in each treatment group.
Table S3. Information of specific m6A methylation positions in Nl-Nip vs. Nip comparison.
Table S4. Motifs of m6A modification sites in each treatment group.
Table S5. Transcripts Per kilobase of exon model per Million mapped reads (TPM) of all genes in each treatment group.
Table S6. Integrated analyses of m6A modification sites, m6A directions, transcriptome expression levels, and transcript regulatory direction.
Table S7. Detailed information of the m6A methylated genes appeared in enriched GO pathways under BPH infestation.
Table S8. Detailed information of the m6A methylated genes appeared in enriched KEGG pathways under BPH infestation.
Table S9.m6A modification profile of the main m6A modification machinery components.
Table S10. Analysis of m6A position and direction of the main m6A modification machinery components genes in Nip and Nl-Nip treatment groups.
Table S11.m6A modification profile of the brown planthopper resistance-, cellulose and hemicellulose synthesis-related genes.
Table S12. Analysis of m6A position and direction of the BPH resistance genes in Nip and Nl-Nip treatment groups.
Table S13. Analysis of m6A position and direction of the cellulose synthesis genes in Nip and Nl-Nip treatment groups.
Table S14. Analysis of m6A position and direction of the hemicellulose synthesis genes in Nip and Nl-Nip treatment groups.
Table S15. m6A modification profile of the JA- and SA-metabolism related genes.
Table S16. Analysis of m6A position and direction of the JA biosynthesis genes in Nip and Nl-Nip treatment groups.
Table S17. Analysis of m6A position and direction of the JA response genes in Nip and Nl-Nip treatment groups.
Table S18. Analysis of m6A position and direction of the SA biosynthesis genes in Nip and Nl-Nip treatment groups.
Table S19. Analysis of m6A position and direction of the SA responsive genes in Nip and Nl-Nip treatment groups.
Table S20.m6A modification profile of rice growth-related genes.
Table S21. Analysis of m6A position and direction of auxin-related genes in Nip and Nl-Nip treatment groups.
Table S22. Analysis of m6A position and direction of gibberellin-related genes in Nip and Nl-Nip treatment groups.
Table S23. Primers used in RT-qPCR qualification of the candidate genes.