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.
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.
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-α.
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.
Full text
Available for:
IZUM, KILJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK
T cell-derived pro-inflammatory cytokines are a major driver of rheumatoid arthritis (RA) pathogenesis. Although these cytokines have traditionally been attributed to CD4 T cells, we have found that ...CD8 T cells are notably abundant in synovium and make more interferon (IFN)-γ and nearly as much tumor necrosis factor (TNF) as their CD4 T cell counterparts. Furthermore, using unbiased high-dimensional single-cell RNA-seq and flow cytometric data, we found that the vast majority of synovial tissue and synovial fluid CD8 T cells belong to an effector CD8 T cell population characterized by high expression of granzyme K (GzmK) and low expression of granzyme B (GzmB) and perforin. Functional experiments demonstrate that these GzmK
GzmB
CD8 T cells are major cytokine producers with low cytotoxic potential. Using T cell receptor repertoire data, we found that CD8 GzmK
GzmB
T cells are clonally expanded in synovial tissues and maintain their granzyme expression and overall cell state in blood, suggesting that they are enriched in tissue but also circulate. Using GzmK and GzmB signatures, we found that GzmK-expressing CD8 T cells were also the major CD8 T cell population in the gut, kidney, and coronavirus disease 2019 (COVID-19) bronchoalveolar lavage fluid, suggesting that they form a core population of tissue-associated T cells across diseases and human tissues. We term this population tissue-enriched expressing GzmK or T
CD8 cells. Armed to produce cytokines in response to both antigen-dependent and antigen-independent stimuli, CD8 T
cells have the potential to drive inflammation.
Abstract
Classical dendritic cells (cDC) are professional antigen-presenting cells (APC) that regulate immunity and tolerance. Neutrophil-derived cells with properties of DCs (nAPC) are observed in ...human diseases and after culture of neutrophils with cytokines. Here we show that FcγR-mediated endocytosis of antibody-antigen complexes or an anti-FcγRIIIB-antigen conjugate converts neutrophils into nAPCs that, in contrast to those generated with cytokines alone, activate T cells to levels observed with cDCs and elicit CD8
+
T cell-dependent anti-tumor immunity in mice. Single cell transcript analyses and validation studies implicate the transcription factor PU.1 in neutrophil to nAPC conversion. In humans, blood nAPC frequency in lupus patients correlates with disease. Moreover, anti-FcγRIIIB-antigen conjugate treatment induces nAPCs that can activate autologous T cells when using neutrophils from individuals with myeloid neoplasms that harbor neoantigens or those vaccinated against bacterial toxins. Thus, anti-FcγRIIIB-antigen conjugate-induced conversion of neutrophils to immunogenic nAPCs may represent a possible immunotherapy for cancer and infectious diseases.
Psoriasis is a chronic, systemic inflammatory condition primarily affecting skin. While the role of the immune compartment (e.g., T cells) is well established, the changes in the skin compartment are ...more poorly understood. Using longitudinal skin biopsies (n = 375) from the "Psoriasis Treatment with Abatacept and Ustekinumab: A Study of Efficacy"(PAUSE) clinical trial (n = 101), we report 953 expression quantitative trait loci (eQTLs). Of those, 116 eQTLs have effect sizes that were modulated by local skin inflammation (eQTL interactions). By examining these eQTL genes (eGenes), we find that most are expressed in the skin tissue compartment, and a subset overlap with the NRF2 pathway. Indeed, the strongest eQTL interaction signal - rs1491377616-LCE3C - links a psoriasis risk locus with a gene specifically expressed in the epidermis. This eQTL study highlights the potential to use biospecimens from clinical trials to discover in vivo eQTL interactions with therapeutically relevant environmental variables.
Rheumatoid arthritis is a prototypical autoimmune disease that causes joint inflammation and destruction
. There is currently no cure for rheumatoid arthritis, and the effectiveness of treatments ...varies across patients, suggesting an undefined pathogenic diversity
. 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.
Full text
Available for:
GEOZS, IJS, IMTLJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBMB, UL, UM, UPUK, ZAGLJ
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.
Full text
Available for:
GEOZS, IJS, IMTLJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBMB, UL, UM, UPUK, ZAGLJ
Objective
Recent studies have uncovered diverse cell types and states in the rheumatoid arthritis (RA) synovium; however, limited data exist correlating these findings with patient‐level clinical ...information. Using the largest cohort to date with clinical and multicell data, we determined associations between RA clinical factors with cell types and states in the RA synovium.
Methods
The Accelerated Medicines Partnership Rheumatoid Arthritis study recruited patients with active RA who were not receiving disease‐modifying antirheumatic drugs (DMARDs) or who had an inadequate response to methotrexate (MTX) or tumor necrosis factor inhibitors. RA clinical factors were systematically collected. Biopsies were performed on an inflamed joint, and tissue were disaggregated and processed with a cellular indexing of transcriptomes and epitopes sequencing pipeline from which the following cell type percentages and cell type abundance phenotypes (CTAPs) were derived: endothelial, fibroblast, and myeloid (EFM); fibroblasts; myeloid; T and B cells; T cells and fibroblasts (TF); and T and myeloid cells. Correlations were measured between RA clinical factors, cell type percentage, and CTAPs.
Results
We studied 72 patients (mean age 57 years, 75% women, 83% seropositive, mean RA duration 6.6 years, mean Disease Activity Score‐28 C‐reactive Protein 3 DAS28‐CRP3 score 4.8). Higher DAS28‐CRP3 correlated with a higher T cell percentage (P < 0.01). Those receiving MTX and not a biologic DMARD (bDMARD) had a higher percentage of B cells versus those receiving no DMARDs (P < 0.01). Most of those receiving bDMARDs were categorized as EFM (57%), whereas none were TF. No significant difference was observed across CTAPs for age, sex, RA disease duration, or DAS28‐CRP3.
Conclusion
In this comprehensive screen of clinical factors, we observed differential associations between DMARDs and cell phenotypes, suggesting that RA therapies, more than other clinical factors, may impact cell type/state in the synovium and ultimately influence response to subsequent therapies.
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SAZU, SBCE, SBMB, UL, UM, UPUK
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.
Full text
Available for:
GEOZS, IJS, IMTLJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBMB, UL, UM, UPUK, ZAGLJ
Delayed detection of LN associates with worse outcomes. There are conflicting recommendations regarding a threshold level of proteinuria at which biopsy will likely yield actionable management. This ...study addressed the association of urine protein:creatinine ratios (UPCR) with clinical characteristics and investigated the incidence of proliferative and membranous histology in patients with a UPCR between 0.5 and 1.
A total of 275 SLE patients (113 first biopsy, 162 repeat) were enrolled in the multicentre multi-ethnic/racial Accelerating Medicines Partnership across 15 US sites at the time of a clinically indicated renal biopsy. Patients were followed for 1 year.
At biopsy, 54 patients had UPCR <1 and 221 had UPCR ≥1. Independent of UPCR or biopsy number, a majority (92%) of patients had class III, IV, V or mixed histology. Moreover, patients with UPCR <1 and class III, IV, V, or mixed had a median activity index of 4.5 and chronicity index of 3, yet 39% of these patients had an inactive sediment. Neither anti-dsDNA nor low complement distinguished class I or II from III, IV, V or mixed in patients with UPCR <1. Of 29 patients with baseline UPCR <1 and class III, IV, V or mixed, 23 (79%) had a UPCR <0.5 at 1 year.
In this prospective study, three-quarters of patients with UPCR <1 had histology showing class III, IV, V or mixed with accompanying activity and chronicity despite an inactive sediment or normal serologies. These data support renal biopsy at thresholds lower than a UPCR of 1.
Rheumatoid arthritis (RA) is an autoimmune disease involving antigen-specific T and B cells. Here, we perform single-cell RNA and repertoire sequencing on paired synovial tissue and blood samples ...from 12 seropositive RA patients. We identify clonally expanded CD4 + T cells, including CCL5+ cells and T peripheral helper (Tph) cells, which show a prominent transcriptomic signature of recent activation and effector function. CD8 + T cells show higher oligoclonality than CD4 + T cells, with the largest synovial clones enriched in GZMK+ cells. CD8 + T cells with possibly virus-reactive TCRs are distributed across transcriptomic clusters. In the B cell compartment, NR4A1+ activated B cells, and plasma cells are enriched in the synovium and demonstrate substantial clonal expansion. We identify synovial plasma cells that share BCRs with synovial ABC, memory, and activated B cells. Receptor-ligand analysis predicted IFNG and TNFRSF members as mediators of synovial Tph-B cell interactions. Together, these results reveal clonal relationships between functionally distinct lymphocyte populations that infiltrate the synovium of patients with RA.Activated B cells and T cells accumulate within joints of patients with rheumatoid arthritis. Here, the authors use single-cell transcriptome and repertoire profiling to identify clonally expanded synovial B cells and T cells and define their phenotypes and predicted cell-cell interactions.