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.

Adaptive plasticity of IL-10+ and IL-35+ Treg cells cooperatively promotes tumor T cell exhaustion.

Regulatory T cells (Treg cells) maintain host self-tolerance but are a major barrier to effective cancer immunotherapy. Treg cells subvert beneficial anti-tumor immunity by modulating inhibitory receptor expression on tumor-infiltrating lymphocytes (TILs); however, the underlying mediators and mechanisms have remained elusive. Here, we found that the cytokines IL-10 and IL-35 (Ebi3-IL-12α heterodimer) were divergently expressed by Treg cell subpopulations in the tumor microenvironment (TME) and cooperatively promoted intratumoral T cell exhaustion by modulating several inhibitory receptor expression and exhaustion-associated transcriptomic signature of CD8+ TILs. While expression of BLIMP1 (encoded by Prdm1) was a common target, IL-10 and IL-35 differentially affected effector T cell versus memory T cell fates, respectively, highlighting their differential, partially overlapping but non-redundant regulation of anti-tumor immunity. Our results reveal previously unappreciated cooperative roles for Treg cell-derived IL-10 and IL-35 in promoting BLIMP1-dependent exhaustion of CD8+ TILs that limits effective anti-tumor immunity.

Immune Landscape of Viral- and Carcinogen-Driven Head and Neck Cancer.

Head and neck squamous cell carcinoma (HNSCC) arises through exposure to environmental carcinogens or malignant transformation by human papillomavirus (HPV). Here, we assessed the transcriptional profiles of 131,224 single cells from peripheral and intra-tumoral immune populations from patients with HPV- and HPV+ HNSCC and healthy donors. Immune cells within tumors of HPV- and HPV+ HNSCC displayed a spectrum of transcriptional signatures, with helper CD4+ T cells and B cells being relatively divergent and CD8+ T cells and CD4+ regulatory T cells being relatively similar. Transcriptional results were contextualized through multispectral immunofluorescence analyses and evaluating putative cell-cell communication based on spatial proximity. These analyses defined a gene expression signature associated with CD4+ T follicular helper cells that is associated with longer progression-free survival in HNSCC patients. The datasets and analytical approaches herein provide a resource for the further study of the impact of immune cells on viral- and carcinogen-induced cancers.

Oncolytic Viruses Engineered to Enforce Leptin Expression Reprogram Tumor-Infiltrating T Cell Metabolism and Promote Tumor Clearance.

Immunotherapy can reinvigorate dormant responses to cancer, but response rates remain low. Oncolytic viruses, which replicate in cancer cells, induce tumor lysis and immune priming, but their immune consequences are unclear. We profiled the infiltrate of aggressive melanomas induced by oncolytic Vaccinia virus using RNA sequencing and found substantial remodeling of the tumor microenvironment, dominated by effector T cell influx. However, responses to oncolytic viruses were incomplete due to metabolic insufficiencies induced by the tumor microenvironment. We identified the adipokine leptin as a potent metabolic reprogramming agent that supported antitumor responses. Leptin metabolically reprogrammed T cells in vitro, and melanoma cells expressing leptin were immunologically controlled in mice. Engineering oncolytic viruses to express leptin in tumor cells induced complete responses in tumor-bearing mice and supported memory development in the tumor infiltrate. Thus, leptin can provide metabolic support to tumor immunity, and oncolytic viruses represent a platform to deliver metabolic therapy.

Treg Cells Promote the SREBP1-Dependent Metabolic Fitness of Tumor-Promoting Macrophages via Repression of CD8+ T Cell-Derived Interferon-γ.

Regulatory T (Treg) cells are crucial for immune homeostasis, but they also contribute to tumor immune evasion by promoting a suppressive tumor microenvironment (TME). Mice with Treg cell-restricted Neuropilin-1 deficiency show tumor resistance while maintaining peripheral immune homeostasis, thereby providing a controlled system to interrogate the impact of intratumoral Treg cells on the TME. Using this and other genetic models, we showed that Treg cells shaped the transcriptional landscape across multiple tumor-infiltrating immune cell types. Treg cells suppressed CD8+ T cell secretion of interferon-γ (IFNγ), which would otherwise block the activation of sterol regulatory element-binding protein 1 (SREBP1)-mediated fatty acid synthesis in immunosuppressive (M2-like) tumor-associated macrophages (TAMs). Thus, Treg cells indirectly but selectively sustained M2-like TAM metabolic fitness, mitochondrial integrity, and survival. SREBP1 inhibition augmented the efficacy of immune checkpoint blockade, suggesting that targeting Treg cells or their modulation of lipid metabolism in M2-like TAMs could improve cancer immunotherapy.

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.

Molecular Changes Implicate Angiogenesis and Arterial Remodeling in Systemic Sclerosis-Associated and Idiopathic Pulmonary Hypertension.

BACKGROUND: Pulmonary hypertension (PH) is a common complication of systemic sclerosis (SSc) and a leading cause of mortality among patients with this disease. PH can also occur as an idiopathic condition (idiopathic pulmonary arterial hypertension). Investigation of transcriptomic alterations in vascular populations is critical to elucidating cellular mechanisms underlying pathobiology of SSc-associated and idiopathic PH. METHODS: We analyzed single-cell RNA sequencing profiles of endothelial and perivascular mesenchymal populations from explanted lung tissue of patients with SSc-associated PH (n=16), idiopathic pulmonary arterial hypertension (n=3), and healthy controls (n=15). Findings were validated by immunofluorescence staining of explanted human lung tissue. RESULTS: Three disease-associated endothelial populations emerged. Two angiogenic endothelial cell (EC) subtypes markedly expanded in SSc-associated PH lungs: tip ECs expressing canonical tip markers PGF and APLN and phalanx ECs expressing genes associated with vascular development, endothelial barrier integrity, and Notch signaling. Gene regulatory network analysis suggested enrichment of Smad1 and PPAR-γ (peroxisome proliferator-activated receptor-γ) regulon activities in these 2 populations, respectively. Mapping of potential ligand-receptor interactions highlighted Notch, apelin-APJ, and angiopoietin-Tie signaling pathways between angiogenic ECs and perivascular cells. Transitional cells, expressing both endothelial and pericyte/smooth muscle cell markers, provided evidence for the presence of endothelial-to-mesenchymal transition. Transcriptional programs associated with arterial endothelial dysfunction implicated VEGF-A (vascular endothelial growth factor-A), TGF-ß1, angiotensin, and TNFSF12/TWEAK in the injury/remodeling phenotype of PH arterial ECs. CONCLUSIONS: These data provide high-resolution insights into the complexity and plasticity of the pulmonary endothelium in SSc-associated PH and idiopathic pulmonary arterial hypertension and provide direct molecular insights into soluble mediators and transcription factors driving PH vasculopathy.

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.

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.

Single-cell transcriptomes and chromatin accessibility of endothelial cells unravel transcription factors associated with dysregulated angiogenesis in systemic sclerosis.

OBJECTIVES: Vasculopathy emerges early in systemic sclerosis (SSc) and links to endothelial cell (EC) injury and angiogenesis. Understanding EC transcriptomes and epigenomes is crucial for unravelling the mechanisms involved. METHODS: Transcriptomes and chromatin accessibility were assessed by single-cell RNA sequencing and single-nucleus transposase-accessible chromatin sequencing. Immunofluorescent staining of skin and proteomics assay were employed to confirm the altered SSc EC phenotypes. Gain-of-function assay was used to evaluate the effects of ETS transcription factors on human dermal ECs (hDECs). RESULTS: Both control and SSc ECs shared transcriptomic signatures of vascular linages (arterial, capillary and venous ECs) and lymphatic ECs. Arterial ECs in SSc showed reduced number and increased expression of genes associated with apoptosis. Two distinct EC subpopulations, tip and proliferating ECs, were markedly upregulated in SSc, indicating enhanced proangiogenic and proliferative activities. Molecular features of aberrant SSc-ECs were associated with disease pathogenesis and clinical traits of SSc, such as skin fibrosis and digital ulcers. Ligand-receptor analysis demonstrated altered intercellular networks of SSc EC subpopulations with perivascular and immune cells. Furthermore, the integration of open chromatin profiles with transcriptomic analysis suggested an increased accessibility of regulatory elements for ETS family transcription factors in SSc ECs. Overexpression of ETS genes in hDECs suggested ELK4, ERF and ETS1 may orchestrate arterial apoptosis and dysregulated angiogenesis in SSc. CONCLUSIONS: This study unveils transcriptional and chromatin alterations in driving endovascular dysregulation in SSc, proposing ELK4, ERF and ETS1 as novel targets in ECs for addressing vascular complications in the condition.

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.

Long non-coding RNA H19X as a regulator of mononuclear cell adhesion to the endothelium in systemic sclerosis.

OBJECTIVE: To define the functional relevance of H19 X-linked co-expressed lncRNA (H19X) in endothelial cell (EC) activation as a key process in systemic sclerosis (SSc) vasculopathy. METHODS: H19X expression in SSc skin biopsies was analyzed from single cell RNA sequencing (scRNA-seq) data. Differential expression and pathway enrichment analysis between cells expressing (H19Xpos) and non expressing H19X (H19Xneg) cells was performed. H19X function was investigated in human dermal microvascular EC (HDMECs) by silencing. H19X and EC adhesion molecules levels were analyzed by RT-qPCR and Western Blot after stimulation with proinflammatory cytokines. Cytoskeletal rearrangements were analyzed by fluorescent staining. Endothelial adhesion was evaluated by co-culture of HDMECs and fluorescent labelled peripheral blood mononuclear cells (PBMCs). Shedding VCAM1 was evaluated by ELISA on HDMEC supernatant. RESULTS: scRNA-seq showed significant upregulation of H19X in SSc compared with healthy EC. In HDMEC, H19X was consistently induced by type I and II interferons. H19X knockdown lead to a significant decrease of the mRNA of several adhesion molecules. Particularly, vascular cell adhesion protein 1 (VCAM1) was significantly reduced at protein and mRNA levels. Co-expression analysis of the scRNA-seq data confirmed a higher expression of VCAM1 in (H19Xpos) EC. EC were also strongly associated with the 'cell adhesion molecule' pathway. Moreover, VCAM1 downstream pathway displayed less activation following H19X knockdown. Contractility of HDMEC, PBMC adhesion to HDMEC and VCAM1 shedding were also reduced following H19X knockdown. CONCLUSIONS: lncRNA H19X may contribute to EC activation in SSc vasculopathy, acting as a regulator of expression of adhesion molecules in EC.

Loss of Protein Kinase D2 Activity Protects Against Bleomycin-Induced Dermal Fibrosis in Mice.

Protein kinase D (PKD) has been linked to inflammatory responses in various pathologic conditions; however, its role in inflammation-induced dermal fibrosis has not been evaluated. In this study, we aimed to investigate the roles and mechanisms of protein kinase D2 (PKD2) in inflammation-induced dermal fibrosis and evaluate the therapeutic potential of PKD inhibitors in this disease. Using homozygous kinase-dead PKD2 knock-in (KI) mice, we examined whether genetic ablation or pharmacologic inhibition of PKD2 activity affected dermal inflammation and fibrosis in a bleomycin (BLM)-induced skin fibrosis model. Our data showed that dermal thickness and collagen fibers were significantly reduced in BLM-treated PKD2 KI mice compared with that in wild-type mice, and so was the expression of α-smooth muscle actin and collagens and the mRNA levels of transforming growth factor-ß1 and interleukin-6 in the KI mice. Corroboratively, pharmacologic inhibition of PKD by CRT0066101 also significantly blocked BLM-induced dermal fibrosis and reduced α-smooth muscle actin, collagen, and interleukin-6 expression. Further analyses indicated that loss of PKD2 activity significantly blocked BLM-induced infiltration of monocytes/macrophages and neutrophils in the dermis. Moreover, using bone marrow-derived macrophages, we demonstrated that PKD activity was required for cytokine production and migration of macrophages. We have further identified Akt as a major downstream target of PKD2 in the early inflammatory phase of the fibrotic process. Taken together, our findings indicate that PKD2 promotes dermal fibrosis via regulating immune cell infiltration, cytokine production, and downstream activation of Akt in lesional skin, and targeted inhibition of PKD2 may benefit the treatment of this condition.

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.

Wnt7a deficit is associated with dysfunctional angiogenesis in pulmonary arterial hypertension.

INTRODUCTION: Pulmonary arterial hypertension (PAH) is characterised by loss of microvessels. The Wnt pathways control pulmonary angiogenesis but their role in PAH is incompletely understood. We hypothesised that Wnt activation in pulmonary microvascular endothelial cells (PMVECs) is required for pulmonary angiogenesis, and its loss contributes to PAH. METHODS: Lung tissue and PMVECs from healthy and PAH patients were screened for Wnt production. Global and endothelial-specific Wnt7a -/- mice were generated and exposed to chronic hypoxia and Sugen-hypoxia (SuHx). RESULTS: Healthy PMVECs demonstrated >6-fold Wnt7a expression during angiogenesis that was absent in PAH PMVECs and lungs. Wnt7a expression correlated with the formation of tip cells, a migratory endothelial phenotype critical for angiogenesis. PAH PMVECs demonstrated reduced vascular endothelial growth factor (VEGF)-induced tip cell formation as evidenced by reduced filopodia formation and motility, which was partially rescued by recombinant Wnt7a. We discovered that Wnt7a promotes VEGF signalling by facilitating Y1175 tyrosine phosphorylation in vascular endothelial growth factor receptor 2 (VEGFR2) through receptor tyrosine kinase-like orphan receptor 2 (ROR2), a Wnt-specific receptor. We found that ROR2 knockdown mimics Wnt7a insufficiency and prevents recovery of tip cell formation with Wnt7a stimulation. While there was no difference between wild-type and endothelial-specific Wnt7a -/- mice under either chronic hypoxia or SuHx, global Wnt7a +/- mice in hypoxia demonstrated higher pulmonary pressures and severe right ventricular and lung vascular remodelling. Similar to PAH, Wnt7a +/- PMVECs exhibited an insufficient angiogenic response to VEGF-A that improved with Wnt7a. CONCLUSIONS: Wnt7a promotes VEGF signalling in lung PMVECs and its loss is associated with an insufficient VEGF-A angiogenic response. We propose that Wnt7a deficiency contributes to progressive small vessel loss in PAH.

Biomarkers of haemodynamic severity of systemic sclerosis-associated pulmonary arterial hypertension by serum proteome analysis.

OBJECTIVES: To mine the serum proteome of patients with systemic sclerosis-associated pulmonary arterial hypertension (SSc-PAH) and to detect biomarkers that may assist in earlier and more effective diagnosis and treatment. METHODS: Patients with limited cutaneous SSc, no extensive interstitial lung disease and no PAH-specific therapy were included. They were classified as cases if they had PAH confirmed by right heart catheterisation (RHC) and serum collected on the same day as RHC; and as controls if they had no clinical evidence of PAH. RESULTS: Patients were mostly middle-aged females with anticentromere-associated SSc. Among 1129 proteins assessed by a high-throughput proteomic assay (SOMAscan), only 2 were differentially expressed and correlated significantly with pulmonary vascular resistance (PVR) in SSc-PAH patients (n=15): chemerin (ρ=0.62, p=0.01) and SET (ρ=0.62, p=0.01). To validate these results, serum levels of chemerin were measured by ELISA in an independent cohort. Chemerin levels were confirmed to be significantly higher (p=0.01) and correlate with PVR (ρ=0.42, p=0.04) in SSc-PAH patients (n=24). Chemerin mRNA expression was detected in fibroblasts, pulmonary artery smooth muscle cells (PA-SMCs)/pericytes and mesothelial cells in SSc-PAH lungs by single-cell RNA-sequencing. Confocal immunofluorescence revealed increased expression of a chemerin receptor, CMKLR1, on SSc-PAH PA-SMCs. SSc-PAH serum seemed to induce higher PA-SMC proliferation than serum from SSc patients without PAH. This difference appeared neutralised when adding the CMKLR1 inhibitor α-NETA. CONCLUSION: Chemerin seems an interesting surrogate biomarker for PVR in SSc-PAH. Increased chemerin serum levels and CMKLR1 expression by PA-SMCs may contribute to SSc-PAH pathogenesis by inducing PA-SMC proliferation.

Impact of enzymatic digestion on single cell suspension yield from peripheral human lung tissue.

An increasing number of translational investigations of lung biology rely on analyzing single cell suspensions obtained from human lungs. To obtain these single cell suspensions, human lungs from biopsies or research-consented organ donors must be subjected to mechanical and enzymatic digestion prior to analysis with either flow cytometry or single cell RNA sequencing. A variety of enzymes have been used to perform tissue digestion, each with potential limitations. To better understand the limitations of each enzymatic digestion protocol and to establish a framework for comparing studies across protocols, we performed five commonly published protocols in parallel from identical samples obtained from 6 human lungs. Following mechanical (gentleMACS™) and enzymatic digestion, we quantified cell count and viability using a Nexcelom Cellometer and determined cell phenotype using multiparameter spectral flow cytometry (Cytek™ Aurora). We found that all protocols were superior in cellular yield and viability when compared to mechanical digestion alone. Protocols high in dispase cleaved immune markers CD4, CD8, CD69, and CD103 and contributed to an increased monocyte to macrophage yield. Similarly, dispase led to a differential epithelial cell yield, with increased TSPN8+ and ITGA6+ epithelial cells and reduced CD66e+ cells. When compared to collagenase D, collagenase P protocols yielded increased AT1 and AT2 cells and decreased endothelial cells. These results provide a framework for selecting an enzymatic digestion protocol best suited to the scientific question and allow for comparison of studies using different protocols.

A data-driven approach finds RNA polymerase III antibody and tendon friction rubs as enrichment tools for early diffuse scleroderma trials.

OBJECTIVE: Clinical trials in early diffuse SSc have consistently shown a placebo group response with a declining modified Rodnan skin score (mRSS), with negative outcomes. Our objective was to identify strategies using clinical characteristics or laboratory values to improve trial design. METHODS: We identified early diffuse SSc patients first seen at the University of Pittsburgh from 1980-2015. Eligible patients had ≥3 visits, with at least two mRSS scores within the first year of follow-up. We performed Kaplan-Meier analyses, group-based trajectory analysis of mRSS scores, followed by multivariable regression analysis and classification tree analysis. We applied the results to the abatacept in early diffuse systemic sclerosis (ASSET) trial outcome data. RESULTS: We identified 403 patients with <18 months, and 514 with <36 months disease duration. The median number of mRSS follow-up scores was 14 (interquartile range 8, 25). All methodologic approaches identified skin thickness progression rate, RNA polymerase III (RNAP3) antibody positivity and presence of tendon friction rubs (TFR) as predictors of mRSS trajectory over 5 years of follow-up, and thereby as potential enrichment variables. When applied to the ASSET data, adjustment for both RNAP3 and TFR demonstrated reduction of the placebo mRSS response, particularly at 6 months. A significant difference in the ACR Composite Response Index in Systemic Sclerosis (CRISS) score was found with adjustment by RNAP3 at 6 months, and TFR or RNAP3 at 12 months. CONCLUSION: Adjustment for both RNAP3 and TFR predicts mRSS trajectory and diminished the mRSS decline in ASSET placebo group, and identified significant differences in CRISS. RNAP3, particularly, is a stratification or enrichment approach to improve early diffuse SSc trial design.

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.