Conclusion
Our study revealed different DEGs, pathways, biological process and PPI networks involved in post-traumatic IVD degeneration. These findings will advance the understanding of the pathogenesis of IVD degeneration, and help to identify novel biomarkers for the disease diagnosis.
Methods
To mimic a traumatic insult, bovine caudal IVDs underwent one strike loading. The control group was cultured under physiological loading. At 24 hours after one strike or physiological loading, RNA was extracted from nucleus pulposus (NP) and annulus fibrosus (AF) tissue. High throughput next generation RNA sequencing was performed to identify differentially expressed genes (DEGs) between the one strike loading group and the control group. Gene Ontology (GO) functional and Kyoto Encyclopedia of Genes and Genomes analyses were performed to analyze DEGs and pathways. Protein-protein interaction (PPI) network was analyzed with cytoscape software. DEGs were verified using qRT-PCR. Degenerated human IVD tissue was collected for immunofluorescence staining to verify the expression of DEGs in human disc tissue.
Results
One strike loading resulted in significant gene expression changes compared with physiological loading. In total 253 DEGs were found in NP tissue and 208 DEGs in AF tissue. Many of the highly dysregulated genes have known functions in disc degeneration and extracellular matrix (ECM) homeostasis. ACTB, ACTG, PFN1, MYL12B in NP tissue and FGF1, SPP1 in AF tissue were verified by qRT-PCR and immunofluorescence imaging. The identified DEGs were involved in focal adhesion, ECM-receptor interaction, PI3K-AKT, and cytokine-cytokine receptor interaction pathways. Three clusters of PPI networks were identified. GO enrichment revealed that these DEGs were mainly involved in inflammatory response, the ECM and growth factor signaling and protein folding biological process.