Tissue-specific enhancer-gene maps from multimodal single-cell data identify causal disease alleles

Translating genome-wide association study (GWAS) loci into causal variants and genes requires accurate cell-type-specific enhancer-gene maps from disease-relevant tissues. Building enhancer-gene maps is essential but challenging with current experimental methods in primary human tissues. Here we developed a nonparametric statistical method, SCENT (single-cell enhancer target gene mapping), that models association between enhancer chromatin accessibility and gene expression in single-cell or nucleus multimodal RNA sequencing and ATAC sequencing data. We applied SCENT to 9 multimodal datasets including >120,000 single cells or nuclei and created 23 cell-type-specific enhancer-gene maps. These maps were highly enriched for causal variants in expression quantitative loci and GWAS for 1,143 diseases and traits. We identified likely causal genes for both common and rare diseases and linked somatic mutation hotspots to target genes. We demonstrate that application of SCENT to multimodal data from disease-relevant human tissue enables the scalable construction of accurate cell-type-specific enhancer-gene maps, essential for defining noncoding variant function.

Mapping the dynamic genetic regulatory architecture of HLA genes at single-cell resolution

The human leukocyte antigen (HLA) locus plays a critical role in complex traits spanning autoimmune and infectious diseases, transplantation and cancer. While coding variation in HLA genes has been extensively documented, regulatory genetic variation modulating HLA expression levels has not been comprehensively investigated. Here we mapped expression quantitative trait loci (eQTLs) for classical HLA genes across 1,073 individuals and 1,131,414 single cells from three tissues. To mitigate technical confounding, we developed scHLApers, a pipeline to accurately quantify single-cell HLA expression using personalized reference genomes. We identified cell-type-specific cis-eQTLs for every classical HLA gene. Modeling eQTLs at single-cell resolution revealed that many eQTL effects are dynamic across cell states even within a cell type. HLA-DQ genes exhibit particularly cell-state-dependent effects within myeloid, B and T cells. For example, a T cell HLA-DQA1 eQTL ( rs3104371 ) is strongest in cytotoxic cells. Dynamic HLA regulation may underlie important interindividual variability in immune responses.

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.

Clonally expanded CD38hi cytotoxic CD8 T cells define the T cell infiltrate in checkpoint inhibitor-associated arthritis

Immune checkpoint inhibitor (ICI) therapies used to treat cancer, such as anti-PD-1 antibodies, can induce autoimmune conditions in some individuals. The T cell mechanisms mediating such iatrogenic autoimmunity and their overlap with spontaneous autoimmune diseases remain unclear. Here, we compared T cells from the joints of 20 patients with an inflammatory arthritis induced by ICI therapy (ICI-arthritis) with two archetypal autoimmune arthritides, rheumatoid arthritis (RA) and psoriatic arthritis (PsA). Single-cell transcriptomic and antigen receptor repertoire analyses highlighted clonal expansion of an activated effector CD8 T cell population in the joints and blood of patients with ICI-arthritis. These cells were identified as CD38hiCD127- CD8 T cells and were uniquely enriched in ICI-arthritis joints compared with RA and PsA and also displayed an elevated interferon signature. In vitro, type I interferon induced CD8 T cells to acquire the ICI-associated CD38hi phenotype and enhanced cytotoxic function. In a cohort of patients with advanced melanoma, ICI therapy markedly expanded circulating CD38hiCD127- T cells, which were frequently bound by the therapeutic anti-PD-1 drug. In patients with ICI-arthritis, drug-bound CD8 T cells in circulation showed marked clonal overlap with drug-bound CD8 T cells from synovial fluid. These results suggest that ICI therapy directly targets CD8 T cells in patients who develop ICI-arthritis and induces an autoimmune pathology that is distinct from prototypical spontaneous autoimmune arthritides.

Drivers of heterogeneity in synovial fibroblasts in rheumatoid arthritis

Inflammation of non-barrier immunologically quiescent tissues is associated with a massive influx of blood-borne innate and adaptive immune cells. Cues from the latter are likely to alter and expand activated states of the resident cells. However, local communications between immigrant and resident cell types in human inflammatory disease remain poorly understood. Here, we explored drivers of fibroblast-like synoviocyte (FLS) heterogeneity in inflamed joints of patients with rheumatoid arthritis using paired single-cell RNA and ATAC sequencing, multiplexed imaging and spatial transcriptomics along with in vitro modeling of cell-extrinsic factor signaling. These analyses suggest that local exposures to myeloid and T cell-derived cytokines, TNF, IFN-γ, IL-1β or lack thereof, drive four distinct FLS states some of which closely resemble fibroblast states in other disease-affected tissues including skin and colon. Our results highlight a role for concurrent, spatially distributed cytokine signaling within the inflamed synovium.

Emerging synovial cell states in rheumatoid arthritis as potential therapeutic targets

PURPOSE OF REVIEW

To summarize recently discovered novel cell states in rheumatoid arthritis (RA) synovium that could have important implications for disease treatment.

RECENT FINDINGS

The use of multiomic technologies, including single-cell and spatial transcriptomics and mass cytometry, has led to the discovery of several novel cell states, which could have important implications for the treatment of RA. These cells can be found in patient blood, synovial fluid, or synovial tissue and span several immune cell subsets as well as stromal cell types. These diverse cell states may represent the targets of current or future therapeutics, while their fluctuations may inform the ideal timing for therapy. Future efforts are needed to implicate how each cell state functions in the pathophysiologic network within affected joints and how medications perturb each cell state and ultimately the tissue.

SUMMARY

Multiomic molecular technologies have afforded the discovery of numerous novel cellular states in RA synovium; the next challenge will be to link these states to pathophysiology and treatment response.

The Chromatin Landscape of Pathogenic Transcriptional Cell States in Rheumatoid Arthritis

Synovial tissue inflammation is the hallmark of rheumatoid arthritis (RA). Recent work has identified prominent pathogenic cell states in inflamed RA synovial tissue, such as T peripheral helper cells; however, the epigenetic regulation of these states has yet to be defined. We measured genome-wide open chromatin at single cell resolution from 30 synovial tissue samples, including 12 samples with transcriptional data in multimodal experiments. We identified 24 chromatin classes and predicted their associated transcription factors, including a CD8+ GZMK+ class associated with EOMES and a lining fibroblast class associated with AP-1. By integrating an RA tissue transcriptional atlas, we found that the chromatin classes represented 'superstates' corresponding to multiple transcriptional cell states. Finally, we demonstrated the utility of this RA tissue chromatin atlas through the associations between disease phenotypes and chromatin class abundance as well as the nomination of classes mediating the effects of putatively causal RA genetic variants.

Clonal associations of lymphocyte subsets and functional states revealed by single cell antigen receptor profiling of T and B cells in rheumatoid arthritis synovium

Rheumatoid arthritis (RA) is an autoimmune disease initiated by antigen-specific T cells and B cells, which promote synovial inflammation through a complex set of interactions with innate immune and stromal cells. To better understand the phenotypes and clonal relationships of synovial T and B cells, we performed single-cell RNA and repertoire sequencing on paired synovial tissue and peripheral blood samples from 12 donors with seropositive RA ranging from early to chronic disease. Paired transcriptomic-repertoire analyses highlighted 3 clonally distinct CD4 T cells populations that were enriched in RA synovium: T peripheral helper (Tph) and T follicular helper (Tfh) cells, CCL5+ T cells, and T regulatory cells (Tregs). Among these cells, Tph cells showed a unique transcriptomic signature of recent T cell receptor (TCR) activation, and clonally expanded Tph cells expressed an elevated transcriptomic effector signature compared to non-expanded Tph cells. CD8 T cells showed higher oligoclonality than CD4 T cells, and the largest CD8 T cell clones in synovium were highly enriched in GZMK+ cells. TCR analyses revealed CD8 T cells with likely viral-reactive TCRs distributed across transcriptomic clusters and definitively identified MAIT cells in synovium, which showed transcriptomic features of TCR activation. Among B cells, non-naive B cells including age-associated B cells (ABC), NR4A1+ activated B cells, and plasma cells, were enriched in synovium and had higher somatic hypermutation rates compared to blood B cells. Synovial B cells demonstrated substantial clonal expansion, with ABC, memory, and activated B cells clonally linked to synovial plasma cells. Together, these results reveal clonal relationships between functionally distinct lymphocyte populations that infiltrate RA synovium.

Mapping the dynamic genetic regulatory architecture of HLA genes at single-cell resolution

The human leukocyte antigen (HLA) locus plays a critical role in complex traits spanning autoimmune and infectious diseases, transplantation, and cancer. While coding variation in HLA genes has been extensively documented, regulatory genetic variation modulating HLA expression levels has not been comprehensively investigated. Here, we mapped expression quantitative trait loci (eQTLs) for classical HLA genes across 1,073 individuals and 1,131,414 single cells from three tissues, using personalized reference genomes to mitigate technical confounding. We identified cell-type-specific cis-eQTLs for every classical HLA gene. Modeling eQTLs at single-cell resolution revealed that many eQTL effects are dynamic across cell states even within a cell type. HLA-DQ genes exhibit particularly cell-state-dependent effects within myeloid, B, and T cells. Dynamic HLA regulation may underlie important interindividual variability in immune responses.

Cytotoxic CD8+ T cells target citrullinated antigens in rheumatoid arthritis

The immune mechanisms that mediate synovitis and joint destruction in rheumatoid arthritis (RA) remain poorly defined. Although increased levels of CD8⁺ T cells have been described in RA, their function in pathogenesis remains unclear. Here we perform single cell transcriptome and T cell receptor (TCR) sequencing of CD8⁺ T cells derived from anti-citrullinated protein antibodies (ACPA)+ RA blood. We identify GZMB⁺CD8⁺ subpopulations containing large clonal lineage expansions that express cytotoxic and tissue homing transcriptional programs, while a GZMK⁺CD8⁺ memory subpopulation comprises smaller clonal expansions that express effector T cell transcriptional programs. We demonstrate RA citrullinated autoantigens presented by MHC class I activate RA blood-derived GZMB⁺CD8⁺ T cells to expand, express cytotoxic mediators, and mediate killing of target cells. We also demonstrate that these clonally expanded GZMB⁺CD8⁺ cells are present in RA synovium. These findings suggest that cytotoxic CD8⁺ T cells targeting citrullinated antigens contribute to synovitis and joint tissue destruction in ACPA+ RA.

Neoplastic synovial lining cells that coexpress podoplanin and CD90 overproduce CSF-1, driving tenosynovial giant cell tumor

Tenosynovial giant cell tumor (TCGT) is a rare neoplasm affecting the synovium of joints, bursae, and tendon sheaths. The overproduction of colony-stimulating factor-1 (CSF-1) by a minority of the tumor population works in a paracrine fashion to drive tumor growth. Pathology of the reactive, monocytic component has been well elucidated, whereas the populations of neoplastic cells and all the sources of CSF-1 overproduction are incompletely characterized. Podoplanin (PDPN), or gp38, is a cell surface glycoprotein that is expressed on fibroblast-like synovial cells and upregulated in rheumatoid arthritis and many cancers; it governs cell mobility, epithelial-mesenchymal transition, and other functions and is associated with lymphangiogenesis and poor prognosis in many solid tumors, which underscores its local and possible systemic effects. We found higher PDPN expression in TGCT than in internal controls of patients' healthy synovium. Flow cytometry partitioned PDPN^high cells into PDPN^high CD90⁺ and PDPN^high CD14⁺ populations. Quantitative real-time polymerase chain reaction analysis of the PDPN^high CD90⁺ cells revealed that CSF-1 expression was 10-fold higher than in PDPN^high CD14⁺ cells. Therefore, we conclude that the lining fibroblast-like synovial cells, which express PDPN^high CD90⁺ , are responsible for the overproduction of CSF-1 and for driving tumor growth.

Dynamic spectrum of ectopic lymphoid B cell activation and hypermutation in the RA synovium characterized by NR4A nuclear receptor expression

Ectopic lymphoid structures (ELS) can develop in rheumatoid arthritis (RA) synovial tissue, but the precise pathways of B cell activation and selection are not well understood. Here, we identify a synovial B cell population characterized by co-expression of a family of orphan nuclear receptors (NR4A1-3), which is highly enriched in RA synovial tissue. A transcriptomic profile of NR4A synovial B cells significantly overlaps with germinal center light zone B cells and an accrual of somatic hypermutation that correlates with loss of naive B cell state. NR4A B cells co-express lymphotoxins α and β and IL-6, supporting functions in ELS promotion. Expanded and shared clones between synovial NR4A B cells and plasma cells and the rapid upregulation with BCR stimulation point to in situ differentiation. Together, we identify a dynamic progression of B cell activation in RA synovial ELS, with NR4A transcription factors having an important role in local adaptive immune responses.

Sequencing of Circulating Microbial Cell-Free DNA Can Identify Pathogens in Periprosthetic Joint Infections

BACKGROUND

Over 1 million Americans undergo joint replacement each year, and approximately 1 in 75 will incur a periprosthetic joint infection. Effective treatment necessitates pathogen identification, yet standard-of-care cultures fail to detect organisms in 10% to 20% of cases and require invasive sampling. We hypothesized that cell-free DNA (cfDNA) fragments from microorganisms in a periprosthetic joint infection can be found in the bloodstream and utilized to accurately identify pathogens via next-generation sequencing.

METHODS

In this prospective observational study performed at a musculoskeletal specialty hospital in the U.S., we enrolled 53 adults with validated hip or knee periprosthetic joint infections. Participants had peripheral blood drawn immediately prior to surgical treatment. Microbial cfDNA from plasma was sequenced and aligned to a genome database with >1,000 microbial species. Intraoperative tissue and synovial fluid cultures were performed per the standard of care. The primary outcome was accuracy in organism identification with use of blood cfDNA sequencing, as measured by agreement with tissue-culture results.

RESULTS

Intraoperative and preoperative joint cultures identified an organism in 46 (87%) of 53 patients. Microbial cfDNA sequencing identified the joint pathogen in 35 cases, including 4 of 7 culture-negative cases (57%). Thus, as an adjunct to cultures, cfDNA sequencing increased pathogen detection from 87% to 94%. The median time to species identification for cases with genus-only culture results was 3 days less than standard-of-care methods. Circulating cfDNA sequencing in 14 cases detected additional microorganisms not grown in cultures. At postoperative encounters, cfDNA sequencing demonstrated no detection or reduced levels of the infectious pathogen.

CONCLUSIONS

Microbial cfDNA from pathogens causing local periprosthetic joint infections can be detected in peripheral blood. These circulating biomarkers can be sequenced from noninvasive venipuncture, providing a novel source for joint pathogen identification. Further development as an adjunct to tissue cultures holds promise to increase the number of cases with accurate pathogen identification and improve time-to-speciation. This test may also offer a novel method to monitor infection clearance during the treatment period.

LEVEL OF EVIDENCE

Diagnostic Level II. See Instructions for Authors for a complete description of levels of evidence.

IFN-γ and TNF-α drive a CXCL10+ CCL2+ macrophage phenotype expanded in severe COVID-19 lungs and inflammatory diseases with tissue inflammation

BACKGROUND

Immunosuppressive and anti-cytokine treatment may have a protective effect for patients with COVID-19. Understanding the immune cell states shared between COVID-19 and other inflammatory diseases with established therapies may help nominate immunomodulatory therapies.

METHODS

To identify cellular phenotypes that may be shared across tissues affected by disparate inflammatory diseases, we developed a meta-analysis and integration pipeline that models and removes the effects of technology, tissue of origin, and donor that confound cell-type identification. Using this approach, we integrated > 300,000 single-cell transcriptomic profiles from COVID-19-affected lungs and tissues from healthy subjects and patients with five inflammatory diseases: rheumatoid arthritis (RA), Crohn's disease (CD), ulcerative colitis (UC), systemic lupus erythematosus (SLE), and interstitial lung disease. We tested the association of shared immune states with severe/inflamed status compared to healthy control using mixed-effects modeling. To define environmental factors within these tissues that shape shared macrophage phenotypes, we stimulated human blood-derived macrophages with defined combinations of inflammatory factors, emphasizing in particular antiviral interferons IFN-beta (IFN-β) and IFN-gamma (IFN-γ), and pro-inflammatory cytokines such as TNF.

RESULTS

We built an immune cell reference consisting of > 300,000 single-cell profiles from 125 healthy or disease-affected donors from COVID-19 and five inflammatory diseases. We observed a CXCL10+ CCL2+ inflammatory macrophage state that is shared and strikingly abundant in severe COVID-19 bronchoalveolar lavage samples, inflamed RA synovium, inflamed CD ileum, and UC colon. These cells exhibited a distinct arrangement of pro-inflammatory and interferon response genes, including elevated levels of CXCL10, CXCL9, CCL2, CCL3, GBP1, STAT1, and IL1B. Further, we found this macrophage phenotype is induced upon co-stimulation by IFN-γ and TNF-α.

CONCLUSIONS

Our integrative analysis identified immune cell states shared across inflamed tissues affected by inflammatory diseases and COVID-19. Our study supports a key role for IFN-γ together with TNF-α in driving an abundant inflammatory macrophage phenotype in severe COVID-19-affected lungs, as well as inflamed RA synovium, CD ileum, and UC colon, which may be targeted by existing immunomodulatory therapies.

Rheumatology in the era of precision medicine: synovial tissue molecular patterns and treatment response in rheumatoid arthritis

PURPOSE OF REVIEW

A critical unmet need in rheumatoid arthritis (RA) is the identification of biomarkers that predict which of the available medications will be most effective for an individual in order to lower disease activity sooner than is afforded by the current treat-to-target approach. Here we will discuss recent reports examining the potential for synovial tissue molecular, cellular, and spatial profiling in defining objective measures of treatment response and therein developing personalized medicine for RA.

RECENT FINDINGS

Recent high-dimensional molecular profiling of RA synovium has provided unprecedented resolution of the cell types and pathways in tissues affected by rheumatic diseases. Heightened attention to tissue architecture is also emerging as a means to classify individual disease variation that may allow patients to be further stratified by therapeutic response. Although this wealth of data may have already pinpointed promising biomarkers, additional studies, likely including tissue-based functional drug response assays, will be required to demonstrate how the complex tissue environment responds.

SUMMARY

Molecular, cellular, and more recently spatial profiling of the RA synovium are uncovering fundamental features of the disease. Current investigations are examining whether this information will provide meaningful biomarkers for individualized medicine in RA.

IFN- γ and TNF- α drive a CXCL10 + CCL2 + macrophage phenotype expanded in severe COVID-19 and other diseases with tissue inflammation

Immunosuppressive and anti-cytokine treatment may have a protective effect for patients with COVID-19. Understanding the immune cell states shared between COVID-19 and other inflammatory diseases with established therapies may help nominate immunomodulatory therapies. Using an integrative strategy, we built a reference by meta-analyzing > 300,000 immune cells from COVID-19 and 5 inflammatory diseases including rheumatoid arthritis (RA), Crohn's disease (CD), ulcerative colitis (UC), lupus, and interstitial lung disease. Our cross-disease analysis revealed that an FCN1 + inflammatory macrophage state is common to COVID-19 bronchoalveolar lavage samples, RA synovium, CD ileum, and UC colon. We also observed that a CXCL10 + CCL2 + inflammatory macrophage state is abundant in severe COVID-19, inflamed CD and RA, and expresses inflammatory genes such as GBP1, STAT1 , and IL1B . We found that the CXCL10 + CCL2 + macrophages are transcriptionally similar to blood-derived macrophages stimulated with TNF- α and IFN- γ ex vivo . Our findings suggest that IFN- γ , alongside TNF- α , might be a key driver of this abundant inflammatory macrophage phenotype in severe COVID-19 and other inflammatory diseases, which may be targeted by existing immunomodulatory therapies.

Rheumatoid Arthritis Flares After Total Hip and Total Knee Arthroplasty: Outcomes at One Year

OBJECTIVE

Most patients with rheumatoid arthritis (RA) undergoing total hip arthroplasty (THA) and total knee arthroplasty (TKA) have active RA and report postoperative flares; whether RA disease activity or flares increase the risk of worse pain and function scores 1 year later is unknown.

METHODS

Patients with RA were enrolled before THA/TKA. Patient-reported outcomes, including the Hip disability and Osteoarthritis Outcome Score (HOOS)/Knee Injury and Osteoarthritis Outcome Score (KOOS) and physician assessments of disease characteristics and activity (Disease Activity Score in 28 joints [DAS28] and Clinical Disease Activity Index), were collected before surgery. Patient-reported outcomes were repeated at 1 year. Postoperative flares were identified using the RA Flare Questionnaire weekly for 6 weeks and were defined by concordance between patient report plus physician assessment. We compared baseline characteristics and HOOS/KOOS scores using 2-sample t-test/Wilcoxon's rank sum test as well as chi-square/Fisher's exact tests. We used multivariate linear and logistic regression to determine the association of baseline characteristics, disease activity, and flares with 1-year outcomes.

RESULTS

One-year HOOS/KOOS scores were available for 122 patients (56 with THA and 66 with TKA). Although HOOS/KOOS pain was worse for patients who experienced a flare within 6 weeks of surgery, absolute improvement was not different. In multivariable models, baseline DAS28 predicted 1-year HOOS/KOOS pain and function; each 1-unit increase in DAS28 worsened 1-year pain by 2.41 (SE 1.05; P = 0.02) and 1-year function by 4.96 (SE 1.17; P = 0.0001). Postoperative flares were not independent risk factors for pain or function scores.

CONCLUSION

Higher disease activity increased the risk of worse pain and function 1 year after arthroplasty, but postoperative flares did not.

Fibroblast subtypes in tissues affected by autoimmunity: with lessons from lymph node fibroblasts

The recent advent of single-cell technologies has fast-tracked the discovery of multiple fibroblast subsets in tissues affected by autoimmune disease. In recent years, interest in lymph node fibroblasts that support and regulate immune cells has also grown, leading to an expanding framework of stromal cell subsets with distinct spatial, transcriptional, and functional characteristics. Inflammation can drive tissue fibroblasts to adopt a lymphoid tissue stromal cell phenotype, suggesting that fibroblasts in diseased tissues can have counterparts in lymphoid tissues. Here, we examine fibroblast subsets in tissues affected by autoimmunity in the context of knowledge gained from studies on lymph node fibroblasts, with the ultimate aim to better understand stromal cell heterogeneity in these immunologically reactive tissues.

Notch signalling drives synovial fibroblast identity and arthritis pathology

The synovium is a mesenchymal tissue composed mainly of fibroblasts with a lining and sublining that surrounds the joints. In rheumatoid arthritis (RA), the synovial tissue undergoes marked hyperplasia, becomes inflamed and invasive and destroys the joint^1,2. Recently, we and others found that a subset of fibroblasts located in the sublining undergoes major expansion in RA and is linked to disease activity^3,4,5. However, the molecular mechanism by which these fibroblasts differentiate and expand in RA remains unknown. Here, we identified a critical role for NOTCH3 signaling in the differentiation of perivascular and sublining CD90(THY1)+ fibroblasts. Using single cell RNA-sequencing and synovial tissue organoids, we found that NOTCH3 signaling drives both transcriptional and spatial gradients in fibroblasts emanating from vascular endothelial cells outward. In active RA, NOTCH3 and NOTCH target genes are markedly upregulated in synovial fibroblasts. Importantly, genetic deletion of Notch3 or monoclonal antibody-blockade of NOTCH3 signaling attenuates inflammation and prevents joint damage in inflammatory arthritis. Our results indicate that synovial fibroblasts exhibit positional identity regulated by endothelium-derived Notch signaling and that this stromal crosstalk pathway underlies inflammation and pathology in inflammatory arthritis.

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 (T_H1)- 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.

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.

Histologic and Transcriptional Evidence of Subclinical Synovial Inflammation in Patients With Rheumatoid Arthritis in Clinical Remission

OBJECTIVE

Patients with rheumatoid arthritis (RA) in clinical remission may have subclinical synovial inflammation. This study was undertaken to determine the proportion of patients with RA in remission or with low disease activity at the time of arthroplasty who had histologic or transcriptional evidence of synovitis, and to identify clinical features that distinguished patients as having subclinical synovitis.

METHODS

We compared Disease Activity Score in 28 joints (DAS28) to synovial histologic features in 135 patients with RA undergoing arthroplasty. We also compared DAS28 scores to RNA-Seq data in a subset of 35 patients.

RESULTS

Fourteen percent of patients met DAS28 criteria for clinical remission (DAS28 <2.6), and another 15% met criteria for low disease activity (DAS28 <3.2). Histologic analysis of synovium revealed synovitis in 27% and 31% of samples from patients in remission and patients with low disease activity, respectively. Patients with low disease activity and synovitis also exhibited increased C-reactive protein (CRP) (P = 0.0006) and increased anti-cyclic citrullinated peptide (anti-CCP) antibody levels (P = 0.03) compared to patients without synovitis. Compared to patients with a "low inflammatory synovium" subtype, 183 genes were differentially expressed in the synovium of patients with subclinical synovitis. The majority of these genes (86%) were also differentially expressed in the synovium of patients with clinically active disease (DAS28 ≥3.2).

CONCLUSION

Thirty-one percent of patients with low clinical disease activity exhibited histologic evidence of subclinical synovitis, which was associated with increased CRP and anti-CCP levels. Our findings suggest that synovial gene expression signatures of clinical synovitis are present in patients with subclinical synovitis.

Insights into rheumatic diseases from next-generation sequencing

Rheumatic diseases have complex aetiologies that are not fully understood, which makes the study of pathogenic mechanisms in these diseases a challenge for researchers. Next-generation sequencing (NGS) and related omics technologies, such as transcriptomics, epigenomics and genomics, provide an unprecedented genome-wide view of gene expression, environmentally responsive epigenetic changes and genetic variation. The integrated application of NGS technologies to samples from carefully phenotyped clinical cohorts of patients has the potential to solve remaining mysteries in the pathogenesis of several rheumatic diseases, to identify new therapeutic targets and to underpin a precision medicine approach to the diagnosis and treatment of rheumatic diseases. This Review provides an overview of the NGS technologies available, showcases important advances in rheumatic disease research already powered by these technologies and highlights NGS approaches that hold particular promise for generating new insights and advancing the field.

HBEGF+ macrophages in rheumatoid arthritis induce fibroblast invasiveness

Macrophages tailor their function according to the signals found in tissue microenvironments, assuming a wide spectrum of phenotypes. A detailed understanding of macrophage phenotypes in human tissues is limited. Using single-cell RNA sequencing, we defined distinct macrophage subsets in the joints of patients with the autoimmune disease rheumatoid arthritis (RA), which affects ~1% of the population. The subset we refer to as HBEGF+ inflammatory macrophages is enriched in RA tissues and is shaped by resident fibroblasts and the cytokine tumor necrosis factor (TNF). These macrophages promoted fibroblast invasiveness in an epidermal growth factor receptor-dependent manner, indicating that intercellular cross-talk in this inflamed setting reshapes both cell types and contributes to fibroblast-mediated joint destruction. In an ex vivo synovial tissue assay, most medications used to treat RA patients targeted HBEGF+ inflammatory macrophages; however, in some cases, medication redirected them into a state that is not expected to resolve inflammation. These data highlight how advances in our understanding of chronically inflamed human tissues and the effects of medications therein can be achieved by studies on local macrophage phenotypes and intercellular interactions.

Increased Staphylococcus aureus Nasal Carriage Rates in Rheumatoid Arthritis Patients on Biologic Therapy

BACKGROUND

Rheumatoid arthritis patients are at increased risk for periprosthetic joint infection after arthroplasty. The reason is multifactorial. Nasal colonization with Staphylococcus aureus is a modifiable risk factor; carriage rates in RA patients are unknown. The goal of this study is to determine the S aureus nasal carriage rates of RA patients on biologics, RA patients on traditional disease-modifying anti-rheumatic drugs (DMARDs), and osteoarthritis.

METHODS

Consecutive patients with RA on biologics (±DMARDs), RA on non-biologic DMARDs, or OA were prospectively enrolled from April 2017 to May 2018. One hundred twenty-three patients were determined necessary per group to show a difference in carriage rates. Patients underwent a nasal swab and answered questions to identify additional risk factors. S aureus positive swabs were further categorized using spa typing. Logistic regression evaluated the association with S aureus colonization between the groups after controlling for known risk factors.

RESULTS

RA patients on biologics, 70% of whom were on DMARDs, had statistically significant increase in S aureus colonization (37%) compared to RA on DMARDs alone (24%), or OA (20%) (P = .01 overall). After controlling for glucocorticoids, antibiotic use, recent hospitalization, and diabetes, RA on biologics had a significant increased risk of S aureus nasal colonization (Odds ratio 1.80, 95% confidence interval 1.00-3.22, P = .047).

CONCLUSION

S aureus colonization risk was increased for RA on biologics compared to RA not on biologics and OA. Nasal S aureus carriage increases the risk of surgical site infection; this modifiable risk factor should be addressed prior to total joint arthroplasty for this higher risk patient group.

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.

Electronic supplementary material

The online version of this article (10.1186/s13075-018-1631-y) contains supplementary material, which is available to authorized users.

Identification of Three Rheumatoid Arthritis Disease Subtypes by Machine Learning Integration of Synovial Histologic Features and RNA Sequencing Data

OBJECTIVE

In this study, we sought to refine histologic scoring of rheumatoid arthritis (RA) synovial tissue by training with gene expression data and machine learning.

METHODS

Twenty histologic features were assessed in 129 synovial tissue samples (n = 123 RA patients and n = 6 osteoarthritis [OA] patients). Consensus clustering was performed on gene expression data from a subset of 45 synovial samples. Support vector machine learning was used to predict gene expression subtypes, using histologic data as the input. Corresponding clinical data were compared across subtypes.

RESULTS

Consensus clustering of gene expression data revealed 3 distinct synovial subtypes, including a high inflammatory subtype characterized by extensive infiltration of leukocytes, a low inflammatory subtype characterized by enrichment in pathways including transforming growth factor β, glycoproteins, and neuronal genes, and a mixed subtype. Machine learning applied to histologic features, with gene expression subtypes serving as labels, generated an algorithm for the scoring of histologic features. Patients with the high inflammatory synovial subtype exhibited higher levels of markers of systemic inflammation and autoantibodies. C-reactive protein (CRP) levels were significantly correlated with the severity of pain in the high inflammatory subgroup but not in the others.

CONCLUSION

Gene expression analysis of RA and OA synovial tissue revealed 3 distinct synovial subtypes. These labels were used to generate a histologic scoring algorithm in which the histologic scores were found to be associated with parameters of systemic inflammation, including the erythrocyte sedimentation rate, CRP level, and autoantibody levels. Comparison of gene expression patterns to clinical features revealed a potentially clinically important distinction: mechanisms of pain may differ in patients with different synovial subtypes.

Flares in Patients with Rheumatoid Arthritis after Total Hip and Total Knee Arthroplasty: Rates, Characteristics, and Risk Factors

OBJECTIVE

Rates of total knee arthroplasty (TKA) and total hip arthroplasty (THA) remain high for patients with rheumatoid arthritis (RA), who are at risk of flaring after surgery. We aimed to describe rates, characteristics, and risk factors of RA flare within 6 weeks of THA and TKA.

METHODS

Patients with RA were recruited prior to elective THA and TKA surgery and prospectively followed. Clinicians evaluated RA clinical characteristics 0-2 weeks before and 6 weeks after surgery. Patients answered questions regarding disease activity including self-reported joint counts and flare status weekly for 6 weeks. Per standard of care, biologics were stopped before surgery, while glucocorticoids and methotrexate (MTX) were typically continued. Multivariable logistic regression was used to identify baseline characteristics associated with postsurgical RA flares.

RESULTS

Of 120 patients, the mean age was 62 years and the median RA duration 14.8 years. Ninety-eight (82%) met 2010/1987 American College of Rheumatology/European League Against Rheumatism criteria, 53 (44%) underwent THA (and the rest TKA), and 61 (51%) were taking biologics. By 6 weeks, 75 (63%) had flared. At baseline, flarers had significantly higher disease activity (as measured by the 28-joint Disease Activity Score), erythrocyte sedimentation rate, C-reactive protein, and pain. Numerically more flarers used biologics, but stopping biologics did not predict flares, and continuing MTX was not protective. A higher baseline disease activity predicted flaring by 6 weeks (OR 2.12, p = 0.02).

CONCLUSION

Flares are frequent in patients with RA undergoing arthroplasty. Higher baseline disease activity significantly increases the risk. Although more patients stopping biologics flared, this did not independently predict flaring. The effect of early postsurgery flares requires further study.

Functionally distinct disease-associated fibroblast subsets in rheumatoid arthritis

Fibroblasts regulate tissue homeostasis, coordinate inflammatory responses, and mediate tissue damage. In rheumatoid arthritis (RA), synovial fibroblasts maintain chronic inflammation which leads to joint destruction. Little is known about fibroblast heterogeneity or if aberrations in fibroblast subsets relate to pathology. Here, we show functional and transcriptional differences between fibroblast subsets from human synovial tissues using bulk transcriptomics of targeted subpopulations and single-cell transcriptomics. We identify seven fibroblast subsets with distinct surface protein phenotypes, and collapse them into three subsets by integrating transcriptomic data. One fibroblast subset, characterized by the expression of proteins podoplanin, THY1 membrane glycoprotein and cadherin-11, but lacking CD34, is threefold expanded in patients with RA relative to patients with osteoarthritis. These fibroblasts localize to the perivascular zone in inflamed synovium, secrete proinflammatory cytokines, are proliferative, and have an in vitro phenotype characteristic of invasive cells. Our strategy may be used as a template to identify pathogenic stromal cellular subsets in other complex diseases.

Synovial fibroblasts are thought to be central mediators of joint destruction in rheumatoid arthritis (RA). Here the authors use single-cell transcriptomics and flow cytometry to identify synovial fibroblast subsets that are expanded and display distinct tissue distribution and function in patients with RA.

Modulation of TNF-induced macrophage polarization by synovial fibroblasts

Mesenchymal stromal cells have emerged as powerful modulators of the immune system. In this study, we explored how the human macrophage response to TNF is regulated by human synovial fibroblasts, the representative stromal cell type in the synovial lining of joints that become activated during inflammatory arthritis. We found that synovial fibroblasts strongly suppressed TNF-mediated induction of an IFN-β autocrine loop and downstream expression of IFN-stimulated genes (ISGs), including chemokines CXCL9 and CXCL10 that are characteristic of classical macrophage activation. TNF induced the production of soluble synovial fibroblast factors that suppressed the macrophage production of IFN-β, and cooperated with TNF to limit the responsiveness of macrophages to IFN-β by suppressing activation of Jak-STAT signaling. Genome-wide transcriptome analysis showed that cocultured synovial fibroblasts modulate the expression of approximately one third of TNF-regulated genes in macrophages, including genes in pathways important for macrophage survival and polarization toward an alternatively activated phenotype. Pathway analysis revealed that gene expression programs regulated by synovial fibroblasts in our coculture system were also regulated in rheumatoid arthritis synovial macrophages, suggesting that these fibroblast-mediated changes may contribute to rheumatoid arthritis pathogenesis. This work furthers our understanding of the interplay between innate immune and stromal cells during an inflammatory response, one that is particularly relevant to inflammatory arthritis. Our findings also identify modulation of macrophage phenotype as a new function for synovial fibroblasts that may prove to be a contributing factor in arthritis pathogenesis.

Regulation of type I interferon responses

Type I interferons (IFNs) activate intracellular antimicrobial programmes and influence the development of innate and adaptive immune responses. Canonical type I IFN signalling activates the Janus kinase (JAK)-signal transducer and activator of transcription (STAT) pathway, leading to transcription of IFN-stimulated genes (ISGs). Host, pathogen and environmental factors regulate the responses of cells to this signalling pathway and thus calibrate host defences while limiting tissue damage and preventing autoimmunity. Here, we summarize the signalling and epigenetic mechanisms that regulate type I IFN-induced STAT activation and ISG transcription and translation. These regulatory mechanisms determine the biological outcomes of type I IFN responses and whether pathogens are cleared effectively or chronic infection or autoimmune disease ensues.

Smyd2 controls cytoplasmic lysine methylation of Hsp90 and myofilament organization

Protein lysine methylation is one of the most widespread post-translational modifications in the nuclei of eukaryotic cells. Methylated lysines on histones and nonhistone proteins promote the formation of protein complexes that control gene expression and DNA replication and repair. In the cytoplasm, however, the role of lysine methylation in protein complex formation is not well established. Here we report that the cytoplasmic protein chaperone Hsp90 is methylated by the lysine methyltransferase Smyd2 in various cell types. In muscle, Hsp90 methylation contributes to the formation of a protein complex containing Smyd2, Hsp90, and the sarcomeric protein titin. Deficiency in Smyd2 results in the loss of Hsp90 methylation, impaired titin stability, and altered muscle function. Collectively, our data reveal a cytoplasmic protein network that employs lysine methylation for the maintenance and function of skeletal muscle.

Regulation of medullary thymic epithelial cell differentiation and function by the signaling protein Sin

Medullary thymic epithelial cells (mTECs) play an important role in T cell tolerance and prevention of autoimmunity. Mice deficient in expression of the signaling protein Sin exhibit exaggerated immune responses and multitissue inflammation. Here, we show that Sin is expressed in the thymic stroma, specifically in mTECs. Sin deficiency led to thymic stroma-dependent autoimmune manifestations shown by radiation chimeras and thymic transplants in nude mice, and associated with defective mTEC-mediated elimination of thymocytes in a T cell receptor transgenic model of negative selection. Lack of Sin expression correlated with a disorganized medullary architecture and fewer functionally mature mTECs under steady-state conditions. Additionally, Sin deficiency inhibited the expansion of mTECs in response to in vivo administration of keratinocyte growth factor (KGF). These results identify Sin as a novel regulator of mTEC development and T cell tolerance, and suggest that Sin is important for homeostatic maintenance of the medullary epithelium in the adult thymus.

Deficiency in expression of the signaling protein Sin/Efs leads to T-lymphocyte activation and mucosal inflammation

Our studies have concentrated on elucidating the role of the signaling protein Sin in T-lymphocyte function. We have previously shown that Sin overexpression inhibits T-lymphocyte development and activation. Here we show that Sin-deficient mice exhibit exaggerated immune responses characterized by enhanced cytokine secretion and T-cell-dependent antibody production. Excessive T-cell responses in young mice correlate with spontaneous development of inflammatory lesions in different organs of aged Sin(-/-) mice, particularly the small intestine. The intestinal inflammation is characterized by T- and B-cell infiltrates in the lamina propria, which correlate with crypt enlargement and marked villus expansion and/or damage. Similar to the human intestinal inflammatory disorder Crohn's disease (CD), and in contrast to most mouse models of mucosal inflammation, inflammatory lesions in the gastrointestinal tract of Sin(-/-) mice are restricted to the small bowel. Taken together, these results suggest that Sin regulates immune system and T-lymphocyte function and that immune system dysfunction in the absence of Sin may underlie the pathogenesis of tissue-specific inflammation and enteropathies such as CD.

The adapter molecule Sin regulates T-cell-receptor-mediated signal transduction by modulating signaling substrate availability

Engagement of the T-cell receptor (TCR) results in the activation of a multitude of signaling events that regulate the function of T lymphocytes. These signaling events are in turn modulated by adapter molecules, which control the final functional output through the formation of multiprotein complexes. In this report, we identified the adapter molecule Sin as a new regulator of T-cell activation. We found that the expression of Sin in transgenic T lymphocytes and Jurkat T cells inhibited interleukin-2 expression and T-cell proliferation. This inhibitory effect was specific and was due to defective phospholipase C-gamma (PLC-gamma) phosphorylation and activation. In contrast to other adapters that become phosphorylated upon TCR stimulation, Sin was constitutively phosphorylated in resting cells by the Src kinase Fyn and bound to signaling intermediates, including PLC-gamma. In stimulated cells, Sin was transiently dephosphorylated, which coincided with transient dissociation of Fyn and PLC-gamma. Downregulation of Sin expression using Sin-specific short interfering RNA oligonucleotides inhibited transcriptional activation in response to TCR stimulation. Our results suggest that endogenous Sin influences T-lymphocyte signaling by sequestering signaling substrates and regulating their availability and/or activity in resting cells, while Sin is required for targeting these intermediates to the TCR for fast signal transmission during stimulation.

Sin: good or bad? A T lymphocyte perspective

Stimulation of T cells through their antigen receptor induces a multitude of signaling networks that regulate T cell activation in the form of cytokine production and T cell proliferation. Multiple signal integration sites exist along these pathways in the form of multiprotein signaling complexes, the formation of which is facilitated by adapter and scaffold molecules. In recent years a number of adapter and scaffold molecules have been described in T cells and shown to play an integral part in T cell function. Among these molecules are proteins that function as positive or negative regulators of T cell activation downstream of the activated T cell receptor (TCR). Here, we discuss the role of a small family of multiadapter proteins on T cell activation, the p130Cas family, with emphasis on one of its members, Sin (Src-interacting protein). Our results suggest that Sin inhibits thymocyte development and T cell activation and is a novel negative regulator of T lymphocyte function.

Defective thymocyte maturation by transgenic expression of a truncated form of the T lymphocyte adapter molecule and Fyn substrate, Sin

Adapter molecules that promote protein-protein interactions play a central role in T lymphocyte differentiation and activation. In this study, we examined the role of the T lymphocyte-expressed adapter protein and Src kinase substrate, Sin, on thymocyte function using transgenic mice expressing an activated, truncated allele of Sin (SinDeltaC). We found that SinDeltaC expression led to reduced numbers of CD4(+) and CD8(+) single-positive cells and reduced thymic cellularity due to increased thymocyte apoptosis. Because the adapter properties of Sin are mediated by tyrosine-based motifs and given that Sin is a substrate for Src tyrosine kinases, we examined the involvement of these kinases in the inhibitory effects of SinDeltaC. We found that in transgenic thymocytes, SinDeltaC was constitutively phosphorylated by the Src kinase Fyn, but not by the related kinase Lck. Using SinDeltaC and fyn(-/-) animals, we also found that the expression of Fyn was required for the inhibitory effect of SinDeltaC on thymocyte apoptosis but not for SinDeltaC-mediated inhibition of T cell maturation. The inhibitory effect of SinDeltaC on thymocyte maturation correlated with defective activation of the mitogen-activated protein kinase extracellular signal-regulated kinase. Our results suggest that the Sin mutant inhibits thymocyte differentiation through Fyn-dependent and -independent mechanisms and that endogenous Sin may be an important regulator of thymocyte development.

Identification of diverse nerve growth factor-regulated genes by serial analysis of gene expression (SAGE) profiling

Neurotrophic factors such as nerve growth factor (NGF) promote a wide variety of responses in neurons, including differentiation, survival, plasticity, and repair. Such actions often require changes in gene expression. To identify the regulated genes and thereby to more fully understand the NGF mechanism, we carried out serial analysis of gene expression (SAGE) profiling of transcripts derived from rat PC12 cells before and after NGF-promoted neuronal differentiation. Multiple criteria supported the reliability of the profile. Approximately 157,000 SAGE tags were analyzed, representing at least 21,000 unique transcripts. Of these, nearly 800 were regulated by 6-fold or more in response to NGF. Approximately 150 of the regulated transcripts have been matched to named genes, the majority of which were not previously known to be NGF-responsive. Functional categorization of the regulated genes provides insight into the complex, integrated mechanism by which NGF promotes its multiple actions. It is anticipated that as genomic sequence information accrues the data derived here will continue to provide information about neurotrophic factor mechanisms.