m6A mRNA methylation analysis provides novel insights into heat stress responses in the liver tissue of sheep

N6-methyladenosine (m6A) mRNA methylation varies in response to stress. However, no map of m6A mRNA methylation has been obtained for sheep, nor is it known what effect this has on regulating heat stress in sheep. Here, we obtained m6A methylation maps of sheep liver tissues with and without heat stress by MeRIP-seq. In total, 8306 m6A peaks associated with 2697 genes were detected in the heat stress group, and 12,958 m6A peaks associated with 5494 genes were detected in the control group. Peaks were mainly enriched in coding regions and near stop codons with classical RRACH motifs. Methylation levels of heat stress and control sheep were higher near stop codons, although methylation was significantly lower in heat stress sheep. GO and KEGG revealed that differential m6A-containing genes were significantly enriched in the stress response and fat metabolism. Our results showed that m6A mRNA methylation modifications regulate heat stress in sheep.