Molecular Portraits of Early Rheumatoid Arthritis Identify Clinical and Treatment Response Phenotypes

Biologics or Janus Kinase Inhibitors in Rheumatoid Arthritis Patients Who are Insufficient Responders to Conventional Anti-Rheumatic Drugs

Rheumatoid arthritis (RA) is a chronic immune-mediated inflammatory disease which can induce progressive disability if not properly treated early. Over the last 20 years, the improvement of knowledge on the pathogenesis of the disease has made available several drugs targeting key elements of the pathogenetic process, which now represent the preferred treatment option after the failure of first-line therapy with conventional drugs such as methotrexate (MTX). To this category of targeted drugs belong anti-cytokine or cell-targeted biological agents and more recently also Janus kinase inhibitors (JAKis). In the absence to date of specific biomarkers to guide the therapeutic choice in the context of true precision medicine, the choice of the first targeted drug after MTX failure is guided by treatment cost (especially after the marketing of biosimilar products) and by the clinical characteristics of the patient (age, sex, comorbidities and compliance) and the disease (presence or absence of autoantibodies and systemic or extra-articular manifestations), which may influence the efficacy and safety profile of the available products. This viewpoint focuses on the decision-making process underlying the personalized approach to RA therapy and will analyse the evidence in the literature supporting the choice of individual products and in particular the differential choice between biological drugs and JAKis.

JAK/STAT signaling in rheumatoid arthritis leukocytes is uncoupled from serum cytokines in a subset of patients

OBJECTIVE

Rheumatoid Arthritis (RA) is a systemic autoimmune disease involving pro-inflammatory cytokines that can be therapeutically targeted by antibodies or kinase inhibitors. Nevertheless, these drugs fail in a subset of patients independent of the abundance of the targeted cytokines. We aim to explore the cellular basis of this phenomenon by analyzing the relation of cytokine abundance and activation of downstream signaling pathways in RA.

METHODS

The study included 62 RA patients and 9 healthy controls (HC). Phosphorylation of STAT 1-6 in various immune cell subsets was determined ex vivo using a novel robust flow cytometry-based protocol. Serum concentrations of IL-6, IL-10, IL-12p70, IL-17 A, interferon gamma, and TNFα in the same samples were measured using highly sensitive single molecule array (SIMOA).

RESULTS

We found an increase in circulating cytokines in RA patients, while STAT activity was lower in RA patients compared to HC. Based on STAT activity we determined three endotypes in active RA patients (cDAI>10, n = 28): 1) those with active STAT5a/b signaling in T cells (n = 7/28), 2) those with a low STAT activity in all assessed cell types (n = 14/28), and 3) those with active STAT1 and STAT3 signaling mainly in myeloid cells (n = 7/28). Integrating intracellular STAT activation and cytokine analysis revealed diminished JAK/STAT signaling in a subset of patients (n = 8/20) despite elevated serum cytokine concentrations.

CONCLUSION

Diminished JAK/STAT signaling in active RA may partly explain unresponsiveness to therapy targeting cytokine signaling. Analysis of JAK/STAT phosphorylation may identify patients at risk for non-response to these therapies.

Deciphering the molecular landscape of rheumatoid arthritis offers new insights into the stratified treatment for the condition

Background

For Rheumatoid Arthritis (RA), a long-term chronic illness, it is essential to identify and describe patient subtypes with comparable goal status and molecular biomarkers. This study aims to develop and validate a new subtyping scheme that integrates genome-scale transcriptomic profiles of RA peripheral blood genes, providing a fresh perspective for stratified treatments.

Methods

We utilized independent microarray datasets of RA peripheral blood mononuclear cells (PBMCs). Up-regulated differentially expressed genes (DEGs) were subjected to functional enrichment analysis. Unsupervised cluster analysis was then employed to identify RA peripheral blood gene expression-driven subtypes. We defined three distinct clustering subtypes based on the identified 404 up-regulated DEGs.

Results

Subtype A, named NE-driving, was enriched in pathways related to neutrophil activation and responses to bacteria. Subtype B, termed interferon-driving (IFN-driving), exhibited abundant B cells and showed increased expression of transcripts involved in IFN signaling and defense responses to viruses. In Subtype C, an enrichment of CD8+ T-cells was found, ultimately defining it as CD8+ T-cells-driving. The RA subtyping scheme was validated using the XGBoost machine learning algorithm. We also evaluated the therapeutic outcomes of biological disease-modifying anti-rheumatic drugs.

Conclusions

The findings provide valuable insights for deep stratification, enabling the design of molecular diagnosis and serving as a reference for stratified therapy in RA patients in the future.

Dynamic Release from Acetalated Dextran Nanoparticles for Precision Therapy of Inflammation

Polymer-based nanoparticles (NPs) that react to altered physiological characteristics have the potential to enhance the delivery of therapeutics to a specific area. These materials can utilize biochemical triggers, such as low pH, which is prone to happen locally in an inflammatory microenvironment due to increased cellular activity. This reduced pH is neutralized when inflammation subsides. For precise delivery of therapeutics to match this dynamic reaction, drug delivery systems (DDS) need to not only release the drug (ON) but also stop the release (OFF) autonomously. In this study, we use a systematic approach to optimize the composition of acetalated dextran (AcDex) NPs to start (ON) and stop (OFF) releasing model cargo, depending on local pH changes. By mixing ratios of AcDex polymers (mixed NPs), we achieved a highly sensitive material that was able to rapidly release cargo when going from pH 7.4 to pH 6.0. At the same time, the mix also offered a stable composition that enabled a rapid ON/OFF/ON/OFF switching within this narrow pH range in only 90 min. These mixed NPs were also sensitive to biological pH changes, with increased release in the presence of inflammatory cells compared to healthy cells. Such precise and controllable characteristics of a DDS position mixed NPs as a potential treatment platform to inhibit disease flare-ups, reducing both systemic and local side effects to offer a superior treatment option for inflammation compared to conventional systems.

The Role of Alpha-Synuclein in Synucleinopathy: Impact on Lipid Regulation at Mitochondria-ER Membranes

The protein alpha-synuclein (αSyn) plays a critical role in the pathogenesis of synucleinopathy, which includes Parkinson's disease and multiple system atrophy, and mounting evidence suggests that lipid dyshomeostasis is a critical phenotype in these neurodegenerative conditions. Previously, we identified that αSyn localizes to mitochondria-associated endoplasmic reticulum membranes (MAMs), temporary functional domains containing proteins that regulate lipid metabolism, including the de novo synthesis of phosphatidylserine. In the present study, we have analyzed the lipid composition of postmortem human samples, focusing on the substantia nigra pars compacta of Parkinson's disease and controls, as well as three less affected brain regions of Parkinson's donors. To further assess synucleinopathy-related lipidome alterations, similar analyses were performed on the striatum of multiple system atrophy cases. Our data show region-and disease-specific changes in the levels of lipid species. Specifically, our data revealed alterations in the levels of specific phosphatidylserine species in brain areas most affected in Parkinson's disease. Some of these alterations, albeit to a lesser degree, are also observed multiples system atrophy. Using induced pluripotent stem cell-derived neurons, we show that αSyn contributes to regulating phosphatidylserine metabolism at MAM domains, and that αSyn dosage parallels the perturbation in phosphatidylserine levels. Our results support the notion that αSyn pathophysiology is linked to the dysregulation of lipid homeostasis, which may contribute to the vulnerability of specific brain regions in synucleinopathy. These findings have significant therapeutic implications.

Phosphoproteomic profiling of early rheumatoid arthritis synovium reveals active signalling pathways and differentiates inflammatory pathotypes

BACKGROUND

Kinases are intracellular signalling mediators and key to sustaining the inflammatory process in rheumatoid arthritis (RA). Oral inhibitors of Janus Kinase family (JAKs) are widely used in RA, while inhibitors of other kinase families e.g. phosphoinositide 3-kinase (PI3K) are under development. Most current biomarker platforms quantify mRNA/protein levels, but give no direct information on whether proteins are active/inactive. Phosphoproteome analysis has the potential to measure specific enzyme activation status at tissue level.

METHODS

We validated the feasibility of phosphoproteome and total proteome analysis on 8 pre-treatment synovial biopsies from treatment-naive RA patients using label-free mass spectrometry, to identify active cell signalling pathways in synovial tissue which might explain failure to respond to RA therapeutics.

RESULTS

Differential expression analysis and functional enrichment revealed clear separation of phosphoproteome and proteome profiles between lymphoid and myeloid RA pathotypes. Abundance of specific phosphosites was associated with the degree of inflammatory state. The lymphoid pathotype was enriched with lymphoproliferative signalling phosphosites, including Mammalian Target Of Rapamycin (MTOR) signalling, whereas the myeloid pathotype was associated with Mitogen-Activated Protein Kinase (MAPK) and CDK mediated signalling. This analysis also highlighted novel kinases not previously linked to RA, such as Protein Kinase, DNA-Activated, Catalytic Subunit (PRKDC) in the myeloid pathotype. Several phosphosites correlated with clinical features, such as Disease-Activity-Score (DAS)-28, suggesting that phosphosite analysis has potential for identifying novel biomarkers at tissue-level of disease severity and prognosis.

CONCLUSIONS

Specific phosphoproteome/proteome signatures delineate RA pathotypes and may have clinical utility for stratifying patients for personalised medicine in RA.

The fibroid phenotype of biological naïve patients with rheumatoid arthritis are less likely to respond to anti-IL-6R treatment

Type III collagen gene expression is upregulated in the synovium of patients with rheumatoid arthritis (RA) presenting the fibroid phenotype. The soluble type III collagen formation biomarker, PRO-C3, is known to measure fibrogenesis in fibrotic diseases. In this exploratory study, we aimed to investigate the association between fibrogenesis (PRO-C3) and the disease- and treatment response in patients with RA. We measured PRO-C3 in subsets of two clinical trials assessing the effect of the anti-interleukin-6 (IL-6) receptor treatment tocilizumab (TCZ) as monotherapy or polytherapy with methotrexate. PRO-C3 levels had weak or very weak correlations with the clinical parameters (Spearman's). However, when the patients were divided into Disease Activity Score-28 groups characterized by the erythrocyte sedimentation rate (DAS28-ESR), there was a statistical difference between the PRO-C3 levels of the different groups (p < 0.05). To determine the response in relation to PRO-C3, a cut-off based on PRO-C3 levels and patients in remission (DAS28-ESR ≤ 2.6) was identified. This showed that a reduction in PRO-C3 after treatment initiation was associated with decreased DAS28-ESR and a higher response rate in patients with low PRO-C3 levels than in those with high PRO-C3 levels. This indicates that a fibrotic component affects the responsiveness of patients.

Notch ligands are biomarkers of anti-TNF response in RA patients

Notch and its ligands play a critical role in rheumatoid arthritis (RA) pathogenesis. Hence, studies were conducted to delineate the functional significance of the Notch pathway in RA synovial tissue (ST) cells and the influence of RA therapies on their expression. Morphological studies reveal that JAG1, DLL4, and Notch1 are highly enriched in RA ST lining and sublining CD68+CD14+ MΦs. JAG1 and DLL4 transcription is jointly upregulated in RA MΦs reprogrammed by TLR4/5 ligation and TNF, whereas Syntenin-1 exposure expands JAG1, DLL4, and Notch1 expression levels in these cells. Single-cell RNA-seq data exhibit that JAG1 and Notch3 are overexpressed on all fibroblast-like synoviocyte (FLS) subpopulations, in parallel, JAG2, DLL1, and Notch1 expression levels are modest on RA FLS and are predominately potentiated by TLR4 ligation. Intriguingly, JAG1, DLL1/4, and Notch1/3 are presented on RA endothelial cells, and their expression is mutually reconfigured by TLR4/5 ligation in the endothelium. Synovial JAG1/JAG2/DLL1 or Notch1/3 transcriptomes were unchanged in patients who received disease-modifying anti-rheumatic drugs (DMARDs) or IL-6R Ab therapy regardless of disease activity score. Uniquely, RA MΦs and endothelial cells rewired by IL-6 displayed DLL4 transcriptional upregulation, and IL-6R antibody treatment disrupted RA ST DLL4 transcription in good responders compared to non-responders or moderate responders. Nevertheless, the JAG1/JAG2/DLL1/DLL4 transcriptome was diminished in anti-TNF good responders with myeloid pathotype and was unaltered in the fibroid pathotype except for DLL4. Taken together, our findings suggest that RA myeloid Notch ligands can serve as markers for anti-TNF responsiveness and trans-activate Notch receptors expressed on RA FLS and/or endothelial cells.

Molecular portrait of chronic joint diseases: Defining endotypes toward personalized medicine

Joint diseases affect hundreds of millions of people worldwide, and their prevalence is constantly increasing. To date, despite recent advances in the development of therapeutic options for most rheumatic conditions, a significant proportion of patients still lack efficient disease management, considerably impacting their quality of life. Through the spectrum of rheumatoid arthritis (RA), psoriatic arthritis (PsA), and osteoarthritis (OA) as quintessential and common rheumatic diseases, this review first provides an overview of their epidemiological and clinical features before exploring how the better definition of clinical phenotypes has helped their clinical management. It then discusses the recent progress in understanding the diversity of endotypes underlying disease phenotypes. Finally, this review highlights the current challenges of implementing molecular endotypes towards the personalized management of RA, PsA and OA patients in the future.

Autoimmunity to stromal-derived autoantigens in rheumatoid ectopic germinal centers exacerbates arthritis and affects clinical response

Ectopic lymphoid structures (ELSs) in the rheumatoid synovial joints sustain autoreactivity against locally expressed autoantigens. We recently identified recombinant monoclonal antibodies (RA-rmAbs) derived from single, locally differentiated rheumatoid arthritis (RA) synovial B cells, which specifically recognize fibroblast-like synoviocytes (FLSs). Here, we aimed to identify the specificity of FLS-derived autoantigens fueling local autoimmunity and the functional role of anti-FLS antibodies in promoting chronic inflammation. A subset of anti-FLS RA-rmAbs reacting with a 60 kDa band from FLS extracts demonstrated specificity for HSP60 and partial cross-reactivity to other stromal autoantigens (i.e., calreticulin/vimentin) but not to citrullinated fibrinogen. Anti-FLS RA-rmAbs, but not anti-neutrophil extracellular traps rmAbs, exhibited pathogenic properties in a mouse model of collagen-induced arthritis. In patients, anti-HSP60 antibodies were preferentially detected in RA versus osteoarthritis (OA) synovial fluid. Synovial HSPD1 and CALR gene expression analyzed using bulk RNA-Seq and GeoMx-DSP closely correlated with the lympho-myeloid RA pathotype, and HSP60 protein expression was predominantly observed around ELS. Moreover, we observed a significant reduction in synovial HSP60 gene expression followed B cell depletion with rituximab that was strongly associated with the treatment response. Overall, we report that synovial stromal-derived autoantigens are targeted by pathogenic autoantibodies and are associated with specific RA pathotypes, with potential value for patient stratification and as predictors of the response to B cell-depleting therapies.

Clinical Phenotypes, Serological Biomarkers, and Synovial Features Defining Seropositive and Seronegative Rheumatoid Arthritis: A Literature Review

Rheumatoid arthritis (RA) is a chronic autoimmune disorder which can lead to long-term joint damage and significantly reduced quality of life if not promptly diagnosed and adequately treated. Despite significant advances in treatment, about 40% of patients with RA do not respond to individual pharmacological agents and up to 20% do not respond to any of the available medications. To address this large unmet clinical need, several recent studies have focussed on an in-depth histological and molecular characterisation of the synovial tissue to drive the application of precision medicine to RA. Currently, RA patients are clinically divided into "seropositive" or "seronegative" RA, depending on the presence of routinely checked antibodies. Recent work has suggested that over the last two decades, long-term outcomes have improved significantly in seropositive RA but not in seronegative RA. Here, we present up-to-date differences in epidemiology, clinical features, and serological biomarkers in seronegative versus seropositive RA and discuss how histological and molecular synovial signatures, revealed by recent large synovial biopsy-based clinical trials, may be exploited to refine the classification of RA patients, especially in the seronegative group.

Synovial fibroblast gene expression is associated with sensory nerve growth and pain in rheumatoid arthritis

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.

Amino Acids Fueling Fibroblast-Like Synoviocyte Activation and Arthritis By Regulating Chemokine Expression and Leukocyte Migration

OBJECTIVE

We analyzed the impact of amino acid (AA) availability on the inflammatory response in arthritis.

METHODS

We stimulated rheumatoid arthritis (RA) fibroblast-like synoviocytes (FLSs) with tumor necrosis factor (TNF) in the presence or absence of proteinogenic AAs and measured their response by QuantSeq 3' messenger RNA sequencing, quantitative polymerase chain reaction, and enzyme-linked immunosorbent assay. Signal transduction events were determined by Western blot. We performed K/BxN serum transfer arthritis in mice receiving a normal and a low-protein diet and analyzed arthritis clinically and histologically.

RESULTS

Deprivation of AAs decreased the expression of a specific subset of genes, including the chemokines CXCL10, CCL2, and CCL5 in TNF-stimulated FLSs. Mechanistically, the presence of AAs was required for the TNF-induced activation of an interferon regulatory factor 1 (IRF1)-STAT1 signaling circuit that drives the expression of chemotactic factors. The expression of IRF1 and the IRF1-dependent gene set in FLSs was highly correlated with the presence of inflammatory cells in human RA, emphasizing the important role of this AA-dependent pathway in inflammatory cell recruitment to the synovial tissue. Finally, we show that mice receiving a low-protein diet expressed less IRF1 in the inflamed synovium and consequently developed reduced clinical and histologic signs of arthritis.

CONCLUSION

AA deprivation reduces the severity of arthritis by suppressing the expression of IRF1-STAT1-driven chemokines, which are crucial for leukocyte recruitment to the arthritic joint. Overall, our study provides novel insights into critical determinants of inflammatory arthritis and may pave the way for dietary intervention trials in RA.

shinyExprPortal: a configurable 'shiny' portal for sharing analysis of molecular expression data

MOTIVATION

The scale of omics research presents many obstacles to full sharing and access to analysis results. Current publication models impose limits on the number of pages and figures, requiring careful preparation and selection of content. At the same time, depositing data in open repositories significantly shifts the burden of access and reproduction to readers, who may include people who are not programmers or analysts.

RESULTS

We introduce shinyExprPortal, an R package that implements omics web portals with minimal coding effort. The portals allow exploration of transcriptomic or proteomic expression data and phenotypes, showcasing results of various types of analysis including differential expression, co-expression and pathways analysis. The integration with bioinformatics workflows enables researchers to focus on their results and share findings using interactive and publication-quality plots.

AVAILABILITY AND IMPLEMENTATION

The shinyExprPortal package is available to download and install from CRAN and https://github.com/C4TB/shinyExprPortal.

Deconvolution of synovial myeloid cell subsets across pathotypes and role of COL3A1+ macrophages in rheumatoid arthritis remission

Background

Monocyte/macrophage (Mo/Mp) is a critical cell population involved in immune modulation of rheumatoid synovitis (RA) across different pathotypes. This study aims to investigate the contribution of Mo/Mp clusters to RA activity, and the biological function of particular subtypes in RA remission.

Methods

We integrated single-cell RNA sequencing datasets from 4 published and 1 in-house studies using Liger selected by comparison. We estimated the abundance of Mo/Mp subtypes in bulk RNA-seq data from the 81 patients of the Pathobiology of Early Arthritis Cohort (PEAC) using deconvolution analysis. Correlations between Mo/Mp subtypes and RA clinical metrics were assessed. A particular cell type was identified using multicolor immunofluorescence and flow cytometry in vivo and successfully induced from a cell line in vitro. Potential immune modulation function of it was performed using immunohistochemical staining, adhesion assay, and RT-qPCR.

Results

We identified 8 Mo/Mp clusters. As a particular subtype among them, COL3A1+ Mp (CD68+, COL3A1+, ACTA2-) enriched in myeloid pathotype and negatively correlated with RA severity metrics in all pathotypes. Flow cytometry and multicolor immunofluorescence evidenced the enrichment and M2-like phenotype of COL3A1+ Mp in the myeloid pathotype. Further assays suggested that COL3A1+ Mp potentially attenuates RA severity via expressing anti-inflammatory cytokines, enhancing Mp adhesion, and forming a physical barrier at the synovial lining.

Conclusion

This study reported unexplored associations between different pathologies and myeloid cell subtypes. We also identified a fibroblast-and-M2-like cluster named COL3A1+ Mp, which potentially contributes to synovial immune homeostasis. Targeting the development of COL3A1+ Mp may hold promise for inducing RA remission.

Synovial 5-Lipoxygenase-Derived Oxylipins Define a Lympho-Myeloid-Enriched Synovium

OBJECTIVE

Oxylipins are bioactive lipids derived from polyunsaturated fatty acids (PUFAs) that modulate inflammation and may remain overexpressed in refractory synovitis. In plasma, they could also be biomarkers of synovial pathology. The aim of this study is to determine if synovial oxylipins in inflamed joints correlate with plasma oxylipins and with synovial histologic patterns.

METHODS

Patients with established rheumatoid or psoriatic arthritis with active disease despite treatment were recruited, and paired synovial tissue (ST) and plasma were collected. Oxylipins were determined by liquid chromatography with tandem mass spectrometry and were classified into groups according to their PUFA precursor and enzyme. The expression of CD20, CD68, CD3, and CD138 was obtained to describe synovial histology. Cell-specific expression of oxylipin-related genes was identified by examining available synovial single-cell RNA sequencing data.

RESULTS

We included a total of 32 ST and 26 paired-plasma samples. A total of 71 oxylipins were identified in ST, but only 24 were identified in plasma. Only levels of 9,10-dihydroxyoctadecenoic acid and tetranor-Prostaglandin FM had a significant positive correlation between plasma and ST. Several oxylipins and oxylipin-related genes were differentially expressed among synovial phenotypes. Specifically, several 5-lipoxygenase (LOX)-derived oxylipins were statistically elevated in the lympho-myeloid phenotype and associated with B cell expression in rheumatoid arthritis samples.

CONCLUSION

The lack of correlation between ST and plasma oxylipins suggests that ST lipid profiling better characterizes active pathways in treated joints. Synovial 5-LOX-derived oxylipins were highly expressed in lympho-myeloid-enriched synovium. Combination therapy with 5-LOX inhibitors to improve refractory inflammation may be needed in patients with this histologic group.

Axl and MerTK regulate synovial inflammation and are modulated by IL-6 inhibition in rheumatoid arthritis

The TAM tyrosine kinases, Axl and MerTK, play an important role in rheumatoid arthritis (RA). Here, using a unique synovial tissue bioresource of patients with RA matched for disease stage and treatment exposure, we assessed how Axl and MerTK relate to synovial histopathology and disease activity, and their topographical expression and longitudinal modulation by targeted treatments. We show that in treatment-naive patients, high AXL levels are associated with pauci-immune histology and low disease activity and inversely correlate with the expression levels of pro-inflammatory genes. We define the location of Axl/MerTK in rheumatoid synovium using immunohistochemistry/fluorescence and digital spatial profiling and show that Axl is preferentially expressed in the lining layer. Moreover, its ectodomain, released in the synovial fluid, is associated with synovial histopathology. We also show that Toll-like-receptor 4-stimulated synovial fibroblasts from patients with RA modulate MerTK shedding by macrophages. Lastly, Axl/MerTK synovial expression is influenced by disease stage and therapeutic intervention, notably by IL-6 inhibition. These findings suggest that Axl/MerTK are a dynamic axis modulated by synovial cellular features, disease stage and treatment.

Patterns and determinants of response to novel therapies in juvenile and adult-onset polyarthritis

Biologic and targeted synthetic DMARDs (b/tsDMARDs) have revolutionized the management of multiple rheumatic inflammatory conditions. Among these, polyarticular JIA (pJIA) and RA display similarities in terms of disease pathophysiology and response pattern to b/tsDMARDs. Indeed, the therapeutic efficacy of novel targeted drugs is variable among individual patients, in both RA and pJIA. The mechanisms and determinants of this heterogeneous response are diverse and complex, such that the development of true 'precision'-medicine strategies has proven highly challenging. In this review, we will discuss pathophysiological, patient-specific, drug-specific and environmental factors contributing to individual therapeutic response in pJIA in comparison with what is known in RA. Although some biomarkers have been identified that stratify with respect to the likelihood of either therapeutic response or non-response, few have proved useful in clinical practice so far, likely due to the complexity of treatment-response mechanisms. Consequently, we propose a pragmatic, patient-centred and clinically based approach, i.e. personalized instead of biomarker-based precision medicine in JIA.

Synovial Tissue Insights into Heterogeneity of Rheumatoid Arthritis

PURPOSE OF REVIEW

Rheumatoid arthritis is one of the most common rheumatic and autoimmune diseases. While it can affect many different organ systems, RA primarily involves inflammation in the synovium, the tissue that lines joints. Patients with RA exhibit significant clinical heterogeneity in terms of presence or absence of autoantibodies, degree of permanent deformities, and most importantly, treatment response. These clinical characteristics point to heterogeneity in the cellular and molecular pathogenesis of RA, an area that several recent studies have begun to address.

RECENT FINDINGS

Single-cell RNA-sequencing initiatives and deeper focused studies have revealed several RA-associated cell populations in synovial tissues, including peripheral helper T cells, autoimmunity-associated B cells (ABCs), and NOTCH3+ sublining fibroblasts. Recent large transcriptional studies and translational clinical trials present frameworks to capture cellular and molecular heterogeneity in RA synovium. Technological developments, such as spatial transcriptomics and machine learning, promise to further elucidate the different types of RA synovitis and the biological mechanisms that characterize them, key elements of precision medicine to optimize patient care and outcomes in RA. This review recaps the findings of those recent studies and puts our current knowledge and future challenges into scientific and clinical perspective.

Pathophysiology and stratification of treatment-resistant rheumatoid arthritis

Early diagnosis and timely therapeutic intervention are clinical challenges of rheumatoid arthritis (RA), especially for treatment-resistant or difficult-to-treat patients. Little is known about the immunological mechanisms involved in refractory RA. In this review, we summarize previous research findings on the immunological mechanisms of treatment-resistant RA. Genetic prediction of treatment-resistant RA is challenging. Patients with and without anti-cyclic citrullinated peptide autoantibodies are considered part of distinct subgroups, especially regarding long-term clinical prognosis and treatment responses. B cells, T cells and other immune cells and fibroblasts are of pathophysiological importance and are associated with treatment responses. Finally, we propose a new hypothesis that stratifies patients with RA into two subgroups with distinct immunological pathologies based on our recent immunomics analysis of RA. One RA subgroup with a favorable prognosis is characterized by increased interferon signaling. Another subgroup with a worse prognosis is characterized by enhanced acquired immune responses. Increases in dendritic cell precursors and diversified autoreactive anti-modified protein antibodies may have pathophysiological roles, especially in the latter subgroup. These findings that improve treatment response predictions might contribute to future precision medicine for RA.

Expression quantitative trait loci analysis in rheumatoid arthritis identifies tissue specific variants associated with severity and outcome

OBJECTIVE

Genome-wide association studies have successfully identified more than 100 loci associated with susceptibility to rheumatoid arthritis (RA). However, our understanding of the functional effects of genetic variants in causing RA and their effects on disease severity and response to treatment remains limited.

METHODS

In this study, we conducted expression quantitative trait locus (eQTL) analysis to dissect the link between genetic variants and gene expression comparing the disease tissue against blood using RNA-Sequencing of synovial biopsies (n=85) and blood samples (n=51) from treatment-naïve patients with RA from the Pathobiology of Early Arthritis Cohort.

RESULTS

This identified 898 eQTL genes in synovium and genes loci in blood, with 232 genes in common to both synovium and blood, although notably many eQTL were tissue specific. Examining the HLA region, we uncovered a specific eQTL at HLA-DPB2 with the critical triad of single-nucleotide polymorphisms (SNPs) rs3128921 driving synovial HLA-DPB2 expression, and both rs3128921 and HLA-DPB2 gene expression correlating with clinical severity and increasing probability of the lympho-myeloid pathotype.

CONCLUSIONS

This analysis highlights the need to explore functional consequences of genetic associations in disease tissue. HLA-DPB2 SNP rs3128921 could potentially be used to stratify patients to more aggressive treatment immediately at diagnosis.

Role of IFN-α in Rheumatoid Arthritis

PURPOSE OF REVIEW

Type 1 interferons (IFN-I) are of increasing interest across a wide range of autoimmune rheumatic diseases. Historically, research into their role in rheumatoid arthritis (RA) has been relatively neglected, but recent work continues to highlight a potential contribution to RA pathophysiology.

RECENT FINDINGS

We emphasise the importance of disease stage when examining IFN-I in RA and provide an overview on how IFN-I may have a direct role on a variety of relevant cellular functions. We explore how clinical trajectory may be influenced by increased IFN-I signalling, and also, the limitations of scores composed of interferon response genes. Relevant environmental triggers and inheritable RA genetic risk relating to IFN-I signalling are explored with emphasis on intriguing data potentially linking IFN-I exposure, epigenetic changes, and disease relevant processes. Whilst these data cumulatively illustrate a likely role for IFN-I in RA, they also highlight the knowledge gaps, particularly in populations at risk for RA, and suggest directions for future research to both better understand IFN-I biology and inform targeted therapeutic strategies.

Tofacitinib restores psoriatic arthritis fibroblast-like synoviocytes function via autophagy and mitochondrial quality control modulation

OBJECTIVES

To evaluate the in vitro effect of tofacitinib on autophagy activity of psoriatic arthritis (PsA) fibroblast-like synoviocytes (FLS), and to confirm its activity on inflammatory and invasive properties of FLS and synovial cells, deepening the impact on mitochondrial function.

METHODS

FLS, peripheral blood mononuclear cells (PBMCs), and synovial cells from active PsA patients were cultured with tofacitinib 1 μM or vehicle control for 24 h. Autophagy was measured by Western blot and by fluorescence microscopy. Chemokines/cytokines released into culture supernatants were quantified by ELISA, while invasive properties of FLS by migration assays. Specific mitochondrial probes were adopted to measure intracellular reactive oxygen species (ROS), mitochondrial potential, morphology, turnover and mitophagy. Oxygen consumption rate (OCR), reflecting oxidative phosphorylation, was quantified using the Seahorse technology. Differences were determined by adopting the non-parametric Wilcoxon signed rank test.

RESULTS

18 patients with moderately-to-severely active PsA were enrolled. Tofacitinib significantly increased the levels of the autophagy markers LC3-II and ATG7 in PsA FLS compared to vehicle control, suggesting an increase in spontaneous autophagy activity; no effect was highlighted in PBMCs and synovial cells cultures. Tofacitinib reduced migration properties of PsA FLS, and reduced MCP-1 and IL-6 release into FLS and synovial cells cultures supernatants. Furthermore, tofacitinib decreased intracellular ROS production, increased basal OCR, ATP production and maximal respiratory capacity, and enhanced mitophagy and mitochondrial turnover.

CONCLUSIONS

The JAK inhibitor tofacitinib reduces the pro-invasive and pro-inflammatory properties of PsA FLS. Autophagy induction and mitochondrial quality control modulation by tofacitinib might contribute to FLS function restoration.

Comparative analysis of late-stage rheumatoid arthritis and osteoarthritis reveals shared histopathological features

OBJECTIVES

Osteoarthritis (OA) is a debilitating and heterogeneous condition, characterized by various levels of articular cartilage degradation, osteophytes formation, and synovial inflammation. Multiple evidences suggest that synovitis may appear early in the disease development and correlates with disease severity and pain, therefore representing a relevant therapeutic target. In a typical synovitis-driven joint disease, namely rheumatoid arthritis (RA), several pathotypes have been described by our group and associated with clinical phenotypes, disease progression, and response to therapy. However, whether these pathotypes can be also observed in the OA synovium is currently unknown.

METHODS

Here, using histological approaches combined with semi-quantitative scoring and quantitative digital image analyses, we comparatively characterize the immune cell infiltration in a large cohort of OA and RA synovial tissue samples collected at the time of total joint replacement.

RESULTS

We demonstrate that OA synovium can be categorized also into three pathotypes and characterized by disease- and stage-specific features. Moreover, we revealed that pathotypes specifically reflect distinct levels of peripheral inflammation.

CONCLUSIONS

In this study, we provide a novel and relevant pathological classification of OA synovial inflammation. Further studies investigating synovial molecular pathology in OA may contribute to the development of disease-modifying therapies.

Predicting best treatment in rheumatoid arthritis

BACKGROUND

Although targeted biological treatments have transformed the outlook for patients with rheumatoid arthritis (RA), 40% of patients show poor clinical response, and there is an imperative to unravel the molecular pathways and mechanisms underlying non-response and disease progression. 5-20% of RA individuals do not respond to all current medications including biologic and targeted therapies, which suggests that distinct pathogenic processes underlie multi-drug refractoriness.

OBJECTIVES

In this brief review we discuss advances from recent studies in precision medicine in rheumatoid arthritis.

METHODS

Bulk RNA-Sequencing of synovial biopsies from RA individuals combined with histology and deep clinical phenotyping has revealed substantial insights into divergent pathogenic pathways which lead to disease progression and illuminated mechanisms underlying failure to response to specific treatments. Biopsy-driven randomised controlled trials, such as R4RA and the forthcoming STRAP trial, have enabled the development of machine learning predictive models for predicting response to different therapies.

RESULTS

In the Pathobiology of Early Arthritis Cohort (PEAC), gene expression analysis showed that individuals could be classified into three gene expression subgroups which correlated with histopathotypes defined by histological markers: pauci-immune fibroid pathotype characterised by fibroblasts and an absence of immune inflammatory cells; diffuse-myeloid pathotype characterised by macrophage influx; and the lympho-myeloid pathotype delineated by the presence of B cells, but typically containing a complex inflammatory infiltrate with ectopic lymphoid structure formation. In the R4RA biopsy-driven randomised controlled trial, patients were randomised to either rituximab or tocilizumab. Comprehensive analysis of synovial biopsies pre/post-treatment identified gene signatures of response associated with pathogenic pathways which could be tracked over time. A group of true refractory patients were identified who had failed anti-TNF prior to the study (it was an entry criterion) and then subsequently failed both trial biologics during the trial. RNA-Seq analysis and digital spatial profiling identified specific cell types including DKK3+ fibroblasts as being associated with the refractory state. We identified machine learning predictive models based on specific gene signatures which were able to predict future response to therapy as well as the refractory state.

CONCLUSIONS

RNA-sequencing of synovial biopsies has enabled substantial progress in understanding disease endotypes in RA and identifying synovial gene signatures which predict prognosis and future response to treatment.

Distinct long-term disease activity trajectories differentiate early on treatment with etanercept in both rheumatoid arthritis and spondylarthritis patients: a prospective cohort study

To characterize disease activity trajectories and compare long-term drug retention between rheumatoid (RA) and spondylarthritis (SpA) patients initiating tumor necrosis factor inhibitor (TNFi) treatment (etanercept). Prospective observational study of RA, axial (AxSpA) and peripheral SpA (PerSpA) patients initiating etanercept during 2004-2020. Kaplan-Meier plots were used for drug retention comparisons and multivariable Cox regression models for predictors of discontinuation. Long-term disease activity trajectories were identified by latent class growth models using DAS28-ESR or ASDAS-CRP as outcome for RA and AxSpA respectively. We assessed 711 patients (450 RA, 178 AxSpA and 83 PerSpA) with a median (IQR) follow-up of 12 (5-32) months. At 5 years, 22%, 30% and 21% of RA, AxSpA and PerSpA patients, respectively, remained on therapy. Etanercept discontinuation was independent of the diagnosis and was predicted by gender and obesity in both RA and SpA groups. Four disease activity (DA) trajectories were identified from 6th month of treatment in both RA and AxSpA. RA patients in remission-low DA groups (33.7%) were younger, had shorter disease duration, fewer comorbidities and lower baseline disease activity compared to moderate (40.6%) & high DA (25.7%) groups. In AxSpA 74% were in inactive-low DA and they were more often males, non-obese and had lower number of comorbidities compared to higher ASDAS-CRP trajectories. In RA and AxSpA patients, disease activity trajectories revealed heterogeneity of TNFi treatment responses and prognosis. Male gender, lower baseline disease activity and fewer comorbidities, characterize a favourable outcome in terms of disease burden accrual and TNFi survival.

High-grade synovitis associates with clinical markers and response to therapy in chronic inflammatory arthritis: post hoc analysis of a synovial biomarkers prospective cohort study

Background

Inflammatory arthritis (IAs), such as rheumatoid arthritis (RA) and psoriatic arthritis (PsA), are characterized by the presence of chronic synovitis. The Krenn's synovitis score (KSS), a simple tool detectable by haematoxylin/eosin staining of synovial biopsy samples, allows the discrimination between high-grade and low-grade synovitis. The aim of this study was to identify the clinical associations of KSS and to evaluate the relationship between high-grade synovitis and treatment response in IA patients.

Methods

Clinical, laboratory and ultrasound data were retrieved from RA and PsA patients recruited in the prospective MATRIX cohort study. Inclusion criteria were age≥18 years, RA or PsA diagnosis, and presence of active disease with eligibility to start/modify therapy. Patients underwent ultrasound-guided synovial biopsy of one of the most involved joints before starting/modifying treatment according to treat-to-target strategy. The samples were analysed by an expert pathologist for KSS calculation. Univariable and multivariable logistic regression analyses were performed to evaluate the relationship between KSS and baseline variables. The association between KSS and treatment response at 24 weeks of follow-up was investigated in univariable logistic regression analysis.

Results

53 patients, 34 RA and 19 PsA, completed 24 weeks of follow-up after synovial biopsy. Patients were either treatment naïve (N=6, 11%), csDMARDs-experienced (N=46, 87%) or b/tsDMARDs-experienced (N=20, 38%). Median KSS was 6.00 (Q1-Q3 4.00-7.00) in RA and 4.00 (3.00-6.00) in PsA (p=0.040), and inflammatory infiltrates score was significantly higher in RA than in PsA patients (median 3.00 vs 2.00, p=0.021). In multivariable analysis, synovial effusion in the biopsied joint (OR 9.26, 95%CI 2.12-53.91) and erythrocyte sedimentation rate (ESR) (OR 1.04, 95%CI 1.01-1.08) associated with high KSS. High-grade synovitis significantly associated with a higher probability of achieving DAS28 remission, ACR20/50 response, and Boolean2.0 remission, independently from diagnosis.

Conclusion

Several markers of pro-inflammatory pathways associated with the presence of high-grade synovitis, and patients with higher KSS shared a higher probability of treatment targets achievement in the follow up. The integration of a simple and feasible tool like KSS in the clinical and prognostic stratification of patients with IA might help in intercepting patients with a disease more prone to respond to available treatment paradigms.

Metabolic reprogramming by Syntenin-1 directs RA FLS and endothelial cell-mediated inflammation and angiogenesis

A novel rheumatoid arthritis (RA) synovial fluid protein, Syntenin-1, and its receptor, Syndecan-1 (SDC-1), are colocalized on RA synovial tissue endothelial cells and fibroblast-like synoviocytes (FLS). Syntenin-1 exacerbates the inflammatory landscape of endothelial cells and RA FLS by upregulating transcription of IRF1/5/7/9, IL-1β, IL-6, and CCL2 through SDC-1 ligation and HIF1α, or mTOR activation. Mechanistically, Syntenin-1 orchestrates RA FLS and endothelial cell invasion via SDC-1 and/or mTOR signaling. In Syntenin-1 reprogrammed endothelial cells, the dynamic expression of metabolic intermediates coincides with escalated glycolysis along with unchanged oxidative factors, AMPK, PGC-1α, citrate, and inactive oxidative phosphorylation. Conversely, RA FLS rewired by Syntenin-1 displayed a modest glycolytic-ATP accompanied by a robust mitochondrial-ATP capacity. The enriched mitochondrial-ATP detected in Syntenin-1 reprogrammed RA FLS was coupled with mitochondrial fusion and fission recapitulated by escalated Mitofusin-2 and DRP1 expression. We found that VEGFR1/2 and Notch1 networks are responsible for the crosstalk between Syntenin-1 rewired endothelial cells and RA FLS, which are also represented in RA explants. Similar to RA explants, morphological and transcriptome studies authenticated the importance of VEGFR1/2, Notch1, RAPTOR, and HIF1α pathways in Syntenin-1 arthritic mice and their obstruction in SDC-1 deficient animals. Consistently, dysregulation of SDC-1, mTOR, and HIF1α negated Syntenin-1 inflammatory phenotype in RA explants, while inhibition of HIF1α impaired synovial angiogenic imprint amplified by Syntenin-1. In conclusion, since the current therapies are ineffective on Syntenin-1 and SDC-1 expression in RA synovial tissue and blood, targeting this pathway and its interconnected metabolic intermediates may provide a novel therapeutic strategy.

The long non-coding RNA HOTAIR contributes to joint-specific gene expression in rheumatoid arthritis

Although patients with rheumatoid arthritis (RA) typically exhibit symmetrical joint involvement, some patients develop alternative disease patterns in response to treatment, suggesting that different molecular mechanism may underlie disease progression depending on joint location. Here, we identify joint-specific changes in RA synovium and synovial fibroblasts (SF) between knee and hand joints. We show that the long non-coding RNA HOTAIR, which is only expressed in knee SF, regulates more than 50% of this site-specific gene expression in SF. HOTAIR is downregulated after stimulation with pro-inflammatory cytokines and is expressed at lower levels in knee samples from patients with RA, compared with osteoarthritis. Knockdown of HOTAIR in knee SF increases PI-Akt signalling and IL-6 production, but reduces Wnt signalling. Silencing HOTAIR inhibits the migratory function of SF, decreases SF-mediated osteoclastogenesis, and increases the recruitment of B cells by SF. We propose that HOTAIR is an important epigenetic factor in joint-specific gene expression in RA.

A Peripheral Blood Signature of Increased Th1 and Myeloid Cells Combined with Serum Inflammatory Mediators Is Associated with Response to Abatacept in Rheumatoid Arthritis Patients

Abatacept (CTLA4-Ig)-a monoclonal antibody which restricts T cell activation-is an effective treatment for rheumatoid arthritis (RA). Nevertheless, only 50% of RA patients attain clinical responses, while predictors of response are rather limited. Herein, we aimed to investigate for early biomarkers of response to abatacept, based on a detailed immunological profiling of peripheral blood (PB) cells and serum proteins. We applied flow cytometry and proteomics analysis on PB immune cells and serum respectively, of RA patients starting abatacept as the first biologic agent. After 6 months of treatment, 34.5% of patients attained response. At baseline, Th1 and FoxP3+ T cell populations were positively correlated with tender joint counts (p-value = 0.047 and p-value = 0.022, respectively). Upon treatment, CTLA4-Ig effectively reduced the percentages of Th1 and Th17 only in responders (p-value = 0.0277 and p-value = 0.0042, respectively). Notably, baseline levels of Th1 and myeloid cell populations were significantly increased in PB of responders compared to non-responders (p-value = 0.009 and p-value = 0.03, respectively). Proteomics analysis revealed that several inflammatory mediators were present in serum of responders before therapy initiation and strikingly 10 amongst 303 serum proteins were associated with clinical responses. Finally, a composite index based on selected baseline cellular and proteomics' analysis could predict response to abatacept with a high sensitivity (90%) and specificity (88.24%).

Disease activity drives transcriptomic heterogeneity in early untreated rheumatoid synovitis

OBJECTIVES

Transcriptomic profiling of synovial tissue from patients with early, untreated rheumatoid arthritis (RA) was used to explore the ability of unbiased, data-driven approaches to define clinically relevant subgroups.

METHODS

RNASeq was performed on 74 samples, with disease activity data collected at inclusion. Principal components analysis (PCA) and unsupervised clustering were used to define patient clusters based on expression of the most variable genes, followed by pathway analysis and inference of relative abundance of immune cell subsets. Histological assessment and multiplex immunofluorescence (for CD45, CD68, CD206) were performed on paraffin sections.

RESULTS

PCA on expression of the (n=894) most variable genes across this series did not divide samples into distinct groups, instead yielding a continuum correlated with baseline disease activity. Two patient clusters (PtC1, n=52; PtC2, n=22) were defined based on expression of these genes. PtC1, with significantly higher disease activity and probability of response to methotrexate therapy, showed upregulation of immune system genes; PtC2 showed upregulation of lipid metabolism genes, described to characterise tissue resident or M2-like macrophages. In keeping with these data, M2-like:M1-like macrophage ratios were inversely correlated with disease activity scores and were associated with lower synovial immune infiltration and the presence of thinner, M2-like macrophage-rich synovial lining layers.

CONCLUSION

In this large series of early, untreated RA, we show that the synovial transcriptome closely mirrors clinical disease activity and correlates with synovial inflammation. Intriguingly, lower inflammation and disease activity are associated with higher ratios of M2:M1 macrophages, particularly striking in the synovial lining layer. This may point to a protective role for tissue resident macrophages in RA.

Single-cell genomics reveals region-specific developmental trajectories underlying neuronal diversity in the human hypothalamus

The development and diversity of neuronal subtypes in the human hypothalamus has been insufficiently characterized. To address this, we integrated transcriptomic data from 241,096 cells (126,840 newly generated) in the prenatal and adult human hypothalamus to reveal a temporal trajectory from proliferative stem cell populations to mature hypothalamic cell types. Iterative clustering of the adult neurons identified 108 robust transcriptionally distinct neuronal subtypes representing 10 hypothalamic nuclei. Pseudotime trajectories provided insights into the genes driving formation of these nuclei. Comparisons to single-cell transcriptomic data from the mouse hypothalamus suggested extensive conservation of neuronal subtypes despite certain differences in species-enriched gene expression. The uniqueness of hypothalamic neuronal lineages was examined developmentally by comparing excitatory lineages present in cortex and inhibitory lineages in ganglionic eminence, revealing both distinct and shared drivers of neuronal maturation across the human forebrain. These results provide a comprehensive transcriptomic view of human hypothalamus development through gestation and adulthood at cellular resolution.

Donor white blood cell differential is the single largest determinant of whole blood gene expression patterns

It has become widely accepted that sample cellular composition is a significant determinant of the gene expression patterns observed in any transcriptomic experiment performed with bulk tissue. Despite this, many investigations currently performed with whole blood do not experimentally account for possible inter-specimen differences in cellularity, and often assume that any observed gene expression differences are a result of true differences in nuclear transcription. In order to determine how confounding of an assumption this may be, in this study, we recruited a large cohort of human donors (n = 138) and used a combination of next generation sequencing and flow cytometry to quantify and compare the underlying contributions of variance in leukocyte counts versus variance in other biological factors to overall variance in whole blood transcript levels. Our results suggest that the combination of donor neutrophil and lymphocyte counts alone are the primary determinants of whole blood transcript levels for up to 75% of the protein-coding genes expressed in peripheral circulation, whereas the other factors such as age, sex, race, ethnicity, and common disease states have comparatively minimal influence. Broadly, this infers that a majority of gene expression differences observed in experiments performed with whole blood are driven by latent differences in leukocyte counts, and that cell count heterogeneity must be accounted for to meaningfully biologically interpret the results.

Deconstruction of rheumatoid arthritis synovium defines inflammatory subtypes

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.

Vitamin B5 supports MYC oncogenic metabolism and tumor progression in breast cancer

Tumors are intrinsically heterogeneous and it is well established that this directs their evolution, hinders their classification and frustrates therapy1-3. Consequently, spatially resolved omics-level analyses are gaining traction4-9. Despite considerable therapeutic interest, tumor metabolism has been lagging behind this development and there is a paucity of data regarding its spatial organization. To address this shortcoming, we set out to study the local metabolic effects of the oncogene c-MYC, a pleiotropic transcription factor that accumulates with tumor progression and influences metabolism10,11. Through correlative mass spectrometry imaging, we show that pantothenic acid (vitamin B5) associates with MYC-high areas within both human and murine mammary tumors, where its conversion to coenzyme A fuels Krebs cycle activity. Mechanistically, we show that this is accomplished by MYC-mediated upregulation of its multivitamin transporter SLC5A6. Notably, we show that SLC5A6 over-expression alone can induce increased cell growth and a shift toward biosynthesis, whereas conversely, dietary restriction of pantothenic acid leads to a reversal of many MYC-mediated metabolic changes and results in hampered tumor growth. Our work thus establishes the availability of vitamins and cofactors as a potential bottleneck in tumor progression, which can be exploited therapeutically. Overall, we show that a spatial understanding of local metabolism facilitates the identification of clinically relevant, tractable metabolic targets.

Distinction and prognosis of early arthritis phenotypes: an analysis in three European cohorts

OBJECTIVES

The objective of this study is to evaluate whether there are differences in the long-term prognosis across various phenotypes of early arthritis (EA).

METHODS

Three EA cohorts (Reade, Etude et Suivi des Polyarthrites Indifférenciées Récentes (ESPOIR) and Early Arthritis Clinic (EAC)) were analysed. Clinical data were collected up to 24 years. Hands and feet radiographs were scored according to the Sharp van der Heijde (SvdH) method. Latent class analysis was applied to determine the EA phenotypes at baseline. Each class received a label reflecting its most prominent features. Prognostic outcomes included Health Assessment Questionnaire (HAQ), Short Form 36 (SF36) and SvdH score. The association between class membership and outcomes over time was tested in multivariable models.

RESULTS

In total, 390 (Reade), 798 (ESPOIR) and 3991 (EAC) patients were analysed separately. Two classes with symmetrical polyarthritis emerged; one of these labelled as autoimmune inflammatory polyarthritis (AIPA), had high likelihood of acute phase reactants (APR) elevation and autoantibody positivity, while the other (mild-inflammatory polyarthritis; MIPA) had not. A third class had oligoarthritis of upper limbs (OAUL) and could be subdivided into autoimmune OAUL and mild-inflammatory OAUL. A fifth class had oligoarthritis of lower limbs. The SvdH scores were worse in patients with APR/autoantibodies (AIPA) than in those without (MIPA). No clinically meaningful differences across classes in HAQ or SF36 over time were found.

CONCLUSION

Radiographic progression over time primarily occurs in EA patients with APR/autoantibodies. The absence of these markers, however, does not necessarily translate into better long-term function and quality of life. Clinicians should not only aim at preventing joint damage, but look beyond structural progression in order to further improve the lives of people with EA.

Stratification of biological therapies by pathobiology in biologic-naive patients with rheumatoid arthritis (STRAP and STRAP-EU): two parallel, open-label, biopsy-driven, randomised trials

BACKGROUND

Despite highly effective targeted therapies for rheumatoid arthritis, about 40% of patients respond poorly, and predictive biomarkers for treatment choices are lacking. We did a biopsy-driven trial to compare the response to rituximab, etanercept, and tocilizumab in biologic-naive patients with rheumatoid arthritis stratified for synovial B cell status.

METHODS

STRAP and STRAP-EU were two parallel, open-label, biopsy-driven, stratified, randomised, phase 3 trials done across 26 university centres in the UK and Europe. Biologic-naive patients aged 18 years or older with rheumatoid arthritis based on American College of Rheumatology (ACR)-European League Against Rheumatism classification criteria and an inadequate response to conventional synthetic disease-modifying antirheumatic drugs (DMARDs) were included. Following ultrasound-guided synovial biopsy, patients were classified as B cell poor or B cell rich according to synovial B cell signatures and randomly assigned (1:1:1) to intravenous rituximab (1000 mg at week 0 and week 2), subcutaneous tocilizumab (162 mg per week), or subcutaneous etanercept (50 mg per week). The primary outcome was the 16-week ACR20 response in the B cell-poor, intention-to-treat population (defined as all randomly assigned patients), with data pooled from the two trials, comparing etanercept and tocilizumab (grouped) versus rituximab. Safety was assessed in all patients who received at least one dose of study drug. These trials are registered with the EU Clinical Trials Register, 2014-003529-16 (STRAP) and 2017-004079-30 (STRAP-EU).

FINDINGS

Between June 8, 2015, and July 4, 2019, 226 patients were randomly assigned to etanercept (n=73), tocilizumab (n=74), and rituximab (n=79). Three patients (one in each group) were excluded after randomisation because they received parenteral steroids in the 4 weeks before recruitment. 168 (75%) of 223 patients in the intention-to-treat population were women and 170 (76%) were White. In the B cell-poor population, ACR20 response at 16 weeks (primary endpoint) showed no significant differences between etanercept and tocilizumab grouped together and rituximab (46 [60%] of 77 patients vs 26 [59%] of 44; odds ratio 1·02 [95% CI 0·47-2·17], p=0·97). No differences were observed for adverse events, including serious adverse events, which occurred in six (6%) of 102 patients in the rituximab group, nine (6%) of 108 patients in the etanercept group, and three (4%) of 73 patients in the tocilizumab group (p=0·53).

INTERPRETATION

In this biologic-naive population of patients with rheumatoid arthrtitis, the dichotomic classification into synovial B cell poor versus rich did not predict treatment response to B cell depletion with rituximab compared with alternative treatment strategies. However, the lack of response to rituximab in patients with a pauci-immune pathotype and the higher risk of structural damage progression in B cell-rich patients treated with rituximab warrant further investigations into the ability of synovial tissue analyses to inform disease pathogenesis and treatment response.

FUNDING

UK Medical Research Council and Versus Arthritis.

SARS-CoV-2 mRNA vaccine-induced immune responses in rheumatoid arthritis

Our objective was to characterize T and B cell responses to vaccination with SARS-CoV-2 antigens in immunocompromised rheumatoid arthritis (RA) patients. In 22 RA patients, clinical and biological variables were analyzed before and 4 weeks after each of 3 messenger RNA (mRNA) vaccine doses and compared with unmatched healthy individuals. Sequentially sampled peripheral blood mononuclear cells and sera were collected to determine immune profiles and to analyze the T cell response to a spike peptide pool and B cell specificity to the receptor-binding domain (RBD). Anti-spike antibodies were detectable in 6 of 22 RA patients after 1 dose of vaccine with increasing titers after each booster dose, although the overall response was lower compared with that in healthy control individuals. Responding patients after the first dose were more likely to have RA antibodies and a higher baseline proportion of circulating follicular B cells. In RA patients, the mRNA vaccine elicited a robust CD4+ T response to a spike peptide pool following the first and second doses. Consistent with the serologies, RBD-specific B cells exhibited a modest increase after the first dose and the second dose resulted in marked increases only in a fraction of the RA patients to both ancestral and omicron RBD. Our results highlight the importance of multidose COVID-19 vaccination in RA patients to develop a protective humoral response. However, these patients rapidly develop specific T CD4+ responses, despite delayed B cell responses.

Cytokine-directed cellular cross-talk imprints synovial pathotypes in rheumatoid arthritis

INTRODUCTION

Structural reorganisation of the synovium with expansion of fibroblast-like synoviocytes (FLS) and influx of immune cells is a hallmark of rheumatoid arthritis (RA). Activated FLS are increasingly recognised as a critical component driving synovial tissue remodelling by interacting with immune cells resulting in distinct synovial pathotypes of RA.

METHODS

Automated high-content fluorescence microscopy of co-cultured cytokine-activated FLS and autologous peripheral CD4+ T cells from patients with RA was established to quantify cell-cell interactions. Phenotypic profiling of cytokine-treated FLS and co-cultured T cells was done by flow cytometry and RNA-Seq, which were integrated with publicly available transcriptomic data from patients with different histological synovial pathotypes. Computational prediction and knock-down experiments were performed in FLS to identify adhesion molecules for cell-cell interaction.

RESULTS

Cytokine stimulation, especially with TNF-α, led to enhanced FLS-T cell interaction resulting in cell-cell contact-dependent activation, proliferation and differentiation of T cells. Signatures of cytokine-activated FLS were significantly enriched in RA synovial tissues defined as lymphoid-rich or leucocyte-rich pathotypes, with the most prominent effects for TNF-α. FLS cytokine signatures correlated with the number of infiltrating CD4+ T cells in synovial tissue of patients with RA. Ligand-receptor pair interaction analysis identified ICAM1 on FLS as an important mediator in TNF-mediated FLS-T cell interaction. Both, ICAM1 and its receptors were overexpressed in TNF-treated FLS and co-cultured T cells. Knock-down of ICAM1 in FLS resulted in reduced TNF-mediated FLS-T cell interaction.

CONCLUSION

Our study highlights the role of cytokine-activated FLS in orchestrating inflammation-associated synovial pathotypes providing novel insights into disease mechanisms of RA.

Rationale, current state and opportunities in combining biologic disease modifying antirheumatic drugs in rheumatoid and psoriatic arthritis

The advent of biologic disease modifying antirheumatic drugs (bDMARDs) has considerably improved patient outcomes in inflammatory arthritis. However, not all patients reach the state of remission, as disease can be resistant even to single cytokine inhibition by bDMARDs. Simultaneous or sequential inhibition of multiple cytokines may be considered in situations where disease control is not adequate under singular inhibition of cytokines. Although there have been some disappointing experiences in the past with combination of bDMARDs, the ongoing improvement of our understanding about inflammatory pathways and the overall better safety understanding of bDMARDs seem to make new biologic treatment combinations possible. This review covers the rationale and current evidence for bDMARDs combination in inflammatory arthritis.

Multi-omics analysis identifies IgG2b class-switching with ALCAM-CD6 co-stimulation in joint-draining lymph nodes during advanced inflammatory-erosive arthritis

Introduction

Defective lymphatic drainage and translocation of B-cells in inflamed (Bin) joint-draining lymph node sinuses are pathogenic phenomena in patients with severe rheumatoid arthritis (RA). However, the molecular mechanisms underlying this lymphatic dysfunction remain poorly understood. Herein, we utilized multi-omic spatial and single-cell transcriptomics to evaluate altered cellular composition (including lymphatic endothelial cells, macrophages, B-cells, and T-cells) in the joint-draining lymph node sinuses and their associated phenotypic changes and cell-cell interactions during RA development using the tumor necrosis factor transgenic (TNF-Tg) mouse model.

Methods

Popliteal lymph nodes (PLNs) from wild-type (n=10) and TNF-Tg male mice with "Early" (5 to 6-months of age; n=6) and "Advanced" (>8-months of age; n=12) arthritis were harvested and processed for spatial transcriptomics. Single-cell RNA sequencing (scRNAseq) was performed in PLNs from the TNF-Tg cohorts (n=6 PLNs pooled/cohort). PLN histopathology and ELISPOT along with ankle histology and micro-CT were evaluated. Histopathology of human lymph nodes and synovia was performed for clinical correlation.

Results

Advanced PLN sinuses exhibited an increased Ighg2b/Ighm expression ratio (Early 0.5 ± 0.1 vs Advanced 1.4 ± 0.5 counts/counts; p<0.001) that significantly correlated with reduced talus bone volumes in the afferent ankle (R2 = 0.54, p<0.001). Integration of single-cell and spatial transcriptomics revealed the increased IgG2b+ plasma cells localized in MARCO+ peri-follicular medullary sinuses. A concomitant decreased Fth1 expression (Early 2.5 ± 0.74 vs Advanced 1.0 ± 0.50 counts, p<0.001) within Advanced PLN sinuses was associated with accumulation of iron-laden Prussian blue positive macrophages in lymph nodes and synovium of Advanced TNF-Tg mice, and further validated in RA clinical samples. T-cells were increased 8-fold in Advanced PLNs, and bioinformatic pathway assessment identified the interaction between ALCAM+ macrophages and CD6+ T-cells as a plausible co-stimulatory mechanism to promote IgG2b class-switching.

Discussion

Collectively, these data support a model of flare in chronic TNF-induced arthritis in which loss of lymphatic flow through affected joint-draining lymph nodes facilitates the interaction between effluxing macrophages and T-cells via ALCAM-CD6 co-stimulation, initiating IgG2b class-switching and plasma cell differentiation of the expanded Bin population. Future work is warranted to investigate immunoglobulin clonality and potential autoimmune consequences, as well as the efficacy of anti-CD6 therapy to prevent these pathogenic events.

Machine Learning Reveals Synovial Fibroblast Genes Associated with Pain Affect Sensory Nerve Growth in Rheumatoid Arthritis

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.

Citrullinated and malondialdehyde-acetaldehyde modified fibrinogen activates macrophages and promotes an aggressive synovial fibroblast phenotype in patients with rheumatoid arthritis

Objective

Post-translational protein modifications with malondialdehyde-acetaldehyde (MAA) and citrulline (CIT) are implicated in the pathogenesis of rheumatoid arthritis (RA). Although precise mechanisms have not been elucidated, macrophage-fibroblast interactions have been proposed to play a central role in the development and progression of RA. The purpose of our study was to evaluate the downstream effects of macrophage released soluble mediators, following stimulation with fibrinogen (FIB) modified antigens, on human fibroblast-like synoviocytes (HFLS).

Methods

PMA-treated U-937 monocytes (Mϕ) and macrophage-differentiated peripheral blood mononuclear cells (MP) were stimulated with FIB, FIB-MAA, FIB-CIT, or FIB-MAA-CIT. HFLS-RA cells were stimulated directly with FIB antigens or with supernatants (SN) from macrophages (Mϕ-SN or MP-SN) stimulated with FIB antigens. Genes associated with an aggressive HFLS phenotype, extracellular matrix proteins, and activated signaling pathways were evaluated.

Results

HFLS-RA cells treated with Mϕ-SNFIB-CIT and Mϕ-SNFIB-MAA-CIT demonstrated significant increases in mRNA expression of genes associated with an aggressive phenotype at 24-h as compared to direct stimulation with the same antigens. Similar results were obtained using MP-SN. Cellular morphology was altered and protein expression of vimentin (p<0.0001 vs. Mϕ-SNFIB) and type II collagen (p<0.0001) were significantly increased in HFLS-RA cells treated with any of the Mϕ-SN generated following stimulation with modified antigens. Phosphorylation of JNK, Erk1/2, and Akt were increased most substantially in HFLS-RA treated with Mϕ-SNFIB-MAA-CIT (p<0.05 vs Mϕ-SNFIB). These and other data suggested the presence of PDGF-BB in Mϕ-SN. Mϕ-SNFIB-MAA-CIT contained the highest concentration of PDGF-BB (p<0.0001 vs. Mϕ-SNFIB) followed by Mϕ-SNFIB-CIT then Mϕ-SNFIB-MAA. HFLS-RA cells treated with PDGF-BB showed similar cellular morphology to the Mϕ-SN generated following stimulation with modified FIB, as well as the increased expression of vimentin, type II collagen, and the phosphorylation of JNK, Erk1/2 and Akt signaling molecules.

Conclusion

Together, these findings support the hypothesis that in response to MAA-modified and/or citrullinated fibrinogen, macrophages release soluble factors including PDGF-BB that induce fibroblast activation and promote an aggressive fibroblast phenotype. These cellular responses were most robust following macrophage activation with dually modified fibrinogen, compared to single modification alone, providing novel insights into the combined role of multiple post-translational protein modifications in the development of RA.

Dietary protein shapes the profile and repertoire of intestinal CD4+ T cells

The intestinal immune system must tolerate food antigens to avoid allergy, a process requiring CD4+ T cells. Combining antigenically defined diets with gnotobiotic models, we show that food and microbiota distinctly influence the profile and T cell receptor repertoire of intestinal CD4+ T cells. Independent of the microbiota, dietary proteins contributed to accumulation and clonal selection of antigen-experienced CD4+ T cells at the intestinal epithelium, imprinting a tissue-specialized transcriptional program including cytotoxic genes on both conventional and regulatory CD4+ T cells (Tregs). This steady state CD4+ T cell response to food was disrupted by inflammatory challenge, and protection against food allergy in this context was associated with Treg clonal expansion and decreased proinflammatory gene expression. Finally, we identified both steady-state epithelium-adapted CD4+ T cells and tolerance-induced Tregs that recognize dietary antigens, suggesting that both cell types may be critical for preventing inappropriate immune responses to food.

Computational approaches in rheumatic diseases - Deciphering complex spatio-temporal cell interactions

Inflammatory arthritis, including rheumatoid (RA), and psoriatic (PsA) arthritis, are clinically and immunologically heterogeneous diseases with no identified cure. Chronic inflammation of the synovial tissue ushers loss of function of the joint that severely impacts the patient's quality of life, eventually leading to disability and life-threatening comorbidities. The pathogenesis of synovial inflammation is the consequence of compounded immune and stromal cell interactions influenced by genetic and environmental factors. Deciphering the complexity of the synovial cellular landscape has accelerated primarily due to the utilisation of bulk and single cell RNA sequencing. Particularly the capacity to generate cell-cell interaction networks could reveal evidence of previously unappreciated processes leading to disease. However, there is currently a lack of universal nomenclature as a result of varied experimental and technological approaches that discombobulates the study of synovial inflammation. While spatial transcriptomic analysis that combines anatomical information with transcriptomic data of synovial tissue biopsies promises to provide more insights into disease pathogenesis, in vitro functional assays with single-cell resolution will be required to validate current bioinformatic applications. In order to provide a comprehensive approach and translate experimental data to clinical practice, a combination of clinical and molecular data with machine learning has the potential to enhance patient stratification and identify individuals at risk of arthritis that would benefit from early therapeutic intervention. This review aims to provide a comprehensive understanding of the effect of computational approaches in deciphering synovial inflammation pathogenesis and discuss the impact that further experimental and novel computational tools may have on therapeutic target identification and drug development.

Complement therapeutics are coming of age in rheumatology

The complement system was described over 100 years ago, and it is well established that activation of this pathway accompanies the great majority of autoimmune and inflammatory diseases. In addition, over three decades of work in murine models of human disease have nearly universally demonstrated that complement activation is upstream of tissue injury and the engagement of pro-inflammatory mechanisms such as the elaboration of cytokines and chemokines, as well as myeloid cell recruitment and activation. With that background, and taking advantage of advances in the development of biologic and small-molecule therapeutics, the creation and clinical evaluation of complement therapeutics is now rapidly expanding. This article reviews the current state of the complement therapeutics field, with a focus on their use in diseases cared for or consulted upon by rheumatologists. Included is an overview of the activation mechanisms and components of the system, in addition to the mechanisms by which the complement system interacts with other immune system constituents. The various therapeutic approaches to modulating the system in rheumatic and autoimmune diseases are reviewed. To understand how best to clinically assess the complement system, methods of its evaluation are described. Finally, next-generation therapeutic and diagnostic advances that can be envisioned for the future are discussed.

A novel 3D spheroid model of rheumatoid arthritis synovial tissue incorporating fibroblasts, endothelial cells, and macrophages

Objective

Rheumatoid Arthritis (RA) is a progressive and systemic autoimmune disorder associated with chronic and destructive joint inflammation. The hallmarks of joint synovial inflammation are cellular proliferation, extensive neoangiogenesis and infiltration of immune cells, including macrophages. In vitro approaches simulating RA synovial tissue are crucial in preclinical and translational research to evaluate novel diagnostic and/or therapeutic markers. Two-dimensional (2D) settings present very limited in vivo physiological proximity as they cannot recapitulate cell-cell and cell-matrix interactions occurring in the three-dimensional (3D) tissue compartment. Here, we present the engineering of a spheroid-based model of RA synovial tissue which mimics 3D interactions between cells and pro-inflammatory mediators present in the inflamed synovium.

Methods

Spheroids were generated by culturing RA fibroblast-like-synoviocytes (RAFLS), human umbilical vein endothelial cells (ECs) and monocyte-derived macrophages in a collagen-based 3D scaffold. The spheroids were cultured in the presence or absence of vascular endothelial growth factor (VEGF) and fibroblast growth factor 2 (bFGF) or RA synovial fluid (SF). Spheroid expansion and cell migration were quantified for all conditions using confocal microscopy and digital image analysis.

Results

A novel approach using machine learning was developed to quantify spheroid outgrowth and used to reexamine the existing spheroid-based model of RA synovial angiogenesis consisting of ECs and RAFLS. A 2-fold increase in the spheroid outgrowth ratio was demonstrated upon VEGF/bFGF stimulation (p<0.05). The addition of macrophages within the spheroid structure (3.75x104 RAFLS, 7.5x104 ECs and 3.0x104 macrophages) resulted in good incorporation of the new cell type. The addition of VEGF/bFGF significantly induced spheroid outgrowth (p<0.05) in the new system. SF stimulation enhanced containment of macrophages within the spheroids.

Conclusion

We present a novel spheroid based model consisting of RAFLS, ECs and macrophages that reflects the RA synovial tissue microenvironment. This model may be used to dissect the role of specific cell types in inflammatory responses in RA, to study specific signaling pathways involved in the disease pathogenesis and examine the effects of novel diagnostic (molecular imaging) and therapeutic compounds, including small molecule inhibitors and biologics.

DNA Methylation Signatures of Response to Conventional Synthetic and Biologic Disease-Modifying Antirheumatic Drugs (DMARDs) in Rheumatoid Arthritis

Rheumatoid arthritis (RA) is a complex condition that displays heterogeneity in disease severity and response to standard treatments between patients. Failure rates for conventional, target synthetic, and biologic disease-modifying rheumatic drugs (DMARDs) are significant. Although there are models for predicting patient response, they have limited accuracy, require replication/validation, or for samples to be obtained through a synovial biopsy. Thus, currently, there are no prediction methods approved for routine clinical use. Previous research has shown that genetics and environmental factors alone cannot explain the differences in response between patients. Recent studies have demonstrated that deoxyribonucleic acid (DNA) methylation plays an important role in the pathogenesis and disease progression of RA. Importantly, specific DNA methylation profiles associated with response to conventional, target synthetic, and biologic DMARDs have been found in the blood of RA patients and could potentially function as predictive biomarkers. This review will summarize and evaluate the evidence for DNA methylation signatures in treatment response mainly in blood but also learn from the progress made in the diseased tissue in cancer in comparison to RA and autoimmune diseases. We will discuss the benefits and challenges of using DNA methylation signatures as predictive markers and the potential for future progress in this area.

BBmix: a Bayesian beta-binomial mixture model for accurate genotyping from RNA-sequencing

MOTIVATION

While many pipelines have been developed for calling genotypes using RNA-sequencing (RNA-Seq) data, they all have adapted DNA genotype callers that do not model biases specific to RNA-Seq such as allele-specific expression (ASE).

RESULTS

Here, we present Bayesian beta-binomial mixture model (BBmix), a Bayesian beta-binomial mixture model that first learns the expected distribution of read counts for each genotype, and then deploys those learned parameters to call genotypes probabilistically. We benchmarked our model on a wide variety of datasets and showed that our method generally performed better than competitors, mainly due to an increase of up to 1.4% in the accuracy of heterozygous calls, which may have a big impact in reducing false positive rate in applications sensitive to genotyping error such as ASE. Moreover, BBmix can be easily incorporated into standard pipelines for calling genotypes. We further show that parameters are generally transferable within datasets, such that a single learning run of less than 1 h is sufficient to call genotypes in a large number of samples.

AVAILABILITY AND IMPLEMENTATION

We implemented BBmix as an R package that is available for free under a GPL-2 licence at https://gitlab.com/evigorito/bbmix and https://cran.r-project.org/package=bbmix with accompanying pipeline at https://gitlab.com/evigorito/bbmix_pipeline.