Coordinated glucocorticoid receptor and MAFB action induces tolerogenesis and epigenome remodeling in dendritic cells.

Glucocorticoids (GCs) exert potent anti-inflammatory effects in immune cells through the glucocorticoid receptor (GR). Dendritic cells (DCs), central actors for coordinating immune responses, acquire tolerogenic properties in response to GCs. Tolerogenic DCs (tolDCs) have emerged as a potential treatment for various inflammatory diseases. To date, the underlying cell type-specific regulatory mechanisms orchestrating GC-mediated acquisition of immunosuppressive properties remain poorly understood. In this study, we investigated the transcriptomic and epigenomic remodeling associated with differentiation to DCs in the presence of GCs. Our analysis demonstrates a major role of MAFB in this process, in synergy with GR. GR and MAFB both interact with methylcytosine dioxygenase TET2 and bind to genomic loci that undergo specific demethylation in tolDCs. We also show that the role of MAFB is more extensive, binding to thousands of genomic loci in tolDCs. Finally, MAFB knockdown erases the tolerogenic properties of tolDCs and reverts the specific DNA demethylation and gene upregulation. The preeminent role of MAFB is also demonstrated in vivo for myeloid cells from synovium in rheumatoid arthritis following GC treatment. Our results imply that, once directly activated by GR, MAFB plays a critical role in orchestrating the epigenomic and transcriptomic remodeling that define the tolerogenic phenotype.

EOMES-positive CD4+ T cells are increased in PTPN22 (1858T) risk allele carriers.

The presence of the PTPN22 risk allele (1858T) is associated with several autoimmune diseases including rheumatoid arthritis (RA). Despite a number of studies exploring the function of PTPN22 in T cells, the exact impact of the PTPN22 risk allele on T-cell function in humans is still unclear. In this study, using RNA sequencing, we show that, upon TCR-activation, naïve human CD4+ T cells homozygous for the PTPN22 risk allele overexpress a set of genes including CFLAR and 4-1BB, which are important for cytotoxic T-cell differentiation. Moreover, the protein expression of the T-box transcription factor Eomesodermin (EOMES) was increased in T cells from healthy donors homozygous for the PTPN22 risk allele and correlated with a decreased number of naïve CD4+ T cells. There was no difference in the frequency of other CD4+ T-cell subsets (Th1, Th17, Tfh, Treg). Finally, an accumulation of EOMES+ CD4+ T cells was observed in synovial fluid of RA patients with a more pronounced production of Perforin-1 in PTPN22 risk allele carriers. Altogether, we propose a novel mechanism of action of PTPN22 risk allele through the generation of cytotoxic CD4+ T cells and identify EOMES+ CD4+ T cells as a relevant T-cell subset in RA pathogenesis.

Conserved 33-kb haplotype in the MHC class III region regulates chronic arthritis.

Genome-wide association studies have revealed many genetic loci associated with complex autoimmune diseases. In rheumatoid arthritis (RA), the MHC gene HLA-DRB1 is the strongest candidate predicting disease development. It has been suggested that other immune-regulating genes in the MHC contribute to the disease risk, but this contribution has been difficult to show because of the strong linkage disequilibrium within the MHC. We isolated genomic regions in the form of congenic fragments in rats to test whether there are additional susceptibility loci in the MHC. By both congenic mapping in inbred strains and SNP typing in wild rats, we identified a conserved, 33-kb large haplotype Ltab-Ncr3 in the MHC-III region, which regulates the onset, severity, and chronicity of arthritis. The Ltab-Ncr3 haplotype consists of five polymorphic immunoregulatory genes: Lta (lymphotoxin-α), Tnf, Ltb (lymphotoxin-ß), Lst1 (leukocyte-specific transcript 1), and Ncr3 (natural cytotoxicity-triggering receptor 3). Significant correlation in the expression of the Ltab-Ncr3 genes suggests that interaction of these genes may be important in keeping these genes clustered together as a conserved haplotype. We studied the arthritis association and the spliceo-transcriptome of four different Ltab-Ncr3 haplotypes and showed that higher Ltb and Ncr3 expression, lower Lst1 expression, and the expression of a shorter splice variant of Lst1 correlate with reduced arthritis severity in rats. Interestingly, patients with mild RA also showed higher NCR3 expression and lower LST1 expression than patients with severe RA. These data demonstrate the importance of a conserved haplotype in the regulation of complex diseases such as arthritis.

T cells are influenced by a long non-coding RNA in the autoimmune associated PTPN2 locus.

Non-coding SNPs in the protein tyrosine phosphatase non-receptor type 2 (PTPN2) locus have been linked with several autoimmune diseases, including rheumatoid arthritis, type I diabetes, and inflammatory bowel disease. However, the functional consequences of these SNPs are poorly characterized. Herein, we show in blood cells that SNPs in the PTPN2 locus are highly correlated with DNA methylation levels at four CpG sites downstream of PTPN2 and expression levels of the long non-coding RNA (lncRNA) LINC01882 downstream of these CpG sites. We observed that LINC01882 is mainly expressed in T cells and that anti-CD3/CD28 activated naïve CD4+ T cells downregulate the expression of LINC01882. RNA sequencing analysis of LINC01882 knockdown in Jurkat T cells, using a combination of antisense oligonucleotides and RNA interference, revealed the upregulation of the transcription factor ZEB1 and kinase MAP2K4, both involved in IL-2 regulation. Overall, our data suggests the involvement of LINC01882 in T cell activation and hints towards an auxiliary role of these non-coding SNPs in autoimmunity associated with the PTPN2 locus.

Molecular maps of synovial cells in inflammatory arthritis using an optimized synovial tissue dissociation protocol.

In this study, we optimized the dissociation of synovial tissue biopsies for single-cell omics studies and created a single-cell atlas of human synovium in inflammatory arthritis. The optimized protocol allowed consistent isolation of highly viable cells from tiny fresh synovial biopsies, minimizing the synovial biopsy drop-out rate. The synovium scRNA-seq atlas contained over 100,000 unsorted synovial cells from 25 synovial tissues affected by inflammatory arthritis, including 16 structural, 11 lymphoid, and 15 myeloid cell clusters. This synovial cell map expanded the diversity of synovial cell types/states, detected synovial neutrophils, and broadened synovial endothelial cell classification. We revealed tissue-resident macrophage subsets with proposed matrix-sensing (FOLR2+COLEC12high) and iron-recycling (LYVE1+SLC40A1+) activities and identified fibroblast subsets with proposed functions in cartilage breakdown (SOD2highSAA1+SAA2+SDC4+) and extracellular matrix remodeling (SERPINE1+COL5A3+LOXL2+). Our study offers an efficient synovium dissociation method and a reference scRNA-seq resource, that advances the current understanding of synovial cell heterogeneity in inflammatory arthritis.

Bromodomain Protein Inhibitors Reorganize the Chromatin of Synovial Fibroblasts.

Bromodomain- and extra-terminal domain (BET) proteins are epigenetic reader proteins that regulate transcription of their target genes by binding to acetylated histone side chains. Small molecule inhibitors, such as I-BET151, have anti-inflammatory properties in fibroblast-like synoviocytes (FLS) and in animal models of arthritis. Here, we investigated whether BET inhibition can also affect the levels of histone modifications, a novel mechanism underlying BET protein inhibition. On the one hand, FLSs were treated with I-BET151 (1 µM) for 24 h in absence and presence of TNF. On the other hand, FLSs were washed with PBS after 48 h of I-BET151 treatment, and the effects were measured 5 days after I-BET151 treatment or after an additional 24 h stimulation with TNF (5 d + 24 h). Mass spectrometry analysis indicated that I-BET151 induced profound changes in histone modifications, with a global reduction in acetylation on different histone side chains 5 days after treatment. We confirmed changes on acetylated histone side chains in independent samples by Western blotting. I-BET151 treatment reduced mean TNF-induced levels of total acetylated histone 3 (acH3), H3K18ac, and H3K27ac. In line with these changes, the TNF-induced expression of BET protein target genes was suppressed 5 d after I-BET151 treatment. Our data indicate that BET inhibitors not only prevent the reading of acetylated histones but directly influence overall chromatin organization, in particular after stimulation with TNF.

Reset of Inflammatory Priming of Joint Tissue and Reduction of the Severity of Arthritis Flares by Bromodomain Inhibition.

OBJECTIVE: We have recently shown that priming of synovial fibroblasts (SFs) drives arthritis flares. Pathogenic priming of SFs is essentially mediated by epigenetic reprogramming. Bromodomain and extraterminal motif (BET) proteins translate epigenetic changes into transcription. Here, we used a BET inhibitor (I-BET151) to target inflammatory tissue priming and to reduce flare severity in a murine experimental arthritis model. METHODS: BALB/c mice were treated by intraperitoneal injection or by local injection in the paw with I-BET151, which blocks the interaction of BET proteins with acetylated histones. We assessed the effects of I-BET151 on acute arthritis and/or inflammatory tissue priming in a model of repeated injections of monosodium urate crystals or zymosan into the mouse paw. I-BET151 was given before arthritis induction, at peak inflammation, or after healing of the first arthritis bout. We performed transcriptomic (RNA-Seq), epigenomic (ATAC-Seq), and functional (invasion, cytokine production, migration, senescence, metabolic flux) analyses of murine and human SFs treated with I-BET151 in vitro or in vivo. RESULTS: Systemic I-BET151 administration did not affect acute inflammation but abolished inflammatory tissue priming and diminished flare severity in both preventive and therapeutic treatment settings. I-BET151 was also effective when applied locally in the joint. BET inhibition also inhibited osteoclast differentiation, while macrophage activation in the joint was not affected. Flare reduction after BET inhibition was mediated, at least in part, by rolling back the primed transcriptional, metabolic, and pathogenic phenotype of SFs. CONCLUSION: Inflammatory tissue priming is dependent on transcriptional regulation by BET proteins, making them promising therapeutic targets for prevention of arthritis flares in previously affected joints.

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.

Five commercially-available antibodies react differentially with allelic forms of human HLA-DR beta chain.

Allelic variants of HLA-DRB1 have been associated with a variety of autoimmune and infectious diseases. Although the precise molecular mechanisms by which HLA-DRB1 alleles predispose to a particular disease are currently unclear, it has been shown that mRNA expression levels of HLA-DRB1 are dependent on the different alleles. We aimed to measure HLA-DR beta chain levels in peripheral blood mononuclear cells of individuals carrying HLA-DRB1*03:01/*04:01 and HLA-DRB1*03:01/*15:01 alleles by western blotting, using five commercially-available HLA-DRB antibodies. We observed highly heterogeneous binding of the tested antibodies to the different allelic forms of the HLA-DR beta chain. Overall, we show that current immunological research that employs available antibodies to detect HLA-DR beta chains is biased towards detection of specific variants of the protein; this may cause significant discrepancy in quantification of protein expression in a heterogeneous human population.

An Optimized Tissue Dissociation Protocol for Single-Cell RNA Sequencing Analysis of Fresh and Cultured Human Skin Biopsies.

We present an optimized dissociation protocol for preparing high-quality skin cell suspensions for in-depth single-cell RNA-sequencing (scRNA-seq) analysis of fresh and cultured human skin. Our protocol enabled the isolation of a consistently high number of highly viable skin cells from small freshly dissociated punch skin biopsies, which we use for scRNA-seq studies. We recapitulated not only the main cell populations of existing single-cell skin atlases, but also identified rare cell populations, such as mast cells. Furthermore, we effectively isolated highly viable single cells from ex vivo cultured skin biopsy fragments and generated a global single-cell map of the explanted human skin. The quality metrics of the generated scRNA-seq datasets were comparable between freshly dissociated and cultured skin. Overall, by enabling efficient cell isolation and comprehensive cell mapping, our skin dissociation-scRNA-seq workflow can greatly facilitate scRNA-seq discoveries across diverse human skin pathologies and ex vivo skin explant experimentations.

Haplotype-Specific Expression Analysis of MHC Class II Genes in Healthy Individuals and Rheumatoid Arthritis Patients.

HLA-DRB1 alleles have been associated with several autoimmune diseases. For anti-citrullinated protein antibody positive rheumatoid arthritis (RA), HLA-DRB1 shared epitope (SE) alleles are the major genetic risk factors. In order to study the genetic regulation of major histocompatibility complex (MHC) Class II gene expression in immune cells, we investigated transcriptomic profiles of a variety of immune cells from healthy individuals carrying different HLA-DRB1 alleles. Sequencing libraries from peripheral blood mononuclear cells, CD4+ T cells, CD8+ T cells, and CD14+ monocytes of 32 genetically pre-selected healthy female individuals were generated, sequenced and reads were aligned to the standard reference. For the MHC region, reads were mapped to available MHC reference haplotypes and AltHapAlignR was used to estimate gene expression. Using this method, HLA-DRB and HLA-DQ were found to be differentially expressed in different immune cells of healthy individuals as well as in whole blood samples of RA patients carrying HLA-DRB1 SE-positive versus SE-negative alleles. In contrast, no genes outside the MHC region were differentially expressed between individuals carrying HLA-DRB1 SE-positive and SE-negative alleles, thus HLA-DRB1 SE alleles have a strong cis effect on gene expression. Altogether, our findings suggest that immune effects associated with different allelic forms of HLA-DR and HLA-DQ may be associated not only with differences in the structure of these proteins, but also with differences in their expression levels.

Functional genomics atlas of synovial fibroblasts defining rheumatoid arthritis heritability.

BACKGROUND: Genome-wide association studies have reported more than 100 risk loci for rheumatoid arthritis (RA). These loci are shown to be enriched in immune cell-specific enhancers, but the analysis so far has excluded stromal cells, such as synovial fibroblasts (FLS), despite their crucial involvement in the pathogenesis of RA. Here we integrate DNA architecture, 3D chromatin interactions, DNA accessibility, and gene expression in FLS, B cells, and T cells with genetic fine mapping of RA loci. RESULTS: We identify putative causal variants, enhancers, genes, and cell types for 30-60% of RA loci and demonstrate that FLS account for up to 24% of RA heritability. TNF stimulation of FLS alters the organization of topologically associating domains, chromatin state, and the expression of putative causal genes such as TNFAIP3 and IFNAR1. Several putative causal genes constitute RA-relevant functional networks in FLS with roles in cellular proliferation and activation. Finally, we demonstrate that risk variants can have joint-specific effects on target gene expression in RA FLS, which may contribute to the development of the characteristic pattern of joint involvement in RA. CONCLUSION: Overall, our research provides the first direct evidence for a causal role of FLS in the genetic susceptibility for RA accounting for up to a quarter of RA heritability.

T-cell transcriptomics from peripheral blood highlights differences between polymyositis and dermatomyositis patients.

BACKGROUND: Polymyositis (PM) and dermatomyositis (DM) are two distinct subgroups of idiopathic inflammatory myopathies, a chronic inflammatory disorder clinically characterized by muscle weakness and inflammatory cell infiltrates in muscle tissue. In PM, a major component of inflammatory cell infiltrates is CD8+ T cells, whereas in DM, CD4+ T cells, plasmacytoid dendritic cells, and B cells predominate. In this study, with the aim to differentiate involvement of CD4+ and CD8+ T-cell subpopulations in myositis subgroups, we investigated transcriptomic profiles of T cells from peripheral blood of patients with myositis. METHODS: Total RNA was extracted from CD4+ T cells (PM = 8 and DM = 7) and CD8+ T cells (PM = 4 and DM = 5) that were isolated from peripheral blood mononuclear cells via positive selection using microbeads. Sequencing libraries were generated using the Illumina TruSeq Stranded Total RNA Kit and sequenced on an Illumina HiSeq 2500 platform, yielding about 50 million paired-end reads per sample. Differential gene expression analyses were conducted using DESeq2. RESULTS: In CD4+ T cells, only two genes, ANKRD55 and S100B, were expressed significantly higher in patients with PM than in patients with DM (false discovery rate [FDR] < 0.05, model adjusted for age, sex, HLA-DRB1*03 status, and RNA integrity number [RIN]). On the contrary, in CD8+ T cells, 176 genes were differentially expressed in patients with PM compared with patients with DM. Of these, 44 genes were expressed significantly higher in CD8+ T cells from patients with PM, and 132 genes were expressed significantly higher in CD8+ T cells from patients with DM (FDR < 0.05, model adjusted for age, sex, and RIN). Gene Ontology analysis showed that genes differentially expressed in CD8+ T cells are involved in lymphocyte migration and regulation of T-cell differentiation. CONCLUSIONS: Our data strongly suggest that CD8+ T cells represent a major divergence between PM and DM patients compared with CD4+ T cells. These alterations in the gene expression in T cells from PM and DM patients might advocate for distinct immune mechanisms in these subphenotypes of myositis.

An antisense RNA capable of modulating the expression of the tumor suppressor microRNA-34a.

The microRNA-34a is a well-studied tumor suppressor microRNA (miRNA) and a direct downstream target of TP53 with roles in several pathways associated with oncogenesis, such as proliferation, cellular growth, and differentiation. Due to its broad tumor suppressive activity, it is not surprising that miR34a expression is altered in a wide variety of solid tumors and hematological malignancies. However, the mechanisms by which miR34a is regulated in these cancers is largely unknown. In this study, we find that a long noncoding RNA transcribed antisense to the miR34a host gene, is critical for miR34a expression and mediation of its cellular functions in multiple types of human cancer. We name this long noncoding RNA lncTAM34a, and characterize its ability to facilitate miR34a expression under different types of cellular stress in both TP53-deficient and wild-type settings.