Identification of Therapeutic Targets and Prognostic Biomarkers Among CXC Chemokines in the Renal Cell Carcinoma Microenvironment

Study on correlation between CXCL13 and prognosis and immune characteristics of ovarian cancer

Ovarian cancer (OC) has a limited immunotherapeutic response; hence, this study aimed to investigate the relationship between CXC-chemokine ligand 13 (CXCL13) expression and overall survival (OS) rate, key immune pathways, degree of immune cell infiltration, and progressive disease (PD)-1 checkpoint blockade. A total of 703 differentially expressed genes were obtained from "The Cancer Genome Atlas" (TCGA) database based on the immune and stromal scores of 379 OC patients for getting the targeted gene CXCL13. The association between CXCL13 and OS in OC patients, biological function annotation of CXCL13, and its correlation with immune components were assessed. The results indicated that upregulated CXCL13 expression was positively correlated with better OC patient prognosis. CXCL13 expression was associated with 6 immune-related pathways, 10 immune cells, and PD-1 expression of OC micro-environment. Moreover, high expression of CXCL13 was related to a better tumor response and more extended tumor-stable stage after PD-1 blocking therapy in IMvigor210. The study concluded that CXCL13 could be a prognostic marker and a potential immunotherapy target for OC patients, especially PD-1 checkpoint blockade.

Predictive value of CXCL1+_FAP+ phenotype in CAFs for distant metastasis and its correlation with PD-L1 expression in locoregionally advanced nasopharyngeal carcinoma patients

OBJECTIVE

There is a lack of effective biomarkers for predicting the distant metastasis in nasopharyngeal carcinoma (NPC). We aimed to explore the expression of FAP+Cancer-associated fibroblasts (CAFs) derived CXCL1 in NPC and its predictive values for distant metastasis and correlation with PD-L1 expression.

MATERIALS AND METHODS

A total of 345 patients with locoregionally advanced NPC were retrospectively enrolled (the training cohort: the validation cohort = 160:185). Co-expression of CXCL1 and FAP and the expression of PD-L1 were detected by multi-immunofluorescence staining and immunohistochemistry, respectively. The primary end-point was distant metastasis-free survival (DMFS). The Kaplan-Meier method was used to calculate the survival. The Cox proportional hazards model was used to assess prognostic risk factors.

RESULTS

A novel CXCL1+_FAP+ phenotype in CAFs was identified in NPC and then used to divide patients into low and high risk groups. Both in the training cohort and validation cohort, patients in the high risk group had poorer DMFS, overall survival (OS), progression-free survival (PFS) and locoregional relapse-free survival (LRFS) than patients in the low risk group. Multivariate analysis revealed CXCL1+_FAP+ phenotype was an independent prognostic factor for DMFS, OS, PFS and LRFS. Further results showed patients in the high risk group had higher PD-L1 expression than those in the low risk group.

CONCLUSION

Our study showed CXCL1+_FAP+ phenotype in CAFs could effectively classified locoregionally advanced NPC patients into different risk groups for distant metastasis and might be a potential biomarker for anti-PD-1/PD-L1 immunotherapy.

CXCL3: A key player in tumor microenvironment and inflammatory diseases

CXCL3 (C-X-C Motif Chemokine 3), a member of the C-X-C chemokine subfamily, operates as a potent chemoattractant for neutrophils, thereby orchestrating the recruitment and migration of leukocytes alongside eliciting an inflammatory response. Recent inquiries have shed light on the pivotal roles of CXCL3 in the context of carcinogenesis. In the tumor microenvironment, CXCL3 emanating from both tumor and stromal cells intricately modulates cellular behaviors through autocrine and paracrine actions, primarily via interaction with its receptor CXCR2. Activation of signaling cascades such as ERK/MAPK, AKT, and JAK2/STAT3 underscores CXCL3's propensity to favor tumorigenic processes. However, CXCL3 exhibits dualistic behaviors, as evidenced by its capacity to exert anti-tumor effects under specific conditions. Additionally, the involvement of CXCL3 extends to inflammatory disorders like eclampsia, obesity, and asthma. This review encapsulates the structural attributes, biological functionalities, and molecular underpinnings of CXCL3 across both tumorigenesis and inflammatory diseases.

Immune analysis of urine and plasma samples from patients with clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is the third most common type of urological malignancy worldwide, and it is associated with a silent progression and late manifestation. Patients with a metastatic form of ccRCC have a poor prognosis; however, when the disease is diagnosed early, it is largely curable. Currently, there are no biomarkers available in clinical practice for ccRCC. Thus, the aim of the present study was to measure 27 biologically relevant cytokines in preoperative and postoperative urine samples, and in preoperative plasma samples from 34 patients with ccRCC, and to evaluate their diagnostic significance. The concentrations of cytokines were assessed by multiplex immune assay. The results showed significantly higher levels of IL-1 receptor antagonist, IL-6, IL-15, chemokine (C-C motif) ligand (CCL)2, CCL3, CCL4, C-X-C motif ligand (CXCL)10, granulocyte-macrophage colony stimulating factor (GM-CSF) and platelet-derived growth factor-BB (PDGF-BB), and lower levels of granulocyte colony stimulating factor (G-CSF) in urine samples from patients prior to surgery compared with those in the controls. Notably, the urine levels of G-CSF, IL-5 and vascular endothelial growth factor differed following tumor removal compared with the preoperative urine levels. In addition, urinary G-CSF, GM-CSF, IL-6, CXCL10, CCL5 and PDGF-BB appeared to be potential markers of tumor grade. Plasma from patients with ccRCC contained significantly higher levels of IL-6 and lower levels of CCL2 than control plasma. In conclusion, the present findings indicated that urinary and circulating cytokines may represent a promising novel tool for the early diagnosis of ccRCC and/or prediction of tumor grade.

The Clinical Significance and Involvement in Molecular Cancer Processes of Chemokine CXCL1 in Selected Tumors

Chemokines play a key role in cancer processes, with CXCL1 being a well-studied example. Due to the lack of a complete summary of CXCL1's role in cancer in the literature, in this study, we examine the significance of CXCL1 in various cancers such as bladder, glioblastoma, hemangioendothelioma, leukemias, Kaposi's sarcoma, lung, osteosarcoma, renal, and skin cancers (malignant melanoma, basal cell carcinoma, and squamous cell carcinoma), along with thyroid cancer. We focus on understanding how CXCL1 is involved in the cancer processes of these specific types of tumors. We look at how CXCL1 affects cancer cells, including their proliferation, migration, EMT, and metastasis. We also explore how CXCL1 influences other cells connected to tumors, like promoting angiogenesis, recruiting neutrophils, and affecting immune cell functions. Additionally, we discuss the clinical aspects by exploring how CXCL1 levels relate to cancer staging, lymph node metastasis, patient outcomes, chemoresistance, and radioresistance.

The immune landscape and prognostic analysis of CXCL8 immune-related genes in cervical squamous cell carcinoma

Cervical squamous cell carcinoma (CESC), one of the most common malignancies in women, imposes a significant burden on women's health worldwide. Despite extensive research, the molecular and pathogenic mechanisms of cervical squamous cell carcinoma and CESC remain unclear. This study aimed to explore the immune-related genes, immune microenvironment infiltration, and prognosis of CESC, providing a theoretical basis for guiding clinical treatment. Initially, by mining four gene sets and immune-related gene sets from public databases, 14 immune-related genes associated with CESC were identified. Through univariate and multivariate COX regression analyses, as well as lasso regression analysis, four CESC-independent prognostic genes were identified, and a prognostic model was constructed, dividing them into high and low-risk groups. The correlation between these genes and immune cells and immune functions were explored through ssGSEA enrichment analysis, revealing a close association between the high-risk group and processes such as angiogenesis and epithelial-mesenchymal transition. Furthermore, using public databases and qRT-PCR experiments, significant differences in CXCL8 expression between normal cervical cells and cervical cancer cells were discovered. Subsequently, a CXCL8 knockdown plasmid was constructed, and the efficiency of CXCL8 knockdown was validated in two CESC cell lines, MEG-01 and HCE-1. Through CCK-8, scratch, and Transwell assays, it was confirmed that CXCL8 knockdown could inhibit the proliferation, invasion, and migration abilities of CESC cells. Targeting CXCL8 holds promise for personalized therapy for CESC, providing a strong theoretical basis for achieving clinical translation.

Comprehensive analysis of m6 A methylome and transcriptome by Nanopore sequencing in clear cell renal carcinoma

N6 -methyladenosine (m6 A) is the most prevalent epigenetic modification on eukaryotic messenger RNAs. Recent studies have focused on elucidating the key role of m6 A modification patterns in tumor progression. However, the relationship between m6 A and transcriptional regulation remains elusive. Nanopore technology enables the quantification of m6 A levels at each genomic site. In this study, a pair of tumor tissues and adjacent normal tissues from clear cell renal cell carcinoma (ccRCC) surgical samples were collected for Nanopore direct RNA sequencing. We identified 9644 genes displaying anomalous m6 A modifications, with 5343 genes upregulated and 4301 genes downregulated. Among these, 5224 genes were regarded as dysregulated genes, encompassing abnormal regulation of both m6 A modification and RNA expression. Gene Set Enrichment Analysis revealed an enrichment of these genes in pathways related to renal system progress and fatty acid metabolic progress. Furthermore, the χ2 test demonstrated a significant association between the levels of m6 A in dysregulated genes and their transcriptional expression levels. Additionally, we identified four obesity-associated genes (FTO, LEPR, ADIPOR2, and NPY5R) among the dysregulated genes. Further analyses using public databases revealed that these four genes were all related to the prognosis and diagnosis of ccRCC. This study introduced the novel approach of employing conjoint analysis of m6 A modification and RNA expression based on Nanopore sequencing to explore potential disease-related genes. Our work demonstrates the feasibility of the application of Nanopore sequencing technology in RNA epigenetic regulation research and identifies new potential therapeutic targets for ccRCC.

Phospholipase A2 enzymes PLA2G2A and PLA2G12B as potential diagnostic and prognostic biomarkers in cholangiocarcinoma

BACKGROUND

Phospholipase A2 (PLA2) enzymes are pivotal in various biological processes, such as lipid mediator production, membrane remodeling, bioenergetics, and maintaining the body surface barrier. Notably, these enzymes play a significant role in the development of diverse tumors.

AIM

To systematically and comprehensively explore the expression of the PLA2 family genes and their potential implications in cholangiocarcinoma (CCA).

METHODS

We conducted an analysis of five CCA datasets from The Cancer Genome Atlas and the Gene Expression Omnibus. The study identified differentially expressed genes between tumor tissues and adjacent normal tissues, with a focus on PLA2G2A and PLA2G12B. Gene Set Enrichment Analysis was utilized to pinpoint associated pathways. Moreover, relevant hub genes and microRNAs for PLA2G2A and PLA2G12B were predicted, and their correlation with the prognosis of CCA was evaluated.

RESULTS

PLA2G2A and PLA2G12B were discerned as differentially expressed in CCA, manifesting significant variations in expression levels in urine and serum between CCA patients and healthy individuals. Elevated expression of PLA2G2A was correlated with poorer overall survival in CCA patients. Additionally, the study delineated pathways and miRNAs associated with these genes.

CONCLUSION

Our findings suggest that PLA2G2A and PLA2G12B may serve as novel potential diagnostic and prognostic markers for CCA. The increased levels of these genes in biological fluids could be employed as non-invasive markers for CCA, and their expression levels are indicative of prognosis, underscoring their potential utility in clinical settings.

Exploring Immune-Related Gene Profiling and Infiltration of Immune Cells in Cervical Squamous Cell Carcinoma and Endocervical Adenocarcinoma

Cervical cancer is a widespread malignancy among women, leading to a substantial global health impact. Despite extensive research, our understanding of the basic molecules and pathogenic processes of cervical squamous cell carcinoma is still insufficient. This investigation aims to uncover immune-related genes linked to CESC and delineate their functions. Leveraging data from the GEO and ImmPort databases, a total of 22 immune-related genes were identified. Multiple tools, including DAVID, the human protein atlas, STRING, GeneMANIA, and TCGA, were employed to delve into the expression and roles of these immune genes in CESC, alongside their connections to the disease's pathological features. Through RT-PCR, the study confirmed notable disparities in CXCL8 and CXCL10 mRNA expression between CESC and normal cervical tissue. The TCGA dataset's immune-related information reinforced the association of CXCL8 and CXCL10 with immune infiltration in CESC. This research sheds light on the potential of CXCL8 and CXCL10 as promising therapeutic targets and essential prognostic factors for individuals diagnosed with CESC.

Circulating inflammatory cytokines in relation to the risk of renal cell carcinoma: A gender-specific two-sample Mendelian randomization study

BACKGROUND

Currently there is no specific molecular biomarker for the diagnosis and treatment of renal cell carcinoma (RCC). Here we performed a gender-specific two-sample Mendelian randomization analysis to systematically assess the effects of circulating cytokines on RCC.

METHODS

We have employed cis-quantitative trait loci as instrumental variables for the protein levels and expression of circulating cytokines. We estimated the causal effects of circulating cytokines on RCC risk in males and females with several Mendelian randomization methods.

RESULTS

We observed a significant causal effect of Eotaxin on the increased risk of RCC in males (Odds ratio [OR] = 2.546, 95% confidence interval [CI] = 1.617-4.010, p value = 5.496 × 10-5), but not in females (OR = 1.352, 95% CI = 0.766-2.388, p value = 0.298). Besides, we also identified several cytokines as potentially associated with RCC in males including RANTES, MCP3, PDGFbb, TRAIL, and several other cytokines as potentially associated with RCC in females including sICAM and SCGFb.

CONCLUSION

Our study highlighted that a higher level of circulating Eotaxin is causally associated with an increased risk of RCC in males but not in females. Further studies are needed to elucidate the exact mechanism and its potential application in the prognosis and treatment of RCC.

Bioinformatics identification and validation of aging‑related molecular subtype and prognostic signature in breast cancer

Patients with metastatic breast cancer have a poor clinical outcome, accounting for more than 90 percent of breast cancer-related deaths. Aging could regulate many biological processes in malignancies by regulating cell senescence. The role of aging has not been fully clarified. Consensus cluster analysis was performed to differentiate The Cancer Genome Atlas (TCGA) breast cancer cases. Least absolute shrinkage and selection operator (LASSO) cox regression analysis was performed to construct an aging-related prognostic signature. A total of 118 differentially expressed aging-related genes (ARGs) was obtained in breast cancer. Consensus clustering analysis identified 3 categories of TCGA-breast cancer with significant difference in prognosis and immune infiltration. We also constructed an aging-related prognostic signature for breast cancer, which had a good performance in predicting the 1-year, 3-year and 5-year OS and disease specific survival (DSS) of breast cancer patients. Further single gene analysis revealed that the expression of PIK3R1 was significantly different in different pT and pN stages of breast cancer. Moreover, low expression of PIK3R1 showed resistance to many drugs based on the data of Genomics of Drug Sensitivity in Cancer (GDSC) and Genomics of Therapeutics Response Portal (CTRP). PIK3R1 played a vital role in many well-known cancer-related pathways. The current study identified 3 clusters of TCGA-breast cancer cases with significant differences in prognosis and immune infiltration. We also constructed an aging-related prognostic signature for breast cancer. However, further in vivo and in vitro studies should be conducted to verify these results.

Integrated bioinformatics analysis for conducting a prognostic model and identifying immunotherapeutic targets in gastric cancer

BACKGROUND

Gastric cancer is the third leading cause of death from cancer worldwide and has a poor prognosis. Practical risk scores and prognostic models for gastric cancer are lacking. While immunotherapy has succeeded in some cancers, few gastric cancer patients benefit from immunotherapy. Immune genes and the tumor microenvironment (TME) are essential for cancer progression and immunotherapy response. However, the roles of immune genes and the tumor microenvironment in immunotherapy remain unclear. The study aimed to construct a prognostic prediction model and identify immunotherapeutic targets for gastric cancer (GC) patients by exploring immune genes and the tumor microenvironment.

RESULTS

An immune-related risk score (IRRS) model, including APOH, RNASE2, F2R, DEFB126, CXCL6, and CXCL3 genes, was constructed for risk stratification. Patients in the low-risk group, which was characterized by elevated tumor mutation burden (TMB) have higher survival rate. The risk level was remarkably correlated with tumor-infiltrating immune cells (TIICs), the immune checkpoint molecule expression, and immunophenoscore (IPS). CXCL3 and CXCL6 were significantly upregulated in gastric cancer tissues compared with normal tissues using the UALCAN database and RT-qPCR. The nomogram showed good calibration and moderate discrimination in predicting overall survival (OS) at 1-, 3-, and 5- year for gastric cancer patients using risk-level and clinical characteristics.

CONCLUSION

Our findings provided a risk stratification and prognosis prediction tool for gastric cancer patients and further the research into immunotherapy in gastric cancer.

Bioinformatics identification and validation of aging-related molecular subtype and prognostic signature in sarcoma

Aging could regulate many biological processes in malignancies by regulating cell senescence. Consensus cluster analysis was conducted to differentiate TCGA sarcoma cases. LASSO cox regression analysis was performed to construct an aging-related prognostic signature. We identified two categories of TCGA-sarcoma with significant difference in prognosis, immune infiltration and chemotherapy and targeted therapy. Moreover, an aging-related prognostic signature was constructed for sarcoma, which had a good performance in predicting the 3-year and 5-year overall survival of sarcoma patients. We also identified a lncRNA MALAT1/miR-508-3p/CCNA2 regulatory axis for sarcoma. This stratification could provide more evidence for estimating prognosis and immunotherapy of sarcoma.

The Role of Cytokines in Cholesterol Accumulation in Cells and Atherosclerosis Progression

Atherosclerosis is the most common cardiovascular disease and is the number one cause of death worldwide. Today, atherosclerosis is a multifactorial chronic inflammatory disease with an autoimmune component, accompanied by the accumulation of cholesterol in the vessel wall and the formation of atherosclerotic plaques, endothelial dysfunction, and chronic inflammation. In the process of accumulation of atherogenic lipids, cells of the immune system, such as monocytes, macrophages, dendritic cells, etc., play an important role, producing and/or activating the production of various cytokines—interferons, interleukins, chemokines. In this review, we have tried to summarize the most important cytokines involved in the processes of atherogenesis.

LRP2 and DOCK8 Are Potential Antigens for mRNA Vaccine Development in Immunologically 'Cold' KIRC Tumours

The administration of mRNA-based tumour vaccines is considered a promising strategy in tumour immunotherapy, although its application against kidney renal clear cell carcinoma (KIRC) is still at its infancy stage. The purpose of this study was to identify potential antigens and to further select suitable patients for vaccination. Gene expression data and clinical information were retrieved from Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases. GEPIA2 was used to evaluate the prognostic value of selected antigens. The relationship of antigens presenting cell infiltration with antigen expression was evaluated by TIMER, and immune subtypes were determined using unsupervised cluster analysis. Tumour antigens LRP2 and DOCK8, which are associated with prognosis and tumour-infiltrating antigen-presenting cells, were identified in KIRC. A total of six immune subtypes were identified, and patients with immune subtype 1-4 (IS1-4) tumours had an immune 'cold' phenotype, a higher tumour mutation burden, and poor survival. Moreover, these immune subtypes showed significant differences in the expression of immune checkpoint and immunogenic cell death modulators. Finally, the immune landscape of KIRC revealed the immune-related cell components in individual patients. This study suggests that LRP2 and DOCK8 are potential KIRC antigens in the development of mRNA vaccines, and patients with immune subtypes IS1-4 are suitable for vaccination.

CXCL1 confers a survival advantage in Kaposi's sarcoma-associated herpesvirus-infected human endothelial cells through STAT3 phosphorylation

Many cytokines produced by Kaposi's sarcoma-associated herpesvirus (KSHV)-infected cells have been shown to participate in the pathogenesis of KSHV. Determination of the exact role of cytokines in Kaposi's sarcoma (KS) pathogenesis is limited, however, by the difficulty to manipulate the target genes in human endothelial cells. In this study, we sought to elucidate the role of cytokines in KSHV-infected human immortalized endothelial cell line (HuARLT cells) by knockout (KO) of the corresponding target genes using the CRISPR/Cas9 system. The cytokine production profile of KSHV-infected HuARLT cells was analyzed using a protein array, and several cytokines were found to be highly upregulated following KSHV infection. This study focused on CXCL1, which was investigated by knocked out in HuARLT cells. KSHV-infected CXCL1 KO cells underwent increased cell death compared to KSHV-infected wild-type (WT) cells and mock-infected CXCL1 KO cells. Lytic replication was not observed in KSHV-infected WT nor CXCL1 KO cells. Phosphorylation of STAT3 was significantly suppressed in KSHV-infected CXCL1 KO cells. Additionally, inhibitors of STAT3 and CXCL1 induced cell death in KSHV-infected endothelial cells. Our results show that CXCL1 production is required for the survival of KSHV-infected endothelial cells, and the CXCL1 to STAT3 phosphorylation signaling pathway may be a therapeutic target for KS.

Bioinformatics analysis of markers based on m6 A related to prognosis combined with immune invasion of renal clear cell carcinoma

The incidence rate of renal cell carcinoma (RCC) is about 3% of all adult cancers. Of these, the Kidney clear cell renal cell carcinoma (KIRC) is the most common type, accounting for about 70%-75% of RCC. KIRC is difficult to be detected in time clinically. KIRC still has no effective treatment at this stage. We combined high-throughput bioinformatics analysis to obtained the structural sequence transcriptome data, relevant clinical information, and m⁶ A gene map of KIRC patients from genomics TCGA database. Pearson's correlation analysis was used to explore m⁶ A related gene long noncoding RNAs (lncRNAs), and then univariate Cox regression analysis was performed to screen the prognostic role of KIRC patients. Lasso-Cox regression was performed to establish the lncRNAs risk model associated with m⁶ A.LINC02154 and AC016773.2, Z98200.2, AL161782.1, EMX2OS, AC021483.2, CD27-AS1, AC006213.3 were iidentif. Compared with the low-risk group, the overall survival of patients in the high-risk group was significantly worse. Analyzing whether there are differences in immune cells between high-risk and low-risk subgroups. There were CD4 memory resting, Monocytes, Macrophages M1, Dendritic cells activated, Mast cells resting, which had higher infiltrations in the low-risk group. We performed Go enrichment analysis, Kyoto Encyclopedia of Genes and Genomes enrichment analysis and gene set enrichment analysis enrichment analysis. Overall, our results suggest that the component of m6A-related lncRNAs in the prognostic signal may be a key mediator in the immune microenvironment of KIRC, which represents a promising therapeutic effect.

Identification and validation of the molecular subtype and prognostic signature for clear cell renal cell carcinoma based on neutrophil extracellular traps

Background: Renal cell carcinoma (RCC) is one of the most common cancers, with an annual incidence of nearly 400,000 cases worldwide. Increasing evidence has also demonstrated the vital role of neutrophil extracellular traps (NETs) in cancer progression and metastatic dissemination. Methods: Consensus cluster analysis was performed to determine the number of ccRCC subtypes. The Kruskal-Wallis test or Student t-test was performed to evaluate the difference of infiltrating immune cell and gene expression in different groups. The Kaplan-Meier (KM) method was used to draw the survival curve. LASSO cox regression analysis was conducted to construct a NET-related prognostic signature. We also constructed a lncRNA-miRNA-mRNA regulatory axis by several miRNA and lncRNA target databases. Results: A total of 23 differentially expressed NET-related genes were obtained in ccRCC. Three clusters of ccRCC cases with significant difference in prognosis, immune infiltration, and chemotherapy and targeted therapy were identified. LASSO Cox regression analysis identified a NET-related prognostic signature including six genes (G0S2, DYSF, MMP9, SLC22A4, SELP, and KCNJ15), and this signature had a good performance in predicting the overall survival of ccRCC patients. The expression of prognostic signature genes was significantly correlated with the pTMN stage, immune infiltration, tumor mutational burdens, microsatellite instability, and drug sensitivity of ccRCC patients. MMP9 was identified as the hub gene. We also identified the lncRNA UBA6-AS1/miR-149-5p/MMP9 regulatory axis for the progression of ccRCC. Conclusion: Collectively, the current study identified three molecular clusters and a prognostic signature for ccRCC based on neutrophil extracellular traps. Integrative transcriptome analyses plus clinical sample validation may facilitate the biomarker discovery and clinical transformation.

CXCL8 Up-Regulated LSECtin through AKT Signal and Correlates with the Immune Microenvironment Modulation in Colon Cancer

Simple Summary

Patients with high expression of CXCL8 are not sensitive to immune checkpoint inhibitors (ICIs) treatment, but the mechanism is unclear. LSECtin is the immune checkpoint ligand of LAG3, and is considered as an important factor of ICIs resistance. This study confirmed the role of CXCL8 and LSECtin in immune microenvironment modulation of colon cancer. The expression of CXCL8 is positively correlated with more than 40 immune checkpoints. CXCL8 could up-regulate LSECtin through AKT signal and promoted the proliferation and invasion ability of colon cancer. These results may be important reasons for the primary drug resistance of ICIs in colon cancer.

Abstract

Background: The role of CXCL8 and LSECtin in colon cancer liver metastasis and immune checkpoint inhibitors (ICIs) treatment effect were widely recognized. However, the regulatory role of CXCL8 on LSECtin is still unclear. Methods: The expression of CXCL8 or LSECtin was analyzed by TCGA database, and verified by GES110225 and clinical samples. The relationship between the expression of CXCL8 or LSECtin and immune cells infiltration, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, Gene Ontology (GO) items, stromal score, Estimation of STromal and Immune cells in MAlignant Tumours (ESTIMAT) immune score, tumor mutation burden (TMB), mismatch repair gene and immune checkpoints expression were analyzed by Spearman. The effects of CXCL8 on LSECtin expression, proliferation, and invasion ability were clarified by recombinant CXCL8 or CXCL8 interfering RNA. Results: In colon cancer, the expression of CXCL8 was higher, but LSECtin was lower than that in normal mucosa. The expression of CXCL8 or LSECtin was significantly positively correlated with immune cells infiltration, stromal score, ESTIMATE immune score, TMB, and immune checkpoints expression. The expression of LSECtin was closely related to the cytokine-cytokine receptor interaction pathway and response of chemokine function, such as CXCL8/CXCR1/2 pathway. There was a significant positive correlation between the expression of CXCL8 and LSECtin in colon cancer. CXCL8 up-regulated LSECtin through AKT signal and promoted the proliferation and invasion ability of colon cancer. Conclusions: CXCL8 up-regulated LSECtin by activating AKT signal and correlated with the immune microenvironment modulation in colon cancer.

CXCL2 acts as a prognostic biomarker and associated with immune infiltrates in stomach adenocarcinoma

BACKGROUND

STAD ranked 5th most common in the incidence of malignant tumors and 3rd most common in the death rate of cancer worldwide. CXC chemokines affect the biological progress of various tumors, resulting in therapeutic failure. The role of CXCL2 in STAD was still a mystery.

METHODS

The expression, prognostic value, and clinical function of CXCL2 were analyzed using several online bioinformatics tools and clinical tissues.

RESULTS

CXCL2 level was significantly upregulated in STAD tissues. Strong correlation was obtained between CXCL2 level and immune cells as well as immune biomarkers. High CXCL2 expression in STAD was correlated with a favorable prognosis. Further analysis revealed that CXCL2, pTNM stage and age were independent factors affecting the prognosis of STAD patients. A predictive nomogram indicated that the calibration plots for the 1-year, 3-year and 5-year OS rates were predicted relatively well compared with an ideal model in the entire cohort. Validation analysis revealed that CXCL2 expression was upregulated in STAD and high CXCL2 level had a better overall survival. CXCL2 was associated with resistance to numerous drugs or small molecules in STAD.

CONCLUSIONS

We identified CXCL2 as a novel therapeutic target and associated with immune infiltration in STAD.

Identification of core genes and pathways between geriatric multimorbidity and renal insufficiency: potential therapeutic agents discovered using bioinformatics analysis

Background

Geriatric people are prone to suffer from multiple chronic diseases, which can directly or indirectly affect renal function. Through bioinformatics analysis, this study aimed to identify key genes and pathways associated with renal insufficiency in patients with geriatric multimorbidity and explore potential drugs against renal insufficiency.

Methods

The text mining tool Pubmed2Ensembl was used to detect genes associated with the keywords including "Geriatric", "Multimorbidity" and "Renal insufficiency". The GeneCodis program was used to specify Gene Ontology (GO) biological process terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Protein–protein interaction (PPI) networks were constructed using STRING and visualized in Cytoscape. Module analysis was performed using CytoHubba and Molecular Complex Detection (MCODE) plugins. GO and KEGG analysis of gene modules was performed using the Database for Annotation, Visualization and Integrated Discover (DAVID) platform database. Genes clustered in salient modules were selected as core genes. Then, the functions and pathways of core genes were visualized using ClueGO and CluePedia. Finally, the drug-gene interaction database was used to explore drug-gene interactions of the core genes to identify drug candidates for renal insufficiency in patients with geriatric multimorbidity.

Results

Through text mining, 351 genes associated with "Geriatric", "Multimorbidity" and "Renal insufficiency" were identified. A PPI network consisting of 216 nodes and 1087 edges was constructed and CytoHubba was used to sequence the genes. Five gene modules were obtained by MCODE analysis. The 26 genes clustered in module1 were selected as core candidate genes primarily associated with renal insufficiency in patients with geriatric multimorbidity. The HIF-1, PI3K-Akt, MAPK, Rap1, and FoxO signaling pathways were enriched. We found that 21 of the 26 selected genes could be targeted by 34 existing drugs.

Conclusion

This study indicated that CST3, SERPINA1, FN1, PF4, IGF1, KNG1, IL6, VEGFA, ALB, TIMP1, TGFB1, HGF, SERPINE1, APOA1, APOB, FGF23, EGF, APOE, VWF, TF, CP, GAS6, APP, IGFBP3, P4HB, and SPP1 were key genes potentially involved with renal insufficiency in patients with geriatric multimorbidity. In addition, 34 drugs were identified as potential agents for the treatment and management of renal insufficiency.

Systematic investigation of the clinical significance and prognostic value of the CBXs in esophageal cancer

Esophageal cancer (ESCA), one of the most aggressive malignant tumors, has been announced to be the ninth most common cancer and the sixth leading cause of cancer-related death in the world. Chromobox family members (CBXs) are important epigenetic regulators which are related with the transcription of target genes. The role of CBXs in carcinomas has been reported in many studies. However, the function and prognostic value of different CBXs in ESCA are still largely unknown. In this article, we first performed differential expression analysis through several methods including Oncomine and Gene Expression Profiling Interactive Analysis. The results led us to determine the differential expression of CBXs in pan-cancer, especially ESCA. Then we evaluated the prognostic value of different CBX messenger RNA (mRNA) expression in patients with ESCA through the Kaplan-Meier plotter and the Human Protein Atlas database. In addition, we used cBioPortal to explore all genetic alterations and mutations in the CBXs in ESCA. Simultaneously, the correlation between its expression and the level of immune infiltration of ESCA was visualized by TIMER. Finally, the biological function of CBXs in ESCA is obtained through Biological Enrichment Analysis including gene ontology and Kyoto Encyclopedia of Genes and Genomes. The expression levels of CBX3/4/5 and CBX8 in ESCA tissues increased significantly and the expression level of CBX7 decreased through differential expression analysis. Additionally, CBX1 is significantly related to the clinical cancer stage and disease-free survival of ESCA patients. The high mRNA expression of CBX4 is related to the short overall survival of patients with esophageal squamous cell carcinoma, and the high mRNA expression of CBX3/7/8 is related to the short overall survival of patients with esophageal adenocarcinoma, indicating that CBX1/3/4/7/8 may be a potential prognostic biomarker for the survival of ESCA patients. Besides, the expression of CBXs is significantly related to the infiltration of a variety of immune cells, including six types of CD4-positive T-lymphocytes, macrophages, neutrophils, bursindependentlymphocyte, CD8-positive T-lymphocytes cells and dendritic cells in ESCA. Moreover, we found that CBXs are mainly associated with the inhibition of cell cycle and apoptosis pathway. Further, enrichment analysis indicated that CBXs and correlated genes were enriched in mismatch repair, DNA replication, cancer pathways, and spliceosomes. Our research may provide new insights into the choice of prognosis biomarkers of the CBXs in ESCA.

A Regulatory Network Analysis of the Importance of USP15 in Breast Cancer Metastasis and Prognosis

Background

Ubiquitin-specific protease15(USP15), is the 16th identified protease in the USP family and is a key protein in tumorigenesis. However, the predictive value and regulatory mechanism of USP15 in breast cancer are unclear.

Methods

The GEPIA, UALCAN, GeneMANIA, and STRING databases were applied to explore the expression of USP15 in breast cancer and associated proteins. In addition, the TIMER database was evaluated for immune infiltration patterns. Moreover, protein immunoblotting assay, cell scratching assay, small compartment invasion assay, 3D stromal gel assay, immunoprecipitation assay, and immunohistochemistry (IHC) were used to USP15 regulatory mechanisms in breast cancer.

Results

In BRCA, several databases, including GEPIA and UALCAN, describe the upregulation of total protein levels and USP15 phosphorylation. In addition, the expression of USP15 was significantly correlated with gender and clinical stage. Overall survival (OS) was lower in patients with high USP15 expression. Functional network analysis showed that USP15 is involved in tumor-associated pathways, DNA replication, and cell cycle signaling through TGFβRI. In addition, USP15 expression was positively correlated with immune infiltration, including immune score, mesenchymal score, and several tumor-infiltrating lymphocytes (TIL). In addition, IHC results further confirmed the high expression of USP15 in breast cancer and its prognostic potential.

Conclusions

These findings demonstrate that high USP15 expression indicates poor prognosis in BRCA and reveal potential regulatory networks and the positive relationship with immune infiltration. Thus, USP15 may be an attractive predictor for BRCA.

Identification of cervical squamous cell carcinoma feature genes and construction of a prognostic model based on immune-related features

As heterogeneity of cervical squamous cell carcinoma (CSCC), prognosis assessment for CSCC patients remain challenging. To develop novel prognostic strategies for CSCC patients, associated biomarkers are urgently needed. This study aimed to cluster CSCC samples from a molecular perspective. CSCC expression data sets were obtained from The Cancer Genome Atlas and based on the accessed expression profile, a co-expression network was constructed with weighted gene co-expression network analysis to form different gene modules. Tumor microenvironment was evaluated using ESTIMATE algorithm, observing that the brown module was highly associated with tumor immunity. CSCC samples were clustered into three subtypes by consensus clustering based on gene expression profiles in the module. Gene set variation analysis showed differences in immune-related pathways among the three subtypes. CIBERSORT and single-sample gene set enrichment analysis analyses showed the difference in immune cell infiltration among subtype groups. Also, Human leukocyte antigen protein expression varied considerably among subtypes. Subsequently, univariate, Lasso and multivariate Cox regression analyses were performed on the genes in the brown module and an 8-gene prognostic model was constructed. Kaplan-Meier analysis illuminated that the low-risk group manifested a favorable prognosis, and receiver operating characteristic curve showed that the model has good predictive performance. qRT-PCR was used to examine the expression status of the prognosis-associated genes. In conclusion, this study identified three types of CSCC from a molecular perspective and established an effective prognostic model for CSCC, which will provide guidance for clinical subtype identification of CSCC and treatment of patients.

A novel chemokine-based signature for prediction of prognosis and therapeutic response in glioma

AIMS

Gliomas are the primary malignant brain tumor and characterized as the striking cellular heterogeneity and intricate tumor microenvironment (TME), where chemokines regulate immune cell trafficking by shaping local networks. This study aimed to construct a chemokine-based gene signature to evaluate the prognosis and therapeutic response in glioma.

METHODS

In this study, 1024 patients (699 from TCGA and 325 from CGGA database) with clinicopathological information and mRNA sequencing data were enrolled. A chemokine gene signature was constructed by combining LASSO and SVM-RFE algorithm. GO, KEGG, and GSVA analyses were performed for function annotations of the chemokine signature. Candidate mRNAs were subsequently verified through qRT-PCR in an independent cohort including 28 glioma samples. Then, through immunohistochemical staining (IHC), we detected the expression of immunosuppressive markers and explore the role of this gene signature in immunotherapy for glioma. Lastly, the Genomics of Drug Sensitivity in Cancer (GDSC) were leveraged to predict the potential drug related to the gene signature in glioma.

RESULTS

A constructed chemokine gene signature was significantly associated with poorer survival, especially in glioblastoma, IDH wildtype. It also played an independent prognostic factor in both datasets. Moreover, biological function annotations of the predictive signature indicated the gene signature was positively associated with immune-relevant pathways, and the immunosuppressive protein expressions (PD-L1, IBA1, TMEM119, CD68, CSF1R, and TGFB1) were enriched in the high-risk group. In an immunotherapy of glioblastoma cohort, we confirmed the chemokine signature showed a good predictor for patients' response. Lastly, we predicted twelve potential agents for glioma patients with higher riskscore.

CONCLUSION

In all, our results highlighted a potential 4-chemokine signature for predicting prognosis in glioma and reflected the intricate immune landscape in glioma. It also threw light on integrating tailored risk stratification with precision therapy for glioblastoma.

Interleukin 20 receptor subunit beta (IL20RB) predicts poor prognosis and regulates immune cell infiltration in clear cell renal cell carcinoma

BACKGROUND AND OBJECTIVE

Emerging evidence has proven the robust role of tumor mutation burden (TMB) and immune cell infiltration (ICI) in cancer immunotherapy. However, the precise effect of TMB and ICI on clear cell renal cell carcinoma (ccRCC) remains elusive and merits further investigation. Therefore, we aim to identify the TMB-related genes in predicting prognosis and to explore the potential mechanisms of the identified Interleukin 20 receptor subunit beta (IL20RB) in ICI in ccRCC.

METHOD

The relative information of patients with ccRCC was obtained from The Cancer Genome Atlas database (TCGA). Immune-related genes were downloaded from the Immunology Database and Analysis Portal database. Cox regression analysis was used to identify prognosis-related immune genes for ccRCC. The relationship of IL20RB expression levels with clinicopathological parameters was analyzed using the "limma" and "survival" packages. Gene Expression Omnibus (GEO) and International Cancer Genome Consortium (ICGC) databases were used as external validation. Quantitative Real-time PCR (qRT-PCR) and western blots were used to validate the expression levels of IL20RB in tumor cells. Cell counting kit-8 (CCK-8) assay and colony formation assay were used to examine the effect of IL20RB on the viability of ccRCC cells. Gene set enrichment analysis (GSEA) was introduced for the analysis of IL20RB-related signaling pathways. Tumor Immune Estimation Resource (TIMER) and Tumor and Immune System Interaction Database (TISIDB) were utilized to determine the correlation of IL20RB expression levels with tumor-infiltrating immune cells (TIICs).

RESULTS

IL20RB was significantly overexpressed in different ccRCC tissues and cells. High IL20RB expression in ccRCC patients was associated with short overall survival, high tumor grade, and advanced TNM stage. After knockdown of IL20RB with small interfering RNA (siRNA) technology, ccRCC cells' proliferation was significantly attenuated. Moreover, overexpression of IL20RB could increase the infiltration level of several immune cells, especially T follicular helper cells (Tfh), and overexpressed Tfh cells were correlated with poor prognosis in ccRCC.

CONCLUSIONS

IL20RB may function as an immune-associated therapeutic target for it determines cancer progression and regulates immune cell infiltration in ccRCC.

Integrin-Linked Kinase Expression Characterizes the Immunosuppressive Tumor Microenvironment in Colorectal Cancer and Regulates PD-L1 Expression and Immune Cell Cytotoxicity

Integrin-linked kinase (ILK) has been implicated as a molecular driver and mediator in both inflammation and tumorigenesis of the colon. However, a role for ILK in the tumor microenvironment (TME) and immune evasion has not been investigated. Here, we show a correlation of ILK expression with the immunosuppressive TME and cancer prognosis. We also uncover a role for ILK in the regulation of programmed death-ligand 1 (PD-L1) expression and immune cell cytotoxicity. Interrogation of web-based data-mining platforms, showed upregulation of ILK expression in tumors and adjacent-non tumor tissue of colorectal cancer (CRC) associated with poor survival and advanced stages. ILK expression was correlated with cancer-associated fibroblast (CAFs) and immunosuppressive cell infiltration including regulatory T cells (Treg) and M2 macrophages (M2) in addition to their gene markers. ILK expression was also significantly correlated with the expression of different cytokines and chemokines. ILK expression showed pronounced association with different important immune checkpoints including PD-L1. Deletion of the ILK gene in PD-L1 positive CRC cell lines using a doxycycline inducible-CRISPR/Cas9, resulted in suppression of both the basal and IFNγ-induced PD-L1 expression via downregulating NF-κB p65. This subsequently sensitized the CRC cells to NK92 immune cell cytotoxicity. These findings suggest that ILK can be used as a biomarker for prognosis and immune cell infiltration in colon cancer. Moreover, ILK could provide a therapeutic target to prevent immune evasion mediated by the expression of PD-L1.

ACE2 maybe serve as a prognostic biomarker in breast invasive carcinoma

BACKGROUND

Breast cancer is a frequently occurring malignant tumor in women. Angiotensin-converting enzyme 2 (ACE2) is widely expressed in most organs; however, the association of ACE2 with prognosis and immune infiltration in breast invasive carcinoma (BRCA) remains elusive.

METHODS

We explored the expression level and prognostic value of ACE2 in patients with BRCA using a series of online bioinformatics analysis databases encompassing Oncomine, UALCAN, Kaplan-Meier plotter, TIMER, LinkedOmics, and GEO. qRT-PCR was performed to verify our findings.

RESULTS

Angiotensin-converting enzyme 2 mRNA and protein expression levels were decreased in BRCA tissues, and patients with low ACE2 expression levels had a poor prognosis. DNA promoter methylation of ACE2 significantly downregulated ACE2 expression in BRCA, while the expression of this protein was positively linked to immune infiltration of B cells, CD8+ and CD4+ T cells, neutrophils, and dendritic cells in BRCA tissues. The high expression level of ACE2 in enriched basophils, CD8+ T cells, and type-2 helper T cells, which showed decreasing levels, indicated a better prognosis for BRCA. Enrichment analyses revealed that NF-κB, IL-17, and TNF signaling pathways were highly correlated to ACE2 in BRCA. Verification study revealed that downregulation of ACE2 was associated with a better prognosis in BRCA. Univariate and multivariate analysis confirmed ACE2 expression and clinical stage as independent prognostic factors for breast cancer.

CONCLUSIONS

Angiotensin-converting enzyme 2 may be a potential prognostic biomarker and target for BRCA. Nevertheless, future investigations are needed for validating our findings and promoting the clinical application of ACE2 in BRCA.

Comprehensive analysis of prognostic value and immune infiltration of CXC chemokines in pancreatic cancer

BACKGROUND

The prognosis of pancreatic cancer is poor, with a 5-year survival rate of less than 10%. Studies have shown that chemokines in the tumour microenvironment are often altered, which is associated with immune infiltration and the prognosis and survival of pancreatic cancer patients.

METHODS

Multiomics and bioinformatics tools were used to clarify CXC chemokine expression and its role in the pancreatic ductal adenocarcinoma (PDAC) immune microenvironment.

RESULTS

Most CXC chemokines were upregulated in pancreatic cancer and correlated with patient prognosis. CXC chemokines can activate cancer-related signalling pathways and affect immune infiltration. Furthermore, most CXC chemokines were significantly correlated with the abundance of macrophages, neutrophils and dendritic cells. CXCL5 was selected as a hub gene, and a variety of immune checkpoints, including PD-1/PD-L1 and CTLA-4, were identified.

CONCLUSION

Our study provides novel insights into CXC chemokine expression and its role in the PDAC immune microenvironment. These results can provide more data about prognostic biomarkers and therapeutic targets of PDAC.

RNA-seq and Mitochondrial DNA Analysis of Adrenal Gland Metastatic Tissue in a Patient with Renal Cell Carcinoma

This study aimed to clarify whether genetic mutations participate in renal cell carcinoma (RCC) metastasis to the adrenal gland (AG). Our study analyzed whole mitochondrial gene and ribonucleic acid sequencing (RNA-seq) data from a male patient in his 60s with metastatic RCC. We confirmed common mutation sites in the mitochondrial gene and carried out Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis using RNA-seq data for RCC and adrenal carcinoma. Furthermore, we confirmed the common mutation sites of mitochondrial genes in which the T3394Y (p.H30Y) site transitioned from histidine (His.; H) to tyrosine (Tyr.; Y) in the NADH dehydrogenase subunit 1 (ND1) gene. The R11,807G (p.T350A) site transitioned from threonine (Thr.; T) to alanine (Ala.; A). Additionally, the G15,438R or A (p.G231D) site transitioned from glycine (Gly.; G) to aspartic acid (Asp.; D) in cytochrome b (CYTB). Furthermore, pathway analysis, using RNA-seq, confirmed the common mutant pathway between RCC and adrenal carcinoma as cytokine–cytokine receptor (CCR) interaction. Confirmation of the original mutation sites suggests that transfer to AG may be related to the CCR interaction. Thus, during metastasis to the AG, mitochondria DNA mutation may represent the initial origin of the metastasis, followed by the likely mutation of the nuclear genes.

The initial expression alterations occurring to transcription factors during the formation of breast cancer: Evidence from bioinformatics

BACKGROUND

Breast cancer (BC) is the leading malignancy among women worldwide.

AIM

This work aimed to present a comprehensively bioinformatic analysis of gene expression profiles and to identify the hub genes during BC tumorigenesis, providing potential biomarkers and targets for the diagnosis and therapy of BC.

MATERIALS & METHODS

In this study, multiple public databases, bioinformatics approaches, and online analytical tools were employed and the real-time reverse transcription polymerase chain reaction was implemented.

RESULTS

First, we identified 10, 107, and 3869 differentially expressed genes (DEGs) from three gene expression datasets (GSE9574, GSE15852, and GSE42568, covering normal, para-cancerous, and BC samples, respectively), and investigated different biological functions and pathways involved. Then, we screened out 8, 16, and 29 module genes from these DEGs, respectively. Next, 10 candidate genes were determined through expression and survival analyses. We noted that seven candidate genes JUN, FOS, FOSB, EGR1, ZFP36, CFD, and PPARG were downregulated in BC compared to normal tissues and lower expressed in aggressive types of BC (basal, HER2+ , and luminal B), TP53 mutation group, younger patients, higher stage BC, and lymph node metastasis BC, while CD27, PSMB9, and SELL were upregulated. The present study discovered that the expression levels of these candidate genes were correlated with the infiltration of immune cells (CD8+ T cell, macrophage, natural killer [NK] cell, and cancer-associated fibroblast) in BC, as well as biomarkers of immune cells and immune checkpoints. We also revealed that promoter methylation, amplification, and deep deletion might contribute to the abnormal expressions of candidate genes. Moreover, we illustrated downstream-targeted genes of JUN, FOS, FOSB, EGR1, and ZFP36 and demonstrated that these targeted genes were involved in "positive regulation of cell death", "pathways in cancer", "PI3K-Akt signaling pathway", and so on.

DISCUSSION & CONCLUSION

We presented differential gene expression profiles among normal, para-cancerous, and BC tissues and further identified candidate genes that might contribute to tumorigenesis and progression of BC, as potential diagnostic and prognostic targets for BC patients.

Comprehensive Analysis of the Glycolysis-Related Gene Prognostic Signature and Immune Infiltration in Endometrial Cancer

Glucose metabolic reprogramming and immune imbalance play important roles in the progression of cancers. The purpose of this study is to develop a glycolysis-related prognostic signature for endometrial cancer (EC) and analyze its relationship with immune function. The mRNA expression profiling of the glycolysis-related genes and clinical data of EC patients were downloaded from The Cancer Genome Atlas (TCGA). We identified a glycolysis-related gene prognostic signature for predicting the prognosis of EC by using The Least Absolute Shrinkage and Selection Operator (LASSO) regression and found the patients in the high-risk group had worse survival prognosis. Multivariate Cox regression analysis showed that the gene signature was an independent prognostic factor for EC. The ROC curve confirmed the accuracy of the prognostic signature (AUC = 0.730). Then, we constructed a nomogram to predict the 1–5 years survival rate of EC patients. The association between the gene signature and immune function was analyzed based on the “ESTIMATE” and “CIBERSORT” algorithm, which showed the immune and ESTIMATE scores of patients in the high-risk group were lower, while the low immune and ESTIMATE scores were associated with a worse prognosis of patients. The imbalance of immune cells was also found in the high-risk group. Further, the protein of CDK1, a gene in the signature, was found to be closely related to prognosis of EC and inhibition of CDK1 could inhibit migration and promote apoptosis of EC cells. This study reveals a link between glycolysis-related gene signature and immunity, and provides personalized therapeutic targets for EC.

Toll-Like Receptors Serve as Biomarkers for Early Diagnosis and Prognosis Assessment of Kidney Renal Clear Cell Carcinoma by Influencing the Immune Microenvironment: Comprehensive Bioinformatics Analysis Combined With Experimental Validation

Background: Toll-like receptors (TLRs) are important initiators of innate and acquired immune responses. However, its role in kidney renal clear cell carcinoma (KIRC) remains unclear.

Methods: TLRs and their relationships with KIRC were studied in detail by ONCOMINE, UALCAN, GEPIA, cBioPortal, GeneMANIA, FunRich, LinkedOmics, TIMER and TRRUST. Moreover, we used clinical samples to verify the expressions of TLR3 and TLR4 in early stage of KIRC by real-time fluorescence quantitative polymerase chain reaction (RT-qPCR), flow cytometry (FC) and immunohistochemistry (IHC).

Results: The expression levels of TLRs in KIRC were generally different compared with adjacent normal tissues. Moreover, the expressions of TLR3 and TLR4 elevated significantly in the early stage of KIRC. Overexpressions of TLR1, TLR3, TLR4 and TLR8 in KIRC patients were associated with longer overall survival (OS), while inhibition of TLR9 expression was related to longer OS. Additionally, overexpressions of TLR1, TLR3 and TLR4 in KIRC patients were associated with longer disease free survival (DFS). There were general genetic alterations and obvious co-expression correlation of TLRs in KIRC. The PPI network between TLRs was rather complex, and the key gene connecting the TLRs interaction was MYD88. The GO analysis and KEGG pathway analysis indicated that TLRs were closely related to adaptive immunity, innate immunity and other immune-related processes. RELA, NFKB1, IRF8, IRF3 and HIF1A were key transcription factors regulating the expressions of TLRs. What’s more, the expression levels of all TLRs in KIRC were positively correlated with the infiltration levels of dendritic cells, macrophages, neutrophils, B cells, CD4⁺ T cells and CD8⁺ T cells. Finally, the results of RT-qPCR, FC and IHC confirmed that TLR3 and TLR4 were significantly elevated in the early stage of KIRC.

Conclusion: The occurrence and development of KIRC are closely related to TLRs, and TLRs have the potential to be early diagnostic biomarkers of KIRC and biomarkers for judging the prognosis and immune status of KIRC. This study may provide new insights into the selection of KIRC immunotherapy targets.

Identification of PDCD1 and PDCD1LG2 as Prognostic Biomarkers and Associated with Immune Infiltration in Hepatocellular Carcinoma

Background

Hepatocellular carcinoma has been identified to be among the most prevalent malignancies in the world and has an unfavorable prognosis. Immune checkpoints perform an essential function in many biological processes and are associated with the survival of cancer patients. The function of immune checkpoints remains unknown.

Methods

We used bioinformatic methods to examine the prognostic value of immune checkpoints and the corresponding link to immune infiltration in HCC. qRT-PCR was used to validate the expression of immune checkpoints and their prognostic significance in HCC.

Results

The level of mRNA of SIGLEC15, PDCD1LG2, LAG3, PDCD1, CTLA4 as well as PDCD1LG2 was increased in HCC tissues as opposed to liver tissues. Immune checkpoints were shown to participate in the activation of the apoptotic pathway in HCC patients. The elevated expression of PDCD1 and PDCD1LG2 were shown to have a favorable recurrence-free survival (RFS), progression-free survival (PFS), disease-specific survival (DSS), and overall survival (OS). PDCD1, PDCD1LG2, and pT stage were independent variables that affect the HCC patients’ prognoses as revealed by the multivariate and univariate analyses. A prediction nomogram indicated that the calibration plots for OS rates over three and five years had a stronger predictive performance in the TCGA HCC cohort in contrast with an ideal model. Positive correlations were observed between the PDCD1 and PDCD1LG2 expression and immune biomarkers, immune cells, chemokine receptors, as well as chemokines.

Conclusion

The present research performed a thorough examination of the prognostic significance of immune checkpoints in HCC and its correlation with immune infiltration, which suggested that PDCD1 and PDCD1LG2 were prognostic biomarkers in HCC and related to the immune infiltration.

A Bioinformatic Analysis of Immune-Related Prognostic Genes in Clear Cell Renal Cell Carcinoma Based on TCGA and GEO Databases

Background

Clear cell renal cell carcinoma (ccRCC) is a commonly occurring tumor. Through a deeper understanding of the immune regulatory mechanisms in the tumor microenvironment, immunotherapy may serve as a potential treatment for cancer patients. This study aimed at identifying the survival-related immune cells and hub genes, which could be potential targets for immunotherapy in ccRCC.

Methods

The gene expression profiles and clinical data of ccRCC patients were extracted from UCSC Xena and Gene Expression Omnibus (GEO) databases. Kaplan–Meier (KM) survival and Least Absolute Shrinkage and Selection Operator (LASSO) regression analyses were utilized to select the survival-related tumor-infiltrating immune cells. Multivariate Cox regression was utilized to develop a signature based on the tumor-infiltrating immune cells (TIICs). Based on the signature, the risk score was calculated, following which the samples were divided into high-risk and low-risk groups. Differentially expressed genes (DEGs) between the two risk groups were identified. Functional enrichment analysis was performed and cytoHubba plug-in of Cytoscape was used to identify the hub genes. Multiple datasets were utilized to validate these hub genes, including the Gene Expression Profiling Interactive Analysis (GEPIA), UALCAN, and the Human Protein Atlas (HPA), and the GEO datasets. Finally, a correlation analysis was performed to evaluate the relationship between the hub genes and TIICs.

Results

Four immune survival-related cells, including T cell CD4 memory-activated, T cell regulatory (Tregs), eosinophils, and mast cell resting were identified. Nine immune-specific hub genes were identified, which included APOE, CASR, CTLA4, CXCL8, EGF, F2, KNG1, MMP9, and IL6. Furthermore, these hub genes were significantly correlated with clinical traits and closely associated with some TIICs.

Conclusion

A total of four survival-related immune cell types and nine hub genes were found to be closely associated with ccRCC. These findings may have implications for the development of novel potential immunotherapeutic targets for ccRCC.

CXCL1: Gene, Promoter, Regulation of Expression, mRNA Stability, Regulation of Activity in the Intercellular Space

CXCL1 is one of the most important chemokines, part of a group of chemotactic cytokines involved in the development of many inflammatory diseases. It activates CXCR2 and, at high levels, CXCR1. The expression of CXCL1 is elevated in inflammatory reactions and also has important functions in physiology, including the induction of angiogenesis and recruitment of neutrophils. Due to a lack of reviews that precisely describe the regulation of CXCL1 expression and function, in this paper, we present the mechanisms of CXCL1 expression regulation with a special focus on cancer. We concentrate on the regulation of CXCL1 expression through the regulation of CXCL1 transcription and mRNA stability, including the involvement of NF-κB, p53, the effect of miRNAs and cytokines such as IFN-γ, IL-1β, IL-17, TGF-β and TNF-α. We also describe the mechanisms regulating CXCL1 activity in the extracellular space, including proteolytic processing, CXCL1 dimerization and the influence of the ACKR1/DARC receptor on CXCL1 localization. Finally, we explain the role of CXCL1 in cancer and possible therapeutic approaches directed against this chemokine.

High Expression of CEMIP Correlates Poor Prognosis and the Tumur Microenvironment in Breast Cancer as a Promisingly Prognostic Biomarker

Cell migration-inducing hyaluronidase 1 (CEMIP), a Wnt-related protein and also known as KIAA1199, is implicated in the process of metastatic colonization in a variety of malignant tumors, including breast cancer (BC), which is one of the most frequently diagnosed tumors in women worldwide. In this study, multiple public databases, online analytical tools, and bioinformatics approaches were applied to explore the expression levels, regulatory mechanisms, and biological functions of CEMIP in BC. We illustrated that CEMIP was highly expressed in various kinds of carcinomas, including BC, especially advanced subtypes, and predicted less favorable prognosis (negatively associated with overall survival) in BC patients, which might be an independent prognostic factor. Then, we revealed that the mutation and high expression of CEMIP might lead to it as an oncogene. We also demonstrated that TP53 mutation, DNA hypo-methylation, and the expression changes of three potential upstream transcription factors (EZH2, EGR1, and JUN) of CEMIP were likely to cause the hyperexpression of CEMIP in BC. Moreover, our findings suggested that CEMIP might exert its carcinogenic roles in the tumor microenvironment via participation in the extracellular matrix formation, increasing cancer-associated fibroblast (CAF), M2 macrophage, and neutrophil infiltration and decreasing CD8⁺ T cell infiltration. In summary, our study provided more solid evidence for CEMIP as a prognostic and metastatic biomarker and a potential therapeutic target in BC. Of course, these findings also need more confirmations of basic experiments and further clinical trials in the future.

A Comprehensive Bioinformatic Analysis of NOTCH Pathway Involvement in Stomach Adenocarcinoma

Background

Activation of NOTCH signaling pathways, which are key regulators of multiple cellular functions, has been frequently implicated in cancer pathogenesis, and NOTCH inhibitors have received much recent focus in the context of cancer therapeutics. However, the role and possible involvement of NOTCH pathways in stomach adenocarcinoma (STAD) are unclear. Here, putative regulatory mechanisms and functions of NOTCH pathways in STAD were investigated.

Methods

Publicly available data from the TCGA-STAD database were utilized to explore the involvement of canonical NOTCH pathways in STAD by analyzing RNA expression levels of NOTCH receptors, ligands, and downstream genes. Statistical analysis of the data pertaining to cancer and noncancerous samples was performed using R software packages and public databases/webservers.

Results

Significant differential gene expression between control and STAD samples was noted for all NOTCH receptors (NOTCH1, 2, 3, and 4), the delta-like NOTCH ligands (DLL-3 and 4), and typical downstream genes (HES1 and HEY1). Four genes (NOTCH1, NOTCH2, NOTCH3, and HEY1) presented prognostic values for the STAD outcome in terms of overall survival. Functional enrichment analysis indicated that NOTCH family genes-strongly correlated genes were mainly enriched in several KEGG signaling pathways such as the PI3K-Akt signaling pathway, human papillomavirus infection, focal adhesion, Rap1 signaling pathway, and ECM-receptor interaction. Gene set enrichment analysis (GSEA) results showed that NOTCH family genes-significantly correlated genes were mainly enriched in four signaling pathways, ECM (extracellular matrix), tumor angiogenesis, inflammatory response, and immune regulation.

Conclusions

NOTCH family genes may play an essential role in the progression of STAD by modulating immune cells and mediating ECM synthesis, angiogenesis, focal adhesion, and PI3K-Akt signaling. Multiple NOTCH family genes are valuable candidate biomarkers or therapeutic targets for the management of STAD.

Prognostic and immune infiltrates for the Chromobox (CBX) protein family in human pancreatic adenocarcinoma

Background

Globally, pancreatic adenocarcinoma (PAAD) is among the most prevalent malignant tumors. The Chromobox (CBX) protein family is a vital component of epigenetic regulatory complexes that have vital biological roles. The biological functions, immune infiltration, expression levels, and the prognostic significance of CBX proteins in PAAD have not yet been established.

Methods

Using bioinformatics tools, such as the Gene Expression Profiling Interactive Analysis (GEPIA), Oncomine, Kaplan-Meier plotter, GeneMANIA, cBioPortal, TIMER and R, we evaluated the prognostic importance, expression levels, gene alterations, risk factors, and immune cell infiltration levels of CBXs in PAAD patients. The expression levels of CBX3 in clinical-pathological samples were also confirmed by immunohistochemistry.

Results

In PAAD tissues, CBX1, CBX3, CBX5, and CBX8 expressions were high. High CBX1, CBX5, CBX6, and CBX7 levels were correlated with tumor stages. Elevated CBX2, CBX6, CBX7, and CBX8 messenger ribonucleic acid (mRNA) levels were markedly correlated with better overall survival (OS) outcomes for pancreatic cancer patients, while elevated CBX3 was markedly correlated with short OS outcomes. In particular, we have validated the differential expression of CBX3 on clinical specimens using immunohistochemical methods. A multivariate logistic analysis revealed that elevated mRNA expression levels of CBX3 and suppressed CBX8 levels were independently associated with short OS outcomes for pancreatic cancer patients. The roles of CBXs and their neighboring proteins are associated with a negatively regulated cell cycle, histone binding, polycomb group protein complexes, and the regulation of pluripotent stem cell signaling pathways. Additionally, CBX levels were found to be markedly associated with immune infiltrates, and found that the immune infiltration score of CBX3 was differentially expressed in cell lines such as CD8 T cells, NK cells, Mast cells and T helper cells.

Conclusions

CBX3/8 is a potential marker for prognostic outcomes in PAAD patients.

Screening of CXC chemokines in the microenvironment of ovarian cancer and the biological function of CXCL10

BACKGROUND

This study aims to screen and identify the biological functions and prognostic value of CXC chemokines in ovarian cancer (OC) through bioinformatics and molecular biology methods, and to provide data support for the selection of biomarkers and prognostic analysis of OC.

METHODS

In this study, GEO, ONCOMINE, GEPIA, cBioPortal, GeneMANIA, Metascape, STRING, TRRUST, and TIMER databases were used to study CXC chemokines. Angiogenesis and T cell killing assay were used to detect the effect of CXCL10 on tumor cell immunity and angiogenesis. Real-time quantitative PCR (qRT-PCR), immunoblotting, and ectopic tumor formation experiments were used to verify the effect of CXCL10 on ovarian cancer tumors.

RESULTS

We found that CXCL1, CXCL10, CXCL11, CXCL13, and CXCL14 were significantly upregulated in OC samples compared with normal tissues. Our data showed that there was a relationship between the expression of CXC chemokines and the infiltration of six types of immune cells significant correlation. In vitro assay confirmed that overexpression of CXCL10 could enhance the killing effect of T cells and inhibit angiogenesis. Further in vivo assay had shown that CXCL10 could affect the progression of ovarian cancer by increasing the expression of cytotoxic T cells and inhibiting angiogenesis.

CONCLUSION

In conclusion, we hope that our data will provide new insights into the development of immunotherapy and the selection of prognostic markers for patients with OC.

Gene signatures with therapeutic value: emerging perspective for personalized immunotherapy in renal cancer

Renal cancer is one of the deadliest urogenital diseases. In recent years, the advent of immunotherapy has led to significant improvement in the management of patients with renal cancer. Although cancer immunotherapy and its combinations had benefited numerous patients, several challenges need to be addressed. Apart from the high costs of treatment, the lack of predictive biomarkers and toxic side-effects have impeded its wider applicability. To address these issues, new biomarkers are required to predict responsiveness and design personalized treatment strategies. Recent advances in the field of single-cell sequencing and multi-dimensional spatial transcriptomics have identified clinically relevant subtypes of renal cancer. Furthermore, there is emerging potential for gene signatures based on immune cells, non-coding RNAs, and pathways such as metabolism and RNA modification. In this review article, we have discussed recent progress in the identification of gene signatures with predictive and prognostic potential in renal cancer.

UGT1A Gene Family Members Serve as Potential Targets and Prognostic Biomarkers for Pancreatic Cancer

Background

Pancreatic cancer (PC) is one of the most common cancers worldwide, with high mortality. The UGT1A gene family plays important roles in pharmacology and toxicology, contributing to interindividual differences in drug disposition. However, mRNA expression and prognostic value of the UGT1A gene family in PC have not been identified.

Methods

Oncomine, GEPIA2, DAVID 6.8, Metascape, Kaplan-Meier plotter, cBioPortal, GeneMANIA, TRRUST v2, TIMER, and R software were used in our study.

Results

The transcriptional levels of UGT1A1/3/6/8/9/10 in PC tissues were significantly higher than those in normal tissues. These results were further validated using five pairs of PC tumor tissues and adjacent nontumor tissues. A significant correlation was found between the expression of UGT1A1/6/10 and the pathological stage of PC. PC patients with lower transcriptional levels of UGT1A1/4/5/6/10 were associated with a better prognosis. The differentially expressed UGT1A gene family functions were primarily related to the glucuronidation pathway, cytokine-cytokine receptor interactions, and the ILK signaling pathway. Our data suggest that HNF1A, AHR, and CDX2 are key transcription factors for the UGT1A gene family. Furthermore, the expression levels of UGT1A1/3/8/9/10 were positively correlated with the activities of tumor-infiltrating immune cells, especially B cells. The expression levels of UGT1A6/9 were negatively correlated with macrophage infiltration levels.

Conclusions

These results suggest that the UGT1A gene family could serve as a potential prognostic biomarker and target for PC. However, future studies are required to validate our findings and promote the clinical utility of the UGT1A gene family in PC.

Comprehensive Analysis of the Prognostic Value and Immune Function of Immune Checkpoints in Stomach Adenocarcinoma

Background

Stomach adenocarcinoma (STAD) is one of the most prevalent malignances and ranks fifth in incidence and third in the cancer-related deaths among all malignances. The prognosis of STAD is poor. Immunotherapy based on immune checkpoint blockade is ever-increasingly suggested as the most promising therapy strategy for STAD. However, the prognosis and therapy value of immune checkpoints in STAD is far from clarified.

Methods

In our study, bioinformatics methods were performed to explore the expression and prognosis value of immune checkpoints in STAD and their association with immune infiltration. qRT-PCR was performed to verify our result.

Results

Most of the immune checkpoints were upregulated in STAD. There were lots of genetic mutations among immune checkpoints in STAD, including missense_mutation, frame_shift_del et al. Interestingly, most of immune checkpoints were associated with drug sensitivity and drug resistance. Moreover, CD274, PVR, LGALS9, ICOSLG and CD70 were associated with the overall survival, post progression survival and first progression in STAD. The univariate and multivariate analysis revealed that CD70, ICOSLG, age, pTNM stage, and radiation therapy were independent factors affecting the prognosis of STAD patients. The expression of ICOSLG and CD70 was correlated with immune cells as well as immune biomarkers, including CD8+ T cells, CD4+ T cells, macrophage, neutrophils and dendritic cells.

Conclusion

All in all, our study performed a comprehensive analysis of the prognostic value and immune function of immune checkpoints in STAD, and our result suggested that immune checkpoint ICOSLG and CD70 serve as prognostic biomarkers and associate with immune infiltration in STAD.

TLR3 Serves as a Prognostic Biomarker and Associates with Immune Infiltration in the Renal Clear Cell Carcinoma Microenvironment

Background

Clear cell renal cancer (KIRC) is one of the most common cancers globally, with a poor prognosis. TLRs play a vital role in anticancer immunity and the regulation of the biological progress of tumour cells. However, the precise role of TLRs in KIRC is still ambiguous.

Methods

Various bioinformatics analysis and clinical validation of tissues were performed to evaluate the prognostic value of TLRs and their correlation with immune infiltration in KIRC.

Results

The expression of TLR2/3/7/8 was increased at both mRNA and protein levels in KIRC. TLRs in KIRC were involved in the activation of apoptosis, EMT, RAS/MAPK, and RTK pathways, as well as the inhibition of the cell cycle and the hormone AR pathway. Drug sensitivity analysis revealed that high expression of TLR3 and low expression of TLR7/9/10 were resistant to most of the small molecules or drugs from CTRP. Enrichment analyses showed that TLRs were mainly involved in innate immune response, toll-like receptor signalling pathway, NF-kappa B signalling pathway, and TNF signalling pathway. Furthermore, a high-level TLR3 expression was associated with a favourable prognosis in KIRC. Validation research further confirmed that TLR3 expression was increased in KIRC tissues, and high TLR3 levels were associated with poor overall survival. Moreover, TLR3 in KIRC showed a positive association with an abundance of immune cells, including B-cells, CD4+ T-cells, CD8+ T-cells, macrophage, neutrophils, and dendritic cells, and the expression of the immune biomarker sets. Several TLR3-associated kinase, miRNA, or transcription factor targets were also identified in KIRC.

Conclusion

Our results indicate that TLR3 serves as a prognostic biomarker and associated with immune infiltration in KIRC. This work lays a foundation for further studies on the role of TLR3 in the carcinogenesis and progression of KIRC.

Exploration of potential therapeutic and prognostic value of CXC chemokines in cervical squamous cell carcinoma and endocervical adenocarcinoma based on bioinformatics analysis

Cervical cancer, as the second most common female malignancy, brings a great health burden to women worldwide. Cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC) are the most common histological subtypes of cervical cancer. CXC chemokines (CXCLs) within the tumor microenvironment can modulate carcinogenesis and progression. The present study aimed to explore the therapeutic and prognostic value of different CXCLs in CESC. ONCOMINE, GEPIA, cBioPortal, TRRUST, GeneMANIA, STRING and TIMER were utilized to explore the expression, mutation and function of CXCLs in CESC, as well as their correlation with pathological and survival features of CESC patients. We found that the mRNA expression levels of CXCL1/8/9/10/11/13/16/17 in CESC were upregulated compared with normal cervical tissues, whereas CXCL12 was downregulated. No significant correlation was found between the expression levels and pathological stage of CESC patients. CESC patients with high expression of CXCL1/2/3/4/5/8 were significantly associated with poor overall survival, additionally, low mRNA level of CXCL3 was associated with better disease-free survival. Besides, a high mutation rate (43%) of CXCLs in CESC was observed. Depicted by co-expression analysis, the expression of CXCL1/2/3/6/8 showed a modest to strong correlation, while that of CXCL9/10/11/13 showed a very strong correlation. Differentially expressed CXCLs primarily functioned in chemokine signaling pathway and inflammation response, such as cell chemotaxis, chemokine activity and chemokine receptor binding. We also found the association of CXCLs with the tumor-infiltration of six types of immune cells (B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils and dendritic cells) in CESC patients. The present study elucidated that CXCLs may have the potential to be novel therapeutic targets and prognosis predictors of CESC patients.

Pan-Cancer Analysis of Prognostic and Immune Infiltrates for CXCs

Simple Summary

CXCs are important genes that regulate inflammation and tumor metastasis. While there are many studies with a focus on individual CXCs, few present a pan-cancer analysis of the whole CXC family. Our results indicate that CXCs are a potential therapeutic target in a variety of tumors and a potential prognostic marker that could improve the survival of cancer patients and the accuracy of prognosis. Meanwhile, we found that CXCs may be involved in diseases caused by intestinal flora.

Abstract

Background: CXCs are important genes that regulate inflammation and tumor metastasis. However, the expression level, prognosis value, and immune infiltration of CXCs in cancers are not clear. Methods: Multiple online datasets were used to analyze the expression, prognosis, and immune regulation of CXCs in this study. Network analysis of the Amadis database and GEO dataset was used to analyze the regulation of intestinal flora on the expression of CXCs. A mouse model was used to verify the fact that intestinal bacterial dysregulation can affect the expression of CXCs. Results: In the three cancers, multiple datasets verified the fact that the mRNA expression of this family was significantly different; the mRNA levels of CXCL3, 8, 9, 10, 14, and 17 were significantly correlated with the prognosis of three cancers. CXCs were correlated with six types of immuno-infiltrating cells in three cancers. Immunohistochemistry of clinical samples confirmed that the expression of CXCL8 and 10 was higher in three cancer tissues. Animal experiments have shown that intestinal flora dysregulation can affect CXCL8 and 10 expressions. Conclusion: Our results further elucidate the function of CXCs in cancers and provide new insights into the prognosis and immune infiltration of breast, colon, and pancreatic cancers, and they suggest that intestinal flora may influence disease progression through CXCs.

Comprehensive Analysis of CD163 as a Prognostic Biomarker and Associated with Immune Infiltration in Glioblastoma Multiforme

Background

Glioblastoma multiforme (GBM) is the most common and aggressive primary malignancy in adults with high aggression. The prognosis of GBM patients is poor. There is a critical need for novel biomarkers for the prognosis and therapy of GBM.

Methods

Differentially expressed genes (DEGs) in GBM were screened using TCGA cohort. Univariate and multivariate Cox regression analyses were performed on DEGs to identify the optimal prognosis-related genes. qRT-PCR was performed to verify the result.

Results

A total of 5216 DEGs, including 2785 upregulated and 2458 downregulated genes, were obtained. Enrichment analysis revealed that these DEGs were mainly involved in the p53 signaling pathway and cell cycle, immune response, and MAPK signaling pathways. Moreover, the top 50 DEGs were associated with drug resistance or drug sensitivity. Prognosis analysis revealed that GBM patients with a high expression of CD163 and CHI3L2 had a poor overall survival, prognosis-free survival, and disease-specific survival. The univariate and multivariate analyses revealed that CD163 and age were independent factors affecting the prognosis of GBM patients. A validation study revealed that CD163 was upregulated in GBM tissues and associated with poor overall survival. Moreover, further analysis revealed that CD163 showed significant correlation with immune cells, immune biomarkers, chemokines, and chemokine receptors. We also identified several CD163-associated kinase, miRNA, and transcription factor targets in GBM, including LCK, miR-483, and ELF1.

Conclusions

In conclusion, our study suggested CD163 as a prognostic biomarker and associated it with immune infiltration in GBM.

Analysis of therapeutic targets and prognostic biomarkers of CXC chemokines in cervical cancer microenvironment

Background

The tumor microenvironment (TME) has received an increasing amount of attention. CXC chemokines can regulate immune cell transport and tumor cell activity to exert anti-tumor immunity. However, studies on the expression and prognosis of CXC chemokines in cervical cancer (CC) are more limited.

Methods

The study investigated the role of CXC chemokines in TME of CC by using public databases. Moreover, quantitative real-time PCR (qRT-PCR) and immunohistochemistry (IHC) of CXC chemokines were performed to further verify.

Results

The transcriptional levels of CXCL1/3/5/6/8/9/10/11/13/16/17 in CC tissues were significantly elevated while the transcriptional levels of CXCL12/14 were significantly reduced. We reached a consistent conclusion that the expression of CXCL9/10/11/13 was verified by quantitative real-time PCR and immunohistochemistry. Moreover, CC patients with low transcriptional levels of CXCL1/2/3/4/5/8 were significantly associated with longer overall survival (OS). The CCL family was related to CXC chemokines neighboring alteration. RELA, NFKB1, LCK and PAK2 were the key transcription factors and kinase targets of CXC chemokines, respectively. We also found there were significant correlations between the expression of CXCL9/10/11 and the infiltration of immune cells (CD8+ T cell, CD4+ T cell, neutrophils and dendritic cells).

Conclusions

In brief, we conducted a comprehensive analysis of CXC chemokines via clinical data and some online public databases. Our results may provide a new idea for the selection of immunotherapeutic targets and prognostic biomarkers for cervical cancer.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02101-9.