Search

Objectives: Autologous chimeric antigen receptor (CAR) T-cell therapy of B-cell malignancies achieves long-term disease remission in a high fraction of patients and has triggered intense research into translating this successful approach into additional cancer types. However, the complex logistics involved in autologous CAR-T manufacturing, the compromised fitness of patient-derived T cells, the high rates of serious toxicities and the overall cost involved with product manufacturing and hospitalisation have driven innovation to overcome such hurdles. One alternative approach is the use of allogeneic natural killer (NK) cells as a source for CAR-NK cell therapy. However, this source has traditionally faced numerous manufacturing challenges., Methods: To address this, we have developed an optimised expansion and transduction protocol for primary human NK cells primed for manufacturing scaling and clinical evaluation. We have performed an in-depth comparison of primary human NK cell sources as a starting material by characterising their phenotype, functionality, expansion potential and transduction efficiency at crucial timepoints of our CAR-NK manufacturing pipeline., Results: We identified adult peripheral blood-derived NK cells to be the superior source for generating a CAR-NK cell product because of a higher maximum yield of CAR-expressing NK cells combined with potent natural, as well as CAR-mediated anti-tumor effector functions., Conclusions: Our optimised manufacturing pipeline dramatically improves lentiviral transduction efficiency of primary human NK cells. We conclude that the exponential expansion pre- and post-transduction and high on-target cytotoxicity make peripheral blood-derived NK cells a feasible and attractive CAR-NK cell product for clinical utility., Competing Interests: JC, RB and IN report employment with oNKo‐Innate. NDH, JR, IN, JC and RB report stock or other ownership in oNKo‐Innate. NDH serves on an advisory board for Bristol Myers Squibb and Syena. CKS and KSS report employment with Kite, a Gilead company, and stock or other ownership in Gilead Sciences. KSS serves on the advisor board of Obsidian Therapeutics., (© 2024 The Authors. Clinical & Translational Immunology published by John Wiley & Sons Australia, Ltd on behalf of Australian and New Zealand Society for Immunology, Inc.)

Ikaros transcription factors are essential for adaptive lymphocyte function, yet their role in innate lymphopoiesis is unknown. Using conditional genetic inactivation, we show that Ikzf1/Ikaros is essential for normal natural killer (NK) cell lymphopoiesis and IKZF1 directly represses Cish, a negative regulator of interleukin-15 receptor resulting in impaired interleukin-15 receptor signaling. Both Bcl2l11 and BIM levels, and intrinsic apoptosis were increased in Ikzf1-null NK cells, which in part accounts for NK lymphopenia as both were restored to normal levels when Ikzf1 and Bcl2l11 were co-deleted. Ikzf1-null NK cells presented extensive transcriptional alterations with reduced AP-1 transcriptional complex expression and increased expression of Ikzf2/Helios and Ikzf3/Aiolos. IKZF1 and IKZF3 directly bound AP-1 family members and deletion of both Ikzf1 and Ikzf3 in NK cells resulted in further reductions in Jun/Fos expression and complete loss of peripheral NK cells. Collectively, we show that Ikaros family members are important regulators of apoptosis, cytokine responsiveness and AP-1 transcriptional activity., (© 2024. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Drugs targeting cyclin-dependent kinases 4 and 6 (CDK4/6) are promising new treatments for melanoma and other solid malignancies. In studies on CDK4/6 inhibitor resistance, protein arginine methyltransferase 5 (PRMT5) regulation of alternative splicing was shown to be an important downstream component of the CDK4/6 pathway. However, the full effects of inhibition of CDK4/6 on splicing events in melanoma and the extent to which they are dependent on PRMT5 has not been established. We performed full-length mRNA sequencing on CHL1 and A375 melanoma cell lines treated with the CDK4/6 inhibitor palbociclib and the PRMT5 inhibitor GSK3326595 and analysed data for differential gene expression and differential pre-mRNA splicing induced by these agents. Changes in gene expression and RNA splicing were more extensive under PRMT5 inhibition than under CDK4/6 inhibition. Although PRMT5 inhibition and CDK4/6 inhibition induced common RNA splicing events and gene expression profiles, the majority of events induced by CDK4/6 inhibition were distinct. Our findings indicate CDK4/6 has the ability to regulate alternative splicing in a manner that is distinct from PRMT5 inhibition, resulting in divergent changes in gene expression under each therapy., Competing Interests: The authors declare that they have no competing interests., (Copyright: © 2023 Chan et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

KAT6A, and its paralog KAT6B, are histone lysine acetyltransferases (HAT) that acetylate histone H3K23 and exert an oncogenic role in several tumor types including breast cancer where KAT6A is frequently amplified/overexpressed. However, pharmacologic targeting of KAT6A to achieve therapeutic benefit has been a challenge. Here we describe identification of a highly potent, selective, and orally bioavailable KAT6A/KAT6B inhibitor CTx-648 (PF-9363), derived from a benzisoxazole series, which demonstrates anti-tumor activity in correlation with H3K23Ac inhibition in KAT6A over-expressing breast cancer. Transcriptional and epigenetic profiling studies show reduced RNA Pol II binding and downregulation of genes involved in estrogen signaling, cell cycle, Myc and stem cell pathways associated with CTx-648 anti-tumor activity in ER-positive (ER+) breast cancer. CTx-648 treatment leads to potent tumor growth inhibition in ER+ breast cancer in vivo models, including models refractory to endocrine therapy, highlighting the potential for targeting KAT6A in ER+ breast cancer., Competing Interests: Declaration of interests S.S., T.A.P., C.Y.C., S.U., J.P., J.C., A.R., N.N., S.G., E.J., O.B., S.K., H.W., Z.W., Y.Z., K.T., L.C., A.M., J.L., P.P.K., S.S., P.F.R., E.C.G., S.Y., R.J., K.A.M., P.B., and H.L. are/were employees of Pfizer and some of the authors are shareholders in Pfizer Inc. L.C. and A.M. are co-inventors on patent WO2022013369. P.S., H.R.L., B.J.M., R.F., C.H., M.C., Y.B., and S.W. are inventors on patent WO2019243491 and have received payments from Monash University derived from distribution of licensing income received from Pfizer. P.P.K., S.S., P.F.R., Y.B., and P.S. are inventors on patent WO2020254946. P.S. is an inventor on patent WO2022013369. T.T., H.R.L and A.K.V. are inventors on patent WO2016198507 A1. The Thomas and Voss laboratories received research funding from the Cooperative Research Centre (CRC) for Cancer Therapeutics (CTx), Australia. T.T. and A.K.V. have received payments from a distribution of the licensing income from Pfizer for serving on an advisory board meeting for Pfizer. M.A.D. has been a member of advisory boards for CTX CRC, GSK, Storm Therapeutics, Celgene and Cambridge Epigenetix. S.-J.D. has been a member of advisory boards for Adela. M.A.D. and S.-J.D. receive research funding from CTx CRC and from the Pfizer Emerging Science Fund. Ian Street is a member of the Scientific Advisory Board of Oncology One Pty Ltd and has received payments from distribution of licensing income from Pfizer. J.N. is currently an adjunct professor at the University of New South Wales, Kensington, NSW, Australia. M.J.D and D.D.B. are currently working at South Australian immunoGENomics Cancer Institute (SAiGENCI), Faculty of Health and Medical Sciences, The University of Adelaide. J.C. is a current employee of oNKo-innate Pty Ltd. J.L. is currently an employee and shareholder of Belharra Therapeutics., (Copyright © 2023 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Epithelial-mesenchymal transition (EMT) is a reversible transcriptional program invoked by cancer cells to drive cancer progression. Transcription factor ZEB1 is a master regulator of EMT, driving disease recurrence in poor-outcome triple negative breast cancers (TNBCs). Here, this work silences ZEB1 in TNBC models by CRISPR/dCas9-mediated epigenetic editing, resulting in highly-specific and nearly complete suppression of ZEB1 in vivo, accompanied by long-lasting tumor inhibition. Integrated "omic" changes promoted by dCas9 linked to the KRAB domain (dCas9-KRAB) enabled the discovery of a ZEB1-dependent-signature of 26 genes differentially-expressed and -methylated, including the reactivation and enhanced chromatin accessibility in cell adhesion loci, outlining epigenetic reprogramming toward a more epithelial state. In the ZEB1 locus transcriptional silencing is associated with induction of locally-spread heterochromatin, significant changes in DNA methylation at specific CpGs, gain of H3K9me3, and a near complete erasure of H3K4me3 in the ZEB1 promoter. Epigenetic shifts induced by ZEB1-silencing are enriched in a subset of human breast tumors, illuminating a clinically-relevant hybrid-like state. Thus, the synthetic epi-silencing of ZEB1 induces stable "lock-in" epigenetic reprogramming of mesenchymal tumors associated with a distinct and stable epigenetic landscape. This work outlines epigenome-engineering approaches for reversing EMT and customizable precision molecular oncology approaches for targeting poor outcome breast cancers., (© 2023 The Authors. Advanced Science published by Wiley-VCH GmbH.)

Background: Epigenetic silencing of tumor suppressor genes (TSGs) is a key feature of oncogenesis in hepatocellular carcinoma (HCC). Liver-targeted delivery of CRISPR-activation (CRISPRa) systems makes it possible to exploit chromatin plasticity, by reprogramming transcriptional dysregulation., Results: Using The Cancer Genome Atlas HCC data, we identify 12 putative TSGs with negative associations between promoter DNA methylation and transcript abundance, with limited genetic alterations. All HCC samples harbor at least one silenced TSG, suggesting that combining a specific panel of genomic targets could maximize efficacy, and potentially improve outcomes as a personalized treatment strategy for HCC patients. Unlike epigenetic modifying drugs lacking locus selectivity, CRISPRa systems enable potent and precise reactivation of at least 4 TSGs tailored to representative HCC lines. Concerted reactivation of HHIP, MT1M, PZP, and TTC36 in Hep3B cells inhibits multiple facets of HCC pathogenesis, such as cell viability, proliferation, and migration., Conclusions: By combining multiple effector domains, we demonstrate the utility of a CRISPRa toolbox of epigenetic effectors and gRNAs for patient-specific treatment of aggressive HCC., (© 2023. The Author(s).)

The advent of massively parallel sequencing revealed extensive transcription beyond protein-coding genes, identifying tens of thousands of long noncoding RNAs (lncRNAs). Selected functional examples raised the possibility that lncRNAs, as a class, may maintain broad regulatory roles. Expression of lncRNAs is strongly linked with adjacent protein-coding gene expression, suggesting potential cis -regulatory functions. A more detailed understanding of these regulatory roles may be obtained through careful examination of the precise timing of lncRNA expression relative to adjacent protein-coding genes. Despite the diversity of reported lncRNA regulatory mechanisms, where causal cis -regulatory relationships exist, lncRNA transcription is expected to precede changes in target gene expression. Using a high temporal resolution RNA-seq time course, we profiled the expression dynamics of several thousand lncRNAs and protein-coding genes in synchronized, transitioning human cells. Our findings reveal that lncRNAs are expressed synchronously with adjacent protein-coding genes. Analysis of lipopolysaccharide-activated mouse dendritic cells revealed the same temporal relationship observed in transitioning human cells. Our findings suggest broad-scale cis -regulatory roles for lncRNAs are not common. The strong association between lncRNAs and adjacent genes may instead indicate an origin as transcriptional by-products from active protein-coding gene promoters and enhancers., (© 2022 Muskovic et al.; Published by Cold Spring Harbor Laboratory Press.)

The development of therapies that target specific disease subtypes has dramatically improved outcomes for patients with breast cancer. However, survival gains have not been uniform across patients, even within a given molecular subtype. Large collections of publicly available drug screening data matched with transcriptomic measurements have facilitated the development of computational models that predict response to therapy. Here, we generated a series of predictive gene signatures to estimate the sensitivity of breast cancer samples to 90 drugs, comprising FDA-approved drugs or compounds in early development. To achieve this, we used a cell line-based drug screen with matched transcriptomic data to derive in silico models that we validated in large independent datasets obtained from cell lines and patient-derived xenograft (PDX) models. Robust computational signatures were obtained for 28 drugs and used to predict drug efficacy in a set of PDX models. We found that our signature for cisplatin can be used to identify tumors that are likely to respond to this drug, even in absence of the BRCA-1 mutation routinely used to select patients for platinum-based therapies. This clinically relevant observation was confirmed in multiple PDXs. Our study foreshadows an effective delivery approach for precision medicine.

It is widely acknowledged that the construction of large-scale dynamic models in systems biology requires complex modelling problems to be broken up into more manageable pieces. To this end, both modelling and software frameworks are required to enable modular modelling. While there has been consistent progress in the development of software tools to enhance model reusability, there has been a relative lack of consideration for how underlying biophysical principles can be applied to this space. Bond graphs combine the aspects of both modularity and physics-based modelling. In this paper, we argue that bond graphs are compatible with recent developments in modularity and abstraction in systems biology, and are thus a desirable framework for constructing large-scale models. We use two examples to illustrate the utility of bond graphs in this context: a model of a mitogen-activated protein kinase (MAPK) cascade to illustrate the reusability of modules and a model of glycolysis to illustrate the ability to modify the model granularity., Competing Interests: The authors have declared that no competing interests exist. Author Edmund Crampin was unable to confirm their authorship contributions. On their behalf, the corresponding author has reported their contributions to the best of their knowledge.

Immunotherapy success in colorectal cancer is mainly limited to patients whose tumors exhibit high microsatellite instability (MSI). However, there is variability in treatment outcomes within this group, which is in part driven by the frequency and characteristics of tumor-infiltrating immune cells. Indeed, the presence of specific infiltrating immune-cell subsets has been shown to correlate with immunotherapy response and is in many cases prognostic of treatment outcome. Tumor-infiltrating lymphocytes (TIL) can undergo distinct differentiation programs, acquiring features of tissue-residency or exhaustion, a process during which T cells upregulate inhibitory receptors, such as PD-1, and lose functionality. Although residency and exhaustion programs of CD8 + T cells are relatively well studied, these programs have only recently been appreciated in CD4 + T cells and remain largely unknown in tumor-infiltrating natural killer (NK) cells. In this study, we used single-cell RNA sequencing (RNA-seq) data to identify signatures of residency and exhaustion in colorectal cancer-infiltrating lymphocytes, including CD8 + , CD4 + , and NK cells. We then tested these signatures in independent single-cell data from tumor and normal tissue-infiltrating immune cells. Furthermore, we used versions of these signatures designed for bulk RNA-seq data to explore tumor-intrinsic mutations associated with residency and exhaustion from TCGA data. Finally, using two independent transcriptomic datasets from patients with colon adenocarcinoma, we showed that combinations of these signatures, in particular combinations of NK-cell activity signatures, together with tumor-associated signatures, such as TGFβ signaling, were associated with distinct survival outcomes in patients with colon adenocarcinoma., (©2021 American Association for Cancer Research.)

Adipogenesis associated Mth938 domain containing (AAMDC) represents an uncharacterized oncogene amplified in aggressive estrogen receptor-positive breast cancers. We uncover that AAMDC regulates the expression of several metabolic enzymes involved in the one-carbon folate and methionine cycles, and lipid metabolism. We show that AAMDC controls PI3K-AKT-mTOR signaling, regulating the translation of ATF4 and MYC and modulating the transcriptional activity of AAMDC-dependent promoters. High AAMDC expression is associated with sensitization to dactolisib and everolimus, and these PI3K-mTOR inhibitors exhibit synergistic interactions with anti-estrogens in IntClust2 models. Ectopic AAMDC expression is sufficient to activate AKT signaling, resulting in estrogen-independent tumor growth. Thus, AAMDC-overexpressing tumors may be sensitive to PI3K-mTORC1 blockers in combination with anti-estrogens. Lastly, we provide evidence that AAMDC can interact with the RabGTPase-activating protein RabGAP1L, and that AAMDC, RabGAP1L, and Rab7a colocalize in endolysosomes. The discovery of the RabGAP1L-AAMDC assembly platform provides insights for the design of selective blockers to target malignancies having the AAMDC amplification.

Gene expression signatures have been critical in defining the molecular phenotypes of cells, tissues, and patient samples. Their most notable and widespread clinical application is stratification of breast cancer patients into molecular (PAM50) subtypes. The cost and relatively large amounts of fresh starting material required for whole-transcriptome sequencing has limited clinical application of thousands of existing gene signatures captured in repositories such as the Molecular Signature Database. We identified genes with stable expression across a range of abundances, and with a preserved relative ordering across thousands of samples, allowing signature scoring and supporting general data normalisation for transcriptomic data. Our new method, stingscore, quantifies and summarises relative expression levels of signature genes from individual samples through the inclusion of these 'stably-expressed genes'. We show that our list of stable genes has better stability across cancer and normal tissue data than previously proposed gene sets. Additionally, we show that signature scores computed from targeted transcript measurements using stingscore can predict docetaxel response in breast cancer patients. This new approach to gene expression signature analysis will facilitate the development of panel-type tests for gene expression signatures, thus supporting clinical translation of the powerful insights gained from cancer transcriptomic studies., (© The Author(s) 2020. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Introduction: SCLC is the most aggressive subtype of lung cancer, and though most patients initially respond to platinum-based chemotherapy, resistance develops rapidly. Immunotherapy holds promise in the treatment of lung cancer; however, patients with SCLC exhibit poor overall responses highlighting the necessity for alternative approaches. Natural killer (NK) cells are an alternative to T cell-based immunotherapies that do not require sensitization to antigens presented on the surface of tumor cells., Methods: We investigated the immunophenotype of human SCLC tumors by both flow cytometry on fresh samples and bioinformatic analysis. Cell lines generated from murine SCLC were transplanted into mice lacking key cytotoxic immune cells. Subcutaneous tumor growth, metastatic dissemination, and activation of CD8 + T and NK cells were evaluated by histology and flow cytometry., Results: Transcriptomic analysis of human SCLC tumors revealed heterogeneous immune checkpoint and cytotoxic signature profiles. Using sophisticated, genetically engineered mouse models, we reported that the absence of NK cells, but not CD8 + T cells, substantially enhanced metastatic dissemination of SCLC tumor cells in vivo. Moreover, hyperactivation of NK cell activity through augmentation of interleukin-15 or transforming growth factor-β signaling pathways ameliorated SCLC metastases, an effect that was enhanced when combined with antiprogrammed cell death-1 therapy., Conclusions: These proof-of-principle findings provide a rationale for exploiting the antitumor functions of NK cells in the treatment of patients with SCLC. Moreover, the distinct immune profiles of SCLC subtypes reveal an unappreciated level of heterogeneity that warrants further investigation in the stratification of patients for immunotherapy., (Copyright © 2020 International Association for the Study of Lung Cancer. Published by Elsevier Inc. All rights reserved.)

Immunotherapy with checkpoint blockade induces rapid and durable immune control of cancer in some patients and has driven a monumental shift in cancer treatment. Neoantigen-specific CD8 + T cells are at the forefront of current immunotherapy strategies, and the majority of drug discovery and clinical trials revolve around further harnessing these immune effectors. Yet the immune system contains a diverse range of antitumour effector cells, and these must function in a coordinated and synergistic manner to overcome the immune-evasion mechanisms used by tumours and achieve complete control with tumour eradication. A key antitumour effector is the natural killer (NK) cells, cytotoxic innate lymphocytes present at high frequency in the circulatory system and identified by their exquisite ability to spontaneously detect and lyse transformed or stressed cells. Emerging data show a role for intratumoural NK cells in driving immunotherapy response and, accordingly, there have been renewed efforts to further elucidate and target the pathways controlling NK cell antitumour function. In this Review, we discuss recent clinical evidence that NK cells are a key immune constituent in the protective antitumour immune response and highlight the major stages of the cancer-NK cell immunity cycle. We also perform a new analysis of publicly available transcriptomic data to provide an overview of the prognostic value of NK cell gene expression in 25 tumour types. Furthermore, we discuss how the role of NK cells evolves with tumour progression, presenting new opportunities to target NK cell function to enhance cancer immunotherapy response rates across a more diverse range of cancers.

Membrane transporters contribute to the regulation of the internal environment of cells by translocating substrates across cell membranes. Like all physical systems, the behaviour of membrane transporters is constrained by the laws of thermodynamics. However, many mathematical models of transporters, especially those incorporated into whole-cell models, are not thermodynamically consistent, leading to unrealistic behaviour. In this paper we use a physics-based modelling framework, in which the transfer of energy is explicitly accounted for, to develop thermodynamically consistent models of transporters. We then apply this methodology to model two specific transporters: the cardiac sarcoplasmic/endoplasmic Ca 2+ ATPase (SERCA) and the cardiac Na + /K + ATPase., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Background: Elucidation of regulatory networks, including identification of regulatory mechanisms specific to a given biological context, is a key aim in systems biology. This has motivated the move from co-expression to differential co-expression analysis and numerous methods have been developed subsequently to address this task; however, evaluation of methods and interpretation of the resulting networks has been hindered by the lack of known context-specific regulatory interactions., Results: In this study, we develop a simulator based on dynamical systems modelling capable of simulating differential co-expression patterns. With the simulator and an evaluation framework, we benchmark and characterise the performance of inference methods. Defining three different levels of "true" networks for each simulation, we show that accurate inference of causation is difficult for all methods, compared to inference of associations. We show that a z-score-based method has the best general performance. Further, analysis of simulation parameters reveals five network and simulation properties that explained the performance of methods. The evaluation framework and inference methods used in this study are available in the dcanr R/Bioconductor package., Conclusions: Our analysis of networks inferred from simulated data show that hub nodes are more likely to be differentially regulated targets than transcription factors. Based on this observation, we propose an interpretation of the inferred differential network that can reconstruct a putative causal network.

Epithelial-mesenchymal transition (EMT) has been a subject of intense scrutiny as it facilitates metastasis and alters drug sensitivity. Although EMT-regulatory roles for numerous miRNAs and transcription factors are known, their functions can be difficult to disentangle, in part due to the difficulty in identifying direct miRNA targets from complex datasets and in deciding how to incorporate 'indirect' miRNA effects that may, or may not, represent biologically relevant information. To better understand how miRNAs exert effects throughout the transcriptome during EMT, we employed Exon-Intron Split Analysis (EISA), a bioinformatic technique that separates transcriptional and post-transcriptional effects through the separate analysis of RNA-Seq reads mapping to exons and introns. We find that in response to the manipulation of miRNAs, a major effect on gene expression is transcriptional. We also find extensive co-ordination of transcriptional and post-transcriptional regulatory mechanisms during both EMT and mesenchymal to epithelial transition (MET) in response to TGF-β or miR-200c respectively. The prominent transcriptional influence of miRNAs was also observed in other datasets where miRNA levels were perturbed. This work cautions against a narrow approach that is limited to the analysis of direct targets, and demonstrates the utility of EISA to examine complex regulatory networks involving both transcriptional and post-transcriptional mechanisms., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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., (Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

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., (©2019 American Association for Cancer Research.)

Advances in RNA sequencing (RNA-seq) technologies that measure the transcriptome of biological samples have revolutionised our ability to understand transcriptional regulatory programs that underpin diseases such as cancer. We recently published singscore - a single sample, rank-based gene set scoring method which quantifies how concordant the transcriptional profile of individual samples are relative to specific gene sets of interest. Here we demonstrate the application of singscore to investigate transcriptional profiles associated with specific mutations or genetic lesions in acute myeloid leukemia. Using matched genomic and transcriptomic data available through the TCGA we show that scoring of appropriate signatures can distinguish samples with corresponding mutations, reflecting the ability of these mutations to drive aberrant transcriptional programs involved in leukemogenesis. We believe the singscore method is particularly useful for studying heterogeneity within a specific subsets of cancers, and as demonstrated, we show the ability of singscore to identify where alternative mutations appear to drive similar transcriptional programs., Competing Interests: No competing interests were disclosed., (Copyright: © 2019 Bhuva DD et al.)

Advances in RNA sequencing (RNA-seq) technologies that measure the transcriptome of biological samples have revolutionised our ability to understand transcriptional regulatory programs that underpin diseases such as cancer. We recently published singscore - a single sample, rank-based gene set scoring method which quantifies how concordant the transcriptional profile of individual samples are relative to specific gene sets of interest. Here we demonstrate the application of singscore to investigate transcriptional profiles associated with specific mutations or genetic lesions in acute myeloid leukemia. Using matched genomic and transcriptomic data available through the TCGA we show that scoring of appropriate signatures can distinguish samples with corresponding mutations, reflecting the ability of these mutations to drive aberrant transcriptional programs involved in leukemogenesis. We believe the singscore method is particularly useful for studying heterogeneity within a specific subsets of cancers, and as demonstrated, we show the ability of singscore to identify where alternative mutations appear to drive similar transcriptional programs., Competing Interests: No competing interests were disclosed., (Copyright: © 2019 Bhuva DD et al.)

Immune 'checkpoint' inhibitors can increase the activity of tumor-resident cytotoxic lymphocytes and have revolutionized cancer treatment. Current therapies block inhibitory pathways in tumor-infiltrating CD8 + T cells and recent studies have shown similar programs in other effector populations such as natural killer (NK) cells. NK cells are critical for immunosurveillance, particularly the control of metastatic cells or hematological cancers. However, how NK cells specifically recognize transformed cells and dominant negative feedback pathways, as well as how tumors escape NK cell control, remains undefined. This review summarizes recent advances that have illuminated inhibitory checkpoints in NK cells, some of which are shared with conventional cytotoxic T lymphocytes. It also outlines emerging approaches aimed at unleashing the potential of NK cells in immunotherapy., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Background: Gene set scoring provides a useful approach for quantifying concordance between sample transcriptomes and selected molecular signatures. Most methods use information from all samples to score an individual sample, leading to unstable scores in small data sets and introducing biases from sample composition (e.g. varying numbers of samples for different cancer subtypes). To address these issues, we have developed a truly single sample scoring method, and associated R/Bioconductor package singscore ( https://bioconductor.org/packages/singscore )., Results: We use multiple cancer data sets to compare singscore against widely-used methods, including GSVA, z-score, PLAGE, and ssGSEA. Our approach does not depend upon background samples and scores are thus stable regardless of the composition and number of samples being scored. In contrast, scores obtained by GSVA, z-score, PLAGE and ssGSEA can be unstable when less data are available (N S  < 25). The singscore method performs as well as the best performing methods in terms of power, recall, false positive rate and computational time, and provides consistently high and balanced performance across all these criteria. To enhance the impact and utility of our method, we have also included a set of functions implementing visual analysis and diagnostics to support the exploration of molecular phenotypes in single samples and across populations of data., Conclusions: The singscore method described here functions independent of sample composition in gene expression data and thus it provides stable scores, which are particularly useful for small data sets or data integration. Singscore performs well across all performance criteria, and includes a suite of powerful visualization functions to assist in the interpretation of results. This method performs as well as or better than other scoring approaches in terms of its power to distinguish samples with distinct biology and its ability to call true differential gene sets between two conditions. These scores can be used for dimensional reduction of transcriptomic data and the phenotypic landscapes obtained by scoring samples against multiple molecular signatures may provide insights for sample stratification.

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., (Copyright © 2018 The Authors. Published by Elsevier Inc. All rights reserved.)