Abstract
Genetically identical female honeybee larvae with different diets develop into sterile workers or fertile queens. It remains unknown whether the reversible RNA N6-methyladenosine (m6A) mark functionally impact this "caste differentiation." Here, we profile the transcriptome-wide m6A methylome of honeybee queen and worker larvae at three instar stages and discover that m6A methylation dynamics are altered by differential feeding. Multiple methylome comparisons show an obvious increase in m6A marks during larval development and reveal a negative correlation between gene expression and m6A methylation. Notably, we find that worker larvae contain more hypermethylated m6A peaks than do queen larvae, and many caste-differentiation-related transcripts are differentially methylated. Chemical suppression of m6A methylation in worker larvae by 3-deazaadenosine (DAA) reduces overall m6A methylation levels and triggers worker larvae to develop queen caste features. Thus, our study demonstrates that m6A functionally impacts caste differentiation and larval development, yet it does not exclude potential contributions from other factors.