Clonally expanded EOMES+ Tr1-like cells in primary and metastatic tumors are associated with disease progression.

Regulatory T (Treg) cells are a barrier for tumor immunity and a target for immunotherapy. Using single-cell transcriptomics, we found that CD4+ T cells infiltrating primary and metastatic colorectal cancer and non-small-cell lung cancer are highly enriched for two subsets of comparable size and suppressor function comprising forkhead box protein P3+ Treg and eomesodermin homolog (EOMES)+ type 1 regulatory T (Tr1)-like cells also expressing granzyme K and chitinase-3-like protein 2. EOMES+ Tr1-like cells, but not Treg cells, were clonally related to effector T cells and were clonally expanded in primary and metastatic tumors, which is consistent with their proliferation and differentiation in situ. Using chitinase-3-like protein 2 as a subset signature, we found that the EOMES+ Tr1-like subset correlates with disease progression but is also associated with response to programmed cell death protein 1-targeted immunotherapy. Collectively, these findings highlight the heterogeneity of Treg cells that accumulate in primary tumors and metastases and identify a new prospective target for cancer immunotherapy.

A Relay Pathway between Arginine and Tryptophan Metabolism Confers Immunosuppressive Properties on Dendritic Cells.

Arginase 1 (Arg1) and indoleamine 2,3-dioxygenase 1 (IDO1) are immunoregulatory enzymes catalyzing the degradation of l-arginine and l-tryptophan, respectively, resulting in local amino acid deprivation. In addition, unlike Arg1, IDO1 is also endowed with non-enzymatic signaling activity in dendritic cells (DCs). Despite considerable knowledge of their individual biology, no integrated functions of Arg1 and IDO1 have been reported yet. We found that IDO1 phosphorylation and consequent activation of IDO1 signaling in DCs was strictly dependent on prior expression of Arg1 and Arg1-dependent production of polyamines. Polyamines, either produced by DCs or released by bystander Arg1+ myeloid-derived suppressor cells, conditioned DCs toward an IDO1-dependent, immunosuppressive phenotype via activation of the Src kinase, which has IDO1-phosphorylating activity. Thus our data indicate that Arg1 and IDO1 are linked by an entwined pathway in immunometabolism and that their joint modulation could represent an important target for effective immunotherapy in several disease settings.

Single-Cell Sequencing of Mouse Heart Immune Infiltrate in Pressure Overload-Driven Heart Failure Reveals Extent of Immune Activation.

BACKGROUND: Inflammation is a key component of cardiac disease, with macrophages and T lymphocytes mediating essential roles in the progression to heart failure. Nonetheless, little insight exists on other immune subsets involved in the cardiotoxic response. METHODS: Here, we used single-cell RNA sequencing to map the cardiac immune composition in the standard murine nonischemic, pressure-overload heart failure model. By focusing our analysis on CD45+ cells, we obtained a higher resolution identification of the immune cell subsets in the heart, at early and late stages of disease and in controls. We then integrated our findings using multiparameter flow cytometry, immunohistochemistry, and tissue clarification immunofluorescence in mouse and human. RESULTS: We found that most major immune cell subpopulations, including macrophages, B cells, T cells and regulatory T cells, dendritic cells, Natural Killer cells, neutrophils, and mast cells are present in both healthy and diseased hearts. Most cell subsets are found within the myocardium, whereas mast cells are found also in the epicardium. Upon induction of pressure overload, immune activation occurs across the entire range of immune cell types. Activation led to upregulation of key subset-specific molecules, such as oncostatin M in proinflammatory macrophages and PD-1 in regulatory T cells, that may help explain clinical findings such as the refractivity of patients with heart failure to anti-tumor necrosis factor therapy and cardiac toxicity during anti-PD-1 cancer immunotherapy, respectively. CONCLUSIONS: Despite the absence of infectious agents or an autoimmune trigger, induction of disease leads to immune activation that involves far more cell types than previously thought, including neutrophils, B cells, Natural Killer cells, and mast cells. This opens up the field of cardioimmunology to further investigation by using toolkits that have already been developed to study the aforementioned immune subsets. The subset-specific molecules that mediate their activation may thus become useful targets for the diagnostics or therapy of heart failure.

Background fluorescence and spreading error are major contributors of variability in high-dimensional flow cytometry data visualization by t-distributed stochastic neighboring embedding.

Multidimensional single-cell analysis requires approaches to visualize complex data in intuitive 2D graphs. In this regard, t-distributed stochastic neighboring embedding (tSNE) is the most popular algorithm for single-cell RNA sequencing and cytometry by time-of-flight (CyTOF), but its application to polychromatic flow cytometry, including the recently developed 30-parameter platform, is still under investigation. We identified differential distribution of background values between samples, generated by either background calculation or spreading error (SE), as a major source of variability in polychromatic flow cytometry data representation by tSNE, ultimately resulting in the identification of erroneous heterogeneity among cell populations. Biexponential transformation of raw data and limiting SE during panel development dramatically improved data visualization. These aspects must be taken into consideration when using computational approaches as discovery tools in large sets of samples from independent experiments or immunomonitoring in clinical trials.

The early expansion of anergic NKG2Apos/CD56dim/CD16neg natural killer represents a therapeutic target in haploidentical hematopoietic stem cell transplantation.

Natural killer cells are the first lymphocyte population to reconstitute early after non-myeloablative and T cell-replete haploidentical hematopoietic stem cell transplantation with post-transplant infusion of cyclophosphamide. The study herein characterizes the transient and predominant expansion starting from the second week following haploidentical hematopoietic stem cell transplantation of a donor-derived unconventional subset of NKp46neg-low/CD56dim/CD16neg natural killer cells expressing remarkably high levels of CD94/NKG2A. Both transcription and phenotypic profiles indicated that unconventional NKp46neg-low/CD56dim/CD16neg cells are a distinct natural killer cell subpopulation with features of late stage differentiation, yet retaining proliferative capability and functional plasticity to generate conventional NKp46pos/CD56bright/CD16neg-low cells in response to interleukin-15 plus interleukin-18. While present at low frequency in healthy donors, unconventional NKp46neg-low/CD56dim/CD16neg cells are greatly expanded in the seven weeks following haploidentical hematopoietic stem cell transplantation, and express high levels of the activating receptors NKG2D and NKp30 as well as of the lytic granules Granzyme-B and Perforin. Nonetheless, NKp46neg-low/CD56dim/CD16neg cells displayed a markedly defective cytotoxicity that could be reversed by blocking the inhibitory receptor CD94/NKG2A. These data open new and important perspectives to better understand the ontogenesis/homeostasis of human natural killer cells and to develop a novel immune-therapeutic approach that targets the inhibitory NKG2A check-point, thus unleashing natural killer cell alloreactivity early after haploidentical hematopoietic stem cell transplantation.

A comparative transcriptomic analysis of astrocytes differentiation from human neural progenitor cells.

Astrocytes are a morphologically and functionally heterogeneous population of cells that play critical roles in neurodevelopment and in the regulation of central nervous system homeostasis. Studies of human astrocytes have been hampered by the lack of specific molecular markers and by the difficulties associated with purifying and culturing astrocytes from adult human brains. Human neural progenitor cells (NPCs) with self-renewal and multipotent properties represent an appealing model system to gain insight into the developmental genetics and function of human astrocytes, but a comprehensive molecular characterization that confirms the validity of this cellular system is still missing. Here we used an unbiased transcriptomic analysis to characterize in vitro culture of human NPCs and to define the gene expression programs activated during the differentiation of these cells into astrocytes using FBS or the combination of CNTF and BMP4. Our results demonstrate that in vitro cultures of human NPCs isolated during the gliogenic phase of neurodevelopment mainly consist of radial glial cells (RGCs) and glia-restricted progenitor cells. In these cells the combination of CNTF and BMP4 activates the JAK/STAT and SMAD signaling cascades, leading to the inhibition of oligodendrocytes lineage commitment and activation of astrocytes differentiation. On the other hand, FBS-derived astrocytes have properties of reactive astrocytes. Our work suggests that in vitro culture of human NPCs represents a valuable cellular system to study human disorders characterized by impairment of astrocytes development and function. Our datasets represent an important resource for researchers studying human astrocytes development and might set the basis for the discovery of novel human-specific astrocyte markers.

CD4+ Regulatory T Cells in Human Cancer: Subsets, Origin, and Molecular Regulation.

CD4+CD25hiFOXP3+ regulatory T cells (Treg) play major roles in the maintenance of immune tolerance, prevention of inflammation, and tissue homeostasis and repair. In contrast with these beneficial roles, Tregs are abundant in virtually all tumors and have been mechanistically linked to disease progression, metastases development, and therapy resistance. Tregs are thus recognized as a major target for cancer immunotherapy. Compared with other sites in the body, tumors harbor hyperactivated Treg subsets whose molecular characteristics are only beginning to be elucidated. Here, we describe current knowledge of intratumoral Tregs and discuss their potential cellular and tissue origin. Furthermore, we describe currently recognized molecular regulators that drive differentiation and maintenance of Tregs in cancer, with a special focus on those signals regulating their chronic immune activation, with relevant implications for cancer progression and therapy.

CRUSTY: a versatile web platform for the rapid analysis and visualization of high-dimensional flow cytometry data.

Flow cytometry (FCM) can investigate dozens of parameters from millions of cells and hundreds of specimens in a short time and at a reasonable cost, but the amount of data that is generated is considerable. Computational approaches are useful to identify novel subpopulations and molecular biomarkers, but generally require deep expertize in bioinformatics and the use of different platforms. To overcome these limitations, we introduce CRUSTY, an interactive, user-friendly webtool incorporating the most popular algorithms for FCM data analysis, and capable of visualizing graphical and tabular results and automatically generating publication-quality figures within minutes. CRUSTY also hosts an interactive interface for the exploration of results in real time. Thus, CRUSTY enables a large number of users to mine complex datasets and reduce the time required for data exploration and interpretation. CRUSTY is accessible at https://crusty.humanitas.it/ .

Single-cell profiling defines the prognostic benefit of CD39high tissue resident memory CD8+ T cells in luminal-like breast cancer.

Luminal-like breast cancer (BC) constitutes the majority of BC subtypes, but, differently from highly aggressive triple negative BC, is poorly infiltrated by the immune system. The quality of the immune infiltrate in luminal-like BCs has been poorly studied, thereby limiting further investigation of immunotherapeutic strategies. By using high-dimensional single-cell technologies, we identify heterogeneous behavior within the tissue-resident memory CD8+ T (Trm) cells infiltrating luminal-like tumors. A subset of CD127- CD39hi Trm cells, preferentially present in the tumor compared to the adjacent normal breast tissue or peripheral blood, retains enhanced degranulation capacity compared to the CD127+ CD39lo Trm counterpart ex vivo, and is specifically associated with positive prognosis. Nevertheless, such prognostic benefit is lost in the presence of highly-suppressive CCR8hi ICOShi IRF4+ effector Tregs. Thus, combinatorial strategies aiming at boosting Trm function and infiltration while relieving from Treg-mediated immunosuppression should be investigated to achieve proper tumor control in luminal-like BCs.

TIM4 expression by dendritic cells mediates uptake of tumor-associated antigens and anti-tumor responses.

Acquisition of cell-associated tumor antigens by type 1 dendritic cells (cDC1) is essential to induce and sustain tumor specific CD8+ T cells via cross-presentation. Here we show that capture and engulfment of cell associated antigens by tissue resident lung cDC1 is inhibited during progression of mouse lung tumors. Mechanistically, loss of phagocytosis is linked to tumor-mediated downregulation of the phosphatidylserine receptor TIM4, that is highly expressed in normal lung resident cDC1. TIM4 receptor blockade and conditional cDC1 deletion impair activation of tumor specific CD8+ T cells and promote tumor progression. In human lung adenocarcinomas, TIM4 transcripts increase the prognostic value of a cDC1 signature and predict responses to PD-1 treatment. Thus, TIM4 on lung resident cDC1 contributes to immune surveillance and its expression is suppressed in advanced tumors.

Cell-Type-Specific Analysis of Molecular Pathology in Autism Identifies Common Genes and Pathways Affected Across Neocortical Regions.

Despite its heterogeneity, autism is characterized by a defined behavioral phenotype, suggesting that the molecular pathology affects specific neural substrates to cause behavioral dysfunction. Previous studies identified genes dysregulated in autism cortex but did not address their cell-type specificity. Moreover, it is unknown whether there is a core of genes dysregulated across multiple neocortical regions. We applied RNA sequencing to postmortem brain tissue samples from autism patients and neurologically normal controls and combined our data with previously published datasets. We then identified genes, pathways, and alternative splicing events which are dysregulated in five neocortical regions in autism. To gain information about cell-type specificity of the dysregulated genes, we analyzed single-nuclei RNA sequencing data of adult human cortex and intersected cell-type-specific gene signatures with genes dysregulated in autism in specific cortical regions. We found that autism-associated gene expression changes across 4 frontal and temporal cortex regions converge on 27 genes related to immune response and enriched in human astrocytes, microglia, and brain endothelium. Shared splicing changes, however, are found in genes predominantly associated with synaptic function and adult interneurons and projection neurons. Finally, we demonstrate that regions of DNA differentially methylated in autism overlap genes associated with development and enriched in human cortical oligodendrocytes. Our study identifies signatures of autism molecular pathology shared across neocortical regions, as well as neural cell types enriched for common dysregulated genes, thus paving way for assessing cell-type-specific mechanisms of autism pathology.

P2X7 Receptor Activity Limits Accumulation of T Cells within Tumors.

Extracellular ATP (eATP) is a signaling molecule that variably affects all cells of the immune system either directly or after hydrolysis to adenosine. Although eATP is virtually absent in the interstitium of normal tissues, it can be present in the hundreds of micromolar range in tumors, a concentration compatible with activation of the ATP-gated ionotropic P2X7 receptor. Here, we show that P2X7 activity in tumor-infiltrating lymphocytes (TIL) induces cellular senescence and limits tumor suppression. P2X7 stimulation affected cell cycling of effector T cells and resulted in generation of mitochondrial reactive oxygen species and p38 MAPK-dependent upregulation of cyclin-dependent kinase inhibitor 1A (Cdkn1a, encoding for p21Waf1/Cip1). Lack of P2X7 promoted a transcriptional signature that correlated with enhanced cytotoxic T-cell response in human solid tumors. In mice, transfer of tumor-specific T cells with deletion of P2rx7 significantly reduced tumor growth and extended survival. Collectively, these findings uncover a purinergic checkpoint that can be targeted to improve the efficacy of cancer immunotherapy strategies. SIGNIFICANCE: These findings suggest that the purinergic checkpoint P2X7 may be targeted to enhance T-cell-mediated cancer immunotherapy and improve T effector cell accumulation in the tumor microenvironment. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/18/3906/F1.large.jpg.

PI3K Inhibitors Curtail MYC-Dependent Mutant p53 Gain-of-Function in Head and Neck Squamous Cell Carcinoma.

PURPOSE: Mutation of TP53 gene is a hallmark of head and neck squamous cell carcinoma (HNSCC) not yet exploited therapeutically. TP53 mutation frequently leads to the synthesis of mutant p53 proteins with gain-of-function activity, associated with radioresistance and high incidence of local recurrences in HNSCC. EXPERIMENTAL DESIGN: Mutant p53-associated functions were investigated through gene set enrichment analysis in the Cancer Genome Atlas cohort of HNSCC and in a panel of 22 HNSCC cell lines. Mutant p53-dependent transcripts were analyzed in HNSCC cell line Cal27, carrying mutant p53H193L; FaDu, carrying p53R248L; and Detroit 562, carrying p53R175H. Drugs impinging on mutant p53-MYC-dependent signature were identified interrogating Connectivity Map (https://clue.io) derived from the Library of Integrated Network-based Cellular Signatures (LINCS) database (http://lincs.hms.harvard.edu/) and analyzed in HNSCC cell lines and patient-derived xenografts (PDX) models. RESULTS: We identified a signature of transcripts directly controlled by gain-of-function mutant p53 protein and prognostic in HNSCC, which is highly enriched of MYC targets. Specifically, both in PDX and cell lines of HNSCC treated with the PI3Kα-selective inhibitor BYL719 (alpelisib) the downregulation of mutant p53/MYC-dependent signature correlates with response to this compound. Mechanistically, mutant p53 favors the binding of MYC to its target promoters and enhances MYC protein stability. Treatment with BYL719 disrupts the interaction of MYC, mutant p53, and YAP proteins with MYC target promoters. Of note, depletion of MYC, mutant p53, or YAP potentiates the effectiveness of BYL719 treatment. CONCLUSIONS: Collectively, the blocking of this transcriptional network is an important determinant for the response to BYL719 in HNSCC.

Global chromatin conformation differences in the Drosophila dosage compensated chromosome X.

In Drosophila melanogaster the single male chromosome X undergoes an average twofold transcriptional upregulation for balancing the transcriptional output between sexes. Previous literature hypothesised that a global change in chromosome structure may accompany this process. However, recent studies based on Hi-C failed to detect these differences. Here we show that global conformational differences are specifically present in the male chromosome X and detectable using Hi-C data on sex-sorted embryos, as well as male and female cell lines, by leveraging custom data analysis solutions. We find the male chromosome X has more mid-/long-range interactions. We also identify differences at structural domain boundaries containing BEAF-32 in conjunction with CP190 or Chromator. Weakening of these domain boundaries in male chromosome X co-localizes with the binding of the dosage compensation complex and its co-factor CLAMP, reported to enhance chromatin accessibility. Together, our data strongly indicate that chromosome X dosage compensation affects global chromosome structure.

Chromosome positioning in interphase nuclei of hematopoietic stem cell and myeloid precursor.

Human myelopoiesis is an intriguing biological process during which multipotent stem cells limit their differentiation potential generating precursors that evolve into terminally differentiated cells. The differentiation process is correlated with differential gene expression and changes in nuclear architecture. In interphase, chromosomes are distinct entities known as chromosome territories and they show a radial localization that could result in a constrain of inter-homologous distance. This element plays a role in genome stability and gene expression. Here, we provide the first experimental evidence of 3D chromosomal arrangement considering two steps of human normal myelopoiesis. Specifically, multicolor 3D-FISH and 3D image analysis revealed that, in both normal human hematopoietic stem cells and myelod precursors CD14-, chromosomal position is correlated with gene density. However, we observed that inter-homologue distances are totally different during differentiation. This could be associated with differential gene expression that we found comparing the two cell types. Our results disclose an unprecedented framework relevant for deciphering the genomic mechanisms at the base of normal human myelopoiesis.

P2RY1/ALK3-Expressing Cells within the Adult Human Exocrine Pancreas Are BMP-7 Expandable and Exhibit Progenitor-like Characteristics.

Treatment of human pancreatic non-endocrine tissue with Bone Morphogenetic Protein 7 (BMP-7) leads to the formation of glucose-responsive ß-like cells. Here, we show that BMP-7 acts on extrainsular cells expressing PDX1 and the BMP receptor activin-like kinase 3 (ALK3/BMPR1A). In vitro lineage tracing indicates that ALK3+ cell populations are multipotent. PDX1+/ALK3+ cells are absent from islets but prominently represented in the major pancreatic ducts and pancreatic duct glands. We identified the purinergic receptor P2Y1 (P2RY1) as a surrogate surface marker for PDX1. Sorted P2RY1+/ALK3bright+ cells form BMP-7-expandable colonies characterized by NKX6.1 and PDX1 expression. Unlike the negative fraction controls, these colonies can be differentiated into multiple pancreatic lineages upon BMP-7 withdrawal. RNA-seq further corroborates the progenitor-like nature of P2RY1+/ALK3bright+ cells and their multilineage differentiation potential. Our studies confirm the existence of progenitor cells in the adult human pancreas and suggest a specific anatomical location within the ductal and glandular networks.

Skeletal muscle characteristics are preserved in hTERT/cdk4 human myogenic cell lines.

BACKGROUND: hTERT/cdk4 immortalized myogenic human cell lines represent an important tool for skeletal muscle research, being used as therapeutically pertinent models of various neuromuscular disorders and in numerous fundamental studies of muscle cell function. However, the cell cycle is linked to other cellular processes such as integrin regulation, the PI3K/Akt pathway, and microtubule stability, raising the question as to whether genetic modification related to the cell cycle results in secondary effects that could undermine the validity of these cell models. RESULTS: Here we subjected five healthy and disease muscle cell isolates to transcriptomic analysis, comparing immortalized lines with their parent primary populations in both differentiated and undifferentiated states, and testing their myogenic character by comparison with non-myogenic (CD56-negative) cells. Principal component analysis of global gene expression showed tight clustering of immortalized myoblasts to their parent primary populations, with clean separation from the non-myogenic reference. Comparison was made to publicly available transcriptomic data from studies of muscle human pathology, cell, and animal models, including to derive a consensus set of genes previously shown to have altered regulation during myoblast differentiation. Hierarchical clustering of samples based on gene expression of this consensus set showed that immortalized lines retained the myogenic expression patterns of their parent primary populations. Of 2784 canonical pathways and gene ontology terms tested by gene set enrichment analysis, none were significantly enriched in immortalized compared to primary cell populations. We observed, at the whole transcriptome level, a strong signature of cell cycle shutdown associated with senescence in one primary myoblast population, whereas its immortalized clone was protected. CONCLUSIONS: Immortalization had no observed effect on the myogenic cascade or on any other cellular processes, and it was protective against the systems level effects of senescence that are observed at higher division counts of primary cells.

Aging: a portrait from gene expression profile in blood cells.

The availability of reliable biomarkers of aging is important not only to monitor the effect of interventions and predict the timing of pathologies associated with aging but also to understand the mechanisms and devise appropriate countermeasures. Blood cells provide an easily available tissue and gene expression profiles from whole blood samples appear to mirror disease states and some aspects of the aging process itself. We report here a microarray analysis of whole blood samples from two cohorts of healthy adult and elderly subjects, aged 43±3 and 68±4 years, respectively, to monitor gene expression changes in the initial phase of the senescence process. A number of significant changes were found in the elderly compared to the adult group, including decreased levels of transcripts coding for components of the mitochondrial respiratory chain, which correlate with a parallel decline in the maximum rate of oxygen consumption (VO2max), as monitored in the same subjects. In addition, blood cells show age-related changes in the expression of several markers of immunosenescence, inflammation and oxidative stress. These findings support the notion that the immune system has a major role in tissue homeostasis and repair, which appears to be impaired since early stages of the aging process.

Genome-Wide Definition of Promoter and Enhancer Usage during Neural Induction of Human Embryonic Stem Cells.

Genome-wide mapping of transcriptional regulatory elements is an essential tool for understanding the molecular events orchestrating self-renewal, commitment and differentiation of stem cells. We combined high-throughput identification of transcription start sites with genome-wide profiling of histones modifications to map active promoters and enhancers in embryonic stem cells (ESCs) induced to neuroepithelial-like stem cells (NESCs). Our analysis showed that most promoters are active in both cell types while approximately half of the enhancers are cell-specific and account for most of the epigenetic changes occurring during neural induction, and most likely for the modulation of the promoters to generate cell-specific gene expression programs. Interestingly, the majority of the promoters activated or up-regulated during neural induction have a "bivalent" histone modification signature in ESCs, suggesting that developmentally-regulated promoters are already poised for transcription in ESCs, which are apparently pre-committed to neuroectodermal differentiation. Overall, our study provides a collection of differentially used enhancers, promoters, transcription starts sites, protein-coding and non-coding RNAs in human ESCs and ESC-derived NESCs, and a broad, genome-wide description of promoter and enhancer usage and of gene expression programs characterizing the transition from a pluripotent to a neural-restricted cell fate.