Network analysis reveals a major role for 14q32 cluster miRNAs in determining transcriptional differences between IGHV-mutated and unmutated CLL.

Chronic lymphocytic leukaemia (CLL) cells can express unmutated (U-CLL) or mutated (M-CLL) immunoglobulin heavy chain (IGHV) genes with differing clinical behaviours, variable B cell receptor (BCR) signalling capacity and distinct transcriptional profiles. As it remains unclear how these differences reflect the tumour cells' innate pre/post germinal centre origin or their BCR signalling competence, we applied mRNA/miRNA sequencing to 38 CLL cases categorised into three subsets by IGHV mutational status and BCR signalling capacity. We identified 492 mRNAs and 38 miRNAs differentially expressed between U-CLL and M-CLL, but only 9 mRNAs and 0 miRNAs associated with BCR competence within M-CLL. Of the IGHV-associated miRNAs, (14/38 (37%)) derived from chr14q32 clusters where all miRNAs were co-expressed with the MEG3 lncRNA from a cancer associated imprinted locus. Integrative analysis of miRNA/mRNA data revealed pronounced regulatory potential for the 14q32 miRNAs, potentially accounting for up to 25% of the IGHV-related transcriptome signature. GAB1, a positive regulator of BCR signalling, was potentially regulated by five 14q32 miRNAs and we confirmed that two of these (miR-409-3p and miR-411-3p) significantly repressed activity of the GAB1 3'UTR. Our analysis demonstrates a potential key role of the 14q32 miRNA locus in the regulation of CLL-related gene regulation.

Evaluating performance of existing computational models in predicting CD8+ T cell pathogenic epitopes and cancer neoantigens.

T cell recognition of a cognate peptide-major histocompatibility complex (pMHC) presented on the surface of infected or malignant cells is of the utmost importance for mediating robust and long-term immune responses. Accurate predictions of cognate pMHC targets for T cell receptors would greatly facilitate identification of vaccine targets for both pathogenic diseases and personalized cancer immunotherapies. Predicting immunogenic peptides therefore has been at the center of intensive research for the past decades but has proven challenging. Although numerous models have been proposed, performance of these models has not been systematically evaluated and their success rate in predicting epitopes in the context of human pathology has not been measured and compared. In this study, we evaluated the performance of several publicly available models, in identifying immunogenic CD8+ T cell targets in the context of pathogens and cancers. We found that for predicting immunogenic peptides from an emerging virus such as severe acute respiratory syndrome coronavirus 2, none of the models perform substantially better than random or offer considerable improvement beyond HLA ligand prediction. We also observed suboptimal performance for predicting cancer neoantigens. Through investigation of potential factors associated with ill performance of models, we highlight several data- and model-associated issues. In particular, we observed that cross-HLA variation in the distribution of immunogenic and non-immunogenic peptides in the training data of the models seems to substantially confound the predictions. We additionally compared key parameters associated with immunogenicity between pathogenic peptides and cancer neoantigens and observed evidence for differences in the thresholds of binding affinity and stability, which suggested the need to modulate different features in identifying immunogenic pathogen versus cancer peptides. Overall, we demonstrate that accurate and reliable predictions of immunogenic CD8+ T cell targets remain unsolved; thus, we hope our work will guide users and model developers regarding potential pitfalls and unsettled questions in existing immunogenicity predictors.

HLA-dependent variation in SARS-CoV-2 CD8 + T cell cross-reactivity with human coronaviruses.

The conditions and extent of cross-protective immunity between severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and common-cold human coronaviruses (HCoVs) remain open despite several reports of pre-existing T cell immunity to SARS-CoV-2 in individuals without prior exposure. Using a pool of functionally evaluated SARS-CoV-2 peptides, we report a map of 126 immunogenic peptides with high similarity to 285 MHC-presented peptides from at least one HCoV. Employing this map of SARS-CoV-2-non-homologous and homologous immunogenic peptides, we observe several immunogenic peptides with high similarity to human proteins, some of which have been reported to have elevated expression in severe COVID-19 patients. After combining our map with SARS-CoV-2-specific TCR repertoire data from COVID-19 patients and healthy controls, we show that public repertoires for the majority of convalescent patients are dominated by TCRs cognate to non-homologous SARS-CoV-2 peptides. We find that for a subset of patients, >50% of their public SARS-CoV-2-specific repertoires consist of TCRs cognate to homologous SARS-CoV-2-HCoV peptides. Further analysis suggests that this skewed distribution of TCRs cognate to homologous or non-homologous peptides in COVID-19 patients is likely to be HLA-dependent. Finally, we provide 10 SARS-CoV-2 peptides with known cognate TCRs that are conserved across multiple coronaviruses and are predicted to be recognized by a high proportion of the global population. These findings may have important implications for COVID-19 heterogeneity, vaccine-induced immune responses, and robustness of immunity to SARS-CoV-2 and its variants.

CTEN Induces Tumour Cell Invasion and Survival and Is Prognostic in Radiotherapy-Treated Head and Neck Cancer.

Head and neck squamous cell carcinoma (HNSCC) is a heterogenous disease treated with surgery and/or (chemo) radiotherapy, but up to 50% of patients with late-stage disease develop locoregional recurrence. Determining the mechanisms underpinning treatment resistance could identify new therapeutic targets and aid treatment selection. C-terminal tensin-like (CTEN) is a member of the tensin family, upregulated in several cancers, although its expression and function in HNSCC are unknown. We found that CTEN is commonly upregulated in HNSCC, particularly HPV-ve tumours. In vitro CTEN was upregulated in HPV-ve (n = 5) and HPV+ve (n = 2) HNSCC cell lines. Stable shRNA knockdown of CTEN in vivo significantly reduced tumour growth (SCC-25), and functional analyses in vitro showed that CTEN promoted tumour cell invasion, colony formation and growth in 3D-culture (SCC-25, Detroit 562). RNA sequencing of SCC-25 cells following CTEN siRNA knockdown identified 349 differentially expressed genes (logFC > 1, p < 0.05). Gene ontology analysis highlighted terms relating to cell locomotion and apoptosis, consistent with in vitro findings. A membrane-based antibody array confirmed that CTEN regulated multiple apoptosis-associated proteins, including HSP60 and cleaved caspase-3. Notably, in a mixed cohort of HPV+ve and HPV-ve HNSCC patients (n = 259), we found a significant, independent negative association of CTEN with prognosis, limited to those patients treated with (chemo)radiotherapy, not surgery, irrespective of human papillomavirus (HPV) status. These data show that CTEN is commonly upregulated in HNSCC and exerts several functional effects. Its potential role in modulating apoptotic response to therapy suggests utility as a predictive biomarker or radio-sensitising target.

Genomic programming of IRF4-expressing human Langerhans cells.

Langerhans cells (LC) can prime tolerogenic as well as immunogenic responses in skin, but the genomic states and transcription factors (TF) regulating these context-specific responses are unclear. Bulk and single-cell transcriptional profiling demonstrates that human migratory LCs are robustly programmed for MHC-I and MHC-II antigen presentation. Chromatin analysis reveals enrichment of ETS-IRF and AP1-IRF composite regulatory elements in antigen-presentation genes, coinciding with expression of the TFs, PU.1, IRF4 and BATF3 but not IRF8. Migration of LCs from the epidermis is accompanied by upregulation of IRF4, antigen processing components and co-stimulatory molecules. TNF stimulation augments LC cross-presentation while attenuating IRF4 expression. CRISPR-mediated editing reveals IRF4 to positively regulate the LC activation programme, but repress NF2EL2 and NF-kB pathway genes that promote responsiveness to oxidative stress and inflammatory cytokines. Thus, IRF4-dependent genomic programming of human migratory LCs appears to enable LC maturation while attenuating excessive inflammatory and immunogenic responses in the epidermis.

Immune subtyping of extranodal NK/T-cell lymphoma: a new biomarker and an immune shift during disease progression.

Extranodal NK/T-cell lymphoma is an aggressive lymphoma that is strongly associated with Epstein-Barr virus infection. Although some extranodal NK/T-cell lymphoma patients have shown responses to immune checkpoint blockade, biomarkers for predicting extranodal NK/T-cell lymphoma patient response to immunotherapy have not yet been defined. To understand the tumor immune microenvironment, we analyzed the expression of 579 immune-related genes and characterized the immune cells using immunohistochemistries and in situ hybridization for EBER. Based on comprehensive analyses, we developed an immune subtyping model that classifies extranodal NK/T-cell lymphoma patients into four tumor immune microenvironment subgroups using three immunohistochemical markers (FoxP3, PD-L1, and CD68). The four tumor immune microenvironment subgroups were named immune tolerance, immune evasion-A, immune evasion-B, and immune silenced. The immune tolerance group was characterized by high-Treg counts and was frequently observed in early stage, and nasal extranodal NK/T-cell lymphoma. The immune evasion group showed high cytotoxic T-cell counts and high PD-L1 expression but low Treg counts. In the immune-silenced group, almost all immune responses were exhausted, most patients were at an advanced stage, and had the poorest disease prognosis among the tumor immune microenvironment subgroups. In some patients (n = 3), a shift in the tumor immune microenvironment subgroup classification was observed in sequential biopsies. The response rate to pembrolizumab, an anti-PD-1 antibody, was 100% (1/1) in the immune tolerance group, 60% (3/5) in the immune evasion group, and 0% (0/5) in the immune-silenced group. We classified extranodal NK/T-cell lymphoma into four tumor immune microenvironment subgroups using a new classification system. In conclusion, we propose that the tumor immune microenvironment of extranodal NK/T-cell lymphoma may change during disease progression and may serve as a useful biomarker for immunotherapy.

Preferential Infiltration of Unique Vγ9Jγ2-Vδ2 T Cells Into Glioblastoma Multiforme.

Glioblastoma multiforme (GBM) is clinically highly aggressive as a result of evolutionary dynamics induced by cross-talk between cancer cells and a heterogeneous group of immune cells in tumor microenvironment. The brain harbors limited numbers of immune cells with few lymphocytes and macrophages; thus, innate-like lymphocytes, such as γδ T cells, have important roles in antitumor immunity. Here, we characterized GBM-infiltrating γδ T cells, which may have roles in regulating the GBM tumor microenvironment and cancer cell gene expression. V(D)J repertoires of tumor-infiltrating and blood-circulating γδ T cells from four patients were analyzed by next-generation sequencing-based T-cell receptor (TCR) sequencing in addition to mutation and immune profiles in four GBM cases. In all tumor tissues, abundant innate and effector/memory lymphocytes were detected, accompanied by large numbers of tumor-associated macrophages and closely located tumor-infiltrating γδ T cells, which appear to have anti-tumor activity. The immune-related gene expression analysis using the TCGA database showed that the signature gene expression extent of γδ T cells were more associated with those of cytotoxic T and Th1 cells and M1 macrophages than those of Th2 cells and M2 macrophages. Although the most abundant γδ T cells were Vγ9Vδ2 T cells in both tumor tissues and blood, the repertoire of intratumoral Vγ9Vδ2 T cells was distinct from that of peripheral blood Vγ9Vδ2 T cells and was dominated by Vγ9Jγ2 sequences, not by canonical Vγ9JγP sequences that are mostly commonly found in blood γδ T cells. Collectively, unique GBM-specific TCR clonotypes were identified by comparing TCR repertoires of peripheral blood and intra-tumoral γδ T cells. These findings will be helpful for the elucidation of tumor-specific antigens and development of anticancer immunotherapies using tumor-infiltrating γδ T cells.

HPV, tumour metabolism and novel target identification in head and neck squamous cell carcinoma.

BACKGROUND: Metabolic changes in tumour cells are used in clinical imaging and may provide potential therapeutic targets. Human papillomavirus (HPV) status is important in classifying head and neck cancers (HNSCC), identifying a distinct clinical phenotype; metabolic differences between these HNSCC subtypes remain poorly understood. METHODS: We used RNA sequencing to classify the metabolic expression profiles of HPV+ve and HPV-ve HNSCC, performed a meta-analysis on FDG-PET imaging characteristics and correlated results with in vitro extracellular flux analysis of HPV-ve and HPV+ve HNSCC cell lines. The monocarboxylic acid transporter-1 (MCT1) was identified as a potential metabolic target and tested in functional assays. RESULTS: Specific metabolic profiles were associated with HPV status, not limited to carbohydrate metabolism. There was dominance of all energy pathways in HPV-negative disease, with elevated expression of genes associated with glycolysis and oxidative phosphorylation. In vitro analysis confirmed comparative increased rates of oxidative phosphorylation and glycolysis in HPV-negative cell lines. PET SUV(max) scores however were unable to reliably differentiate between HPV-positive and HPV-negative tumours. MCT1 expression was significantly increased in HPV-negative tumours, and inhibition suppressed tumour cell invasion, colony formation and promoted radiosensitivity. CONCLUSION: HPV-positive and negative HNSCC have different metabolic profiles which may have potential therapeutic applications.

Head and Neck Squamous Cell Carcinomas Are Characterized by a Stable Immune Signature Within the Primary Tumor Over Time and Space.

Purpose: Genetic and morphologic heterogeneity is well-documented in solid cancers. Immune cells are also variably distributed within the tumor; this heterogeneity is difficult to assess in small biopsies, and may confound our understanding of the determinants of successful immunotherapy. We examined the transcriptomic variability of the immunologic signature in head and neck squamous cell carcinoma (HNSCC) within individual tumors using transcriptomic and IHC assessments.Experimental Design: Forty-four tumor biopsies from 16 HNSCC patients, taken at diagnosis and later at resection, were analyzed using RNA-sequencing. Variance filtering was used to identify the top 4,000 most variable genes. Principal component analysis, hierarchical clustering, and correlation analysis were performed. Gene expression of CD8A was correlated to IHC analysis.Results: Analysis of immunologic gene expression was highly consistent in replicates from the same cancer. Across the cohort, samples from the same patient were most similar to each other, both spatially (at diagnosis) and, notably, over time (diagnostic biopsy compared with resection); comparison of global gene expression by hierarchical clustering (P ≤ 0.0001) and correlation analysis [median intrapatient r = 0.82; median interpatient r = 0.63]. CD8A gene transcript counts were highly correlated with CD8 T-cell counts by IHC (r = 0.82).Conclusions: Our data demonstrate that in HNSCC the global tumor and adaptive immune signatures are stable between discrete parts of the same tumor and also at different timepoints. This suggests that immunologic heterogeneity may not be a key reason for failure of immunotherapy and underpins the use of transcriptomics for immunologic evaluation of novel agents in HNSCC patients. Clin Cancer Res; 23(24); 7641-9. ©2017 AACR.

Induction of fibroblast senescence generates a non-fibrogenic myofibroblast phenotype that differentially impacts on cancer prognosis.

Cancer-associated fibroblasts (CAF) remain a poorly characterized, heterogeneous cell population. Here we characterized two previously described tumor-promoting CAF sub-types, smooth muscle actin (SMA)-positive myofibroblasts and senescent fibroblasts, identifying a novel link between the two. Analysis of CAF cultured ex vivo, showed that senescent CAF are predominantly SMA-positive; this was confirmed by immunochemistry in head & neck (HNSCC) and esophageal (EAC) cancers. In vitro, we found that fibroblasts induced to senesce develop molecular, ultrastructural and contractile features typical of myofibroblasts and this is dependent on canonical TGF-ß signaling. Similar to TGF-ß1-generated myofibroblasts, these cells secrete soluble factors that promote tumor cell motility. However, RNA-sequencing revealed significant transcriptomic differences between the two SMA-positive CAF groups, particularly in genes associated with extracellular matrix (ECM) deposition and organization, which differentially promote tumor cell invasion. Notably, second harmonic generation imaging and bioinformatic analysis of SMA-positive human HNSCC and EAC showed that collagen fiber organization correlates with poor prognosis, indicating that heterogeneity within the SMA-positive CAF population differentially impacts on survival. These results show that non-fibrogenic, SMA-positive myofibroblasts can be directly generated through induction of fibroblast senescence and suggest that senescence and myofibroblast differentiation are closely linked processes.

Gene expression analysis of TIL rich HPV-driven head and neck tumors reveals a distinct B-cell signature when compared to HPV independent tumors.

Human papilloma virus (HPV)-associated head and neck squamous cell carcinoma (HNSCC) has a better prognosis than it's HPV negative (HPV(-)) counterpart. This may be due to the higher numbers of tumor-infiltrating lymphocytes (TILs) in HPV positive (HPV(+)) tumors. RNA-Sequencing (RNA-Seq) was used to evaluate whether the differences in clinical behaviour simply reflect a numerical difference in TILs or whether there is a fundamental behavioural difference between TILs in these two settings. Thirty-nine HNSCC tumors were scored for TIL density by immunohistochemistry. After the removal of 16 TILlow tumors, RNA-Seq analysis was performed on 23 TILhigh/med tumors (HPV(+) n=10 and HPV(-) n=13). Using EdgeR, differentially expressed genes (DEG) were identified. Immune subset analysis was performed using Functional Analysis of Individual RNA-Seq/ Microarray Expression (FAIME) and immune gene RNA transcript count analysis. In total, 1,634 DEGs were identified, with a dominant immune signature observed in HPV(+) tumors. After normalizing the expression profiles to account for differences in B- and T-cell number, 437 significantly DEGs remained. A B-cell associated signature distinguished HPV(+) from HPV(-) tumors, and included the DEGs CD200, GGA2, ADAM28, STAG3, SPIB, VCAM1, BCL2 and ICOSLG; the immune signal relative to T-cells was qualitatively similar between TILs of both tumor cohorts. Our findings were validated and confirmed in two independent cohorts using TCGA data and tumor-infiltrating B-cells from additional HPV(+) HNSCC patients. A B-cell associated signal segregated tumors relative to HPV status. Our data suggests that the role of B-cells in the adaptive immune response to HPV(+) HNSCC requires re-assessment.

Upregulated Glucose Metabolism Correlates Inversely with CD8+ T-cell Infiltration and Survival in Squamous Cell Carcinoma.

Antibodies that block T-cell-regulatory checkpoints have recently emerged as a transformative approach to cancer treatment. However, the clinical efficacy of checkpoint blockade depends upon inherent tumor immunogenicity, with variation in infiltrating T cells contributing to differences in objective response rates. Here, we sought to understand the molecular correlates of tumor-infiltrating T lymphocytes (TIL) in squamous cell carcinoma (SCC), using a systems biologic approach to integrate publicly available omics datasets with histopathologic features. We provide evidence that links TIL abundance and therapeutic outcome to the regulation of tumor glycolysis by EGFR and HIF, both of which are attractive molecular targets for use in combination with immunotherapeutics. Cancer Res; 76(14); 4136-48. ©2016 AACR.