Analysis of gene expression profiles at different stages during preadipocyte differentiation in rabbits

Excessive accumulation of fat is harmful to human health. The preadipocyte differentiation is a critical process of fat development. Studying the expression profiles of genes related to preadipocyte differentiation contributes to understanding of the mechanism of fat accumulation. Despite being considered an ideal animal model for studying adipogenesis, little is known about the gene expression profiles at different stages during preadipocyte differentiation in rabbits. In the present study, rabbit preadipocytes were cultured in vitro and induced for differentiation, and gene expression profiles of adipocytes collected at days 0, 3, and 9 of differentiation were analyzed by RNA-seq. We identified 1352 differentially expressed genes (DEGs) when comparing day 3 with day 0 and identified 888 DEGs when comparing day 9 with day 3. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that the PPAR signaling pathway and PI3K-Akt signaling pathway were significantly enriched by the DEGs that up-regulated within the period of day 0 - day 3, and the GO terms and KEGG pathways that were associated with cell cycle were enriched by the DEGs that up-regulated within the period of day 3 - day 9. The DEGs that specifically up-regulated within the period of day 0 - day 3 might play roles in the cytoplasm, and the DEGs that specifically up-regulated within the period of day 3 - day 9 might act in the nucleus. The protein-protein interaction (PPI) network constructed by DEGs showed that hub node genes might modulate rabbit preadipocyte differentiation via regulating cell cycle.