Abstract:Objective: To explore the clinical significance of differential expression of serum peripheral micrornas(miRNAs)in spinal tuberculosis. Methods: From May 2020 to May 2023, 50 patients with spinal tuberculosis (experimental group) and 30 healthy people (control group) admitted to the Department of Spinal Orthopaedics of our hospital were taken as the research objects. MiRCURYTM LNA array (v.18.0) was used to screen the expression profile of miRNAs in 10 experimental patients and 10 healthy controls. The expression levels of 10 selected miRNAs were verified by RT-qCPR in the research cohort. Hsa-miR-4708-3p, hsa-miR-5100 and hsa-miR-493-3p were selected to construct the spinal tuberculosis diagnosis model (STB model) based on miRNAs. Diagnostic value of STB model in spinal tuberculosis by analyzing receiver operating characteristic (ROC) curve. Results: Compared with healthy people, 288 kinds of miRNAs were found in the serum of patients with spinal tuberculosis(123 kinds were up-regulated and 165 kinds were down-regulated) (multiple change>2, P<0.05). RT-qPCR analysis confirmed that compared with the control group, the levels of hsa-miR-4708-3p and hsa-miR-5100 in serum samples of the experimental group increased significantly(P<0.05), and the level of hsa-miR-493-3p decreased significantly(P<0.05). Based on the expression levels of hsa-miR-4708-3p, hsa-miR-5100 and hsa-miR-493-3p as variables, a multivariate logistic regression model was established as follows: 0.582+1.055×(hsa-miR-4708-3p)+0.811× (hsa-miR-5100)-0.970×(hsa-miR-493-3p). Compared with the control group, the STB model score in the experimental group was significantly higher (0.32±0.04 vs 0.65±0.08, t=6.09, P<0.001). When the critical value was set to 0.46, the area under ROC curve (AUC) of STB model in identifying spinal tuberculosis was 0.911, with sensitivity of 88.9% and specificity of 77.0%. The STB model scores of patients with spinal tuberculosis with abnormal CRP and erythrocyte sedimentation rate were significantly higher than those with normal CRP and erythrocyte sedimentation rate (P<0.05). The score of STB model was positively correlated with CRP and erythrocyte sedimentation rate (r=0.633, 0.499, all P<0.001). Conclusion: Hsa-miR-4708-3p, hsa-miR-5100 and hsa-miR-493-3p have stable differential expressions in peripheral blood of patients with spinal tuberculosis, and the diagnostic model containing these three miRNAs has good diagnostic performance for spinal tuberculosis.
[1] 毕圣贤,胡锡健,张辉国.2014年至2018年中国大陆结核病疫情的时空分布特征及其影响因素[J].中华疾病控制杂志,2021,25(7):758-762,811. [2] Pattnaik B,Patnaik N,Mittal S,et al.Micro RNAs as potential biomarkers in tuberculosis:a systematic review[J].Noncoding RNA Res,2022,7(1):16-26. [3] Kaushik A C,Wu Q,Lin L,et al.Exosomal ncRNAs profiling of mycobacterial infection identified miRNA-185-5p as a novel biomarker for tuberculosis[J].Brief Bioinform,2021,22(6):210. [4] Riviere E,Heupink T H,Ismail N,et al.Capacity building for whole genome sequencing of mycobacterium tuberculosis and bioinformatics in high TB burden countries[J].Brief Bioinform,2021,22(4):246. [5] Wang L,Xiong Y,Fu B,et al.MicroRNAs as immune regulators and biomarkers in tuberculosis[J].Front Immunol,2022,13:1027472. [6] Alipoor S D,Tabarsi P,Varahram M,et al.Serum exosomal miRNAs are associated with active pulmonary tuberculosis[J].Dis Markers,2019,2019:1907426. [7] Gong Z,Han S,Liang T,et al.Mycobacterium tuberculosis effector PPE36 attenuates host cytokine storm damage via inhibiting macrophage M1 polarization[J].Cell Physiol,2021,236(11):7405-7420. [8] Yang C,Shi Z,Hu J,et al.miRNA-155 expression and role in pathogenesis in spinal tuberculosis-induced intervertebral disc destruction[J].Exp Ther Med,2019,17(4):3239-3246. [9] Shams R,Saberi S,Zali M,et al.Identification of potential microRNA panels for pancreatic cancer diagnosis using microarray datasets and bioinformatics methods[J].Sci Rep,2020,10(1):7559. [10] Zhang H,Wang J,Wang Y,et al.Long non-coding LEF1-AS1 sponge miR-5100 regulates apoptosis and autophagy in gastric cancer cells via the miR-5100/DEK/AMPK-mTOR axis[J].Int Mol Sci,2022,23(9):4787. [11] Kisan A,Chhabra R.Modulation of gene expression by YTH domain family (YTHDF) proteins in human physiology and pathology[J].Cell Physiol,2023,238(1):5-31.