Identification of gene regulatory networks affected across drug-resistant epilepsies.

Epilepsy is a chronic and heterogenous disease characterized by recurrent unprovoked seizures, that are commonly resistant to antiseizure medications. This study applies a transcriptome network-based approach across epilepsies aiming to improve understanding of molecular disease pathobiology, recognize affected biological mechanisms and apply causal reasoning to identify therapeutic hypotheses. This study included the most common drug-resistant epilepsies (DREs), such as temporal lobe epilepsy with hippocampal sclerosis (TLE-HS), and mTOR pathway-related malformations of cortical development (mTORopathies). This systematic comparison characterized the global molecular signature of epilepsies, elucidating the key underlying mechanisms of disease pathology including neurotransmission and synaptic plasticity, brain extracellular matrix and energy metabolism. In addition, specific dysregulations in neuroinflammation and oligodendrocyte function were observed in TLE-HS and mTORopathies, respectively. The aforementioned mechanisms are proposed as molecular hallmarks of DRE with the identified upstream regulators offering opportunities for drug-target discovery and development.

Thioredoxin is a metabolic rheostat controlling regulatory B cells.

Metabolic programming is important for B cell fate, but the bioenergetic requirement for regulatory B (Breg) cell differentiation and function is unknown. Here we show that Breg cell differentiation, unlike non-Breg cells, relies on mitochondrial electron transport and homeostatic levels of reactive oxygen species (ROS). Single-cell RNA sequencing analysis revealed that TXN, encoding the metabolic redox protein thioredoxin (Trx), is highly expressed by Breg cells, unlike Trx inhibitor TXNIP which was downregulated. Pharmacological inhibition or gene silencing of TXN resulted in mitochondrial membrane depolarization and increased ROS levels, selectively suppressing Breg cell differentiation and function while favoring pro-inflammatory B cell differentiation. Patients with systemic lupus erythematosus (SLE), characterized by Breg cell deficiencies, present with B cell mitochondrial membrane depolarization, elevated ROS and fewer Trx+ B cells. Exogenous Trx stimulation restored Breg cells and mitochondrial membrane polarization in SLE B cells to healthy B cell levels, indicating Trx insufficiency underlies Breg cell impairment in patients with SLE.

Differential regulation of IL-17A and IL-17F via STAT5 contributes to psoriatic disease.

BACKGROUND: IL-17A plays a pivotal pathogenic role in several immune-mediated inflammatory diseases. Despite sharing 50% sequence homology with IL-17A, the role of IL-17F remains less clear. Clinical findings suggest that dual inhibition of IL-17A and IL-17F in psoriatic disease is more efficacious than IL-17A inhibition alone, positing a pathogenic role for IL-17F. OBJECTIVE: We characterized the regulation of IL-17A and IL-17F in psoriatic disease. METHODS: Using both in vitro systems and lesional skin tissue from patients, we interrogated the chromosomal, transcriptional, and protein expression landscape of IL-17A+ and IL-17F+ TH17 cells. Alongside established assays such as single-cell RNA sequencing, we developed a novel cytokine-capture technique that was combined with chromatin immunoprecipitation sequencing and RNA sequencing. RESULTS: We confirm a preferential elevation of IL-17F over IL-17A in psoriatic disease and show that expression of each isoform predominantly occurs in distinct cell populations. The expression of both IL-17A and IL-17F exhibited a high degree of plasticity, with the balance between the 2 isoforms influenced by proinflammatory signaling and by anti-inflammatory drugs such as methylprednisolone. This plasticity was reflected in a broad H3K4me3 region at the IL17A-F locus, while opposing effects of STAT5/IL-2 signaling were observed for each of the 2 genes. Functionally, higher IL17F expression was linked to greater cell proliferation. CONCLUSION: There are key differences in the regulation of IL-17A and IL-17F in psoriatic disease, leading to distinct inflammatory cell populations. As such, we propose that both IL-17A and IL-17F neutralization may be required to maximally inhibit IL-17-driven pathology.

DNA methylation analysis identifies key transcription factors involved in mesenchymal stem cell osteogenic differentiation.

BACKGROUND: Knowledge about regulating transcription factors (TFs) for osteoblastogenesis from mesenchymal stem cells (MSCs) is limited. Therefore, we investigated the relationship between genomic regions subject to DNA-methylation changes during osteoblastogenesis and the TFs known to directly interact with these regulatory regions. RESULTS: The genome-wide DNA-methylation signature of MSCs differentiated to osteoblasts and adipocytes was determined using the Illumina HumanMethylation450 BeadChip array. During adipogenesis no CpGs passed our test for significant methylation changes. Oppositely, during osteoblastogenesis we identified 2462 differently significantly methylated CpGs (adj. p < 0.05). These resided outside of CpGs islands and were significantly enriched in enhancer regions. We confirmed the correlation between DNA-methylation and gene expression. Accordingly, we developed a bioinformatic tool to analyse differentially methylated regions and the TFs interacting with them. By overlaying our osteoblastogenesis differentially methylated regions with ENCODE TF ChIP-seq data we obtained a set of candidate TFs associated to DNA-methylation changes. Among them, ZEB1 TF was highly related with DNA-methylation. Using RNA interference, we confirmed that ZEB1, and ZEB2, played a key role in adipogenesis and osteoblastogenesis processes. For clinical relevance, ZEB1 mRNA expression in human bone samples was evaluated. This expression positively correlated with weight, body mass index, and PPARγ expression. CONCLUSIONS: In this work we describe an osteoblastogenesis-associated DNA-methylation profile and, using these data, validate a novel computational tool to identify key TFs associated to age-related disease processes. By means of this tool we identified and confirmed ZEB TFs as mediators involved in the MSCs differentiation to osteoblasts and adipocytes, and obesity-related bone adiposity.

DNA methylation analysis identifies key transcription factors involved in mesenchymal stem cell osteogenic differentiation

BACKGROUND: Knowledge about regulating transcription factors (TFs) for osteoblastogenesis from mesenchymal stem cells (MSCs) is limited. Therefore, we investigated the relationship between genomic regions subject to DNA-methylation changes during osteoblastogenesis and the TFs known to directly interact with these regulatory regions. RESULTS: The genome-wide DNA-methylation signature of MSCs differentiated to osteoblasts and adipocytes was determined using the Illumina HumanMethylation450 BeadChip array. During adipogenesis no CpGs passed our test for significant methylation changes. Oppositely, during osteoblastogenesis we identified 2462 differently significantly methylated CpGs (adj. p < 0.05). These resided outside of CpGs islands and were significantly enriched in enhancer regions. We confirmed the correlation between DNA-methylation and gene expression. Accordingly, we developed a bioinformatic tool to analyse differentially methylated regions and the TFs interacting with them. By overlaying our osteoblastogenesis differentially methylated regions with ENCODE TF ChIP-seq data we obtained a set of candidate TFs associated to DNA-methylation changes. Among them, ZEB1 TF was highly related with DNA-methylation. Using RNA interference, we confirmed that ZEB1, and ZEB2, played a key role in adipogenesis and osteoblastogenesis processes. For clinical relevance, ZEB1 mRNA expression in human bone samples was evaluated. This expression positively correlated with weight, body mass index, and PPARγ expression. CONCLUSIONS: In this work we describe an osteoblastogenesis-associated DNA-methylation profile and, using these data, validate a novel computational tool to identify key TFs associated to age-related disease processes. By means of this tool we identified and confirmed ZEB TFs as mediators involved in the MSCs differentiation to osteoblasts and adipocytes, and obesity-related bone adiposity.

Regulation of microRNA-221, -222, -21 and -27 in articular cartilage subjected to abnormal compressive forces.

KEY POINTS: microRNAs (miRs) are small non-coding molecules that regulate post-transcriptional target gene expression. miRs are involved in regulating cellular activities in response to mechanical loading in all physiological systems, although it is largely unknown whether this response differs with increasing magnitudes of load. miR-221, miR-222, miR-21-5p and miR-27a-5p were significantly increased in ex vivo cartilage explants subjected to increasing load magnitude and in in vivo joint cartilage exposed to abnormal loading. TIMP3 and CPEB3 are putative miR targets in chondrocytes Identification of mechanically regulated miRs that have potential to impact on tissue homeostasis provides a mechanism by which load-induced tissue behaviour is regulated, in both health and pathology, in all physiological systems. ABSTRACT: MicroRNAs (miRs) are small non-coding molecules that regulate post-transcriptional target gene expression and are involved in mechano-regulation of cellular activities in all physiological systems. It is unknown whether such epigenetic mechanisms are regulated in response to increasing magnitudes of load. The present study investigated mechano-regulation of miRs in articular cartilage subjected to 'physiological' and 'non-physiological' compressive loads in vitro as a model system and validated findings in an in vivo model of abnormal joint loading. Bovine full-depth articular cartilage explants were loaded to 2.5 MPa (physiological) or 7 MPa (non-physiological) (1 Hz, 15 min) and mechanically-regulated miRs identified using next generation sequencing and verified using a quantitative PCR. Downstream targets were verified using miR-specific mimics or inhibitors in conjunction with 3'-UTR luciferase activity assays. A subset of miRs were mechanically-regulated in ex vivo cartilage explants and in vivo joint cartilage. miR-221, miR-222, miR-21-5p and miR-27a-5p were increased and miR-483 levels decreased with increasing load magnitude. Tissue inhibitor of metalloproteinase 3 (TIMP3) and cytoplasmic polyadenylation element binding protein 3 (CPEB3) were identified as putative downstream targets. Our data confirm miR-221 and -222 mechano-regulation and demonstrates novel mechano-regulation of miR-21-5p and miR-27a-5p in ex vivo and in vivo cartilage loading models. TIMP3 and CPEB3 are putative miR targets in chondrocytes. Identification of specific miRs that are regulated by increasing load magnitude, as well as their potential to impact on tissue homeostasis, has direct relevance to other mechano-sensitive physiological systems and provides a mechanism by which load-induced tissue behaviour is regulated, in both health and pathology.

Early-onset autoimmunity associated with SOCS1 haploinsufficiency.

Autoimmunity can occur when a checkpoint of self-tolerance fails. The study of familial autoimmune diseases can reveal pathophysiological mechanisms involved in more common autoimmune diseases. Here, by whole-exome/genome sequencing we identify heterozygous, autosomal-dominant, germline loss-of-function mutations in the SOCS1 gene in ten patients from five unrelated families with early onset autoimmune manifestations. The intracellular protein SOCS1 is known to downregulate cytokine signaling by inhibiting the JAK-STAT pathway. Accordingly, patient-derived lymphocytes exhibit increased STAT activation in vitro in response to interferon-γ, IL-2 and IL-4 that is reverted by the JAK1/JAK2 inhibitor ruxolitinib. This effect is associated with a series of in vitro and in vivo immune abnormalities consistent with lymphocyte hyperactivity. Hence, SOCS1 haploinsufficiency causes a dominantly inherited predisposition to early onset autoimmune diseases related to cytokine hypersensitivity of immune cells.

microRNA-seq of cartilage reveals an overabundance of miR-140-3p which contains functional isomiRs.

miR-140 is selectively expressed in cartilage. Deletion of the entire Mir140 locus in mice results in growth retardation and early-onset osteoarthritis-like pathology; however, the relative contribution of miR-140-5p or miR-140-3p to the phenotype remains to be determined. An unbiased small RNA sequencing approach identified miR-140-3p as significantly more abundant (>10-fold) than miR-140-5p in human cartilage. Analysis of these data identified multiple miR-140-3p isomiRs differing from the miRBase annotation at both the 5' and 3' end, with >99% having one of two seed sequences (5' bases 2-8). Canonical (miR-140-3p.2) and shifted (miR-140-3p.1) seed isomiRs were overexpressed in chondrocytes and transcriptomics performed to identify targets. miR-140-3p.1 and miR-140-3p.2 significantly down-regulated 694 and 238 genes, respectively, of which only 162 genes were commonly down-regulated. IsomiR targets were validated using 3'UTR luciferase assays. miR-140-3p.1 targets were enriched within up-regulated genes in rib chondrocytes of Mir140-null mice and within down-regulated genes during human chondrogenesis. Finally, through imputing the expression of miR-140 from the expression of the host gene WWP2 in 124 previously published data sets, an inverse correlation with miR-140-3p.1 predicted targets was identified. Together these data suggest the novel seed containing isomiR miR-140-3p.1 is more functional than original consensus miR-140-3p seed containing isomiR.

Isolation and next generation sequencing of archival formalin-fixed DNA.

DNA from archived organs is presumed unsuitable for genomic studies because of excessive formalin-fixation. As next generation sequencing (NGS) requires short DNA fragments, and Uracil-N-glycosylase (UNG) can be used to overcome deamination, there has been renewed interest in the possibility of genomic studies using these collections. We describe a novel method of DNA extraction capable of providing PCR amplicons of at least 400 bp length from such excessively formalin-fixed human tissues. When compared with a leading commercial formalin-fixed DNA extraction kit, our method produced greater yields of DNA and reduced sequence variations. Analysis of PCR products using bacterial sub-cloning and Sanger sequencing from UNG-treated DNA unexpectedly revealed increased sequence variations, compared with untreated samples. Finally, whole exome NGS was performed on a myocardial sample fixed in formalin for 2 years and compared with lymphocyte-derived DNA (as a gold standard) from the same patient. Despite the reduction in the number and quality of reads in the formalin-fixed DNA, we were able to show that bioinformatic processing by joint calling and variant quality score recalibration (VQSR) increased the sensitivity four-fold to 56% and doubled specificity to 68% when compared with a standard hard-filtering approach. Thus, high-quality DNA can be extracted from excessively formalin-fixed tissues and bioinformatic processing can optimise sensitivity and specificity of results. Sequencing of several sub-cloned amplicons is an important methodological step in assessing DNA quality.

Lymphocyte DNA methylation mediates genetic risk at shared immune-mediated disease loci.

BACKGROUND: Defining regulatory mechanisms through which noncoding risk variants influence the cell-mediated pathogenesis of immune-mediated disease (IMD) has emerged as a priority in the post-genome-wide association study era. OBJECTIVES: With a focus on rheumatoid arthritis, we sought new insight into genetic mechanisms of adaptive immune dysregulation to help prioritize molecular pathways for targeting in this and related immune pathologies. METHODS: Whole-genome methylation and transcriptional data from isolated CD4+ T cells and B cells of more than 100 genotyped and phenotyped patients with inflammatory arthritis, all of whom were naive to immunomodulatory treatments, were obtained. Analysis integrated these comprehensive data with genome-wide association study findings across IMDs and other publicly available resources. RESULTS: We provide strong evidence that disease-associated DNA variants regulate cis-CpG methylation in CD4+ T and/or B cells at 37% RA loci. Using paired, cell-specific transcriptomic data and causal inference testing, we identify examples where site-specific DNA methylation in turn mediates gene expression, including FCRL3 in both cell types and ORMDL3/GSDMB, IL6ST/ANKRD55, and JAZF1 in CD4+ T cells. A number of genes regulated in this way highlight mechanisms common to RA and other IMDs including multiple sclerosis and asthma, in turn distinguishing them from osteoarthritis, a primarily degenerative disease. Finally, we corroborate the observed effects experimentally. CONCLUSIONS: Our observations highlight important mechanisms of genetic risk in RA and the wider context of immune dysregulation. They confirm the utility of DNA methylation profiling as a tool for causal gene prioritization and, potentially, therapeutic targeting in complex IMD.

DNA hypomethylation during MSC chondrogenesis occurs predominantly at enhancer regions.

Regulation of transcription occurs in a cell type specific manner orchestrated by epigenetic mechanisms including DNA methylation. Methylation changes may also play a key role in lineage specification during stem cell differentiation. To further our understanding of epigenetic regulation in chondrocytes we characterised the DNA methylation changes during chondrogenesis of mesenchymal stem cells (MSCs) by Infinium 450 K methylation array. Significant DNA hypomethylation was identified during chondrogenic differentiation including changes at many key cartilage gene loci. Integration with chondrogenesis gene expression data revealed an enrichment of significant CpGs in upregulated genes, while characterisation of significant CpG loci indicated their predominant localisation to enhancer regions. Comparison with methylation profiles of other tissues, including healthy and diseased adult cartilage, identified chondrocyte-specific regions of hypomethylation and the overlap with differentially methylated CpGs in osteoarthritis. Taken together we have associated DNA methylation levels with the chondrocyte phenotype. The consequences of which has potential to improve cartilage generation for tissue engineering purposes and also to provide context for observed methylation changes in cartilage diseases such as osteoarthritis.

Severe type I interferonopathy and unrestrained interferon signaling due to a homozygous germline mutation in STAT2.

Excessive type I interferon (IFNα/ß) activity is implicated in a spectrum of human disease, yet its direct role remains to be conclusively proven. We investigated two siblings with severe early-onset autoinflammatory disease and an elevated IFN signature. Whole-exome sequencing revealed a shared homozygous missense Arg148Trp variant in STAT2, a transcription factor that functions exclusively downstream of innate IFNs. Cells bearing STAT2R148W in homozygosity (but not heterozygosity) were hypersensitive to IFNα/ß, which manifest as prolonged Janus kinase-signal transducers and activators of transcription (STAT) signaling and transcriptional activation. We show that this gain of IFN activity results from the failure of mutant STAT2R148W to interact with ubiquitin-specific protease 18, a key STAT2-dependent negative regulator of IFNα/ß signaling. These observations reveal an essential in vivo function of STAT2 in the regulation of human IFNα/ß signaling, providing concrete evidence of the serious pathological consequences of unrestrained IFNα/ß activity and supporting efforts to target this pathway therapeutically in IFN-associated disease.

RNASeq_similarity_matrix: visually identify sample mix-ups in RNASeq data using a 'genomic' sequence similarity matrix.

MOTIVATION: Mistakes in linking a patient's biological samples with their phenotype data can confound RNA-Seq studies. The current method for avoiding such sample mixups is to test for inconsistencies between biological data and known phenotype data such as sex. However, in DNA studies a common QC step is to check for unexpected relatedness between samples. Here, we extend this method to RNA-Seq, which allows the detection of duplicated samples without relying on identifying inconsistencies with phenotype data. SUMMARY: We present RNASeq_similarity_matrix: an automated tool to generate a sequence similarity matrix from RNA-Seq data, which can be used to visually identify sample mix-ups. This is particularly useful when a study contains multiple samples from the same individual, but can also detect contamination in studies with only one sample per individual. AVAILABILITY: RNASeq_similarity_matrix has been made available as a documented GPL licensed Docker image on www.github.com/nicokist/RNASeq_similarity_matrix. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Predicting drug-free remission in rheumatoid arthritis: A prospective interventional cohort study.

BACKGROUND: Many patients with rheumatoid arthritis (RA) achieve disease remission with modern treatment strategies. However, having achieved this state, there are no tests that predict when withdrawal of therapy will result in drug-free remission rather than flare. We aimed to identify predictors of drug-free remission in RA. METHODS: The Biomarkers of Remission in Rheumatoid Arthritis (BioRRA) Study was a unique, prospective, interventional cohort study of complete and abrupt cessation of conventional synthetic disease-modifying anti-rheumatic drugs (DMARDs). Patients with RA of at least 12 months duration and in clinical and ultrasound remission discontinued DMARDs and were monitored for six months. The primary outcome was time-to-flare, defined as disease activity score in 28 joints with C-reactive protein (DAS28-CRP) ≥ 2.4. Baseline clinical and ultrasound measures, circulating inflammatory biomarkers, and peripheral CD4+ T cell gene expression were assessed for their ability to predict time-to-flare and flare/remission status by Cox regression and receiver-operating characteristic (ROC) analysis respectively. RESULTS: 23/44 (52%) eligible patients experienced an arthritis flare after a median (IQR) of 48 (31.5-86.5) days following DMARD cessation. A composite score incorporating five baseline variables (three transcripts [FAM102B, ENSG00000228010, ENSG00000227070], one cytokine [interleukin-27], one clinical [Boolean remission]) differentiated future flare from drug-free remission with an area under the ROC curve of 0.96 (95% CI 0.91-1.00), sensitivity 0.91 (0.78-1.00) and specificity 0.95 (0.84-1.00). CONCLUSION: We provide proof-of-concept evidence for predictors of drug-free remission in RA. If validated, these biomarkers could help to personalize immunosuppressant withdrawal: a therapy paradigm shift with ensuing patient and economic benefits.

IL-6 Mediated Transcriptional Programming of Naïve CD4+ T Cells in Early Rheumatoid Arthritis Drives Dysregulated Effector Function.

Objective: We have previously shown that increased circulating interleukin-6 (IL-6) results in enhanced CD4+ T cell signaling via signal transduction and activator of transcription-3 (STAT3) in early rheumatoid arthritis (RA). We tested the hypothesis that transcriptional "imprinting" of T-cells by this mechanism skews downstream effector responses, reinforcing immune dysregulation at a critical, but targetable, disease phase. Methods: We modeled naïve CD4+ T cell exposure to pathophysiological concentrations of IL-6 in vitro, assessing the dynamic transcriptional and functional consequences for downstream effector cells utilizing microarray and flow cytometry. Fresh blood from treatment-naïve early arthritis patients was phenotyped in parallel for comparison. Results: T cell sensitivity to IL-6 was most marked in the naïve subset, and related to gp130 rather than IL-6R expression. Exposure of healthy naïve CD4+ T cells to IL-6 induced the same STAT3 target genes as previously seen to discriminate RA patients from disease controls. After TCR stimulation IL-6 pre-exposed cells exhibited enhanced proliferative capacity, activation, and a propensity toward Th1 differentiation, compared to non-exposed cells. An entirely analogous phenotype was observed in early RA compared to control CD4+ T cells. Conclusions: Sustained IL-6 exposure at a critical point in the natural history of RA "primes" the adaptive immune system to respond aberrantly to TCR stimulation, potentiating disease induction with implications for the optimal timing of targeted therapy.

Prioritization of PLEC and GRINA as Osteoarthritis Risk Genes Through the Identification and Characterization of Novel Methylation Quantitative Trait Loci.

OBJECTIVE: To identify methylation quantitative trait loci (mQTLs) correlating with osteoarthritis (OA) risk alleles and to undertake mechanistic characterization as a means of target gene prioritization. METHODS: We used genome-wide genotyping and cartilage DNA methylation array data in a discovery screen of novel OA risk loci. This was followed by methylation, gene expression analysis, and genotyping studies in additional cartilage samples, accompanied by in silico analyses. RESULTS: We identified 4 novel OA mQTLs. The most significant mQTL contained 9 CpG sites where methylation correlated with OA risk genotype, with 5 of the CpG sites having P values <1 × 10-10 . The 9 CpG sites reside in an interval of only 7.7 kb within the PLEC gene and form 2 distinct clusters. We were able to prioritize PLEC and the adjacent gene GRINA as independent targets of the OA risk. We identified PLEC and GRINA expression QTLs operating in cartilage, as well as methylation-expression QTLs operating on the 2 genes. GRINA and PLEC also demonstrated differential expression between OA hip and non-OA hip cartilage. CONCLUSION: PLEC encodes plectin, a cytoskeletal protein that maintains tissue integrity by regulating intracellular signaling in response to mechanical stimuli. GRINA encodes the ionotropic glutamate receptor TMBIM3 (transmembrane BAX inhibitor 1 motif-containing protein family member 3), which regulates cell survival. Based on our results, we hypothesize that in a joint predisposed to OA, expression of these genes alters in order to combat aberrant biomechanics, and that this is epigenetically regulated. However, carriage of the OA risk-conferring allele at this locus hinders this response and contributes to disease development.

Expression of STAT3-regulated genes in circulating CD4+ T cells discriminates rheumatoid arthritis independently of clinical parameters in early arthritis.

OBJECTIVES: Dysregulated signal transduction and activator of transcription-3 (STAT3) signalling in CD4+ T cells has been proposed as an early pathophysiological event in RA. We sought further evidence for this observation, and to determine its clinical relevance. METHODS: Microarray technology was used to measure gene expression in purified peripheral blood CD4+ T cells from treatment-naïve RA patients and disease controls newly recruited from an early arthritis clinic. Analysis focused on 12 previously proposed transcripts, and concurrent STAT3 pathway activation was determined in the same cells by flow cytometry. A pooled analysis of previous and current gene expression findings incorporated detailed clinical parameters and employed multivariate analysis. RESULTS: In an independent cohort of 161 patients, expression of 11 of 12 proposed signature genes differed significantly between RA patients and controls, robustly validating the earlier findings. Differential regulation was most pronounced for the STAT3 target genes PIM1, BCL3 and SOCS3 (>1.3-fold difference; P < 0.005), each of whose expression correlated strongly with paired intracellular phospho-STAT3. In a meta-analysis of 279 patients the same three genes accounted for the majority of the signature's ability to discriminate RA patients, which was found to be independent of age, joint involvement or acute phase response. CONCLUSION: The STAT3-mediated dysregulation of BCL3, SOCS3 and PIM1 in circulating CD4+ T cells is a discriminatory feature of early RA that occurs independently of acute phase response. The mechanistic and functional implications of this observation at a cellular level warrant clarification.

Symptom-based stratification of patients with primary Sjögren's syndrome: multi-dimensional characterisation of international observational cohorts and reanalyses of randomised clinical trials.

BACKGROUND: Heterogeneity is a major obstacle to developing effective treatments for patients with primary Sjögren's syndrome. We aimed to develop a robust method for stratification, exploiting heterogeneity in patient-reported symptoms, and to relate these differences to pathobiology and therapeutic response. METHODS: We did hierarchical cluster analysis using five common symptoms associated with primary Sjögren's syndrome (pain, fatigue, dryness, anxiety, and depression), followed by multinomial logistic regression to identify subgroups in the UK Primary Sjögren's Syndrome Registry (UKPSSR). We assessed clinical and biological differences between these subgroups, including transcriptional differences in peripheral blood. Patients from two independent validation cohorts in Norway and France were used to confirm patient stratification. Data from two phase 3 clinical trials were similarly stratified to assess the differences between subgroups in treatment response to hydroxychloroquine and rituximab. FINDINGS: In the UKPSSR cohort (n=608), we identified four subgroups: Low symptom burden (LSB), high symptom burden (HSB), dryness dominant with fatigue (DDF), and pain dominant with fatigue (PDF). Significant differences in peripheral blood lymphocyte counts, anti-SSA and anti-SSB antibody positivity, as well as serum IgG, κ-free light chain, ß2-microglobulin, and CXCL13 concentrations were observed between these subgroups, along with differentially expressed transcriptomic modules in peripheral blood. Similar findings were observed in the independent validation cohorts (n=396). Reanalysis of trial data stratifying patients into these subgroups suggested a treatment effect with hydroxychloroquine in the HSB subgroup and with rituximab in the DDF subgroup compared with placebo. INTERPRETATION: Stratification on the basis of patient-reported symptoms of patients with primary Sjögren's syndrome revealed distinct pathobiological endotypes with distinct responses to immunomodulatory treatments. Our data have important implications for clinical management, trial design, and therapeutic development. Similar stratification approaches might be useful for patients with other chronic immune-mediated diseases. FUNDING: UK Medical Research Council, British Sjogren's Syndrome Association, French Ministry of Health, Arthritis Research UK, Foundation for Research in Rheumatology. VIDEO ABSTRACT.