Aims. To explore the novel molecular mechanisms of histone deacetylase 4 (HDAC4) in chondrocytes via RNA sequencing (RNA-seq) analysis. Methods. Empty adenovirus (EP) and a HDAC4 overexpression adenovirus were transfected into cultured human chondrocytes. The cell survival rate was examined by real-time cell analysis (RTCA) and EdU and flow cytometry assays. Cell biofunction was detected by Western blotting. The expression profiles of messenger RNAs (mRNAs) in the EP and HDAC4 transfection groups were assessed using whole-transcriptome sequencing (RNA-seq). Volcano plot, Gene Ontology, and pathway analyses were performed to identify differentially expressed genes (DEGs). For verification of the results, the A289E/S246/467/632 A sites of HDAC4 were mutated to enhance the function of HDAC4 by increasing HDAC4 expression in the nucleus. RNA-seq was performed to identify the molecular mechanism of HDAC4 in chondrocytes. Finally, the top ten DEGs associated with ribosomes were verified by quantitative polymerase chain reaction (QPCR) in chondrocytes, and the top gene was verified both in vitro and in vivo. Results. HDAC4 markedly improved the survival rate and biofunction of chondrocytes. RNA-seq analysis of the EP and HDAC4 groups showed that HDAC4 induced 2,668 significant gene expression changes in chondrocytes (1,483 genes upregulated and 1,185 genes downregulated, p < 0.05), and ribosomes exhibited especially large increases. The results were confirmed by RNA-seq of the EP versus mutated HDAC4 groups and the validations in vitro and in vivo. Conclusion. The enhanced ribosome pathway plays a key role in the mechanism by which HDAC4 improves the survival rate and biofunction of chondrocytes. Cite this article: Bone Joint Res 2023;12(7):433–446

Objectives. The aim of this study was to provide a comprehensive understanding of alterations in messenger RNAs (mRNAs), long noncoding RNAs (lncRNAs), and circular RNAs (circRNAs) in cartilage affected by osteoarthritis (OA). Methods. The expression profiles of mRNAs, lncRNAs, and circRNAs in OA cartilage were assessed using whole-transcriptome sequencing. Bioinformatics analyses included prediction and reannotation of novel lncRNAs and circRNAs, their classification, and their placement into subgroups. Gene ontology and pathway analysis were performed to identify differentially expressed genes (DEGs), differentially expressed lncRNAs (DELs), and differentially expressed circRNAs (DECs). We focused on the overlap of DEGs and targets of DELs previously identified in seven high-throughput studies. The top ten DELs were verified by quantitative reverse transcriptase polymerase chain reaction (qRT-PCR) in articular chondrocytes, both in vitro and in vivo. Results. In total, 739 mRNAs, 1152 lncRNAs, and 42 circRNAs were found to be differentially expressed in OA cartilage tissue. Among these, we identified 18 overlapping DEGs and targets of DELs, and the top ten DELs were screened by expression profile analysis as candidate OA-related genes. WISP2, ATF3, and CHI3L1 were significantly increased in both normal versus OA tissues and normal versus interleukin (IL)-1β-induced OA-like cell models, while ADAM12, PRELP, and ASPN were shown to be significantly decreased. Among the identified DELs, we observed higher expression of ENST00000453554 and MSTRG.99593.3, and lower expression of MSTRG.44186.2 and NONHSAT186094.1 in normal versus OA cells and tissues. Conclusion. This study revealed expression patterns of coding and noncoding RNAs in OA cartilage, which added sets of genes and noncoding RNAs to the list of candidate diagnostic biomarkers and therapeutic agents for OA patients. Cite this article: H. Li, H. H. Yang, Z. G. Sun, H. B. Tang, J. K. Min. Whole-transcriptome sequencing of knee joint cartilage from osteoarthritis patients. Bone Joint Res 2019;8:290–303. DOI: 10.1302/2046-3758.87.BJR-2018-0297.R1