Screening for potential genes associated with bone overgrowth after mid-shaft femur fracture in a rat model

BACKGROUND We investigated the underlying molecular mechanisms of bone overgrowth after femoral fracture by using high-throughput bioinformatics approaches. METHODS The gene expression profile of GSE3298 (accession number) was obtained from the Gene Expression Omnibus database. Sixteen femoral growth plate samples, including nine samples without fracture and seven fracture samples for seven time points, were used for analysis. The Limma package was applied to identify differentially expressed genes (DEGs) between fractured and intact samples. The DAVID online tool was used for Gene ontology functional and pathway enrichment analysis. A protein-protein interaction (PPI) network established by String software was used to identify interactions between significant DEGs, and network modules were detected using plug-in MCODE. Additionally, a transcription regulatory network was constructed based on the ENCODE Project and PPI network. RESULTS A total of 680 DEGs were screened in fractured femoral growth plate samples compared with controls, including 238 up- and 442 down-regulated genes. These DEGs were significantly involved in the calcium signaling pathway and cancer pathway. A PPI network was constructed with 167 nodes and 233 edges, and module analysis demonstrated that CCL2, CSF2, NOS2, and DLC1 may stimulate bone overgrowth after femoral fracture via anti-apoptosis-related functions. A transcription regulatory network was constructed with 387 interacting pairs, and overlapping nodes were significantly enriched in intracellular signaling cascade and regulation of cell proliferation, among others. CONCLUSIONS Bone overgrowth was associated with changes in the expression of identified DEGs such as CCL2, NOS2, CSF2, and DLC1 in the femoral head. They may be important in regulating bone overgrowth via the anti-apoptosis of osteoblasts.