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.