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.

SARS-CoV-2 infection elucidates unique features of pregnancy-specific immunity.

Pregnancy is a risk factor for increased severity of SARS-CoV-2 and other respiratory infections. The mechanisms underlying this risk have not been well-established, partly due to a limited understanding of how pregnancy shapes immune responses. To gain insight into the role of pregnancy in modulating immune responses at steady state and upon perturbation, we collected peripheral blood mononuclear cells (PBMC), plasma, and stool from 226 women, including 152 pregnant individuals (n = 96 with SARS-CoV-2 infection and n = 56 healthy controls) and 74 non-pregnant women (n = 55 with SARS-CoV-2 and n = 19 healthy controls). We found that SARS-CoV-2 infection was associated with altered T cell responses in pregnant compared to non-pregnant women. Differences included a lower percentage of memory T cells, a distinct clonal expansion of CD4-expressing CD8 + T cells, and the enhanced expression of T cell exhaustion markers, such as programmed cell death-1 (PD-1) and T cell immunoglobulin and mucin domain-3 (Tim-3), in pregnant women. We identified additional evidence of immune dysfunction in severely and critically ill pregnant women, including a lack of expected elevation in regulatory T cell (Treg) levels, diminished interferon responses, and profound suppression of monocyte function. Consistent with earlier data, we found maternal obesity was also associated with altered immune responses to SARS-CoV-2 infection, including enhanced production of inflammatory cytokines by T cells. Certain gut bacterial species were altered in pregnancy and upon SARS-CoV-2 infection in pregnant individuals compared to non-pregnant women. Shifts in cytokine and chemokine levels were also identified in the sera of pregnant individuals, most notably a robust increase of interleukin-27 (IL-27), a cytokine known to drive T cell exhaustion, in the pregnant uninfected control group compared to all non-pregnant groups. IL-27 levels were also significantly higher in uninfected pregnant controls compared to pregnant SARS-CoV-2-infected individuals. Using two different preclinical mouse models of inflammation-induced fetal demise and respiratory influenza viral infection, we found that enhanced IL-27 protects developing fetuses from maternal inflammation but renders adult female mice vulnerable to viral infection. These combined findings from human and murine studies reveal nuanced pregnancy-associated immune responses, suggesting mechanisms underlying the increased susceptibility of pregnant individuals to viral respiratory infections.

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.

Immune Responses in Checkpoint Myocarditis Across Heart, Blood, and Tumor.

Immune checkpoint inhibitors (ICIs) are widely used anti-cancer therapies that can cause morbid and potentially fatal immune-related adverse events (irAEs). ICI-related myocarditis (irMyocarditis) is uncommon but has the highest mortality of any irAE. The pathogenesis of irMyocarditis and its relationship to anti-tumor immunity remain poorly understood. We sought to define immune responses in heart, tumor, and blood during irMyocarditis and identify biomarkers of clinical severity by leveraging single-cell (sc)RNA-seq coupled with T cell receptor (TCR) sequencing, microscopy, and proteomics analysis of 28 irMyocarditis patients and 23 controls. Our analysis of 284,360 cells from heart and blood specimens identified cytotoxic T cells, inflammatory macrophages, conventional dendritic cells (cDCs), and fibroblasts enriched in irMyocarditis heart tissue. Additionally, potentially targetable, pro-inflammatory transcriptional programs were upregulated across multiple cell types. TCR clones enriched in heart and paired tumor tissue were largely non-overlapping, suggesting distinct T cell responses within these tissues. We also identify the presence of cardiac-expanded TCRs in a circulating, cycling CD8 T cell population as a novel peripheral biomarker of fatality. Collectively, these findings highlight critical biology driving irMyocarditis and putative biomarkers for therapeutic intervention.

A human model of asthma exacerbation reveals transcriptional programs and cell circuits specific to allergic asthma.

Asthma is a chronic disease most commonly associated with allergy and type 2 inflammation. However, the mechanisms that link airway inflammation to the structural changes that define asthma are incompletely understood. Using a human model of allergen-induced asthma exacerbation, we compared the lower airway mucosa in allergic asthmatics and allergic non-asthmatic controls using single-cell RNA sequencing. In response to allergen, the asthmatic airway epithelium was highly dynamic and up-regulated genes involved in matrix degradation, mucus metaplasia, and glycolysis while failing to induce injury-repair and antioxidant pathways observed in controls. IL9-expressing pathogenic TH2 cells were specific to asthmatic airways and were only observed after allergen challenge. Additionally, conventional type 2 dendritic cells (DC2 that express CD1C) and CCR2-expressing monocyte-derived cells (MCs) were uniquely enriched in asthmatics after allergen, with up-regulation of genes that sustain type 2 inflammation and promote pathologic airway remodeling. In contrast, allergic controls were enriched for macrophage-like MCs that up-regulated tissue repair programs after allergen challenge, suggesting that these populations may protect against asthmatic airway remodeling. Cellular interaction analyses revealed a TH2-mononuclear phagocyte-basal cell interactome unique to asthmatics. These pathogenic cellular circuits were characterized by type 2 programming of immune and structural cells and additional pathways that may sustain and amplify type 2 signals, including TNF family signaling, altered cellular metabolism, failure to engage antioxidant responses, and loss of growth factor signaling. Our findings therefore suggest that pathogenic effector circuits and the absence of proresolution programs drive structural airway disease in response to type 2 inflammation.

Single-cell transcriptomics in cancer: computational challenges and opportunities.

Intratumor heterogeneity is a common characteristic across diverse cancer types and presents challenges to current standards of treatment. Advancements in high-throughput sequencing and imaging technologies provide opportunities to identify and characterize these aspects of heterogeneity. Notably, transcriptomic profiling at a single-cell resolution enables quantitative measurements of the molecular activity that underlies the phenotypic diversity of cells within a tumor. Such high-dimensional data require computational analysis to extract relevant biological insights about the cell types and states that drive cancer development, pathogenesis, and clinical outcomes. In this review, we highlight emerging themes in the computational analysis of single-cell transcriptomics data and their applications to cancer research. We focus on downstream analytical challenges relevant to cancer research, including how to computationally perform unified analysis across many patients and disease states, distinguish neoplastic from nonneoplastic cells, infer communication with the tumor microenvironment, and delineate tumoral and microenvironmental evolution with trajectory and RNA velocity analysis. We include discussions of challenges and opportunities for future computational methodological advancements necessary to realize the translational potential of single-cell transcriptomic profiling in cancer.

Synoviocyte-targeted therapy synergizes with TNF inhibition in arthritis reversal.

Fibroblast-like synoviocytes (FLS) are joint-lining cells that promote rheumatoid arthritis (RA) pathology. Current disease-modifying antirheumatic agents (DMARDs) operate through systemic immunosuppression. FLS-targeted approaches could potentially be combined with DMARDs to improve control of RA without increasing immunosuppression. Here, we assessed the potential of immunoglobulin-like domains 1 and 2 (Ig1&2), a decoy protein that activates the receptor tyrosine phosphatase sigma (PTPRS) on FLS, for RA therapy. We report that PTPRS expression is enriched in synovial lining RA FLS and that Ig1&2 reduces migration of RA but not osteoarthritis FLS. Administration of an Fc-fusion Ig1&2 attenuated arthritis in mice without affecting innate or adaptive immunity. Furthermore, PTPRS was down-regulated in FLS by tumor necrosis factor (TNF) via a phosphatidylinositol 3-kinase-mediated pathway, and TNF inhibition enhanced PTPRS expression in arthritic joints. Combination of ineffective doses of TNF inhibitor and Fc-Ig1&2 reversed arthritis in mice, providing an example of synergy between FLS-targeted and immunosuppressive DMARD therapies.

A positively selected FBN1 missense variant reduces height in Peruvian individuals.

On average, Peruvian individuals are among the shortest in the world1. Here we show that Native American ancestry is associated with reduced height in an ethnically diverse group of Peruvian individuals, and identify a population-specific, missense variant in the FBN1 gene (E1297G) that is significantly associated with lower height. Each copy of the minor allele (frequency of 4.7%) reduces height by 2.2 cm (4.4 cm in homozygous individuals). To our knowledge, this is the largest effect size known for a common height-associated variant. FBN1 encodes the extracellular matrix protein fibrillin 1, which is a major structural component of microfibrils. We observed less densely packed fibrillin-1-rich microfibrils with irregular edges in the skin of individuals who were homozygous for G1297 compared with individuals who were homozygous for E1297. Moreover, we show that the E1297G locus is under positive selection in non-African populations, and that the E1297 variant shows subtle evidence of positive selection specifically within the Peruvian population. This variant is also significantly more frequent in coastal Peruvian populations than in populations from the Andes or the Amazon, which suggests that short stature might be the result of adaptation to factors that are associated with the coastal environment in Peru.

Using genetics to prioritize diagnoses for rheumatology outpatients with inflammatory arthritis.

It is challenging to quickly diagnose slowly progressing diseases. To prioritize multiple related diagnoses, we developed G-PROB (Genetic Probability tool) to calculate the probability of different diseases for a patient using genetic risk scores. We tested G-PROB for inflammatory arthritis-causing diseases (rheumatoid arthritis, systemic lupus erythematosus, spondyloarthropathy, psoriatic arthritis, and gout). After validating on simulated data, we tested G-PROB in three cohorts: 1211 patients identified by International Classification of Diseases (ICD) codes within the eMERGE database, 245 patients identified through ICD codes and medical record review within the Partners Biobank, and 243 patients first presenting with unexplained inflammatory arthritis and with final diagnoses by record review within the Partners Biobank. Calibration of G-probabilities with disease status was high, with regression coefficients from 0.90 to 1.08 (1.00 is ideal). G-probabilities discriminated true diagnoses across the three cohorts with pooled areas under the curve (95% CI) of 0.69 (0.67 to 0.71), 0.81 (0.76 to 0.84), and 0.84 (0.81 to 0.86), respectively. For all patients, at least one disease could be ruled out, and in 45% of patients, a likely diagnosis was identified with a 64% positive predictive value. In 35% of cases, the clinician's initial diagnosis was incorrect. Initial clinical diagnosis explained 39% of the variance in final disease, which improved to 51% (P < 0.0001) after adding G-probabilities. Converting genotype information before a clinical visit into an interpretable probability value for five different inflammatory arthritides could potentially be used to improve the diagnostic efficiency of rheumatic diseases in clinical practice.

CUX1 and IκBζ (NFKBIZ) mediate the synergistic inflammatory response to TNF and IL-17A in stromal fibroblasts.

The role of stromal fibroblasts in chronic inflammation is unfolding. In rheumatoid arthritis, leukocyte-derived cytokines TNF and IL-17A work together, activating fibroblasts to become a dominant source of the hallmark cytokine IL-6. However, IL-17A alone has minimal effect on fibroblasts. To identify key mediators of the synergistic response to TNF and IL-17A in human synovial fibroblasts, we performed time series, dose-response, and gene-silencing transcriptomics experiments. Here we show that in combination with TNF, IL-17A selectively induces a specific set of genes mediated by factors including cut-like homeobox 1 (CUX1) and IκBζ (NFKBIZ). In the promoters of CXCL1, CXCL2, and CXCL3, we found a putative CUX1-NF-κB binding motif not found elsewhere in the genome. CUX1 and NF-κB p65 mediate transcription of these genes independent of LIFR, STAT3, STAT4, and ELF3. Transcription of NFKBIZ, encoding the atypical IκB factor IκBζ, is IL-17A dose-dependent, and IκBζ only mediates the transcriptional response to TNF and IL-17A, but not to TNF alone. In fibroblasts, IL-17A response depends on CUX1 and IκBζ to engage the NF-κB complex to produce chemoattractants for neutrophil and monocyte recruitment.

Fast, sensitive and accurate integration of single-cell data with Harmony.

The emerging diversity of single-cell RNA-seq datasets allows for the full transcriptional characterization of cell types across a wide variety of biological and clinical conditions. However, it is challenging to analyze them together, particularly when datasets are assayed with different technologies, because biological and technical differences are interspersed. We present Harmony (https://github.com/immunogenomics/harmony), an algorithm that projects cells into a shared embedding in which cells group by cell type rather than dataset-specific conditions. Harmony simultaneously accounts for multiple experimental and biological factors. In six analyses, we demonstrate the superior performance of Harmony to previously published algorithms while requiring fewer computational resources. Harmony enables the integration of ~106 cells on a personal computer. We apply Harmony to peripheral blood mononuclear cells from datasets with large experimental differences, five studies of pancreatic islet cells, mouse embryogenesis datasets and the integration of scRNA-seq with spatial transcriptomics data.

Author Correction: Tubular cell and keratinocyte single-cell transcriptomics applied to lupus nephritis reveal type I IFN and fibrosis relevant pathways.

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

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.

Tubular cell and keratinocyte single-cell transcriptomics applied to lupus nephritis reveal type I IFN and fibrosis relevant pathways.

The molecular and cellular processes that lead to renal damage and to the heterogeneity of lupus nephritis (LN) are not well understood. We applied single-cell RNA sequencing (scRNA-seq) to renal biopsies from patients with LN and evaluated skin biopsies as a potential source of diagnostic and prognostic markers of renal disease. Type I interferon (IFN)-response signatures in tubular cells and keratinocytes distinguished patients with LN from healthy control subjects. Moreover, a high IFN-response signature and fibrotic signature in tubular cells were each associated with failure to respond to treatment. Analysis of tubular cells from patients with proliferative, membranous and mixed LN indicated pathways relevant to inflammation and fibrosis, which offer insight into their histologic differences. In summary, we applied scRNA-seq to LN to deconstruct its heterogeneity and identify novel targets for personalized approaches to therapy.

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.

Lymphocyte innateness defined by transcriptional states reflects a balance between proliferation and effector functions.

How innate T cells (ITC), including invariant natural killer T (iNKT) cells, mucosal-associated invariant T (MAIT) cells, and γδ T cells, maintain a poised effector state has been unclear. Here we address this question using low-input and single-cell RNA-seq of human lymphocyte populations. Unbiased transcriptomic analyses uncover a continuous 'innateness gradient', with adaptive T cells at one end, followed by MAIT, iNKT, γδ T and natural killer cells at the other end. Single-cell RNA-seq reveals four broad states of innateness, and heterogeneity within canonical innate and adaptive populations. Transcriptional and functional data show that innateness is characterized by pre-formed mRNA encoding effector functions, but impaired proliferation marked by decreased baseline expression of ribosomal genes. Together, our data shed new light on the poised state of ITC, in which innateness is defined by a transcriptionally-orchestrated trade-off between rapid cell growth and rapid effector function.

Mixed-effects association of single cells identifies an expanded effector CD4+ T cell subset in rheumatoid arthritis.

High-dimensional single-cell analyses have improved the ability to resolve complex mixtures of cells from human disease samples; however, identifying disease-associated cell types or cell states in patient samples remains challenging because of technical and interindividual variation. Here, we present mixed-effects modeling of associations of single cells (MASC), a reverse single-cell association strategy for testing whether case-control status influences the membership of single cells in any of multiple cellular subsets while accounting for technical confounders and biological variation. Applying MASC to mass cytometry analyses of CD4+ T cells from the blood of rheumatoid arthritis (RA) patients and controls revealed a significantly expanded population of CD4+ T cells, identified as CD27- HLA-DR+ effector memory cells, in RA patients (odds ratio, 1.7; P = 1.1 × 10-3). The frequency of CD27- HLA-DR+ cells was similarly elevated in blood samples from a second RA patient cohort, and CD27- HLA-DR+ cell frequency decreased in RA patients who responded to immunosuppressive therapy. Mass cytometry and flow cytometry analyses indicated that CD27- HLA-DR+ cells were associated with RA (meta-analysis P = 2.3 × 10-4). Compared to peripheral blood, synovial fluid and synovial tissue samples from RA patients contained about fivefold higher frequencies of CD27- HLA-DR+ cells, which comprised ~10% of synovial CD4+ T cells. CD27- HLA-DR+ cells expressed a distinctive effector memory transcriptomic program with T helper 1 (TH1)- and cytotoxicity-associated features and produced abundant interferon-γ (IFN-γ) and granzyme A protein upon stimulation. We propose that MASC is a broadly applicable method to identify disease-associated cell populations in high-dimensional single-cell data.

Discovering in vivo cytokine-eQTL interactions from a lupus clinical trial.

BACKGROUND: Cytokines are critical to human disease and are attractive therapeutic targets given their widespread influence on gene regulation and transcription. Defining the downstream regulatory mechanisms influenced by cytokines is central to defining drug and disease mechanisms. One promising strategy is to use interactions between expression quantitative trait loci (eQTLs) and cytokine levels to define target genes and mechanisms. RESULTS: In a clinical trial for anti-IL-6 in patients with systemic lupus erythematosus, we measure interferon (IFN) status, anti-IL-6 drug exposure, and whole blood genome-wide gene expression at three time points. We show that repeat transcriptomic measurements increases the number of cis eQTLs identified compared to using a single time point. We observe a statistically significant enrichment of in vivo eQTL interactions with IFN status and anti-IL-6 drug exposure and find many novel interactions that have not been previously described. Finally, we find transcription factor binding motifs interrupted by eQTL interaction SNPs, which point to key regulatory mediators of these environmental stimuli and therefore potential therapeutic targets for autoimmune diseases. In particular, genes with IFN interactions are enriched for ISRE binding site motifs, while those with anti-IL-6 interactions are enriched for IRF4 motifs. CONCLUSIONS: This study highlights the potential to exploit clinical trial data to discover in vivo eQTL interactions with therapeutically relevant environmental variables.

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.