Fig. 1
RNA methylation and expression variability is inversely correlated. a Density distribution of m6A in the control and heat-stressed leaves and flowers of Arabidopsis. UTR and CDS are marked as dark and light grey boxes, respectively. CS, control sample; HS, heat-stressed sample. b Consensus sequence motif of m6A-modified sites in the non-stressed leaf sample. The top 1000 most m6A-enriched regions were used for the analysis. c Number of genes with m6A enrichment score greater than 0. Enrichment of m6A was determined by log2-converted fold change of m6A-immunoprecipitated to input RNA levels. d, e Number of HVGs (d) and LVGs (e) in leaves and flowers subjected to the control and heat stress treatments. HVGs and LVGs are genes with log2CV2 greater than 5 and lower than − 10, respectively. CV2 refers to coefficient of variation of a gene which was determined by FPKMs of RNA-seq performed for three biological replications. f Fraction of m6A-modified transcripts in quartiles of gene expression variability. Methylated transcript was as defined in c