Fig. 5

Relationship among m1A, m5C, m6A, m7G and miRNA-binding sites within 3’UTRs. A Association between 3’UTR methylation abundance and miRNA targeting. Methylation abundances of the targeted 3’UTRs for the 50 most/least abundant miRNAs in brain are shown using a bar chart, revealing a significantly greater percentage of m6A peaks within target mRNA 3‘UTRs for the 50 weakly expressed brain miRNAs compared with the 50 most abundant brain miRNAs (*p < 0.05, Wilcoxon test). The error bars indicate the highest and lowest values, and the box boundaries denote the first quartile, median, and third quartile (Public data: SRX685388, SRX685392). B (Upper) Venn diagram showing the overlap of genes comprising different kinds of methylation in their 3’UTR. (Middle) Bar chart displaying the number of M1a, m5C, m6A and m7G peaks within the 3’UTRs of methylation genes. (Lower) The number of 3’UTRs containing one, two, three or four kinds of methylation modification. C The enrichment scores of m1A, m5C, m6A and m7G of the five genes containing all four kinds of modifications are shown in the radar chart (uvrag: UV radiation resistance-associated gene; rab33ba: rab3b, member RAS oncogene family a; akap17a: A kinase (PRKA) anchor protein 17A). D Distribution of methylation peaks and miRNA target sites within 3’ UTRs. The frequency of methylation peaks under normoxia (red) or hypoxia (blue) and miRNA target sites (black) along the length of 3′ UTRs is shown. The relativity of methylation peak sites and miRNA target sites was calculated and displayed using the Pearson correlation coefficient