Figure 3
The evaluation of RNA editing sites by distribution in different cell types and genomic elements and their dynamic changes during hematopoiesis. (A) The number of RNA editing sites in each HSPC population. (B) Pearson correlation coefficient between the total reads count and the number of editing sites in each HSPC population. (C) Pearson correlation coefficient between the cell number and the number of editing sites in each HSPC population. The number of RNA editing sites were independent of cell number and sequencing depth. (D) Bar plot showing the proportion of editing sites within and outside Alus in HSPC populations. Most RNA editing sites located in Alu genomics element. Data present mean ± s.e.m. P-value was calculated by T-Test. ****P < 0.0001. (E) Bar plot showing the distribution of editing sites across different genomic elements in HSPC populations. Data present mean ± s.e.m. (F) A heatmap showing confidence scores (P_edit) calculated by RED-ML of all editing sites. The color indicates P_edit (from 0 to 1) for a given site (row) in a population (column). The P_edit equal to 0 means this site is not edited in this population whose coverage is greater than 30×. The right listed the most significant enriched GO pathways in each pattern and the common pathways in all patterns. These GO terms were analysed with genes of these RNA editing sites annotated by ANNOVAR. (G) A heatmap showing the most significant enriched GO pathways in each module related to hematopoiesis. The genes annotated by ANNOVAR were used to do GO analysis.
