Intramuscular preadipocyte differentiation plays a critical role in bovine intramuscular fat (IMF) deposition. However, the roles of different RNAs, including mRNAs, circRNAs, lncRNAs and miRNAs, in regulating the adipogenic differentiation of intramuscular preadipocytes remain largely unclear.
In the present study, a whole transcriptome sequencing and analysis, including the analysis of mRNAs, circRNAs, lncRNAs and miRNAs, during different differentiation stages (0, 3, 6, and 9 d) of intramuscular preadipocytes from Qinchuan cattle was performed. All samples were prepared with 3 biological replicates. Here, a total of 27,153 mRNAs, 14,070 circRNAs, 7035 lncRNAs, and 427 miRNAs were annotated. Among them, we identified 4848 differentially expressed mRNAs (DEMs), 181 DE circRNAs (DECs), 501 DE lncRNAs (DELs) and 77 DE miRNAs (DEmiRs) between 0 d and other differentiation days (3, 6, and 9 d). GO and KEGG functional enrichment analyses showed that these differentially expressed genes were mainly enriched in cell differentiation, fat metabolism and adipogenesis-related pathways. Furthermore, weighted gene coexpression network analysis (WGCNA) and co-expression network analysis screened out multiple important mRNAs, circRNAs and lncRNAs related to intramuscular adipogenesis. Based on the competing endogenous RNA (ceRNA) regulatory mechanism, we finally identified 24 potential ceRNA networks and 31 potential key genes, including FOXO1/miR-330/circRNA2018/MSTRG.20301, GPAM/miR-27b/ciRNA489 and SESN3/miR-433/circRNA2627MSTRG.20342.
This study provides new insights into the differential expression patterns of different transcript types (i.e., mRNAs, circRNAs, lncRNAs and miRNAs) in intramuscular preadipocyte differentiation. Our findings provide data support for studying the molecular mechanism of key mRNAs and noncoding RNAs in IMF deposition, and provide new candidate markers for the molecular breeding of beef cattle.