Identification of DNA repair-related genes predicting pathogenesis and prognosis for liver cancer

Fatty acid binding protein 5 is a novel therapeutic target for hepatocellular carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) is an aggressive subtype of liver cancer and is one of the most common cancers with high mortality worldwide. Reprogrammed lipid metabolism plays crucial roles in HCC cancer cell survival, growth, and evolution. Emerging evidence suggests the importance of fatty acid binding proteins (FABPs) in contribution to cancer progression and metastasis; however, how these FABPs are dysregulated in cancer cells, especially in HCC, and the roles of FABPs in cancer progression have not been well defined.

AIM

To understand the genetic alterations and expression of FABPs and their associated cancer hallmarks and oncogenes in contributing to cancer malignancies.

METHODS

We used The Cancer Genome Atlas datasets of pan cancer and liver hepatocellular carcinoma (LIHC) as well as patient cohorts with other cancer types in this study. We investigated genetic alterations of FABPs in various cancer types. mRNA expression was used to determine if FABPs are abnormally expressed in tumor tissues compared to non-tumor controls and to investigate whether their expression correlates with patient clinical outcome, enriched cancer hallmarks and oncogenes previously reported for patients with HCC. We determined the protein levels of FABP5 and its correlated genes in two HCC cell lines and assessed the potential of FABP5 inhibition in treating HCC cells.

RESULTS

We discovered that a gene cluster including five FABP family members (FABP4, FABP5, FABP8, FABP9 and FABP12) is frequently co-amplified in cancer. Amplification, in fact, is the most common genetic alteration for FABPs, leading to overexpression of FABPs. FABP5 showed the greatest differential mRNA expression comparing tumor with non-tumor tissues. High FABP5 expression correlates well with worse patient outcomes (P < 0.05). FABP5 expression highly correlates with enrichment of G2M checkpoint (r = 0.33, P = 1.1e-10), TP53 signaling pathway (r = 0.22, P = 1.7e-5) and many genes in the gene sets such as CDK1 (r = 0.56, P = 0), CDK4 (r = 0.49, P = 0), and TP53 (r = 0.22, P = 1.6e-5). Furthermore, FABP5 also correlates well with two co-expressed oncogenes PLK1 and BIRC5 in pan cancer especially in LIHC patients (r = 0.58, P = 0; r = 0.58, P = 0; respectively). FABP5high Huh7 cells also expressed higher protein levels of p53, BIRC5, CDK1, CDK2, and CDK4 than FABP5low HepG2 cells. FABP5 inhibition more potently inhibited the tumor cell growth in Huh7 cells than in HepG2 cells.

CONCLUSION

We discovered that FABP5 gene is frequently amplified in cancer, especially in HCC, leading to its significant elevated expression in HCC. Its high expression correlates well with worse patient outcome, enriched cancer hallmarks and oncogenes in HCC. FABP5 inhibition impaired the cell viability of FABP5high Huh7 cells. All these support that FABP5 is a novel therapeutic target for treating FABP5high HCC.

Prognostic and immunological roles of heat shock protein A4 in lung adenocarcinoma

BACKGROUND

Heat shock protein A4 (HSPA4) belongs to molecular chaperone protein family which plays important roles within variable cellular activities, including cancer initiation and progression. However, the prognostic and immunological significance of HSPA4 in lung adenocarcinoma (LUAD) has not been revealed yet.

AIM

To explore the prognostic and immunological roles of HSPA4 to identify a novel prognostic biomarker and therapeutic target for LUAD.

METHODS

We assessed the prognostic and immunological significance of HSPA4 in LUAD using data from The Cancer Genome Atlas database. The association between HSPA4 expression and clinical-pathological features was assessed through Kruskal-Wallis and Wilcoxon signed-rank test. Univariate/multivariate Cox regression analyses and Kaplan-Meier curves were employed to evaluate prognostic factors, including HSPA4, in LUAD. Gene set enrichment analysis (GSEA) was conducted to identify the key signaling pathways associated with HSPA4. The correlation between HSPA4 expression and cancer immune infiltration was evaluated using single-sample gene set enrichment analysis (ssGSEA).

RESULTS

Overexpressing HSPA4 was significantly related to advanced pathologic TNM stage, advanced pathologic stage, progression disease status of primary therapy outcome and female subgroups with LUAD. In addition, increased HSPA4 expression was found to be related to worse disease-specific survival and overall survival. GSEA analysis indicated a significant correlation between HSPA4 and cell cycle regulation and immune response, particularly through diminishing the function of cytotoxicity cells and CD8 T cells. The ssGSEA algorithm showed a positive correlation between HSPA4 expression and infiltrating levels of Th2 cells, while a negative correlation was observed with cytotoxic cell infiltration levels.

CONCLUSION

Our findings indicate HSPA4 is related to prognosis and immune cell infiltrates and may act as a novel prognostic biomarker and therapeutic target for LUAD.

A pan-cancer analysis of pituitary tumor-transforming 3, pseudogene

BACKGROUND

Although evidence regarding pituitary tumor-transforming 3, pseudogene (PTTG3P) involvement in human cancers has been acquired via human and animal model-based molecular studies, there is a lack of pan-cancer analysis of this gene in human tumors.

METHODS

Tumor-causing effects of PTTG3P in 24 human tumors were explored using The Cancer Genome Atlas (TCGA) datasets from different bioinformatics databases and applying in silico tools such as The University of ALabama at Birmingham CANcer (UALCAN), Human Protein Atlas (HPA), Kaplan Meier (KM) plotter, cBioPortal, Search Tool for the Retrieval of Interacting Genes/Proteins (STRING), Cytoscape, Database for Annotation, Visualization, and Integrated Discovery (DAVID), Tumor IMmune Estimation Resource (TIMER), and Comparative Toxicogenomics Database (CTD). Then, via in vitro experiments, including RNA sequencing (RNA-seq) and targeted bisulfite sequencing (bisulfite-seq), expression and promoter methylation levels of PTTG3P were verified in cell lines.

RESULTS

The PTTG3P expression was overexpressed across 23 malignancies and its overexpression was further found significantly effecting the overall survival (OS) durations of the esophageal carcinoma (ESCA) and head and neck cancer (HNSC) patients. This important information helps us to understand that PTTG3P plays a significant role in the development and progression of ESCA and HNSC. As for PTTG3P functional mechanisms, this gene along with its other binding partners was significantly concentrated in "Oocyte meiosis", "Cell cycle", "Ubiquitin mediated proteolysis", and "Progesterone-mediated oocyte maturation". Moreover, ESCA and HNSC tissues having the higher expression of PTTG3P were found to have lower promoter methylation levels of PTTG3P and higher CD8+ T immune cells level. Additionally, PTTG3P expression-regulatory drugs were also explored in the current manuscript for designing appropriate treatment strategies for ESCA and HNSC with respect to PTTG3P expression.

CONCLUSION

Our pan-cancer based findings provided a comprehensive account of the oncogenic role and utilization of PTTG3P as a novel molecular biomarker of ESCA and HNSC.

Knockdown of RFC4 inhibits the cell proliferation of nasopharyngeal carcinoma in vitro and in vivo

Nasopharyngeal carcinoma (NPC) is a malignant tumor that mainly occurs in East and Southeast Asia. Although patients benefit from the main NPC treatments (e.g., radiotherapy and concurrent chemotherapy), persistent and recurrent diseases still occur in some NPC patients. Therefore, investigating the pathogenesis of NPC is of great clinical significance. In the present study, replication factor c subunit 4 (RFC4) is a key potential target involved in NPC progression via bioinformatics analysis. Furthermore, the expression and mechanism of RFC4 in NPC were investigated in vitro and in vivo. Our results revealed that RFC4 was more elevated in NPC tumor tissues than in normal tissues. RFC4 knockdown induced G2/M cell cycle arrest and inhibited NPC cell proliferation in vitro and in vivo. Interestingly, HOXA10 was confirmed as a downstream target of RFC4, and the overexpression of HOXA10 attenuated the silencing of RFC4-induced cell proliferation, colony formation inhibition, and cell cycle arrest. For the first time, this study reveals that RFC4 is required for NPC cell proliferation and may play a pivotal role in NPC tumorigenesis.

Up-Frameshift Suppressor 3 as a prognostic biomarker and correlated with immune infiltrates: A pan-cancer analysis

BACKGROUND

The mRNA expression of protein Up-Frameshift Suppressor 3 Homolog B (UPF3B) differ in different tumors. However, the clinical relevance of UPF3B in cancer patients, such as with prognosis, tumor stage, and levels of tumor-infiltrating immune cells remain unclear.

METHODS

We performed bioinformatics analysis of UPF3B with The Cancer Genome Atlas (TCGA) database (https://xenabrowser.net) and TIMER2.0 (Tumor Immune Estimation Resource 2.0, http://timer.comp-genomics.org/). UPF3B expression in 33 cancers versus counterpart normal tissues was analyzed using TCGA pan-cancer data. The influence of UPF3B in long-term prognosis was evaluated using Kaplan-Meier method, and the associations between UPF3B transcription levels and immune-related gene expression, immune cell infiltration, tumor microenvironment (TME) score are analyzed by spearman correlation analysis. Enrichment analysis of UPF3B was conducted using the R package "clusterProfiler."

RESULTS

The transcriptional level of UPF3B was dysregulated in the human pan-cancer dataset. A significant correlation was found between the expression of UPF3B and the pathological stage of Esophageal Carcinoma (ESCA), Kidney Chromophobe (KIHC), Liver Hepatocellular Carcinoma (LIHC), and Skin Cutaneous Melanoma (SKCM). Multiple cancer types with high transcriptional levels of UPF3B were associated with a significantly worse prognosis. The functions of expressed UPF3B gene are primarily related to ubiquitin mediated proteolysis, cell cycle, and mRNA surveillance pathway. Our results also show that immune cells infiltration and immunosuppressive markers such as CTLA-4, PD-1 and PD-L1 significantly correlate with UPF3B expression.

CONCLUSIONS

In the present study, we synthetically explored the expression status and prognostic significance of UPF3B, and the relationship with clinic characters and immune microenvironment across cancers. Our results may provide novel insights for UPF3B as an immunotherapeutic target and valuable prognostic biomarker in various malignant tumor.

The pattern of expression and prognostic value of key regulators for m7G RNA methylation in hepatocellular carcinoma

N7-methylguanosine (m7G) modification on internal RNA positions plays a vital role in several biological processes. Recent research shows m7G modification is associated with multiple cancers. However, in hepatocellular carcinoma (HCC), its implications remain to be determined. In this place, we need to interrogate the mRNA patterns for 29 key regulators of m7G RNA modification and assess their prognostic value in HCC. Initial, the details from The Cancer Genome Atlas (TCGA) database concerning transcribed gene data and clinical information of HCC patients were inspected systematically. Second, according to the mRNA profiles of 29 m7G RNA methylation regulators, two clusters (named 1 and 2, respectively) were identified by consensus clustering. Furthermore, robust risk signature for seven m7G RNA modification regulators was constructed. Last, we used the Gene Expression Omnibus (GEO) dataset to validate the prognostic associations of the seven-gene risk signature. We figured out that 24/29 key regulators of m7G RNA modification varied remarkably in their grades of expression between the HCC and the adjacent tumor control tissues. Cluster one compared with cluster two had a substandard prognosis and was also positively correlated with T classification (T), pathological stage, and vital status (fustat) significantly. Consensus clustering results suggested the expression pattern of m7G RNA modification regulators was correlated with the malignancy of HCC strongly. In addition, cluster one was extensively enriched in metabolic-related pathways. Seven optimal genes (METTL1, WDR4, NSUN2, EIF4E, EIF4E2, NCBP1, and NCBP2) were selected to establish the risk model for HCC. Indicating by further analyses and validation, the prognostic model has fine anticipating command and this probability signature might be a self supporting presage factor for HCC. Finally, a new prognostic nomogram based on age, gender, pathological stage, histological grade, and prospects were established to forecast the prognosis of HCC patients accurately. In essence, we detected association of HCC severity and expression levels of m7G RNA modification regulators, and developed a risk score model for predicting prognosis of HCC patients’ progression.

DNA Damage Repair Classifier Defines Distinct Groups in Hepatocellular Carcinoma

Simple Summary

DNA repair pathways have been implicated in hepatocellular carcinoma outcomes. We found that hepatocellular carcinomas (HCC) could be separated into two groups (high and low) based on the overall expression of genes involved in DNA repair. Among the low repair group, there were three subgroups, one of which shared features of the high repair group. Given the important role of liver in metabolism and detoxification and its regenerative capacity, proliferation and DNA damage responses are critical in subdividing major biological categories of liver tumors. High repair samples showed more proliferative and regenerative signatures and had poorer outcomes versus the low repair that were more associated with the genes involved in normal liver biology. These biological groups suggest that dysregulation in endogenous liver processes promotes a pro-tumorigenic microenvironment that may facilitate tumor progression or identify tumors that require more substantial clinical intervention.

Abstract

DNA repair pathways have been associated with variability in hepatocellular carcinoma (HCC) clinical outcomes, but the mechanism through which DNA repair varies as a function of liver regeneration and other HCC characteristics is poorly understood. We curated a panel of 199 genes representing 15 DNA repair pathways to identify DNA repair expression classes and evaluate their associations with liver features and clinicopathologic variables in The Cancer Genome Atlas (TCGA) HCC study. We identified two groups in HCC, defined by low or high expression across all DNA repair pathways. The low-repair group had lower grade and retained the expression of classical liver markers, whereas the high-repair group had more clinically aggressive features, increased p53 mutant-like gene expression, and high liver regenerative gene expression. These pronounced features overshadowed the variation in the low-repair subset, but when considered separately, the low-repair samples included three subgroups: L1, L2, and L3. L3 had high DNA repair expression with worse progression-free (HR 1.24, 95% CI 0.81–1.91) and overall (HR 1.63, 95% CI 0.98–2.71) survival. High-repair outcomes were also significantly worse compared with the L1 and L2 groups. HCCs vary in DNA repair expression, and a subset of tumors with high regeneration profoundly disrupts liver biology and poor prognosis.

A Panel of E2F Target Gene Signature Predicting the Prognosis of Hepatocellular Carcinoma

Hepatocellular carcinoma is one of the most malignant tumors, and the therapeutic effects of traditional treatments are poor. It is urgent to explore and identify new biomarkers and therapeutic targets to develop novel treatments which are individualized and effective. Three hallmarks, including E2F targets, G2M checkpoint and DNA repair, were collected by GSEA analysis. The panel of E2F-related gene signature consisted of five genes: HN1, KIF4A, CDCA3, CDCA8 and SSRP1. They had various mutation rates ranging from 0.8 to 5% in hepatocellular carcinoma, and patients with gene mutation had poorer prognosis. Among these genes, HN1 has the greatest mutation rate, and SSRP1 has the greatest impact on the model with a B (COX) value of 0.8842. Patients with higher expression of these genes had poorer prognosis. Kaplan-Meier curves in stratified survival analysis confirmed that patients with high risk scores had poor prognosis (p < 0.05). The results of univariate and multivariate COX survival analysis showed that risk score was closely related to the overall survival of patients with hepatocellular carcinoma. For clinical validation, we found that all the genes in the model were upregulated in hepatocellular carcinoma tissues compared to normal liver tissues, which was consistent with the previous results we obtained. Our study demonstrated that a panel of E2F target genes signature including five genes could predict the prognosis of hepatocellular carcinoma. This panel gene signature can facilitate the development of individualized and effective treatment for hepatocellular carcinoma.

CircSCAP interacts with SF3A3 to inhibit the malignance of non-small cell lung cancer by activating p53 signaling

BACKGROUND

Circular RNA (circRNA) has been recently identified as a critical regulator during carcinogenesis. However, the biological function and potential underlying mechanisms of circRNAs in lung cancer remain to be further elucidated.

METHODS

Here, we first evaluated the differentially expressed circRNAs between tumor and the matched adjacent nontumor tissues (3 pairs) of lung cancer patients via circRNA microarray. The expression of top five dysregulated circRNAs were tested in lung cancer cell lines and the circSCAP with concordant alteration in microarray data and cell lines was selected for further investigation. Then we validated the expression level of circSCAP in tumor and corresponding adjacent tissues (161 pairs) from a lung cancer cohort by RT-PCR analysis followed by correlation and prognosis analysis between circSCAP and clinical characteristics. Non-small cell lung cancer (NSCLC) accounts for the majority of lung cancer diagnosis (about 80% in the cohort used in this study). Therefore, we focused the role of circSCAP in NSCLC in the present study. In vitro and in vivo assays were performed to study the biological function of circSCAP in NSCLC. Biotin-labeled RNA pulldown and RNA immunoprecipitation (RIP) assays were carried out to identify the proteins directly interacting with circSCAP. The molecular mechanism of circSCAP-driven tumor suppression was demonstrated by immunoblotting, immunoprecipitation and luciferase reporter assays. In vitro and in vivo rescue experiments were conducted to verify the role of the circSCAP/SF3A3/p53 signaling axis in NSCLC.

RESULTS

We screened the expression profiles of human circRNAs in lung cancer tissues and found that hsa_circ_0065214 (termed as circSCAP) was significantly decreased. Kaplan-Meier analysis showed that patients with low level of circSCAP had a significantly poor prognosis. Gain- and loss-of-function experiments suggested that circSCAP played an important role in NSCLC cell proliferation, cell migration and apoptosis. Mechanistically, circSCAP directly binds to the SF3A3 protein, facilitating the reduction of SF3A3 by promoting its ubiquitin-proteasome-mediated degradation, which enhances the expression of MDM4-S to finally activate its downstream p53 signaling.

CONCLUSION

These findings illustrate a novel circSCAP/SF3A3/p53 signaling axis involved in suppressing the malignance of NSCLC and provide a promising target for NSCLC prognosis prediction and treatment.

Identification of DNA Repair-Related Genes Predicting Clinical Outcome for Thyroid Cancer

Recent studies have demonstrated the utility and superiority of DNA repair-related genes as novel biomarkers for cancer diagnosis, prognosis, and therapy. Here, we aimed to screen the potential survival-related DNA repair-related genes in thyroid cancer (TC). TCGA datasets were utilized to analyze the differentially expressed DNA repair-related genes between TC and nontumor tissues. The K–M approach and univariate analysis were employed to screen survival-related genes. RT-PCR was employed to examine the expression of DNA repair-related genes in TC samples and matched noncancer samples. CCK-8 analyses were used to determine cellular proliferation. Herein, our team discovered that the expression of four DNA repair-related genes was remarkably upregulated in TC samples in contrast to noncancer samples. Survival assays identified 14 DNA repair-related genes. In our cohort, we observed that the expression of TAF13 and DCTN4 was distinctly elevated in TC specimens in contrast to nontumor specimens. Moreover, knockdown of TAF13 and DCTN4 was observed to inhibit the TC cellular proliferation. Overall, the upregulation of TAF13 and DCTN4 is related to decreased overall survival in TC patients. Therefore, the assessment of TAF13 and DCTN4 expression may be useful for predicting prognosis in these patients.

DNA Repair and Replication-Related Gene Signature Based on Tumor Mutation Burden Reveals Prognostic and Immunotherapy Response in Gastric Cancer

The genomic variant features (mutations, deletions, structural variants, etc.) within gastric cancer impact its evolution and immunogenicity. The tumor has developed several coping strategies to respond to these changes by DNA repair and replication (DRR). However, the intrinsic relationship between the associated DRR-related genes and gastric cancer progression remained unknown. This study selected DRR-related genes with tumor mutation burden based on the TCGA (The Cancer Genome Atlas) database of gastric cancer transcriptome and mutation data. The prognosis model of seven genes (LAMA2, CREB3L3, SELP, ABCC9, CYP1B1, CDH2, and GAMT) was constructed by a univariate and LASSO regression analysis and divided into high-risk and low-risk groups with the median risk score. Survival analysis showed that overall survival (OS) was lower in the high-risk group than that in the low-risk group. Moreover, patients with gastric cancer in the high-risk group have worse survival in different subgroups, including age, gender, histological grade, and TNM stage. The nomogram that included risk scores for DRR-related genes could accurately foresee OS of patients with gastric cancer. Interestingly, the tumor mutation burden score was higher in the low-risk group than that in the high-risk group, and the risk score for DRR-related genes was negatively correlated with tumor mutation burden in gastric cancer. Next, we further combined the risk score and tumor mutation burden to evaluate the prognosis of gastric cancer patients. The low-risk cohort had a better prognosis than the high-risk cohort in the high tumor mutation burden subgroup. The number of mutation types in the high-risk group was lower than that in the low-risk group. In the immune microenvironment of gastric cancer, more naïve B cells, memory resting CD4+ T cells, Treg cells, monocytes cells, and resting mast cells were infiltrated in the high-risk group. At last, PD-L1 and IAP expressions were negatively correlated with the risk scores; patients with gastric cancer in the low-risk group showed better immunotherapy outcomes than those in the high-risk group. Overall, the DRR-related gene signature based on tumor mutation burden is a novel biomarker for prognostic and immunotherapy response in patients with gastric cancer.

ZWINT is a Promising Therapeutic Biomarker Associated with the Immune Microenvironment of Hepatocellular Carcinoma

Background

The prognosis of patients with advanced hepatocellular carcinoma (HCC) is still poor, effective therapeutic targets are needed. ZW10 interacting kinetochore protein (Zwint) is an essential component of the mitotic spindle checkpoint and is upregulated in cancers. Disappointing, the role of ZWINT in HCC has not been fully illuminated.

Methods

Multiple tools, including TIMER2.0, Oncomine, GEPIA2, UALCAN, LinkedOmics, Kaplan-Meier Plotter, cBioPortal, and MethSurv, etc. were applied to comprehensively analyze the expression, genetic alternations, clinicopathological relevance, prognostic value, and DNA methylation of ZWINT, along with its correlations with immune infiltration in HCC. Besides, gene set enrichment analysis (GSEA) and protein-protein interaction (PPI) analysis were performed for the correlated genes of ZWINT, closely interconnected clusters and hub proteins in the PPI network were discovered to learn the underlying biological mechanisms.

Results

We found ZWINT was significantly upregulated in diverse cancers including HCC, compared with the corresponding normal controls. ZWINT upregulation was significantly associated with unfavorable clinicopathological features and survivals of HCC patients. Genetic alternations of ZWINT frequently occurred, which were linked to worse outcomes of HCC patients. The results of GSEA displayed ZWINT and its correlated genes might be components of condensed chromosomes and spindles, which participated in biological processes and signaling pathways involving DNA replication, cytokinesis, and cell cycle checkpoint, etc. Three highly interconnected clusters and 10 hub proteins were identified from the PPI network constructed with the correlated genes of ZWINT. Moreover, ZWINT expression was found positively correlated with infiltration levels of various immune cells, especially myeloid-derived suppressor cells.

Conclusion

This study demonstrated ZWINT might be a promising unfavorable prognostic biomarker and a therapeutic target of HCC, which could regulate HCC progression through cell division and immunosuppression.