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.

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.

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.

Outlier detection at the transcriptome-proteome interface.

BACKGROUND: In high-throughput experimental biology, it is widely acknowledged that while expression levels measured at the levels of transcriptome and the corresponding proteome do not, in general, correlate well, messenger RNA levels are used as convenient proxies for protein levels. Our interest is in developing data-driven computational models that can bridge the gap between these two levels of measurement at which different mechanisms of regulation may act on different molecular species causing any observed lack of correlations. To this end, we build data-driven predictors of protein levels using mRNA levels and known proxies of translation efficiencies as covariates. Previous work showed that in such a setting, outliers with respect to the model are reliable candidates for post-translational regulation. RESULTS: Here, we present and compare two novel formulations of deriving a protein concentration predictor from which outliers may be extracted in a systematic manner. The first approach, outlier rejecting regression, allows explicit specification of a certain fraction of the data as outliers. In a regression setting, this is a non-convex optimization problem which we solve by deriving a difference of convex functions algorithm (DCA). With post-translationally regulated proteins, one expects their concentrations to be affected primarily by disruption of protein stability. Our second algorithm exploits this observation by minimizing an asymmetric loss using quantile regression and extracts outlier proteins whose measured concentrations are lower than what a genome-wide regression would predict. We validate the two approaches on a dataset of yeast transcriptome and proteome. Functional annotation check on detected outliers demonstrate that the methods are able to identify post-translationally regulated genes with high statistical confidence.

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.

Group A Streptococcus induces CD1a-autoreactive T cells and promotes psoriatic inflammation.

Group A Streptococcus (GAS) infection is associated with multiple clinical sequelae, including different subtypes of psoriasis. Such post-streptococcal disorders have been long known but are largely unexplained. CD1a is expressed at constitutively high levels by Langerhans cells and presents lipid antigens to T cells, but the potential relevance to GAS infection has not been studied. Here, we investigated whether GAS-responsive CD1a-restricted T cells contribute to the pathogenesis of psoriasis. Healthy individuals had high frequencies of circulating and cutaneous GAS-responsive CD4+ and CD8+ T cells with rapid effector functions, including the production of interleukin-22 (IL-22). Human skin and blood single-cell CITE-seq analyses of IL-22-producing T cells showed a type 17 signature with proliferative potential, whereas IFN-γ-producing T cells displayed cytotoxic T lymphocyte characteristics. Furthermore, individuals with psoriasis had significantly higher frequencies of circulating GAS-reactive T cells, enriched for markers of activation, cytolytic potential, and tissue association. In addition to responding to GAS, subsets of expanded GAS-reactive T cell clones/lines were found to be autoreactive, which included the recognition of the self-lipid antigen lysophosphatidylcholine. CD8+ T cell clones/lines produced cytolytic mediators and lysed infected CD1a-expressing cells. Furthermore, we established cutaneous models of GAS infection in a humanized CD1a transgenic mouse model and identified enhanced and prolonged local and systemic inflammation, with resolution through a psoriasis-like phenotype. Together, these findings link GAS infection to the CD1a pathway and show that GAS infection promotes the proliferation and activation of CD1a-autoreactive T cells, with relevance to post-streptococcal disease, including the pathogenesis and treatment of psoriasis.

Single cell spatial analysis reveals inflammatory foci of immature neutrophil and CD8 T cells in COVID-19 lungs.

Single cell spatial interrogation of the immune-structural interactions in COVID -19 lungs is challenging, mainly because of the marked cellular infiltrate and architecturally distorted microstructure. To address this, we develop a suite of mathematical tools to search for statistically significant co-locations amongst immune and structural cells identified using 37-plex imaging mass cytometry. This unbiased method reveals a cellular map interleaved with an inflammatory network of immature neutrophils, cytotoxic CD8 T cells, megakaryocytes and monocytes co-located with regenerating alveolar progenitors and endothelium. Of note, a highly active cluster of immature neutrophils and CD8 T cells, is found spatially linked with alveolar progenitor cells, and temporally with the diffuse alveolar damage stage. These findings offer further insights into how immune cells interact in the lungs of severe COVID-19 disease. We provide our pipeline [Spatial Omics Oxford Pipeline (SpOOx)] and visual-analytical tool, Multi-Dimensional Viewer (MDV) software, as a resource for spatial analysis.

Potential Biomarker Identification by RNA-Seq Analysis in Antibiotic-Related Drug Reaction with Eosinophilia and Systemic Symptoms (DRESS): A Pilot Study.

One of the most severe forms of cutaneous adverse drug reactions is "drug reaction with eosinophilia and systemic symptoms" (DRESS), hence subsequent avoidance of the causal drug is imperative. However, attribution of drug culpability in DRESS is challenging and standard skin allergy tests are not recommended due to patient safety reasons. Whilst incidence of DRESS is relatively low, between 1:1000 and 1:10 000 drug exposures, antibiotics are a commoner cause of DRESS and absence of confirmatory diagnostic test can result in unnecessary avoidance of efficacious treatment. We therefore sought to identify potential biomarkers for development of a diagnostic test in antibiotic-associated DRESS. Peripheral blood mononuclear cells from a "discovery" cohort (n = 5) challenged to causative antibiotic or control were analyzed for transcriptomic profile. A panel of genes was then tested in a validation cohort (n = 6) and compared with tolerant controls and other inflammatory conditions which can clinically mimic DRESS. A scoring system to identify presence of drug hypersensitivity was developed based on gene expression alterations of this panel. The DRESS transcriptomic panel identified antibiotic-DRESS cases in a validation cohort but was not altered in other inflammatory conditions. Machine learning or differential expression selection of a biomarker panel consisting of 6 genes (STAC, GPR183, CD40, CISH, CD4, and CCL8) showed high sensitivity and specificity (100% and 85.7%-100%, respectively) for identification of the culprit drug in these cohorts of antibiotic-associated DRESS. Further work is required to determine whether the same panel can be repeated for larger cohorts, different medications, and other T-cell-mediated drug hypersensitivity reactions.

Dengue virus co-opts innate type 2 pathways to escape early control of viral replication.

Mast cell products and high levels of type 2 cytokines are associated with severe dengue disease. Group 2 innate lymphoid cells (ILC2) are type-2 cytokine-producing cells that are activated by epithelial cytokines and mast cell-derived lipid mediators. Through ex vivo RNAseq analysis, we observed that ILC2 are activated during acute dengue viral infection, and show an impaired type I-IFN signature in severe disease. We observed that circulating ILC2 are permissive for dengue virus infection in vivo and in vitro, particularly when activated through prostaglandin D2 (PGD2). ILC2 underwent productive dengue virus infection, which was inhibited through CRTH2 antagonism. Furthermore, exogenous IFN-ß induced expression of type I-IFN responsive anti-viral genes by ILC2. PGD2 downregulated type I-IFN responsive gene and protein expression; and urinary prostaglandin D2 metabolite levels were elevated in severe dengue. Moreover, supernatants from activated ILC2 enhanced monocyte infection in a GM-CSF and mannan-dependent manner. Our results indicate that dengue virus co-opts an innate type 2 environment to escape early type I-IFN control and facilitate viral dissemination. PGD2 downregulates type I-IFN induced anti-viral responses in ILC2. CRTH2 antagonism may be a therapeutic strategy for dengue-associated disease.

Constraints on HIV-1 diversity from protein structure.

The high rate of HIV-1 evolution contributes to immune escape, enables the virus to escape drug therapy, and may underlie the difficulty of producing an effective vaccine. Identifying constraints on HIV evolution is therefore of prime importance. To investigate this problem, we examined the relationships between sequence diversity, selection, and protein structure. We found that while there was an increase in sequence diversity over time, this variation had a tendency to be limited to specific structural regions. When individual sites were analyzed, there was, in contrast, substantial and widespread evolutionary constraint over gag and env. This constraint was present even in the highly variable envelope proteins. The evolutionary significance of an individual site is indicated by the change in selection pressure along the time course: increasing entropy indicates that the site is evolving predominantly in a more "clock"-like manner, low entropy values with no increase indicate a high degree of constraint, and high entropy values indicate a lack of constraint. Few sites display high degrees of turnover. Mapping these sites onto the three-dimensional protein structure, we found a significant difference between evolutionary rates for regions buried in the core of the protein and those on the surface. This constraint did not change over the time period analyzed and was not subtype dependent, as similar results were found for subtypes B and C. This link between sequence and structure not only demonstrates the limits of recent HIV-1 evolution but also highlights the origins of evolutionary constraint on viral change.

Micro RNA Targets in HIV Latency: Insights into Novel Layers of Latency Control.

Despite the considerable progress that has been made in identifying cellular factors and pathways that contribute to establishment and maintenance of the latent HIV reservoir, it remains the major obstacle to eradicating this virus. Most recently, noncoding genes have been implicated in regulation of HIV expression. In this study, small RNA sequencing was used to profile expression of microRNAs (miRNAs) in a primary CD4+ T cell in vitro model of HIV latency. Previously, we have shown that protein-coding genes dysregulated in this model were enriched for the p53 signaling pathway, which was confirmed experimentally. We further found a link between p53 signaling and dysregulated long noncoding RNAs. In this study, we hypothesized that miRNAs may provide an additional level of regulation of the p53 signaling pathway during HIV latency. Twenty-six miRNAs were identified to be dysregulated in our latency model. A subset of these miRNAs was validated by real-time quantitative polymerase chain reaction. Predicted messenger RNA (mRNA) targets and cellular pathways enriched for mRNA targets were identified using several analytical methods. Our analyses showed that many protein-coding genes and pathways targeted by dysregulated miRNAs have relevance to regulation of HIV expression or establishment of HIV latency. The p53 signaling pathway was found among pathways that were targeted by dysregulated miRNAs at a greater level than expected by chance. This study provides a mechanistic insight into regulation of the p53 pathway through miRNAs that may contribute to the establishment of latency.

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.

Upregulation of epithelial metallothioneins by metal-rich ultrafine particulate matter from an underground railway.

Airborne particulate matter (PM) is a leading cause of mortality and morbidity. However, understanding of the range and mechanisms of effects of PM components is poor. PM generated in underground railways is rich in metals, especially iron. In the ultrafine (UFPM; <0.1 µm diameter) fraction, the combination of small size and metal enrichment poses an unknown health risk. This study aimed to analyse transcriptomic responses to underground UFPM in primary bronchial epithelial cells (PBECs), a key site of PM deposition. The oxidation state of iron in UFPM from an underground station was determined by X-ray absorption near edge structure (XANES) spectroscopy. Antioxidant response was assayed using a reporter cell line transfected with an antioxidant response element (ARE)-luciferase construct. Differentiated PBECs were exposed to UFPM for 6 h or 24 h for RNA-Seq and RT-qPCR analysis. XANES showed predominance of redox-active Fe3O4, with ROS generation confirmed by induction of ARE-luciferase expression. 6 h exposure of PBECs to UFPM identified 52 differentially expressed genes (DEGs), especially associated with epithelial maintenance, whereas 24 h exposure yielded 23 DEGs, particularly involved with redox homeostasis and metal binding. At both timepoints, there was upregulation of members of the metallothionein family, low molecular weight proteins with antioxidant activity whose main function is binding and homeostasis of zinc and copper ions, but not iron ions. This upregulation was partially inhibited by metal chelation or ROS scavenging. These data suggest differential regulation of responses to metal-rich UFPM depending on exposure period, and highlight novel pathways and markers of PM exposure, with the role of metallothioneins warranting further investigation.

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.

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.

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.

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.

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.

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.