Systemic vaccination induces CD8+ T cells and remodels the tumor microenvironment.

Therapeutic cancer vaccines are designed to increase tumor-specific T cell immunity. However, suppressive mechanisms within the tumor microenvironment (TME) may limit T cell function. Here, we assessed how the route of vaccination alters intratumoral myeloid cells. Using a self-assembling nanoparticle vaccine that links tumor antigen peptides to a Toll-like receptor 7/8 agonist (SNP-7/8a), we treated tumor-bearing mice subcutaneously (SNP-SC) or intravenously (SNP-IV). Both routes generated antigen-specific CD8+ T cells that infiltrated tumors. However, only SNP-IV mediated tumor regression, dependent on systemic type I interferon at the time of boost. Single-cell RNA-sequencing revealed that intratumoral monocytes expressing an immunoregulatory gene signature (Chil3, Anxa2, Wfdc17) were reduced after SNP-IV boost. In humans, the Chil3+ monocyte gene signature is enriched in CD16- monocytes and associated with worse outcomes. Our results show that the generation of tumor-specific CD8+ T cells combined with remodeling of the TME is a promising approach for tumor immunotherapy.

Gremlin 1+ fibroblastic niche maintains dendritic cell homeostasis in lymphoid tissues.

Fibroblastic reticular cells (FRCs) are specialized stromal cells that define tissue architecture and regulate lymphocyte compartmentalization, homeostasis, and innate and adaptive immunity in secondary lymphoid organs (SLOs). In the present study, we used single-cell RNA sequencing (scRNA-seq) of human and mouse lymph nodes (LNs) to identify a subset of T cell-zone FRCs defined by the expression of Gremlin1 (Grem1) in both species. Grem1-CreERT2 knock-in mice enabled localization, multi-omics characterization and genetic depletion of Grem1+ FRCs. Grem1+ FRCs primarily localize at T-B cell junctions of SLOs, neighboring pre-dendritic cells and conventional dendritic cells (cDCs). As such, their depletion resulted in preferential loss and decreased homeostatic proliferation and survival of resident cDCs and compromised T cell immunity. Trajectory analysis of human LN scRNA-seq data revealed expression similarities to murine FRCs, with GREM1+ cells marking the endpoint of both trajectories. These findings illuminate a new Grem1+ fibroblastic niche in LNs that functions to maintain the homeostasis of lymphoid tissue-resident cDCs.

LRRC15+ myofibroblasts dictate the stromal setpoint to suppress tumour immunity.

Recent single-cell studies of cancer in both mice and humans have identified the emergence of a myofibroblast population specifically marked by the highly restricted leucine-rich-repeat-containing protein 15 (LRRC15)1-3. However, the molecular signals that underlie the development of LRRC15+ cancer-associated fibroblasts (CAFs) and their direct impact on anti-tumour immunity are uncharacterized. Here in mouse models of pancreatic cancer, we provide in vivo genetic evidence that TGFß receptor type 2 signalling in healthy dermatopontin+ universal fibroblasts is essential for the development of cancer-associated LRRC15+ myofibroblasts. This axis also predominantly drives fibroblast lineage diversity in human cancers. Using newly developed Lrrc15-diphtheria toxin receptor knock-in mice to selectively deplete LRRC15+ CAFs, we show that depletion of this population markedly reduces the total tumour fibroblast content. Moreover, the CAF composition is recalibrated towards universal fibroblasts. This relieves direct suppression of tumour-infiltrating CD8+ T cells to enhance their effector function and augments tumour regression in response to anti-PDL1 immune checkpoint blockade. Collectively, these findings demonstrate that TGFß-dependent LRRC15+ CAFs dictate the tumour-fibroblast setpoint to promote tumour growth. These cells also directly suppress CD8+ T cell function and limit responsiveness to checkpoint blockade. Development of treatments that restore the homeostatic fibroblast setpoint by reducing the population of pro-disease LRRC15+ myofibroblasts may improve patient survival and response to immunotherapy.

Cross-tissue organization of the fibroblast lineage.

Fibroblasts are non-haematopoietic structural cells that define the architecture of organs, support the homeostasis of tissue-resident cells and have key roles in fibrosis, cancer, autoimmunity and wound healing1. Recent studies have described fibroblast heterogeneity within individual tissues1. However, the field lacks a characterization of fibroblasts at single-cell resolution across tissues in healthy and diseased organs. Here we constructed fibroblast atlases by integrating single-cell transcriptomic data from about 230,000 fibroblasts across 17 tissues, 50 datasets, 11 disease states and 2 species. Mouse fibroblast atlases and a DptIRESCreERT2 knock-in mouse identified two universal fibroblast transcriptional subtypes across tissues. Our analysis suggests that these cells can serve as a reservoir that can yield specialized fibroblasts across a broad range of steady-state tissues and activated fibroblasts in disease. Comparison to an atlas of human fibroblasts from perturbed states showed that fibroblast transcriptional states are conserved between mice and humans, including universal fibroblasts and activated phenotypes associated with pathogenicity in human cancer, fibrosis, arthritis and inflammation. In summary, a cross-species and pan-tissue approach to transcriptomics at single-cell resolution has identified key organizing principles of the fibroblast lineage in health and disease.

A bird's eye view of fibroblast heterogeneity: A pan-disease, pan-cancer perspective.

Fibroblasts, custodians of tissue architecture and function, are no longer considered a monolithic entity across tissues and disease indications. Recent advances in single-cell technologies provide an unrestricted, high-resolution view of fibroblast heterogeneity that exists within and across tissues. In this review, we summarize a compendium of single-cell transcriptomic studies and provide a comprehensive accounting of fibroblast subsets, many of which have been described to occupy specific niches in tissues at homeostatic and pathologic states. Understanding this heterogeneity is particularly important in the context of cancer, as the diverse cancer-associated fibroblast (CAF) phenotypes in the tumor microenvironment (TME) are directly impacted by the expression phenotypes of their predecessors. Relationships between these heterogeneous populations often accompany and influence response to therapy in cancer and fibrosis. We further highlight the importance of integrating single-cell studies to deduce common fibroblast phenotypes across disease states, which will facilitate the identification of common signaling pathways, gene regulatory programs, and cell surface markers that are going to advance drug discovery and targeting.

MHC class I H2-Kb negatively regulates neural progenitor cell proliferation by inhibiting FGFR signaling.

Proteins of the major histocompatibility complex class I (MHC I), predominantly known for antigen presentation in the immune system, have recently been shown to be necessary for developmental neural refinement and adult synaptic plasticity. However, their roles in nonneuronal cell populations in the brain remain largely unexplored. Here, we identify classical MHC I molecule H2-Kb as a negative regulator of proliferation in neural stem and progenitor cells (NSPCs). Using genetic knockout mouse models and in vivo viral-mediated RNA interference (RNAi) and overexpression, we delineate a role for H2-Kb in negatively regulating NSPC proliferation and adult hippocampal neurogenesis. Transcriptomic analysis of H2-Kb knockout NSPCs, in combination with in vitro RNAi, overexpression, and pharmacological approaches, further revealed that H2-Kb inhibits cell proliferation by dampening signaling pathways downstream of fibroblast growth factor receptor 1 (Fgfr1). These findings identify H2-Kb as a critical regulator of cell proliferation through the modulation of growth factor signaling.

A membrane protein display platform for receptor interactome discovery.

Cell surface receptors are critical for cell signaling and constitute a quarter of all human genes. Despite their importance and abundance, receptor interaction networks remain understudied because of difficulties associated with maintaining membrane proteins in their native conformation and their typically weak interactions. To overcome these challenges, we developed an extracellular vesicle-based method for membrane protein display that enables purification-free and high-throughput detection of receptor-ligand interactions in membranes. We demonstrate that this platform is broadly applicable to a variety of membrane proteins, enabling enhanced detection of extracellular interactions over a wide range of binding affinities. We were able to recapitulate and expand the interactome for prominent members of the B7 family of immunoregulatory proteins such as PD-L1/CD274 and B7-H3/CD276. Moreover, when applied to the orphan cancer-associated fibroblast protein, LRRC15, we identified a membrane-dependent interaction with the tumor stroma marker TEM1/CD248. Furthermore, this platform enabled profiling of cellular receptors for target-expressing as well as endogenous extracellular vesicles. Overall, this study presents a sensitive and easy to use screening platform that bypasses membrane protein purification and enables characterization of interactomes for any cell surface-expressed target of interest in its native state.

miRNA-independent function of long noncoding pri-miRNA loci.

Among the large, diverse set of mammalian long noncoding RNAs (lncRNAs), long noncoding primary microRNAs (lnc-pri-miRNAs) are those that host miRNAs. Whether lnc-pri-miRNA loci have important biological function independent of their cognate miRNAs is poorly understood. From a genome-scale lncRNA screen, lnc-pri-miRNA loci were enriched for function in cell proliferation, and in glioblastoma (i.e., GBM) cells with DGCR8 or DROSHA knockdown, lnc-pri-miRNA screen hits still regulated cell growth. To molecularly dissect the function of a lnc-pri-miRNA locus, we studied LOC646329 (also known as MIR29HG), which hosts the miR-29a/b1 cluster. In GBM cells, LOC646329 knockdown reduced miR-29a/b1 levels, and these cells exhibited decreased growth. However, genetic deletion of the miR-29a/b1 cluster (LOC646329-miR29Δ) did not decrease cell growth, while knockdown of LOC646329-miR29Δ transcripts reduced cell proliferation. The miR-29a/b1-independent activity of LOC646329 corresponded to enhancer-like activation of a neighboring oncogene (MKLN1), regulating cell propagation. The LOC646329 locus interacts with the MKLN1 promoter, and antisense oligonucleotide knockdown of the lncRNA disrupts these interactions and reduces the enhancer-like activity. More broadly, analysis of genome-wide data from multiple human cell types showed that lnc-pri-miRNA loci are significantly enriched for DNA looping interactions with gene promoters as well as genomic and epigenetic characteristics of transcriptional enhancers. Functional studies of additional lnc-pri-miRNA loci demonstrated cognate miRNA-independent enhancer-like activity. Together, these data demonstrate that lnc-pri-miRNA loci can regulate cell biology via both miRNA-dependent and miRNA-independent mechanisms.

CONICS integrates scRNA-seq with DNA sequencing to map gene expression to tumor sub-clones.

Motivation: Single-cell RNA-sequencing (scRNA-seq) has enabled studies of tissue composition at unprecedented resolution. However, the application of scRNA-seq to clinical cancer samples has been limited, partly due to a lack of scRNA-seq algorithms that integrate genomic mutation data. Results: To address this, we present. CONICS: COpy-Number analysis In single-Cell RNA-Sequencing. CONICS is a software tool for mapping gene expression from scRNA-seq to tumor clones and phylogenies, with routines enabling: the quantitation of copy-number alterations in scRNA-seq, robust separation of neoplastic cells from tumor-infiltrating stroma, inter-clone differential-expression analysis and intra-clone co-expression analysis. Availability and implementation: CONICS is written in Python and R, and is available from https://github.com/diazlab/CONICS. Supplementary information: Supplementary data are available at Bioinformatics online.

Transcriptomic Analyses of Adipocyte Differentiation From Human Mesenchymal Stromal-Cells (MSC).

Adipogenesis is a physiological process required for fat-tissue development, mainly involved in regulating the organism energetic-state. Abnormal distribution-changes and dysfunctions in such tissue are associated to different pathologies. Adipocytes are generated from progenitor cells, via a complex differentiating process not yet well understood. Therefore, we investigated differential mRNA and miRNA expression patterns of human mesenchymal stromal-cells (MSC) induced and not induced to differentiate into adipocytes by next (second)-generation sequencing. A total of 2,866 differentially expressed genes (101 encoding miRNA) were identified, with 705 (46 encoding miRNA) being upregulated in adipogenesis. They were related to different pathways, including PPARG, lipid, carbohydrate and energy metabolism, redox, membrane-organelle biosynthesis, and endocrine system. Downregulated genes were related to extracellular matrix and cell migration, proliferation, and differentiation. Analyses of mRNA-miRNA interaction showed that repressed miRNA-encoding genes can act downregulating PPARG-related genes; mostly the PPARG activator (PPARGC1A). Induced miRNA-encoding genes regulate downregulated genes related to TGFB1. These results shed new light to understand adipose-tissue differentiation and physiology, increasing our knowledge about pathologies like obesity, type-2 diabetes and osteoporosis. J. Cell. Physiol. 232: 771-784, 2017. © 2016 Wiley Periodicals, Inc.

Single-cell sequencing maps gene expression to mutational phylogenies in PDGF- and EGF-driven gliomas.

Glioblastoma multiforme (GBM) is the most common and aggressive type of primary brain tumor. Epidermal growth factor (EGF) and platelet-derived growth factor (PDGF) receptors are frequently amplified and/or possess gain-of-function mutations in GBM However, clinical trials of tyrosine-kinase inhibitors have shown disappointing efficacy, in part due to intra-tumor heterogeneity. To assess the effect of clonal heterogeneity on gene expression, we derived an approach to map single-cell expression profiles to sequentially acquired mutations identified from exome sequencing. Using 288 single cells, we constructed high-resolution phylogenies of EGF-driven and PDGF-driven GBMs, modeling transcriptional kinetics during tumor evolution. Descending the phylogenetic tree of a PDGF-driven tumor corresponded to a progressive induction of an oligodendrocyte progenitor-like cell type, expressing pro-angiogenic factors. In contrast, phylogenetic analysis of an EGFR-amplified tumor showed an up-regulation of pro-invasive genes. An in-frame deletion in a specific dimerization domain of PDGF receptor correlates with an up-regulation of growth pathways in a proneural GBM and enhances proliferation when ectopically expressed in glioma cell lines. In-frame deletions in this domain are frequent in public GBM data.

Systematic screen of chemotherapeutics in Drosophila stem cell tumors.

Here we report the development of an in vivo system to study the interaction of stem cells with drugs using a tumor model in the adult Drosophila intestine. Strikingly, we find that some Food and Drug Administration-approved chemotherapeutics that can inhibit the growth of Drosophila tumor stem cells can paradoxically promote the hyperproliferation of their wild-type counterparts. These results reveal an unanticipated side effect on stem cells that may contribute to tumor recurrence. We propose that the same side effect may occur in humans based on our finding that it is driven in Drosophila by the evolutionarily conserved Janus kinase-signal transducers and activators of transcription (JAK-STAT) pathway. An immediate implication of our findings is that supplementing traditional chemotherapeutics with anti-inflammatories may reduce tumor recurrence.

Principles of miRNA-mRNA interactions: beyond sequence complementarity.

MicroRNAs (miRNAs) are small non-coding RNAs that post-transcriptionally regulate gene expression by altering the translation efficiency and/or stability of targeted mRNAs. In vertebrates, more than 50% of all protein-coding RNAs are assumed to be subject to miRNA-mediated control, but current high-throughput methods that reliably measure miRNA-mRNA interactions either require prior knowledge of target mRNAs or elaborate preparation procedures. Consequently, experimentally validated interactions are relatively rare. Furthermore, in silico prediction based on sequence complementarity of miRNAs and their corresponding target sites suffers from extremely high false positive rates. Apparently, sequence complementarity alone is often insufficient to reflect the complex post-transcriptional regulation of mRNAs by miRNAs, which is especially true for animals. Therefore, combined analysis of small non-coding and protein-coding RNAs is indispensable to better understand and predict the complex dynamics of miRNA-regulated gene expression. Single-nucleotide polymorphisms (SNPs) and alternative polyadenylation (APA) can affect miRNA binding of a given transcript from different individuals and tissues, and especially APA is currently emerging as a major factor that contributes to variations in miRNA-mRNA interplay in animals. In this review, we focus on the influence of APA and SNPs on miRNA-mediated gene regulation and discuss the computational approaches that take these mechanisms into account.

Dual 3'Seq using deepSuperSAGE uncovers transcriptomes of interacting Salmonella enterica Typhimurium and human host cells.

BACKGROUND: The interaction of eukaryotic host and prokaryotic pathogen cells is linked to specific changes in the cellular proteome, and consequently to infection-related gene expression patterns of the involved cells. To simultaneously assess the transcriptomes of both organisms during their interaction we developed dual 3'Seq, a tag-based sequencing protocol that allows for exact quantification of differentially expressed transcripts in interacting pro- and eukaryotic cells without prior fixation or physical disruption of the interaction. RESULTS: Human epithelial cells were infected with Salmonella enterica Typhimurium as a model system for invasion of the intestinal epithelium, and the transcriptional response of the infected host cells together with the differential expression of invading and intracellular pathogen cells was determined by dual 3'Seq coupled with the next-generation sequencing-based transcriptome profiling technique deepSuperSAGE (deep Serial Analysis of Gene Expression). Annotation to reference transcriptomes comprising the operon structure of the employed S. enterica Typhimurium strain allowed for in silico separation of the interacting cells including quantification of polycistronic RNAs. Eighty-nine percent of the known loci are found to be transcribed in prokaryotic cells prior or subsequent to infection of the host, while 75% of all protein-coding loci are represented in the polyadenylated transcriptomes of human host cells. CONCLUSIONS: Dual 3'Seq was alternatively coupled to MACE (Massive Analysis of cDNA ends) to assess the advantages and drawbacks of a library preparation procedure that allows for sequencing of longer fragments. Additionally, the identified expression patterns of both organisms were validated by qRT-PCR using three independent biological replicates, which confirmed that RELB along with NFKB1 and NFKB2 are involved in the initial immune response of epithelial cells after infection with S. enterica Typhimurium.

Identification of novel small ncRNAs in pollen of tomato.

BACKGROUND: The unprecedented role of sncRNAs in the regulation of pollen biogenesis on both transcriptional and epigenetic levels has been experimentally proven. However, little is known about their global regulation, especially under stress conditions. We used tomato pollen in order to identify pollen stage-specific sncRNAs and their target mRNAs. We further deployed elevated temperatures to discern stress responsive sncRNAs. For this purpose high throughput sncRNA-sequencing as well as Massive Analysis of cDNA Ends (MACE) were performed for three-replicated sncRNAs libraries derived from tomato tetrad, post-meiotic, and mature pollen under control and heat stress conditions. RESULTS: Using the omiRas analysis pipeline we identified known and predicted novel miRNAs as well as sncRNAs from other classes, responsive or not to heat. Differential expression analysis revealed that post-meiotic and mature pollen react most strongly by regulation of the expression of coding and non-coding genomic regions in response to heat. To gain insight to the function of these miRNAs, we predicted targets and annotated them to Gene Ontology terms. This approach revealed that most of them belong to protein binding, transcription, and Serine/Threonine kinase activity GO categories. Beside miRNAs, we observed differential expression of both tRNAs and snoRNAs in tetrad, post-meiotic, and mature pollen when comparing normal and heat stress conditions. CONCLUSIONS: Thus, we describe a global spectrum of sncRNAs expressed in pollen as well as unveiled those which are regulated at specific time-points during pollen biogenesis. We integrated the small RNAs into the regulatory network of tomato heat stress response in pollen.

Next-generation sequencing reveals novel differentially regulated mRNAs, lncRNAs, miRNAs, sdRNAs and a piRNA in pancreatic cancer.

BACKGROUND: Previous studies identified microRNAs (miRNAs) and messenger RNAs with significantly different expression between normal pancreas and pancreatic cancer (PDAC) tissues. Due to technological limitations of microarrays and real-time PCR systems these studies focused on a fixed set of targets. Expression of other RNA classes such as long intergenic non-coding RNAs or sno-derived RNAs has rarely been examined in pancreatic cancer. Here, we analysed the coding and non-coding transcriptome of six PDAC and five control tissues using next-generation sequencing. RESULTS: Besides the confirmation of several deregulated mRNAs and miRNAs, miRNAs without previous implication in PDAC were detected: miR-802, miR-2114 or miR-561. SnoRNA-derived RNAs (e.g. sno-HBII-296B) and piR-017061, a piwi-interacting RNA, were found to be differentially expressed between PDAC and control tissues. In silico target analysis of miR-802 revealed potential binding sites in the 3' UTR of TCF4, encoding a transcription factor that controls Wnt signalling genes. Overexpression of miR-802 in MiaPaCa pancreatic cancer cells reduced TCF4 protein levels. Using Massive Analysis of cDNA Ends (MACE) we identified differential expression of 43 lincRNAs, long intergenic non-coding RNAs, e.g. LINC00261 and LINC00152 as well as several natural antisense transcripts like HNF1A-AS1 and AFAP1-AS1. Differential expression was confirmed by qPCR on the mRNA/miRNA/lincRNA level and by immunohistochemistry on the protein level. CONCLUSIONS: Here, we report a novel lncRNA, sncRNA and mRNA signature of PDAC. In silico prediction of ncRNA targets allowed for assigning potential functions to differentially regulated RNAs.

The Other Itch That Rashes: a Clinical and Therapeutic Approach to Pruritus and Skin Picking Disorders.

Patients with skin picking disorders (SPDs) have historically been an under recognized and under treated group. Originally classified an impulse control disorder, skin picking disorder is now considered under the obsessive-compulsive disorder (OCD) and obsessive-compulsive-related disorder (OCRD) spectrum. Workup of SPD concentrates on differentiating this condition from primary skin disorders, systemic conditions associated with pruritus, and psychocutaneous syndromes. It is important to first address any underlying pruritic disorders (if present) that may be a trigger to pick the skin. First-line treatment for SPD is cognitive behavioral therapy. When necessary, pharmacologic therapy can be helpful as well.

omiRas: a Web server for differential expression analysis of miRNAs derived from small RNA-Seq data.

SUMMARY: Small RNA deep sequencing is widely used to characterize non-coding RNAs (ncRNAs) differentially expressed between two conditions, e.g. healthy and diseased individuals and to reveal insights into molecular mechanisms underlying condition-specific phenotypic traits. The ncRNAome is composed of a multitude of RNAs, such as transfer RNA, small nucleolar RNA and microRNA (miRNA), to name few. Here we present omiRas, a Web server for the annotation, comparison and visualization of interaction networks of ncRNAs derived from next-generation sequencing experiments of two different conditions. The Web tool allows the user to submit raw sequencing data and results are presented as: (i) static annotation results including length distribution, mapping statistics, alignments and quantification tables for each library as well as lists of differentially expressed ncRNAs between conditions and (ii) an interactive network visualization of user-selected miRNAs and their target genes based on the combination of several miRNA-mRNA interaction databases. AVAILABILITY AND IMPLEMENTATION: The omiRas Web server is implemented in Python, PostgreSQL, R and can be accessed at: http://tools.genxpro.net/omiras/.

One Mutation to Rule Them All: Mutant KRAS Controls Tumor Intrinsic and Microenvironment Signaling.

In a recent study published in Cancer Discovery, Hsu and colleagues employ an elegant combination of single-cell and bulk RNA sequencing experiments from mouse and human colorectal cancer samples, patient-derived organoids, two-dimensional in vitro systems, and in vivo validation in genetically engineered colorectal cancer mouse models to investigate how mutant KRAS (KRAS*) impacts the tumor microenvironment. They identify a molecular signaling cascade downstream of KRAS* that activates a specific program of lipid-rich cancer-associated fibroblasts, promoting tumor angiogenesis and progression. These findings may lead to new therapeutic strategies for patients with colorectal cancer with KRAS*.

The Exonuclease TREX1 Constitutes an Innate Immune Checkpoint Limiting cGAS/STING-Mediated Antitumor Immunity.

The DNA exonuclease three-prime repair exonuclease 1 (TREX1) is critical for preventing autoimmunity in mice and humans by degrading endogenous cytosolic DNA, which otherwise triggers activation of the innate cGAS/STING pathway leading to the production of type I IFNs. As tumor cells are prone to aberrant cytosolic DNA accumulation, we hypothesized that they are critically dependent on TREX1 activity to limit their immunogenicity. Here, we show that in tumor cells, TREX1 restricts spontaneous activation of the cGAS/STING pathway, and the subsequent induction of a type I IFN response. As a result, TREX1 deficiency compromised in vivo tumor growth in mice. This delay in tumor growth depended on a functional immune system, systemic type I IFN signaling, and tumor-intrinsic cGAS expression. Mechanistically, we show that tumor TREX1 loss drove activation of CD8+ T cells and NK cells, prevented CD8+ T-cell exhaustion, and remodeled an immunosuppressive myeloid compartment. Consequently, TREX1 deficiency combined with T-cell-directed immune checkpoint blockade. Collectively, we conclude that TREX1 is essential to limit tumor immunogenicity, and that targeting this innate immune checkpoint remodels the tumor microenvironment and enhances antitumor immunity by itself and in combination with T-cell-targeted therapies. See related article by Toufektchan et al., p. 673.