Background N6-methyladenosine (m6A) is the most common RNA modification in eukaryotes and has been implicated as a novel epigenetic marker that is involved in various biological processes. The pattern and functional dissection of m6A in the regulation of several major human viral diseases have already been reported. However, the patterns and functions of m6A distribution in plant disease bursting remain largely unknown. Results We analyse the high-quality m6A methylomes in rice plants infected with two devastating viruses. We find that the m6A methylation is mainly associated with genes that are not actively expressed in virus-infected rice plants. We also detect different m6A peak distributions on the same gene, which may contribute to different antiviral modes between rice stripe virus or rice black-stripe dwarf virus infection. Interestingly, we observe increased levels of m6A methylation in rice plant response to virus infection. Several antiviral pathway-related genes, such as RNA silencing-, resistance-, and fundamental antiviral phytohormone metabolic-related genes, are also m6A methylated. The level of m6A methylation is tightly associated with its relative expression levels. Conclusions We revealed the dynamics of m6A modification during the interaction between rice and viruses, which may act as a main regulatory strategy in gene expression. Our investigations highlight the significance of m6A modifications in interactions between plant and viruses, especially in regulating the expression of genes involved in key pathways.

No abstract

No abstract

Broad bean wilt virus 2 (BBWV-2), which belongs to the genus Fabavirus of the family Secoviridae, is an important pathogen that causes damage to broad bean, pepper, yam, spinach and other economically important ornamental and horticultural crops worldwide. Previously, only limited reports have shown the genetic variation of BBWV2. Meanwhile, the detailed evolutionary changes, synonymous codon usage bias and host adaptation of this virus are largely unclear. Here, we performed comprehensive analyses of the phylodynamics, reassortment, composition bias and codon usage pattern of BBWV2 using forty-two complete genome sequences of BBWV-2 isolates together with two other full-length RNA1 sequences and six full-length RNA2 sequences. Both recombination and reassortment had a significant influence on the genomic evolution of BBWV2. Through phylogenetic analysis we detected three and four lineages based on the ORF1 and ORF2 nonrecombinant sequences, respectively. The evolutionary rates of the two BBWV2 ORF coding sequences were 8.895 × 10−4 and 4.560 × 10−4 subs/site/year, respectively. We found a relatively conserved and stable genomic composition with a lower codon usage choice in the two BBWV2 protein coding sequences. ENC-plot and neutrality plot analyses showed that natural selection is the key factor shaping the codon usage pattern of BBWV2. Strong correlations between BBWV2 and broad bean and pepper were observed from similarity index (SiD), codon adaptation index (CAI) and relative codon deoptimization index (RCDI) analyses. Our study is the first to evaluate the phylodynamics, codon usage patterns and adaptive evolution of a fabavirus, and our results may be useful for the understanding of the origin of this virus.