Meat quality is one of the most important economic traits in pig breeding and production. Intramuscular fat (IMF) is a major factor that improves meat quality. To better understand the alternative splicing (AS) events underlying meat quality, long-read isoform sequencing (Iso-seq) was used to identify differential (D)AS events between the longissimus thoracis (LT) and semitendinosus (ST), which differ in IMF content, together with short-read RNA-seq. Through Iso-seq analysis, we identified a total of 56,789 novel transcripts covering protein-coding genes, lncRNA, and fusion transcripts that were not previously annotated in pigs. We also identified 456,965 AS events, among which 3930 were DAS events, corresponding to 2364 unique genes. Through integrative analysis of Iso-seq and RNA-seq, we identified 1174 differentially expressed genes (DEGs), among which 122 were DAS genes, i.e., DE-DAS genes. There are 12 overlapped pathways between the top 20 DEGs and DE-DAS genes, as revealed by KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis, indicating that DE-DAS genes play important roles in the differential phenotype of LT and ST. Further analysis showed that upregulated DE-DAS genes are more important than downregulated ones in IMF deposition. Fatty acid degradation and the PPAR (peroxisome proliferator-activated receptor) signaling pathway were found to be the most important pathways regulating the differential fat deposition of the two muscles. The results update the existing porcine genome annotations and provide data for the in-depth exploration of the mechanisms underlying meat quality and IMF deposition.

The cAMP response element-binding protein (CREB), a basic leucine zipper transcription factor, is involved in the activation of numerous genes in a variety of cell types. The CREB gene is rich in alternative splicing (AS) events. However, studies on the AS of CREB genes in pigs are limited, and few reports have compared the roles of isoforms in activating gene expression. Here, five AS transcripts, V1–5, were characterized by RT-PCR and two, V3 and V5, were new identifications. Both V1 and V2 have all the functional domains of the CREB protein, with similar tissue expression profiles and mRNA stability, suggesting that they have similar roles. The transcriptional transactivation activities of four isoforms encoding complete polypeptides were analyzed on the expression of the B-cell CLL/lymphoma 2-like protein 2 and the poly (A)-binding protein, nuclear 1 genes with a dual-luciferase reporter system, and differential activities were observed. Both V1 and V2 have promoting effects, but their roles are gene-specific. V3 has no effect on the promoter of the two genes, while V4 functions as a repressor. The mechanisms underlying the differential roles of V1 and V2 were analyzed with RNA-seq, and the genes specifically regulated by V1 and V2 were identified. These results will contribute to further revealing the role of CREB and to analyzing the significance of AS in genes.