SLIDE: Significant Latent Factor Interaction Discovery and Exploration across biological domains.

Modern multiomic technologies can generate deep multiscale profiles. However, differences in data modalities, multicollinearity of the data, and large numbers of irrelevant features make analyses and integration of high-dimensional omic datasets challenging. Here we present Significant Latent Factor Interaction Discovery and Exploration (SLIDE), a first-in-class interpretable machine learning technique for identifying significant interacting latent factors underlying outcomes of interest from high-dimensional omic datasets. SLIDE makes no assumptions regarding data-generating mechanisms, comes with theoretical guarantees regarding identifiability of the latent factors/corresponding inference, and has rigorous false discovery rate control. Using SLIDE on single-cell and spatial omic datasets, we uncovered significant interacting latent factors underlying a range of molecular, cellular and organismal phenotypes. SLIDE outperforms/performs at least as well as a wide range of state-of-the-art approaches, including other latent factor approaches. More importantly, it provides biological inference beyond prediction that other methods do not afford. Thus, SLIDE is a versatile engine for biological discovery from modern multiomic datasets.

Increased CD8+ tissue resident memory T cells, regulatory T cells and activated natural killer cells in systemic sclerosis lungs.

OBJECTIVE: Multiple observations indicate a role for lymphocytes in driving autoimmunity in SSc. While T and NK cells have been studied in SSc whole blood and bronchoalveolar lavage fluid, their role remains unclear, partly because no studies have analysed these cell types in SSc-interstitial lung disease (ILD) lung tissue. This research aimed to identify and analyse the lymphoid subpopulations in SSc-ILD lung explants. METHODS: Lymphoid populations from 13 SSc-ILD and 6 healthy control (HC) lung explants were analysed using Seurat following single-cell RNA sequencing. Lymphoid clusters were identified by their differential gene expression. Absolute cell numbers and cell proportions in each cluster were compared between cohorts. Additional analyses were performed using pathway analysis, pseudotime and cell ligand-receptor interactions. RESULTS: Activated CD16+ NK cells, CD8+ tissue resident memory T cells and Treg cells were proportionately higher in SSc-ILD compared with HC lungs. Activated CD16+ NK cells in SSc-ILD showed upregulated granzyme B, IFN-γ and CD226. Amphiregulin, highly upregulated by NK cells, was predicted to interact with epidermal growth factor receptor on several bronchial epithelial cell populations. Shifts in CD8+ T cell populations indicated a transition from resting to effector to tissue resident phenotypes in SSc-ILD. CONCLUSIONS: SSc-ILD lungs show activated lymphoid populations. Activated cytotoxic NK cells suggest they may kill alveolar epithelial cells, while their expression of amphiregulin suggests they may also induce bronchial epithelial cell hyperplasia. CD8+ T cells in SSc-ILD appear to transition from resting to the tissue resident memory phenotype.

Increased expression of CXCL6 in secretory cells drives fibroblast collagen synthesis and is associated with increased mortality in idiopathic pulmonary fibrosis.

RATIONALE: Recent data suggest that the localisation of airway epithelial cells in the distal lung in idiopathic pulmonary fibrosis (IPF) may drive pathology. We set out to discover whether chemokines expressed in these ectopic airway epithelial cells may contribute to the pathogenesis of IPF. METHODS: We analysed whole lung and single-cell transcriptomic data obtained from patients with IPF. In addition, we measured chemokine levels in blood, bronchoalveolar lavage (BAL) of IPF patients and air-liquid interface cultures. We employed ex vivo donor and IPF lung fibroblasts and an animal model of pulmonary fibrosis to test the effects of chemokine signalling on fibroblast function. RESULTS: By analysis of whole-lung transcriptomics, protein and BAL, we discovered that CXCL6 (a member of the interleukin-8 family) was increased in patients with IPF. Elevated CXCL6 levels in the BAL of two cohorts of patients with IPF were associated with poor survival (hazard ratio of death or progression 1.89, 95% CI 1.16-3.08; n=179, p=0.01). By immunostaining and single-cell RNA sequencing, CXCL6 was detected in secretory cells. Administration of mCXCL5 (LIX, murine CXCL6 homologue) to mice increased collagen synthesis with and without bleomycin. CXCL6 increased collagen I levels in donor and IPF fibroblasts 4.4-fold and 1.7-fold, respectively. Both silencing of and chemical inhibition of CXCR1/2 blocked the effects of CXCL6 on collagen, while overexpression of CXCR2 increased collagen I levels 4.5-fold in IPF fibroblasts. CONCLUSIONS: CXCL6 is expressed in ectopic airway epithelial cells. Elevated levels of CXCL6 are associated with IPF mortality. CXCL6-driven collagen synthesis represents a functional consequence of ectopic localisation of airway epithelial cells in IPF.

Fibroblast Subpopulations in Systemic Sclerosis: Functional Implications of Individual Subpopulations and Correlations with Clinical Features.

Fibroblasts constitute a heterogeneous population of cells. In this study, we integrated single-cell RNA-sequencing and bulk RNA-sequencing data as well as clinical information to study the role of individual fibroblast populations in systemic sclerosis (SSc). SSc skin demonstrated an increased abundance of COMP+, COL11A1+, MYOC+, CCL19+, SFRP4/SFRP2+, and PRSS23/SFRP2+ fibroblasts signatures and decreased proportions of CXCL12+ and PI16+ fibroblast signatures in the Prospective Registry of Early Systemic Sclerosis and Genetics versus Environment in Scleroderma Outcome Study cohorts. Numerical differences were confirmed by multicolor immunofluorescence for selected fibroblast populations. COMP+, COL11A1+, SFRP4/SFRP2+, PRSS23/SFRP2+, and PI16+ fibroblasts were similarly altered between normal wound healing and patients with SSc. The proportions of profibrotic COMP+, COL11A1+, SFRP4/SFRP2+, and PRSS23/SFRP2+ and proinflammatory CCL19+ fibroblast signatures were positively correlated with clinical and histopathological parameters of skin fibrosis, whereas signatures of CXCL12+ and PI16+ fibroblasts were inversely correlated. Incorporating the proportions of COMP+, COL11A1+, SFRP4/SFRP2+, and PRSS23/SFRP2+ fibroblast signatures into machine learning models improved the classification of patients with SSc into those with progressive versus stable skin fibrosis. In summary, the profound imbalance of fibroblast subpopulations in SSc may drive the progression of skin fibrosis. Specific targeting of disease-relevant fibroblast populations may offer opportunities for the treatment of SSc and other fibrotic diseases.

Antigen-driven T cell-macrophage interactions mediate the interface between innate and adaptive immunity in histidyl-tRNA synthetase-induced myositis.

Introduction: Previous work in humans has demonstrated that both innate and adaptive immune signaling pathways contribute to the pathogenesis of idiopathic inflammatory myopathy (IIM), a systemic autoimmune disease targeting muscle as well as extra-muscular organs. To better define interactive signaling networks in IIM, we characterized the cellular phenotype and transcriptomic profiles of muscle-infiltrating cells in our established murine model of histidyl-tRNA synthetase (HRS)-induced myositis. Methods: Myositis was induced in wild type (WT) and various congenic/mutant strains of C57BL/6 mice through intramuscular immunization with recombinant HRS. Histopathological, immunohistochemical, flow cytometric, and transcriptomic assessments were used to characterize the functional relationship between muscle-infiltrating cell populations in these strains lacking different components of innate and/or adaptive immune signaling. Results: RAG1 KO mice developed markedly reduced muscle inflammation relative to WT mice, demonstrating a key requirement for T cells in driving HRS-induced myositis. While the reduction of mononuclear cell infiltrates in CD4-Cre.MyD88fl/fl conditional knockout mice and OT-II TCR transgenic mice highlighted roles for both innate and TCR-mediated/adaptive immune signaling in T cells, diminished inflammation in Lyz2-Cre.MyD88fl/fl conditional knockout mice underscored the importance of macrophage/myeloid cell populations in supporting T cell infiltration. Single cell RNA sequencing-based clustering of muscle-infiltrating subpopulations and associated pathway analyses showed that perturbations of T cell signaling/function alter the distribution and phenotype of macrophages, fibroblasts, and other non-lymphoid cell populations contributing to HRS-induced myositis. Discussion: Overall, HRS-induced myositis reflects the complex interplay between multiple cell types that collectively drive a TH1-predominant, pro-inflammatory tissue phenotype requiring antigen-mediated activation of both MyD88- and TCR-dependent T cell signaling pathways.

The mycosis fungoides cutaneous microenvironment shapes dysfunctional cell trafficking, antitumor immunity, matrix interactions, and angiogenesis.

Malignant T lymphocyte proliferation in mycosis fungoides (MF) is largely restricted to the skin, implying that malignant cells are dependent on their specific cutaneous tumor microenvironment (TME), including interactions with non-malignant immune and stromal cells, cytokines, and other immunomodulatory factors. To explore these interactions, we performed a comprehensive transcriptome analysis of the TME in advanced-stage MF skin tumors by single-cell RNA sequencing. Our analysis identified cell-type compositions, cellular functions, and cell-to-cell interactions in the MF TME that were distinct from those from healthy skin and benign dermatoses. While patterns of gene expression were common among patient samples, high transcriptional diversity was also observed in immune and stromal cells, with dynamic interactions and crosstalk between these cells and malignant T lymphocytes. This heterogeneity mapped to processes such as cell trafficking, matrix interactions, angiogenesis, immune functions, and metabolism that affect cancer cell growth, migration, and invasion, as well as antitumor immunity. By comprehensively characterizing the transcriptomes of immune and stromal cells within the cutaneous microenvironment of individual MF tumors, we have identified patterns of dysfunction common to all tumors that represent a resource for identifying candidates with therapeutic potential as well as patient-specific heterogeneity that has important implications for personalized disease management.

Transcriptional changes of the aging lung.

Aging is a natural process associated with declined organ function and higher susceptibility to developing chronic diseases. A systemic single-cell type-based study provides a unique opportunity to understand the mechanisms behind age-related pathologies. Here, we use single-cell gene expression analysis comparing healthy young and aged human lungs from nonsmoker donors to investigate age-related transcriptional changes. Our data suggest that aging has a heterogenous effect on lung cells, as some populations are more transcriptionally dynamic while others remain stable in aged individuals. We found that monocytes and alveolar macrophages were the most transcriptionally affected populations. These changes were related to inflammation and regulation of the immune response. Additionally, we calculated the LungAge score, which reveals the diversity of lung cell types during aging. Changes in DNA damage repair, fatty acid metabolism, and inflammation are essential for age prediction. Finally, we quantified the senescence score in aged lungs and found that the more biased cells toward senescence are immune and progenitor cells. Our study provides a comprehensive and systemic analysis of the molecular signatures of lung aging. Our LungAge signature can be used to predict molecular signatures of physiological aging and to detect common signatures of age-related lung diseases.

End-Stage Idiopathic Pulmonary Fibrosis Lung Microenvironment Promotes Impaired NK Activity.

Idiopathic pulmonary fibrosis (IPF) is a fibrotic age-related chronic lung disease characterized by the accumulation of senescent cells. Whether impaired immune response is responsible for the accumulation of senescent cells in the IPF lung remains unknown. In this study, we characterized the NK phenotype in IPF lungs via flow cytometry using 5-dodecanoylaminofluorescein di-ß-d-galactopyranoside, markers of tissue residence, and chemokine receptors. The effect of the lung microenvironment was evaluated using lung fibroblast (LF) conditioned media (CM), and the bleomycin-induced pulmonary fibrosis mouse model was used to assess the in vivo relationship between NK cells and the accumulation of senescent cells. We found that NK cells from the lower lobe of IPF patients exhibited immune-senescent and impaired CD57-NKG2A+ phenotype. We also observed that culture of NK cells from healthy donors in CM from IPF lower lobe lung fibroblasts induced a senescent-like phenotype and impaired cytotoxic capacity. There is an impaired NK recruitment by LF, and NKs presented decreased migration toward their CM. In addition, NK cell-depleted mice treated with bleomycin showed increased collagen deposition and accumulation of different populations of senescent cells compared with controls. The IPF lung microenvironment induces a dysfunctional NK phenotype limiting the clearance of lung senescent cells and the resolution of lung fibrosis. We propose that impaired NK activity could be one of the mechanisms responsible for perpetuating the accumulation of senescent cells in IPF lungs.

An integrated cell atlas of the lung in health and disease.

Single-cell technologies have transformed our understanding of human tissues. Yet, studies typically capture only a limited number of donors and disagree on cell type definitions. Integrating many single-cell datasets can address these limitations of individual studies and capture the variability present in the population. Here we present the integrated Human Lung Cell Atlas (HLCA), combining 49 datasets of the human respiratory system into a single atlas spanning over 2.4 million cells from 486 individuals. The HLCA presents a consensus cell type re-annotation with matching marker genes, including annotations of rare and previously undescribed cell types. Leveraging the number and diversity of individuals in the HLCA, we identify gene modules that are associated with demographic covariates such as age, sex and body mass index, as well as gene modules changing expression along the proximal-to-distal axis of the bronchial tree. Mapping new data to the HLCA enables rapid data annotation and interpretation. Using the HLCA as a reference for the study of disease, we identify shared cell states across multiple lung diseases, including SPP1+ profibrotic monocyte-derived macrophages in COVID-19, pulmonary fibrosis and lung carcinoma. Overall, the HLCA serves as an example for the development and use of large-scale, cross-dataset organ atlases within the Human Cell Atlas.

Single-Cell Transcriptome Analysis Identifies Subclusters with Inflammatory Fibroblast Responses in Localized Scleroderma.

Localized scleroderma (LS) is an autoimmune disease with both inflammatory and fibrotic components causing an abnormal deposition of collagen in the skin and underlying tissue, often leading to disfigurement and disability. Much of its pathophysiology is extrapolated from systemic sclerosis (SSc) since the histopathology findings in the skin are nearly identical. However, LS is critically understudied. Single-cell RNA sequencing (scRNA seq) technology provides a novel way to obtain detailed information at the individual cellular level, overcoming this barrier. Here, we analyzed the affected skin of 14 patients with LS (pediatric and adult) and 14 healthy controls. Fibroblast populations were the focus, since they are the main drivers of fibrosis in SSc. We identified 12 fibroblast subclusters in LS, which overall had an inflammatory gene expression (IFN and HLA-associated genes). A myofibroblast-like cluster (SFRP4/PRSS23) was more prevalent in LS subjects and shared many upregulated genes expressed in SSc-associated myofibroblasts, though it also had strong expression of CXCL9/10/11, known CXCR3 ligands. A CXCL2/IRF1 cluster identified was unique to LS, with a robust inflammatory gene signature, including IL-6, and according to cell communication analysis are influenced by macrophages. In summary, potential disease-propagating fibroblasts and associated gene signatures were identified in LS skin via scRNA seq.

Single-nucleus chromatin accessibility identifies a critical role for TWIST1 in idiopathic pulmonary fibrosis myofibroblast activity.

BACKGROUND: In idiopathic pulmonary fibrosis (IPF), myofibroblasts are key effectors of fibrosis and architectural distortion by excessive deposition of extracellular matrix and their acquired contractile capacity. Single-cell RNA-sequencing (scRNA-seq) has precisely defined the IPF myofibroblast transcriptome, but identifying critical transcription factor activity by this approach is imprecise. METHODS: We performed single-nucleus assay for transposase-accessible chromatin sequencing on explanted lungs from patients with IPF (n=3) and donor controls (n=2) and integrated this with a larger scRNA-seq dataset (10 IPF, eight controls) to identify differentially accessible chromatin regions and enriched transcription factor motifs within lung cell populations. We performed RNA-sequencing on pulmonary fibroblasts of bleomycin-injured Twist1-overexpressing COL1A2 Cre-ER mice to examine alterations in fibrosis-relevant pathways following Twist1 overexpression in collagen-producing cells. RESULTS: TWIST1, and other E-box transcription factor motifs, were significantly enriched in open chromatin of IPF myofibroblasts compared to both IPF nonmyogenic (log2 fold change (FC) 8.909, adjusted p-value 1.82×10-35) and control fibroblasts (log2FC 8.975, adjusted p-value 3.72×10-28). TWIST1 expression was selectively upregulated in IPF myofibroblasts (log2FC 3.136, adjusted p-value 1.41×10- 24), with two regions of TWIST1 having significantly increased accessibility in IPF myofibroblasts. Overexpression of Twist1 in COL1A2-expressing fibroblasts of bleomycin-injured mice resulted in increased collagen synthesis and upregulation of genes with enriched chromatin accessibility in IPF myofibroblasts. CONCLUSIONS: Our studies utilising human multiomic single-cell analyses combined with in vivo murine disease models confirm a critical regulatory function for TWIST1 in IPF myofibroblast activity in the fibrotic lung. Understanding the global process of opening TWIST1 and other E-box transcription factor motifs that govern myofibroblast differentiation may identify new therapeutic interventions for fibrotic pulmonary diseases.

Cell Type-Specific Biomarkers of Systemic Sclerosis Disease Severity Capture Cell-Intrinsic and Cell-Extrinsic Circuits.

OBJECTIVE: Systemic sclerosis (SSc) is a multifactorial autoimmune fibrotic disorder involving complex rewiring of cell-intrinsic and cell-extrinsic signaling coexpression networks involving a range of cell types. However, the rewired circuits as well as corresponding cell-cell interactions remain poorly understood. To address this, we used a predictive machine learning framework to analyze single-cell RNA-sequencing data from 24 SSc patients across the severity spectrum as quantified by the modified Rodnan skin score (MRSS). METHODS: We used a least absolute shrinkage and selection operator (LASSO)-based predictive machine learning approach on the single-cell RNA-sequencing data set to identify predictive biomarkers of SSc severity, both across and within cell types. The use of L1 regularization helps prevent overfitting on high-dimensional data. Correlation network analyses were coupled to the LASSO model to identify cell-intrinsic and cell-extrinsic co-correlates of the identified biomarkers of SSc severity. RESULTS: We found that the uncovered cell type-specific predictive biomarkers of MRSS included previously implicated genes in fibroblast and myeloid cell subsets (e.g., SFPR2+ fibroblasts and monocytes), as well as novel gene biomarkers of MRSS, especially in keratinocytes. Correlation network analyses revealed novel cross-talk between immune pathways and implicated keratinocytes in addition to fibroblast and myeloid cells as key cell types involved in SSc pathogenesis. We then validated the uncovered association of key gene expression and protein markers in keratinocytes, KRT6A and S100A8, with SSc skin disease severity. CONCLUSION: Our global systems analyses reveal previously uncharacterized cell-intrinsic and cell-extrinsic signaling coexpression networks underlying SSc severity that involve keratinocytes, myeloid cells, and fibroblasts.

Early events marking lung fibroblast transition to profibrotic state in idiopathic pulmonary fibrosis.

BACKGROUND: Idiopathic Pulmonary Fibrosis (IPF) is an age-associated progressive lung disease with accumulation of scar tissue impairing gas exchange. Previous high-throughput studies elucidated the role of cellular heterogeneity and molecular pathways in advanced disease. However, critical pathogenic pathways occurring in the transition of fibroblasts from normal to profibrotic have been largely overlooked. METHODS: We used single cell transcriptomics (scRNA-seq) from lungs of healthy controls and IPF patients (lower and upper lobes). We identified fibroblast subclusters, genes and pathways associated with early disease. Immunofluorescence assays validated the role of MOXD1 early in fibrosis. RESULTS: We identified four distinct fibroblast subgroups, including one marking the normal-to-profibrotic state transition. Our results show for the first time that global downregulation of ribosomal proteins and significant upregulation of the majority of copper-binding proteins, including MOXD1, mark the IPF transition. We find no significant differences in gene expression in IPF upper and lower lobe samples, which were selected to have low and high degree of fibrosis, respectively. CONCLUSIONS: Early events during IPF onset in fibroblasts include dysregulation of ribosomal and copper-binding proteins. Fibroblasts in early stage IPF may have already acquired a profibrotic phenotype while hallmarks of advanced disease, including fibroblast foci and honeycomb formation, are still not evident. The new transitional fibroblasts we discover could prove very important for studying the role of fibroblast plasticity in disease progression and help develop early diagnosis tools and therapeutic interventions targeting earlier disease states.

An individualized causal framework for learning intercellular communication networks that define microenvironments of individual tumors.

Cells within a tumor microenvironment (TME) dynamically communicate and influence each other's cellular states through an intercellular communication network (ICN). In cancers, intercellular communications underlie immune evasion mechanisms of individual tumors. We developed an individualized causal analysis framework for discovering tumor specific ICNs. Using head and neck squamous cell carcinoma (HNSCC) tumors as a testbed, we first mined single-cell RNA-sequencing data to discover gene expression modules (GEMs) that reflect the states of transcriptomic processes within tumor and stromal single cells. By deconvoluting bulk transcriptomes of HNSCC tumors profiled by The Cancer Genome Atlas (TCGA), we estimated the activation states of these transcriptomic processes in individual tumors. Finally, we applied individualized causal network learning to discover an ICN within each tumor. Our results show that cellular states of cells in TMEs are coordinated through ICNs that enable multi-way communications among epithelial, fibroblast, endothelial, and immune cells. Further analyses of individual ICNs revealed structural patterns that were shared across subsets of tumors, leading to the discovery of 4 different subtypes of networks that underlie disparate TMEs of HNSCC. Patients with distinct TMEs exhibited significantly different clinical outcomes. Our results show that the capability of estimating individual ICNs reveals heterogeneity of ICNs and sheds light on the importance of intercellular communication in impacting disease development and progression.

Tofacitinib blocks IFN-regulated biomarker genes in skin fibroblasts and keratinocytes in a systemic sclerosis trial.

BACKGROUNDSystemic sclerosis (SSc) is an autoimmune, connective tissue disease characterized by vasculopathy and fibrosis of the skin and internal organs.METHODSWe randomized 15 participants with early diffuse cutaneous SSc to tofacitinib 5 mg twice a day or matching placebo in a phase I/II double-blind, placebo-controlled trial. The primary outcome measure was safety and tolerability at or before week 24. To understand the changes in gene expression associated with tofacitinib treatment in each skin cell population, we compared single-cell gene expression in punch skin biopsies obtained at baseline and 6 weeks following the initiation of treatment.RESULTSTofacitinib was well tolerated; no participants experienced grade 3 or higher adverse events before or at week 24. Trends in efficacy outcome measures favored tofacitnib. Baseline gene expression in fibroblast and keratinocyte subpopulations indicated IFN-activated gene expression. Tofacitinib inhibited IFN-regulated gene expression in SFRP2/DPP4 fibroblasts (progenitors of myofibroblasts) and in MYOC/CCL19, representing adventitial fibroblasts (P < 0.05), as well as in the basal and keratinized layers of the epidermis. Gene expression in macrophages and DCs indicated inhibition of STAT3 by tofacitinib (P < 0.05). No clinically meaningful inhibition of T cells and endothelial cells in the skin tissue was observed.CONCLUSIONThese results indicate that mesenchymal and epithelial cells of a target organ in SSc, not the infiltrating lymphocytes, may be the primary focus for therapeutic effects of a Janus kinase inhibitor.TRIAL REGISTRATIONClinicalTrials.gov NCT03274076.FUNDINGPfizer, NIH/National Institute of Arthritis and Musculoskeletal and Skin Diseases (NIAMS) R01 AR070470, NIH/NIAMS K24 AR063120, Taubman Medical Research Institute and NIH P30 AR075043, and NIH/NIAMS K01 AR072129.

Epigenetic Regulation of Profibrotic Macrophages in Systemic Sclerosis-Associated Interstitial Lung Disease.

OBJECTIVE: Systemic sclerosis-associated interstitial lung disease (SSc-ILD) is the leading cause of death in patients with SSc with unclear pathogenesis and limited treatment options. Evidence strongly supports an important role for profibrotic secreted phosphoprotein 1 (SPP1)-expressing macrophages in SSc-ILD. This study was undertaken to define the transcriptome and chromatin structural changes of SPP1 SSc-ILD macrophages in order to better understand their role in promoting fibrosis and to identify transcription factors associated with open chromatin driving their altered phenotype. METHODS: We performed single-cell RNA sequencing (scRNA-Seq) on 11 explanted SSc-ILD and healthy control lung samples, as well as single-cell assay for transposase-accessible chromatin sequencing on 5 lung samples to define altered chromatin accessibility of SPP1 macrophages. We predicted transcription factors regulating SPP1 macrophages using single-cell regulatory network inference and clustering (SCENIC) and determined transcription factor binding sites associated with global alterations in SPP1 chromatin accessibility using Signac/Seurat. RESULTS: We identified distinct macrophage subpopulations using scRNA-Seq analysis in healthy and SSc-ILD lungs and assessed gene expression changes during the change of healthy control macrophages into SPP1 macrophages. Analysis of open chromatin validated SCENIC predictions, indicating that microphthalmia-associated transcription factor, transcription factor EB, activating transcription factor 6, sterol regulatory element binding transcription factor 1, basic helix-loop-helix family member E40, Kruppel-like factor 6, ETS variant transcription factor 5, and/or members of the activator protein 1 family of transcription factors regulate SPP1 macrophage differentiation. CONCLUSION: Our findings shed light on the underlying changes in chromatin structure and transcription factor regulation of profibrotic SPP1 macrophages in SSc-ILD. Similar alterations in SPP1 macrophages may underpin fibrosis in other organs involved in SSc and point to novel targets for the treatment of SSc-ILD, specifically targeting profibrotic macrophages.

Autoreactive CD8+ T cells are restrained by an exhaustion-like program that is maintained by LAG3.

Impaired chronic viral and tumor clearance has been attributed to CD8+ T cell exhaustion, a differentiation state in which T cells have reduced and altered effector function that can be partially reversed upon blockade of inhibitory receptors. The role of the exhaustion program and transcriptional networks that control CD8+ T cell function and fate in autoimmunity is not clear. Here we show that intra-islet CD8+ T cells phenotypically, transcriptionally, epigenetically and metabolically possess features of canonically exhausted T cells, yet maintain important differences. This 'restrained' phenotype can be perturbed and disease accelerated by CD8+ T cell-restricted deletion of the inhibitory receptor lymphocyte activating gene 3 (LAG3). Mechanistically, LAG3-deficient CD8+ T cells have enhanced effector-like functions, trafficking to the islets, and have a diminished exhausted phenotype, highlighting a physiological role for an exhaustion program in limiting autoimmunity and implicating LAG3 as a target for autoimmune therapy.

Modulation of tissue resident memory T cells by glucocorticoids after acute cellular rejection in lung transplantation.

Acute cellular rejection is common after lung transplantation and is associated with an increased risk of early chronic rejection. We present combined single-cell RNA and TCR sequencing on recipient-derived T cells obtained from the bronchoalveolar lavage of three lung transplant recipients with rejection and compare them with T cells obtained from the same patients after treatment of rejection with high-dose systemic glucocorticoids. At the time of rejection, we found an oligoclonal expansion of cytotoxic CD8+ T cells that all persisted as tissue resident memory T cells after successful treatment. Persisting CD8+ allograft-resident T cells have reduced gene expression for cytotoxic mediators after therapy with glucocorticoids but accumulate around airways. This clonal expansion is discordant with circulating T cell clonal expansion at the time of rejection, suggesting in situ expansion. We thus highlight the accumulation of cytotoxic, recipient-derived tissue resident memory T cells within the lung allograft that persist despite the administration of high-dose systemic glucocorticoids. The long-term clinical consequences of this persistence have yet to be characterized.

Expansion of Fcγ Receptor IIIa-Positive Macrophages, Ficolin 1-Positive Monocyte-Derived Dendritic Cells, and Plasmacytoid Dendritic Cells Associated With Severe Skin Disease in Systemic Sclerosis.

OBJECTIVE: In this study, we sought a comprehensive understanding of myeloid cell types driving fibrosis in diffuse cutaneous systemic sclerosis (dcSSc) skin. METHODS: We analyzed the transcriptomes of 2,465 myeloid cells from skin biopsy specimens from 12 dcSSc patients and 10 healthy control subjects using single-cell RNA sequencing. Monocyte-derived dendritic cells (mo-DCs) were assessed using immunohistochemical staining and immunofluorescence analyses targeting ficolin-1 (FCN-1). RESULTS: A t-distributed stochastic neighbor embedding analysis of single-cell transcriptome data revealed 12 myeloid cell clusters, 9 of which paralleled previously described healthy control macrophage/DC clusters, and 3 of which were dcSSc-specific myeloid cell clusters. One SSc-associated macrophage cluster, highly expressing Fcγ receptor IIIA, was suggested on pseudotime analysis to be derived from normal CCR1+ and MARCO+ macrophages. A second SSc-associated myeloid population highly expressed monocyte markers FCN-1, epiregulin, S100A8, and S100A9, but was closely related to type 2 conventional DCs on pseudotime analysis and identified as mo-DCs. Mo-DCs were associated with more severe skin disease. Proliferating macrophages and plasmacytoid DCs were detected almost exclusively in dcSSc skin, the latter clustering with B cells and apparently derived from lymphoid progenitors. CONCLUSION: Transcriptional signatures in these and other myeloid populations indicate innate immune system activation, possibly through Toll-like receptors and highly up-regulated chemokines. However, the appearance and activation of myeloid cells varies between patients, indicating potential differences in the underlying pathogenesis and/or temporal disease activity in dcSSc.