Transcriptome analysis of buffalo granulosa cells in three dimensional culture systems

Hanging drop (HD) three-dimensional (3D) culture model for buffalo granulosa cells (GC) was reported to mimic the preovulatory stage of ovarian follicles in our previous study. To further verify its reliability, the present study attempted a comparative transcriptome profile of buffalo GC freshly isolated from ovarian follicles (<8 mm diameter) (FC) and their cultures in normal culture dish (ND or 2D), polyHEMA coated dish (PH) and HD culture systems (3D). Out of 223 significantly (-log2 fold change: >3; p < .0005; false discovery rate [FDR]: <0.1) differentially expressed genes (SDEGs) among different culture systems, 137 were found unannotated, and 94, 29, and 66 were exclusively expressed in FC, PH, and HD, respectively. However, on eliminating the fixed points of p values and FDR from the entire raw data, only 11 genes related to long noncoding RNA, 12 genes related to luteinization, and 3 genes related to follicular maturation were exclusively expressed in FC, PH, and HD culture systems, respectively. The quantitative real time-PCR validation and the next generation sequencing data had more than 90% correlation. Bioinformatics analyses of the exclusively expressed SDEG revealed that the freshly aspirated GCs were a true representative of GCs from small follicles (<8 mm diameter), the GC spheroids under PH maintained mitochondrial function, and those cultured in HD system for 6 days simulated the inflammatory milieu required for ovulation. Therefore, the comparative transcriptome profile also reinforced that HD culture system is better in vitro culture method than the other methods analyzed in this study for buffalo GC.