Identification of key genes associated with the effect of osmotic stimuli on intervertebral discs using microarray analysis.

The present study aimed to explore the effect of osmotic stimuli on intervertebral discs (IVDs) using microarray analysis. Gene expression dataset GSE1648 was downloaded from the Gene Expression Omnibus database. There were 11 IVD cell samples in this dataset, which included 4 hyperosmotic stimuli samples, 3 hypoosmotic stimuli samples and 4 isosmotic stimuli samples. The differentially expressed genes (DEGs) in hyperosmotic or hypoosmotic IVD cells (designated DEGs-hyper or DEGs-hypo) were identified, compared with isosmotic cells, using the limma package of R software. The Database for Annotation, Visualization and Integrated Discovery was used to perform a Gene Ontology (GO) term enrichment analysis for the DEG sets. Protein-protein interaction (PPI) network and microRNA (miRNA) gene-regulatory network data for the DEG sets were obtained using the Human Protein Reference Database (HPRD) and the TargetScan database, respectively, and these networks were constructed and visualized using Cytoscape software. There was a total of 43 DEGs in DEGs-hyper and 9 in DEGs-hypo. Analysis of DEGs-hyper revealed that 41 GO terms were significantly enriched. In total, 376 pairs and 382 nodes were involved in the PPI network, and 1,314 miRNA-gene pairs and 422 nodes were contained in the miRNA-gene-regulated network. The results of the present study indicated that potential target genes (including NCOA3, SOS1, XPO1, ZBTB18, EFNB2 and SOBP) may be involved in the effect of osmotic stimuli on IVD, and the biological processes of apoptosis and cell death may be associated with the effect of high osmolality on IVD disease. The potential targets identified in the present study are more reliable than those identified by previous studies.