Progranulin Deficiency Promotes Circuit-Specific Synaptic Pruning by Microglia via Complement Activation.

Microglia maintain homeostasis in the brain, but whether aberrant microglial activation can cause neurodegeneration remains controversial. Here, we use transcriptome profiling to demonstrate that deficiency in frontotemporal dementia (FTD) gene progranulin (Grn) leads to an age-dependent, progressive upregulation of lysosomal and innate immunity genes, increased complement production, and enhanced synaptic pruning in microglia. During aging, Grn(-/-) mice show profound microglia infiltration and preferential elimination of inhibitory synapses in the ventral thalamus, which lead to hyperexcitability in the thalamocortical circuits and obsessive-compulsive disorder (OCD)-like grooming behaviors. Remarkably, deleting C1qa gene significantly reduces synaptic pruning by Grn(-/-) microglia and mitigates neurodegeneration, behavioral phenotypes, and premature mortality in Grn(-/-) mice. Together, our results uncover a previously unrecognized role of progranulin in suppressing aberrant microglia activation during aging. These results represent an important conceptual advance that complement activation and microglia-mediated synaptic pruning are major drivers, rather than consequences, of neurodegeneration caused by progranulin deficiency.

Neutrophils prime a long-lived effector macrophage phenotype that mediates accelerated helminth expulsion.

We examined the role of innate cells in acquired resistance to the natural murine parasitic nematode, Nippostrongylus brasiliensis. Macrophages obtained from lungs as late as 45 d after N. brasiliensis inoculation were able to transfer accelerated parasite clearance to naive recipients. Primed macrophages adhered to larvae in vitro and triggered increased mortality of parasites. Neutrophil depletion in primed mice abrogated the protective effects of transferred macrophages and inhibited their in vitro binding to larvae. Neutrophils in parasite-infected mice showed a distinct transcriptional profile and promoted alternatively activated M2 macrophage polarization through secretory factors including IL-13. Differentially activated neutrophils in the context of a type 2 immune response therefore prime a long-lived effector macrophage phenotype that directly mediates rapid nematode damage and clearance.

An End-to-End Platform for Digital Pathology Using Hyperspectral Autofluorescence Microscopy and Deep Learning-Based Virtual Histology.

Conventional histopathology involves expensive and labor-intensive processes that often consume tissue samples, rendering them unavailable for other analyses. We present a novel end-to-end workflow for pathology powered by hyperspectral microscopy and deep learning. First, we developed a custom hyperspectral microscope to nondestructively image the autofluorescence of unstained tissue sections. We then trained a deep learning model to use autofluorescence to generate virtual histologic stains, which avoids the cost and variability of chemical staining procedures and conserves tissue samples. We showed that the virtual images reproduce the histologic features present in the real-stained images using a randomized nonalcoholic steatohepatitis (NASH) scoring comparison study, where both real and virtual stains are scored by pathologists (D.T., A.D.B., R.K.P.). The test showed moderate-to-good concordance between pathologists' scoring on corresponding real and virtual stains. Finally, we developed deep learning-based models for automated NASH Clinical Research Network score prediction. We showed that the end-to-end automated pathology platform is comparable with an independent panel of pathologists for NASH Clinical Research Network scoring when evaluated against the expert pathologist consensus scores. This study provides proof of concept for this virtual staining strategy, which could improve cost, efficiency, and reliability in pathology and enable novel approaches to spatial biology research.

Molecular definition of the identity and activation of natural killer cells.

Using whole-genome microarray data sets of the Immunological Genome Project, we demonstrate a closer transcriptional relationship between NK cells and T cells than between any other leukocytes, distinguished by their shared expression of genes encoding molecules with similar signaling functions. Whereas resting NK cells are known to share expression of a few genes with cytotoxic CD8(+) T cells, our transcriptome-wide analysis demonstrates that the commonalities extend to hundreds of genes, many encoding molecules with unknown functions. Resting NK cells demonstrate a 'preprimed' state compared with naive T cells, which allows NK cells to respond more rapidly to viral infection. Collectively, our data provide a global context for known and previously unknown molecular aspects of NK cell identity and function by delineating the genome-wide repertoire of gene expression of NK cells in various states.

Multi-omic profiling reveals early immunological indicators for identifying COVID-19 Progressors.

We sought to better understand the immune response during the immediate post-diagnosis phase of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) by identifying molecular associations with longitudinal disease outcomes. Multi-omic analyses identified differences in immune cell composition, cytokine levels, and cell subset-specific transcriptomic and epigenomic signatures between individuals on a more serious disease trajectory (Progressors) as compared to those on a milder course (Non-progressors). Higher levels of multiple cytokines were observed in Progressors, with IL-6 showing the largest difference. Blood monocyte cell subsets were also skewed, showing a comparative decrease in non-classical CD14-CD16+ and intermediate CD14+CD16+ monocytes. In lymphocytes, the CD8+ T effector memory cells displayed a gene expression signature consistent with stronger T cell activation in Progressors. These early stage observations could serve as the basis for the development of prognostic biomarkers of disease risk and interventional strategies to improve the management of severe COVID-19. BACKGROUND: Much of the literature on immune response post-SARS-CoV-2 infection has been in the acute and post-acute phases of infection. TRANSLATIONAL SIGNIFICANCE: We found differences at early time points of infection in approximately 160 participants. We compared multi-omic signatures in immune cells between individuals progressing to needing more significant medical intervention and non-progressors. We observed widespread evidence of a state of increased inflammation associated with progression, supported by a range of epigenomic, transcriptomic, and proteomic signatures. The signatures we identified support other findings at later time points and serve as the basis for prognostic biomarker development or to inform interventional strategies.

The SLC36 transporter Pathetic is required for extreme dendrite growth in Drosophila sensory neurons.

Dendrites exhibit enormous diversity in form and can differ in size by several orders of magnitude even in a single animal. However, whether neurons with large dendrite arbors have specialized mechanisms to support their growth demands is unknown. To address this question, we conducted a genetic screen for mutations that differentially affected growth in neurons with different-sized dendrite arbors. From this screen, we identified a mutant that selectively affects dendrite growth in neurons with large dendrite arbors without affecting dendrite growth in neurons with small dendrite arbors or the animal overall. This mutant disrupts a putative amino acid transporter, Pathetic (Path), that localizes to the cell surface and endolysosomal compartments in neurons. Although Path is broadly expressed in neurons and nonneuronal cells, mutation of path impinges on nutrient responses and protein homeostasis specifically in neurons with large dendrite arbors but not in other cells. Altogether, our results demonstrate that specialized molecular mechanisms exist to support growth demands in neurons with large dendrite arbors and define Path as a founding member of this growth program.

A20 upregulation during treated HIV disease is associated with intestinal epithelial cell recovery and function.

TRIAL REGISTRATION: ClinicalTrials.gov Clinical Trial NCT00594880.

Correction: Macrophage Colony Stimulating Factor Derived from CD4+ T Cells Contributes to Control of a Blood-Borne Infection.

[This corrects the article DOI: 10.1371/journal.ppat.1006046.].

Correction: Cathelicidin Insufficiency in Patients with Fatal Leptospirosis.

[This corrects the article DOI: 10.1371/journal.ppat.1005943.].

Empirical assessment of the impact of sample number and read depth on RNA-Seq analysis workflow performance.

BACKGROUND: RNA-Sequencing analysis methods are rapidly evolving, and the tool choice for each step of one common workflow, differential expression analysis, which includes read alignment, expression modeling, and differentially expressed gene identification, has a dramatic impact on performance characteristics. Although a number of workflows are emerging as high performers that are robust to diverse input types, the relative performance characteristics of these workflows when either read depth or sample number is limited-a common occurrence in real-world practice-remain unexplored. RESULTS: Here, we evaluate the impact of varying read depth and sample number on the performance of differential gene expression identification workflows, as measured by precision, or the fraction of genes correctly identified as differentially expressed, and by recall, or the fraction of differentially expressed genes identified. We focus our analysis on 30 high-performing workflows, systematically varying the read depth and number of biological replicates of patient monocyte samples provided as input. We find that, in general for most workflows, read depth has little effect on workflow performance when held above two million reads per sample, with reduced workflow performance below this threshold. The greatest impact of decreased sample number is seen below seven samples per group, when more heterogeneity in workflow performance is observed. The choice of differential expression identification tool, in particular, has a large impact on the response to limited inputs. CONCLUSIONS: Among the tested workflows, the recall/precision balance remains relatively stable at a range of read depths and sample numbers, although some workflows are more sensitive to input restriction. At ranges typically recommended for biological studies, performance is more greatly impacted by the number of biological replicates than by read depth. Caution should be used when selecting analysis workflows and interpreting results from low sample number experiments, as all workflows exhibit poorer performance at lower sample numbers near typically reported values, with variable impact on recall versus precision. These analyses highlight the performance characteristics of common differential gene expression workflows at varying read depths and sample numbers, and provide empirical guidance in experimental and analytical design.

Macrophage Colony Stimulating Factor Derived from CD4+ T Cells Contributes to Control of a Blood-Borne Infection.

Dynamic regulation of leukocyte population size and activation state is crucial for an effective immune response. In malaria, Plasmodium parasites elicit robust host expansion of macrophages and monocytes, but the underlying mechanisms remain unclear. Here we show that myeloid expansion during P. chabaudi infection is dependent upon both CD4+ T cells and the cytokine Macrophage Colony Stimulating Factor (MCSF). Single-cell RNA-Seq analysis on antigen-experienced T cells revealed robust expression of Csf1, the gene encoding MCSF, in a sub-population of CD4+ T cells with distinct transcriptional and surface phenotypes. Selective deletion of Csf1 in CD4+ cells during P. chabaudi infection diminished proliferation and activation of certain myeloid subsets, most notably lymph node-resident CD169+ macrophages, and resulted in increased parasite burden and impaired recovery of infected mice. Depletion of CD169+ macrophages during infection also led to increased parasitemia and significant host mortality, confirming a previously unappreciated role for these cells in control of P. chabaudi. This work establishes the CD4+ T cell as a physiologically relevant source of MCSF in vivo; probes the complexity of the CD4+ T cell response during type 1 infection; and delineates a novel mechanism by which T helper cells regulate myeloid cells to limit growth of a blood-borne intracellular pathogen.

Cathelicidin Insufficiency in Patients with Fatal Leptospirosis.

Leptospirosis causes significant morbidity and mortality worldwide; however, the role of the host immune response in disease progression and high case fatality (>10-50%) is poorly understood. We conducted a multi-parameter investigation of patients with acute leptospirosis to identify mechanisms associated with case fatality. Whole blood transcriptional profiling of 16 hospitalized Brazilian patients with acute leptospirosis (13 survivors, 3 deceased) revealed fatal cases had lower expression of the antimicrobial peptide, cathelicidin, and chemokines, but more abundant pro-inflammatory cytokine receptors. In contrast, survivors generated strong adaptive immune signatures, including transcripts relevant to antigen presentation and immunoglobulin production. In an independent cohort (23 survivors, 22 deceased), fatal cases had higher bacterial loads (P = 0.0004) and lower anti-Leptospira antibody titers (P = 0.02) at the time of hospitalization, independent of the duration of illness. Low serum cathelicidin and RANTES levels during acute illness were independent risk factors for higher bacterial loads (P = 0.005) and death (P = 0.04), respectively. To investigate the mechanism of cathelicidin in patients surviving acute disease, we administered LL-37, the active peptide of cathelicidin, in a hamster model of lethal leptospirosis and found it significantly decreased bacterial loads and increased survival. Our findings indicate that the host immune response plays a central role in severe leptospirosis disease progression. While drawn from a limited study size, significant conclusions include that poor clinical outcomes are associated with high systemic bacterial loads, and a decreased antibody response. Furthermore, our data identified a key role for the antimicrobial peptide, cathelicidin, in mounting an effective bactericidal response against the pathogen, which represents a valuable new therapeutic approach for leptospirosis.

Empirical assessment of analysis workflows for differential expression analysis of human samples using RNA-Seq.

BACKGROUND: RNA-Seq has supplanted microarrays as the preferred method of transcriptome-wide identification of differentially expressed genes. However, RNA-Seq analysis is still rapidly evolving, with a large number of tools available for each of the three major processing steps: read alignment, expression modeling, and identification of differentially expressed genes. Although some studies have benchmarked these tools against gold standard gene expression sets, few have evaluated their performance in concert with one another. Additionally, there is a general lack of testing of such tools on real-world, physiologically relevant datasets, which often possess qualities not reflected in tightly controlled reference RNA samples or synthetic datasets. RESULTS: Here, we evaluate 219 combinatorial implementations of the most commonly used analysis tools for their impact on differential gene expression analysis by RNA-Seq. A test dataset was generated using highly purified human classical and nonclassical monocyte subsets from a clinical cohort, allowing us to evaluate the performance of 495 unique workflows, when accounting for differences in expression units and gene- versus transcript-level estimation. We find that the choice of methodologies leads to wide variation in the number of genes called significant, as well as in performance as gauged by precision and recall, calculated by comparing our RNA-Seq results to those from four previously published microarray and BeadChip analyses of the same cell populations. The method of differential gene expression identification exhibited the strongest impact on performance, with smaller impacts from the choice of read aligner and expression modeler. Many workflows were found to exhibit similar overall performance, but with differences in their calibration, with some biased toward higher precision and others toward higher recall. CONCLUSIONS: There is significant heterogeneity in the performance of RNA-Seq workflows to identify differentially expressed genes. Among the higher performing workflows, different workflows exhibit a precision/recall tradeoff, and the ultimate choice of workflow should take into consideration how the results will be used in subsequent applications. Our analyses highlight the performance characteristics of these workflows, and the data generated in this study could also serve as a useful resource for future development of software for RNA-Seq analysis.

The microRNA bantam regulates a developmental transition in epithelial cells that restricts sensory dendrite growth.

As animals grow, many early born structures grow by cell expansion rather than cell addition; thus growth of distinct structures must be coordinated to maintain proportionality. This phenomenon is particularly widespread in the nervous system, with dendrite arbors of many neurons expanding in concert with their substrate to sustain connectivity and maintain receptive field coverage as animals grow. After rapidly growing to establish body wall coverage, dendrites of Drosophila class IV dendrite arborization (C4da) neurons grow synchronously with their substrate, the body wall epithelium, providing a system to study how proportionality is maintained during animal growth. Here, we show that the microRNA bantam (ban) ensures coordinated growth of C4da dendrites and the epithelium through regulation of epithelial endoreplication, a modified cell cycle that entails genome amplification without cell division. In Drosophila larvae, epithelial endoreplication leads to progressive changes in dendrite-extracellular matrix (ECM) and dendrite-epithelium contacts, coupling dendrite/substrate expansion and restricting dendrite growth beyond established boundaries. Moreover, changes in epithelial expression of cell adhesion molecules, including the beta-integrin myospheroid (mys), accompany this developmental transition. Finally, endoreplication and the accompanying changes in epithelial mys expression are required to constrain late-stage dendrite growth and structural plasticity. Hence, modulating epithelium-ECM attachment probably influences substrate permissivity for dendrite growth and contributes to the dendrite-substrate coupling that ensures proportional expansion of the two cell types.

FCRL5 Delineates Functionally Impaired Memory B Cells Associated with Plasmodium falciparum Exposure.

Exposure to Plasmodium falciparum is associated with circulating "atypical" memory B cells (atMBCs), which appear similar to dysfunctional B cells found in HIV-infected individuals. Functional analysis of atMBCs has been limited, with one report suggesting these cells are not dysfunctional but produce protective antibodies. To better understand the function of malaria-associated atMBCs, we performed global transcriptome analysis of these cells, obtained from individuals living in an area of high malaria endemicity in Uganda. Comparison of gene expression data suggested down-modulation of B cell receptor signaling and apoptosis in atMBCs compared to classical MBCs. Additionally, in contrast to previous reports, we found upregulation of Fc receptor-like 5 (FCRL5), but not FCRL4, on atMBCs. Atypical MBCs were poor spontaneous producers of antibody ex vivo, and higher surface expression of FCRL5 defined a distinct subset of atMBCs compromised in its ability to produce antibody upon stimulation. Moreover, higher levels of P. falciparum exposure were associated with increased frequencies of FCRL5+ atMBCs. Together, our findings suggest that FCLR5+ identifies a functionally distinct, and perhaps dysfunctional, subset of MBCs in individuals exposed to P. falciparum.

JUNB is a key transcriptional modulator of macrophage activation.

Activated macrophages are crucial for restriction of microbial infection but may also promote inflammatory pathology in a wide range of both infectious and sterile conditions. The pathways that regulate macrophage activation are therefore of great interest. Recent studies in silico have putatively identified key transcription factors that may control macrophage activation, but experimental validation is lacking. In this study, we generated a macrophage regulatory network from publicly available microarray data, employing steps to enrich for physiologically relevant interactions. Our analysis predicted a novel relationship between the AP-1 family transcription factor Junb and the gene Il1b, encoding the pyrogen IL-1ß, which macrophages express upon activation by inflammatory stimuli. Previously, Junb has been characterized primarily as a negative regulator of the cell cycle, whereas AP-1 activity in myeloid inflammatory responses has largely been attributed to c-Jun. We confirmed experimentally that Junb is required for full expression of Il1b, and of additional genes involved in classical inflammation, in macrophages treated with LPS and other immunostimulatory molecules. Furthermore, Junb modulates expression of canonical markers of alternative activation in macrophages treated with IL-4. Our results demonstrate that JUNB is a significant modulator of both classical and alternative macrophage activation. Further, this finding provides experimental validation for our network modeling approach, which will facilitate the future use of gene expression data from open databases to reveal novel, physiologically relevant regulatory relationships.

(+)-SJ733, a clinical candidate for malaria that acts through ATP4 to induce rapid host-mediated clearance of Plasmodium.

Drug discovery for malaria has been transformed in the last 5 years by the discovery of many new lead compounds identified by phenotypic screening. The process of developing these compounds as drug leads and studying the cellular responses they induce is revealing new targets that regulate key processes in the Plasmodium parasites that cause malaria. We disclose herein that the clinical candidate (+)-SJ733 acts upon one of these targets, ATP4. ATP4 is thought to be a cation-transporting ATPase responsible for maintaining low intracellular Na(+) levels in the parasite. Treatment of parasitized erythrocytes with (+)-SJ733 in vitro caused a rapid perturbation of Na(+) homeostasis in the parasite. This perturbation was followed by profound physical changes in the infected cells, including increased membrane rigidity and externalization of phosphatidylserine, consistent with eryptosis (erythrocyte suicide) or senescence. These changes are proposed to underpin the rapid (+)-SJ733-induced clearance of parasites seen in vivo. Plasmodium falciparum ATPase 4 (pfatp4) mutations that confer resistance to (+)-SJ733 carry a high fitness cost. The speed with which (+)-SJ733 kills parasites and the high fitness cost associated with resistance-conferring mutations appear to slow and suppress the selection of highly drug-resistant mutants in vivo. Together, our data suggest that inhibitors of PfATP4 have highly attractive features for fast-acting antimalarials to be used in the global eradication campaign.

Trimming of sequence reads alters RNA-Seq gene expression estimates.

BACKGROUND: High-throughput RNA-Sequencing (RNA-Seq) has become the preferred technique for studying gene expression differences between biological samples and for discovering novel isoforms, though the techniques to analyze the resulting data are still immature. One pre-processing step that is widely but heterogeneously applied is trimming, in which low quality bases, identified by the probability that they are called incorrectly, are removed. However, the impact of trimming on subsequent alignment to a genome could influence downstream analyses including gene expression estimation; we hypothesized that this might occur in an inconsistent manner across different genes, resulting in differential bias. RESULTS: To assess the effects of trimming on gene expression, we generated RNA-Seq data sets from four samples of larval Drosophila melanogaster sensory neurons, and used three trimming algorithms--SolexaQA, Trimmomatic, and ConDeTri-to perform quality-based trimming across a wide range of stringencies. After aligning the reads to the D. melanogaster genome with TopHat2, we used Cuffdiff2 to compare the original, untrimmed gene expression estimates to those following trimming. With the most aggressive trimming parameters, over ten percent of genes had significant changes in their estimated expression levels. This trend was seen with two additional RNA-Seq data sets and with alternative differential expression analysis pipelines. We found that the majority of the expression changes could be mitigated by imposing a minimum length filter following trimming, suggesting that the differential gene expression was primarily being driven by spurious mapping of short reads. Slight differences with the untrimmed data set remained after length filtering, which were associated with genes with low exon numbers and high GC content. Finally, an analysis of paired RNA-seq/microarray data sets suggests that no or modest trimming results in the most biologically accurate gene expression estimates. CONCLUSIONS: We find that aggressive quality-based trimming has a large impact on the apparent makeup of RNA-Seq-based gene expression estimates, and that short reads can have a particularly strong impact. We conclude that implementation of trimming in RNA-Seq analysis workflows warrants caution, and if used, should be used in conjunction with a minimum read length filter to minimize the introduction of unpredictable changes in expression estimates.

Brain trauma elicits non-canonical macrophage activation states.

BACKGROUND: Macrophage polarization programs, commonly referred to as "classical" and "alternative" activation, are widely considered as distinct states that are exclusive of one another and are associated with different functions such as inflammation and wound healing, respectively. In a number of disease contexts, such as traumatic brain injury (TBI), macrophage polarization influences the extent of pathogenesis, and efforts are underway to eliminate pathogenic subsets. However, previous studies have not distinguished whether the simultaneous presence of both classical and alternative activation signatures represents the admixture of differentially polarized macrophages or if they have adopted a unique state characterized by components of both classical and alternative activation. METHODS: We analyzed the gene expression profiles of individual monocyte-derived brain macrophages responding to TBI using single-cell RNA sequencing. RNA flow cytometry was used as another single-cell analysis technique to validate the single-cell RNA sequencing results. RESULTS: The analysis of signature polarization genes by single-cell RNA sequencing revealed the presence of diverse activation states, including M(IL4), M(IL10), and M(LPS, IFNγ). However, the expression of a given polarization marker was no more likely than at random to predict simultaneous expression or repression of markers of another polarization program within the same cell, suggesting a lack of exclusivity in macrophage polarization states in vivo in TBI. Also unexpectedly, individual TBI macrophages simultaneously expressed high levels of signature polarization genes across two or three different polarization states and in several distinct and seemingly incompatible combinations. CONCLUSIONS: Single-cell gene expression profiling demonstrated that monocytic macrophages in TBI are not comprised of distinctly polarized subsets but are uniquely and broadly activated. TBI macrophage activation in vivo is deeply complex, with individual cells concurrently adopting both inflammatory and reparative features with a lack of exclusivity. These data provide physiologically relevant evidence that the early macrophage response to TBI is comprised of novel activation states that are discordant with the current paradigm of macrophage polarization-a key consideration for therapeutic modulation.

Distinct functional programming of human fetal and adult monocytes.

Preterm birth affects 1 out of 9 infants in the United States and is the leading cause of long-term neurologic handicap and infant mortality, accounting for 35% of all infant deaths in 2008. Although cytokines including interferon-γ (IFN-γ), interleukin-10 (IL-10), IL-6, and IL-1 are produced in response to in utero infection and are strongly associated with preterm labor, little is known about how human fetal immune cells respond to these cytokines. We demonstrate that fetal and adult CD14(+)CD16(-) classical monocytes are distinct in terms of basal transcriptional profiles and in phosphorylation of signal transducers and activators of transcription (STATs) in response to cytokines. Fetal monocytes phosphorylate canonical and noncanonical STATs and respond more strongly to IFN-γ, IL-6, and IL-4 than adult monocytes. We demonstrate a higher ratio of SOCS3 to IL-6 receptor in adult monocytes than in fetal monocytes, potentially explaining differences in STAT phosphorylation. Additionally, IFN-γ signaling results in upregulation of antigen presentation and costimulatory machinery in adult, but not fetal, monocytes. These findings represent the first evidence that primary human fetal and adult monocytes are functionally distinct, potentially explaining how these cells respond differentially to cytokines implicated in development, in utero infections, and the pathogenesis of preterm labor.