ABSTRACT
Rheumatoid arthritis (RA) is a complex immune-mediated inflammatory disorder in which patients suffer from inflammatory-erosive arthritis. Recent advances on histopathology heterogeneity of RA synovial tissue revealed three distinct phenotypes based on cellular composition (pauci-immune, diffuse and lymphoid), suggesting that distinct etiologies warrant specific targeted therapy which motivates a need for cost effective phenotyping tools in preclinical and clinical settings. To this end, we developed an automated multi-scale computational pathotyping (AMSCP) pipeline for both human and mouse synovial tissue with two distinct components that can be leveraged together or independently: (1) segmentation of different tissue types to characterize tissue-level changes, and (2) cell type classification within each tissue compartment that assesses change across disease states. Here, we demonstrate the efficacy, efficiency, and robustness of the AMSCP pipeline as well as the ability to discover novel phenotypes. Taken together, we find AMSCP to be a valuable cost-effective method for both pre-clinical and clinical research.
Subject(s)
Arthritis, Rheumatoid , Synovial Membrane , Humans , Synovial Membrane/pathology , Synovial Membrane/immunology , Animals , Arthritis, Rheumatoid/pathology , Arthritis, Rheumatoid/immunology , Mice , Phenotype , Computational Biology/methods , Inflammation/pathologyABSTRACT
It has been presumed that rheumatoid arthritis (RA) joint pain is related to inflammation in the synovium; however, recent studies reveal that pain scores in patients do not correlate with synovial inflammation. We developed a machine-learning approach (graph-based gene expression module identification or GbGMI) to identify an 815-gene expression module associated with pain in synovial biopsy samples from patients with established RA who had limited synovial inflammation at arthroplasty. We then validated this finding in an independent cohort of synovial biopsy samples from patients who had early untreated RA with little inflammation. Single-cell RNA sequencing analyses indicated that most of these 815 genes were most robustly expressed by lining layer synovial fibroblasts. Receptor-ligand interaction analysis predicted cross-talk between human lining layer fibroblasts and human dorsal root ganglion neurons expressing calcitonin gene-related peptide (CGRP+). Both RA synovial fibroblast culture supernatant and netrin-4, which is abundantly expressed by lining fibroblasts and was within the GbGMI-identified pain-associated gene module, increased the branching of pain-sensitive murine CGRP+ dorsal root ganglion neurons in vitro. Imaging of solvent-cleared synovial tissue with little inflammation from humans with RA revealed CGRP+ pain-sensing neurons encasing blood vessels growing into synovial hypertrophic papilla. Together, these findings support a model whereby synovial lining fibroblasts express genes associated with pain that enhance the growth of pain-sensing neurons into regions of synovial hypertrophy in RA.
Subject(s)
Arthritis, Rheumatoid , Calcitonin Gene-Related Peptide , Humans , Mice , Animals , Calcitonin Gene-Related Peptide/genetics , Calcitonin Gene-Related Peptide/metabolism , Arthritis, Rheumatoid/complications , Arthritis, Rheumatoid/genetics , Arthritis, Rheumatoid/metabolism , Synovial Membrane/pathology , Inflammation/pathology , Fibroblasts/pathology , Pain/metabolism , Gene Expression , Cells, CulturedABSTRACT
Rheumatoid arthritis is a prototypical autoimmune disease that causes joint inflammation and destruction1. There is currently no cure for rheumatoid arthritis, and the effectiveness of treatments varies across patients, suggesting an undefined pathogenic diversity1,2. Here, to deconstruct the cell states and pathways that characterize this pathogenic heterogeneity, we profiled the full spectrum of cells in inflamed synovium from patients with rheumatoid arthritis. We used multi-modal single-cell RNA-sequencing and surface protein data coupled with histology of synovial tissue from 79 donors to build single-cell atlas of rheumatoid arthritis synovial tissue that includes more than 314,000 cells. We stratified tissues into six groups, referred to as cell-type abundance phenotypes (CTAPs), each characterized by selectively enriched cell states. These CTAPs demonstrate the diversity of synovial inflammation in rheumatoid arthritis, ranging from samples enriched for T and B cells to those largely lacking lymphocytes. Disease-relevant cell states, cytokines, risk genes, histology and serology metrics are associated with particular CTAPs. CTAPs are dynamic and can predict treatment response, highlighting the clinical utility of classifying rheumatoid arthritis synovial phenotypes. This comprehensive atlas and molecular, tissue-based stratification of rheumatoid arthritis synovial tissue reveal new insights into rheumatoid arthritis pathology and heterogeneity that could inform novel targeted treatments.
Subject(s)
Arthritis, Rheumatoid , Humans , Arthritis, Rheumatoid/complications , Arthritis, Rheumatoid/genetics , Arthritis, Rheumatoid/immunology , Arthritis, Rheumatoid/pathology , Cytokines/metabolism , Inflammation/complications , Inflammation/genetics , Inflammation/immunology , Inflammation/pathology , Synovial Membrane/pathology , T-Lymphocytes/immunology , B-Lymphocytes/immunology , Genetic Predisposition to Disease/genetics , Phenotype , Single-Cell Gene Expression AnalysisABSTRACT
It has been presumed that rheumatoid arthritis (RA) joint pain is related to inflammation in the synovium; however, recent studies reveal that pain scores in patients do not correlate with synovial inflammation. We identified a module of 815 genes associated with pain, using a novel machine learning approach, Graph-based Gene expression Module Identification (GbGMI), in samples from patients with longstanding RA, but limited synovial inflammation at arthroplasty, and validated this finding in an independent cohort of synovial biopsy samples from early, untreated RA patients. Single-cell RNA-seq analyses indicated these genes were most robustly expressed by lining layer fibroblasts and receptor-ligand interaction analysis predicted robust lining layer fibroblast crosstalk with pain sensitive CGRP+ dorsal root ganglion sensory neurons. Netrin-4, which is abundantly expressed by lining fibroblasts and associated with pain, significantly increased the branching of pain-sensitive CGRP+ neurons in vitro . We conclude GbGMI is a useful method for identifying a module of genes that associate with a clinical feature of interest. Using this approach, we find that Netrin-4 is produced by synovial fibroblasts in the absence of inflammation and can enhance the outgrowth of CGRP+ pain sensitive nerve fibers. One Sentence Summary: Machine Learning reveals synovial fibroblast genes related to pain affect sensory nerve growth in Rheumatoid Arthritis addresses unmet clinical need.
ABSTRACT
OBJECTIVE: We sought to develop computer vision methods to quantify aggregates of cells in synovial tissue and compare these with clinical and gene expression parameters. METHODS: We assembled a computer vision pipeline to quantify five features encompassing synovial cell density and aggregates and compared these with pathologist scores, disease classification, autoantibody status, and RNA expression in a cohort of 156 patients with rheumatoid arthritis (RA) and 149 patients with osteoarthritis (OA). RESULTS: All five features were associated with pathologist scores of synovial lymphocytic inflammation (P < 0.0001). Three features that related to the cells per unit of tissue were significantly increased in patients with both seronegative and seropositive RA compared with those with OA; on the other hand, aggregate features (number and diameter) were significantly increased in seropositive, but not seronegative, RA compared with OA. Aggregate diameter was associated with the gene expression of immunoglobulin heavy-chain genes in the synovial tissue. Compared with blood, synovial immunoglobulin isotypes were skewed from IGHM and IGHD to IGHG3 and IGHG1. Further, patients with RA with high levels of lymphocytic infiltrates in the synovium demonstrated parallel skewing in their blood with a relative decrease in IGHGM (P < 0.002) and IGHD (P < 0.03) and an increase in class-switched immunoglobulin genes IGHG3 (P < 0.03) and IGHG1 (P < 0.002). CONCLUSION: High-resolution automated identification and quantification of synovial immune cell aggregates uncovered skewing in the synovium from naïve IGHD and IGHM to memory IGHG3 and IGHG1 and revealed that this process is reflected in the blood of patients with high inflammatory synovium.
Subject(s)
Arthritis, Rheumatoid , Osteoarthritis , Humans , Arthritis, Rheumatoid/genetics , Synovial Membrane/metabolism , Osteoarthritis/genetics , Autoantibodies/metabolism , Inflammation/metabolismABSTRACT
BACKGROUND: We sought to identify features that distinguish osteoarthritis (OA) and rheumatoid arthritis (RA) hematoxylin and eosin (H&E)-stained synovial tissue samples. METHODS: We compared fourteen pathologist-scored histology features and computer vision-quantified cell density (147 OA and 60 RA patients) in H&E-stained synovial tissue samples from total knee replacement (TKR) explants. A random forest model was trained using disease state (OA vs RA) as a classifier and histology features and/or computer vision-quantified cell density as inputs. RESULTS: Synovium from OA patients had increased mast cells and fibrosis (p < 0.001), while synovium from RA patients exhibited increased lymphocytic inflammation, lining hyperplasia, neutrophils, detritus, plasma cells, binucleate plasma cells, sub-lining giant cells, fibrin (all p < 0.001), Russell bodies (p = 0.019), and synovial lining giant cells (p = 0.003). Fourteen pathologist-scored features allowed for discrimination between OA and RA, producing a micro-averaged area under the receiver operating curve (micro-AUC) of 0.85±0.06. This discriminatory ability was comparable to that of computer vision cell density alone (micro-AUC = 0.87±0.04). Combining the pathologist scores with the cell density metric improved the discriminatory power of the model (micro-AUC = 0.92±0.06). The optimal cell density threshold to distinguish OA from RA synovium was 3400 cells/mm2, which yielded a sensitivity of 0.82 and specificity of 0.82. CONCLUSIONS: H&E-stained images of TKR explant synovium can be correctly classified as OA or RA in 82% of samples. Cell density greater than 3400 cells/mm2 and the presence of mast cells and fibrosis are the most important features for making this distinction.
Subject(s)
Arthritis, Rheumatoid , Osteoarthritis , Humans , Inflammation , Osteoarthritis/diagnosis , Arthritis, Rheumatoid/diagnosis , Synovial Membrane , Machine LearningABSTRACT
Periodontal disease is more common in individuals with rheumatoid arthritis (RA) who have detectable anti-citrullinated protein antibodies (ACPAs), implicating oral mucosal inflammation in RA pathogenesis. Here, we performed paired analysis of human and bacterial transcriptomics in longitudinal blood samples from RA patients. We found that patients with RA and periodontal disease experienced repeated oral bacteremias associated with transcriptional signatures of ISG15+HLADRhi and CD48highS100A2pos monocytes, recently identified in inflamed RA synovia and blood of those with RA flares. The oral bacteria observed transiently in blood were broadly citrullinated in the mouth, and their in situ citrullinated epitopes were targeted by extensively somatically hypermutated ACPAs encoded by RA blood plasmablasts. Together, these results suggest that (i) periodontal disease results in repeated breaches of the oral mucosa that release citrullinated oral bacteria into circulation, which (ii) activate inflammatory monocyte subsets that are observed in inflamed RA synovia and blood of RA patients with flares and (iii) activate ACPA B cells, thereby promoting affinity maturation and epitope spreading to citrullinated human antigens.
Subject(s)
Arthritis, Rheumatoid , Periodontal Diseases , Humans , Autoantibodies , Mouth Mucosa , Antibody Formation , Epitopes , BacteriaABSTRACT
To detail the unmet clinical and scientific needs in the field of rheumatology. After a 2-year hiatus due to the SARS-CoV-2 pandemic, the 22nd annual international Advances in Targeted Therapies meeting brought together more than 100 leading basic scientists and clinical researchers in rheumatology, immunology, epidemiology, molecular biology and other specialties. Breakout sessions were convened with experts in five rheumatological disease-specific groups including: rheumatoid arthritis (RA), psoriatic arthritis, axial spondyloarthritis, systemic lupus erythematosus and connective tissue diseases (CTDs). In each group, experts were asked to identify and prioritise current unmet needs in clinical and translational research, as well as highlight recent progress in meeting formerly identified unmet needs. Clinical trial design innovation was emphasised across all disease states. Within RA, developing therapies and trials for refractory disease patients remained among the most important identified unmet needs and within lupus and spondyloarthritis the need to account for disease endotypes was highlighted. The RA group also identified the need to better understand the natural history of RA, pre-RA states and the need ultimately for precision medicine. In CTD generally, experts focused on the need to better identify molecular, cellular and clinical signals of early and undifferentiated disease in order to identify novel drug targets. There remains a strong need to develop therapies and therapeutic strategies for those with treatment-refractory disease. Increasingly it is clear that we need to better understand the natural history of these diseases, including their 'predisease' states, and identify molecular signatures, including at a tissue level, which can facilitate disease diagnosis and treatment. As these unmet needs in the field of rheumatic diseases have been identified based on consensus of expert clinicians and scientists in the field, this document may serve individual researchers, institutions and industry to help prioritise their scientific activities.
Subject(s)
Arthritis, Psoriatic , Arthritis, Rheumatoid , COVID-19 , Rheumatic Diseases , Rheumatology , Humans , SARS-CoV-2 , Rheumatic Diseases/drug therapy , Rheumatic Diseases/diagnosis , Arthritis, Rheumatoid/drug therapy , Arthritis, Psoriatic/drug therapyABSTRACT
Background: Total hip arthroplasty (THA) and total knee arthroplasty (TKA) are cost-effective procedures that decrease pain and improve health-related quality of life for patients with advanced symptomatic arthritis, including rheumatoid arthritis (RA). Patients with RA have a longer length of stay (LOS) after THA or TKA than patients with osteoarthritis, yet the factors contributing to LOS have not been investigated. Purpose: We sought to identify the factors contributing to LOS for patients with RA undergoing THA and TKA at a single tertiary care orthopedic specialty hospital. Methods: We retrospectively reviewed data from a prospectively collected cohort of 252 RA patients undergoing either THA or TKA. Demographics, RA characteristics, medications, serologies, and disease activity were collected preoperatively. Linear regression was performed to explore the relationship between LOS (log-transformed) and possible predictors. A multivariate model was constructed through backward selection using significant predictors from a univariate analysis. Results: Of the 252 patients with RA, 83% were women; they had a median disease duration of 14 years and moderate disease activity at the time of arthroplasty. We had LOS data on 240 (95%) of the cases. The mean LOS was 3.4 ± 1.5 days. The multivariate analysis revealed a longer LOS for RA patients who underwent TKA versus THA, were women versus men, required a blood transfusion, and took preoperative opioids. Conclusion: Our retrospective study found that increased postoperative LOS in RA patients undergoing THA or TKA was associated with factors both non-modifiable (type of surgery, sex) and modifiable (postoperative blood transfusion, preoperative opioid use). These findings suggest that preoperative optimization of the patient with RA might focus on improving anemia and reducing opioid use in efforts to shorten LOS. More rigorous study is warranted.
ABSTRACT
Ectopic lymphoid structures (ELS) can develop in rheumatoid arthritis (RA) synovial tissue, but the precise pathways of B cell activation and selection are not well understood. Here, we identify a synovial B cell population characterized by co-expression of a family of orphan nuclear receptors (NR4A1-3), which is highly enriched in RA synovial tissue. A transcriptomic profile of NR4A synovial B cells significantly overlaps with germinal center light zone B cells and an accrual of somatic hypermutation that correlates with loss of naive B cell state. NR4A B cells co-express lymphotoxins α and ß and IL-6, supporting functions in ELS promotion. Expanded and shared clones between synovial NR4A B cells and plasma cells and the rapid upregulation with BCR stimulation point to in situ differentiation. Together, we identify a dynamic progression of B cell activation in RA synovial ELS, with NR4A transcription factors having an important role in local adaptive immune responses.
Subject(s)
Arthritis, Rheumatoid , Synovial Membrane , B-Lymphocytes , Humans , Plasma Cells/metabolism , Receptors, Cytoplasmic and Nuclear/metabolism , Synovial Membrane/metabolismABSTRACT
OBJECTIVE: We quantified inflammatory burden in rheumatoid arthritis (RA) synovial tissue by using computer vision to automate the process of counting individual nuclei in hematoxylin and eosin images. METHODS: We adapted and applied computer vision algorithms to quantify nuclei density (count of nuclei per unit area of tissue) on synovial tissue from arthroplasty samples. A pathologist validated algorithm results by labeling nuclei in synovial images that were mislabeled or missed by the algorithm. Nuclei density was compared with other measures of RA inflammation such as semiquantitative histology scores, gene-expression data, and clinical measures of disease activity. RESULTS: The algorithm detected a median of 112,657 (range 8,160-821,717) nuclei per synovial sample. Based on pathologist-validated results, the sensitivity and specificity of the algorithm was 97% and 100%, respectively. The mean nuclei density calculated by the algorithm was significantly higher (P < 0.05) in synovium with increased histology scores for lymphocytic inflammation, plasma cells, and lining hyperplasia. Analysis of RNA sequencing identified 915 significantly differentially expressed genes in correlation with nuclei density (false discovery rate is less than 0.05). Mean nuclei density was significantly higher (P < 0.05) in patients with elevated levels of C-reactive protein, erythrocyte sedimentation rate, rheumatoid factor, and cyclized citrullinated protein antibody. CONCLUSION: Nuclei density is a robust measurement of inflammatory burden in RA and correlates with multiple orthogonal measurements of inflammation.
ABSTRACT
To address the need for simple, safe, sensitive, and scalable SARS-CoV-2 tests, we validated and implemented a PCR test that uses a saliva collection kit use at home. Individuals self-collected 300 µl saliva in vials containing Darnell Rockefeller University Laboratory (DRUL) buffer and extracted RNA was assayed by RT-PCR (the DRUL saliva assay). The limit of detection was confirmed to be 1 viral copy/µl in 20 of 20 replicate extractions. Viral RNA was stable in DRUL buffer at room temperature up to seven days after sample collection, and safety studies demonstrated that DRUL buffer immediately inactivated virus at concentrations up to 2.75x106 PFU/ml. Results from SARS-CoV-2 positive nasopharyngeal (NP) swab samples collected in viral transport media and assayed with a standard FDA Emergency Use Authorization (EUA) test were highly correlated with samples placed in DRUL buffer. Direct comparison of results from 162 individuals tested by FDA EUA oropharyngeal (OP) or NP swabs with co-collected saliva samples identified four otherwise unidentified positive cases in DRUL buffer. Over six months, we collected 3,724 samples from individuals ranging from 3 months to 92 years of age. This included collecting weekly samples over 10 weeks from teachers, children, and parents from a pre-school program, which allowed its safe reopening while at-risk pods were quarantined. In sum, we validated a simple, sensitive, stable, and safe PCR-based test using a self-collected saliva sample as a valuable tool for clinical diagnosis and screening at workplaces and schools.
Subject(s)
COVID-19 Nucleic Acid Testing , COVID-19 , SARS-CoV-2 , Saliva/virology , Schools , Specimen Handling , COVID-19/diagnosis , COVID-19/genetics , Child , Female , Humans , MaleABSTRACT
OBJECTIVE: We sought to determine histologic and gene expression features of clinical improvement in early diffuse cutaneous systemic sclerosis (dcSSc; scleroderma). METHODS: Fifty-eight forearm biopsies were evaluated from 26 individuals with dcSSc in two clinical trials. Histologic/immunophenotypic assessments of global severity, alpha-smooth muscle actin (aSMA), CD34, collagen, inflammatory infiltrate, follicles and thickness were compared with gene expression and clinical data. Support vector machine learning was performed using scleroderma gene expression subset (normal-like, fibroproliferative, inflammatory) as classifiers and histology scores as inputs. Comparison of w-vector mean absolute weights was used to identify histologic features most predictive of gene expression subset. We then tested for differential gene expression according to histologic severity and compared those with clinical improvement (according to the Combined Response Index in Systemic Sclerosis). RESULTS: aSMA was highest and CD34 lowest in samples with highest local Modified Rodnan Skin Score. CD34 and aSMA changed significantly from baseline to 52 weeks in clinical improvers. CD34 and aSMA were the strongest predictors of gene expression subset, with highest CD34 staining in the normal-like subset (p<0.001) and highest aSMA staining in the inflammatory subset (p=0.016). Analysis of gene expression according to CD34 and aSMA binarised scores identified a 47-gene fibroblast polarisation signature that decreases over time only in improvers (vs non-improvers). Pathway analysis of these genes identified gene expression signatures of inflammatory fibroblasts. CONCLUSION: CD34 and aSMA stains describe distinct fibroblast polarisation states, are associated with gene expression subsets and clinical assessments, and may be useful biomarkers of clinical severity and improvement in dcSSc.
Subject(s)
Fibroblasts/physiology , Machine Learning , Scleroderma, Diffuse/genetics , Severity of Illness Index , Actins/metabolism , Adult , Antigens, CD34/metabolism , Clinical Trials as Topic , Collagen/metabolism , Female , Forearm , Gene Expression , Humans , Male , Middle Aged , Skin/metabolismABSTRACT
BACKGROUND: Rheumatoid arthritis, like many inflammatory diseases, is characterized by episodes of quiescence and exacerbation (flares). The molecular events leading to flares are unknown. METHODS: We established a clinical and technical protocol for repeated home collection of blood in patients with rheumatoid arthritis to allow for longitudinal RNA sequencing (RNA-seq). Specimens were obtained from 364 time points during eight flares over a period of 4 years in our index patient, as well as from 235 time points during flares in three additional patients. We identified transcripts that were differentially expressed before flares and compared these with data from synovial single-cell RNA-seq. Flow cytometry and sorted-blood-cell RNA-seq in additional patients were used to validate the findings. RESULTS: Consistent changes were observed in blood transcriptional profiles 1 to 2 weeks before a rheumatoid arthritis flare. B-cell activation was followed by expansion of circulating CD45-CD31-PDPN+ preinflammatory mesenchymal, or PRIME, cells in the blood from patients with rheumatoid arthritis; these cells shared features of inflammatory synovial fibroblasts. Levels of circulating PRIME cells decreased during flares in all 4 patients, and flow cytometry and sorted-cell RNA-seq confirmed the presence of PRIME cells in 19 additional patients with rheumatoid arthritis. CONCLUSIONS: Longitudinal genomic analysis of rheumatoid arthritis flares revealed PRIME cells in the blood during the period before a flare and suggested a model in which these cells become activated by B cells in the weeks before a flare and subsequently migrate out of the blood into the synovium. (Funded by the National Institutes of Health and others.).