Full-Length Transcriptome and Gene Expression Analysis of Different Ovis aries Adipose Tissues Reveals Transcript Variants Involved in Lipid Biosynthesis

Animals (Basel). 2023 Dec 19;14(1):7. doi: 10.3390/ani14010007.

Abstract

Sheep have historically been bred globally as a vital food source. To explore the transcriptome of adipose tissue and investigate key genes regulating adipose metabolism in sheep, adipose tissue samples were obtained from F1 Dorper × Hu sheep. High-throughput sequencing libraries for second- and third-generation sequencing were constructed using extracted total RNA. Functional annotation of differentially expressed genes and isoforms facilitated the identification of key regulatory genes and isoforms associated with sheep fat metabolism. SMRT-seq generated 919,259 high-accuracy cDNA sequences after filtering. Full-length sequences were corrected using RNA-seq sequences, and 699,680 high-quality full-length non-chimeric (FLNC) reads were obtained. Upon evaluating the ratio of total lengths based on FLNC sequencing, it was determined that 36,909 out of 56,316 multiple-exon isoforms met the criteria for full-length status. This indicates the identification of 330,375 full-length FLNC transcripts among the 370,114 multiple-exon FLNC transcripts. By comparing the reference genomes, 60,276 loci and 111,302 isoforms were identified. In addition, 43,423 new genes and 44,563 new isoforms were identified. The results identified 185 (3198), 394 (3592), and 83 (3286) differentially expressed genes (transcripts) between tail and subcutaneous, tail and visceral, and subcutaneous and visceral adipose tissues, respectively. Functional annotation and pathway analysis revealed the following observations. (1) Among the differentially expressed genes (DEGs) of TF and SF tissues, the downregulation of ACADL, ACSL6, and NC_056060.1.2536 was observed in SF, while FFAR4 exhibited upregulation. (2) Among the DEGs of TF and VF tissues, expressions of ACADL, ACSL6, COL1A1, COL1A2, and SCD were downregulated in VF, with upregulation of FFAR4. (3) Among SF and VF expressions of COL1A1, COL1A2, and NC_056060.1.2536 were downregulated in VF. Specific differentially expressed genes (ACADL, ACSL6, COL1A1, COL1A2, FFAR4, NC_056060.1.2536, and SCD) and transcripts (NC_056066.1.1866.16 and NC_056066.1.1866.22) were identified as relevant to fat metabolism. These results provide a dataset for further verification of the regulatory pathway associated with fat metabolism in sheep.

Keywords: RNA sequencing; adipose tissue; full-length transcriptome; new transcript; sheep.