Single-cell transcriptomic analyses reveal distinct immune cell contributions to epithelial barrier dysfunction in checkpoint inhibitor colitis.

Immune checkpoint inhibitor (ICI) therapy has revolutionized oncology, but treatments are limited by immune-related adverse events, including checkpoint inhibitor colitis (irColitis). Little is understood about the pathogenic mechanisms driving irColitis, which does not readily occur in model organisms, such as mice. To define molecular drivers of irColitis, we used single-cell multi-omics to profile approximately 300,000 cells from the colon mucosa and blood of 13 patients with cancer who developed irColitis (nine on anti-PD-1 or anti-CTLA-4 monotherapy and four on dual ICI therapy; most patients had skin or lung cancer), eight controls on ICI therapy and eight healthy controls. Patients with irColitis showed expanded mucosal Tregs, ITGAEHi CD8 tissue-resident memory T cells expressing CXCL13 and Th17 gene programs and recirculating ITGB2Hi CD8 T cells. Cytotoxic GNLYHi CD4 T cells, recirculating ITGB2Hi CD8 T cells and endothelial cells expressing hypoxia gene programs were further expanded in colitis associated with anti-PD-1/CTLA-4 therapy compared to anti-PD-1 therapy. Luminal epithelial cells in patients with irColitis expressed PCSK9, PD-L1 and interferon-induced signatures associated with apoptosis, increased cell turnover and malabsorption. Together, these data suggest roles for circulating T cells and epithelial-immune crosstalk critical to PD-1/CTLA-4-dependent tolerance and barrier function and identify potential therapeutic targets for irColitis.

Blood immunophenotyping identifies distinct kidney histopathology and outcomes in patients with lupus nephritis.

Lupus nephritis (LN) is a frequent manifestation of systemic lupus erythematosus, and fewer than half of patients achieve complete renal response with standard immunosuppressants. Identifying non-invasive, blood-based pathologic immune alterations associated with renal injury could aid therapeutic decisions. Here, we used mass cytometry immunophenotyping of peripheral blood mononuclear cells in 145 patients with biopsy-proven LN and 40 healthy controls to evaluate the heterogeneity of immune activation in patients with LN and to identify correlates of renal parameters and treatment response. Unbiased analysis identified 3 immunologically distinct groups of patients with LN that were associated with different patterns of histopathology, renal cell infiltrates, urine proteomic profiles, and treatment response at one year. Patients with enriched circulating granzyme B+ T cells at baseline showed more severe disease and increased numbers of activated CD8 T cells in the kidney, yet they had the highest likelihood of treatment response. A second group characterized primarily by a high type I interferon signature had a lower likelihood of response to therapy, while a third group appeared immunologically inactive by immunophenotyping at enrollment but with chronic renal injuries. Main immune profiles could be distilled down to 5 simple cytometric parameters that recapitulate several of the associations, highlighting the potential for blood immune profiling to translate to clinically useful non-invasive metrics to assess immune-mediated disease in LN.

Cytotoxic CD4+ T cells eliminate senescent cells by targeting cytomegalovirus antigen.

Senescent cell accumulation has been implicated in the pathogenesis of aging-associated diseases, including cancer. The mechanism that prevents the accumulation of senescent cells in aging human organs is unclear. Here, we demonstrate that a virus-immune axis controls the senescent fibroblast accumulation in the human skin. Senescent fibroblasts increased in old skin compared with young skin. However, they did not increase with advancing age in the elderly. Increased CXCL9 and cytotoxic CD4+ T cells (CD4 CTLs) recruitment were significantly associated with reduced senescent fibroblasts in the old skin. Senescent fibroblasts expressed human leukocyte antigen class II (HLA-II) and human cytomegalovirus glycoprotein B (HCMV-gB), becoming direct CD4 CTL targets. Skin-resident CD4 CTLs eliminated HCMV-gB+ senescent fibroblasts in an HLA-II-dependent manner, and HCMV-gB activated CD4 CTLs from the human skin. Collectively, our findings demonstrate HCMV reactivation in senescent cells, which CD4 CTLs can directly eliminate through the recognition of the HCMV-gB antigen.

Systematic identification of genomic elements that regulate FCGR2A expression and harbor variants linked with autoimmune disease.

BACKGROUND: FCGR2A binds antibody-antigen complexes to regulate the abundance of circulating and deposited complexes along with downstream immune and autoimmune responses. Although the abundance of FCRG2A may be critical in immune-mediated diseases, little is known about whether its surface expression is regulated through cis genomic elements and non-coding variants. In the current study, we aimed to characterize the regulation of FCGR2A expression, the impact of genetic variation and its association with autoimmune disease. METHODS: We applied CRISPR-based interference and editing to scrutinize 1.7 Mb of open chromatin surrounding the FCGR2A gene to identify regulatory elements. Relevant transcription factors (TFs) binding to these regions were defined through public databases. Genetic variants affecting regulation were identified using luciferase reporter assays and were verified in a cohort of 1996 genotyped healthy individuals using flow cytometry. RESULTS: We identified a complex proximal region and five distal enhancers regulating FCGR2A. The proximal region split into subregions upstream and downstream of the transcription start site, was enriched in binding of inflammation-regulated TFs, and harbored a variant associated with FCGR2A expression in primary myeloid cells. One distal enhancer region was occupied by CCCTC-binding factor (CTCF) whose binding site was disrupted by a rare genetic variant, altering gene expression. CONCLUSIONS: The FCGR2A gene is regulated by multiple proximal and distal genomic regions, with links to autoimmune disease. These findings may open up novel therapeutic avenues where fine-tuning of FCGR2A levels may constitute a part of treatment strategies for immune-mediated diseases.

Differentiation of exhausted CD8+ T cells after termination of chronic antigen stimulation stops short of achieving functional T cell memory.

T cell exhaustion is associated with failure to clear chronic infections and malignant cells. Defining the molecular mechanisms of T cell exhaustion and reinvigoration is essential to improving immunotherapeutic modalities. Here we confirmed pervasive phenotypic, functional and transcriptional differences between memory and exhausted antigen-specific CD8+ T cells in human hepatitis C virus (HCV) infection before and after treatment. After viral cure, phenotypic changes in clonally stable exhausted T cell populations suggested differentiation toward a memory-like profile. However, functionally, the cells showed little improvement, and critical transcriptional regulators remained in the exhaustion state. Notably, T cells from chronic HCV infection that were exposed to antigen for less time because of viral escape mutations were functionally and transcriptionally more similar to memory T cells from spontaneously resolved HCV infection. Thus, the duration of T cell stimulation impacts exhaustion recovery, with antigen removal after long-term exhaustion being insufficient for the development of functional T cell memory.

Genome-wide enhancer maps link risk variants to disease genes.

Genome-wide association studies (GWAS) have identified thousands of noncoding loci that are associated with human diseases and complex traits, each of which could reveal insights into the mechanisms of disease1. Many of the underlying causal variants may affect enhancers2,3, but we lack accurate maps of enhancers and their target genes to interpret such variants. We recently developed the activity-by-contact (ABC) model to predict which enhancers regulate which genes and validated the model using CRISPR perturbations in several cell types4. Here we apply this ABC model to create enhancer-gene maps in 131 human cell types and tissues, and use these maps to interpret the functions of GWAS variants. Across 72 diseases and complex traits, ABC links 5,036 GWAS signals to 2,249 unique genes, including a class of 577 genes that appear to influence multiple phenotypes through variants in enhancers that act in different cell types. In inflammatory bowel disease (IBD), causal variants are enriched in predicted enhancers by more than 20-fold in particular cell types such as dendritic cells, and ABC achieves higher precision than other regulatory methods at connecting noncoding variants to target genes. These variant-to-function maps reveal an enhancer that contains an IBD risk variant and that regulates the expression of PPIF to alter the membrane potential of mitochondria in macrophages. Our study reveals principles of genome regulation, identifies genes that affect IBD and provides a resource and generalizable strategy to connect risk variants of common diseases to their molecular and cellular functions.

Differentiation of exhausted CD8 T cells after termination of chronic antigen stimulation stops short of achieving functional T cell memory

T cell exhaustion is associated with failure to clear chronic infections and malignant cells. Defining the molecular mechanisms of T cell exhaustion and reinvigoration is essential to improving immunotherapeutic modalities. Analysis of antigen-specific CD8+ T cells before and after antigen removal in human hepatitis C virus (HCV) infection confirmed pervasive phenotypic, functional, and transcriptional differences between exhausted and memory CD8+ T cells. After viral cure, we observed broad phenotypic and transcriptional changes in clonally stable exhausted T-cell populations suggesting differentiation towards a memory-like profile. However, functionally, the cells showed little improvement and critical transcriptional regulators remained in the exhaustion state. Notably, T cells from chronic HCV infection that were exposed to antigen for shorter periods of time because of viral escape mutations were functionally and transcriptionally more similar to memory T cells from spontaneously resolved acute HCV infection. Thus, duration of T cell stimulation impacts the ability to recover from exhaustion, as antigen removal after long-term T cell exhaustion is insufficient for the development of key T cell memory characteristics.

Loss of the Nuclear Protein RTF2 Enhances Influenza Virus Replication.

While hundreds of genes are induced by type I interferons, their roles in restricting the influenza virus life cycle remain mostly unknown. Using a loss-of-function CRISPR screen in cells prestimulated with interferon beta (IFN-ß), we identified a small number of factors required for restricting influenza A virus replication. In addition to known components of the interferon signaling pathway, we found that replication termination factor 2 (RTF2) restricts influenza virus at the nuclear stage (and perhaps other stages) of the viral life cycle, based on several lines of evidence. First, a deficiency in RTF2 leads to higher levels of viral primary transcription, even in the presence of cycloheximide to block genome replication and secondary transcription. Second, cells that lack RTF2 have enhanced activity of a viral reporter that depends solely on four viral proteins that carry out replication and transcription in the nucleus. Third, when the RTF2 protein is mislocalized outside the nucleus, it is not able to restrict replication. Finally, the absence of RTF2 leads not only to enhanced viral transcription but also to reduced expression of antiviral factors in response to interferon. RTF2 thus inhibits primary influenza virus transcription, likely acts in the nucleus, and contributes to the upregulation of antiviral effectors in response to type I interferons.IMPORTANCE Viral infection triggers the secretion of type I interferons, which in turn induce the expression of hundreds of antiviral genes. However, the roles of these induced genes in controlling viral infections remain largely unknown, limiting our ability to develop host-based antiviral therapeutics against pathogenic viruses, such as influenza virus. Here, we performed a loss-of-function genetic CRISPR screen in cells prestimulated with type I interferon to identify antiviral genes that restrict influenza A virus replication. Besides finding key components of the interferon signaling pathway, we discovered a new restriction factor, RTF2, which acts in the nucleus, restricts influenza virus transcription, and contributes to the interferon-induced upregulation of known restriction factors. Our work contributes to the field of antiviral immunology by discovering and characterizing a novel restriction factor of influenza virus and may ultimately be useful for understanding how to control a virus that causes significant morbidity and mortality worldwide.

Disabled Homolog 2 Controls Prometastatic Activity of Tumor-Associated Macrophages.

Tumor-associated macrophages (TAM) are regulators of extracellular matrix (ECM) remodeling and metastatic progression, the main cause of cancer-associated death. We found that disabled homolog 2 mitogen-responsive phosphoprotein (DAB2) is highly expressed in tumor-infiltrating TAMs and that its genetic ablation significantly impairs lung metastasis formation. DAB2-expressing TAMs, mainly localized along the tumor-invasive front, participate in integrin recycling, ECM remodeling, and directional migration in a tridimensional matrix. DAB2+ macrophages escort the invasive dissemination of cancer cells by a mechanosensing pathway requiring the transcription factor YAP. In human lobular breast and gastric carcinomas, DAB2+ TAMs correlated with a poor clinical outcome, identifying DAB2 as potential prognostic biomarker for stratification of patients with cancer. DAB2 is therefore central for the prometastatic activity of TAMs. SIGNIFICANCE: DAB2 expression in macrophages is essential for metastasis formation but not primary tumor growth. Mechanosensing cues, activating the complex YAP-TAZ, regulate DAB2 in macrophages, which in turn controls integrin recycling and ECM remodeling in 3-D tissue matrix. The presence of DAB2+ TAMs in patients with cancer correlates with worse prognosis.This article is highlighted in the In This Issue feature, p. 1611.

SLAMF6​ deficiency augments tumor killing and skews toward an effector phenotype revealing it as a novel T cell checkpoint.

SLAMF6 is a homotypic receptor of the Ig-superfamily whose exact role in immune modulation has remained elusive. Its constitutive expression on resting and activated T cells precludes it from being a bona fide exhaustion marker. By breeding Pmel-1 mice with SLAMF6 -/- mice, we generated donors for T cells lacking SLAMF6 and expressing a transgenic TCR for gp100-melanoma antigen. Activated Pmel-1xSLAMF6 -/- CD8+ T cells displayed improved polyfunctionality and strong tumor cytolysis. T-bet was the dominant transcription factor in Pmel-1 x SLAMF6 -/- cells, and upon activation, they acquired an effector-memory phenotype. Adoptive transfer of Pmel-1 x SLAMF6 -/- T cells to melanoma-bearing mice resulted in lasting tumor regression in contrast to temporary responses achieved with Pmel-1 T cells. LAG-3 expression was elevated in the SLAMF6 -/- cells, and the addition of the LAG-3-blocking antibody to the adoptive transfer protocol improved the SLAMF6 -/- T cells and expedited the antitumor response even further. The results from this study support the notion that SLAMF6 is an inhibitory immune receptor whose absence enables powerful CD8+ T cells to eradicate tumors.

The immune system helps to protect our bodies from illnesses and infections. Immunotherapies are medicines designed to treat diseases, such as cancer, by boosting the immune system against the condition. This is a powerful approach but so far immunotherapies have only had partial success and there is a need for further improvements. One protein called SLAMF6 is found on cells from the immune system that attack and kill cancer cells. Immunotherapies that suppress SLAMF6 on immune cells called killer T cells could increase immune system activity helping to treat cancers, particularly melanoma skin cancers. So far the potential for SLAMF6 as a target for immunotherapy has not been fully explored. Hajaj et al. created mice with killer T cells that recognized skin cancer cells and lacked SLAMF6. These modified cells were better at fighting cancer, producing more anti-cancer chemicals called cytokines and killing more cancer cells. The modified cells had a lasting effect on tumors and helped the mice to live longer. The effects could be further boosted by treating the mice in combination with other immunotherapies. SLAMF6 is a possible new target for skin cancer immunotherapy that could help more people to live longer following cancer diagnosis. The next step is to create a drug to target SLAMF6 in humans and to test it in clinical trials.

An immune-cell signature of bacterial sepsis.

Dysregulation of the immune response to bacterial infection can lead to sepsis, a condition with high mortality. Multiple whole-blood gene-expression studies have defined sepsis-associated molecular signatures, but have not resolved changes in transcriptional states of specific cell types. Here, we used single-cell RNA-sequencing to profile the blood of people with sepsis (n = 29) across three clinical cohorts with corresponding controls (n = 36). We profiled total peripheral blood mononuclear cells (PBMCs, 106,545 cells) and dendritic cells (19,806 cells) across all subjects and, on the basis of clustering of their gene-expression profiles, defined 16 immune-cell states. We identified a unique CD14+ monocyte state that is expanded in people with sepsis and validated its power in distinguishing these individuals from controls using public transcriptomic data from subjects with different disease etiologies and from multiple geographic locations (18 cohorts, n = 1,467 subjects). We identified a panel of surface markers for isolation and quantification of the monocyte state and characterized its epigenomic and functional phenotypes, and propose a model for its induction from human bone marrow. This study demonstrates the utility of single-cell genomics in discovering disease-associated cytologic signatures and provides insight into the cellular basis of immune dysregulation in bacterial sepsis.

Genome-wide CRISPR screen identifies host dependency factors for influenza A virus infection.

Host dependency factors that are required for influenza A virus infection may serve as therapeutic targets as the virus is less likely to bypass them under drug-mediated selection pressure. Previous attempts to identify host factors have produced largely divergent results, with few overlapping hits across different studies. Here, we perform a genome-wide CRISPR/Cas9 screen and devise a new approach, meta-analysis by information content (MAIC) to systematically combine our results with prior evidence for influenza host factors. MAIC out-performs other meta-analysis methods when using our CRISPR screen as validation data. We validate the host factors, WDR7, CCDC115 and TMEM199, demonstrating that these genes are essential for viral entry and regulation of V-type ATPase assembly. We also find that CMTR1, a human mRNA cap methyltransferase, is required for efficient viral cap snatching and regulation of a cell autonomous immune response, and provides synergistic protection with the influenza endonuclease inhibitor Xofluza.

A large peptidome dataset improves HLA class I epitope prediction across most of the human population.

Prediction of HLA epitopes is important for the development of cancer immunotherapies and vaccines. However, current prediction algorithms have limited predictive power, in part because they were not trained on high-quality epitope datasets covering a broad range of HLA alleles. To enable prediction of endogenous HLA class I-associated peptides across a large fraction of the human population, we used mass spectrometry to profile >185,000 peptides eluted from 95 HLA-A, -B, -C and -G mono-allelic cell lines. We identified canonical peptide motifs per HLA allele, unique and shared binding submotifs across alleles and distinct motifs associated with different peptide lengths. By integrating these data with transcript abundance and peptide processing, we developed HLAthena, providing allele-and-length-specific and pan-allele-pan-length prediction models for endogenous peptide presentation. These models predicted endogenous HLA class I-associated ligands with 1.5-fold improvement in positive predictive value compared with existing tools and correctly identified >75% of HLA-bound peptides that were observed experimentally in 11 patient-derived tumor cell lines.

Publisher Correction: The immune cell landscape in kidneys of patients with lupus nephritis.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The immune cell landscape in kidneys of patients with lupus nephritis.

Lupus nephritis is a potentially fatal autoimmune disease for which the current treatment is ineffective and often toxic. To develop mechanistic hypotheses of disease, we analyzed kidney samples from patients with lupus nephritis and from healthy control subjects using single-cell RNA sequencing. Our analysis revealed 21 subsets of leukocytes active in disease, including multiple populations of myeloid cells, T cells, natural killer cells and B cells that demonstrated both pro-inflammatory responses and inflammation-resolving responses. We found evidence of local activation of B cells correlated with an age-associated B-cell signature and evidence of progressive stages of monocyte differentiation within the kidney. A clear interferon response was observed in most cells. Two chemokine receptors, CXCR4 and CX3CR1, were broadly expressed, implying a potentially central role in cell trafficking. Gene expression of immune cells in urine and kidney was highly correlated, which would suggest that urine might serve as a surrogate for kidney biopsies.

Defining inflammatory cell states in rheumatoid arthritis joint synovial tissues by integrating single-cell transcriptomics and mass cytometry.

To define the cell populations that drive joint inflammation in rheumatoid arthritis (RA), we applied single-cell RNA sequencing (scRNA-seq), mass cytometry, bulk RNA sequencing (RNA-seq) and flow cytometry to T cells, B cells, monocytes, and fibroblasts from 51 samples of synovial tissue from patients with RA or osteoarthritis (OA). Utilizing an integrated strategy based on canonical correlation analysis of 5,265 scRNA-seq profiles, we identified 18 unique cell populations. Combining mass cytometry and transcriptomics revealed cell states expanded in RA synovia: THY1(CD90)+HLA-DRAhi sublining fibroblasts, IL1B+ pro-inflammatory monocytes, ITGAX+TBX21+ autoimmune-associated B cells and PDCD1+ peripheral helper T (TPH) cells and follicular helper T (TFH) cells. We defined distinct subsets of CD8+ T cells characterized by GZMK+, GZMB+, and GNLY+ phenotypes. We mapped inflammatory mediators to their source cell populations; for example, we attributed IL6 expression to THY1+HLA-DRAhi fibroblasts and IL1B production to pro-inflammatory monocytes. These populations are potentially key mediators of RA pathogenesis.

Multiplexed enrichment and genomic profiling of peripheral blood cells reveal subset-specific immune signatures.

Specialized immune cell subsets are involved in autoimmune disease, cancer immunity, and infectious disease through a diverse range of functions mediated by overlapping pathways and signals. However, subset-specific responses may not be detectable in analyses of whole blood samples, and no efficient approach for profiling cell subsets at high throughput from small samples is available. We present a low-input microfluidic system for sorting immune cells into subsets and profiling their gene expression. We validate the system's technical performance against standard subset isolation and library construction protocols and demonstrate the importance of subset-specific profiling through in vitro stimulation experiments. We show the ability of this integrated platform to identify subset-specific disease signatures by profiling four immune cell subsets in blood from patients with systemic lupus erythematosus (SLE) and matched control subjects. The platform has the potential to make multiplexed subset-specific analysis routine in many research laboratories and clinical settings.

Extranuclear DNA accumulates in aged cells and contributes to senescence and inflammation.

Systemic inflammation is central to aging-related conditions. However, the intrinsic factors that induce inflammation are not well understood. We previously identified a cell-autonomous pathway through which damaged nuclear DNA is trafficked to the cytosol where it activates innate cytosolic DNA sensors that trigger inflammation. These results led us to hypothesize that DNA released after cumulative damage contributes to persistent inflammation in aging cells through a similar mechanism. Consistent with this notion, we found that older cells harbored higher levels of extranuclear DNA compared to younger cells. Extranuclear DNA was exported by a leptomycin B-sensitive process, degraded through the autophagosome-lysosomal pathway and triggered innate immune responses through the DNA-sensing cGAS-STING pathway. Patient cells from the aging diseases ataxia and progeria also displayed extranuclear DNA accumulation, increased pIRF3 and pTBK1, and STING-dependent p16 expression. Removing extranuclear DNA in old cells using DNASE2A reduced innate immune responses and senescence-associated (SA) ß-gal enzyme activity. Cells and tissues of Dnase2a-/- mice with defective DNA degradation exhibited slower growth, higher activity of ß-gal, or increased expression of HP-1ß and p16 proteins, while Dnase2a-/- ;Sting-/- cells and tissues were rescued from these phenotypes, supporting a role for extranuclear DNA in senescence. We hypothesize a direct role for excess DNA in aging-related inflammation and in replicative senescence, and propose DNA degradation as a therapeutic approach to remove intrinsic DNA and revert inflammation associated with aging.

A secreted PD-L1 splice variant that covalently dimerizes and mediates immunosuppression.

Targeting immune checkpoint pathways, such as programmed death ligand-1 (PD-L1, also known as CD274 or B7-H1) or its receptor programmed cell death-1 (PD-1) has shown improved survival for patients with numerous types of cancers, not limited to lung cancer, melanoma, renal cell carcinoma, and Hodgkin lymphoma. PD-L1 is a co-inhibitory molecule whose expression on the surface of tumor cells is associated with worse prognosis in many tumors. Here we describe a splice variant (secPD-L1) that does not splice into the transmembrane domain, but instead produces a secreted form of PD-L1 that has a unique 18 amino acid tail containing a cysteine that allows it to homodimerize and more effectively inhibit lymphocyte function than monomeric soluble PD-L1. We show that recombinant secPD-L1 can dimerize and inhibit T-cell proliferation and IFN-gamma production in vitro. The secPD-L1 variant is expressed by malignant cells in vitro that also express high levels of full-length PD-L1. Transcriptomic analysis of gene expression across The Cancer Genome Atlas found the strongest association of secPD-L1 with full-length PD-L1, but also with subsets of immunologic genes, such as in myeloid-derived suppressor cells. Moreover, the splice variant is also expressed in normal tissues and within normal peripheral blood cells it is preferentially expressed in activated myeloid cells. This is the first report of a form of secreted PD-L1 that homodimerizes and is functionally active. SecPD-L1 may function as a paracrine negative immune regulator within the tumor, since secPD-L1 does not require a cell-to-cell interaction to mediate its inhibitory effect.

Methods for high-dimensional analysis of cells dissociated from cryopreserved synovial tissue

BACKGROUND: Detailed molecular analyses of cells from rheumatoid arthritis (RA) synovium hold promise in identifying cellular phenotypes that drive tissue pathology and joint damage. The Accelerating Medicines Partnership RA/SLE Network aims to deconstruct autoimmune pathology by examining cells within target tissues through multiple high-dimensional assays. Robust standardized protocols need to be developed before cellular phenotypes at a single cell level can be effectively compared across patient samples. METHODS: Multiple clinical sites collected cryopreserved synovial tissue fragments from arthroplasty and synovial biopsy in a 10% DMSO solution. Mechanical and enzymatic dissociation parameters were optimized for viable cell extraction and surface protein preservation for cell sorting and mass cytometry, as well as for reproducibility in RNA sequencing (RNA-seq). Cryopreserved synovial samples were collectively analyzed at a central processing site by a custom-designed and validated 35-marker mass cytometry panel. In parallel, each sample was flow sorted into fibroblast, T-cell, B-cell, and macrophage suspensions for bulk population RNA-seq and plate-based single-cell CEL-Seq2 RNA-seq. RESULTS: Upon dissociation, cryopreserved synovial tissue fragments yielded a high frequency of viable cells, comparable to samples undergoing immediate processing. Optimization of synovial tissue dissociation across six clinical collection sites with ~ 30 arthroplasty and ~ 20 biopsy samples yielded a consensus digestion protocol using 100 µg/ml of Liberase™ TL enzyme preparation. This protocol yielded immune and stromal cell lineages with preserved surface markers and minimized variability across replicate RNA-seq transcriptomes. Mass cytometry analysis of cells from cryopreserved synovium distinguished diverse fibroblast phenotypes, distinct populations of memory B cells and antibody-secreting cells, and multiple CD4+ and CD8+ T-cell activation states. Bulk RNA-seq of sorted cell populations demonstrated robust separation of synovial lymphocytes, fibroblasts, and macrophages. Single-cell RNA-seq produced transcriptomes of over 1000 genes/cell, including transcripts encoding characteristic lineage markers identified. CONCLUSIONS: We have established a robust protocol to acquire viable cells from cryopreserved synovial tissue with intact transcriptomes and cell surface phenotypes. A centralized pipeline to generate multiple high-dimensional analyses of synovial tissue samples collected across a collaborative network was developed. Integrated analysis of such datasets from large patient cohorts may help define molecular heterogeneity within RA pathology and identify new therapeutic targets and biomarkers.