In silico analysis of serum miRNA profiles in seronegative and seropositive rheumatoid arthritis patients by small RNA sequencing

Abstract

Rheumatoid arthritis (RA) is a refractory autoimmune disease, affecting about 1% of the world’s population. RA is divided into seronegative RA and seropositive RA. However, biomarkers for discriminating between seronegative and seropositive RA have not been reported. In this study, we profiled serum miRNAs in seronegative RA patients (N-RA), seropositive RA patients (P-RA) and healthy controls (HC) by small RNA sequencing. Results indicated that compared with HC group, there were one up-regulated and four downregulated miRNAs in N-RA group (fold change ≥ 2 and P value < 0.05); compared with P-RA group, there were two up-regulated and four downregulated miRNAs in N-RA group; compared with HC group, there were three up-regulated and four downregulated miRNAs in P-RA group. Among them, the level of hsa-miR-362-5p in N-RA group was up-regulated compared with that in HC group and P-RA group, and the level of hsa-miR-6855-5p and hsa-miR-187-3p in P-RA group was upregulated compared with that in N-RA group and HC group. Validation by qPCR confirmed that serum hsa-miR-362-5p level was elevated in N-RA group. Subsequently, by analyzing the target genes using RNAhybrid, PITA, Miranda and TargetScan and functions of differential miRNAs utilizing Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG), we found that the target genes and molecular pathways regulated by miRNAs in seronegative RA and seropositive RA were roughly the same, and miRNAs in these two diseases may participate in the occurrence and development of diseases by regulating the immune system. In conclusion, this study revealed the profiles of serum miRNAs in seronegative and seropositive RA patients for the first time, providing potential biomarkers and targets for the diagnosis and treatment of seronegative and seropositive RA.