Stricturing Crohn's Disease Single-Cell RNA Sequencing Reveals Fibroblast Heterogeneity and Intercellular Interactions.

BACKGROUND & AIMS: Fibroblasts play a key role in stricture formation in Crohn's disease (CD) but understanding its pathogenesis requires a systems-level investigation to uncover new treatment targets. We studied full-thickness CD tissues to characterize fibroblast heterogeneity and function by generating the first single-cell RNA sequencing (scRNAseq) atlas of strictured bowel and providing proof of principle for therapeutic target validation. METHODS: We performed scRNAseq of 13 fresh full-thickness CD resections containing noninvolved, inflamed nonstrictured, and strictured segments as well as 7 normal non-CD bowel segments. Each segment was separated into mucosa/submucosa or muscularis propria and analyzed separately for a total of 99 tissue samples and 409,001 cells. We validated cadherin-11 (CDH11) as a potential therapeutic target by using whole tissues, isolated intestinal cells, NanoString nCounter, next-generation sequencing, proteomics, and animal models. RESULTS: Our integrated dataset revealed fibroblast heterogeneity in strictured CD with the majority of stricture-selective changes detected in the mucosa/submucosa, but not the muscle layer. Cell-cell interaction modeling revealed CXCL14+ as well as MMP/WNT5A+ fibroblasts displaying a central signaling role in CD strictures. CDH11, a fibroblast cell-cell adhesion molecule, was broadly expressed and up-regulated, and its profibrotic function was validated using NanoString nCounter, RNA sequencing, tissue target expression, in vitro gain- and loss-of-function experiments, proteomics, and knock-out and antibody-mediated CDH11 blockade in experimental colitis. CONCLUSIONS: A full-thickness bowel scRNAseq atlas revealed previously unrecognized fibroblast heterogeneity and interactions in CD strictures and CDH11 was validated as a potential therapeutic target. These results provide a new resource for a better understanding of CD stricture formation and open potential therapeutic developments. This work has been posted as a preprint on Biorxiv under doi: 10.1101/2023.04.03.534781.

Immune Complex-Driven Generation of Human Macrophages with Anti-Inflammatory and Growth-Promoting Activity.

To maintain homeostasis, macrophages must be capable of assuming either an inflammatory or an anti-inflammatory phenotype. To better understand the latter, we stimulated human macrophages in vitro with TLR ligands in the presence of high-density immune complexes (IC). This combination of stimuli resulted in a broad suppression of inflammatory mediators and an upregulation of molecules involved in tissue remodeling and angiogenesis. Transcriptomic analysis of TLR stimulation in the presence of IC predicted the downstream activation of AKT and the inhibition of GSK3. Consequently, we pretreated LPS-stimulated human macrophages with small molecule inhibitors of GSK3 to partially phenocopy the regulatory effects of stimulation in the presence of IC. The upregulation of DC-STAMP and matrix metalloproteases was observed on these cells and may represent potential biomarkers for this regulatory activation state. To demonstrate the presence of these anti-inflammatory, growth-promoting macrophages in a human infectious disease, biopsies from patients with leprosy (Hanseniasis) were analyzed. The lepromatous form of this disease is characterized by hypergammaglobulinemia and defective cell-mediated immunity. Lesions in lepromatous leprosy contained macrophages with a regulatory phenotype expressing higher levels of DC-STAMP and lower levels of IL-12, relative to macrophages in tuberculoid leprosy lesions. Therefore, we propose that increased signaling by FcγR cross-linking on TLR-stimulated macrophages can paradoxically promote the resolution of inflammation and initiate processes critical to tissue growth and repair. It can also contribute to infectious disease progression.

Macrophages and the Recovery from Acute and Chronic Inflammation.

In recent years, researchers have devoted much attention to the diverse roles of macrophages and their contributions to tissue development, wound healing, and angiogenesis. What should not be lost in the discussions regarding the diverse biology of these cells is that when perturbed, macrophages are the primary contributors to potentially pathological inflammatory processes. Macrophages stand poised to rapidly produce large amounts of inflammatory cytokines in response to danger signals. The production of these cytokines can initiate a cascade of inflammatory mediator release that can lead to wholesale tissue destruction. The destructive inflammatory capability of macrophages is amplified by exposure to exogenous interferon-γ, which prolongs and heightens inflammatory responses. In simple terms, macrophages can thus be viewed as incendiary devices with hair triggers waiting to detonate. We have begun to ask questions about how these cells can be regulated to mitigate the collateral destruction associated with macrophage activation.

The influence of HLA genetic variation on plasma protein expression.

Genetic variation in the human leukocyte antigen (HLA) loci is associated with risk of immune-mediated diseases, but the molecular effects of HLA polymorphism are unclear. Here we examined the effects of HLA genetic variation on the expression of 2940 plasma proteins across 45,330 Europeans in the UK Biobank, with replication analyses across multiple ancestry groups. We detected 504 proteins affected by HLA variants (HLA-pQTL), including widespread trans effects by autoimmune disease risk alleles. More than 80% of the HLA-pQTL fine-mapped to amino acid positions in the peptide binding groove. HLA-I and II affected proteins expressed in similar cell types but in different pathways of both adaptive and innate immunity. Finally, we investigated potential HLA-pQTL effects on disease by integrating HLA-pQTL with fine-mapped HLA-disease signals in the UK Biobank. Our data reveal the diverse effects of HLA genetic variation and aid the interpretation of associations between HLA alleles and immune-mediated diseases.

Single-cell isolation from full-thickness human intestinal tissue resections for single-cell RNA sequencing.

Single-cell isolation techniques allow the investigation of physical and functional relationships between individual cells within a complex cell population. Here, we present a protocol for single-cell isolation from full-thickness intestinal tissue resections. We describe steps for pre-processing specimens, isolation of lamina propria and muscular layers, and red blood cell lysis. We then detail fixation of isolated cells and assessment of cell quality. The resulting cell suspension can be subjected to RNA sequencing on the 10× Chromium platform. For complete details on the use and execution of this protocol, please refer to Mukherjee et al.1.

Identification of a broadly fibrogenic macrophage subset induced by type 3 inflammation.

Macrophages are central orchestrators of the tissue response to injury, with distinct macrophage activation states playing key roles in fibrosis progression and resolution. Identifying key macrophage populations found in human fibrotic tissues could lead to new treatments for fibrosis. Here, we used human liver and lung single-cell RNA sequencing datasets to identify a subset of CD9+TREM2+ macrophages that express SPP1, GPNMB, FABP5, and CD63. In both human and murine hepatic and pulmonary fibrosis, these macrophages were enriched at the outside edges of scarring and adjacent to activated mesenchymal cells. Neutrophils expressing MMP9, which participates in the activation of TGF-ß1, and the type 3 cytokines GM-CSF and IL-17A coclustered with these macrophages. In vitro, GM-CSF, IL-17A, and TGF-ß1 drive the differentiation of human monocytes into macrophages expressing scar-associated markers. Such differentiated cells could degrade collagen IV but not collagen I and promote TGF-ß1-induced collagen I deposition by activated mesenchymal cells. In murine models blocking GM-CSF, IL-17A or TGF-ß1 reduced scar-associated macrophage expansion and hepatic or pulmonary fibrosis. Our work identifies a highly specific macrophage population to which we assign a profibrotic role across species and tissues. It further provides a strategy for unbiased discovery, triage, and preclinical validation of therapeutic targets based on this fibrogenic macrophage population.

Distinct use of super-enhancer elements controls cell type-specific CD25 transcription and function.

The IL-2 receptor α chain (IL-2Rα/CD25) is constitutively expressed on double-negative (DN2/DN3 thymocytes and regulatory T cells (Tregs) but induced by IL-2 on T and natural killer (NK) cells, with Il2ra expression regulated by a STAT5-dependent super-enhancer. We investigated CD25 regulation and function using a series of mice with deletions spanning STAT5-binding elements. Deleting the upstream super-enhancer region mainly affected constitutive CD25 expression on DN2/DN3 thymocytes and Tregs, with these mice developing autoimmune alopecia, whereas deleting an intronic region decreased IL-2-induced CD25 on peripheral T and NK cells. Thus, distinct super-enhancer elements preferentially control constitutive versus inducible expression in a cell type-specific manner. The mediator-1 coactivator colocalized with specific STAT5-binding sites. Moreover, both upstream and intronic regions had extensive chromatin interactions, and deletion of either region altered the super-enhancer structure in mature T cells. These results demonstrate differential functions for distinct super-enhancer elements, thereby indicating previously unknown ways to manipulate CD25 expression in a cell type-specific fashion.

Stricturing Crohn's disease single-cell RNA sequencing reveals fibroblast heterogeneity and intercellular interactions.

Background: Fibroblasts play a key role in stricture formation in Crohn's disease (CD) but understanding it's pathogenesis requires a systems-level investigation to uncover new treatment targets. We studied full thickness CD tissues to characterize fibroblast heterogeneity and function by generating the first single cell RNA sequencing (scRNAseq) atlas of strictured bowel and providing proof of principle for therapeutic target validation. Methods: We performed scRNAseq of 13 fresh full thickness CD resections containing non-involved, inflamed non-strictured, and strictured segments as well as 7 normal non-CD bowel segments. Each segment was separated into mucosa/submucosa or muscularis propria and analyzed separately for a total of 99 tissue samples and 409,001 cells. We validated cadherin-11 (CDH11) as a potential therapeutic target by using whole tissues, isolated intestinal cells, NanoString nCounter, next generation sequencing, proteomics and animal models. Results: Our integrated dataset revealed fibroblast heterogeneity in strictured CD with the majority of stricture-selective changes detected in the mucosa/submucosa, but not the muscle layer. Cell-cell interaction modeling revealed CXCL14+ as well as MMP/WNT5A+ fibroblasts displaying a central signaling role in CD strictures. CDH11, a fibroblast cell-cell adhesion molecule, was broadly expressed and upregulated, and its pro-fibrotic function was validated by NanoString nCounter, RNA sequencing, tissue target expression, in vitro gain- and loss-of-function experiments, proteomics, and two animal models of experimental colitis. Conclusion: A full-thickness bowel scRNAseq atlas revealed previously unrecognized fibroblast heterogeneity and interactions in CD strictures and CDH11 was validated as a potential therapeutic target. These results provide a new resource for a better understanding of CD stricture formation and opens potential therapeutic developments.

IFNγ regulates NAD+ metabolism to promote the respiratory burst in human monocytes.

Interferon γ (IFNγ) is an essential and pleiotropic activator of human monocytes, but little is known about the changes in cellular metabolism required for IFNγ-induced activation. We sought to elucidate the mechanisms by which IFNγ reprograms monocyte metabolism to support its immunologic activities. We found that IFNγ increased oxygen consumption rates (OCR) in monocytes, indicative of reactive oxygen species generation by both mitochondria and nicotinamide adenine dinucleotide phosphate (NADPH) oxidase. Transcriptional profiling revealed that this oxidative phenotype was driven by IFNγ-induced reprogramming of NAD+ metabolism, which is dependent on nicotinamide phosphoribosyltransferase (NAMPT)-mediated NAD+ salvage to generate NADH and NADPH for oxidation by mitochondrial complex I and NADPH oxidase, respectively. Consistent with this pathway, monocytes from patients with gain-of-function mutations in STAT1 demonstrated higher-than-normal OCR, whereas chemical or genetic disruption of mitochondrial complex I (rotenone treatment or Leigh syndrome patient monocytes) or NADPH oxidase (diphenyleneiodonium treatment or chronic granulomatous disease [CGD] patient monocytes) reduced OCR. Interestingly, inhibition of NAMPT in healthy monocytes completely abrogated the IFNγ-induced oxygen consumption, comparable to levels observed in CGD monocytes. These data identify an IFNγ-induced, NAMPT-dependent, NAD+ salvage pathway that is critical for IFNγ activation of human monocytes.

Single cell sequencing identifies clonally expanded synovial CD4+ TPH cells expressing GPR56 in rheumatoid arthritis.

Rheumatoid arthritis (RA) is an autoimmune disease affecting synovial joints where different CD4+ T cell subsets may contribute to pathology. Here, we perform single cell sequencing on synovial CD4+ T cells from anti-citrullinated protein antibodies (ACPA)+ and ACPA- RA patients and identify two peripheral helper T cell (TPH) states and a cytotoxic CD4+ T cell subset. We show that the adhesion G-protein coupled receptor 56 (GPR56) delineates synovial CXCL13high TPH CD4+ T cells expressing LAG-3 and the tissue-resident memory receptors CXCR6 and CD69. In ACPA- SF, TPH cells display lower levels of GPR56 and LAG-3. Further, most expanded T cell clones in the joint are within CXCL13high TPH CD4+ T cells. Finally, RNA-velocity analyses suggest a common differentiation pathway between the two TPH clusters and effector CD4+ T cells. Our study provides comprehensive immunoprofiling of the synovial CD4+ T cell subsets in ACPA+ and ACPA- RA.

Humoral immunity in leishmaniasis - Prevention or promotion of parasite growth?

Leishmaniasis can present as a "spectrum" of clinical outcomes. There is evidence that these divergent clinical outcomes are attributable to genetic differences in the human host [1] as well the species of infecting parasite [2]. The spectrum of disease has largely been described by defining the polar opposites of T cell immune responses. In the mouse model, a TH1 immune response is associated with low numbers of Leishmania parasites in lesions, whereas a TH2 immune response has been associated with unrestricted parasite growth. In the present work, we revisit leishmaniasis and seek to better define the clinical spectrum as a function of divergent humoral immune responses. We describe examples in human, canine, and even some murine models of leishmaniasis that reveal a direct correlation between high anti-parasite antibody responses and unrestricted parasite growth. Therefore, we propose that the spectral nature of this disease may be due to quantitative and qualitative differences in the antibodies that are produced during disease. In human visceral leishmaniasis, a decrease in anti-parasite antibody levels may actually predict disease resolution. Thus, rather than defining this disease as a simple TH1/TH2 dichotomy, we propose that clinical leishmaniasis depends on the degree of humoral immunity, with high IgG predicting parasite persistence. These observations have obvious implications for vaccine development in leishmaniasis, and they may extend to other diseases caused by intracellular pathogens.

Host and parasite responses in human diffuse cutaneous leishmaniasis caused by L. amazonensis.

Diffuse cutaneous leishmaniasis (DCL) is a rare form of leishmaniasis where parasites grow uncontrolled in diffuse lesions across the skin. Meta-transcriptomic analysis of biopsies from DCL patients infected with Leishmania amazonensis demonstrated an infiltration of atypical B cells producing a surprising preponderance of the IgG4 isotype. DCL lesions contained minimal CD8+ T cell transcripts and no evidence of persistent TH2 responses. Whereas localized disease exhibited activated (so-called M1) macrophage presence, transcripts in DCL suggested a regulatory macrophage (R-Mϕ) phenotype with higher levels of ABCB5, DCSTAMP, SPP1, SLAMF9, PPARG, MMPs, and TM4SF19. The high levels of parasite transcripts in DCL and the remarkable uniformity among patients afforded a unique opportunity to study parasite gene expression in this disease. Patterns of parasite gene expression in DCL more closely resembled in vitro parasite growth in resting macrophages, in the absence of T cells. In contrast, parasite gene expression in LCL revealed 336 parasite genes that were differently upregulated, relative to DCL and in vitro macrophage growth, and these transcripts may represent transcripts that are produced by the parasite in response to host immune pressure.

Correction: Meta-transcriptome Profiling of the Human-Leishmania braziliensis Cutaneous Lesion.

[This corrects the article DOI: 10.1371/journal.pntd.0004992.].

Meta-transcriptome Profiling of the Human-Leishmania braziliensis Cutaneous Lesion.

Host and parasite gene expression in skin biopsies from Leishmania braziliensis-infected patients were simultaneously analyzed using high throughput RNA-sequencing. Biopsies were taken from 8 patients with early cutaneous leishmaniasis and 17 patients with late cutaneous leishmaniasis. Although parasite DNA was found in all patient lesions at the time of biopsy, the patients could be stratified into two groups: one lacking detectable parasite transcripts (PTNeg) in lesions, and another in which parasite transcripts were readily detected (PTPos). These groups exhibited substantial differences in host responses to infection. PTPos biopsies contained an unexpected increase in B lymphocyte-specific and immunoglobulin transcripts in the lesions, and an upregulation of immune inhibitory molecules. Biopsies without detectable parasite transcripts showed decreased evidence for B cell activation, but increased expression of antimicrobial genes and genes encoding skin barrier functions. The composition and abundance of L. braziliensis transcripts in PTPos lesions were surprisingly conserved among all six patients, with minimal meaningful differences between lesions from patients with early and late cutaneous leishmaniasis. The most abundant parasite transcripts expressed in lesions were distinct from transcripts expressed in vitro in human macrophage cultures infected with L. amazonensis or L. major. Therefore in vitro gene expression in macrophage monolayers may not be a strong predictor of gene expression in lesions. Some of the most highly expressed in vivo transcripts encoded amastin-like proteins, hypothetical genes, putative parasite virulence factors, as well as histones and tubulin. In summary, RNA sequencing allowed us to simultaneously analyze human and L. braziliensis transcriptomes in lesions of infected patients, and identify unexpected differences in host immune responses which correlated with active transcription of parasite genes.

Extensive rewiring of epithelial-stromal co-expression networks in breast cancer.

BACKGROUND: Epithelial-stromal crosstalk plays a critical role in invasive breast cancer pathogenesis; however, little is known on a systems level about how epithelial-stromal interactions evolve during carcinogenesis. RESULTS: We develop a framework for building genome-wide epithelial-stromal co-expression networks composed of pairwise co-expression relationships between mRNA levels of genes expressed in the epithelium and stroma across a population of patients. We apply this method to laser capture micro-dissection expression profiling datasets in the setting of breast carcinogenesis. Our analysis shows that epithelial-stromal co-expression networks undergo extensive rewiring during carcinogenesis, with the emergence of distinct network hubs in normal breast, and estrogen receptor-positive and estrogen receptor-negative invasive breast cancer, and the emergence of distinct patterns of functional network enrichment. In contrast to normal breast, the strongest epithelial-stromal co-expression relationships in invasive breast cancer mostly represent self-loops, in which the same gene is co-expressed in epithelial and stromal regions. We validate this observation using an independent laser capture micro-dissection dataset and confirm that self-loop interactions are significantly increased in cancer by performing computational image analysis of epithelial and stromal protein expression using images from the Human Protein Atlas. CONCLUSIONS: Epithelial-stromal co-expression network analysis represents a new approach for systems-level analyses of spatially localized transcriptomic data. The analysis provides new biological insights into the rewiring of epithelial-stromal co-expression networks and the emergence of epithelial-stromal co-expression self-loops in breast cancer. The approach may facilitate the development of new diagnostics and therapeutics targeting epithelial-stromal interactions in cancer.