Abstract
Fast-twitch and slow-twitch muscle fibers not only differ in metabolic characteristics and physiological functions but also significantly influence the texture of livestock meat. RNA editing represents an important post-transcriptional regulatory process that can influence both gene expression and the resulting protein function. However, studies on RNA editing events in yak muscle remain limited. This study systematically identified RNA editing events in yak biceps femoris (BF, n = 3) and obliquus externus abdominis (OEA, n = 3) using transcriptomic data, discovering 17,713 unique editing sites, most located in non-coding regions. Within coding regions, 3350 sites were detected, with 1195 resulting in non-synonymous amino acid substitutions. Further analysis revealed that 785 sites potentially affected miRNA binding sites, suggesting RNA editing may participate in miRNA-mediated gene regulation. Tukey's post hoc test (p < 0.05) identified 242 sites (involving 170 genes) with significantly different editing levels between BF and OEA. KEGG pathway analysis indicated that genes with differential RNA editing were predominantly associated with pathways involved in muscle fiber type transitions, including the MAPK and calcium signaling pathways. Collectively, this study maps the RNA editing landscape in yak muscle tissue and identifies distinct, fiber-type-specific RNA editing patterns between oxidative and glycolytic muscle fibers, including differences in editing levels and site distributions, supporting a potential association between RNA editing and muscle fiber type transformation.