Gene co-expression and alternative splicing analysis of key metabolic tissues to unravel the regulatory signatures of fatty acid composition in cattle.

Increasing the healthy/unhealthy fatty acid (FA) ratio in meat is one of the urgent tasks required to address consumer concerns. However, the regulatory mechanisms ultimately resulting in FA profiles vary among animals and remain largely unknown. In this study, using ~1.2 Tb high-quality RNA-Seq based transcriptomic data of 188 samples from four key metabolic tissues (rumen, liver, muscle, and backfat) together with the contents of 49 FAs in backfat, the molecular regulatory mechanisms of these tissues contributing to FA formation in cattle were explored. Using this large dataset, the alternative splicing (AS) events, one of the transcriptional regulatory mechanisms in four tissues were identified. The highly conserved and absent AS events were detected in rumen tissue, which may contribute to its functional differences compared with the other three tissues. In addition, the healthy/unhealthy FA ratio related AS events, differential expressed (DE) genes, co-expressed genes, and their functions in four tissues were analyzed. Eight key genes were identified from the integrated analysis of DE, co-expressed, and AS genes between animals with high and low healthy/unhealthy FA ratios. This study provides an applicable pipeline for AS events based on comprehensive RNA-Seq analysis and improves our understanding of the regulatory mechanism of FAs in beef cattle.