CXCL9: evidence and contradictions for its role in tumor progression

Chemokines and Cytokines in Immunotherapy of Melanoma and Other Tumors: From Biomarkers to Therapeutic Targets

Chemokines and cytokines represent an emerging field of immunotherapy research. They are responsible for the crosstalk and chemoattraction of immune cells and tumor cells. For instance, CXCL9/10/11 chemoattract effector CD8+ T cells to the tumor microenvironment, making an argument for their promising role as biomarkers for a favorable outcome. The cytokine Interleukin-15 (IL-15) can promote the chemokine expression of CXCR3 ligands but also XCL1, contributing to an important DC-T cell interaction. Recruited cytotoxic T cells can be clonally expanded by IL-2. Delivering or inducing these chemokines and cytokines can result in tumor shrinkage and might synergize with immune checkpoint inhibition. In addition, blocking specific chemokine and cytokine receptors such as CCR2, CCR4 or Il-6R can reduce the recruitment of tumor-associated macrophages (TAMs), myeloid-derived suppressor cells (MDSCs) or regulatory T cells (Tregs). Efforts to target these chemokines and cytokines have the potential to personalize cancer immunotherapy further and address patients that are not yet responsive because of immune cell exclusion. Targeting cytokines such as IL-6 and IL-15 is currently being evaluated in clinical trials in combination with immune checkpoint-blocking antibodies for the treatment of metastatic melanoma. The improved overall survival of melanoma patients might outweigh potential risks such as autoimmunity. However, off-target toxicity needs to be elucidated.

Tuberculous pleural effusion-induced Arg-1+ macrophage polarization contributes to lung cancer progression via autophagy signaling

BACKGROUND

The association between tuberculous fibrosis and lung cancer development has been reported by some epidemiological and experimental studies; however, its underlying mechanisms remain unclear, and the role of macrophage (MФ) polarization in cancer progression is unknown. The aim of the present study was to investigate the role of M2 Arg-1+ MФ in tuberculous pleurisy-assisted tumorigenicity in vitro and in vivo.

METHODS

The interactions between tuberculous pleural effusion (TPE)-induced M2 Arg-1+ MФ and A549 lung cancer cells were evaluated. A murine model injected with cancer cells 2 weeks after Mycobacterium bovis bacillus Calmette-Guérin pleural infection was used to validate the involvement of tuberculous fibrosis to tumor invasion.

RESULTS

Increased CXCL9 and CXCL10 levels of TPE induced M2 Arg-1+ MФ polarization of murine bone marrow-derived MФ. TPE-induced M2 Arg-1+ MФ polarization facilitated lung cancer proliferation via autophagy signaling and E-cadherin signaling in vitro. An inhibitor of arginase-1 targeting M2 Arg-1+ MФ both in vitro and in vivo significantly reduced tuberculous fibrosis-induced metastatic potential of lung cancer and decreased autophagy signaling and E-cadherin expression.

CONCLUSION

Tuberculous pleural fibrosis induces M2 Arg-1+ polarization, and M2 Arg-1+ MФ contribute to lung cancer metastasis via autophagy and E-cadherin signaling. Therefore, M2 Arg-1+ tumor associated MФ may be a novel therapeutic target for tuberculous fibrosis-induced lung cancer progression.

Identification of PANoptosis-related subtypes, construction of a prognosis signature, and tumor microenvironment landscape of hepatocellular carcinoma using bioinformatic analysis and experimental verification

Background

Hepatocellular carcinoma (HCC) is one of the most lethal malignancies worldwide. PANoptosis is a recently unveiled programmed cell death pathway, Nonetheless, the precise implications of PANoptosis within the context of HCC remain incompletely elucidated.

Methods

We conducted a comprehensive bioinformatics analysis to evaluate both the expression and mutation patterns of PANoptosis-related genes (PRGs). We categorized HCC into two clusters and identified differentially expressed PANoptosis-related genes (DEPRGs). Next, a PANoptosis risk model was constructed using LASSO and multivariate Cox regression analyses. The relationship between PRGs, risk genes, the risk model, and the immune microenvironment was studies. In addition, drug sensitivity between high- and low-risk groups was examined. The expression profiles of these four risk genes were elucidate by qRT-PCR or immunohistochemical (IHC). Furthermore, the effect of CTSC knock down on HCC cell behavior was verified using in vitro experiments.

Results

We constructed a prognostic signature of four DEPRGs (CTSC, CDCA8, G6PD, and CXCL9). Receiver operating characteristic curve analyses underscored the superior prognostic capacity of this signature in assessing the outcomes of HCC patients. Subsequently, patients were stratified based on their risk scores, which revealed that the low-risk group had better prognosis than those in the high-risk group. High-risk group displayed a lower Stromal Score, Immune Score, ESTIMATE score, and higher cancer stem cell content, tumor mutation burden (TMB) values. Furthermore, a correlation was noted between the risk model and the sensitivity to 56 chemotherapeutic agents, as well as immunotherapy efficacy, in patient with. These findings provide valuable guidance for personalized clinical treatment strategies. The qRT-PCR analysis revealed that upregulated expression of CTSC, CDCA8, and G6PD, whereas downregulated expression of CXCL9 in HCC compared with adjacent tumor tissue and normal liver cell lines. The knockdown of CTSC significantly reduced both HCC cell proliferation and migration.

Conclusion

Our study underscores the promise of PANoptosis-based molecular clustering and prognostic signatures in predicting patient survival and discerning the intricacies of the tumor microenvironment within the context of HCC. These insights hold the potential to advance our comprehension of the therapeutic contribution of PANoptosis plays in HCC and pave the way for generating more efficacious treatment strategies.

Integration of single-cell RNA-seq and bulk RNA-seq data to construct and validate a cancer-associated fibroblast-related prognostic signature for patients with ovarian cancer

BACKGROUND

To establish a prognostic risk profile for ovarian cancer (OC) patients based on cancer-associated fibroblasts (CAFs) and gain a comprehensive understanding of their role in OC progression, prognosis, and therapeutic efficacy.

METHODS

Data on OC single-cell RNA sequencing (scRNA-seq) and total RNA-seq were collected from the GEO and TCGA databases. Seurat R program was used to analyze scRNA-seq data and identify CAFs clusters corresponding to CAFs markers. Differential expression analysis was performed on the TCGA dataset to identify prognostic genes. A CAF-associated risk signature was designed using Lasso regression and combined with clinicopathological variables to develop a nomogram. Functional enrichment and the immune landscape were also analyzed.

RESULTS

Five CAFs clusters were identified in OC using scRNA-seq data, and 2 were significantly associated with OC prognosis. Seven genes were selected to develop a CAF-based risk signature, primarily associated with 28 pathways. The signature was a key independent predictor of OC prognosis and relevant in predicting the results of immunotherapy interventions. A novel nomogram combining CAF-based risk and disease stage was developed to predict OC prognosis.

CONCLUSION

The study highlights the importance of CAFs in OC progression and suggests potential for innovative treatment strategies. A CAF-based risk signature provides a highly accurate prediction of the prognosis of OC patients, and the developed nomogram shows promising results in predicting the OC prognosis.

Decreased expression of TXNIP is associated with poor prognosis and immune infiltration in kidney renal clear cell carcinoma

The most prevalent and insidious type of kidney cancer is kidney clear cell carcinoma (KIRC). Thioredoxin-interacting protein (TXNIP) encodes a thioredoxin-binding protein involved in cellular energy metabolism, redox homeostasis, apoptosis induction and inflammatory responses. However, the relationship between TXNIP, immune infiltration and its prognostic value in KIRC remains unclear. Thus, the present study evaluated the potential for TXNIP as a prognostic marker in patients with KIRC. Data from The Cancer Genome Atlas were used to assess relative mRNA expression levels of TXNIP in different types of cancer. The protein expression levels of TXNIP were evaluated using the Human Protein Atlas. Enrichment analysis of genes co-expressed with TXNIP was performed to assess relevant biological processes that TXNIP may be involved in. CIBERSORT was used to predict the infiltration of 21 tumor-infiltrating immune cells (TIICs). Univariate and multivariate Cox regression analyses were used to assess the relationship between TXNIP expression and prognosis. Single-cell RNA-sequencing datasets were used to evaluate the mRNA expression levels of TXNIP in certain immune cells in KIRC. The CellMiner database was used to analyze the relationship between TXNIP mRNA expression and drug sensitivity in KIRC. The results from the present study demonstrated that TXNIP expression was significantly decreased in KIRC tissue compared with that in normal tissue, as confirmed by western blotting and reverse transcription-quantitative PCR. In addition, downregulated TXNIP expression was significantly associated with poor prognosis, a high histological grade and an advanced stage. The Cell Counting Kit-8 assay demonstrated that TXNIP overexpression significantly suppressed tumor cell proliferation. Univariate and multivariate Cox regression analyses indicated that TXNIP served as a separate prognostic factor in KIRC. Moreover, TXNIP expression was significantly correlated with the accumulation of several TIICs and its overexpression significantly downregulated the mRNA expression levels of CD25 and cytotoxic T-lymphocyte-associated protein 4, immune cell surface markers in CD4+ T lymphocytes. In conclusion, TXNIP may be used as a possible biomarker to assess unfavorable prognostic outcomes and identify immunotherapy targets in KIRC.

Single-cell colocalization analysis using a deep generative model

Analyzing colocalization of single cells with heterogeneous molecular phenotypes is essential for understanding cell-cell interactions, and cellular responses to external stimuli and their biological functions in diseases and tissues. However, existing computational methodologies identified the colocalization patterns between predefined cell populations, which can obscure the molecular signatures arising from intercellular communication. Here, we introduce DeepCOLOR, a computational framework based on a deep generative model that recovers intercellular colocalization networks with single-cell resolution by the integration of single-cell and spatial transcriptomes. Along with colocalized population detection accuracy that is superior to existing methods in simulated dataset, DeepCOLOR identified plausible cell-cell interaction candidates between colocalized single cells and segregated cell populations defined by the colocalization relationships in mouse brain tissues, human squamous cell carcinoma samples, and human lung tissues infected with SARS-CoV-2. DeepCOLOR is applicable to studying cell-cell interactions behind various spatial niches. A record of this paper's transparent peer review process is included in the supplemental information.

Blockade of C5a receptor unleashes tumor-associated macrophage antitumor response and enhances CXCL9-dependent CD8+ T cell activity

Macrophages play a crucial role in shaping the immune state within the tumor microenvironment (TME) and are often influenced by tumors to hinder antitumor immunity. However, the underlying mechanisms are still elusive. Here, we observed abnormal expression of complement 5a receptor (C5aR) in human ovarian cancer (OC), and identified high levels of C5aR expression on tumor-associated macrophages (TAMs), which led to the polarization of TAMs toward an immunosuppressive phenotype. C5aR knockout or inhibitor treatment restored TAM antitumor response and attenuated tumor progression. Mechanistically, C5aR deficiency reprogrammed macrophages from a protumor state to an antitumor state, associating with the upregulation of immune response and stimulation pathways, which in turn resulted in the enhanced antitumor response of cytotoxic T cells in a manner dependent on chemokine (C-X-C motif) ligand 9 (CXCL9). The pharmacological inhibition of C5aR also improved the efficacy of immune checkpoint blockade therapy. In patients, C5aR expression associated with CXCL9 production and infiltration of CD8+ T cells, and a high C5aR level predicted poor clinical outcomes and worse benefits from anti-PD-1 therapy. Thus, our study sheds light on the mechanisms underlying the modulation of TAM antitumor immune response by the C5a-C5aR axis and highlights the potential of targeting C5aR for clinical applications.

C-X-C Motif Chemokine Ligand 9 Correlates with Favorable Prognosis in Triple-Negative Breast Cancer by Promoting Immune Cell Infiltration

C-X-C motif chemokine ligand 9 (CXCL9) plays an important role in antitumor immunity through the recruitment, proliferation, and activation of immune cells (IC). Here, we evaluated the expression patterns of CXCL9 and programmed death-ligand 1 (PD-L1) in a cohort of 268 patients with triple-negative breast cancer (TNBC) by tissue microarray (TMA). The correlations between CXCL9 expression in ICs or tumor cells (TC) and clinicopathologic parameters, PD-L1 expression, tumor-infiltrating lymphocytes (TIL) and survival were analyzed in this cohort (n = 268). In addition, we analyzed a TNBC dataset (n = 138) from The Cancer Genome Atlas (TCGA) to identify correlation between CXCL9 expression and other immune gene expression, immune infiltration, and prognosis. The results of the TMA cohort (n = 268) showed that CXCL9 was expressed in 80.6% cases, with elevated expression levels in ICs relative to in TCs (median: 1% vs. 0%). CXCL9 expressed in ≥1% of ICs was categorized as the CXCL9-IC-positive group. CXCL9-IC expression was strongly and positively correlated with the PD-L1 expression, CD3+ TILs, CD4+ TILs, CD8+ TILs, and CD19+ TILs (all P < 0.0001). Survival analyses showed that the CXCL9-IC-positive group demonstrated prolonged disease-free survival (P = 0.038) and overall survival (P = 0.023) compared with the negative group. The analyses from TCGA cohort (n = 138) showed that elevated CXCL9 expression correlated with increased infiltration of B cells, macrophages, natural killer cells, monocytes and increased expression of immune checkpoint molecules and other CXCL family members, including CXCL10 and CXCL11. These findings confirm the regulatory role of CXCL9 in antitumor immunity and suggest a potential role in treatments involving immune checkpoint blockade.

Effect of TRAF6-knockout on gene expression and lncRNA expression in Epithelioma papulosum cyprini (EPC) cells

TRAF6 is a key immune gene that plays a significant role in toll-like receptor signal transduction and activates downstream immune genes involved in antiviral immunity in fish. To explore the role of TRAF6 in Epithelioma papulosum cyprini (EPC) cells, we knocked out the TRAF6 gene using the Clustered Regularly Interspaced Short Palindromic Repeats-Cas9 (CRISPR-Cas9) technique and then analyzed the transcriptomes of the knockout cells. In this study, we identified that 232 transcripts were differentially expressed in naive cells. Using the pipeline, we identified 381 novel lncRNAs in EPC cells, 23 of which were differentially expressed. Gene Ontology enrichment analysis demonstrated that differentially expressed genes (DEG) are implicated in various immune processes, such as neutrophil chemotaxis and mitogen-activated protein kinase binding. In addition, the KEGG pathway analysis revealed enrichment in immune-related pathways (Interleukin-17 signaling pathway, cytokine-cytokine receptor interaction, and TNF signaling pathway). Furthermore, the target genes of the differentially expressed lncRNAs were implicated in the negative regulation of interleukin-6 and tumor necrosis factor production. These results indicate that lncRNAs and protein-coding genes participate in the regulation of immune and metabolic processes in fish.

Establishment of SLC7A11-knockout mouse and its preliminary investigation in melanoma

Solute carrier family 7 member 11 (SLC7A11)/xCT is an amino acid transporter that mediates the cystine uptake and glutamate export, participates in several malignant tumors' progression. However, the role of SLC7A11 on the occurrence and development of melanoma still remains unclear. Here, the transcribed mRNA encoding for Cas9 and sgRNA targeting SLC7A11 in vitro were microinjected into zygotes, to establish the SLC7A11 knockout (KO) mice (SLC7A11-/-). Further, we conducted melanoma-bearing mice using the metastatic melanoma cell line (B16-F10) to observe the melanoma development. There was no off-target in KO mice detected by T7E1 cleavage assay. The results showed that the tumor volume of KO mice was significantly lower than that of SLC7A11+/+ (WT) mice at 8d, 10d, 12d, 14d, and 16d (P < 0.05). The tumors of WT appeared to more disorganized morphology, more unbalanced nuclear-cytoplasmic ratio, less defined boundary, and increased tumor necrosis. And after SLC7A11 deletion, the expression of CXCL9 and TLR6 were significantly up-regulated, and that of NOS2 and CCL8 were significantly down-regulated (P < 0.01). Additionally, Ki67 immunostaining revealed lower proliferating cells in the tumors of SLC7A11 KO mice compared to WT mice. In summary, the deletion of SLC7A11 significantly inhibited the development of melanoma. Our results provide direct evidence to identify SLC7A11 as a novel target for molecular therapy and prognosis judgment of melanoma.

SIMVI reveals intrinsic and spatial-induced states in spatial omics data

Spatial omics analyze gene expression and interaction dynamics in relation to tissue structure and function. However, existing methods cannot model the intrinsic and spatial-induced variation in spatial omics data, thus failing to identify true spatial interaction effects. Here, we present Spatial Interaction Modeling using Variational Inference (SIMVI), an annotation-free framework that disentangles cell intrinsic and spatial-induced latent variables for modeling gene expression in spatial omics data. SIMVI enables novel downstream analyses, such as clustering and differential expression analysis based on disentangled representations, spatial effect (SE) identification, SE interpretation, and transfer learning on new measurements / modalities. We benchmarked SIMVI on both simulated and real datasets and show that SIMVI uniquely generates highly accurate SE inferences in synthetic datasets and unveils intrinsic variation in complex real datasets. We applied SIMVI to spatial omics data from diverse platforms and tissues (MERFISH human cortex / mouse liver, Slide-seqv2 mouse hippocampus, Spatial-ATAC-RNA-seq) and revealed various region-specific and cell-type-specific spatial interactions. In addition, our experiments on MERFISH human cortex and spatial-ATAC-RNA-seq showcased SIMVI's power in identifying SEs for new samples / modalities. Finally, we applied SIMVI on a newly collected CosMx melanoma dataset. Using SIMVI, we identified immune cells associated with spatial-dependent interactions and revealed the underlying spatial variations associated with patient outcomes.

Targeting CXCL9/10/11-CXCR3 axis: an important component of tumor-promoting and antitumor immunity

Chemokines are chemotactic-competent molecules composed of a family of small cytokines, playing a key role in regulating tumor progression. The roles of chemokines in antitumor immune responses are of great interest. CXCL9, CXCL10, and CXCL11 are important members of chemokines. It has been widely investigated that these three chemokines can bind to their common receptor CXCR3 and regulate the differentiation, migration, and tumor infiltration of immune cells, directly or indirectly affecting tumor growth and metastasis. Here, we summarize the mechanism of how the CXCL9/10/11-CXCR3 axis affects the tumor microenvironment, and list the latest researches to find out how this axis predicts the prognosis of different cancers. In addition, immunotherapy improves the survival of tumor patients, but some patients show drug resistance. Studies have found that the regulation of CXCL9/10/11-CXCR3 on the tumor microenvironment is involved in the process of changing immunotherapy resistance. Here we also describe new approaches to restoring sensitivity to immune checkpoint inhibitors through the CXCL9/10/11-CXCR3 axis.

Colon cancer transcriptome

Over the last four decades, methodological innovations have continuously changed transcriptome profiling. It is now feasible to sequence and quantify the transcriptional outputs of individual cells or thousands of samples using RNA sequencing (RNA-seq). These transcriptomes serve as a connection between cellular behaviors and their underlying molecular mechanisms, such as mutations. This relationship, in the context of cancer, provides a chance to unravel tumor complexity and heterogeneity and uncover novel biomarkers or treatment options. Since colon cancer is one of the most frequent malignancies, its prognosis and diagnosis seem to be critical. The transcriptome technology is developing for an earlier and more accurate diagnosis of cancer which can provide better protectivity and prognostic utility to medical teams and patients. A transcriptome is a whole set of expressed coding and non-coding RNAs in an individual or cell population. The cancer transcriptome includes RNA-based changes. The combined genome and transcriptome of a patient may provide a comprehensive picture of their cancer, and this information is beginning to affect treatment decision-making in real-time. A full assessment of the transcriptome of colon (colorectal) cancer has been assessed in this review paper based on risk factors such as age, obesity, gender, alcohol use, race, and also different stages of cancer, as well as non-coding RNAs like circRNAs, miRNAs, lncRNAs, and siRNAs. Similarly, they have been examined independently in the transcriptome study of colon cancer.

Comparative Transcriptomic Analysis of Archival Human Vestibular Schwannoma Tissue from Patients with and without Tinnitus

Vestibular schwannoma (VS) is an intracranial tumor that commonly presents with tinnitus and hearing loss. To uncover the molecular mechanisms underlying VS-associated tinnitus, we applied next-generation sequencing (Illumina HiSeq) to formalin-fixed paraffin-embedded archival VS samples from nine patients with tinnitus (VS-Tin) and seven patients without tinnitus (VS-NoTin). Bioinformatic analysis was used to detect differentially expressed genes (DEG; i.e., ≥two-fold change [FC]) while correcting for multiple comparisons. Using RNA-seq analysis, VS-Tin had significantly lower expression of GFAP (logFC = −3.04), APLNR (logFC = −2.95), PREX2 (logFC = −1.44), and PLVAP (logFC = −1.04; all p < 0.01) vs. VS-NoTin. These trends were validated by using real-time RT-qPCR. At the protein level, immunohistochemistry revealed a trend for less PREX2 and apelin expression and greater expression of NLRP3 inflammasome and CD68-positive macrophages in VS-Tin than in VS-NoTin, suggesting the activation of inflammatory processes in VS-Tin. Functional enrichment analysis revealed that the top three protein categories—glycoproteins, signal peptides, and secreted proteins—were significantly enriched in VS-Tin in comparison with VS-NoTin. In a gene set enrichment analysis, the top pathway was allograft rejection, an inflammatory pathway that includes the MMP9, CXCL9, IL16, PF4, ITK, and ACVR2A genes. Future studies are needed to examine the importance of these candidates and of inflammation in VS-associated tinnitus.

EBV-positive pyothorax-associated lymphoma expresses CXCL9 and CXCL10 chemokines that attract cytotoxic lymphocytes via CXCR3

Epstein-Barr virus (EBV)-positive diffuse large B-cell lymphoma associated with chronic inflammation (DLBCL-CI) develops in the setting of long-standing inflammation. This type of lymphoma may have specific expression profiles of chemokines involved in the pathogenesis of DLBCL-CI. EBV-positive pyothorax-associated lymphoma (PAL) is a prototype of DLBCL-CI and represents a valuable model for the study of this disease category. Using a panel of PAL cell lines, we found that PAL cells expressed and secreted C-X-C motif chemokine ligands 9 and 10 (CXCL9 and CXCL10), the ligands of CXCR3, in contrast to EBV-negative DLBCL cell lines, which did not. Culture supernatants from PAL cell lines attracted CXCR3-expressing CD4⁺ T cells, CD8⁺ T cells, and CD56⁺ natural killer cells from human peripheral blood mononuclear cells. PAL cells injected into mice also attracted CXCR3-positive cytotoxic lymphocytes that expressed interferon-γ. The expression of CXCL9 and CXCL10 was detected in PAL tumor biopsy samples from patients, and CXCR3-positive lymphocytes were abundant in the tissue samples. Collectively, these findings suggest that CXCL9 and CXCL10 are produced by PAL cells and can elicit cytotoxic responses via CXCR3. This chemokine system is also likely to contribute to tissue necrosis, which is a signature histological feature of DLBCL-CI. Further studies are warranted to determine whether the CXCL9-CXCL10/CXCR3 axis exerts antitumor effects in DLBCL-CI.

CXCL9 correlates with antitumor immunity and is predictive of a favorable prognosis in uterine corpus endometrial carcinoma

Background

The C-X-C motif chemokine ligand-9 (CXCL9) is related to the progression of multiple neoplasms. Yet, its biological functions in uterine corpus endometrioid carcinoma (UCEC) remain shrouded in confusion. Here, we assessed the prognostic significance and potential mechanism of CXCL9 in UCEC.

Methods

Firstly, bioinformatics analysis of the public cancer database, including the Cancer Genome Atlas / the Genotype-Tissue Expression project (TCGA+ GTEx, n=552) and Gene Expression Omnibus (GEO): GSE63678 (n=7), were utilized for the CXCL9 expression-related analysis in UCEC. Then, the survival analysis of TCGA-UCEC was performed. Futher, the gene set enrichment analysis (GSEA) was carried out to reveal the potential molecular signaling pathway in UCEC associated with CXCL9 expression. Moreover, the immunohistochemistry (IHC) assay of our validation cohort (n=124) from human specimens were used to demonstrate the latent significance of CXCL9 in UCEC.

Results

The bioinformatics analysis suggested that CXCL9 expression was significantly upregulated in UCEC patients; and hyper-expression of CXCL9 was related to prolonged survival. the GSEA enrichment analysis showed various immune response-related pathways, including T/NK cell, lymphocyte activation, cytokine-cytokine receptor interaction network, and chemokine signaling pathway, mediated by CXCL9. In addition, the cytotoxic molecules (IFNG, SLAMF7, JCHAIN, NKG7, GBP5, LYZ, GZMA, GZMB, and TNF3F9) and the immunosuppressive genes (including PD-L1) were positively related to the expression of CXCL9. Further, the IHC assay indicated that the CXCL9 protein expression was mainly located in intertumoral and significantly upregulated in the UCEC patients; UCEC with high intertumoral CXCL9 cell abundance harbored an improved prognosis; a higher ratio of anti-tumor immune cells (CD4⁺, CD8⁺, and CD56⁺ cell) and PD-L1 was found in UCEC with CXCL9 high expression.

Conclusion

Overexpressed CXCL9 correlates with antitumor immunity and is predictive of a favorable prognosis in UCEC. It hinted that CXCL9 may serve as an independent prognostic biomarker or therapeutic target in UCEC patients, which augmented anti-tumor immune effects to furnish survival benefits.

Identification of Comorbidities, Genomic Associations, and Molecular Mechanisms for COVID-19 Using Bioinformatics Approaches

Several studies have been done to identify comorbidities of COVID-19. In this work, we developed an analytical bioinformatics framework to reveal COVID-19 comorbidities, their genomic associations, and molecular mechanisms accomplishing transcriptomic analyses of the RNA-seq datasets provided by the Gene Expression Omnibus (GEO) database, where normal and infected tissues were evaluated. Using the framework, we identified 27 COVID-19 correlated diseases out of 7,092 collected diseases. Analyzing clinical and epidemiological research, we noticed that our identified 27 diseases are associated with COVID-19, where hypertension, diabetes, obesity, and lung cancer are observed several times in COVID-19 patients. Therefore, we selected the above four diseases and performed assorted analyses to demonstrate the association between COVID-19 and hypertension, diabetes, obesity, and lung cancer as comorbidities. We investigated genomic associations with the cross-comparative analysis and Jaccard's similarity index, identifying shared differentially expressed genes (DEGs) and linking DEGs of COVID-19 and the comorbidities, in which we identified hypertension as the most associated illness. We also revealed molecular mechanisms by identifying statistically significant ten pathways and ten ontologies. Moreover, to understand cellular physiology, we did protein-protein interaction (PPI) analyses among the comorbidities and COVID-19. We also used the degree centrality method and identified ten biomarker hub proteins (IL1B, CXCL8, FN1, MMP9, CXCL10, IL1A, IRF7, VWF, CXCL9, and ISG15) that associate COVID-19 with the comorbidities. Finally, we validated our findings by searching the published literature. Thus, our analytical approach elicited interconnections between COVID-19 and the aforementioned comorbidities in terms of remarkable DEGs, pathways, ontologies, PPI, and biomarker hub proteins.

Physiology of chemokines in the cancer microenvironment

Chemokines are chemotactic cytokines whose canonical functions govern movement of receptor-expressing cells along chemical gradients. Chemokines are a physiological system that is finely tuned by ligand and receptor expression, ligand or receptor oligomerization, redundancy, expression of atypical receptors, and non-GPCR binding partners that cumulatively influence discrete pharmacological signaling responses and cellular functions. In cancer, chemokines play paradoxical roles in both the directed emigration of metastatic, receptor-expressing cancer cells out of the tumor as well as immigration of tumor-infiltrating immune cells that culminate in a tumor-unique immune microenvironment. In the age of precision oncology, strategies to effectively harness the power of immunotherapy requires consideration of chemokine gradients within the unique spatial topography and temporal influences with heterogeneous tumors. In this article, we review current literature on the diversity of chemokine ligands and their cellular receptors that detect and process chemotactic gradients and illustrate how differences between ligand recognition and receptor activation influence the signaling machinery that drives cellular movement into and out of the tumor microenvironment. Facets of chemokine physiology across discrete cancer immune phenotypes are contrasted to existing chemokine-centered therapies in cancer.

Tumor Microenvironment before and after Chemoradiation in Locally Advanced Rectal Cancer: Beyond PD-L1

BACKGROUND

In locally advanced rectal cancer treatment, neoadjuvant concurrent chemoradiation therapy (cCRT) is the standard of care. The tumor microenvironment (TME) is a complex entity comprising of tumor cells, immune cells and surrounding stroma and is closely associated with tumor growth and survival, response to antitumor therapies and also resistance to treatment. We aimed to assess the change in biomarkers associated with TME following standard neoadjuvant cCRT in rectal cancer.

METHODS

We accessed archival tissue from rectal cancer patients treated with neoadjuvant cCRT at Allegheny Health Network (AHN) facilities over the past 14 years. Pre-treatment and post-treatment biopsies were assayed for PD-L1, CD8+ T-cells, CXCL9, TIM-3, IDO-1, IFN-G, IL17RE, LAG-3, and OX40 in 41 patients.

RESULTS

We found statistically significant upregulation in multiple biomarkers namely CD8, IL17RE, LAG3 and OX40 post neoadjuvant cCRT and a trend towards upregulation, although not statistically significant, in biomarkers PD-L1, CXCL9, TIM-3, IDO-1 and IFN-G expression.

CONCLUSIONS

This provides a glimpse into the TME before and after neoadjuvant cCRT. We suggest that the biomarkers noted to be upregulated could be used for designing appropriate clinical trials and development of therapeutic targeted drug therapy in an effort to achieve better response to neoadjuvant therapy, increasing clinical and pathological complete response rates and improved overall outcomes.

The role of CXCR3 and its ligands in cancer

Chemokines are a class of small cytokines or signaling proteins that are secreted by cells. Owing to their ability to induce directional chemotaxis of nearby responding cells, they are called chemotactic cytokines. Chemokines and chemokine receptors have now been shown to influence many cellular functions, including survival, adhesion, invasion, and proliferation, and regulate chemokine levels. Most malignant tumors express one or more chemokine receptors. The CXC subgroup of chemokine receptors, CXCR3, is mainly expressed on the surface of activated T cells, B cells, and natural killer cells, and plays an essential role in infection, autoimmune diseases, and tumor immunity by binding to specific receptors on target cell membranes to induce targeted migration and immune responses. It is vital to treat infections, autoimmune diseases, and tumors. CXCR3 and its ligands, CXCL9, CXCL10, and CXCL11, are closely associated with the development and progression of many tumors. With the elucidation of its mechanism of action, CXCR3 is expected to become a new indicator for evaluating the prognosis of patients with tumors and a new target for clinical tumor immunotherapy. This article reviews the significance and mechanism of action of the chemokine receptor CXCR3 and its specific ligands in tumor development.

Chimeric antigen receptor T cells applied to solid tumors

Chimeric antigen receptor T cell (CAR-T) therapy is novel tumor immunotherapy that enables autologous T to express synthetic receptors to specifically recognize the surface tumor-associated antigens for exerting subsequent antitumor effects, and eliminating the resistance, metastases and recurrence of cancer. Although CAR T cells have exhibited success in eradicating hematologic malignancies, their applications to solid tumors has not yet been achieved due to obstacles such as the immune-suppressor tumor microenvironment and lack of tumor specific target antigens. In this review, we presented advancements in the development of CAR T cell therapy in solid tumors, and offered a brief summary of the challenges, as well as novel engineering and pharmaceutical interventions to overcome these barriers. Looking forward, we discussed the latest studies which are expected to reach the clinicals in the next few years, including CRISPR screens-based CAR modification and CAR T cells driven from progenitor-like T cells. Collectively, this review may inspire researchers and clinicians to develop clinical available strategies of CAR T cell therapies in solid tumor.

CXCL9-11 chemokines and CXCR3 receptor in teleost fish species

The coordinated migration of immune cells from lymphoid organs to in or out of the bloodstream, and towards the site of infection or tissue damage is fundamental for an efficient innate and adaptive immune response. Interestingly, an essential part of this movement is mediated by chemoattractant cytokines called chemokines. Although the nature and function of chemokines and their receptors are well documented in mammals, much research is needed to accomplish a similar level of understanding of the role of chemokines in fish immunity. The first chemokine gene identified in teleosts (rainbow trout, Oncorhynchus mykiss) was CK1 in 1998. Since then, the identification of fish chemokine orthologue genes and characterization of their role has been more complex than expected, primarily because of the whole genome duplication processes occurring in fish, and because chemokines evolve faster than other immune genes. Some of the most studied chemokines are CXCL9, CXCL10, CXCL11, and the CXCR3 receptor, all involved in T cell migration and in the induction of the T helper 1 (Th1) immune response. Data from the zebrafish and rainbow trout CXCL9-11/CXCR3 axis suggest that these chemokines and the receptor arose early in evolution and must be present in most teleost fish. However, the pieces of knowledge also indicate that different numbers of gene copies can be present in different species, with distinct regulatory expression mechanisms and probably, also with different roles, as the differential expression in fish tissues suggest. Here, we revised the current knowledge of the CXCL9-11/CXCR3 axis in teleost fishes, identifying the gaps in knowledge, and raising some hypotheses for the role of CXCL9, CXCL10 CXCL11, and CXCR3 receptor axis in fish, which can encourage further studies in the field.

NK and cells with NK-like activities in cancer immunotherapy-clinical perspectives

Natural killer (NK) cells are lymphoid cells of innate immunity that take important roles in immune surveillance. NK cells are considered as a bridge between innate and adaptive immunity, and their infiltration into tumor area is related positively with prolonged patient survival. They are defined as CD16⁺ CD56⁺ CD3^- cells in clinic. NK cells promote cytolytic effects on target cells and induce their apoptosis. Loss of NK cell cytotoxic activity and reduction in the number of activating receptors are the current issues for application of such cells in cellular immunotherapy, which resulted in the diminished long-term effects. The focus of this review is to discuss about the activity of NK cells and cells with NK-like activity including natural killer T (NKT), cytokine-induced killer (CIK) and lymphokine-activated killer (LAK) cells in immunotherapy of human solid cancers.

Chemokines and chemokine receptors in colorectal cancer; multifarious roles and clinical impact

Colorectal cancer (CRC) is considered the second cause of cancer death worldwide. The early diagnosis plays a key role in patient prognosis and subsequently overall survival. Similar to several types of cancer, colorectal cancer is also characterised by drug resistance and heterogeneity that contribute to its complexity -especially at advanced stages. However, despite the extensive research related to the identification of biomarkers associated to early diagnosis, accurate prognosis and the management of CRC patients, little progress has been made thus far. Therefore, the mortality rates, especially at advanced stages, remain high. A large family of chemoattractant cytokines called chemokines are known for their significant role in inflammation and immunity. Chemokines released by the different tumorous cells play a key role in increasing the complexity of the tumour's microenvironment. The current review investigates the role of chemokines and chemokine receptors in colorectal cancer and their potential as clinical molecular signatures that could be effectively used as a personalised therapeutic approach. We discussed how chemokine and chemokine receptors regulate the microenvironment and lead to heterogeneity in CRC. An important aspect of chemokines is their role in drug resistance which has been extensively discussed. This review also provides an overview of the current advances in the search for chemokines and chemokine receptors in CRC.

Longitudinal plasma proteomic profiling of patients with non-small cell lung cancer undergoing immune checkpoint blockade

Background

Immune checkpoint inhibitors (ICIs) have revolutionized the cancer therapy landscape due to long-term benefits in patients with advanced metastatic disease. However, robust predictive biomarkers for response are still lacking and treatment resistance is not fully understood.

Methods

We profiled approximately 800 pre-treatment and on-treatment plasma proteins from 143 ICI-treated patients with non-small cell lung cancer (NSCLC) using ELISA-based arrays. Different clinical parameters were collected from the patients including specific mutations, smoking habits, and body mass index, among others. Machine learning algorithms were used to identify a predictive signature for response. Bioinformatics tools were used for the identification of patient subtypes and analysis of differentially expressed proteins and pathways in each response group.

Results

We identified a predictive signature for response to treatment comprizing two proteins (CXCL8 and CXCL10) and two clinical parameters (age and sex). Bioinformatic analysis of the proteomic profiles identified three distinct patient clusters that correlated with multiple parameters such as response, sex and TNM (tumors, nodes, and metastasis) staging. Patients who did not benefit from ICI therapy exhibited significantly higher plasma levels of several proteins on-treatment, and enrichment in neutrophil-related proteins.

Conclusions

Our study reveals potential biomarkers in blood plasma for predicting response to ICI therapy in patients with NSCLC and sheds light on mechanisms underlying therapy resistance.

Recombinant BCG-LTAK63 Vaccine Candidate for Tuberculosis Induces an Inflammatory Profile in Human Macrophages

Tuberculosis (TB) is one of the top 10 leading causes of death worldwide. The recombinant BCG strain expressing the genetically detoxified A subunit of the thermolabile toxin from Escherichia coli (LTAK63) adjuvant (rBCG-LTAK63) has previously been shown to confer superior protection and immunogenicity compared to BCG in a murine TB infection model. To further investigate the immunological mechanisms induced by rBCG-LTAK63, we evaluated the immune responses induced by rBCG-LTAK63, BCG, and Mycobacterium tuberculosis (Mtb) H37Rv strains in experimental infections of primary human M1 and M2 macrophages at the transcriptomic and cytokine secretion levels. The rBCG-LTAK63-infected M1 macrophages more profoundly upregulated interferon-inducible genes such as IFIT3, OAS3, and antimicrobial gene CXCL9 compared to BCG, and induced higher levels of inflammatory cytokines such as IL-12(p70), TNF-β, and IL-15. The rBCG-LTAK63-infected M2 macrophages more extensively upregulated transcripts of inflammation-related genes, TAP1, GBP1, SLAMF7, TNIP1, and IL6, and induced higher levels of cytokines related to inflammation and tissue repair, MCP-3 and EGF, as compared to BCG. Thus, our data revealed an important signature of immune responses induced in human macrophages by rBCG-LTAK63 associated with increased inflammation, activation, and tissue repair, which may be correlated with a protective immune response against TB.

Pathophysiological role of chemokine MIG/CXCL9 in the development of atopic dermatitis

Background. In patients with atopic dermatitis, the persistence of microorganisms on the skin surface is high, which can enhance the expression of MIG/CXCL9, exacerbating inflammation and activating keratinocyte apoptosis, however, the dynamics of this chemokine in atopic dermatitis has not been studied.

The aim. To study the concentration of chemokine MIG/CXCL9 in the dynamics of atopic dermatitis and determine its role in the pathogenesis of dermatosis.

Materials and methods. The study included 80 patients aged 13 to 44 years with limited and widespread atopic dermatitis and 30 practically healthy volunteers. The therapy of patients, the collection of biological material were carried out in the Regional Dermatovenerologic Dispensary in Chita, laboratory tests were performed at the Chita State Medical Academy in the period from 2018 to 2020. The MIG/CXCL9 level was studied during exacerbation and remission of the disease in blood serum and skin exudate by flow cytofluorimetry using LEGENDplex Human Proinflammatory Chemokine Panel (BioLegend, USA). In healthy volunteers, skin exudate sampling was carried out by the “skin window” method. For statistical processing, the Microsoft Excel (Microsoft Corp., USA) application software package and SPSS Statistics 25.0 (IBM Corp., USA) were used.

Results. The concentration of chemokine MIG/CXCL9 in the skin exudate is greater than in the blood serum. With a limited form of the disease in adolescents, the level of MIG/CXCL9 in the skin exudate is 8.1 times higher than the control values, with a common form – 9.3 times. In adults with advanced atopic dermatitis, the concentration of chemokine IL/CXCL9 in the skin exudate is 20.8 times higher than the values of the control group.

Conclusion. In atopic dermatitis, the level of chemokine MIG/CXCL9 is higher in the cutaneous pathological process. In the pathogenesis of the disease, MIG/ CXCL9 inhibits collagen synthesis and promotes apoptosis of keratinocytes, followed by the formation of hyperreactivity of the skin, its dryness and peeling.

Altered Vaginal Microbiota Composition Correlates With Human Papillomavirus and Mucosal Immune Responses in Women With Symptomatic Cervical Ectopy

Cervical ectopy is a benign condition of the lower genital tract that is frequently detected in women of reproductive age. Although cervical ectopy is regarded as a physiological condition, some women experience symptoms such as leucorrhoea, persistent bleeding and recurrent vaginal infections that require medical intervention. Cervical ectopy has not been linked to cervical cancer, but it is thought to facilitate the acquisition of sexually transmitted diseases (STDs), like Human Papillomavirus (HPV) infection, as it provides a favorable microenvironment for virus infection and dissemination. We and others have described the presence of oncogenic HPV types in women with symptomatic cervical ectopy. The relevance of this finding and the impact of symptomatic cervical ectopy on the cervicovaginal microenvironment (vaginal microbiota, immune and inflammatory responses) are currently unknown. To shed some light into the interplay between HPV, the vaginal microbiota and mucosal immune and inflammatory responses in the context of this condition, we enrolled 156 women with symptomatic cervical ectopy and determined the presence of HPV using a type-specific multiplex genotyping assay. Overall, HPV was detected in 54.48% women, oncogenic HPV types were found in more than 90% of HPV-positive cases. The most prevalent HPV types were HPV16 (29.4%), HPV31 (21.17%) and HPV18 (15.29%). Next, we evaluated the vaginal microbial composition and diversity by 16S rDNA sequencing, and quantified levels of cytokines and chemokines by flow cytometry using bead-based multiplex assays in a sub-cohort of 63 women. IL-21 and CXCL9 were significantly upregulated in HPV-positive women (p=0.0002 and p=0.013, respectively). Women with symptomatic cervical ectopy and HPV infection had increased diversity (p<0.001), and their vaginal microbiota was enriched in bacterial vaginosis-associated anaerobes (Sneathia, Shuttleworthia, Prevotella, and Atopobium) and depleted in Lactobacillus spp. Furthermore, the vaginal microbiota of women with symptomatic cervical ectopy and HPV infection correlated with vaginal inflammation (IL-1β, rho=0.56, p=0.0004) and increased mucosal homeostatic response (IL-22, rho=0.60, p=0.0001). Taken together, our results suggest that HPV infection and dysbiotic vaginal communities could favor a vaginal microenvironment that might delay the recovery of the cervical epithelium in women with symptomatic cervical ectopy and favor STDs acquisition.

Differences in Circulating Extracellular Vesicle and Soluble Cytokines in Older Versus Younger Breast Cancer Patients With Distinct Symptom Profiles

Because extracellular vesicle (EV)-associated cytokines, both encapsulated and surface bound, have been associated with symptom severity, and may vary over the lifespan, they may be potential biomarkers to uncover underlying mechanisms of various conditions. This study evaluated the associations of soluble and EV-associated cytokine concentrations with distinct symptom profiles reported by 290 women with breast cancer prior to surgery. Patients were classified into older (≥60 years, n = 93) and younger (< 60 years, n = 197) cohorts within two previously identified distinct symptom severity profiles, that included pain, depressive symptoms, sleep disturbance, and fatigue (i.e., High Fatigue Low Pain and All Low). EVs were extracted using ExoQuick. Cytokine concentrations were determined using Luminex multiplex assay. Mann Whitney U test evaluated the differences in EV and soluble cytokine levels between symptom classes and between and within the older and younger cohorts adjusting for Karnofsky Performance Status (KPS) score, body mass index (BMI), and stage of disease. Partial correlation analyses were run between symptom severity scores and cytokine concentrations. Results of this study suggest that levels of cytokine concentrations differ between EV and soluble fractions. Several EV and soluble pro-inflammatory cytokines had positive associations with depressive symptoms and fatigue within both age cohorts and symptom profiles. In addition, in the older cohort with High Fatigue Low Pain symptom profile, EV GM-CSF concentrations were higher compared to the All Low symptom profile (p < 0.05). Albeit limited by a small sample size, these exploratory analyses provide new information on the association between cytokines and symptom profiles of older and younger cohorts. Of note, unique EV-associated cytokines were found in older patients and in specific symptom classes. These results suggest that EVs may be potential biomarker discovery tools. Understanding the mechanisms that underlie distinct symptom class profiles categorized by age may inform intervention trials and offer precision medicine approaches.

Extracellular vesicles from triple negative breast cancer promote pro-inflammatory macrophages associated with better clinical outcome

Our work uncovers mechanisms by which tumor cells impact tumor-associated macrophages in human triple-negative breast cancer. Via extracellular vesicles (EVs), these tumors promote macrophages with proinflammatory features, correlated with better clinical outcome. Our results suggest exploration of these EVs as tools, alone, or in combination with other therapies, to promote a favorable environment for the generation of anti-tumor immune responses.

Tumor associated macrophages (TAMs), which differentiate from circulating monocytes, are pervasive across human cancers and comprise heterogeneous populations. The contribution of tumor-derived signals to TAM heterogeneity is not well understood. In particular, tumors release both soluble factors and extracellular vesicles (EVs), whose respective impact on TAM precursors may be different. Here, we show that triple negative breast cancer cells (TNBCs) release EVs and soluble molecules promoting monocyte differentiation toward distinct macrophage fates. EVs specifically promoted proinflammatory macrophages bearing an interferon response signature. The combination in TNBC EVs of surface CSF-1 promoting survival and cargoes promoting cGAS/STING or other activation pathways led to differentiation of this particular macrophage subset. Notably, macrophages expressing the EV-induced signature were found among patients’ TAMs. Furthermore, higher expression of this signature was associated with T cell infiltration and extended patient survival. Together, this data indicates that TNBC-released CSF-1-bearing EVs promote a tumor immune microenvironment associated with a better prognosis in TNBC patients.

CXCL9-modified CAR T cells improve immune cell infiltration and antitumor efficacy

Chimeric antigen receptor (CAR) T cells remain unsatisfactory in treating solid tumors. The frequency of tumor-infiltrating T cells is closely related to the good prognosis of patients. Augmenting T cell accumulation in the tumor microenvironment is essential for tumor clearance. To overcome insufficient immune cell infiltration, innovative CAR designs need to be developed immediately. CXCL9 plays a pivotal role in regulating T cell migration and inhibiting tumor angiogenesis. Therefore, we engineered CAR T cells expressing CXCL9 (CART-CXCL9). The addition of CXCL9 enhanced cytokine secretion and cytotoxicity of CAR T cells and endowed CAR T cells with the ability to recruit activated T cells and antiangiogenic effect. In tumor-bearing mice, CART-CXCL9 cells attracted more T cell trafficking to the tumor site and inhibited angiogenesis than conventional CAR T cells. Additionally, CART-CXCL9 cell therapy slowed tumor growth and prolonged mouse survival, displaying superior antitumor activity. Briefly, modifying CAR T cells to express CXCL9 could effectively improve CAR T cell efficacy against solid tumors.

Identifying Potential Biomarkers of Prognostic Value in Colorectal Cancer via Tumor Microenvironment Data Mining

Colorectal cancer (CRC) is a common cancer that has increased rapidly worldwide in the past decades with a relatively high mortality rate. An increasing body of evidence has highlighted the importance of infiltrating immune and stromal cells in CRC. In this study, based on gene expression data of CRC patients in TCGA database we evaluated immune and stromal scores in tumor microenvironment using ESTIMATE method. Results showed there was potential correlation between these scores and the prognosis, and that patients with higher immune score and lower stromal score had longer survival time. We found that immune score was correlated with clinical characteristics including tumor location, tumor stage, and survival time. Specifically, the right-sided colon cancer had markedly elevated immune score, compared to left-sided colon cancer and rectal cancer. These results might be useful for understanding tumor microenvironment in colorectal cancer. Through the differential analysis we got a list of genes significantly associated with immune and stromal scores. Gene Set Enrichment and protein-protein interaction network analysis were used to further illustrate these differentially expressed genes. Finally, 15 hub genes were identified, and three (CXCL9, CXCL10 and SELL) of them were validated with favorable outcomes in CRC patients. Our result suggested that these tumor microenvironment related genes might be potential biomarkers for the prognosis of CRC.

Immune Score-based Molecular Subtypes and Signature Associated with Clinical Outcome in Hepatoblastoma

Background: This study aimed to identify genes related to the immune score of hepatoblastoma, examine the characteristics of the immune microenvironment of hepatoblastoma, and construct a risk scoring system for predicting the prognosis of hepatoblastoma. Methods: Through using the gene chip data of patients with hepatoblastoma with survival data in the ArrayExpress and GEO databases, the immune score of hepatoblastoma was calculated by the ESITIMATE algorithm, and the prognostic value of immune score in patients with hepatoblastoma was studied by the survival analysis. Genes related to the immune score were identified by the WGCNA algorithm. According to these genes, patients with hepatoblastoma were clustered unsupervised. Finally, the risk scoring system was constructed according to the immune score-related genes. Results: The immune score calculated by the ESTIMATE algorithm had a good prognostic value in patients with hepatoblastoma. Patients with high immune scores had better OS than those with low immune scores (P < 0.001). A total of 146 immune score-related genes were identified by WGCNA analysis, and univariate COX regression analysis indicated that 59 of the genes had prognostic value. According to the unsupervised clustering results of the 146 immune score-related genes, patients with hepatoblastoma could be divided into two subtypes with different prognoses, namely molecular subtype 1 and subtype 2, with molecular subtype 1 having a better prognosis. The immunocyte infiltration analysis results showed that the difference between the two subtypes was mainly in activated CD4 T cells, activated dendritic cells, CD56 bright natural killer cells, the macrophage, and regulatory T cells. According to the immune score-related genes, a risk scoring system was constructed based on a five-gene signature. After the cut-off value was determined, patients with hepatoblastoma were divided into a high-risk group and a low-risk group. The prognosis of the two groups was different. Conclusions: The immune score has a good prognostic value in patients with hepatoblastoma. Based on the different expression patterns of immune score-related genes, hepatoblastoma can be divided into two different prognostic molecular subtypes, showing different immunocyte infiltration patterns. The established risk scoring system based on a five-gene signature has a good predictive value in patients with hepatoblastoma.

Moving towards the Future of Radio-Immunotherapy: Could We “Tailor” the Abscopal Effect on Head and Neck Cancer Patients?

The abscopal effect (AbE) is defined as radiation-induced shrinkage of distant, non-treated, neoplastic lesions and it is considered the best clinical picture of the efficient immune stimulation by irradiation. The first report about abscopal tumor regression upon radiotherapy dates back to the beginning of the 20th century. The growing preclinical and clinical synergism between radiation and immunotherapy gave birth the purpose to more easily reproduce the abscopal effect, nevertheless, it is still rare in clinical practice. In this review we summarize immunological modulation of radiotherapy, focusing on the well-balanced equilibrium of tumor microenvironment and how radio-immunotherapy combinations can perturb it, with particular attention on head and neck squamous cell cancer. Finally, we investigate future perspectives, with the aim to “tailor” the abscopal effect to the patient.

Oxidized carbon black nanoparticles induce endothelial damage through C-X-C chemokine receptor 3-mediated pathway

Oxidation of engineered nanomaterials during application in various industrial sectors can alter their toxicity. Oxidized nanomaterials also have widespread industrial and biomedical applications. In this study, we evaluated the cardiopulmonary hazard posed by these nanomaterials using oxidized carbon black (CB) nanoparticles (CBox) as a model particle. Particle surface chemistry was characterized by X-ray photo electron spectroscopy (XPS) and Fourier-transform infrared spectroscopy (FTIR). Colloidal characterization and in vitro dosimetry modeling (particle kinetics, fate and transport modeling) were performed. Lung inflammation was assessed following oropharyngeal aspiration of CB or oxidized CBox particles (20 μg per mouse) in C57BL/6J mice. Toxicity and functional assays were also performed on murine macrophage (RAW 264.7) and endothelial cell lines (C166) with and without pharmacological inhibitors. Oxidant generation was assessed by electron paramagnetic resonance spectroscopy (EPR) and via flow cytometry. Endothelial toxicity was evaluated by quantifying pro-inflammatory mRNA expression, monolayer permeability, and wound closure. XPS and FTIR spectra indicated surface modifications, the appearance of new functionalities, and greater oxidative potential (both acellular and in vitro) of CBox particles. Treatment with CBox demonstrated greater in vivo inflammatory potentials (lavage neutrophil counts, secreted cytokine, and lung tissue mRNA expression) and air-blood barrier disruption (lavage proteins). Oxidant-dependent pro-inflammatory signaling in macrophages led to the production of CXCR3 ligands (CXCL9,10,11). Conditioned medium from CBox-treated macrophages induced significant elevation in endothelial cell pro-inflammatory mRNA expression, enhanced monolayer permeability and impairment of scratch healing in CXCR3 dependent manner. In summary, this study mechanistically demonstrated an increased biological potency of CBox particles and established the role of macrophage-released chemical mediators in endothelial damage.

CSF1R Inhibition Combined with GM-CSF Reprograms Macrophages and Disrupts Protumoral Interplays with AML Cells

Relapse is a major issue in acute myeloid leukemia (AML) and while the contribution of gene mutations in developing drug resistance is well established, little is known on the role of macrophages (MΦs) in an AML cell microenvironment. We examined whether myeloblasts could educate MΦs to adopt a protumoral orientation supporting myeloblast survival and resistance to therapy. Flow cytometry analyses demonstrated that M2-like CD163⁺ MΦs are abundantly present, at diagnosis, in the bone marrow of AML patients. We showed that myeloblasts, or their conditioned medium, polarize monocytes to M2-like CD163⁺ MΦs, induce the secretion of many protumoral factors, and promote myeloblast survival and proliferation as long as close intercellular contacts are maintained. Importantly, pharmacologic inhibition of the CSF1 receptor (CSF1R), in the presence of GM-CSF, reprogrammed MΦ polarization to an M1-like orientation, induced the secretion of soluble factors with antitumoral activities, reduced protumoral agonists, and promoted the apoptosis of myeloblasts interacting with MΦs. Furthermore, myeloblasts, which became resistant to venetoclax or midostaurin during their interplay with protumoral CD163⁺ MΦs, regained sensitivity to these targeted therapies following CSF1R inhibition in the presence of GM-CSF. These data reveal a crucial role of CD163⁺ MΦ interactions with myeloblasts that promote myeloblast survival and identify CSF1R inhibition as a novel target for AML therapy.

Irradiation Mediates IFNα and CXCL9 Expression in Non-Small Cell Lung Cancer to Stimulate CD8+ T Cells Activity and Migration toward Tumors

Irradiation-broken DNA fragments increase type I interferon and chemokines secretion in tumor cells. Since radiotherapy may augment tumor immunotherapy, we hypothesize that the chemokines increased by irradiation could recruit CD8⁺ T cells to suppress tumor proliferation. This study intended to unveil the secreted factors activating and recruiting CD8⁺ T cells in non-small-cell lung cancer (NSCLC). EGFR-positive A549 was selected and treated by X-irradiation (IR) to identify the overexpression of chemokines associated to CD8⁺ T cell cytotoxicity and recruitment. A transwell assay with Alexa 488-labeled CD8⁺ T cells was used to evaluate CD8⁺ T cell motility in vitro. A nuclear imaging platform by In¹¹¹-labeled nivolumab was used to track CD8⁺ T cells homing to tumors in vivo. The activation markers GZMB, PRF-1, and IFNγ, migration marker CD183 (CXCR3), and inhibitory marker CD274 (PD-1), were measured and compared in CD8⁺ T cells with A549 co-cultured, chemokines treated, and patients with late-stage lung cancer. We found that IR not only suppressed A549 proliferation but also induced IFNα and CXCL9 expression (p < 0.05). IFNα majorly increased IFNγ levels in CD8⁺ T cells (p < 0.05) and synergistically with CXCL9 enhanced CD8⁺ T cell migration in vitro (p < 0.05). We found that CXCR3 and PD-1 were down-regulated and up-regulated, respectively, in the peripheral blood CD8⁺ T cells in patients with lung cancer (n = 4 vs. healthy n = 3, both p < 0.05), which exhibited reduction of cell motility (p < 0.05). The in vivo nuclear imaging data indicated highly CD8⁺ T cells migrated to A549-induced tumors. In addition, we demonstrated that healthy PBMCs significantly suppressed the parallel tumor growth (p < 0.05) and the radioresistant tumor growth in the tumor xenograft mice (p < 0.05), but PBMCs from patients with lung cancer had lost the anti-tumor capacity. We demonstrated that IR induced IFNα and CXCL9 expression in A549 cells, leading to CD8⁺ T cell migration. This study unveiled a potential mechanism for radiotherapy to activate and recruit CD8⁺ T cells to suppress lung tumors.

The Right Partner in Crime: Unlocking the Potential of the Anti-EGFR Antibody Cetuximab via Combination With Natural Killer Cell Chartering Immunotherapeutic Strategies

Cetuximab has an established role in the treatment of patients with recurrent/metastatic colorectal cancer and head and neck squamous cell cancer (HNSCC). However, the long-term effectiveness of cetuximab has been limited by the development of acquired resistance, leading to tumor relapse. By contrast, immunotherapies can elicit long-term tumor regression, but the overall response rates are much more limited. In addition to epidermal growth factor (EGFR) inhibition, cetuximab can activate natural killer (NK) cells to induce antibody-dependent cellular cytotoxicity (ADCC). In view of the above, there is an unmet need for the majority of patients that are treated with both monotherapy cetuximab and immunotherapy. Accumulated evidence from (pre-)clinical studies suggests that targeted therapies can have synergistic antitumor effects through combination with immunotherapy. However, further optimizations, aimed towards illuminating the multifaceted interplay, are required to avoid toxicity and to achieve better therapeutic effectiveness. The current review summarizes existing (pre-)clinical evidence to provide a rationale supporting the use of combined cetuximab and immunotherapy approaches in patients with different types of cancer.

Cytotoxic T-Cell Trafficking Chemokine Profiles Correlate With Defined Mucosal Microbial Communities in Colorectal Cancer

The involvement of gut microbiota in T-cell trafficking into tumor tissue of colorectal cancer (CRC) remains to be further elucidated. The current study aimed to evaluate the expression of major cytotoxic T-cell trafficking chemokines (CTTCs) and chemokine-associated microbiota profiles in both tumor and adjacent normal tissues during CRC progression. We analyzed the expression of chemokine C-X-C motif ligands 9, 10, and 11 (CXCL9, CXCL10, and CXCL11), and C-C motif ligand 5 (CCL5), characterized gut mucosa-associated microbiota (MAM), and investigated their correlations in CRC patients. Our results showed that the expression of CXCL9, CXCL10, and CXCL11 was significantly higher in tumor than in adjacent normal tissues in 136 CRC patients. Notably, the high expression of CXCL9 in tumor tissues was associated with enhanced CD8⁺ T-cell infiltration and improved survival. Moreover, the MAM in tumor tissues showed reduction of microbial diversity and increase of oral bacteria. Microbial network analysis identified differences in microbial composition and structure between tumor and adjacent normal tissues. In addition, stronger associations between oral bacteria and other gut microbes were observed. Furthermore, the correlation analysis between the defined MAM and individual CTTCs showed that the CTTCs' correlated operational taxonomic units (OTUs) in tumor and adjacent normal tissues rarely overlap with each other. Notably, all the enriched OTUs were positively correlated with the CTTCs in either tumor or adjacent normal tissues. Our findings demonstrated stronger interactions between oral bacteria and gut microbes, and a shifted correlation pattern between MAM and major CTTCs in tumor tissues, underlining possible mechanisms of gut microbiota-host interaction in CRC.

Pan-cancer analysis of longitudinal metastatic tumors reveals genomic alterations and immune landscape dynamics associated with pembrolizumab sensitivity

Serial circulating tumor DNA (ctDNA) monitoring is emerging as a non-invasive strategy to predict and monitor immune checkpoint blockade (ICB) therapeutic efficacy across cancer types. Yet, limited data exist to show the relationship between ctDNA dynamics and tumor genome and immune microenvironment in patients receiving ICB. Here, we present an in-depth analysis of clinical, whole-exome, transcriptome, and ctDNA profiles of 73 patients with advanced solid tumors, across 30 cancer types, from a phase II basket clinical trial of pembrolizumab (NCT02644369) and report changes in genomic and immune landscapes (primary outcomes). Patients stratified by ctDNA and tumor burden dynamics correspond with survival and clinical benefit. High mutation burden, high expression of immune signatures, and mutations in BRCA2 are associated with pembrolizumab molecular sensitivity, while abundant copy-number alterations and B2M loss-of-heterozygosity corresponded with resistance. Upon treatment, induction of genes expressed by T cell, B cell, and myeloid cell populations are consistent with sensitivity and resistance. We identified the upregulated expression of PLA2G2D, an immune-regulating phospholipase, as a potential biomarker of adaptive resistance to ICB. Together, these findings provide insights into the diversity of immunogenomic mechanisms that underpin pembrolizumab outcomes.

Interpretable systems biomarkers predict response to immune-checkpoint inhibitors

Cancer cells can leverage several cell-intrinsic and -extrinsic mechanisms to escape immune system recognition. The inherent complexity of the tumor microenvironment, with its multicellular and dynamic nature, poses great challenges for the extraction of biomarkers of immune response and immunotherapy efficacy. Here, we use RNA-sequencing (RNA-seq) data combined with different sources of prior knowledge to derive system-based signatures of the tumor microenvironment, quantifying immune-cell composition and intra- and intercellular communications. We applied multi-task learning to these signatures to predict different hallmarks of immune responses and derive cancer-type-specific models based on interpretable systems biomarkers. By applying our models to independent RNA-seq data from cancer patients treated with PD-1/PD-L1 inhibitors, we demonstrated that our method to Estimate Systems Immune Response (EaSIeR) accurately predicts therapeutic outcome. We anticipate that EaSIeR will be a valuable tool to provide a holistic description of immune responses in complex and dynamic systems such as tumors using available RNA-seq data.

Toxoplasma gondii could have a possible role in the cancer mechanism by modulating the host's cell response

Toxoplasma gondii, which manipulates many signaling pathways to achieve persistence in host cells, is intimately linked to immune and inflammation responses. However, there is still lack of information about the impact of T. gondii on cellular and immune responses. This study was designed to seek the impact of T. gondii infection causing life-long inflammation in brain, on cancer mechanism. To identify molecular effects of the T. gondii and understand the association between the functional perturbations occurring during infection and cancer development, the transcriptomic datasets obtained mice infected with T. gondii were downloaded from GEO. The differentially expressed genes (DEGs) were identified and functional enrichment analysis was performed using IPA platform, then all results were evaluated with comparison analyses. Subsequently, a T. gondii infection model with human neuroepithelioma cell culture was performed in order to validate top DEGs participated in common networks/pathways in cancer mechanism. Transcriptomic analyses of infected mice and in vitro cell culture model revealed a strong immune response and inflammation occurred by parasite-induced damage and parasite-associated immunopathology in host cell and tissue. T. gondii infection could modulate certain signaling pathways of host, which were also common to those perturbed in carcinogenesis. Interestingly, the network analysis of the data sets predicted an activation in development of solid cancer vice versa inhibition in hematological cancer during T. gondii infection. Parasite might also control the tumor growth due to its potent immune-stimulant effects. As result, T. gondii infection generating a continual inflammation in tissues might potentially contribute to cancer development by regulating critical host signaling pathways or reveal an anti-tumoral activity.

Analysis of Immune-Related Signatures Related to CD4+ T Cell Infiltration With Gene Co-Expression Network in Pancreatic Adenocarcinoma

Background

Pancreatic ductal adenocarcinoma (PDAC) is one of the most invasive solid malignancies. Immunotherapy and targeted therapy confirmed an existing certain curative effect in treating PDAC. The aim of this study was to develop an immune-related molecular marker to enhance the ability to predict Stages III and IV PDAC patients.

Method

In this study, weighted gene co-expression network (WGCNA) analysis and a deconvolution algorithm (CIBERSORT) that evaluated the cellular constituent of immune cells were used to evaluate PDAC expression data from the GEO (Gene Expression Omnibus) datasets, and identify modules related to CD4+ T cells. LASSO Cox regression analysis and Kaplan-Meier curve were applied to select and build prognostic multi-gene signature in TCGA Stages III and IV PDAC patients (N = 126). This was followed by independent Stages III and IV validation of the gene signature in the International Cancer Genome Consortium (ICGC, N = 62) and the Fudan University Shanghai Cancer Center (FUSCC, N = 42) cohort. Inherited germline mutations and tumor immunity exploration were applied to elucidate the molecular mechanisms in PDAC. Univariate and Multivariate Cox regression analyses were applied to verify the independent prognostic factors. Finally, a prognostic nomogram was created according to the TCGA-PDAC dataset.

Results

A four-gene signature comprising NAPSB, ZNF831, CXCL9 and PYHIN1 was established to predict overall survival of PDAC. This signature also robustly predicted survival in two independent validation cohorts. The four-gene signature could divide patients into high and low-risk groups with disparity overall survival verified by a Log-rank test. Expression of four genes positively correlated with immunosuppression activity (PD-L1 and PD1). Immune-related genes nomogram and corresponding calibration curves showed significant performance for predicting 3-year survival in TCGA-PDAC dataset.

Conclusion

We constructed a novel four-gene signature to predict the prognosis of Stages III and IV PDAC patients by applying WGCNA and CIBERSORT algorithm scoring to transcriptome data different from traditional methods of filtrating for differential genes in cancer and healthy tissues. The findings may provide reference to predict survival and was beneficial to individualized management for advanced PDAC patients.

Infiltrative tumour growth pattern correlates with poor outcome in oesophageal cancer

Objective

Oesophageal cancer (OEC) is an aggressive disease with a poor survival rate. Prognostic markers are thus urgently needed. Due to the demonstrated prognostic value of histopathological growth pattern (HGP) in other cancers, we performed a retrospective assessment of HGP in patients suffering from invasive OEC.

Design

A first cohort composed of 89 treatment-naïve operated patients with OEC from The Cancer Genome Atlas (TCGA) public database was constituted, from which H&E images and RNA-sequencing data were retrieved. Next, a second cohort composed of 99 patients with OEC treated and operated in a Belgian hospital was established. H&E-stained sections and extracted tumorous RNA were obtained from the samples. HGP were assessed on H&E slides as infiltrative (IGP) or expansive (EGP). TCGA RNA-sequencing data were analysed through the gene set enrichment analysis and Cytoscape softwares. Real-time quantitative PCR (qPCR) experiments were performed to assess gene expression in the Belgian cohort.

Results

IGP patients displayed a grim prognosis compared with EGP patients, while IGP was found as associated with numerous lymphovascular emboli and perinervous infiltrations. Analyses of the TCGA expression data showed that angiogenesis, epithelial-to-mesenchymal transition (EMT) and inflammation were significantly upregulated in IGP compared with EGP samples. qPCR experiments of three genes appearing as highly upregulated in each pathway showed no difference in expression according to the HGP.

Conclusion

The current study demonstrates the poor prognostic value carried by IGP in OC and suggests angiogenesis, EMT and inflammation as key carcinogenetic pathways upregulated in this pattern.

Construction of a five-gene prognostic model based on immune-related genes for the prediction of survival in pancreatic cancer

Purpose: To identify differentially expressed immune-related genes (DEIRGs) and construct a model with survival-related DEIRGs for evaluating the prognosis of patients with pancreatic cancer (PC).

Methods: Six microarray gene expression datasets of PC from the Gene Expression Omnibus (GEO) and Immunology Database and Analysis Portal (ImmPort) were used to identify DEIRGs. RNA sequencing and clinical data from The Cancer Genome Atlas Program-Pancreatic Adenocarcinoma (TCGA-PAAD) database were used to establish the prognostic model. Univariate, least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression analyses were applied to determine the final variables of the prognostic model. The median risk score was used as the cut-off value to classify samples into low- and high-risk groups. The prognostic model was further validated using an internal validation set of TCGA and an external validation set of GSE62452.

Results: In total, 142 DEIRGs were identified from six GEO datasets, 47 were survival-related DEIRGs. A prognostic model comprising five genes (i.e., ERAP2, CXCL9, AREG, DKK1, and IL20RB) was established. High-risk patients had poor survival compared with low-risk patients. The 1-, 2-, 3-year area under the receiver operating characteristic (ROC) curve of the model reached 0.85, 0.87, and 0.93, respectively. Additionally, the prognostic model reflected the infiltration of neutrophils and dendritic cells. The expression of most characteristic immune checkpoints was significantly higher in the high-risk group versus the low-risk group.

Conclusions: The five-gene prognostic model showed reliably predictive accuracy. This model may provide useful information for immunotherapy and facilitate personalized monitoring for patients with PC.

A comparison expression analysis of CXCR4, CXCL9 and Caspase-9 in dermal vascular endothelial cells between keloids and normal skin on chemotaxis and apoptosis

This present study was designed to explore key biological characteristics and biomarkers associated with dermal vascular endothelial cells of keloids. GSE121618 dataset was downloaded in the Gene Expression Omnibus (GEO) Database, including the KECs group and NVECs group. Through GEO2R, we have screened the differentially expressed genes (DEGs) and performed gene ontology (GO), Gene Set Enrichment Analysis (GSEA), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Then, we constructed a protein-protein interaction (PPI) network and analyzed hub genes via the Search Tool for the Retrieval of Interacting Genes (STRING) Online Database and Cytoscape software. Furthermore, experiments were performed to validate the expression of selected genes, including H&E staining, immunohistochemical staining, Western blot, and RT-qPCR. A total of 1040 DEGs were selected with GEO2R online tools. Most of the enriched pathways and processes focus on cell migration, tube development, chemotaxis, cell motility, and regulation of apoptosis. With the assistance of STRING and Cytoscape, hub genes were selected. In our validation experiments of RT-qPCR, the mRNA expression of selected genes has significant differences between different groups in tissue and cell experiments. As was shown in immunohistochemical staining, the proteins of CXCR4, CXCL9, and Caspase-9 had higher expression levels in tissue samples of the Keloid group than the Normal skin group. Western blot and RT-qPCR in dermal vascular endothelial cell experiments were consistent with the aforementioned results. This study has provided a deeper analysis of the pathogenesis of dermal vascular endothelial cells in keloids. Genes of CXCR4, CXCL9, and Caspase-9 may influence the processes of inflammatory responses and vascular endothelial cell apoptosis to exert crucial effects in the development of keloids. Abbreviations: GEO: gene expression omnibus; DEGs: differentially expressed genes; KVECs: keloid vascular endothelial cells; NVECs: normal skin vascular endothelial cells; GO: gene ontology; KEGG: Kyoto encyclopedia of genes and genomes; PPI: protein protein interaction; BP: biological process; CC: cellular component; MF: molecular function; GSEA: gene set enrichment analysis; STRING: search tool for the retrieval of interacting genes; MCODE: molecular complex detection.

Construction of TME and Identification of crosstalk between malignant cells and macrophages by SPP1 in hepatocellular carcinoma

Liver cancer accounts for 6% of all malignancies causing death worldwide, and hepatocellular carcinoma (HCC) is the most common histological type. HCC is a heterogeneous cancer, but how the tumour microenvironment (TME) of HCC contributes to the progression of HCC remains unclear. In this study, we investigated the immune microenvironment by multiomics analysis. The tumour immune infiltration characteristics of HCC were determined at the genomic, epigenetic, bulk transcriptome and single-cell levels by data from The Cancer Genome Atlas portal and the Gene Expression Omnibus (GEO). An epigenetic immune-related scoring system (EIRS) was developed to stratify patients with poor prognosis. SPP1, one gene in the EIRS system, was identified as an immune-related predictor of poor survival in HCC patients. Through receptor-ligand pair analysis in single-cell RNA-seq, SPP1 was indicated to mediate the crosstalk between HCC cells and macrophages via SPP1-CD44 and SPP1-PTGER4 association. In vitro experiments further validate SPP1 can trigger the polarization of macrophages to M2-phenotype tumour-associated macrophages (TAMs).

Pparγ1 Facilitates ErbB2-Mammary Adenocarcinoma in Mice

HER2, which is associated with clinically aggressive disease, is overexpressed in 15-20% of breast cancers (BC). The host immune system participates in the therapeutic response of HER2⁺ breast cancer. Identifying genetic programs that participate in ErbB2-induced tumors may provide the rational basis for co-extinction therapeutic approaches. Peroxisome proliferator-activated receptor γ (PPARγ), which is expressed in a variety of malignancies, governs biological functions through transcriptional programs. Herein, genetic deletion of endogenous Pparγ1 restrained mammary tumor progression, lipogenesis, and induced local mammary tumor macrophage infiltration, without affecting other tissue hematopoietic stem cell pools. Endogenous Pparγ1 induced expression of both an EphA2-Amphiregulin and an inflammatory INFγ and Cxcl5 signaling module, that was recapitulated in human breast cancer. Pparγ1 bound directly to growth promoting and proinflammatory target genes in the context of chromatin. We conclude Pparγ1 promotes ErbB2-induced tumor growth and inflammation and represents a relevant target for therapeutic coextinction. Herein, endogenous Pparγ1 promoted ErbB2-mediated mammary tumor onset and progression. PPARγ1 increased expression of an EGF-EphA2 receptor tyrosine kinase module and a cytokine/chemokine 1 transcriptional module. The induction of a pro-tumorigenic inflammatory state by Pparγ1 may provide the rationale for complementary coextinction programs in ErbB2 tumors.

Integrated Genomic and Transcriptomic Analysis reveals key genes for predicting dual-phenotype Hepatocellular Carcinoma Prognosis

Dual-phenotype hepatocellular carcinoma (DPHCC) expresses both hepatocyte and cholangiocyte markers, and is characterized by high recurrence and low survival rates. The underlying molecular mechanisms of DPHCC pathogenesis are unclear.

We performed whole exome sequencing and RNA sequencing of three subtypes of HCC (10 DPHCC, 10 CK19-positive HCC, and 14 CK19-negative HCC), followed by integrated bioinformatics analysis, including somatic mutation analysis, mutation signal analysis, differential gene expression analysis, and pathway enrichment analysis. Cox proportional hazard regression analyses were applied for exploring survival related characteristics.

We found that mutated genes in DPHCC patients were associated with carcinogenesis and immunity, and the up-regulated genes were mainly enriched in transcription-related and cancer-related pathways, and the down-regulated genes were mainly enriched in immune-related pathways. CXCL9 was selected as the hub gene, which is associated with immune cells and survival prognosis. Our results showed that low CXCL9 expression was significantly associated with poor prognosis, and its expression was significantly reduced in DPHCC samples.

In conclusion, we explored the molecular mechanisms governing DPHCC development and progression and identified CXCL9, which influences the immune microenvironment and prognosis of DPHCC and might be new clinically significant biomarkers for predicting prognosis.