The Endogenous Plasma Small RNAome of Rheumatoid Arthritis

The intricacies of isomiRs: from classification to clinical relevance

MicroRNAs (miRNAs) and isoforms of their archetype, called isomiRs, regulate gene expression via complementary base-pair binding to messenger RNAs (mRNAs). The partially evolutionarily conserved isomiR sequence variations are differentially expressed among tissues, populations, and genders, and between healthy and diseased states. Aiming towards the clinical use of isomiRs as diagnostic biomarkers and for therapeutic purposes, several challenges need to be addressed, including (i) clarification of isomiR definition, (ii) improved annotation in databases with new standardization (such as the mirGFF3 format), and (iii) improved methods of isomiR detection, functional verification, and in silico analysis. In this review we discuss the respective challenges, and highlight the opportunities for clinical use of isomiRs, especially in the light of increasing amounts of next-generation sequencing (NGS) data.

Anti-Inflammatory Effect of High-Density Lipoprotein Blunted by Delivery of Altered MicroRNA Cargo in Patients With Rheumatoid Arthritis

OBJECTIVE

High-density lipoprotein (HDL) has well-characterized anti-atherogenic cholesterol efflux and antioxidant functions. Another function of HDL uncharacterized in rheumatoid arthritis (RA) is its ability to transport microRNAs (miRNAs) between cells and thus alter cellular function. The study's purpose was to determine if HDL-miRNA cargo is altered and affects inflammation in RA.

METHODS

HDL-microRNAs were characterized in 30 RA and 30 control participants by next generation sequencing and quantitative polymerase chain reaction. The most abundant differentially expressed miRNA was evaluated further. The function of miR-1246 was assessed by miRNA mimics, antagomiRs, small interfering RNA knockdown, and luciferase assays. Monocyte-derived macrophages were treated with miR-1246-loaded HDL and unmodified HDL from RA and control participants to measure delivery of miR-1246 and its effect on interleukin-6 (IL-6).

RESULTS

The most abundant miRNA on HDL was miR-1246; it was significantly enriched two-fold on HDL from RA versus control participants. HDL-mediated miR-1246 delivery to macrophages significantly increased IL6 expression 43-fold. miR-1246 delivery significantly decreased DUSP3 1.5-fold and DUSP3 small interfering RNA knockdown increased macrophage IL6 expression. Luciferase assay indicated DUSP3 is a direct target of miR-1246. Unmodified HDL from RA delivered 1.6-fold more miR-1246 versus control participant HDL. Unmodified HDL from both RA and control participants attenuated activated macrophage IL6 expression, but this effect was significantly blunted in RA so that IL6 expression was 3.4-fold higher after RA versus control HDL treatment.

CONCLUSION

HDL-miR-1246 was increased in RA versus control participants and delivery of miR-1246 to macrophages increased IL-6 expression by targeting DUSP3. The altered HDL-miRNA cargo in RA blunted HDL's anti-inflammatory effect.

Multi-omics analysis reveals COVID-19 vaccine induced attenuation of inflammatory responses during breakthrough disease

The immune mechanisms mediating COVID-19 vaccine attenuation of COVID-19 remain undescribed. We conducted comprehensive analyses detailing immune responses to SARS-CoV-2 virus in blood post-vaccination with ChAdOx1 nCoV-19 or a placebo. Samples from randomised placebo-controlled trials (NCT04324606 and NCT04400838) were taken at baseline, onset of COVID-19-like symptoms, and 7 days later, confirming COVID-19 using nucleic amplification test (NAAT test) via real-time PCR (RT-PCR). Serum cytokines were measured with multiplexed immunoassays. The transcriptome was analysed with long, short and small RNA sequencing. We found attenuation of RNA inflammatory signatures in ChAdOx1 nCoV-19 compared with placebo vaccinees and reduced levels of serum proteins associated with COVID-19 severity. KREMEN1, a putative alternative SARS-CoV-2 receptor, was downregulated in placebo compared with ChAdOx1 nCoV-19 vaccinees. Vaccination ameliorates reductions in cell counts across leukocyte populations and platelets noted at COVID-19 onset, without inducing potentially deleterious Th2-skewed immune responses. Multi-omics integration links a global reduction in miRNA expression at COVID-19 onset to increased pro-inflammatory responses at the mRNA level. This study reveals insights into the role of COVID-19 vaccines in mitigating disease severity by abrogating pro-inflammatory responses associated with severe COVID-19, affirming vaccine-mediated benefit in breakthrough infection, and highlighting the importance of clinically relevant endpoints in vaccine evaluation.

LncRNA expression profiling in exosomes derived from synovial fluid of patients with rheumatoid arthritis

OBJECTIVE

To identify the long non-coding RNA (lncRNA) expression profiling in exosomes derived from synovial fluid of rheumatoid arthritis (RA) patients, and carry out bioinformatics analysis on target genes of differentially expressed lncRNAs.

METHODS

Exosomes were isolated from synovial fluid via ultracentrifugation. RNAs were extracted from exosomes by using HiPure Liquid RNA/miRNA kits, followed by lncRNA sequencing. Differentially expressed lncRNAs in RA were screened, and bioinformatics analysis of their target genes was carried out. qRT-PCR was used to verify the lncRNA expression levels.

RESULTS

Compared with osteoarthritis (OA), 347 lncRNAs were found differentially expressed in RA. Compared with gout, 805 lncRNAs were found differentially expressed in RA. Compared with both OA and gout, 85 lncRNAs were found specially expressed in RA (65 were upregulated (including ENST00000433825.1)). Functional analysis of target genes of the specially expressed lncRNAs revealed significant enrichment of "autophagy" and "mTOR signaling pathway". The qRT-PCR results indicated that ENST00000433825.1 was highly expressed in RA, compared with both OA and gout (P < 0.05), which matched the lncRNA sequencing results. Correlation analysis showed that the level of ENST00000433825.1 in RA patients was significantly and positively correlated with the level of C-reactive protein (CRP) (P < 0.001).

CONCLUSIONS

The lncRNA expression profiling in exosomes derived from synovial fluid of RA was significantly different from OA and gout. ENST00000433825.1 was highly and uniquely expressed in RA and significantly and positively correlated with CRP, which might provide a diagnostic and therapeutic biomarker for RA.

Emerging roles of SnoRNAs in the pathogenesis and treatment of autoimmune disorders

SnoRNAs (small non-coding RNAs) have recently gained prominence in autoimmune diseases, revealing their crucial role in modulating the immune response and contributing to disease pathogenesis. Initially known for their involvement in ribosomal RNA processing and modification, molecular biology and genomics advancements have uncovered their broader impact on cellular function, especially in autoimmune disorders. Autoimmune diseases represent conditions characterized by the immune system's erroneous attacks on self-tissues, encompassing disorders like systemic lupus erythematosus, rheumatoid arthritis, and multiple sclerosis. The complex etiology of these conditions involves a delicate interplay of genetic and environmental factors. Emerging evidence suggests that snoRNAs initially recognized for their housekeeping roles, extend their influence on immune regulation through diverse mechanisms. SnoRNAs have been implicated in epigenetic modification, directly affecting the gene expression profiles of immune cells. Their ability to guide site-specific changes on ribosomal RNAs and other non-coding RNAs can significantly influence the translation of proteins involved in immune response pathways. Moreover, snoRNAs interact with key immune-related proteins, modulating their functions and subsequently impacting immune cell development, activation, and tolerance. Dysregulation of snoRNA expression has been observed in various autoimmune diseases, underscoring their potential as biomarkers for disease diagnosis, prognosis, and therapeutic targets. Manipulating snoRNA expression or activity is a promising therapeutic intervention avenue, offering the potential for personalized treatment strategies in autoimmune diseases. However, there remains a need for comprehensive research efforts to elucidate the precise molecular mechanisms underlying snoRNA-mediated immune modulation. Further investigations in this domain are essential to unravel the potential of snoRNAs in autoimmune disorders.

Interrelated grid of non-coding RNA: An important aspect in Rheumatoid Arthritis pathogenesis

Inflammation and autoimmunity are the root cause of rheumatoid arthritis, a destructive disease of joints. Multiple biomolecules are involved in the pathogenesis of RA and are related to various events of molecular biology. RNA is a versatile biomolecule, playing numerous roles at structural, functional, and regulatory stages to maintain cellular homeostasis. The involvement of RNA (coding/non-coding) in disease development and progression has left a wide whole to fill with newer approaches. Non-coding RNAs belong to the housekeeping and regulatory categories and both have their specific roles, and their alteration causes specific implications in disease pathogenesis. Housekeeping RNAs, rRNA, tRNA and regulatory RNA, micro-RNA, circular RNA, piRNA and long non-coding RNA were found to be important regulators of inflammation. They work at the pre-and post-transcriptional levels and were found to be more intriguing to study their regulatory impact on disease pathogenesis. The review addresses a question on how the non-coding RNA gets involved in early RA pathogenesis and can be utilized to know their targets to understand the disease better and make way towards the unresolved mystery of RA development.

iSnoDi-LSGT: identifying snoRNA-disease associations based on local similarity constraints and global topological constraints

Growing evidence proves that small nucleolar RNAs (snoRNAs) have important functions in various biological processes, the malfunction of which leads to the emergence and development of complex diseases. However, identifying snoRNA-disease associations is an ongoing challenging task due to the considerable time- and money-consuming biological experiments. Therefore, it is urgent to design efficient and economical methods for the identification of snoRNA-disease associations. In this regard, we propose a computational method named iSnoDi-LSGT, which utilizes snoRNA sequence similarity and disease similarity as local similarity constraints. The iSnoDi-LSGT predictor further employs network embedding technology to extract topological features of snoRNAs and diseases, based on which snoRNA topological similarity and disease topological similarity are calculated as global topological constraints. To the best of our knowledge, the iSnoDi-LSGT is the first computational method for snoRNA-disease association identification. The experimental results indicate that the iSnoDi-LSGT predictor can effectively predict unknown snoRNA-disease associations. The web server of the iSnoDi-LSGT predictor is freely available at http://bliulab.net/iSnoDi-LSGT.

Integration of small RNAs from plasma and cerebrospinal fluid for classification of multiple sclerosis

Multiple Sclerosis (MS) is an autoimmune, neurological disease, commonly presenting with a relapsing-remitting form, that later converts to a secondary progressive stage, referred to as RRMS and SPMS, respectively. Early treatment slows disease progression, hence, accurate and early diagnosis is crucial. Recent advances in large-scale data processing and analysis have progressed molecular biomarker development. Here, we focus on small RNA data derived from cell-free cerebrospinal fluid (CSF), cerebrospinal fluid cells, plasma and peripheral blood mononuclear cells as well as CSF cell methylome data, from people with RRMS (n = 20), clinically/radiologically isolated syndrome (CIS/RIS, n = 2) and neurological disease controls (n = 14). We applied multiple co-inertia analysis (MCIA), an unsupervised and thereby unbiased, multivariate method for simultaneous data integration and found that the top latent variable classifies RRMS status with an Area Under the Receiver Operating Characteristics (AUROC) score of 0.82. Variable selection based on Lasso regression reduced features to 44, derived from the small RNAs from plasma (20), CSF cells (8) and cell-free CSF (16), with a marginal reduction in AUROC to 0.79. Samples from SPMS patients (n = 6) were subsequently projected on the latent space and differed significantly from RRMS and controls. On contrary, we found no differences between relapse and remission or between inflammatory and non-inflammatory disease controls, suggesting that the latent variable is not prone to inflammatory signals alone, but could be MS-specific. Hence, we here showcase that integration of small RNAs from plasma and CSF can be utilized to distinguish RRMS from SPMS and neurological disease controls.

Changes in the Small Noncoding RNAome During M1 and M2 Macrophage Polarization

Macrophages belong to a special phagocytic subgroup of human leukocytes and are one of the important cells of the human immune system. Small noncoding RNAs are a group of small RNA molecules that can be transcribed without the ability to encode proteins but could play a specific function in cells. SncRNAs mainly include microRNAs (miRNAs) and piwi-interacting RNAs (piRNAs), small nucleolar RNAs (snoRNAs), small nuclear RNAs (snRNAs) and repeat RNAs. We used high-throughput sequencing analysis and qPCR to detect the expression changes of the small noncoding RNAome during macrophage polarization. Our results showed that 84 miRNAs and 47 miRNAs with were downregulated during M1 macrophage polarization and that 11 miRNAs were upregulated and 19 miRNAs were downregulated during M2 macrophage polarization. MiR-novel-3-nature and miR-27b-5p could promote expression of TNF-α which was marker gene of M1 macrophages. The piRNA analysis results showed that 69 piRNAs were upregulated and 61 piRNAs were downregulated during M1 macrophage polarization and that 3 piRNAs were upregulated and 10 piRNAs were downregulated during M2 macrophage polarization. DQ551351 and DQ551308 could promote the mRNA expression of TNF-α and DQ551351overexpression promoted the antitumor activity of M1 macrophages. SnoRNA results showed that 62 snoRNAs were upregulated and 59 snoRNAs were downregulated during M1 macrophage polarization, whereas 6 snoRNAs were upregulated and 10 snoRNAs were downregulated during M2 macrophage polarization. Overexpression of snoRNA ENSMUST00000158683.2 could inhibit expression of TNF-α. For snRNA, we found that 12 snRNAs were upregulated and 15 snRNAs were downregulated during M1 macrophage polarization and that 2 snRNAs were upregulated during M2 macrophage polarization. ENSMUSG00000096786 could promote expression of IL-1 and iNOS and ENSMUSG00000096786 overexpression promoted the antitumor activity of M1 macrophages. Analysis of repeat RNAs showed that 7 repeat RNAs were upregulated and 9 repeat RNAs were downregulated during M1 macrophage polarization and that 2 repeat RNAs were downregulated during M2 macrophage polarization. We first reported the expression changes of piRNA, snoRNA, snRNA and repeat RNA during macrophage polarization, and preliminarily confirmed that piRNA, snoRNA and snRNA can regulate the function of macrophages.

MicroRNA-Mediated Epigenetic Regulation of Rheumatoid Arthritis Susceptibility and Pathogenesis

MicroRNAs (miRNAs) play crucial roles in regulating the transcriptome and development of rheumatoid arthritis (RA). Currently, a comprehensive map illustrating how miRNAs regulate transcripts, pathways, immune system differentiation, and their interactions with terminal cells such as fibroblast-like synoviocytes (FLS), immune-cells, osteoblasts, and osteoclasts are still laking. In this review, we summarize the roles of miRNAs in the susceptibility, pathogenesis, diagnosis, therapeutic intervention, and prognosis of RA. Numerous miRNAs are abnormally expressed in cells involved in RA and regulate target genes and pathways, including NF-κB, Fas-FasL, JAK-STAT, and mTOR pathways. We outline how functional genetic variants of miR-499 and miR-146a partly explain susceptibility to RA. By regulating gene expression, miRNAs affect T cell differentiation into diverse cell types, including Th17 and Treg cells, thus constituting promising gene therapy targets to modulate the immune system in RA. We summarize the diagnostic and prognostic potential of blood-circulating and cell-free miRNAs, highlighting the opportunity to combine these miRNAs with antibodies to cyclic citrullinated peptide (ACCP) to allow accurate diagnosis and prognosis, particularly for seronegative patients. Furthermore, we review the evidence implicating miRNAs as promising biomarkers of efficiency and response of, and resistance to, disease-modifying anti-rheumatic drugs and immunotherapy. Finally, we discuss the autotherapeutic effect of miRNA intervention as a step toward the development of miRNA-based anti-RA drugs. Collectively, the current evidence supports miRNAs as interesting targets to better understand the pathogenetic mechanisms of RA and design more efficient therapeutic interventions.

The non-coding RNA interactome in joint health and disease

Non-coding RNAs have distinct regulatory roles in the pathogenesis of joint diseases including osteoarthritis (OA) and rheumatoid arthritis (RA). As the amount of high-throughput profiling studies and mechanistic investigations of microRNAs, long non-coding RNAs and circular RNAs in joint tissues and biofluids has increased, data have emerged that suggest complex interactions among non-coding RNAs that are often overlooked as critical regulators of gene expression. Identifying these non-coding RNAs and their interactions is useful for understanding both joint health and disease. Non-coding RNAs regulate signalling pathways and biological processes that are important for normal joint development but, when dysregulated, can contribute to disease. The specific expression profiles of non-coding RNAs in various disease states support their roles as promising candidate biomarkers, mediators of pathogenic mechanisms and potential therapeutic targets. This Review synthesizes literature published in the past 2 years on the role of non-coding RNAs in OA and RA with a focus on inflammation, cell death, cell proliferation and extracellular matrix dysregulation. Research to date makes it apparent that 'non-coding' does not mean 'non-essential' and that non-coding RNAs are important parts of a complex interactome that underlies OA and RA.

Possible Roles of tRNA Fragments, as New Regulatory ncRNAs, in the Pathogenesis of Rheumatoid Arthritis

Understanding the pathophysiology of rheumatoid arthritis (RA) has led to the successful development of molecule-targeted drugs for the treatment of RA. However, some RA patients are refractory to these treatments, suggesting that the pathological mechanism of the disease is not entirely understood. Genome and transcriptome analysis is essential for understanding the unknown pathophysiology of human diseases. Rapid and more comprehensive gene analysis technologies have revealed notable changes in the expression of coding RNA and non-coding RNA in RA patients. This review focuses on the current state of non-coding RNA research in relation to RA, especially on tRNA fragments. Interestingly, it has been found that tRNA fragments repress translation and are antiapoptotic. The association between tRNA fragments and various diseases has been studied, and this article reviews the possible role of tRNA fragments in RA.

Plasma miRNAs improve the prediction of coronary atherosclerosis in patients with rheumatoid arthritis

OBJECTIVE

MicroRNAs (miRNAs) regulate gene expression and are disease biomarkers. Rheumatoid arthritis (RA) patients have accelerated atherosclerosis leading to excess cardiovascular morbidity and mortality, but traditional risk factors for cardiovascular risk stratification are inadequate. In the general population, miRNAs improve cardiovascular risk estimation beyond traditional risk factors. Our objective was to develop a miRNA panel that predicts coronary atherosclerosis in RA patients.

METHODS

Plasma small RNA next-generation sequencing (NGS) was performed on 161 RA patients whose Agatston scores for coronary artery calcium were previously measured. Random forest analysis of plasma NGS miRNA expression was used to determine which miRNAs best differentiated between those patients with and without coronary artery calcium. Top predictive miRNAs were assayed by quantitative PCR (qPCR). Elastic net regression was used to develop the most parsimonious models with qPCR-measured miRNA concentrations and clinical variables (age, sex, ACC/AHA 10-year risk score, DAS28 score, and diabetes) separately to predict the presence of coronary artery calcium and high coronary artery calcium. C-statistics were used to assess performance model performance.

RESULTS

The top miRNAs which differentiated those with and without coronary atherosclerosis based on random forest analysis included let-7c-5p, miR-30e-5p, miR-30c-5p, miR-4446-3p, miR-126-5p, miR-3168, miR-425-5p, miR-126-3p, miR-30a-5p, and miR-125a-5p. For coronary artery calcium prediction, addition of all miRNAs except miR-126-3p to clinical factors improved the c-statistic modestly from 0.86 to 0.87. For high coronary artery calcium prediction, addition of all miRNAs except miR-30c-5p to clinical factors improved the c-statistic from 0.75 to 0.80.

CONCLUSION

A plasma miRNA panel improved the prediction of high coronary artery calcium beyond traditional risk factors and RA disease activity. Further evaluation of the miRNA panel for prediction of coronary events in RA is necessary. Key Point • A plasma microRNA panel including let-7c-5p, miR-30a-5p, miR-30e-5p, miR-125a-5p, miR-126-3p, miR-126-5p, miR-425-5p, miR-3168, and miR-4446-3p improved the prediction of high coronary artery calcium beyond clinical factors in patients with rheumatoid arthritis.

The Expression of Non-Coding RNAs and Their Target Molecules in Rheumatoid Arthritis: A Molecular Basis for Rheumatoid Pathogenesis and Its Potential Clinical Applications

Rheumatoid arthritis (RA) is a typical autoimmune-mediated rheumatic disease presenting as a chronic synovitis in the joint. The chronic synovial inflammation is characterized by hyper-vascularity and extravasation of various immune-related cells to form lymphoid aggregates where an intimate cross-talk among innate and adaptive immune cells takes place. These interactions facilitate production of abundant proinflammatory cytokines, chemokines and growth factors for the proliferation/maturation/differentiation of B lymphocytes to become plasma cells. Finally, the autoantibodies against denatured immunoglobulin G (rheumatoid factors), EB virus nuclear antigens (EBNAs) and citrullinated protein (ACPAs) are produced to trigger the development of RA. Furthermore, it is documented that gene mutations, abnormal epigenetic regulation of peptidylarginine deiminase genes 2 and 4 (PADI2 and PADI4), and thereby the induced autoantibodies against PAD2 and PAD4 are implicated in ACPA production in RA patients. The aberrant expressions of non-coding RNAs (ncRNAs) including microRNAs (miRs) and long non-coding RNAs (lncRNAs) in the immune system undoubtedly derange the mRNA expressions of cytokines/chemokines/growth factors. In the present review, we will discuss in detail the expression of these ncRNAs and their target molecules participating in developing RA, and the potential biomarkers for the disease, its diagnosis, cardiovascular complications and therapeutic response. Finally, we propose some prospective investigations for unraveling the conundrums of rheumatoid pathogenesis.

Circulating microbial small RNAs are altered in patients with rheumatoid arthritis

OBJECTIVES

To determine if plasma microbial small RNAs (sRNAs) are altered in patients with rheumatoid arthritis (RA) compared with control subjects, associated with RA disease-related features, and altered by disease-modifying antirheumatic drugs (DMARDs).

METHODS

sRNA sequencing was performed on plasma from 165 patients with RA and 90 matched controls and a separate cohort of 70 patients with RA before and after starting a DMARD. Genome alignments for RA-associated bacteria, representative bacterial and fungal human microbiome genomes and environmental bacteria were performed. Microbial genome counts and individual sRNAs were compared across groups and correlated with disease features. False discovery rate was set at 0.05.

RESULTS

Genome counts of Lactobacillus salivarius, Anaerobaculum hydrogeniformans, Staphylococcus epidermidis, Staphylococcus aureus, Paenisporosarcina spp, Facklamia hominis, Sphingobacterium spiritivorum, Lentibacillus amyloliquefaciens, Geobacillus spp, and Pseudomonas fluorescens were significantly decreased in the plasma of RA compared with control subjects. Three microbial transfer RNA-derived sRNAs were increased in RA versus controls and inversely associated with disease activity. Higher total microbial sRNA reads were associated with lower disease activity in RA. Baseline total microbial sRNAs were threefold higher among patients who improved with DMARD versus those who did not but did not change significantly after 6 months of treatment.

CONCLUSION

Plasma microbial sRNA composition is altered in RA versus control subjects and associated with some measures of RA disease activity. DMARD treatment does not alter microbial sRNA abundance or composition, but increased abundance of microbial sRNAs at baseline was associated with disease activity improvement at 6 months.

Small Non-Coding RNAome of Ageing Chondrocytes

Ageing is a leading risk factor predisposing cartilage to osteoarthritis. However, little research has been conducted on the effect of ageing on the expression of small non-coding RNAs (sncRNAs). RNA from young and old chondrocytes from macroscopically normal equine metacarpophalangeal joints was extracted and subjected to small RNA sequencing (RNA-seq). Differential expression analysis was performed in R using package DESeq2. For transfer RNA (tRNA) fragment analysis, tRNA reads were aligned to horse tRNA sequences using Bowtie2 version 2.2.5. Selected microRNA (miRNAs or miRs) and small nucleolar RNA (snoRNA) findings were validated using real-time quantitative Polymerase Chain Reaction (qRT-PCR) in an extended cohort of equine chondrocytes. tRNA fragments were further investigated in low- and high-grade OA human cartilage tissue. In total, 83 sncRNAs were differentially expressed between young and old equine chondrocytes, including miRNAs, snoRNAs, small nuclear RNAs (snRNAs), and tRNAs. qRT-PCR analysis confirmed findings. tRNA fragment analysis revealed that tRNA halves (tiRNAs), tiRNA-5035-GluCTC and tiRNA-5031-GluCTC-1 were reduced in both high grade OA human cartilage and old equine chondrocytes. For the first time, we have measured the effect of ageing on the expression of sncRNAs in equine chondrocytes. Changes were detected in a number of different sncRNA species. This study supports a role for sncRNAs in ageing cartilage and their potential involvement in age-related cartilage diseases.