Search

Innate T cells, including CD1d-restricted invariant natural killer T (iNKT) cells, are characterized by their rapid activation in response to non-peptide antigens, such as lipids. While the transcriptional profiles of naive, effector and memory adaptive T cells have been well studied, less is known about transcriptional regulation of different iNKT cell activation states. Here, using single cell RNA-sequencing, we performed longitudinal profiling of activated iNKT cells, generating a transcriptomic atlas of iNKT cell activation states. We found that transcriptional signatures of activation are highly conserved among heterogeneous iNKT cell populations, including NKT1, NKT2 and NKT17 subsets, and human iNKT cells. Strikingly, we found that regulatory iNKT cells, such as adipose iNKT cells, undergo blunted activation, and display constitutive enrichment of memory-like cMAF+ and KLRG1+ populations. Moreover, we identify a conserved cMAF-associated transcriptional network among NKT10 cells, providing novel insights into the biology of regulatory and antigen experienced iNKT cells.

Immune checkpoint inhibitor (ICI) therapies that promote T cell activation have improved outcomes for advanced malignancies yet also elicit harmful autoimmune reactions. The T cell mechanisms mediating these iatrogenic autoimmune events remain unclear. Here we assayed T cells from joints of patients affected by ICI-induced inflammatory arthritis (ICI-arthritis), which can present clinically indistinguishable from rheumatoid arthritis (RA). Compared to the autoimmune arthritides RA and psoriatic arthritis (PsA), ICI-arthritis joints contained an expanded CD38hi CD127− CD8+ T cell subset that displays cytotoxic, effector, and interferon (IFN) response signatures. The abundance of CD38hi CD8 T cells in ICI-arthritis resulted from a limited number of clones that could be found proliferating in the joint. Exposure of synovial T cells to Type I IFN, more so than IFN-γ, induces the CD38hi cytotoxic phenotype. Relative to other CD8+ T cell subsets in the joints, the CD38hi population is distinct from a dysfunctional population and clonally most related to TCF7+ memory populations. Examination of synovial tissue from bilateral knee arthroplasty demonstrated considerable sharing of TCR clonotypes in the CD38hi CD8 T cell fraction from both knees. These results define a distinct CD8 T cell subset that may be directly activated by ICI therapy and mediate a tissue-specific autoimmune cellular reaction in patient joints.

SummaryPro-inflammatory fibroblasts are critical to pathogenesis in rheumatoid arthritis, inflammatory bowel disease, interstitial lung disease, and Sjögren’s syndrome, and represent a novel therapeutic target for chronic inflammatory disease. However, the heterogeneity of fibroblast phenotypes, exacerbated by the lack of a common cross-tissue taxonomy, has limited the understanding of which pathways are shared by multiple diseases. To investigate, we profiled patient-derived fibroblasts from inflamed and non-inflamed synovium, intestine, lung, and salivary glands with single-cell RNA-sequencing. We integrated all fibroblasts into a multi-tissue atlas to characterize shared and tissue-specific phenotypes. Two shared clusters, CXCL10+CCL19+ immune-interacting and SPARC+COL3A1+ vascular-interacting fibroblasts were expanded in all inflamed tissues and additionally mapped to dermal analogues in a public atopic dermatitis atlas. We further confirmed these human pro-inflammatory fibroblasts in animal models of lung, joint, and intestinal inflammation. This work represents the first cross-tissue, single-cell fibroblast atlas revealing shared pathogenic activation states across four chronic inflammatory diseases.

Macrophages regulate protective immune responses to infectious microbes, but aberrant macrophage activation frequently drives pathological inflammation. To identify regulators of vigorous macrophage activation, we analyzed RNA-seq data from synovial macrophages and identified SLAMF7 as a receptor associated with a super-activated macrophage state in rheumatoid arthritis. We implicated IFN-γ as a key regulator of SLAMF7 expression. Engaging this receptor drove an exuberant wave of inflammatory cytokine expression, and induction of TNF-α following SLAMF7 engagement amplified inflammation through an autocrine signaling loop. We observed SLAMF7-induced gene programs not only in macrophages from rheumatoid arthritis patients, but in gut macrophages from active Crohn’s disease patients and lung macrophages from severe COVID-19 patients. This suggests a central role for SLAMF7 in macrophage super-activation with broad implications in pathology.

Current prevailing knowledge on EMT (epithelial mesenchymal transition) deems epithelial cells acquire the characters of mesenchymal cells to be capable of invading and metastasizing on their own. One of the signature events of EMT is called “cadherin switch”, e.g. the epithelial E-cadherin switching to the mesenchymal Cadherin-11. Here, we report the critical events after EMT that cancer cells utilize cadherin-11 to hijack the endogenous cadherin-11 positive fibroblasts. Numerous 3-D cell invasion assays with high-content live cell imaging methods reveal that cadherin-11 positive cancer cells adhere to and migrate back and forth dynamically on the cell bodies of fibroblasts. By adhering to fibroblasts for co-invasion through 3-D matrices, cancer cells acquire higher invasion speed and velocity, as well as significantly elevated invasion persistence, which are exclusive characteristics of fibroblast invasion. Silencing cadherin-11 in cancer cells or in fibroblasts, or in both, significantly decouples such physical co-invasion. Additional bioinformatics studies and PDX (patient derived xenograft) studies link such cadherin-11 mediated cancer hijacking fibroblasts to the clinical cancer progression in human such as triple-negative breast cancer patients. Further animal studies confirm cadherin-11 mediates cancer hijacking fibroblasts in vivo and promotes significant solid tumor progression and distant metastasis. Moreover, overexpression of cadherin-11 strikingly protects 4T1-luc cells from implant rejection against firefly luciferase in immunocompetent mice. Overall, our findings report and characterize the critical post-EMT event of cancer hijacking fibroblasts in cancer progression and suggest cadherin-11 can be a therapeutic target for solid tumors with stroma. Our studies hence provide significant updates on the “EMT” theory that EMT cancer cells can hijack fibroblasts to achieve full mesenchymal behaviors in vivo for efficient homing, growth, metastasis and evasion of immune surveillance. Our studies also reveal that cadherin-11 is the key molecule that helps link cancer cells to stromal fibroblasts in the “Seed & Soil” theory.Graphical Abstract

Macrophages tailor their function to the signals found in tissue microenvironments, taking on a wide spectrum of phenotypes. In human tissues, a detailed understanding of macrophage phenotypes is limited. Using single-cell RNA-sequencing, we define 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 shaped by resident fibroblasts and the cytokine TNF. These macrophages promote fibroblast invasiveness in an EGF receptor dependent manner, indicating that inflammatory intercellular crosstalk reshapes both cell types and contributes to fibroblast-mediated joint destruction. In anex vivotissue assay, the HBEGF+inflammatory macrophage is targeted by several anti-inflammatory RA medications, however, COX inhibition redirects it towards a different inflammatory phenotype that is also expected to perpetuate pathology. These data highlight advances in understanding the pathophysiology and drug mechanisms in chronic inflammatory disorders can be achieved by focusing on macrophage phenotypes in the context of complex interactions in human tissues.One Sentence SummaryA newly identified human macrophage phenotype from patients with the autoimmune condition RA is found to promote joint tissue invasiveness and demonstrates variable sensitivities to anti-inflammatory medications used to treat the disease.

The rapidly emerging diversity of single cell RNAseq datasets allows us to characterize the transcriptional behavior of cell types across a wide variety of biological and clinical conditions. With this comprehensive breadth comes a major analytical challenge. The same cell type across tissues, from different donors, or in different disease states, may appear to express different genes. A joint analysis of multiple datasets requires the integration of cells across diverse conditions. This is particularly challenging when datasets are assayed with different technologies in which real biological differences are interspersed with technical differences. We present Harmony, an algorithm that projects cells into a shared embedding in which cells group by cell type rather than dataset-specific conditions. Unlike available single-cell integration methods, Harmony can simultaneously account for multiple experimental and biological factors. We develop objective metrics to evaluate the quality of data integration. In four separate analyses, we demonstrate the superior performance of Harmony to four single-cell-specific integration algorithms. Moreover, we show that Harmony requires dramatically fewer computational resources. It is the only available algorithm that makes the integration of ∼ 106 cells feasible on a personal computer. We demonstrate that Harmony identifies both broad populations and fine-grained subpopulations of PBMCs from datasets with large experimental differences. In a meta-analysis of 14,746 cells from 5 studies of human pancreatic islet cells, Harmony accounts for variation among technologies and donors to successfully align several rare subpopulations. In the resulting integrated embedding, we identify a previously unidentified population of potentially dysfunctional alpha islet cells, enriched for genes active in the Endoplasmic Reticulum (ER) stress response. The abundance of these alpha cells correlates across donors with the proportion of dysfunctional beta cells also enriched in ER stress response genes. Harmony is a fast and flexible general purpose integration algorithm that enables the identification of shared fine-grained subpopulations across a variety of experimental and biological conditions.

SUMMARYThe identification of lymphocyte subsets with non-overlapping effector functions has been pivotal to the development of targeted therapies in immune mediated inflammatory diseases (IMIDs). However it remains unclear whether fibroblast subclasses with non-overlapping functions also exist and are responsible for the wide variety of tissue driven processes observed in IMIDs such as inflammation and damage. Here we identify and describe the biology of distinct subsets of fibroblasts responsible for mediating either inflammation or tissue damage in arthritis. We show that deletion of FAPα+ synovial cells suppressed both inflammation and bone erosions in murine models of resolving and persistent arthritis. Single cell transcriptional analysis identified two distinct fibroblast subsets: FAPα+ THY1+ immune effector fibroblasts located in the synovial sub-lining, and FAPα+ THY1- destructive fibroblasts restricted to the synovial lining. When adoptively transferred into the joint, FAP α+ THY1- fibroblasts selectively mediate bone and cartilage damage with little effect on inflammation whereas transfer of FAP α+ THY1+ fibroblasts resulted in a more severe and persistent inflammatory arthritis, with minimal effect on bone and cartilage. Our findings describing anatomically discrete, functionally distinct fibroblast subsets with non-overlapping functions have important implications for cell based therapies aimed at modulating inflammation and tissue damage.

Lupus nephritis is a potentially fatal autoimmune disease, whose current treatment is ineffective and often toxic. To gain insights into disease mechanisms, we analyzed kidney samples from lupus nephritis patients and healthy controls using single-cell RNA-seq. Our analysis revealed 21 subsets of leukocytes active in disease, including multiple populations of myeloid, T, NK and B cells, demonstrating both pro-inflammatory and resolving responses. We found evidence of local activation of B cells correlated with an age-associated B cell signature, and 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, pointing to potential therapeutic targets. Gene expression of immune cells in urine and kidney was highly correlated, suggesting urine may be a surrogate for kidney biopsies. Our results provide a first comprehensive view of the complex network of leukocytes active in lupus nephritis kidneys.

To define the cell populations in rheumatoid arthritis (RA) driving joint inflammation, we applied single-cell RNA-seq (scRNA-seq), mass cytometry, bulk RNA-seq, and flow cytometry to sorted T cells, B cells, monocytes, and fibroblasts from 51 synovial tissue RA and osteoarthritis (OA) patient samples. Utilizing an integrated computational strategy based on canonical correlation analysis to 5,452 scRNA-seq profiles, we identified 18 unique cell populations. Combining mass cytometry and transcriptomics together revealed cell states expanded in RA synovia: THY1+HLAhigh sublining fibroblasts (OR=33.8), IL1B+ pro-inflammatory monocytes (OR=7.8), CD11c+T-bet+ autoimmune-associated B cells (OR=5.7), and PD-1+Tph/Tfh (OR=3.0). We also defined CD8+ T cell subsets characterized by GZMK+, GZMB+, and GNLY+ expression. Using bulk and single-cell data, we mapped inflammatory mediators to source cell populations, for example attributing IL6 production to THY1+HLAhigh fibroblasts and naïve B cells, and IL1B to pro-inflammatory monocytes. These populations are potentially key mediators of RA pathogenesis.

Innate T cells (ITCs), including invariant natural killer T (iNKT) cells, mucosal-associated invariant T (MAIT) cells, and γδ T cell populations, use conserved antigen receptors generated by somatic recombination to respond to non-peptide antigens in an innate-like manner. Understanding where these cells fit in the scheme of immunity has been a puzzle since their discovery. Here, immunophenotyping of 101 individuals revealed that these populations account for as much as 25% of peripheral human T cells. To better understand these cells, we generated detailed gene expression profiles using low-input RNA-seq and confirmed key findings through protein-level and functional validation. Unbiased transcriptomic analyses revealed a continuous ‘innateness gradient’ with adaptive T cells at one end followed by MAIT, iNKT, Vδ1+T, Vδ2+T, and natural killer cells at the other end. Innateness was characterized by decreased expression of translational machinery genes and reduced proliferative potential, which allowed for prioritization of effector functions, including rapid cytokine and chemokine production, and cytotoxicity. Thus, global transcriptional programs uncovered rapid proliferation and rapid effector functions as competing goals that respectively define adaptive and innate states.One Sentence SummaryAdaptive and innate T cells align along a continuous innateness gradient, reflecting a trade-off between effector function and proliferative capacity.

OBJECTIVEThere is a critical need to define the cells that mediate tissue damage in lupus nephritis. Here we aimed to establish a protocol to preserve lupus nephritis kidney biopsies and urine cell samples obtained at multiple clinical sites for subsequent isolation and transcriptomic analysis of single cells.METHODSFresh and cryopreserved kidney tissue from tumor nephrectomies and lupus nephritis kidney biopsies were disaggregated by enzymatic digestion. Cell yields and cell composition were assessed by flow cytometry. Transcriptomes of leukocytes and epithelial cells were evaluated by low-input and single cell RNA-seq.RESULTSCryopreserved kidney tissue from tumor nephrectomies and lupus nephritis biopsies can be thawed and dissociated to yield intact, viable leukocytes and epithelial cells. Cryopreservation of intact kidney tissue provides higher epithelial cell yields compared to cryopreservation of single cell suspensions from dissociated kidneys. Cell yields and flow cytometric cell phenotypes are comparable between cryopreserved kidney samples and paired kidney samples shipped overnight on wet ice. High-quality single cell and low-input transcriptomic data were generated from leukocytes from both cryopreserved lupus nephritis kidney biopsies and urine, as well as from a subset of kidney epithelial cells.CONCLUSIONThe AMP RA/SLE cryopreserved tissue analysis pipeline provides a method for centralized processing of lupus nephritis kidney biopsies and urine samples to generate robust transcriptomic analyses in multi-center studies.

High dimensional single-cell analyses have dramatically 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 due to technical and inter-individual variation. Here we present Mixed effects modeling of Associations of Single Cells (MASC), a novel 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 confounds and biological variation. Applying MASC to mass cytometry analyses of CD4+ T cells from 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 (OR = 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 respond to immunosuppressive therapy. Compared to peripheral blood, synovial fluid and synovial tissue samples from RA patients contained ∼5-fold higher frequencies of CD27- HLA-DR+ cells, which comprised ∼10% of synovial CD4+ T cells. We find that CD27- HLA-DR+ cells are abundant producers of IFN-γ and also express perforin and granzyme A at elevated levels. Thus MASC identified the expansion of a unique Th1 skewed effector T cell population with cytotoxic capacity in RA. We propose that MASC is a broadly applicable method to identify disease-associated cell populations in high-dimensional single cell data.One Sentence SummaryMixed-effects regression of single cells identifies a cytotoxic Th1-like CD4+ T cell subset while accounting for inter-individual and technical variation.

Fibroblasts mediate normal tissue matrix remodeling, but they can cause fibrosis or tissue destruction following chronic inflammation. In rheumatoid arthritis (RA), synovial fibroblasts expand, degrade cartilage, and drive joint inflammation. Little is known about fibroblast heterogeneity or if aberrations in fibroblast subsets relate to disease pathology. Here, we used an integrative strategy, including bulk transcriptomics on targeted subpopulations and unbiased single-cell transcriptomics, to analyze fibroblasts from synovial tissues. We identify 7 phenotypic fibroblast subsets with distinct surface protein phenotypes, and these collapsed into 3 subsets based on transcriptomics data. One subset expressing PDPN, THY1, but lacking CD34 was 3-fold expanded in RA relative to osteoarthritis (P=0.007); most of these cells expressed CDH11. The subsets were found to differ in expression of cytokines and matrix metalloproteinases, localization in synovial microanatomy, and in response to TNF. Our approach provides a template to identify pathogenic stromal cellular subsets in complex diseases.