Deterministic reprogramming of neutrophils within tumors.

Neutrophils are increasingly recognized as key players in the tumor immune response and are associated with poor clinical outcomes. Despite recent advances characterizing the diversity of neutrophil states in cancer, common trajectories and mechanisms governing the ontogeny and relationship between these neutrophil states remain undefined. Here, we demonstrate that immature and mature neutrophils that enter tumors undergo irreversible epigenetic, transcriptional, and proteomic modifications to converge into a distinct, terminally differentiated dcTRAIL-R1+ state. Reprogrammed dcTRAIL-R1+ neutrophils predominantly localize to a glycolytic and hypoxic niche at the tumor core and exert pro-angiogenic function that favors tumor growth. We found similar trajectories in neutrophils across multiple tumor types and in humans, suggesting that targeting this program may provide a means of enhancing certain cancer immunotherapies.

Harnessing Single-Cell RNA Sequencing to Identify Dendritic Cell Types, Characterize Their Biological States, and Infer Their Activation Trajectory.

Dendritic cells (DCs) orchestrate innate and adaptive immunity, by translating the sensing of distinct danger signals into the induction of different effector lymphocyte responses, to induce the defense mechanisms the best suited to face the threat. Hence, DCs are very plastic, which results from two key characteristics. First, DCs encompass distinct cell types specialized in different functions. Second, each DC type can undergo different activation states, fine-tuning its functions depending on its tissue microenvironment and the pathophysiological context, by adapting the output signals it delivers to the input signals it receives. Hence, to better understand DC biology and harness it in the clinic, we must determine which combinations of DC types and activation states mediate which functions and how.To decipher the nature, functions, and regulation of DC types and their physiological activation states, one of the methods that can be harnessed most successfully is ex vivo single-cell RNA sequencing (scRNAseq). However, for new users of this approach, determining which analytics strategy and computational tools to choose can be quite challenging, considering the rapid evolution and broad burgeoning in the field. In addition, awareness must be raised on the need for specific, robust, and tractable strategies to annotate cells for cell type identity and activation states. It is also important to emphasize the necessity of examining whether similar cell activation trajectories are inferred by using different, complementary methods. In this chapter, we take these issues into account for providing a pipeline for scRNAseq analysis and illustrating it with a tutorial reanalyzing a public dataset of mononuclear phagocytes isolated from the lungs of naïve or tumor-bearing mice. We describe this pipeline step-by-step, including data quality controls, dimensionality reduction, cell clustering, cell cluster annotation, inference of the cell activation trajectories, and investigation of the underpinning molecular regulation. It is accompanied with a more complete tutorial on GitHub. We hope that this method will be helpful for both wet lab and bioinformatics researchers interested in harnessing scRNAseq data for deciphering the biology of DCs or other cell types and that it will contribute to establishing high standards in the field.

Conditional Knockout of Hypoxia-Inducible Factor 1-Alpha in Tumor-Infiltrating Neutrophils Protects against Pancreatic Ductal Adenocarcinoma.

Large numbers of neutrophils infiltrate tumors and comprise a notable component of the inflammatory tumor microenvironment. While it is established that tumor cells exhibit the Warburg effect for energy production, the contribution of the neutrophil metabolic state to tumorigenesis is unknown. Here, we investigated whether neutrophil infiltration and metabolic status promotes tumor progression in an orthotopic mouse model of pancreatic ductal adenocarcinoma (PDAC). We observed a large increase in the proportion of neutrophils in the blood and tumor upon orthotopic transplantation. Intriguingly, these tumor-infiltrating neutrophils up-regulated glycolytic factors and hypoxia-inducible factor 1-alpha (HIF-1α) expression compared to neutrophils from the bone marrow and blood of the same mouse. This enhanced glycolytic signature was also observed in human PDAC tissue samples. Strikingly, neutrophil-specific deletion of HIF-1α (HIF-1αΔNφ) significantly reduced tumor burden and improved overall survival in orthotopic transplanted mice, by converting the pro-tumorigenic neutrophil phenotype to an anti-tumorigenic phenotype. This outcome was associated with elevated reactive oxygen species production and activated natural killer cells and CD8+ cytotoxic T cells compared to littermate control mice. These data suggest a role for HIF-1α in neutrophil metabolism, which could be exploited as a target for metabolic modulation in cancer.

Transcriptomics of rhinovirus persistence reveals sustained expression of RIG-I and interferon-stimulated genes in nasal epithelial cells in vitro.

BACKGROUND: Human rhinoviruses (HRVs) are frequently associated with asthma exacerbations, and have been found in the airways of asthmatic patients. While HRV-induced acute infection is well-documented, it is less clear whether the nasal epithelium sustains prolonged HRV infections along with the associated activation of host immune responses. OBJECTIVE: To investigate sustainably regulated host responses of human nasal epithelial cells (hNECs) during HRV persistence. METHODS: Using a time-course study, HRV16 persistence and viral replication dynamics were established using an in vitro infection model of hNECs. RNA sequencing was performed on hNECs in the early and late stages of infection at 3 and 14 days post-infection (dpi), respectively. The functional enrichment of differentially expressed genes (DEGs) was evaluated using gene ontology (GO) and Ingenuity pathway analysis. RESULTS: HRV RNA and protein expression persisted throughout prolonged infections, even after decreased production of infectious virus progeny. GO analysis of unique DEGs indicated altered regulation of pathways related to ciliary function and airway remodeling at 3 dpi and serine-type endopeptidase activity at 14 dpi. The functional enrichment of shared DEGs between the two time-points was related to interferon (IFN) and cytoplasmic pattern recognition receptor (PRR) signaling pathways. Validation of the sustained regulation of candidate genes confirmed the persistent expression of RIG-I and revealed its close co-regulation with interferon-stimulated genes (ISGs) during HRV persistence. CONCLUSIONS: The persistence of HRV RNA does not necessarily indicate an active infection during prolonged infection. The sustained expression of RIG-I and ISGs in response to viral RNA persistence highlights the importance of assessing how immune-activating host factors can change during active HRV infection and the immune regulation that persists thereafter.

Circulating CD1c+ myeloid dendritic cells are potential precursors to LCH lesion CD1a+CD207+ cells.

Langerhans cell histiocytosis (LCH) is a myeloproliferative disorder that is characterized by the inflammatory lesions with pathogenic CD1a+CD207+ dendritic cells (DCs). BRAFV600E and other somatic activating MAPK gene mutations have been identified in differentiating bone marrow and blood myeloid cells, but the origin of the LCH lesion CD1a+CD207+ DCs and mechanisms of lesion formation remain incompletely defined. To identify candidate LCH CD1a+CD207+ DC precursor populations, gene-expression profiles of LCH lesion CD1a+CD207+ DCs were first compared with established gene signatures from human myeloid cell subpopulations. Interestingly, the CD1c+ myeloid DC (mDC) gene signature was most enriched in the LCH CD1a+CD207+ DC transcriptome. Additionally, the BRAFV600E allele was not only localized to CD1a+CD207- DCs and CD1a+CD207+ DCs, but it was also identified in CD1c+ mDCs in LCH lesions. Transcriptomes of CD1a+CD207- DCs were nearly indistinguishable from CD1a+CD207+ DCs (both CD1a+CD207low and CD1a+CD207high subpopulations). Transcription profiles of LCH lesion CD1a+CD207+ DCs and peripheral blood CD1c+ mDCs from healthy donors were compared to identify potential LCH DC-specific biomarkers: HLA-DQB2 expression was significantly increased in LCH lesion CD1a+CD207+ DCs compared with circulating CD1c+ mDCs from healthy donors. HLA-DQB2 antigen was identified on LCH lesion CD1a+CD207- DCs and CD1a+CD207+ DCs as well as on CD1c+(CD1a+CD207-) mDCs, but it was not identified in any other lesion myeloid subpopulations. HLA-DQB2 expression was specific to peripheral blood of patients with BRAFV600E+ peripheral blood mononuclear cells, and HLA-DQB2+CD1c+ blood cells were highly enriched for the BRAFV600E in these patients. These data support a model in which blood CD1c+HLA-DQB2+ mDCs with activated ERK migrate to lesion sites where they differentiate into pathogenic CD1a+CD207+ DCs.

Cell-autonomous FLT3L shedding via ADAM10 mediates conventional dendritic cell development in mouse spleen.

Conventional dendritic cells (cDCs) derive from bone marrow (BM) precursors that undergo cascades of developmental programs to terminally differentiate in peripheral tissues. Pre-cDC1s and pre-cDC2s commit in the BM to each differentiate into CD8α+/CD103+ cDC1s and CD11b+ cDC2s, respectively. Although both cDCs rely on the cytokine FLT3L during development, mechanisms that ensure cDC accessibility to FLT3L have yet to be elucidated. Here, we generated mice that lacked a disintegrin and metalloproteinase (ADAM) 10 in DCs (Itgax-cre × Adam10-fl/fl; ADAM10∆DC) and found that ADAM10 deletion markedly impacted splenic cDC2 development. Pre-cDC2s accumulated in the spleen with transcriptomic alterations that reflected their inability to differentiate and exhibited abrupt failure to survive as terminally differentiated cDC2s. Induced ADAM10 ablation also led to the reduction of terminally differentiated cDC2s, and restoration of Notch signaling, a major pathway downstream of ADAM10, only modestly rescued them. ADAM10∆DC BM failed to generate cDC2s in BM chimeric mice with or without cotransferred ADAM10-sufficient BM, indicating that cDC2 development required cell-autonomous ADAM10. We determined cDC2s to be sources of soluble FLT3L, as supported by decreased serum FLT3L concentration and the retention of membrane-bound FLT3L on cDC2 surfaces in ADAM10∆DC mice, and by demonstrating the release of soluble FLT3L by cDC2 in ex vivo culture supernatants. Through in vitro studies utilizing murine embryonic fibroblasts, we determined FLT3L to be a substrate for ADAM10. These data collectively reveal cDC2s as FLT3L sources and highlight a cell-autonomous mechanism that may enhance FLT3L accessibility for cDC2 development and survival.

Two distinct interstitial macrophage populations coexist across tissues in specific subtissular niches.

Macrophages are a heterogeneous cell population involved in tissue homeostasis, inflammation, and various pathologies. Although the major tissue-resident macrophage populations have been extensively studied, interstitial macrophages (IMs) residing within the tissue parenchyma remain poorly defined. Here we studied IMs from murine lung, fat, heart, and dermis. We identified two independent IM subpopulations that are conserved across tissues: Lyve1loMHCIIhiCX3CR1hi (Lyve1loMHCIIhi) and Lyve1hiMHCIIloCX3CR1lo (Lyve1hiMHCIIlo) monocyte-derived IMs, with distinct gene expression profiles, phenotypes, functions, and localizations. Using a new mouse model of inducible macrophage depletion (Slco2b1 flox/DTR), we found that the absence of Lyve1hiMHCIIlo IMs exacerbated experimental lung fibrosis. Thus, we demonstrate that two independent populations of IMs coexist across tissues and exhibit conserved niche-dependent functional programming.

Organ-Specific Fate, Recruitment, and Refilling Dynamics of Tissue-Resident Macrophages during Blood-Stage Malaria.

Inflammation-induced disappearance of tissue-resident macrophages represents a key pathogen defense mechanism. Using a model of systemic blood-stage malaria, we studied the dynamics of tissue-resident macrophages in multiple organs to determine how they are depleted and refilled during the course of disease. We show that Plasmodium infection results in a transient loss of embryonically established resident macrophages prior to the parasitemia peak. Fate-mapping analysis reveals that inflammatory monocytes contribute to the repopulation of the emptied niches of splenic red pulp macrophages and hepatic Kupffer cells, while lung alveolar macrophages refill their niche predominantly through self-renewal. Interestingly, the local microenvironment of the spleen and liver can "imprint" the molecular characteristics of fetal-derived macrophages on newly differentiated bone marrow-derived immigrants with remarkably similar gene expression profiles and turnover kinetics. Thus, the mononuclear phagocytic system has developed distinct but effective tissue-specific strategies to replenish emptied niches to guarantee the functional integrity of the system.

Bystander CD8+ T cells are abundant and phenotypically distinct in human tumour infiltrates.

Various forms of immunotherapy, such as checkpoint blockade immunotherapy, are proving to be effective at restoring T cell-mediated immune responses that can lead to marked and sustained clinical responses, but only in some patients and cancer types1-4. Patients and tumours may respond unpredictably to immunotherapy partly owing to heterogeneity of the immune composition and phenotypic profiles of tumour-infiltrating lymphocytes (TILs) within individual tumours and between patients5,6. Although there is evidence that tumour-mutation-derived neoantigen-specific T cells play a role in tumour control2,4,7-10, in most cases the antigen specificities of phenotypically diverse tumour-infiltrating T cells are largely unknown. Here we show that human lung and colorectal cancer CD8+ TILs can not only be specific for tumour antigens (for example, neoantigens), but also recognize a wide range of epitopes unrelated to cancer (such as those from Epstein-Barr virus, human cytomegalovirus or influenza virus). We found that these bystander CD8+ TILs have diverse phenotypes that overlap with tumour-specific cells, but lack CD39 expression. In colorectal and lung tumours, the absence of CD39 in CD8+ TILs defines populations that lack hallmarks of chronic antigen stimulation at the tumour site, supporting their classification as bystanders. Expression of CD39 varied markedly between patients, with some patients having predominantly CD39- CD8+ TILs. Furthermore, frequencies of CD39 expression among CD8+ TILs correlated with several important clinical parameters, such as the mutation status of lung tumour epidermal growth factor receptors. Our results demonstrate that not all tumour-infiltrating T cells are specific for tumour antigens, and suggest that measuring CD39 expression could be a straightforward way to quantify or isolate bystander T cells.

The tumour microenvironment creates a niche for the self-renewal of tumour-promoting macrophages in colon adenoma.

Circulating CCR2+ monocytes are crucial for maintaining the adult tissue-resident F4/80hiMHCIIhi macrophage pool in the intestinal lamina propria. Here we show that a subpopulation of CCR2-independent F4/80hiMHCIIlow macrophages, which are the most abundant F4/80hi cells in neonates, gradually decline in number in adulthood; these macrophages likely represent the fetal contribution to F4/80hi cells. In colon adenomas of ApcMin/+ mice, F4/80hiMHCIIlow macrophages are not only preserved, but become the dominant subpopulation among tumour-resident macrophages during tumour progression. Furthermore, these pro-tumoural F4/80hiMHCIIlow and F4/80hiMHCIIhi macrophages can self-renew in the tumour and maintain their numbers mostly independent from bone marrow contribution. Analyses of colon adenomas indicate that CSF1 may be a key facilitator of macrophage self-renewal. In summary, the tumour microenvironment creates an isolated niche for tissue-resident macrophages that favours macrophage survival and self-renewal.

Human fetal dendritic cells promote prenatal T-cell immune suppression through arginase-2.

During gestation the developing human fetus is exposed to a diverse range of potentially immune-stimulatory molecules including semi-allogeneic antigens from maternal cells, substances from ingested amniotic fluid, food antigens, and microbes. Yet the capacity of the fetal immune system, including antigen-presenting cells, to detect and respond to such stimuli remains unclear. In particular, dendritic cells, which are crucial for effective immunity and tolerance, remain poorly characterized in the developing fetus. Here we show that subsets of antigen-presenting cells can be identified in fetal tissues and are related to adult populations of antigen-presenting cells. Similar to adult dendritic cells, fetal dendritic cells migrate to lymph nodes and respond to toll-like receptor ligation; however, they differ markedly in their response to allogeneic antigens, strongly promoting regulatory T-cell induction and inhibiting T-cell tumour-necrosis factor-α production through arginase-2 activity. Our results reveal a previously unappreciated role of dendritic cells within the developing fetus and indicate that they mediate homeostatic immune-suppressive responses during gestation.

E2F1 Orchestrates Transcriptomics and Oxidative Metabolism in Wharton's Jelly-Derived Mesenchymal Stem Cells from Growth-Restricted Infants.

Wharton's jelly-derived Mesenchymal Stem Cells (MSCs) isolated from newborns with intrauterine fetal growth restriction were previously shown to exert anabolic features including insulin hypersensitivity. Here, we extend these observations and demonstrate that MSCs from small for gestational age (SGA) individuals have decreased mitochondrial oxygen consumption rates. Comparing normally grown and SGA MSCs using next generation sequencing studies, we measured global transcriptomic and epigenetic profiles and identified E2F1 as an over-expressed transcription factor regulating oxidative metabolism in the SGA group. We further show that E2F1 regulates the differential transcriptome found in SGA derived MSCs and is associated with the activating histone marks H3K27ac and H3K4me3. One of the key genes regulated by E2F1 was found to be the fatty acid elongase ELOVL2, a gene involved in the endogenous synthesis of docosahexaenoic acid (DHA). Finally, we shed light on how the E2F1-ELOVL2 pathway may alter oxidative respiration in the SGA condition by contributing to the maintenance of cellular metabolic homeostasis.

Transcriptional and functional characterization of CD137L-dendritic cells identifies a novel dendritic cell phenotype.

The importance of monocyte-derived dendritic cells (DCs) is evidenced by the fact that they are essential for the elimination of pathogens. Although in vitro DCs can be generated by treatment of monocytes with GM-CSF and IL-4, it is unknown what stimuli induce differentiation of DCs in vivo. CD137L-DCs are human monocyte-derived DC that are generated by CD137 ligand (CD137L) signaling. We demonstrate that the gene signature of in vitro generated CD137L-DCs is most similar to those of GM-CSF and IL-4-generated immature DCs and of macrophages. This is reminiscent of in vivo inflammatory DC which also have been reported to share gene signatures with monocyte-derived DCs and macrophages. Performing direct comparison of deposited human gene expression data with a CD137L-DC dataset revealed a significant enrichment of CD137L-DC signature genes in inflammatory in vivo DCs. In addition, surface marker expression and cytokine secretion by CD137L-DCs resemble closely those of inflammatory DCs. Further, CD137L-DCs express high levels of adhesion molecules, display strong attachment, and employ the adhesion molecule ALCAM to stimulate T cell proliferation. This study characterizes the gene expression profile of CD137L-DCs, and identifies significant similarities of CD137L-DCs with in vivo inflammatory monocyte-derived DCs and macrophages.

Human dermal CD14⁺ cells are a transient population of monocyte-derived macrophages.

Dendritic cells (DCs), monocytes, and macrophages are leukocytes with critical roles in immunity and tolerance. The DC network is evolutionarily conserved; the homologs of human tissue CD141(hi)XCR1⁺ CLEC9A⁺ DCs and CD1c⁺ DCs are murine CD103⁺ DCs and CD64⁻ CD11b⁺ DCs. In addition, human tissues also contain CD14⁺ cells, currently designated as DCs, with an as-yet unknown murine counterpart. Here we have demonstrated that human dermal CD14⁺ cells are a tissue-resident population of monocyte-derived macrophages with a short half-life of <6 days. The decline and reconstitution kinetics of human blood CD14⁺ monocytes and dermal CD14⁺ cells in vivo supported their precursor-progeny relationship. The murine homologs of human dermal CD14⁺ cells are CD11b⁺ CD64⁺ monocyte-derived macrophages. Human and mouse monocytes and macrophages were defined by highly conserved gene transcripts, which were distinct from DCs. The demonstration of monocyte-derived macrophages in the steady state in human tissue supports a conserved organization of human and mouse mononuclear phagocyte system.

Human regulatory B cells combine phenotypic and genetic hallmarks with a distinct differentiation fate.

Regulatory B cells (B-reg) produce IL-10 and suppress inflammation in both mice and humans, but limited data on the phenotype and function of these cells have precluded detailed assessment of their contribution to host immunity. In this article, we report that human B-reg cannot be defined based on a phenotype composed of conventional B cell markers, and that IL-10 production can be elicited in both the CD27(+) memory population and naive B cell subset after only a brief stimulation in vitro. We therefore sought to obtain a better definition of IL-10-producing human B-regs using a multiparameter analysis of B cell phenotype, function, and gene expression profile. Exposure to CpG and anti-Ig are the most potent stimuli for IL-10 secretion in human B cells, but microarray analysis revealed that human B cells cotreated with these reagents resulted in only ∼0.7% of genes being differentially expressed between IL-10(+) and IL-10(-) cells. Instead, connectivity map analysis revealed that IL-10-secreting B cells are those undergoing specific differentiation toward a germinal center fate, and we identified a CD11c(+) B cell subset that was not capable of producing IL-10 even under optimal conditions. Our findings will assist in the identification of a broader range of human pro-B-reg populations that may represent novel targets for immunotherapy.

Multiple SVM-RFE for gene selection in cancer classification with expression data.

This paper proposes a new feature selection method that uses a backward elimination procedure similar to that implemented in support vector machine recursive feature elimination (SVM-RFE). Unlike the SVM-RFE method, at each step, the proposed approach computes the feature ranking score from a statistical analysis of weight vectors of multiple linear SVMs trained on subsamples of the original training data. We tested the proposed method on four gene expression datasets for cancer classification. The results show that the proposed feature selection method selects better gene subsets than the original SVM-RFE and improves the classification accuracy. A Gene Ontology-based similarity assessment indicates that the selected subsets are functionally diverse, further validating our gene selection method. This investigation also suggests that, for gene expression-based cancer classification, average test error from multiple partitions of training and test sets can be recommended as a reference of performance quality.

Proteomic cancer classification with mass spectrometry data.

The ultimate goal of cancer proteomics is to adapt proteomic technologies for routine use in clinical laboratories for the purpose of diagnostic and prognostic classification of disease states, as well as in evaluating drug toxicity and efficacy. Analysis of tumor-specific proteomic profiles may also allow better understanding of tumor development and the identification of novel targets for cancer therapy. The biological variability among patient samples as well as the huge dynamic range of biomarker concentrations are currently the main challenges facing efforts to deduce diagnostic patterns that are unique to specific disease states. While several strategies exist to address this problem, we focus here on cancer classification using mass spectrometry (MS) for proteomic profiling and biomarker identification. Recent advances in MS technology are starting to enable high-throughput profiling of the protein content of complex samples. For cancer classification, the protein samples from cancer patients and noncancer patients or from different cancer stages are analyzed through MS instruments and the MS patterns are used to build a diagnostic classifier. To illustrate the importance of feature selection in cancer classification, we present a method based on support vector machine-recursive feature elimination (SVM-RFE), demonstrated on two cancer datasets from ovarian and lung cancer.