Abstract
To elucidate the currently unknown mechanisms responsible for the diverse biological aspects between two-dimensional (2D) and three-dimensional (3D) cultured 3T3-L1 preadipocytes, RNA-sequencing analyses were performed. During a 7-day culture period, 2D- and 3D-cultured 3T3-L1 cells were subjected to lipid staining by BODIPY, qPCR for adipogenesis related genes, including peroxisome proliferator-activated receptor γ (Pparγ), CCAAT/enhancer-binding protein alpha (Cebpa), Ap2 (fatty acid-binding protein 4; Fabp4), leptin, and AdipoQ (adiponectin), and RNA-sequencing analysis. Differentially expressed genes (DEGs) were detected by next-generation RNA sequencing (RNA-seq) and validated by a quantitative reverse transcription-polymerase chain reaction (qRT-PCR). Bioinformatic analyses were performed on DEGs using a Gene Ontology (GO) enrichment analysis and an Ingenuity Pathway Analysis (IPA). Significant spontaneous adipogenesis was observed in 3D 3T3-L1 spheroids, but not in 2D-cultured cells. The mRNA expression of Pparγ, Cebpa, and Ap2 among the five genes tested were significantly higher in 3D spheroids than in 2D-cultured cells, thus providing support for this conclusion. RNA analysis demonstrated that a total of 826 upregulated and 725 downregulated genes were identified as DEGs. GO enrichment analysis and IPA found 50 possible upstream regulators, and among these, 6 regulators-transforming growth factor β1 (TGFβ1), signal transducer and activator of transcription 3 (STAT3), interleukin 6 (IL6), angiotensinogen (AGT), FOS, and MYC-were, in fact, significantly upregulated. Further analyses of these regulators by causal networks of the top 14 predicted diseases and functions networks (IPA network score indicated more than 30), suggesting that STAT3 was the most critical upstream regulator. The findings presented herein suggest that STAT3 has a critical role in regulating the unique biological properties of 3D spheroids that are produced from 3T3-L1 preadipocytes.