A Novel MIP-1-Expressing Macrophage Subtype in BAL Fluid from Healthy Volunteers.

Tissue availability remains an important limitation of single-cell genomic technologies for investigating cellular heterogeneity in human health and disease. BAL represents a minimally invasive approach to assessing an individual's lung cellular environment for diagnosis and research. However, the lack of high-quality, healthy lung reference data is a major obstacle to using single-cell approaches to study a plethora of lung diseases. Here, we performed single-cell RNA sequencing on over 40,000 cells isolated from the BAL of four healthy volunteers. Of the six cell types or lineages we identified, macrophages were consistently the most numerous across individuals. Our analysis confirmed the expression of marker genes defining cell types despite background signals because of the ambient RNA found in many single-cell studies. We assessed the variability of gene expression across macrophages and defined a distinct subpopulation of cells expressing a set of genes associated with Macrophage Inflammatory Protein 1 (MIP-1). RNA in situ hybridization and reanalysis of published lung single-cell data validated the presence of this macrophage subpopulation. Thus, our study characterizes lung macrophage heterogeneity in healthy individuals and provides a valuable resource for future studies to understand the lung environment in health and disease.

Microglia Adopt Longitudinal Transcriptional Changes After Traumatic Brain Injury.

BACKGROUND: Traumatic brain injury (TBI) is an under-recognized public health threat. Even mild brain injuries can lead to long-term neurologic impairment. Microglia play a fundamental role in the development and progression of this ensuing neurologic impairment. Despite this, a microglia-specific injury signature has yet to be identified. We hypothesized that TBI would lead to long-term changes in the transcriptional profile of microglial pathways associated with the development of subsequent neurologic impairment. MATERIALS AND METHODS: Male C57BL/6 mice underwent TBI via a controlled cortical impact and were followed longitudinally. FACSorted microglia from TBI mice were subjected to Quantiseq 3'-biased RNA sequencing at 7, 30, and 90 d after TBI. K-means clustering on 396 differentially expressed genes was performed, and gene ontology enrichment analysis was used to determine corresponding enriched processes. RESULTS: Differentially expressed genes in microglia exhibited four main patterns of expression over the course of TBI. In particular, we identified four gene clusters which corresponded to the host defense response, synaptic plasticity, lipid remodeling, and membrane polarization. CONCLUSIONS: Transcriptional profiling within individual populations of microglia after TBI remains a critical unmet research need within the field of TBI. This focused study identified several physiologic processes within microglia that may be associated with development of long-term neurologic impairment after TBI. These data demonstrate the capability of longitudinal transcriptional profiling to uncover potential cell-specific targets for the treatment of TBI.

Lipocalin-2 is a pathogenic determinant and biomarker of neuropsychiatric lupus.

Neuropsychiatric manifestations in lupus (NPSLE) affect ∼20-40% of patients. In the central nervous system, lipocalin-2 (LCN2) can promote injury through mechanisms directly linked to NPSLE, including brain barrier disruption, neurotoxicity, and glial activation. Since LCN2 is elevated in lupus and has been implicated in neuroinflammation, we investigated whether LCN2 is required for the pathogenesis of NPSLE. Here, we investigated the effects of LCN2 deficiency on the development of neurobehavioral deficits in the B6.Sle1.Sle3 (Sle1,3) mouse lupus model. Sle1,3 mice exhibited depression-like behavior and impaired spatial and recognition memory, and these deficits were attenuated in Sle1,3-LCN2KO mice. Whole-brain flow cytometry showed a significant increase in brain infiltrating leukocytes in Sle1,3 mice that was not reduced by LCN2 deficiency. RNA sequencing on sorted microglia revealed that several genes differentially expressed between B6 and Sle1,3 mice were regulated by LCN2, and that these genes are key mediators of the neuroinflammatory cascade. Importantly, LCN2 is upregulated in the cerebrospinal fluid of NPSLE patients across 2 different ethnicities. Our findings establish the Sle1,3 strain as an NPSLE model, demonstrate that LCN2 is a major regulator of the detrimental neuroimmune response in NPSLE, and identify CSF LCN2 as a novel biomarker for NPSLE.

Retinal microglia express more MHC class I and promote greater T-cell-driven inflammation than brain microglia.

Introduction: Macrophage function is determined by microenvironment and origin. Brain and retinal microglia are both derived from yolk sac progenitors, yet their microenvironments differ. Utilizing single-cell RNA sequencing (scRNA-seq) data from mice, we tested the hypothesis that retinal and brain microglia exhibit distinct transcriptional profiles due to their unique microenvironments. Methods: Eyes and brains from 2-4 month wildtype mice were combined (20 eyes; 3 brains) to yield one biologically diverse sample per organ. Each tissue was digested into single cell suspensions, enriched for immune cells, and sorted for scRNA-seq. Analysis was performed in Seurat v3 including clustering, integration, and differential expression. Multi-parameter flow cytometry was used for validation of scRNA-seq results. Lymphocytic choriomeningitis virus (LCMV) Clone 13, which produces a systemic, chronic, and neurotropic infection, was used to validate scRNA-seq and flow cytometry results in vivo. Results: Cluster analysis of integrated gene expression data from eye and brain identified 6 Tmem119 + P2ry12 + microglial clusters. Differential expression analysis revealed that eye microglia were enriched for more pro-inflammatory processes including antigen processing via MHC class I (14.0-fold, H2-D1 and H2-K1) and positive regulation of T-cell immunity (8.4-fold) compared to brain microglia. Multi-parameter flow cytometry confirmed that retinal microglia expressed 3.2-fold greater H2-Db and 263.3-fold more H2-Kb than brain microglia. On Day 13 and 29 after LCMV infection, CD8+ T-cell density was greater in the retina than the brain. Discussion: Our data demonstrate that the microenvironment of retina and brain differs, resulting in microglia-specific gene expression changes. Specifically, retinal microglia express greater MHC class I by scRNA-seq and multi-parameter flow cytometry, resulting in a possibly enhanced capability to stimulate CD8+ T-cell inflammation during LCMV infection. These results may explain tissue-specific differences between retina and brain during systemic viral infections and CD8+ T-cell driven autoimmune disease.

MAGNET: A web-based application for gene set enrichment analysis using macrophage data sets.

Characterization of gene lists obtained from high-throughput genomic experiments is an essential task to uncover the underlying biological insights. A common strategy is to perform enrichment analyses that utilize standardized biological annotations, such as GO and KEGG pathways, which attempt to encompass all domains of biology. However, this approach provides generalized, static results that may fail to capture subtleties associated with research questions within a specific domain. Thus, there is a need for an application that can provide precise, relevant results by leveraging the latest research. We have therefore developed an interactive web application, Macrophage Annotation of Gene Network Enrichment Tool (MAGNET), for performing enrichment analyses on gene sets that are specifically relevant to macrophages. Using the hypergeometric distribution, MAGNET assesses the significance of overlapping genes with annotations that were curated from published manuscripts and data repositories. We implemented numerous features that enhance utility and user-friendliness, such as the simultaneous testing of multiple gene sets, different visualization options, option to upload custom datasets, and downloadable outputs. Here, we use three example studies compared against our current database of ten publications on mouse macrophages to demonstrate that MAGNET provides relevant and unique results that complement conventional enrichment analysis tools. Although specific to macrophage datasets, we envision MAGNET will catalyze developments of similar applications in other domains of interest. MAGNET can be freely accessed at the URL https://magnet-winterlab.herokuapp.com. Website implemented in Python and PostgreSQL, with all major browsers supported. The source code is available at https://github.com/sychen9584/MAGNET.

Single-cell sequencing of a novel model of neonatal bile duct ligation in mice identifies macrophage heterogeneity in obstructive cholestasis.

Macrophages (MΦ) play a role in neonatal etiologies of obstructive cholestasis, however, the role for precise MΦ subsets remains poorly defined. We developed a neonatal murine model of bile duct ligation (BDL) to characterize etiology-specific differences in neonatal cholestatic MΦ polarization. Neonatal BDL surgery was performed on female BALB/c mice at 10 days of life (DOL) with sham laparotomy as controls. Comparison was made to the Rhesus Rotavirus (RRV)-induced murine model of biliary atresia (BA). Evaluation of changes at day 7 after surgery (BDL and sham groups) and murine BA (DOL14) included laboratory data, histology (H&E, anti-CD45 and anti-CK19 staining), flow cytometry of MΦ subsets by MHCII and Ly6c expression, and single cell RNA-sequencing (scRNA-seq) analysis. Neonatal BDL achieved a 90% survival rate; mice had elevated bile acids, bilirubin, and alanine aminotransferase (ALT) versus controls (p < 0.05 for all). Histology demonstrated hepatocellular injury, CD45+ portal infiltrate, and CK19+ bile duct proliferation in neonatal BDL. Comparison to murine BA showed increased ALT in neonatal BDL despite no difference in histology Ishak score. Neonatal BDL had significantly lower MHCII-Ly6c+ MΦ versus murine BA, however, scRNA-seq identified greater etiology-specific MΦ heterogeneity with increased endocytosis in neonatal BDL MΦ versus cellular killing in murine BA MΦ. We generated an innovative murine model of neonatal obstructive cholestasis with low mortality. This model enabled comparison to murine BA to define etiology-specific cholestatic MΦ function. Further comparisons to human data may enable development of immune modulatory therapies to improve patient outcomes.

Three Distinct Transcriptional Profiles of Monocytes Associate with Disease Activity in Scleroderma Patients.

OBJECTIVE: Patients with diffuse cutaneous systemic sclerosis (dcSSc) display a complex clinical phenotype. Transcriptional profiling of whole blood or tissue from patients are affected by changes in cellular composition that drive gene expression and an inability to detect minority cell populations. We undertook this study to focus on the 2 main subtypes of circulating monocytes, classical monocytes (CMs) and nonclassical monocytes (NCMs) as a biomarker of SSc disease severity. METHODS: SSc patients were recruited from the Prospective Registry for Early Systemic Sclerosis. Clinical data were collected, as well as peripheral blood for isolation of CMs and NCMs. Age-, sex-, and race-matched healthy volunteers were recruited as controls. Bulk macrophages were isolated from the skin in a separate cohort. All samples were assayed by RNA sequencing (RNA-seq). RESULTS: We used an unbiased approach to cluster patients into 3 groups (groups A-C) based on the transcriptional signatures of CMs relative to controls. Each group maintained their characteristic transcriptional signature in NCMs. Genes up-regulated in group C demonstrated the highest expression compared to the other groups in SSc skin macrophages, relative to controls. Patients from groups B and C exhibited worse lung function than group A, although there was no difference in SSc skin disease at baseline, relative to controls. We validated our approach by applying our group classifications to published bulk monocyte RNA-seq data from SSc patients, and we found that patients without skin disease were most likely to be classified as group A. CONCLUSION: We are the first to show that transcriptional signatures of CMs and NCMs can be used to unbiasedly stratify SSc patients and correlate with disease activity outcome measures.

Tissue-resident, extravascular Ly6c- monocytes are critical for inflammation in the synovium.

Monocytes are abundant immune cells that infiltrate inflamed organs. However, the majority of monocyte studies focus on circulating cells, rather than those in tissue. Here, we identify and characterize an intravascular synovial monocyte population resembling circulating non-classical monocytes and an extravascular tissue-resident monocyte-lineage cell (TR-MC) population distinct in surface marker and transcriptional profile from circulating monocytes, dendritic cells, and tissue macrophages that are conserved in rheumatoid arthritis (RA) patients. TR-MCs are independent of NR4A1 and CCR2, long lived, and embryonically derived. TR-MCs undergo increased proliferation and reverse diapedesis dependent on LFA1 in response to arthrogenic stimuli and are required for the development of RA-like disease. Moreover, pathways that are activated in TR-MCs at the peak of arthritis overlap with those that are downregulated in LFA1-/- TR-MCs. These findings show a facet of mononuclear cell biology that could be imperative to understanding tissue-resident myeloid cell function in RA.

Stromal Transcription Factor 21 Regulates Development of the Renal Stroma via Interaction with Wnt/ß-Catenin Signaling.

Background: Kidney formation requires coordinated interactions between multiple cell types. Input from the interstitial progenitor cells is implicated in multiple aspects of kidney development. We previously reported that transcription factor 21 (Tcf21) is required for ureteric bud branching. Here, we show that Tcf21 in Foxd1+ interstitial progenitors regulates stromal formation and differentiation via interaction with ß-catenin. Methods: We utilized the Foxd1Cre;Tcf21f/f murine kidney for morphologic analysis. We used the murine clonal mesenchymal cell lines MK3/M15 to study Tcf21 interaction with Wnt/ß-catenin. Results: Absence of Tcf21 from Foxd1+ stromal progenitors caused a decrease in stromal cell proliferation, leading to marked reduction of the medullary stromal space. Lack of Tcf21 in the Foxd1+ stromal cells also led to defective differentiation of interstitial cells to smooth-muscle cells, perivascular pericytes, and mesangial cells. Foxd1Cre;Tcf21f/f kidney showed an abnormal pattern of the renal vascular tree. The stroma of Foxd1Cre;Tcf21f/f kidney demonstrated marked reduction in ß-catenin protein expression compared with wild type. Tcf21 was bound to ß-catenin both upon ß-catenin stabilization and at basal state as demonstrated by immunoprecipitation in vitro. In MK3/M15 metanephric mesenchymal cells, Tcf21 enhanced TCF/LEF promoter activity upon ß-catenin stabilization, whereas DNA-binding deficient mutated Tcf21 did not enhance TCF/LEF promoter activity. Kidney explants of Foxd1Cre;Tcf21f/f showed low mRNA expression of stromal Wnt target genes. Treatment of the explants with CHIR, a Wnt ligand mimetic, restored Wnt target gene expression. Here, we also corroborated previous evidence that normal development of the kidney stroma is required for normal development of the Six2+ nephron progenitor cells, loop of Henle, and the collecting ducts. Conclusions: These findings suggest that stromal Tcf21 facilitates medullary stroma development by enhancing Wnt/ß-catenin signaling and promotes stromal cell proliferation and differentiation. Stromal Tcf21 is also required for the development of the adjacent nephron epithelia.

Transcriptional profiling of pediatric cholestatic livers identifies three distinct macrophage populations.

BACKGROUND & AIMS: Limited understanding of the role for specific macrophage subsets in the pathogenesis of cholestatic liver injury is a barrier to advancing medical therapy. Macrophages have previously been implicated in both the mal-adaptive and protective responses in obstructive cholestasis. Recently two macrophage subsets were identified in non-diseased human liver; however, no studies to date fully define the heterogeneous macrophage subsets during the pathogenesis of cholestasis. Here, we aim to further characterize the transcriptional profile of macrophages in pediatric cholestatic liver disease. METHODS: We isolated live hepatic immune cells from patients with biliary atresia (BA), Alagille syndrome (ALGS), and non-cholestatic pediatric liver by fluorescence activated cell sorting. Through single-cell RNA sequencing analysis and immunofluorescence, we characterized cholestatic macrophages. We next compared the transcriptional profile of pediatric cholestatic and non-cholestatic macrophage populations to previously published data on normal adult hepatic macrophages. RESULTS: We identified 3 distinct macrophage populations across cholestatic liver samples and annotated them as lipid-associated macrophages, monocyte-like macrophages, and adaptive macrophages based on their transcriptional profile. Immunofluorescence of liver tissue using markers for each subset confirmed their presence across BA (n = 6) and ALGS (n = 6) patients. Cholestatic macrophages demonstrated reduced expression of immune regulatory genes as compared to normal hepatic macrophages and were distinct from macrophage populations defined in either healthy adult or pediatric non-cholestatic liver. CONCLUSIONS: We are the first to perform single-cell RNA sequencing on human pediatric cholestatic liver and identified three macrophage subsets with distinct transcriptional signatures from healthy liver macrophages. Further analyses will identify similarities and differences in these macrophage sub-populations across etiologies of cholestatic liver disease. Taken together, these findings may allow for future development of targeted therapeutic strategies to reprogram macrophages to an immune regulatory phenotype and reduce cholestatic liver injury.

A Novel Microglia-Specific Transcriptional Signature Correlates With Behavioral Deficits in Neuropsychiatric Lupus.

Neuropsychiatric symptoms of systemic lupus erythematosus (NP-SLE) affect over one-half of SLE patients, yet underlying mechanisms remain largely unknown. We demonstrate that SLE-prone mice (CReCOM) develop NP-SLE, including behavioral deficits prior to systemic autoimmunity, reduced brain volumes, decreased vascular integrity, and brain-infiltrating leukocytes. NP-SLE microglia exhibit numerical expansion, increased synaptic uptake, and a more metabolically active phenotype. Microglia from multiple SLE-prone models express a "NP-SLE signature" unrelated to type I interferon. Rather, the signature is associated with lipid metabolism, scavenger receptor activity and downregulation of inflammatory and chemotaxis processes, suggesting a more regulatory, anti-inflammatory profile. NP-SLE microglia also express genes associated with disease-associated microglia (DAM), a subset of microglia thought to be instrumental in neurodegenerative diseases. Further, expression of "NP-SLE" and "DAM" signatures correlate with the severity of behavioral deficits in young SLE-prone mice prior to overt systemic disease. Our data are the first to demonstrate the predictive value of our newly identified microglia-specific "NP-SLE" and "DAM" signatures as a surrogate for NP-SLE clinical outcomes and suggests that microglia-intrinsic defects precede contributions from systemic SLE for neuropsychiatric manifestations.

Neuropsychiatric Systemic Lupus Erythematosus Is Dependent on Sphingosine-1-Phosphate Signaling.

About 40% of patients with systemic lupus erythematosus experience diffuse neuropsychiatric manifestations, including impaired cognition and depression. Although the pathogenesis of diffuse neuropsychiatric SLE (NPSLE) is not fully understood, loss of brain barrier integrity, autoreactive antibodies, and pro-inflammatory cytokines are major contributors to disease development. Fingolimod, a sphingosine-1-phosphate (S1P) receptor modulator, prevents lymphocyte egress from lymphoid organs through functional antagonism of S1P receptors. In addition to reducing the circulation of autoreactive lymphocytes, fingolimod has direct neuroprotective effects such as preserving brain barrier integrity and decreasing pro-inflammatory cytokine secretion by astrocytes and microglia. Given these effects, we hypothesized that fingolimod would attenuate neurobehavioral deficits in MRL-lpr/lpr (MRL/lpr) mice, a validated neuropsychiatric lupus model. Fingolimod treatment was initiated after the onset of disease, and mice were assessed for alterations in cognitive function and emotionality. We found that fingolimod significantly attenuated spatial memory deficits and depression-like behavior in MRL/lpr mice. Immunofluorescent staining demonstrated a dramatic lessening of brain T cell and macrophage infiltration, and a significant reduction in cortical leakage of serum albumin, in fingolimod treated mice. Astrocytes and endothelial cells from treated mice exhibited reduced expression of inflammatory genes, while microglia showed differential regulation of key immune pathways. Notably, cytokine levels within the cortex and hippocampus were not appreciably decreased with fingolimod despite the improved neurobehavioral profile. Furthermore, despite a reduction in splenomegaly, lymphadenopathy, and circulating autoantibody titers, IgG deposition within the brain was unaffected by treatment. These findings suggest that fingolimod mediates attenuation of NPSLE through a mechanism that is not dependent on reduction of autoantibodies or cytokines, and highlight modulation of the S1P signaling pathway as a novel therapeutic target in lupus involving the central nervous system.