35In livestock, glycolytic potential (GP) is a critical indicator for evaluating the meat 36 quality. To date, two major genes protein kinase AMP-activated γ3 non-catalytic 37 subunit gene (PRKAG3) and phosphorylase kinase catalytic subunit gamma 38 1(PHKG1), and corresponding cause mutations influencing GP have been confirmed 39 in pigs. Therefore, the aim of this study to identify the novel candidate genes and 40 variations related to GP-related traits using a four-hybrid pig model [Pietrain (P)× 41 Duroc (D)] ×[(Landrace) ×(Yorkshire)]. We totally constructed six RNA-seq libraries 42 using longissimus dorsi (LD) muscles, and each library contained two higher GP (H) 43 or two lower GP (L) individuals. A total of 525, 698 and 135 differentially expressed 44 genes (DEGs) were identified between H11 vs L11, H9 vs L9, and H5 vs L5 groups 45 using PossionDis method, respectively. Notably, we found 97 non-redundant DEGs 46 were mapped to GP related QTLs from three paired comparison groups. Moreover, 69 47 DEGs were identified between H (H11, H9 and H5) and L (L11, L9 and L5) groups 48 using NOIseq method. Additionally, 1,076 potential specific SNPs were figured out 49 between H and L groups, and approximately 40 large Indels with a length ≥ 5bp were 50 5identified in each sequencing library. In conclusion, our data provide foundation for 51 further confirming the key genes and the functional mutations affecting GP-related 52 traits in pigs, and also pave the way for elucidating the underling molecular regulatory 53 mechanisms of glycogen metabolism in future study. Moreover, this study might 54 provide valuable information for study on human glycogen storage diseases.