Natural killer (NK) cell activity is essential for initiating antitumor responses and may be linked to immunotherapy success. NK cells and other innate immune components could be exploitable for cancer treatment, which drives the need for tools and methods that identify therapeutic avenues. Here, we extend our gene-set scoring method singscore to investigate NK cell infiltration by applying RNA-seq analysis to samples from bulk tumors. Computational methods have been developed for the deconvolution of immune cell types within solid tumors. We have taken the NK cell gene signatures from several such tools, then curated the gene list using a comparative analysis of tumors and immune cell types. Using a gene-set scoring method to investigate RNA-seq data from The Cancer Genome Atlas (TCGA), we show that patients with metastatic cutaneous melanoma have an improved survival rate if their tumor shows evidence of NK cell infiltration. Furthermore, these survival effects are enhanced in tumors that show higher expression of genes that encode NK cell stimuli such as the cytokine IL15. Using this signature, we then examine transcriptomic data to identify tumor and stromal components that may influence the penetrance of NK cells into solid tumors. Our results provide evidence that NK cells play a role in the regulation of human tumors and highlight potential survival effects associated with increased NK cell activity. Our computational analysis identifies putative gene targets that may be of therapeutic value for boosting NK cell antitumor immunity.

Recent advances in single-cell technologies have enabled high-throughput molecular profiling of cells across modalities and locations. Single-cell transcriptomics data can now be complemented by chromatin accessibility, surface protein expression, adaptive immune receptor repertoire profiling and spatial information. The increasing availability of single-cell data across modalities has motivated the development of novel computational methods to help analysts derive biological insights. As the field grows, it becomes increasingly difficult to navigate the vast landscape of tools and analysis steps. Here, we summarize independent benchmarking studies of unimodal and multimodal single-cell analysis across modalities to suggest comprehensive best-practice workflows for the most common analysis steps. Where independent benchmarks are not available, we review and contrast popular methods. Our article serves as an entry point for novices in the field of single-cell (multi-)omic analysis and guides advanced users to the most recent best practices.

Group 2 innate lymphoid cells (ILC2) are essential to maintain tissue homeostasis. In cancer, ILC2 can harbor both pro- and anti-tumorigenic functions but we know very little about their underlying mechanisms, nor whether they could be clinically relevant or targeted to improve patient outcomes. Here, we found that high ILC2 infiltration in human melanoma was associated with a good clinical prognosis. ILC2 are critical producers of the cytokine granulocyte-macrophage colony-stimulating factor (GM-CSF) which coordinate the recruitment and activation of eosinophils to enhance anti-tumor responses. Tumor-infiltrating ILC2 expressed programmed cell death protein-1 (PD-1), which limited their intratumoral accumulation, proliferation and anti-tumor effector functions. This inhibition could be overcome in vivo by combining IL-33-driven ILC2 activation with PD-1 blockade to significantly increase anti-tumor responses. Together, our results identified ILC2 as a critical immune cell type involved in melanoma immunity and revealed a potential synergistic approach to harness ILC2 function for anti-tumor immunotherapies.

MicroRNAs (miRNAs) are important post-transcriptional regulators of gene expression, functioning in part by facilitating the degradation of target mRNAs. They have an established role in controlling epithelial-mesenchymal transition (EMT), a reversible phenotypic program underlying normal and pathological processes. Many studies demonstrate the role of individual miRNAs using overexpression at levels greatly exceeding physiological abundance. This can influence transcripts with relatively poor targeting and may in part explain why over 130 different miRNAs are directly implicated as EMT regulators. Analyzing a human mammary cell model of EMT we found evidence that a set of miRNAs, including the miR-200 and miR-182/183 family members, co-operate in post-transcriptional regulation, both reinforcing and buffering transcriptional output. Investigating this, we demonstrate that combinatorial treatment altered cellular phenotype with miRNA concentrations much closer to endogenous levels and with less off-target effects. This suggests that co-operative targeting by miRNAs is important for their physiological function and future work classifying miRNAs should consider such combinatorial effects.

Most cancer deaths are due to metastasis, and epithelial-to-mesenchymal transition (EMT) plays a central role in driving cancer cell metastasis. EMT is induced by different stimuli, leading to different signaling patterns and therapeutic responses. TGFβ is one of the best-studied drivers of EMT, and many drugs are available to target this signaling pathway. A comprehensive bioinformatics approach was employed to derive a signature for TGFβ-induced EMT which can be used to score TGFβ-driven EMT in cells and clinical specimens. Considering this signature in pan-cancer cell and tumor datasets, a number of cell lines (including basal B breast cancer and cancers of the central nervous system) show evidence for TGFβ-driven EMT and carry a low mutational burden across the TGFβ signaling pathway. Furthermore, significant variation is observed in the response of high scoring cell lines to some common cancer drugs. Finally, this signature was applied to pan-cancer data from The Cancer Genome Atlas to identify tumor types with evidence of TGFβ-induced EMT. Tumor types with high scores showed significantly lower survival rates than those with low scores and also carry a lower mutational burden in the TGFβ pathway. The current transcriptomic signature demonstrates reproducible results across independent cell line and cancer datasets and identifies samples with strong mesenchymal phenotypes likely to be driven by TGFβ. The TGFβ-induced EMT signature may be useful to identify patients with mesenchymal-like tumors who could benefit from targeted therapeutics to inhibit promesenchymal TGFβ signaling and disrupt the metastatic cascade. .

Despite increasing recognition of the importance of GM-CSF in autoimmune disease, it remains unclear how GM-CSF is regulated at sites of tissue inflammation. Using GM-CSF fate reporter mice, we show that synovial NK cells produce GM-CSF in autoantibody-mediated inflammatory arthritis. Synovial NK cells promote a neutrophilic inflammatory cell infiltrate, and persistent arthritis, via GM-CSF production, as deletion of NK cells, or specific ablation of GM-CSF production in NK cells, abrogated disease. Synovial NK cell production of GM-CSF is IL-18–dependent. Furthermore, we show that cytokine-inducible SH2-containing protein (CIS) is crucial in limiting GM-CSF signaling not only during inflammatory arthritis but also in experimental allergic encephalomyelitis (EAE), a murine model of multiple sclerosis. Thus, a cellular cascade of synovial macrophages, NK cells, and neutrophils mediates persistent joint inflammation via production of IL-18 and GM-CSF. Endogenous CIS provides a key brake on signaling through the GM-CSF receptor. These findings shed new light on GM-CSF biology in sterile tissue inflammation and identify several potential therapeutic targets.

Interleukin (IL)-17–producing CD8+ T (Tc17) cells have emerged as key players in host-microbiota interactions, infection, and cancer. The factors that drive their development, in contrast to interferon (IFN)-γ–producing effector CD8+ T cells, are not clear. Here we demonstrate that the transcription factor TCF-1 (Tcf7) regulates CD8+ T cell fate decisions in double-positive (DP) thymocytes through the sequential suppression of MAF and RORγt, in parallel with TCF-1–driven modulation of chromatin state. Ablation of TCF-1 resulted in enhanced Tc17 cell development and exposed a gene set signature to drive tissue repair and lipid metabolism, which was distinct from other CD8+ T cell subsets. IL-17–producing CD8+ T cells isolated from healthy humans were also distinct from CD8+IL-17− T cells and enriched in pathways driven by MAF and RORγt. Overall, our study reveals how TCF-1 exerts central control of T cell differentiation in the thymus by normally repressing Tc17 differentiation and promoting an effector fate outcome.

Natural killer (NK) cells are innate lymphocytes that play a major role in immunosurveillance against tumor initiation and metastatic spread. The signals and checkpoints that regulate NK cell fitness and function in the tumor microenvironment are not well defined. Transforming growth factor–β (TGF-β) is a suppressor of NK cells that inhibits interleukin-15 (IL-15)–dependent signaling events and increases the abundance of receptors that promote tissue residency. Here, we showed that NK cells express the type I activin receptor ALK4, which, upon binding to its ligand activin-A, phosphorylated SMAD2/3 to suppress IL-15–mediated NK cell metabolism. Activin-A impaired human and mouse NK cell proliferation and reduced the production of granzyme B to impair tumor killing. Similar to TGF-β, activin-A also induced SMAD2/3 phosphorylation and stimulated NK cells to increase their cell surface expression of several markers of ILC1 cells. Activin-A also induced these changes in TGF-β receptor–deficient NK cells, suggesting that activin-A and TGF-β stimulate independent pathways that drive SMAD2/3-mediated NK cell suppression. Last, inhibition of activin-A by follistatin substantially slowed orthotopic melanoma growth in mice. These data highlight the relevance of examining TGF-β–independent SMAD2/3 signaling mechanisms as a therapeutic axis to relieve NK cell suppression and promote antitumor immunity.