Search

Abstract Background Synovial fibroblasts (SFs) are specialized cells of the synovium that provide nutrients and lubricants for the proper function of diarthrodial joints. Recent evidence appreciates the contribution of SF heterogeneity in arthritic pathologies. However, the normal SF profiles and the molecular networks that govern the transition from homeostatic to arthritic SF heterogeneity remain poorly defined. Methods We applied a combined analysis of single-cell (sc) transcriptomes and epigenomes (scRNA-seq and scATAC-seq) to SFs derived from naïve and hTNFtg mice (mice that overexpress human TNF, a murine model for rheumatoid arthritis), by employing the Seurat and ArchR packages. To identify the cellular differentiation lineages, we conducted velocity and trajectory analysis by combining state-of-the-art algorithms including scVelo, Slingshot, and PAGA. We integrated the transcriptomic and epigenomic data to infer gene regulatory networks using ArchR and custom-implemented algorithms. We performed a canonical correlation analysis-based integration of murine data with publicly available datasets from SFs of rheumatoid arthritis patients and sought to identify conserved gene regulatory networks by utilizing the SCENIC algorithm in the human arthritic scRNA-seq atlas. Results By comparing SFs from healthy and hTNFtg mice, we revealed seven homeostatic and two disease-specific subsets of SFs. In healthy synovium, SFs function towards chondro- and osteogenesis, tissue repair, and immune surveillance. The development of arthritis leads to shrinkage of homeostatic SFs and favors the emergence of SF profiles marked by Dkk3 and Lrrc15 expression, functioning towards enhanced inflammatory responses and matrix catabolic processes. Lineage inference analysis indicated that specific Thy1+ SFs at the root of trajectories lead to the intermediate Thy1+/Dkk3+/Lrrc15+ SF states and culminate in a destructive and inflammatory Thy1− SF identity. We further uncovered epigenetically primed gene programs driving the expansion of these arthritic SFs, regulated by NFkB and new candidates, such as Runx1. Cross-species analysis of human/mouse arthritic SF data determined conserved regulatory and transcriptional networks. Conclusions We revealed a dynamic SF landscape from health to arthritis providing a functional genomic blueprint to understand the joint pathophysiology and highlight the fibroblast-oriented therapeutic targets for combating chronic inflammatory and destructive arthritic disease.

Abstract Background The transmembrane-TNF transgenic mouse, TgA86, has been shown to develop spontaneously peripheral arthritis with signs of axial involvement. To assess similarity to human spondyloarthritis, we performed detailed characterization of the axial, peripheral, and comorbid pathologies of this model. Methods TgA86 bone pathologies were assessed at different ages using CT imaging of the spine, tail vertebrae, and hind limbs and characterized in detail by histopathological and immunohistochemical analysis. Cardiac function was examined by echocardiography and electrocardiography and bone structural parameters by μCT analysis. The response of TgA86 mice to either early or late anti-TNF treatment was evaluated clinically, histopathologically, and by μCT analysis. Results TgA86 mice developed with 100% penetrance spontaneous axial and peripheral pathology which progressed with time and manifested as reduced body weight and body length, kyphosis, tail bendings, as well as swollen and distorted hind joints. Whole-body CT analysis at advanced ages revealed bone erosions of sacral and caudal vertebrae as well as of sacroiliac joints and hind limbs and, also, new ectopic bone formation and eventually vertebral fusion. The pathology of these mice highly resembled that of SpA patients, as it evolved through an early inflammatory phase, evident as enthesitis and synovitis in the affected joints, characterized by mesenchymal cell accumulation, and neutrophilic infiltration. Subsequently, regression of inflammation was accompanied by ectopic bone formation, leading to ankylosis. In addition, both systemic bone loss and comorbid heart valve pathology were evident. Importantly, early anti-TNF treatment, similar to clinical treatment protocols, significantly reduced the inflammatory phase of both the axial and peripheral pathology of TgA86 mice. Conclusions The TgA86 mice develop a spontaneous peripheral and axial biphasic pathology accompanied by comorbid heart valvular dysfunction and osteoporosis, overall reproducing the progression of pathognomonic features of human spondyloarthritis. Therefore, the TgA86 mouse represents a valuable model for deciphering the role of transmembrane TNF in the pathogenic mechanisms of spondyloarthritis and for assessing the efficacy of human therapeutics targeting different phases of the disease.