Immune Landscape of Viral- and Carcinogen-Driven Head and Neck Cancer

The diversification of methods for studying cell-cell interactions and communication

No cell lives in a vacuum, and the molecular interactions between cells define most phenotypes. Transcriptomics provides rich information to infer cell-cell interactions and communication, thus accelerating the discovery of the roles of cells within their communities. Such research relies heavily on algorithms that infer which cells are interacting and the ligands and receptors involved. Specific pressures on different research niches are driving the evolution of next-generation computational tools, enabling new conceptual opportunities and technological advances. More sophisticated algorithms now account for the heterogeneity and spatial organization of cells, multiple ligand types and intracellular signalling events, and enable the use of larger and more complex datasets, including single-cell and spatial transcriptomics. Similarly, new high-throughput experimental methods are increasing the number and resolution of interactions that can be analysed simultaneously. Here, we explore recent progress in cell-cell interaction research and highlight the diversification of the next generation of tools, which have yielded a rich ecosystem of tools for different applications and are enabling invaluable discoveries.

RNA-binding protein HuR reprograms immune T cells and promotes oral squamous cell carcinoma

Hu Antigen R, also known as ELAVL1 (HuR), is a key posttranscriptional regulator in eukaryotic cells. HuR overexpression promotes several malignancies, including head and neck squamous cell carcinoma (HNSCC). However, its immune dysfunction-associated tumorigenesis pathways remain unknown. We examined HuR's effects on oral malignancies and immune cell function in vitro and in vivo using oral carcinoma cells and transgenic HuR knockout (KO) mice. CRISPR/Cas9-mediated HuR deletion in mice syngeneic oral cancer cells eliminated colony formation and tumor development. HuR-KO tumors had a lower tumor volume, fewer CD4+CD25+FoxP3+ regulatory T cells, and more CD8+ T cells, suggesting that HuR may suppress the immune response during oral cancer progression. In contrast, HuR KO oral epithelial tissues are resistant to 4NQO-induced oral malignancies compared to control tumor-bearing mice. HuR KO mice showed fewer Tregs and greater IFN levels than WT tumor-bearing mice, suggesting anticancer activity. Finally, the HuR inhibitor pyrvinium pamoate lowers tumor burden by enhancing CD8+ infiltration at the expense of CD4+, suggesting anticancer benefits. Thus, HuR-dependent oral neoplasia relies on immunological dysfunction, suggesting that decreasing HuR may boost antitumor potential and offer a novel HNSCC therapy.

Metabolic pathway-based subtypes associate glycan biosynthesis and treatment response in head and neck cancer

Head and Neck Squamous Cell Carcinoma (HNSCC) is a heterogeneous malignancy that remains a significant challenge in clinical management due to frequent treatment failures and pronounced therapy resistance. While metabolic dysregulation appears to be a critical factor in this scenario, comprehensive analyses of the metabolic HNSCC landscape and its impact on clinical outcomes are lacking. This study utilized transcriptomic data from four independent clinical cohorts to investigate metabolic heterogeneity in HNSCC and define metabolic pathway-based subtypes (MPS). In HPV-negative HNSCCs, MPS1 and MPS2 were identified, while MPS3 was enriched in HPV-positive cases. MPS classification was associated with clinical outcome post adjuvant radio(chemo)therapy, with MPS1 consistently exhibiting the highest risk of therapeutic failure. MPS1 was uniquely characterized by upregulation of glycan (particularly chondroitin/dermatan sulfate) metabolism genes. Immunohistochemistry and pilot mass spectrometry imaging analyses confirmed this at metabolite level. The histological context and single-cell RNA sequencing data identified the malignant cells as key contributors. Globally, MPS1 was distinguished by a unique transcriptomic landscape associated with increased disease aggressiveness, featuring motifs related to epithelial-mesenchymal transition, immune signaling, cancer stemness, tumor microenvironment assembly, and oncogenic signaling. This translated into a distinct histological appearance marked by extensive extracellular matrix remodeling, abundant spindle-shaped cancer-associated fibroblasts, and intimately intertwined populations of malignant and stromal cells. Proof-of-concept data from orthotopic xenotransplants replicated the MPS phenotypes on the histological and transcriptome levels. In summary, this study introduces a metabolic pathway-based classification of HNSCC, pinpointing glycan metabolism-enriched MPS1 as the most challenging subgroup that necessitates alternative therapeutic strategies.

The diversity of inhibitory receptor co-expression patterns of exhausted CD8+ T cells in oropharyngeal carcinoma

Exhausted CD8+ T cells (Texs) are characterized by the expression of various inhibitory receptors (IRs), whereas the functional attributes of these co-expressed IRs remain limited. Here, we systematically characterized the diversity of IR co-expression patterns in Texs from both human oropharyngeal squamous cell carcinoma (OPSCC) tissues and syngeneic OPSCC model. Nearly 60% of the Texs population co-expressed two or more IRs, and the number of co-expressed IRs was positively associated with superior exhaustion and cytotoxicity phenotypes. In OPSCC patients, programmed cell death-1 (PD-1) blockade significantly enhanced PDCD1-based co-expression with other IR genes, whereas dual blockades of PD-1 and cytotoxic T lymphocyte-associated protein 4 (CTLA-4) significantly upregulated CTLA4-based co-expression with other IR genes. Collectively, our findings demonstrate that highly diverse IR co-expression is a leading feature of Texs and represents their functional states, which might provide essential clues for the rational selection of immune checkpoint inhibitors in treating OPSCC.

Distinct molecular profiles drive multifaceted characteristics of colorectal cancer metastatic seeds

Metastasis of primary tumors remains a challenge for early diagnosis and prevention. The cellular properties and molecular drivers of metastatically competent clones within primary tumors remain unclear. Here, we generated 10-16 single cell-derived lines from each of three colorectal cancer (CRC) tumors to identify and characterize metastatic seeds. We found that intrinsic factors conferred clones with distinct metastatic potential and cellular communication capabilities, determining organ-specific metastasis. Poorly differentiated or highly metastatic clones, rather than drug-resistant clones, exhibited poor clinical prognostic impact. Personalized genetic alterations, instead of mutation burden, determined the occurrence of metastatic potential during clonal evolution. Additionally, we developed a gene signature for capturing metastatic potential of primary CRC tumors and demonstrated a strategy for identifying metastatic drivers using isogenic clones with distinct metastatic potential in primary tumors. This study provides insight into the origin and mechanisms of metastasis and will help develop potential anti-metastatic therapeutic targets for CRC patients.

Association of immune cells and the risk of esophageal cancer: A Mendelian randomization study in a East Asian population

Immunotherapy has been used in esophageal cancer (EC), but the causal relationship between EC and immune cells is not clear. Although the cellular phenotype has been reported as a biomarker for immunotherapy, the biomarker studies for immunotherapy in EC still face great challenges. Comprehensive 2-sample Mendelian randomization (MR) analysis was performed to determine the causal association between immune cell signatures and EC in this study. Based on publicly available genetic data, we explored causal associations between 731 immune cell signatures and EC risk. EC had no statistically significant effect on immunophenotypes. Nine immunophenotype types were positively associated with the risk of EC: CD20-%B cell, CD20% lymphocytes, CD25 on IgD- CD27-, CD25 on IgD+ CD24+, CD27 on IgD+ CD24+, CD28+ CD45RA- CD8br AC, CD3 on TD CD8br, IgD-CD38dim%B cells, and Mo MDSC AC. In addition, a total of 15 immunophenotypes were identified as causally associated with EC. IgD+ CD38- %B cell, IgD- CD24- %lymphocyte, CD19 on IgD- CD38dim, CD20 on IgD+ CD24+, CD62L-myeloid DC AC, CD4+ AC, Lymphocyte %leukocyte, CD3 on HLA-DR+ T cell, CD3 on CD45RA- CD4+, HVEM on naive CD4+ AC, HVEM on CD45RA- CD4+, CD4 on TD CD4+, CD4 on CD4 Treg, and CD4 on CD39+ resting Treg, and CD4 on activated & secreting Treg. Our study has demonstrated the close connection between immune cells and EC by genetic means, thus providing guidance for future clinical research.

Neoadjuvant immunotherapy for head and neck squamous cell carcinoma

The neoadjuvant immunotherapy approach marks a significant shift in the treatment paradigm of potentially curable HNSCC. Here, current therapies, despite being highly individualized and advanced, often fall short in achieving satisfactory long-term survival rates and are frequently associated with substantial morbidity.The primary advantage of this approach lies in its potential to intensify and enhance treatment regimens, offering a distinct modality that complements the existing triad of surgery, radiotherapy, and chemotherapy. Checkpoint inhibitors have been at the forefront of this evolution. Demonstrating moderate yet significant survival benefits in the recurrent-metastatic setting with a relatively better safety profile compared to conventional treatments, these agents hold promise when considered for earlier stages of HNSCC.On the other hand, a significant potential benefit of introducing immunotherapy in the neoadjuvant phase is the possibility of treatment de-escalation. By reducing the tumor burden before surgery, this strategy could lead to less invasive surgical interventions. The prospect of organ-sparing protocols becomes a realistic and highly valued goal in this context. Further, the early application of immunotherapy might catalyze a more effective and durable immune response. The induction of an immune memory may potentially lead to a more effective surveillance of residual disease, decreasing the rates of local, regional, and distant recurrences, thereby enhancing overall and recurrence-free survival.However, neoadjuvant immunotherapy is not without its challenges. One of the primary concerns is the safety and adverse events profile. While data suggest that adverse events are relatively rare and manageable, the long-term safety profile in the neoadjuvant setting, especially in the context of curative intent, remains a subject for ongoing research. Another unsolved issue lies in the accurate assessment of treatment response. The discrepancy between radiographic assessment using RECIST criteria and histological findings has been noted, indicating a gap in current imaging techniques' ability to accurately reflect the true efficacy of immunotherapy. This gap underscores the necessity for improved imaging methodologies and the development of new radiologic and pathologic criteria tailored to evaluate the response to immunotherapy accurately.Treatment combinations and timing represent another layer of complexity. There is a vast array of possibilities in combining immunotherapy agents with conventional chemotherapy, targeted therapy, radiation, and other experimental treatments. Determining the optimal treatment regimen for individual patients becomes an intricate task, especially when comparing small, single-arm, non-randomized trials with varying regimens and outcome measures.Moreover, one needs to consider the importance of pre- and intraoperative decision-making in the context of neoadjuvant immunotherapy. As experience with this treatment paradigm grows, there is potential for more tailored surgical approaches based on the patient's remaining disease post-neoadjuvant treatment. This consideration is particularly relevant in extensive surgeries, where organ-sparing protocols could be evaluated.In practical terms, the multi-modal nature of this treatment strategy introduces complexities, especially outside clinical trial settings. Patients face challenges in navigating the treatment landscape, which involves coordination across multiple medical disciplines, highlighting the necessity for streamlined care pathways at specialized centers to facilitate effective treatment management if the neoadjuvant approach is introduced to the real-world.These potential harms and open questions underscore the critical need for meticulously designed clinical trials and correlational studies to ensure patient safety and efficacy. Only these can ensure that this new treatment approach is introduced in a safe way and fulfils the promise it theoretically holds.

Evolving Tumor Characteristics and Smart Nanodrugs for Tumor Immunotherapy

Typical physiological characteristics of tumors, such as weak acidity, low oxygen content, and upregulation of certain enzymes in the tumor microenvironment (TME), provide survival advantages when exposed to targeted attacks by drugs and responsive nanomedicines. Consequently, cancer treatment has significantly progressed in recent years. However, the evolution and adaptation of tumor characteristics still pose many challenges for current treatment methods. Therefore, efficient and precise cancer treatments require an understanding of the heterogeneity degree of various factors in cancer cells during tumor evolution to exploit the typical TME characteristics and manage the mutation process. The highly heterogeneous tumor and infiltrating stromal cells, immune cells, and extracellular components collectively form a unique TME, which plays a crucial role in tumor malignancy, including proliferation, invasion, metastasis, and immune escape. Therefore, the development of new treatment methods that can adapt to the evolutionary characteristics of tumors has become an intense focus in current cancer treatment research. This paper explores the latest understanding of cancer evolution, focusing on how tumors use new antigens to shape their "new faces"; how immune system cells, such as cytotoxic T cells, regulatory T cells, macrophages, and natural killer cells, help tumors become "invisible", that is, immune escape; whether the diverse cancer-associated fibroblasts provide support and coordination for tumors; and whether it is possible to attack tumors in reverse. This paper discusses the limitations of targeted therapy driven by tumor evolution factors and explores future strategies and the potential of intelligent nanomedicines, including the systematic coordination of tumor evolution factors and adaptive methods, to meet this therapeutic challenge.

Single Cell Atlas: a single-cell multi-omics human cell encyclopedia

Single-cell sequencing datasets are key in biology and medicine for unraveling insights into heterogeneous cell populations with unprecedented resolution. Here, we construct a single-cell multi-omics map of human tissues through in-depth characterizations of datasets from five single-cell omics, spatial transcriptomics, and two bulk omics across 125 healthy adult and fetal tissues. We construct its complement web-based platform, the Single Cell Atlas (SCA, www.singlecellatlas.org ), to enable vast interactive data exploration of deep multi-omics signatures across human fetal and adult tissues. The atlas resources and database queries aspire to serve as a one-stop, comprehensive, and time-effective resource for various omics studies.

Single-Cell Informatics for Tumor Microenvironment and Immunotherapy

Cancer comprises malignant cells surrounded by the tumor microenvironment (TME), a dynamic ecosystem composed of heterogeneous cell populations that exert unique influences on tumor development. The immune community within the TME plays a substantial role in tumorigenesis and tumor evolution. The innate and adaptive immune cells "talk" to the tumor through ligand-receptor interactions and signaling molecules, forming a complex communication network to influence the cellular and molecular basis of cancer. Such intricate intratumoral immune composition and interactions foster the application of immunotherapies, which empower the immune system against cancer to elicit durable long-term responses in cancer patients. Single-cell technologies have allowed for the dissection and characterization of the TME to an unprecedented level, while recent advancements in bioinformatics tools have expanded the horizon and depth of high-dimensional single-cell data analysis. This review will unravel the intertwined networks between malignancy and immunity, explore the utilization of computational tools for a deeper understanding of tumor-immune communications, and discuss the application of these approaches to aid in diagnosis or treatment decision making in the clinical setting, as well as the current challenges faced by the researchers with their potential future improvements.

Heterogeneous characterization of neutrophilic cells in head and neck cancers

BACKGROUND

Neutrophilic cells are among the most abundant immune populations within the head and neck tumor microenvironment (TME) and harbor multiple mechanisms of immunosuppression. Despite these important features, neutrophilic cells may be underrepresented in contemporary studies that aim to comprehensively characterize the immune landscape of the TME due to discrepancies in tissue processing and analysis techniques. Here, we review the role of pathologically activated neutrophilic cells within the TME and pitfalls of various approaches used to study their frequency and function in clinical samples.

METHODS

The literature was identified by searching PubMed for "immune landscape" and "tumor immune microenvironment" in combination with keywords describing solid tumor malignancies. Key publications that assessed the immune composition of solid tumors derived from human specimens were included. The tumor and blood processing methodologies in each study were reviewed in depth and correlated with the reported abundance of neutrophilic cells.

RESULTS

Neutrophilic cells do not survive cryopreservation, and many studies fail to identify and study neutrophilic cell populations due to cryopreservation of clinical samples for practical reasons. Additional single-cell transcriptomic studies filter out neutrophilic cells due to low transcriptional counts.

CONCLUSIONS

This report can help readers critically interpret studies aiming to comprehensively study the immune TME that fail to identify and characterize neutrophilic cells.

Exploring dendritic cell subtypes in cancer immunotherapy: unraveling the role of mature regulatory dendritic cells

Dendritic cells (DCs) possess a specialized function in presenting antigens and play pivotal roles in both innate and adaptive immune responses. Their ability to cross-present antigens from tumor cells to naïve T cells is instrumental in generating specific T-cell-mediated antitumor responses, crucial for controlling tumor growth and preventing tumor cell dissemination. However, within a tumor immune microenvironment (TIME), the functions of DCs can be significantly compromised. This review focuses on the profile, function, and activation of DCs, leveraging recent studies that reveal insights into their phenotype acquisition, transcriptional state, and functional programs through single-cell RNA sequence (scRNA-seq) analysis. Additionally, the therapeutic potential of DC-mediated tumor antigen sensing in priming antitumor immunity is discussed.

CHD4 acts as a prognostic factor and drives radioresistance in HPV negative HNSCC

Despite great efforts in improving existing therapies, the outcome of patients with advanced radioresistant HPV-negative head and neck squamous cell carcinoma (HNSCC) remains poor. The chromatin remodeler Chromodomain helicase DNA binding protein 4 (CHD4) is involved in different DNA-repair mechanisms, but the role and potential in HNSCC has not been explored yet. In the present study, we evaluated the prognostic significance of CHD4 expression using in silico analysis of the pan-cancer dataset. Furthermore, we established a monoclonal HNSCC CHD4 knockdown cell clone utilizing the CRISPR/Cas9 system. Effects of lower CHD4 expression on radiosensitivity after increasing doses of ionizing radiation were characterized using clonogenic assays and cell numbers. The in silico analysis revealed that high CHD4 expression is associated with significant poorer overall survival of HPV-negative HNSCC patients. Additionally, the knockdown of CHD4 significantly increased the radiosensitivity of HNSCC cells. Therefore, CHD4 might be involved in promoting radioresistance in hard-to-treat HPV-negative HNSCC entities. We conclude that CHD4 could serve as a prognostic factor in HPV-negative HNSCC tumors and is a potential target protein overcoming radioresistance in HNSCC. Our results and the newly established cell clone laid the foundation to further characterize the underlying mechanisms and ultimately use CHD4 in HNSCC therapies.

Tumor-specific T cells in head and neck cancer have rescuable functionality and can be identified through single-cell co-culture

BACKGROUND

Human papillomavirus (HPV)-negative head and neck squamous cell carcinoma (HNSCC) remains a treatment-resistance disease with limited response to immunotherapy. While T cells in HNSCC are known to display phenotypic dysfunction, whether they retain rescuable functional capacity and tumor-killing capability remains unclear.

METHODS

To investigate the functionality and tumor-specificity of tumor-infiltrating lymphocytes (TILs) across HNSCCs, malignant cell lines and TILs were derived from 31 HPV-negative HNSCCs at the time of standard surgical resection. T cell functional capacity was evaluated through ex vivo expansion, immunophenotyping, and IsoLight single-cell proteomics. Tumor-specificity was investigated through both bulk and single-cell tumor-TIL co-culture.

RESULTS

TILs could be successfully generated from 24 patients (77%), including both previously untreated and radiation recurrent HNSCCs. We demonstrate that across HNSCCs, TILs express multiple exhaustion markers but maintain a predominantly effector memory phenotype. After ex vivo expansion, TILs retain immunogenic functionality even from radiation-resistant, exhausted, and T cell-depleted disease. We further demonstrate tumor-specificity of T cells across HNSCC patients through patient-matched malignant cell-T cell co-culture. Finally, we use optofluidic technology to establish an autologous single tumor cell-single T cell co-culture platform for HNSCC. Cells derived from three HNSCC patients underwent single-cell co-culture which enabled identification and visualization of individual tumor-killing TILs in real-time in all patients.

CONCLUSIONS

These studies show that cancer-specific T cells exist across HNSCC patients with rescuable immunogenicity and can be identified on a single-cell level. These data lay the foundation for development of patient-specific T cell immunotherapies in HNSCC.

Dendritic cells as orchestrators of anticancer immunity and immunotherapy

Dendritic cells (DCs) are a heterogeneous group of antigen-presenting innate immune cells that regulate adaptive immunity, including against cancer. Therefore, understanding the precise activities of DCs in tumours and patients with cancer is important. The classification of DC subsets has historically been based on ontogeny; however, single-cell analyses are now additionally revealing a diversity of functional states of DCs in cancer. DCs can promote the activation of potent antitumour T cells and immune responses via numerous mechanisms, although they can also be hijacked by tumour-mediated factors to contribute to immune tolerance and cancer progression. Consequently, DC activities are often key determinants of the efficacy of immunotherapies, including immune-checkpoint inhibitors. Potentiating the antitumour functions of DCs or using them as tools to orchestrate short-term and long-term anticancer immunity has immense but as-yet underexploited therapeutic potential. In this Review, we outline the nature and emerging complexity of DC states as well as their functions in regulating adaptive immunity across different cancer types. We also describe how DCs are required for the success of current immunotherapies and explore the inherent potential of targeting DCs for cancer therapy. We focus on novel insights on DCs derived from patients with different cancers, single-cell studies of DCs and their relevance to therapeutic strategies.

IL7 in combination with radiotherapy stimulates a memory T-cell response to improve outcomes in HNSCC models

Clinically approved head and neck squamous cell carcinoma (HNSCC) immunotherapies manipulate the immune checkpoint blockade (ICB) axis but have had limited success outside of recurrent/metastatic disease. Interleukin-7 (IL7) has been shown to be essential for effector T-cell survival, activation, and proliferation. Here, we show that IL7 in combination with radiotherapy (RT) is effective in activating CD8 + T-cells for reducing tumor growth. Our studies were conducted using both human papillomavirus related and unrelated orthotopic HNSCC murine models. Immune populations from the tumor, draining lymph nodes, and blood were compared between treatment groups and controls using flow cytometry, proteomics, immunofluorescence staining, and RNA sequencing. Treatment with RT and IL7 (RT + IL7) resulted in significant tumor growth reduction, high CD8 T-cell tumor infiltration, and increased proliferation of T-cell progenitors in the bone marrow. IL7 also expanded a memory-like subpopulation of CD8 T-cells. These results indicate that IL7 in combination with RT can serve as an effective immunotherapy strategy outside of the conventional ICB axis to drive the antitumor activity of CD8 T-cells.

Modulator of TMB-associated immune infiltration (MOTIF) predicts immunotherapy response and guides combination therapy

Patients with high tumor mutational burden (TMB) levels do not consistently respond to immune checkpoint inhibitors (ICIs), possibly because a high TMB level does not necessarily result in adequate infiltration of CD8+ T cells. Using bulk ribonucleic acid sequencing (RNA-seq) data from 9311 tumor samples across 30 cancer types, we developed a novel tool called the modulator of TMB-associated immune infiltration (MOTIF), which comprises genes that can determine the extent of CD8+ T cell infiltration prompted by a certain TMB level. We confirmed that MOTIF can accurately reflect the integrity and defects of the cancer-immunity cycle. By analyzing 84 human single-cell RNA-seq datasets from 32 types of solid tumors, we revealed that MOTIF can provide insights into the diverse roles of various cell types in the modulation of CD8+ T cell infiltration. Using pretreatment RNA-seq data from 13 ICI-treated cohorts, we validated the use of MOTIF in predicting CD8+ T cell infiltration and ICI efficacy. Among the components of MOTIF, we identified EMC3 as a negative regulator of CD8+ T cell infiltration, which was validated via in vivo studies. Additionally, MOTIF provided guidance for the potential combinations of programmed death 1 blockade with certain immunostimulatory drugs to facilitate CD8+ T cell infiltration and improve ICI efficacy.

SEGCECO: Subgraph Embedding of Gene expression matrix for prediction of CEll-cell COmmunication

Recent advances in single-cell RNA sequencing technology have eased analyses of signaling networks of cells. Recently, cell-cell interaction has been studied based on various link prediction approaches on graph-structured data. These approaches have assumptions about the likelihood of node interaction, thus showing high performance for only some specific networks. Subgraph-based methods have solved this problem and outperformed other approaches by extracting local subgraphs from a given network. In this work, we present a novel method, called Subgraph Embedding of Gene expression matrix for prediction of CEll-cell COmmunication (SEGCECO), which uses an attributed graph convolutional neural network to predict cell-cell communication from single-cell RNA-seq data. SEGCECO captures the latent and explicit attributes of undirected, attributed graphs constructed from the gene expression profile of individual cells. High-dimensional and sparse single-cell RNA-seq data make converting the data into a graphical format a daunting task. We successfully overcome this limitation by applying SoptSC, a similarity-based optimization method in which the cell-cell communication network is built using a cell-cell similarity matrix which is learned from gene expression data. We performed experiments on six datasets extracted from the human and mouse pancreas tissue. Our comparative analysis shows that SEGCECO outperforms latent feature-based approaches, and the state-of-the-art method for link prediction, WLNM, with 0.99 ROC and 99% prediction accuracy. The datasets can be found at https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE84133 and the code is publicly available at Github https://github.com/sheenahora/SEGCECO and Code Ocean https://codeocean.com/capsule/8244724/tree.

Single cell deciphering of progression trajectories of the tumor ecosystem in head and neck cancer

Head and neck squamous cell carcinoma is the sixth most common cancer worldwide and has high heterogeneity and unsatisfactory outcomes. To better characterize the tumor progression trajectory, we perform single-cell RNA sequencing of normal tissue, precancerous tissue, early-stage, advanced-stage cancer tissue, lymph node, and recurrent tumors tissue samples. We identify the transcriptional development trajectory of malignant epithelial cells and a tumorigenic epithelial subcluster regulated by TFDP1. Furthermore, we find that the infiltration of POSTN+ fibroblasts and SPP1+ macrophages gradually increases with tumor progression; their interaction or interaction with malignant cells also gradually increase to shape the desmoplastic microenvironment and reprogram malignant cells to promote tumor progression. Additionally, we demonstrate that during lymph node metastasis, exhausted CD8+ T cells with high CXCL13 expression strongly interact with tumor cells to acquire more aggressive phenotypes of extranodal expansion. Finally, we delineate the distinct features of malignant epithelial cells in primary and recurrent tumors, providing a theoretical foundation for the precise selection of targeted therapy for tumors at different stages. In summary, the current study offers a comprehensive landscape and deep insight into epithelial and microenvironmental reprogramming throughout initiation, progression, lymph node metastasis and recurrence of head and neck squamous cell carcinoma.

Reciprocal inhibition between TP63 and STAT1 regulates anti-tumor immune response through interferon-γ signaling in squamous cancer

Squamous cell carcinomas (SCCs) are common and aggressive malignancies. Immune check point blockade (ICB) therapy using PD-1/PD-L1 antibodies has been approved in several types of advanced SCCs. However, low response rate and treatment resistance are common. Improving the efficacy of ICB therapy requires better understanding of the mechanism of immune evasion. Here, we identify that the SCC-master transcription factor TP63 suppresses interferon-γ (IFNγ) signaling. TP63 inhibition leads to increased CD8+ T cell infiltration and heighten tumor killing in in vivo syngeneic mouse model and ex vivo co-culture system, respectively. Moreover, expression of TP63 is negatively correlated with CD8+ T cell infiltration and activation in patients with SCC. Silencing of TP63 enhances the anti-tumor efficacy of PD-1 blockade by promoting CD8+ T cell infiltration and functionality. Mechanistically, TP63 and STAT1 mutually suppress each other to regulate the IFNγ signaling by co-occupying and co-regulating their own promoters and enhancers. Together, our findings elucidate a tumor-extrinsic function of TP63 in promoting immune evasion of SCC cells. Over-expression of TP63 may serve as a biomarker predicting the outcome of SCC patients treated with ICB therapy, and targeting TP63/STAT/IFNγ axis may enhance the efficacy of ICB therapy for this deadly cancer.

Combined IL6 and CCR2 blockade potentiates antitumor activity of NK cells in HPV-negative head and neck cancer

BACKGROUND

While T cell-activating immunotherapies against recurrent head and neck squamous cell carcinoma (HNSCC) have shown impressive results in clinical trials, they are often ineffective in the majority of patients. NK cells are potential targets for immunotherapeutic intervention; however, the setback in monalizumab-based therapy in HNSCC highlights the need for an alternative treatment to enhance their antitumor activity.

METHODS

Single-cell RNA sequencing (scRNA-seq) and TCGA HNSCC datasets were used to identify key molecular alterations in NK cells. Representative HPV-positive ( +) and HPV-negative ( -) HNSCC cell lines and orthotopic mouse models were used to validate the bioinformatic findings. Changes in immune cells were examined by flow cytometry and immunofluorescence.

RESULTS

Through integration of scRNA-seq data with TCGA data, we found that the impact of IL6/IL6R and CCL2/CCR2 signaling pathways on evasion of immune attack by NK cells is more pronounced in the HPV - HNSCC cohort compared to the HPV + HNSCC cohort. In orthotopic mouse models, blocking IL6 with a neutralizing antibody suppressed HPV - but not HPV + tumors, which was accompanied by increased tumor infiltration and proliferation of CD161+ NK cells. Notably, combining the CCR2 chemokine receptor antagonist RS504393 with IL6 blockade resulted in a more pronounced antitumor effect that was associated with more activated intratumoral NK cells in HPV - HNSCC compared to either agent alone.

CONCLUSIONS

These findings demonstrate that dual blockade of IL6 and CCR2 pathways effectively enhances the antitumor activity of NK cells in HPV-negative HNSCC, providing a novel strategy for treating this type of cancer.

The prognostic value of LAYN in HPV-related head and neck squamous cell carcinoma and its influence on immune cell infiltration

BACKGROUND

HPV-positive head and neck squamous cell carcinoma (HNSCC) exhibits different characteristics from HPV-negative tumors in terms of tumor development, clinical features, treatment response, and prognosis. Layilin (LAYN), which contains homology with C-type lectins, plays a critical role in tumorigenesis and cancer progression. However, the prognostic value of LAYN and the relationship between LAYN and immune infiltration levels in HPV-related HNSCC patients still require a comprehensive understanding. Herein, we aimed to assess the prognostic value of LAYN and to investigate its underlying immunological function in HPV-related HNSCC.

METHODS

Through various bioinformatics methods, we analyzed the data from The Cancer Genome Atlas (TCGA), Tumor Immune Estimation Resource (TIMER) and Gene Expression Profiling Interactive Analysis (GEPIA) databases to explore the potential underlying oncogenic impression of LAYN, including the relevance of LAYN to survival outcomes, clinicopathological factors, immune cell infiltration, and immune marker sets in HPV-related HNSCC. The expression levels of LAYN and HPV were also verified in HNSCC patient tissues.

RESULTS

LAYN was differentially expressed in a variety of tumors. The expression of LAYN in HNSCC was significantly higher than that in adjacent normal tissues (P < 0.0001), and high expression of LAYN was correlated with poor overall survival (OS) in HNSCC patients (Hazard Ratio (HR) = 1.3, P = 0.035). Moreover, LAYN expression level in HPV-positive HNSCC patients was significantly lower than that in HPV-negative patients, with HPV-positive HNSCC patients displaying a trend of favorable prognosis. In addition, the relationship between LAYN expression and immune infiltration levels in HPV-positive HNSCC group was less tightly correlated than that in HPV-negative HNSCC group, and there was a strong relationship between LAYN expression and markers of M2 macrophage (P < 0.001) and exhausted T cells (P < 0.05) in HPV-negative HNSCC. Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis suggested that LAYN potentially influenced tumor progression through HPV infection and other cancer-related pathways.

CONCLUSIONS

LAYN might contribute to tumorigenesis via its positive correlation with immune checkpoint molecules and tumor-associated macrophages (TAMs). Our study might provide a novel prognostic biomarker and latent therapeutic target for the treatment of HPV-related HNSCC.

B cell responses and antibody-based therapeutic perspectives in human cancers

BACKGROUND

Immuno-oncology has been focused on T cell-centric approaches until the field recently started appreciating the importance of tumor-reactive antibody production by tumor-infiltrating plasma B cells, and the necessity of developing novel therapeutic antibodies for the treatment of different cancers.

RECENT FINDINGS

B lymphocytes often infiltrate solid tumors and the extent of B cell infiltration normally correlates with stronger T cell responses while generating humoral responses against malignant progression by producing tumor antigens-reactive antibodies that bind and coat the tumor cells and promote cytotoxic effector mechanisms, reiterating the fact that the adaptive immune system works by coordinated humoral and cellular immune responses. Isotypes, magnitude, and the effector functions of antibodies produced by the B cells within the tumor environment differ among cancer types. Interestingly, apart from binding with specific tumor antigens, antibodies produced by tumor-infiltrating B cells could bind to some non-specific receptors, peculiarly expressed by cancer cells. Antibody-based immunotherapies have revolutionized the modalities of cancer treatment across the world but are still limited against hematological malignancies and a few types of solid tumor cancers with a restricted number of targets, which necessitates the expansion of the field to have newer effective targeted antibody therapeutics.

CONCLUSION

Here, we discuss about recent understanding of the protective spontaneous antitumor humoral responses in human cancers, with an emphasis on the advancement and future perspectives of antibody-based immunotherapies in cancer.

Navigating the landscapes of spatial transcriptomics: How computational methods guide the way

Spatially resolved transcriptomics has been dramatically transforming biological and medical research in various fields. It enables transcriptome profiling at single-cell, multi-cellular, or sub-cellular resolution, while retaining the information of geometric localizations of cells in complex tissues. The coupling of cell spatial information and its molecular characteristics generates a novel multi-modal high-throughput data source, which poses new challenges for the development of analytical methods for data-mining. Spatial transcriptomic data are often highly complex, noisy, and biased, presenting a series of difficulties, many unresolved, for data analysis and generation of biological insights. In addition, to keep pace with the ever-evolving spatial transcriptomic experimental technologies, the existing analytical theories and tools need to be updated and reformed accordingly. In this review, we provide an overview and discussion of the current computational approaches for mining of spatial transcriptomics data. Future directions and perspectives of methodology design are proposed to stimulate further discussions and advances in new analytical models and algorithms. This article is categorized under: RNA Methods > RNA Analyses in Cells RNA Evolution and Genomics > Computational Analyses of RNA RNA Export and Localization > RNA Localization.

Identifying potential ligand-receptor interactions based on gradient boosted neural network and interpretable boosting machine for intercellular communication analysis

Cell-cell communication is essential to many key biological processes. Intercellular communication is generally mediated by ligand-receptor interactions (LRIs). Thus, building a comprehensive and high-quality LRI resource can significantly improve intercellular communication analysis. Meantime, due to lack of a "gold standard" dataset, it remains a challenge to evaluate LRI-mediated intercellular communication results. Here, we introduce CellGiQ, a high-confident LRI prediction framework for intercellular communication analysis. Highly confident LRIs are first inferred by LRI feature extraction with BioTriangle, LRI selection using LightGBM, and LRI classification based on ensemble of gradient boosted neural network and interpretable boosting machine. Subsequently, known and identified high-confident LRIs are filtered by combining single-cell RNA sequencing (scRNA-seq) data and further applied to intercellular communication inference through a quartile scoring strategy. To validation the predictions, CellGiQ exploited several evaluation strategies: using AUC and AUPR, it surpassed six competing LRI prediction models on four LRI datasets; through Venn diagrams and molecular docking, its predicted LRIs were validated by five other popular intercellular communication inference methods; based on the overlapping LRIs, it computed high Jaccard index with six other state-of-the-art intercellular communication prediction tools within human HNSCC tissues; by comparing with classical models and literature retrieve, its inferred HNSCC-related intercellular communication results was further validated. The novelty of this study is to identify high-confident LRIs based on machine learning as well as design several LRI validation ways, providing reference for computational LRI prediction. CellGiQ provides an open-source and useful tool to decompose LRI-mediated intercellular communication at single cell resolution. CellGiQ is freely available at https://github.com/plhhnu/CellGiQ.

Development of a novel head and neck squamous cell carcinoma prognostic signature by bulk/single-cell sequencing data integration

BACKGROUND

Identifying cell subpopulations conferring unfavorable prognosis in cancer holds clinical significance. Here, we sought to identify prognostic cell subsets and develop a novel, prognostic signature for head neck squamous cell carcinoma (HNSCC).

METHODS

Highly prognostic cell subpopulations in HNSCC were identified by integrating single-cell and bulk transcriptomic datasets. The prognostic signature and nomogram were developed by least absolute shrinkage and selection operator and multivariate Cox regression analyses based on significantly upregulated genes in this specific cell subpopulation, respectively. The qRT-PCR experiments were utilized for independent validation in our patient cohort.

RESULTS

A specific cancer cell subset associated with unfavorable prognoses was identified. Functional dissections revealed that its transcriptional programs were significantly enriched in E2F, epithelial-to-mesenchymal transition, and glycolysis. A novel prognostic signature comprising six genes was developed and further validated. Risk scores based on qRT-PCR data robustly stratified patients into subgroups with distinct prognoses. A nomogram integrated from this signature and clinical stage had superior performance.

CONCLUSION

Our model derived from integrative analyses of single-cell and bulk RNA-sequencing is a novel, robust prognostic biomarker for HNSCC.

Spatial Transcriptomic and Metabolomic Landscapes of Oral Submucous Fibrosis-Derived Oral Squamous Cell Carcinoma and its Tumor Microenvironment

In South and Southeast Asia, the habit of chewing betel nuts is prevalent, which leads to oral submucous fibrosis (OSF). OSF is a well-established precancerous lesion, and a portion of OSF cases eventually progress to oral squamous cell carcinoma (OSCC). However, the specific molecular mechanisms underlying the malignant transformation of OSCC from OSF are poorly understood. In this study, the leading-edge techniques of Spatial Transcriptomics (ST) and Spatial Metabolomics (SM) are integrated to obtain spatial location information of cancer cells, fibroblasts, and immune cells, as well as the transcriptomic and metabolomic landscapes in OSF-derived OSCC tissues. This work reveals for the first time that some OSF-derived OSCC cells undergo partial epithelial-mesenchymal transition (pEMT) within the in situ carcinoma (ISC) region, eventually acquiring fibroblast-like phenotypes and participating in collagen deposition. Complex interactions among epithelial cells, fibroblasts, and immune cells in the tumor microenvironment are demonstrated. Most importantly, significant metabolic reprogramming in OSF-derived OSCC, including abnormal polyamine metabolism, potentially playing a pivotal role in promoting tumorigenesis and immune evasion is discovered. The ST and SM data in this study shed new light on deciphering the mechanisms of OSF-derived OSCC. The work also offers invaluable clues for the prevention and treatment of OSCC.

Promoting cardiomyocyte proliferation for myocardial regeneration in large mammals

From molecular and cellular perspectives, heart failure is caused by the loss of cardiomyocytes-the fundamental contractile units of the heart. Because mammalian cardiomyocytes exit the cell cycle shortly after birth, the cardiomyocyte damage induced by myocardial infarction (MI) typically leads to dilatation of the left ventricle (LV) and often progresses to heart failure. However, recent findings indicate that the hearts of neonatal pigs completely regenerated the cardiomyocytes that were lost to MI when the injury occurred on postnatal day 1 (P1). This recovery was accompanied by increases in the expression of markers for cell-cycle activity in cardiomyocytes. These results suggest that the repair process was driven by cardiomyocyte proliferation. This review summarizes findings from recent studies that found evidence of cardiomyocyte proliferation in 1) the uninjured hearts of newborn pigs on P1, 2) neonatal pig hearts after myocardial injury on P1, and 3) the hearts of pigs that underwent apical resection surgery (AR) on P1 followed by MI on postnatal day 28 (P28). Analyses of cardiomyocyte single-nucleus RNA sequencing data collected from the hearts of animals in these three experimental groups, their corresponding control groups, and fetal pigs suggested that although the check-point regulators and other molecules that direct cardiomyocyte cell-cycle progression and proliferation in fetal, newborn, and postnatal pigs were identical, the mechanisms that activated cardiomyocyte proliferation in response to injury may differ from those that regulate cardiomyocyte proliferation during development.

Identification of cuproptosis-related lncRNAs with the significance in prognosis and immunotherapy of oral squamous cell carcinoma

Cuproptosis, a recently characterized programmed cell death mechanism, has emerged as a potential contributor to tumorigenesis, metastasis, and immune modulation. Long non-coding RNAs (lncRNAs) have demonstrated diverse regulatory roles in cancer and hold promise as biomarkers. However, the involvement and prognostic significance of cuproptosis-related lncRNAs (CRLs) in oral squamous cell carcinoma (OSCC) remain poorly understood. Based on TCGA-OSCC data, we integrated single-sample gene set enrichment analysis (ssGSEA), the LASSO algorithm, and the tumor immune dysfunction and exclusion (TIDE) algorithm. We identified 11 CRLs through differential expression, Spearman correlation, and univariate Cox regression analyses. Two distinct CRL-related subtypes were unveiled, delineating divergent survival patterns, tumor microenvironments (TME), and mutation profiles. A robust CRL-based signature (including AC107027.3, AC008011.2, MYOSLID, AC005785.1, AC019080.5, AC020558.2, AC025265.1, FAM27E3, and LINC02367) prognosticated OSCC outcomes, immunotherapy responses, and anti-tumor strategies. Superior predictive power compared to other lncRNA models was demonstrated. Functional assessments confirmed the influence of FAM27E3, LINC02367, and MYOSLID knockdown on OSCC cell behaviors. Remarkably, the CRLs-based signature maintained stability across OSCC patient subgroups, underscoring its clinical potential for survival prediction. This study elucidates CRLs' roles in TME of OSCC and establishes a potential signature for precision therapy.

Integrated analysis of single-cell RNA-seq and bulk RNA-seq unravels the molecular feature of M2 macrophages of head and neck squamous cell carcinoma

The connection between head and neck squamous cell carcinoma (HNSC) and M2 tumour-associated macrophages is not yet fully understood. We gathered gene expression profiles and clinical data from HNSC patients in the TCGA database. Using Consensus Clustering, we categorized these patients into M2 macrophage-related clusters. We developed a M2 macrophage-related signature (MRS) through statistical analyses. Additionally, we assessed gene expression in HNSC cells using single-cell sequencing data (GSE139324). We identified three distinct M2 macrophage-related clusters in HNSC, each with different prognostic outcomes and immune characteristics. Patients with different MRS profiles exhibited variations in immune infiltration, genetic mutations and prognosis. FCGR2A may play a role in creating an immunosuppressive tumour microenvironment and could potentially serve as a therapeutic target for HNSC. Our study demonstrated that M2 macrophage-related genes significantly impact the development and progression of HNSC. The M2 macrophage-related model offered a more comprehensive assessment of HNSC patient prognosis, genetic mutations and immune features. FCGR2A was implicated in immunosuppressive microenvironments and may hold promise for the development of novel immunotherapeutic strategies for HNSC.

MDIC3: Matrix decomposition to infer cell-cell communication

Crosstalk among cells is vital for maintaining the biological function and intactness of systems. Most existing methods for investigating cell-cell communications are based on ligand-receptor (L-R) expression, and they focus on the study between two cells. Thus, the final communication inference results are particularly sensitive to the completeness and accuracy of the prior biological knowledge. Because existing L-R research focuses mainly on humans, most existing methods can only examine cell-cell communication for humans. As far as we know, there is currently no effective method to overcome this species limitation. Here, we propose MDIC3 (matrix decomposition to infer cell-cell communication), an unsupervised tool to investigate cell-cell communication in any species, and the results are not limited by specific L-R pairs or signaling pathways. By comparing it with existing methods for the inference of cell-cell communication, MDIC3 obtained better performance in both humans and mice.

Single-Cell Transcriptomics Reveals the Heterogeneity of the Immune Landscape of IDH-Wild-Type High-Grade Gliomas

Isocitrate dehydrogenase (IDH)-wild-type (WT) high-grade gliomas, especially glioblastomas, are highly aggressive and have an immunosuppressive tumor microenvironment. Although tumor-infiltrating immune cells are known to play a critical role in glioma genesis, their heterogeneity and intercellular interactions remain poorly understood. In this study, we constructed a single-cell transcriptome landscape of immune cells from tumor tissue and matching peripheral blood mononuclear cells (PBMC) from IDH-WT high-grade glioma patients. Our analysis identified two subsets of tumor-associated macrophages (TAM) in tumors with the highest protumorigenesis signatures, highlighting their potential role in glioma progression. We also investigated the T-cell trajectory and identified the aryl hydrocarbon receptor (AHR) as a regulator of T-cell dysfunction, providing a potential target for glioma immunotherapy. We further demonstrated that knockout of AHR decreased chimeric antigen receptor (CAR) T-cell exhaustion and improved CAR T-cell antitumor efficacy both in vitro and in vivo. Finally, we explored intercellular communication mediated by ligand-receptor interactions within the tumor microenvironment and PBMCs and revealed the unique cellular interactions present in the tumor microenvironment. Taken together, our study provides a comprehensive immune landscape of IDH-WT high-grade gliomas and offers potential drug targets for glioma immunotherapy.

Pan-cancer scRNA-seq analysis reveals immunological and diagnostic significance of the peripheral blood mononuclear cells

Peripheral blood mononuclear cells (PBMCs) reflect systemic immune response during cancer progression. However, a comprehensive understanding of the composition and function of PBMCs in cancer patients is lacking, and the potential of these features to assist cancer diagnosis is also unclear. Here, the compositional and status differences between cancer patients and healthy donors in PBMCs were investigated by single-cell RNA sequencing (scRNA-seq), involving 262,025 PBMCs from 68 cancer samples and 14 healthy samples. We observed an enhanced activation and differentiation of most immune subsets in cancer patients, along with reduction of naïve T cells, expansion of macrophages, impairment of NK cells and myeloid cells, as well as tumor promotion and immunosuppression. Based on characteristics including differential cell type abundances and/or hub genes identified from weight gene co-expression network analysis (WGCNA) modules of each major cell type, we applied logistic regression to construct cancer diagnosis models. Furthermore, we found that the above models can distinguish cancer patients and healthy donors with high sensitivity. Our study provided new insights into using the features of PBMCs in non-invasive cancer diagnosis.

Natural killer cell-based signature: Prognostic analysis in head and neck squamous cell carcinoma

BACKGROUND

Head and neck squamous cell carcinoma (HNSC) is a challenging cancer with significant clinical implications. Natural killer (NK) cells have emerged as important players in tumor immunosurveillance, yet their role and potential as prognostic biomarkers in HNSC remain unclear.

METHODS

Quantitative analysis using multiple algorithms identified FCRL1, KIR3DL2 and ZNF541 as molecules significantly associated with local NK cell infiltration and patient survival. A prognostic model based on these molecules demonstrated robust predictive performance.

RESULTS

Analysis of high- and low-risk patient groups revealed distinct differences in the tumor microenvironment, indicating an inhibitory immune microenvironment in high-risk patients. Notably, low-risk patients exhibited potential sensitivity to immunotherapy and showed favorable responses to specific drugs such as axitinib, methotrexate, rapamycin and vorinostat. NK cells, important effectors of the innate immune response, were found to play a crucial role in HNSC immunity. The present study provides valuable insights into the correlation between FCRL1, KIR3DL2, ZNF541 and NK cell infiltration, paving the way for future investigations into their roles in HNSC. Activation of NOTCH signaling, MYC targets, DNA repair, E2F targets, epithelial-mesenchymal transition, G2M checkpoint and mitotic spindle pathways in high-risk patients suggests their involvement in disease progression and poor prognosis.

CONCLUSIONS

The present study reveals the significance of NK cells in HNSC and their potential as prognostic biomarkers. The CFKZ score offers a promising approach for predicting patient outcomes and guiding personalized treatment decisions in HNSC. These findings contribute to our understanding of HNSC immunobiology and hold implications for precision medicine in HNSC management.

Immune escape of head and neck cancer mediated by the impaired MHC-I antigen presentation pathway

Tumor immune evasion is a hallmark of Head and Neck Cancers. The advent of immune checkpoint inhibitors (ICIs) in the first-line setting has transformed the management of these tumors. Unfortunately, the response rate of Head and Neck Squamous Cell Carcinomas (HNSCC) to ICIs is below 15%, regardless of the human papillomavirus (HPV) status, which might be partially related with impaired antigen presentation machinery (APM). Mechanistically, HNSCC cells are usually defective in the expression of MHC-I associated APM, while this transcriptional pathway is critical for the activation of tumor-killing effector T-cells. To specifically illuminate the phenomenon and seek for therapeutic strategies, this review summarizes the most recently identified role of genetic and functional dysregulation of the MHC-I pathway, specifically through changes at the genetic, epigenetic, post-transcriptional, and post-translational levels, which substantially contributes to HNSCC immune escape and ICI resistance. Several treatment modalities can be potentially exploited to restore APM signaling in tumors, which improves anti-tumor immunity through the activation of interferons, vaccines or rimantadine against HPV and the inhibition of EGFR, SHP-2, PI3K and MEK. Additionally, the combinatorial use of radiotherapy or cytotoxic agents with ICIs can synergize to potentiate APM signaling. Future directions would include further dissection of MHC-I related APM signaling in HNSCC and whether reversing this inhibition in combination with ICIs would elicit a more robust immune response leading to improved response rates in HNSCC. Therapeutic approaches to restore the MHC-I antigen presentation machinery in Head and Neck Cancer. (Red color texts represent the according strategies and the outcomes).

HPV-related diseases in male patients: an underestimated conundrum

PURPOSE

Human papillomavirus (HPV) infection is the most common sexually transmitted disease, in males and females worldwide. While the role of HPV in female diseases is well known and largely studied, males have negligibly been included in these programs, also because the proportion of women suffering and dying from HPV-related diseases is much larger than men. The aim of this review is to focus on HPV-related diseases in male patients.

METHODS

We performed a literature analysis on the electronic database PubMed. We considered randomized trials, observational and retrospective studies, original articles having as topic the relationship between HPV male infection and the following items: oral, anal penile cancers, warts, condylomas, male infertility, altered sperm parameters, anti-sperm antibodies (ASA). We also included experimental in vitro studies focused on the effects of HPV infection on oocyte fertilization, blastocyst development, and trophoblastic cell invasiveness. In addition, studies describing the adjuvant administration of the HPV vaccination as a possible strategy to promote HPV clearance from semen in infected males were included.

RESULTS

Regarding head and neck HPV-related diseases, the most important non-neoplastic disease is recurrent respiratory papillomatosis (RRP). Regarding neoplastic diseases, the proportion of head and neck cancers attributable to HPV has increased dramatically worldwide. In addition, nowadays, it is thought that half of head and neck squamous cell carcinomas (HNSCCs) cases in the United States are caused by infection with high-risk HPV. HPV is noteworthy in andrological practice too. It was described as having a high HPV prevalence, ranging between 50 and 70%, in male penile shaft, glans penis/coronal sulcus, semen as well as in scrotal, perianal, and anal regions. Moreover, in male patients, HPV infection has been associated, among other diseases, with penile cancers. HPV semen infection has been reported in about 10% in men from the general population and about 16% in men with unexplained infertility, although these data seem widely underestimated according to clinical experience. In particular, HPV semen infection seems to be most related to asthenozoospermia and to anti-sperm antibodies (ASAs).

CONCLUSIONS

HPV infection represents a health problem with a detrimental social and public impact. Despite this evidence, little has been done to date to widely promote vaccination among young males.

Spatial transcriptomics reveals that metabolic characteristics define the tumor immunosuppression microenvironment via iCAF transformation in oral squamous cell carcinoma

Tumor progression is closely related to tumor tissue metabolism and reshaping of the microenvironment. Oral squamous cell carcinoma (OSCC), a representative hypoxic tumor, has a heterogeneous internal metabolic environment. To clarify the relationship between different metabolic regions and the tumor immune microenvironment (TME) in OSCC, Single cell (SC) and spatial transcriptomics (ST) sequencing of OSCC tissues were performed. The proportion of TME in the ST data was obtained through SPOTlight deconvolution using SC and GSE103322 data. The metabolic activity of each spot was calculated using scMetabolism, and k-means clustering was used to classify all spots into hyper-, normal-, or hypometabolic regions. CD4T cell infiltration and TGF-β expression is higher in the hypermetabolic regions than in the others. Through CellPhoneDB and NicheNet cell-cell communication analysis, it was found that in the hypermetabolic region, fibroblasts can utilize the lactate produced by glycolysis of epithelial cells to transform into inflammatory cancer-associated fibroblasts (iCAFs), and the increased expression of HIF1A in iCAFs promotes the transcriptional expression of CXCL12. The secretion of CXCL12 recruits regulatory T cells (Tregs), leading to Treg infiltration and increased TGF-β secretion in the microenvironment and promotes the formation of a tumor immunosuppressive microenvironment. This study delineates the coordinate work axis of epithelial cells-iCAFs-Tregs in OSCC using SC, ST and TCGA bulk data, and highlights potential targets for therapy.

Association between spatial distribution of leukocyte subsets and clinical presentation of head and neck squamous cell carcinoma

Background

Interactions between tumor cells and cells in the microenvironment contribute to tumor development and metastasis. The spatial arrangement of individual cells in relation to each other influences the likelihood of whether and how these cells interact with each other.

Methods

This study investigated the effect of spatial distribution on the function of leukocyte subsets in the microenvironment of human head and neck squamous cell carcinoma (HNSCC) using multiplex immunohistochemistry (IHC). Leukocyte subsets were further classified based on analysis of two previously published HNSCC single-cell RNA datasets and flow cytometry (FC).

Results

IHC revealed distinct distribution patterns of leukocytes differentiated by CD68 and CD163. While CD68hiCD163lo and CD68hiCD163hi cells accumulated near tumor sites, CD68loCD163hi cells were more evenly distributed in the tumor stroma. PD-L1hi and PD-1hi cells accumulated predominantly around tumor sites. High cell density of PD-L1hi CD68hiCD163hi cells or PD-1hi T cells near the tumor site correlated with improved survival. FC and single cell RNA revealed high variability within the CD68/CD163 subsets. CD68hiCD163lo and CD68hiCD163hi cells were predominantly macrophages (MΦ), whereas CD68loCD163hi cells appeared to be predominantly dendritic cells (DCs). Differentiation based on CD64, CD80, CD163, and CD206 revealed that TAM in HNSCC occupy a broad spectrum within the classical M1/M2 polarization. Notably, the MΦ subsets expressed predominantly CD206 and little CD80. The opposite was observed in the DC subsets.

Conclusion

The distribution patterns and their distinct interactions via the PD-L1/PD-1 pathway suggest divergent roles of CD68/CD163 subsets in the HNSCC microenvironment. PD-L1/PD-1 interactions appear to occur primarily between specific cell types close to the tumor site. Whether PD-L1/PD-1 interactions have a positive or negative impact on patient survival appears to depend on both the spatial localization and the entity of the interacting cells. Co-expression of other markers, particularly CD80 and CD206, supports the hypothesis that CD68/CD163 IHC subsets have distinct functions. These results highlight the association between spatial leukocyte distribution patterns and the clinical presentation of HNSCC.

B cells in head and neck squamous cell carcinoma: current opinion and novel therapy

Head and neck squamous cell carcinoma (HNSCC) is a common malignant tumour. Despite advancements in surgery, radiotherapy and chemotherapy, which have improved the prognosis of most patients, a subset of patients with poor prognoses still exist due to loss of surgical opportunities, postoperative recurrence, and metastasis, among other reasons. The tumour microenvironment (TME) is a complex organization composed of tumour, stromal, and endothelial cells. Communication and interaction between tumours and immune cells within the TME are increasingly being recognized as pivotal in inhibiting or promoting tumour development. Previous studies on T cells in the TME of HNSCC have yielded novel therapeutic possibilities. However, the function of B cells, another adaptive immune cell type, in the TME of HNSCC patients has yet to be determined. Recent studies have revealed various distinct subtypes of B cells and tertiary lymphoid structures (TLSs) in the TME of HNSCC patients, which are believed to impact the efficacy of immune checkpoint inhibitors (ICIs). Therefore, this paper focuses on B cells in the TME to explore potential directions for future immunotherapy for HNSCC.

Lymph node and tumor-associated PD-L1+ macrophages antagonize dendritic cell vaccines by suppressing CD8+ T cells

Current immunotherapies provide limited benefits against T cell-depleted tumors, calling for therapeutic innovation. Using multi-omics integration of cancer patient data, we predict a type I interferon (IFN) responseHIGH state of dendritic cell (DC) vaccines, with efficacious clinical impact. However, preclinical DC vaccines recapitulating this state by combining immunogenic cancer cell death with induction of type I IFN responses fail to regress mouse tumors lacking T cell infiltrates. Here, in lymph nodes (LNs), instead of activating CD4+/CD8+ T cells, DCs stimulate immunosuppressive programmed death-ligand 1-positive (PD-L1+) LN-associated macrophages (LAMs). Moreover, DC vaccines also stimulate PD-L1+ tumor-associated macrophages (TAMs). This creates two anatomically distinct niches of PD-L1+ macrophages that suppress CD8+ T cells. Accordingly, a combination of PD-L1 blockade with DC vaccines achieves significant tumor regression by depleting PD-L1+ macrophages, suppressing myeloid inflammation, and de-inhibiting effector/stem-like memory T cells. Importantly, clinical DC vaccines also potentiate T cell-suppressive PD-L1+ TAMs in glioblastoma patients. We propose that a multimodal immunotherapy and vaccination regimen is mandatory to overcome T cell-depleted tumors.

Comparative analysis of immune infiltrates in head and neck cancers across anatomical sites

BACKGROUND

The response rate to immune checkpoint inhibitors targeting programmed cell death 1 (PD-1) receptor is 13%-18% for patients with recurrent or metastatic head and neck squamous cell carcinoma (HNSCC). Detailed understanding of the tumor immune microenvironment (TIME) is crucial in order to explain and improve this response rate. HNSCCs arise at various anatomical locations including the oral cavity, hypopharynx, larynx and oropharynx. Studies directly comparing immune infiltration between anatomical sites are scarce. Since the distinct locations could drive deviating microenvironments, we questioned whether the immune composition varies across these HNSCC sites.

METHODS

Here, we characterized the TIME of 76 fresh tumor specimens using flow cytometry and performed single-cell RNA-sequencing on nine head and neck tumor samples.

RESULTS

We found major differences in the composition of the TIME between patients. When comparing anatomical sites: tumors originating from the oral cavity had higher T cell infiltrates than tumors from other anatomical sites. The percentage of tumor-infiltrating T-lymphocytes positive for the immune checkpoint PD-1 varied considerably between patients, with the highest fraction of PD-1+ T cells found in larynx squamous cell carcinomas (SCCs). While we had hypothesized that the anatomical sites of tumor origin would drive sample clustering, our data showed that the type of TIME was more dominant and was particularly driven by the fraction of T cells positive for PD-1. Moreover, a high proportion of PD-1+ CD8+ T cells associated with an improved overall survival. Using single-cell RNA-sequencing, we observed that PD-1 expression was highest in the CD8-ENTPD1 tissue resident memory T cell/exhausted T cell and CD4-CXCL13 type 1 T helper cell clusters.

CONCLUSIONS

We found that oral cavity SCCs had the highest frequencies of T cells. We also observed considerable interpatient heterogeneity for PD-1 on T cells, with noticeably higher frequencies of PD-1+ CD4+ T helper cells in larynx SCCs. Within the entire cohort, a higher fraction of CD8+ T cells positive for PD-1 was linked to improved overall survival. Whether the fraction of PD-1+ T cells within the TIME enables immune checkpoint inhibitor response prediction for patients with head and neck cancer remains to be determined.

SCAR: Single-cell and Spatially-resolved Cancer Resources

Advances in sequencing and imaging technologies offer a unique opportunity to unravel cell heterogeneity and develop new immunotherapy strategies for cancer research. There is an urgent need for a resource that effectively integrates a vast amount of transcriptomic profiling data to comprehensively explore cancer tissue heterogeneity and the tumor microenvironment. In this context, we developed the Single-cell and Spatially-resolved Cancer Resources (SCAR) database, a combined tumor spatial and single-cell transcriptomic platform, which is freely accessible at http://8.142.154.29/SCAR2023 or http://scaratlas.com. SCAR contains spatial transcriptomic data from 21 tumor tissues and single-cell transcriptomic data from 11 301 352 cells encompassing 395 cancer subtypes and covering a wide variety of tissues, organoids, and cell lines. This resource offers diverse functional modules to address key cancer research questions at multiple levels, including the screening of tumor cell types, metabolic features, cell communication and gene expression patterns within the tumor microenvironment. Moreover, SCAR enables the analysis of biomarker expression patterns and cell developmental trajectories. SCAR also provides a comprehensive analysis of multi-dimensional datasets based on 34 state-of-the-art omics techniques, serving as an essential tool for in-depth mining and understanding of cell heterogeneity and spatial location. The implications of this resource extend to both cancer biology research and cancer immunotherapy development.

CellCommuNet: an atlas of cell-cell communication networks from single-cell RNA sequencing of human and mouse tissues in normal and disease states

Cell-cell communication, as a basic feature of multicellular organisms, is crucial for maintaining the biological functions and microenvironmental homeostasis of cells, organs, and whole organisms. Alterations in cell-cell communication contribute to many diseases, including cancers. Single-cell RNA sequencing (scRNA-seq) provides a powerful method for studying cell-cell communication by enabling the analysis of ligand-receptor interactions. Here, we introduce CellCommuNet (http://www.inbirg.com/cellcommunet/), a comprehensive data resource for exploring cell-cell communication networks in scRNA-seq data from human and mouse tissues in normal and disease states. CellCommuNet currently includes 376 single datasets from multiple sources, and 118 comparison datasets between disease and normal samples originating from the same study. CellCommuNet provides information on the strength of communication between cells and related signalling pathways and facilitates the exploration of differences in cell-cell communication between healthy and disease states. Users can also search for specific signalling pathways, ligand-receptor pairs, and cell types of interest. CellCommuNet provides interactive graphics illustrating cell-cell communication in different states, enabling differential analysis of communication strength between disease and control samples. This comprehensive database aims to be a valuable resource for biologists studying cell-cell communication networks.

Heterogeneity and subtypes of CD4+ regulatory T cells: implications for tumor therapy

In the conventional view, CD4+ regulatory T cell (Treg) represents a subset of lymphocytes that involve the perception and negative regulation of the immune response. CD4+Treg plays an important role in the maintenance of immune homeostasis and immune tolerance. However, recent studies have revealed that CD4+Treg do not suppress the immune response in some diseases, but promote inflammatory injury or inhibit tissue remodeling, suggesting the functional heterogeneity of CD4+Treg. Their involvement in tumor pathogenesis is more complex than previously understood. This article reviews the relevant research on the heterogeneity of CD4+Treg, subtype classification, and their relationship with tumor therapy.

Deciphering head and neck cancer microenvironment: Single-cell and spatial transcriptomics reveals human papillomavirus-associated differences

Human papillomavirus (HPV) is a major causative factor of head and neck squamous cell carcinoma (HNSCC), and the incidence of HPV- associated HNSCC is increasing. The role of tumor microenvironment in viral infection and metastasis needs to be explored further. We studied the molecular characteristics of primary tumors (PTs) and lymph node metastatic tumors (LNMTs) by stratifying them based on their HPV status. Eight samples for single-cell RNA profiling and six samples for spatial transcriptomics (ST), composed of matched primary tumors (PT) and lymph node metastases (LNMT), were collected from both HPV- negative (HPV- ) and HPV-positive (HPV+ ) patients. Using the 10x Genomics Visium platform, integrative analyses with single-cell RNA sequencing were performed. Intracellular and intercellular alterations were analyzed, and the findings were confirmed using experimental validation and publicly available data set. The HPV+ tissues were composed of a substantial amount of lymphoid cells regardless of the presence or absence of metastasis, whereas the HPV- tissue exhibited remarkable changes in the number of macrophages and plasma cells, particularly in the LNMT. From both single-cell RNA and ST data set, we discovered a central gene, pyruvate kinase muscle isoform 1/2 (PKM2), which is closely associated with the stemness of cancer stem cell-like populations in LNMT of HPV- tissue. The consistent expression was observed in HPV- HNSCC cell line and the knockdown of PKM2 weakened spheroid formation ability. Furthermore, we found an ectopic lymphoid structure morphology and clinical effects of the structure in ST slide of the HPV+ patients and verified their presence in tumor tissue using immunohistochemistry. Finally, the ephrin-A (EPHA2) pathway was detected as important signals in angiogenesis for HPV- patients from single-cell RNA and ST profiles, and knockdown of EPHA2 declined the cell migration. Our study described the distinct cellular composition and molecular alterations in primary and metastatic sites in HNSCC patients based on their HPV status. These results provide insights into HNSCC biology in the context of HPV infection and its potential clinical implications.

Mapping the landscape of HPV integration and characterising virus and host genome interactions in HPV-positive oropharyngeal squamous cell carcinoma

BACKGROUND

Human papillomavirus (HPV) integration into the host genome is an important factor in HPV(+)OPSCC carcinogenesis, in conjunction with HPV oncoproteins E6/E7. However, a well-studied investigation about virus-host interaction still needs to be completed. Our objective is to characterise HPV integration to investigate potential mechanisms of tumourigenesis independent of E6/E7 oncoproteins.

MATERIALS AND METHODS

High-throughput viral integration detection was performed on 109 HPV(+)OPSCC tumours with relevant clinicopathological information. Of these tumours, 38 tumours underwent targeted gene sequencing, 29 underwent whole exome sequencing and 26 underwent RNA sequencing.

RESULTS

HPV integration was detected in 94% of tumours (with a mean integration count of 337). Tumours occurring at the tonsil/oropharyngeal wall that exhibit higher PD-L1 expression demonstrated increased integration sites (p = .024). HPV exhibited a propensity for integration at genomic sites located within specific fragile sites (FRA19A) or genes associated with functional roles such as cell proliferation and differentiation (PTEN, AR), immune evasion (CD274) and glycoprotein biosynthesis process (FUT8). The viral oncogenes E2, E4, E6 and E7 tended to remain intact. HPV fragments displayed enrichment within host copy number variation (CNV) regions. However, insertions into genes related to altered homologous recombination repair were infrequent. Genes with integration had distinct expression levels. Fifty-nine genes whose expression level was affected by viral integration were identified, for example, EPHB1, which was reported to be involved in cellular protein metabolic process.

CONCLUSIONS

HPV can promote oncogenesis through recurrent integration into functional host genome regions, leading to subsequent genomic aberrations and gene expression disruption. This study characterises viral integrations and virus-host interactions, enhancing our understanding of HPV-related carcinogenesis mechanisms.

Immune landscape in molecular subtypes of human papillomavirus-negative head and neck cancer

Head and neck squamous cell carcinomas (HNSCC) remain a poorly understood disease clinically and immunologically. HPV is a known risk factor of HNSCC associated with better outcome, whereas HPV-negative HNSCC are more heterogeneous in outcome. Gene expression signatures have been developed to classify HNSCC into four molecular subtypes (classical, basal, mesenchymal, and atypical). However, the molecular underpinnings of treatment response and the immune landscape for these molecular subtypes are largely unknown. Herein, we described a comprehensive immune landscape analysis in three independent HNSCC cohorts (>700 patients) using transcriptomics data. We assigned the HPV- HNSCC patients into these four molecular subtypes and characterized the tumor microenvironment using deconvolution method. We determined that atypical and mesenchymal subtypes have greater immune enrichment and exhibit a T-cell exhaustion phenotype, compared to classical and basal subtypes. Further analyses revealed different B cell maturation and antibody isotypes enrichment patterns, and distinct immune microenvironment crosstalk in the atypical and mesenchymal subtypes. Taken together, our study suggests that treatments that enhances B cell activity may benefit patients with HNSCC of the atypical subtypes. The rationale can be utilized in the design of future precision immunotherapy trials based on the molecular subtypes of HPV- HNSCC.

Identification of TTC21A as a Potential Prognostic Marker in Head and Neck Squamous Cell Carcinoma: In Silico Analysis

BACKGROUND/AIM

Tetratricopeptide repeat domain 21A (TTC21A) plays a crucial role in ciliary function and has been associated with various pathogenic processes, including carcinogenesis. However, its role in head and neck squamous cell carcinoma (HNSCC) has not been elucidated.

MATERIALS AND METHODS

Based on the sequencing and microarray data of HNSCC from publicly available databases, the expression of TTC21A was compared between different subgroups based on clinical and molecular parameters. The survival analysis and regression analysis were conducted using the Kaplan-Meier method and the Cox method, respectively. Functional analysis was performed by the Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Ontology (GO), and gene set enrichment analysis (GSEA) tools. Immune infiltration analysis was performed based on the expression of TTC21A.

RESULTS

TTC21A decreased in tumor tissues and was associated with N stage, histologic grade, HPV infection, and TP53 mutation in HNSCC. TTC21A was an independent indicator of overall survival for patients with HNSCC. A high level of TTC21A expression indicated a favorable prognosis. The TTC21A expression level was involved with immune-related signaling regulation, immune-related gene expression, and immune cell infiltration. TTC21A expression was potent in predicting immunotherapeutic benefits.

CONCLUSION

TTC21A, as a potential predictor of favorable outcomes and immunotherapy response for HNSCC, is related to immune-related signaling regulation, immune-related gene expression, and immune cell infiltration.