Comprehensive Analysis of the Expression and Prognosis for GBPs in Head and neck squamous cell carcinoma

Protein structural context of cancer mutations reveals molecular mechanisms and candidate driver genes

Advances in protein structure determination and modeling allow us to study the structural context of human genetic variants on an unprecedented scale. Here, we analyze millions of cancer-associated missense mutations based on their structural locations and predicted perturbative effects. By considering the collective properties of mutations at the level of individual proteins, we identify distinct patterns associated with tumor suppressors and oncogenes. Tumor suppressors are enriched in structurally damaging mutations, consistent with loss-of-function mechanisms, while oncogene mutations tend to be structurally mild, reflecting selection for gain-of-function driver mutations and against loss-of-function mutations. Although oncogenes are difficult to distinguish from genes with no role in cancer using only structural damage, we find that the three-dimensional clustering of mutations is highly predictive. These observations allow us to identify candidate driver genes and speculate about their molecular roles, which we expect will have general utility in the analysis of cancer sequencing data.

Guanylate binding protein 5 is an immune-related biomarker of oral squamous cell carcinoma: A retrospective prognostic study with bioinformatic analysis

BACKGROUND

Cancer utilizes immunosuppressive mechanisms to create a tumor microenvironment favorable for its progression. The purpose of this study is to histologically characterize the immunological properties of the tumor microenvironment of oral squamous cell carcinoma (OSCC) and identify key molecules involved in the immunological microenvironment and patient prognosis.

METHODS

First, overlapping differentially expressed genes (DEGs) were screened from OSCC transcriptome data in public databases. Correlation analysis of DEGs with known immune-related genes identified genes involved in the immune microenvironment of OSCC. Next, stromal patterns of tumor were classified and immunohistochemical staining was performed for immune cell markers (CD3, CD4, Foxp3, CD8, CD20, CD68, and CD163), programmed death-ligand 1 (PD-L1), and guanylate binding protein 5 (GBP5) in resected specimens obtained from 110 patients with OSCC who underwent resection. Correlations between each factor and their prognostic impact were analyzed.

RESULTS

Among the novel OSCC-specific immune-related genes screened (including ADAMDEC1, CXCL9, CXCL13, DPT, GBP5, IDO1, and PLA2G7), GBP5 was selected as the target gene. Histopathologic analysis showed that multiple T-cell subsets and CD20-positive cells were less common in the advanced stages, whereas CD163-positive cells were more common in advanced stages. The immature type in the stromal pattern category was associated with less immune cell infiltration, lower expression of PD-L1 in immune cells, lower expression of GBP5 in the stroma, and shorter overall survival and recurrence-free survival. Expression of GBP5 in the tumor and stroma correlated with immune cell infiltration of tumors and PD-L1 expression in tumor and immune cells. Patients with low tumor GBP5 expression and high stromal expression had significantly longer overall survival and recurrence-free survival.

CONCLUSIONS

The stromal pattern category may reflect both invasive and immunomodulatory potentials of cancer-associated fibroblasts in OSCC. GBP5 has been suggested as a potential biomarker to predict the prognosis and therapeutic efficacy of immune checkpoint inhibitors.

Identification of MUC1-C as a Target for Suppressing Progression of Head and Neck Squamous Cell Carcinomas

The MUC1-C protein is aberrantly expressed in adenocarcinomas of epithelial barrier tissues and contributes to their progression. Less is known about involvement of MUC1-C in the pathogenesis of squamous cell carcinomas (SCC). Here, we report that the MUC1 gene is upregulated in advanced head and neck SCCs (HNSCC). Studies of HNSCC cell lines demonstrate that the MUC1-C subunit regulates expression of (i) RIG-I and MDA5 pattern recognition receptors, (ii) STAT1 and IFN regulatory factors, and (iii) downstream IFN-stimulated genes. MUC1-C integrates chronic activation of the STAT1 inflammatory pathway with induction of the ∆Np63 and SOX2 genes that are aberrantly expressed in HNSCCs. In extending those dependencies, we demonstrate that MUC1-C is necessary for NOTCH3 expression, self-renewal capacity, and tumorigenicity. The findings that MUC1 associates with ∆Np63, SOX2 and NOTCH3 expression by single-cell RNA sequencing analysis further indicate that MUC1-C drives the HNSCC stem cell state and is a target for suppressing HNSCC progression.

SIGNIFICANCE

This work reports a previously unrecognized role for MUC1-C in driving STAT1-mediated chronic inflammation with the progression of HNSCC and identifies MUC1-C as a druggable target for advanced HNSCC treatment.

Up-regulation of extracellular-matrix and inflammation related genes in oral squamous cell carcinoma

OBJECTIVE

Oral squamous cell carcinoma (OSCC) is the most prevalent malignancy with late-presentation, site-specific heterogeneity, and high-propensity for recurrence/metastasis that has shown rise in mortality. Lately, research emphasize on dynamic interactions between tumor-cells and extracellular-matrix components within tumor-microenvironment that promote tissue integrity loss and carcinogenesis. Therefore, OSCC clinical-management is still challenging.

DESIGN

Present study validated clinical utility of a 13 gene-panel in two chief sub-sites of OSCC: Buccal mucosa squamous cell carcinoma (BMSCC) (N = 50) and Tongue squamous cell carcinoma (TSCC) (N = 52) using qRT-PCR. Principal component analysis and binary logistic regression analysis were applied to acquire definite multi gene models. Protein expression analysis was employed using the Human Protein Atlas, UALCAN and TIMER 2.0 databases to explore potential correlation between immune cells and gene-panels.

RESULTS

Significant up-regulation of CXCL8, CXCL10, FN1, GBP1, IFIT3, ISG15, MMP1, MMP3, MMP10, PLAU, SERPINE1 and SPP1 except OASL was observed in OSCC tissue in comparison of absolute normal controls. Although, this gene-panel could potentially discriminate OSCC tissues from absolute normal controls as solitarily diagnostic and/or predictive biomarkers, models generated also showed substantial discriminating efficacy. Eight-genes were found to be significantly associated with poor-prognosis on clinico-pathological association. Protein-expression confirmed overexpression of gene-panel and added advantage of being secretory-protein. Importantly, up-regulated genes in our study showed significant relation with immune-cells infiltration suggesting their contribution in immune-escape.

CONCLUSION

Thus, we propose that the 13 gene-panel could pave the way to effective and personalized clinical-management of OSCC in terms of diagnostic and prognostic measures and thereby as therapeutic targets.

Smoking-associated gene expression alterations in nasal epithelium reveal immune impairment linked to lung cancer risk

BACKGROUND

Lung cancer is the leading cause of cancer-related death in the world. In contrast to many other cancers, a direct connection to modifiable lifestyle risk in the form of tobacco smoke has long been established. More than 50% of all smoking-related lung cancers occur in former smokers, 40% of which occur more than 15 years after smoking cessation. Despite extensive research, the molecular processes for persistent lung cancer risk remain unclear. We thus set out to examine whether risk stratification in the clinic and in the general population can be improved upon by the addition of genetic data and to explore the mechanisms of the persisting risk in former smokers.

METHODS

We analysed transcriptomic data from accessible airway tissues of 487 subjects, including healthy volunteers and clinic patients of different smoking statuses. We developed a computational model to assess smoking-associated gene expression changes and their reversibility after smoking is stopped, comparing healthy subjects to clinic patients with and without lung cancer.

RESULTS

We find persistent smoking-associated immune alterations to be a hallmark of the clinic patients. Integrating previous GWAS data using a transcriptional network approach, we demonstrate that the same immune- and interferon-related pathways are strongly enriched for genes linked to known genetic risk factors, demonstrating a causal relationship between immune alteration and lung cancer risk. Finally, we used accessible airway transcriptomic data to derive a non-invasive lung cancer risk classifier.

CONCLUSIONS

Our results provide initial evidence for germline-mediated personalized smoke injury response and risk in the general population, with potential implications for managing long-term lung cancer incidence and mortality.

Guanylate binding protein 5 triggers NF-κB activation to foster radioresistance, metastatic progression and PD-L1 expression in oral squamous cell carcinoma

Radioresistance and metastasis are critical issues in managing oral squamous cell carcinoma (OSCC). Although immune checkpoint inhibitors (ICIs) has been recommended to treat OSCC, lacking useful biomarkers limited their anti-cancer effectiveness. We found that guanylate binding protein 5 (GBP5) is upregulated in primary tumors and associates with radioresistance in OSCC. GBP5 expression causally associated with cellular radioresistance and migration ability in the OSCC cell variants. GBP5 upregulation was examined to be correlated with NF-κB activation and programmed cell death-ligand 1 (PD-L1) elevation in OSCC samples. GBP5 knockdown was mitigated, but overexpression enhanced, NF-κB activity and PD-L1 expression in the OSCC cells. NF-κB inhibition by SN50 dramatically suppressed the GBP5-forested irradiation resistance, cellular migration ability and PD-L1 expression in OSCC cells. Importantly, GBP5 upregulation predicted a favorable outcome in cancer patients received ICI treatment. Our findings provide GBP5 as a useful biomarker to predict the anti-OSCC effectiveness of irradiation and ICIs.

Function and mechanism of GBP1 in the development and progression of cervical cancer

Guanylate binding protein 1 (GBP1) is the most concerned member of the GBP family, which has a series of effects such as anti-infection and anti-angiogenesis. Its role in malignant tumors including cervical cancer is still controversial. We aim to explore the effects of GBP1 on cervical cancer through bioinformatics and related experiments. In this study, we first found that GBP1 was generally expressed in cervical cancer in various online databases and was closely related to immune invasion. Secondly, we used multicolor immunofluorescence technology to verify the expression of GBP1 in cervical cancer tissues and its relationship with immune invasion, and explored its relationship with the prognosis of patients with cervical cancer. Knockdown and overexpression assays of GBP1 in vitro were used to prove GBP1 as a potential oncogene of cervical cancer, and its carcinogenicity was verified by in vivo experiment. In order to explore the potential mechanism of GBP1 in promoting cancer, RNA-seq was performed on GBP1 overexpression and knockdown expression cell lines, and GBP1 knockdown and overexpression were found to be associated with many RNA alternative splicing events, suggesting that GBP1 maybe a RNA binding protein (RBP) which affect the biological characteristics of cervical cancer cells through the alternative splicing pathway. However, the later RNA binding protein immunoprecipitation (RIP) assay proved that GBP1 was not a direct alternative splicing factor, while the co-immunoprecipitation (CoIP)-mass spectroscopy (MS) assay combined with protein protein interaction (PPI) analysis proved that 8 alternative splicing factors including Heterogeneous Nuclear Ribonucleoprotein K (HNRNPK) were interacting proteins of GBP1. Combined with the existing reports and the results of RNA-seq alternative splicing analysis, it is speculated that GBP1 may regulate the alternative splicing of CD44 protein by binding to interacting protein-HNRNPK, and thus play a role in promoting cancer in cervical cancer.

Expression and role of FKBPL in lung adenocarcinoma

Dysregulated expression of FK506-binding protein like (FKBPL) has been demonstrated to play crucial roles in tumour development. However, the role of FKBPL in lung adenocarcinoma (ADC) remains unclear. Using immunohistochemical staining, we showed that FKBPL expression was significantly lower in lung ADC than the normal tissues (P < 0.0001). Patients with well or moderately differentiated tumours have higher FKBPL expression compared with patients with poor differentiated tumours (P = 0.037). However, no significant associations were found between FKBPL expression and other clinicopathological variables (P > 0.05 for all). Cox univariate analysis showed that high FKBPL expression was correlated with prolonged overall survival (OS) (P = 0.010). Kaplan-Meier analysis further confirmed that the FKBPL-low group showed a significantly shorter OS than the FKBPL-high group (P = 0.0081). FKBPL expression was not shown as an independent prognostic factor for OS in the multivariate analysis (P = 0.063). Moreover, our study demonstrated that FKBPL could suppress the proliferation of lung ADC cells by delaying cell cycle G1/S phase transition. In addition, FKBPL resulted in increased apoptosis in lung ADC cells. Using the Human Apoptosis Array Kit, we observed that overexpression of FKBPL in lung ADC A549 cells significantly decreased the anti-apoptotic proteins, including heat shock protein 32 (HSP32), heat shock protein 27 (HSP27), and paraoxonase-2 (PON2). FKBPL depletion significantly attenuated the pro-apoptotic protein, phospho-p53 (S46), in lung ADC H1975 cells. These new findings provide an experimental basis for further theoretical investigation of lung ADC.

A Genomic Analysis of Esophageal Squamous Cell Carcinoma in Eastern Africa

BACKGROUND

Esophageal squamous cell carcinoma (ESCC) comprises 90% of all esophageal cancer cases globally and is the most common histology in low-resource settings. Eastern Africa has a disproportionately high incidence of ESCC.

METHODS

We describe the genomic profiles of 61 ESCC cases from Tanzania and compare them to profiles from an existing cohort of ESCC cases from Malawi. We also provide a comparison to ESCC tumors in The Cancer Genome Atlas (TCGA).

RESULTS

We observed substantial transcriptional overlap with other squamous histologies via comparison with TCGA PanCan dataset. DNA analysis revealed known mutational patterns, both genome-wide as well as in genes known to be commonly mutated in ESCC. TP53 mutations were the most common somatic mutation in tumors from both Tanzania and Malawi but were detected at lower frequencies than previously reported in ESCC cases from other settings. In a combined analysis, two unique transcriptional clusters were identified: a proliferative/epithelial cluster and an invasive/migrative/mesenchymal cluster. Mutational signature analysis of the Tanzanian cohort revealed common signatures associated with aging and cytidine deaminase activity (APOBEC) and an absence of signature 29, which was previously reported in the Malawi cohort.

CONCLUSIONS

This study defines the molecular characteristics of ESCC in Tanzania, and enriches the Eastern African dataset, with findings of overall similarities but also some heterogeneity across two unique sites.

IMPACT

Despite a high burden of ESCC in Eastern Africa, investigations into the genomics in this region are nascent. This represents the largest comprehensive genomic analysis ESCC from sub-Saharan Africa to date.

Guanylate-binding proteins signature predicts favorable prognosis, immune-hot microenvironment, and immunotherapy response in hepatocellular carcinoma

BACKGROUND

The role of guanylate-binding proteins (GBPs) in various cancers has been elucidated recently. However, our knowledge of the clinical relevance and biological characteristics of GBPs in hepatocellular carcinoma (HCC) remains limited.

METHODS

A total of 955 HCC patients were enrolled from five independent public HCC cohorts. The role of GBP molecules in HCC was preliminarily investigated, and a GBP family signature, termed GBPs-score, was constructed by principal component analysis to combine the GBP molecule values. We revealed the effects of GBP genes and GBPs-score in HCC via well-established bioinformatics methods and validated GBP1-5 experimentally in a tissue microarray (TMA) cohort.

RESULTS

GBPs molecules were closely associated with the prognosis of patients with HCC, and a high GBPs-score highly inferred a favorable survival outcome. We also revealed high GBPs-score was related to anti-tumor immunity, the immune-hot tumor microenvironment (TME), and immunotherapy response. Among the GBPs members, GBP1-5 rather than GBP6/7 may be dominant in these fields. The TMA analysis based on immunohistochemistry showed positive correlations between GBP1-5 and the immune-hot TME with abundant infiltration of CD8+ T cells in HCC.

CONCLUSIONS

Our integrative study revealed the genetic and immunologic characterizations of GBPs in HCC and highlighted their potential values as promising biomarkers for prognosis and immunotherapy.

Construction of a four-mRNA prognostic signature with its ceRNA network in CESC

Cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC) tumorigenesis involves a combination of multiple genetic alteration processes. Constructing a survival-associated competing endogenous RNA (ceRNA) network and a multi-mRNA-based prognostic signature model can help us better understand the complexity and genetic characteristics of CESC. In this study, the RNA-seq data and clinical information of CESC patients were downloaded from The Cancer Genome Atlas. Differentially expressed mRNAs, lncRNAs and miRNAs were identified with the edgeR R package. A four-mRNA prognostic signature was developed by multivariate Cox regression analysis. Kaplan–Meier survival with the log-rank tests was performed to assess survival rates. The relationships between overall survival (OS) and clinical parameters were evaluated by Cox regression analysis. A survival-associated ceRNA network was constructed with the multiMiR package and miRcode database. Kyoto encyclopedia of genes and genomes (KEGG) analysis and gene ontology analyses were used to identify the functional role of the ceRNA network in the prognosis of CESC. A total of 298 differentially expressed mRNAs, 8 miRNAs, and 29 lncRNAs were significantly associated with the prognosis of CESC. A prognostic signature model based on 4 mRNAs (OPN3, DAAM2, HENMT1, and CAVIN3) was developed, and the prognostic ability of this signature was indicated by the AUC of 0.726. Patients in the high-risk group exhibited significantly worse OS. The KEGG pathways, TGF-β and Cell adhesion molecules, were significantly enriched. In this study, a CESC-associated ceRNA network was constructed, and a multi-mRNA-based prognostic model for CESC was developed based on the ceRNA network, providing a new perspective for cancer pathogenesis research.

Survival-related genes are diversified across cancers but generally enriched in cancer hallmark pathways

Background

Pan-cancer studies have disclosed many commonalities and differences in mutations, copy number variations, and gene expression alterations among cancers. Some of these features are significantly associated with clinical outcomes, and many prognosis-predictive biomarkers or biosignatures have been proposed for specific cancer types. Here, we systematically explored the biological functions and the distribution of survival-related genes (SRGs) across cancers.

Results

We carried out two different statistical survival models on the mRNA expression profiles in 33 cancer types from TCGA. We identified SRGs in each cancer type based on the Cox proportional hazards model and the log-rank test. We found a large difference in the number of SRGs among different cancer types, and most of the identified SRGs were specific to a particular cancer type. While these SRGs were unique to each cancer type, they were found mostly enriched in cancer hallmark pathways, e.g., cell proliferation, cell differentiation, DNA metabolism, and RNA metabolism. We also analyzed the association between cancer driver genes and SRGs and did not find significant over-representation amongst most cancers.

Conclusions

In summary, our work identified all the SRGs for 33 cancer types from TCGA. In addition, the pan-cancer analysis revealed the similarities and the differences in the biological functions of SRGs across cancers. Given the potential of SRGs in clinical utility, our results can serve as a resource for basic research and biotech applications.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08581-x.

Guanylate-Binding Protein 1 as a Potential Predictor of Immunotherapy: A Pan-Cancer Analysis

Background: Mainstream application of cancer immunotherapy is hampered by the low response rate of most cancer patients. A novel immunotherapeutic target or a biomarker predicting response to immunotherapy needs to be developed. Guanylate-binding protein 1 (GBP1) is an interferon (IFN)-inducible guanosine triphosphatases (GTPases) involving inflammation and infection. However, the immunological effects of GBP1 in pan-cancer patients are still obscure.

Methods: Using large-scale public data, we delineated the landscape of GBP1 across 33 cancer types. The correlation between GBP1 expression or mutation and immune cell infiltration was estimated by ESTIMATE, TIMER, xCell, and quanTIseq algorithms. GBP1-related genes and proteins were subjected to function enrichment analysis. Clustering analysis explored the relationship between GBP1 expression and anti-tumor immune phenotypes. We assessed the patient’s response to immunotherapy using the tumor immune dysfunction and exclusion (TIDE) score and immunophenoscore (IPS). Furthermore, we validated the predictive power of GBP1 expression in four independent immunotherapy cohorts.

Results: GBP1 was differentially expressed in tumors and normal tissues in multiple cancer types. Distinct correlations existed between GBP1 expression and prognosis in cancer patients. GBP1 expression and mutation were positively associated with immune cell infiltration. Function enrichment analysis showed that GBP1-related genes were enriched in immune-related pathways. Positive correlations were also observed between GBP1 expression and the expression of immune checkpoints, as well as tumor mutation burden (TMB). Pan-cancer patients with higher GBP1 expression were more inclined to display “hot” anti-tumor immune phenotypes and had lower TIDE scores and higher immunophenoscore, suggesting that these patients had better responses to immunotherapy. Patients with higher GBP1 expression exhibited improved overall survival and clinical benefits in immunotherapy cohorts, including the Gide et al. cohort [area under the curve (AUC): 0.813], the IMvigor210 cohort (AUC: 0.607), the Lauss et al. cohort (AUC: 0.740), and the Kim et al. cohort (AUC: 0.793).

Conclusion: This study provides comprehensive insights into the role of GBP1 in a pan-cancer manner. We identify GBP1 expression as a predictive biomarker for immunotherapy, potentially enabling more precise and personalized immunotherapeutic strategies in the future.

Multi-Omics Analysis of the Effects of Smoking on Human Tumors

Comprehensive studies on cancer patients with different smoking histories, including non-smokers, former smokers, and current smokers, remain elusive. Therefore, we conducted a multi-omics analysis to explore the effect of smoking history on cancer patients. Patients with smoking history were screened from The Cancer Genome Atlas database, and their multi-omics data and clinical information were downloaded. A total of 2,317 patients were included in this study, whereby current smokers presented the worst prognosis, followed by former smokers, while non-smokers showed the best prognosis. More importantly, smoking history was an independent prognosis factor. Patients with different smoking histories exhibited different immune content, and former smokers had the highest immune cells and tumor immune microenvironment. Smokers are under a higher incidence of genomic instability that can be reversed following smoking cessation in some changes. We also noted that smoking reduced the sensitivity of patients to chemotherapeutic drugs, whereas smoking cessation can reverse the situation. Competing endogenous RNA network revealed that mir-193b-3p, mir-301b, mir-205-5p, mir-132-3p, mir-212-3p, mir-1271-5p, and mir-137 may contribute significantly in tobacco-mediated tumor formation. We identified 11 methylation driver genes (including EIF5A2, GBP6, HGD, HS6ST1, ITGA5, NR2F2, PLS1, PPP1R18, PTHLH, SLC6A15, and YEATS2), and methylation modifications of some of these genes have not been reported to be associated with tumors. We constructed a 46-gene model that predicted overall survival with good predictive power. We next drew nomograms of each cancer type. Interestingly, calibration diagrams and concordance indexes are verified that the nomograms were highly accurate for the prognosis of patients. Meanwhile, we found that the 46-gene model has good applicability to the overall survival as well as to disease-specific survival and progression-free intervals. The results of this research provide new and valuable insights for the diagnosis, treatment, and follow-up of cancer patients with different smoking histories.

Clinical Significance and the Role of Guanylate-Binding Protein 5 in Oral Squamous Cell Carcinoma

Guanylate binding protein 5 (GBP5) is the interferon (IFN)-inducible subfamily of guanosine triphosphatases (GTPases) and is involved in pathogen defense. However, the role played by GBP5 in cancer development, especially in oral squamous cell carcinoma (OSCC), is still unknown. Herein, next-generation sequencing analysis showed that the gene expression levels of GBP5 were significantly higher in OSCC tissues compared with those found in corresponding tumor adjacent normal tissues (CTAN) from two pairs of OSCC patients. Higher gene expression levels of GBP5 were also found in tumor tissues of 23 buccal mucosal squamous cell carcinoma (BMSCC)/14 tongue squamous cell carcinoma (TSCC) patients and 30 oral cancer patients from The Cancer Genome Atlas (TCGA) database compared with those in CTAN tissues. Immunohistochemical results showed that protein expression levels of GBP5 were also higher in the tumor tissues of 353 OSCC patients including 117 BMSCC, 187 TSCC, and 49 lip squamous cell carcinoma patients. Moreover, TCGA database analysis indicated that high gene expression levels of GBP5 were associated with poor overall survival in oral cancer patients with moderate/poor cell differentiation, and associated with poor disease-free survival in oral cancer patients with moderate/poor cell differentiation and lymph node metastasis. Furthermore, GBP5-knockdowned cells exhibited decreased cell growth, arrest at G1 phase, and decreased invasion/migration. The gene expression of markers for epithelial-mesenchymal transition and cancer stemness was also reduced in GBP5-silenced oral cancer cells. Taken together, GBP5 might be a potential biomarker and therapeutic target for OSCC patients, especially for those with poor cell differentiation and lymph node metastasis.

Comprehensive analysis of the expression and prognosis for TNFAIPs in head and neck cancer

Head and neck cancer (HNC) tumorigenesis involves a combination of multiple genetic alteration processes. Tumour necrosis factor-alpha-induced proteins (TNFAIPs) are involved in tumour development and progression, but few studies have been conducted on these factors in HNC. We aimed to analyse TNFAIPs and assess their potential as prognostic biomarkers and therapeutic targets using the Oncomine, UALCAN, Human Protein Atlas, LinkedOmics, cBioPortal, GeneMANIA, Enrichr, and Tumor IMmune Estimation Resource databases. We found that the transcript levels of TNFAIP1, TNFAIP3, EFNA1, TNFAIP6 and TNFAIP8 were increased, while those of TNFAIP8L3 and STEAP4 were reduced in HNC tissues versus normal tissues. The EFNA1, TNFAIP8 and TNFAIP8L3 expression levels were significantly correlated with the pathological stage. In HNC patients, high PTX3 and TNFAIP6 transcript levels were significantly associated with shorter overall survival (OS). Moreover, genetic alterations in TNFAIP1, TNFAIP6, and STEAP4 resulted in poorer disease-free survival, progression-free survival, and OS, respectively. TNFAIPs may mediate HNC tumorigenesis by regulating PI3K-Akt, Ras and other signalling pathways. TNFAIPs are also closely correlated with the infiltration of immune cells, including B cells, CD8⁺ T cells, CD4⁺ T cells, etc. The data above indicate that TNFAIPs may be potential biomarkers and therapeutic targets for HNC.

Identification of Prognostic Biomarkers and Correlation With Immune Infiltrates in Hepatocellular Carcinoma Based on a Competing Endogenous RNA Network

Background

Hepatocellular carcinoma (HCC) is one of the most common malignant tumors worldwide. Recently, competing endogenous RNAs (ceRNA) have revealed a significant role in the progression of HCC. Herein, we aimed to construct a ceRNA network to identify potential biomarkers and illustrate its correlation with immune infiltration in HCC.

Methods

RNA sequencing data and clinical traits of HCC patients were downloaded from TCGA. The limma R package was used to identify differentially expressed (DE) RNAs. The predicted prognostic model was established using univariate and multivariate Cox regression. A K-M curve, TISIDB and GEPIA website were utilized for survival analysis. Functional annotation was determined using Enrichr and Reactome. Protein-to-protein network analysis was implemented using SRTNG and Cytoscape. Hub gene expression was validated by quantitative polymerase chain reaction, Oncomine and the Hunan Protein Atlas database. Immune infiltration was analyzed by TIMMER, and Drugbank was exploited to identify bioactive compounds.

Results

The predicted model that was established revealed significant efficacy with 3- and 5-years of the area under ROC at 0.804 and 0.744, respectively. Eleven DEmiRNAs were screened out by a K-M survival analysis. Then, we constructed a ceRNA network, including 56 DElncRNAs, 6 DEmiRNAs, and 28 DEmRNAs. The 28 DEmRNAs were enriched in cancer-related pathways, for example, the TNF signaling pathway. Moreover, six hub genes, CEP55, DEPDC1, KIF23, CLSPN, MYBL2, and RACGAP1, were all overexpressed in HCC tissues and independently correlated with survival rate. Furthermore, expression of hub genes was related to immune cell infiltration in HCC, including B cells, CD8⁺ T cells, CD4⁺ T cells, monocytes, macrophages, neutrophils, and dendritic cells.

Conclusion

The findings from this study demonstrate that CEP55, DEPDC1, KIF23, CLSPN, MYBL2, and RACGAP1 are closely associated with prognosis and immune infiltration, representing potential therapeutic targets or prognostic biomarkers in HCC.