DNMIVD: DNA methylation interactive visualization database

HLF is a promising prognostic, immunological, and therapeutic biomarker in human tumors

Despite past research linking HLF mutations to cancer development, no pan-cancer analyses of HLF have been published. As a result, we utilized multiple databases to illustrate the potential roles of HLF in diverse types of cancers. Several databases were used to assess HLF expression in the TCGA cancer samples. Additional assessments were undertaken to investigate the relationship between HLF and overall survival, immune cell infiltration, genetic alterations, promoter methylation, and protein-protein interaction. HLF's putative roles and the relationship between HLF expression and drug reactivity were investigated. HLF expression was shown to be lower in tumor tissues from a variety of malignancies when compared to normal tissues. There was a substantial link found between HLF expression and patient survival, genetic mutations, and immunological infiltration. HLF influenced the pathways of apoptosis, cell cycle, EMT, and PI3K/AKT signaling. Abnormal expression of HLF lowered sensitivity to numerous anti-tumor drugs and small compounds. According to our findings, reduced HLF expression drives cancer growth, and it has the potential to be identified as a vital biomarker for use in prognosis, immunotherapy, and targeted treatment of a range of malignancies.

Oncogenic and immunological role of EDIL3 in human tumours: From pan-cancer analysis to validation in gastric cancer

Background

Epidermal growth factor-like repeats and discoidin I-like domains 3 (EDIL3) is a secreted extracellular matrix protein implicated in diverse physiological and pathological processes including embryonic development, angiogenesis, and anti-inflammatory responses. Recent reports have indicated that EDIL3 play critical roles in carcinogenesis and progression of many cancers. Herein, we performed a pan-cancer investigation to study the potential functions of EDIL3 in various cancers and experimentally validate its function in gastric cancer (GC).

Methods

We analysed EDIL3 expression profiles in different tumours using The Cancer Genome Atlas database. The Kaplan-Meier Plotter was used to investigate the prognostic value of EDIL3, while receiver operating characteristic curve was performed to analyze its diagnostic efficacy. Several bioinformatics tools were used to study the association between EDIL3 and promoter methylation, gene enrichment analysis, immune infiltration, immune-related genes, and drug sensitivity. Molecular biology experiments were conducted to validate the tumorigenic effects of EDIL3.

Results

EDIL3 is variably expressed in different cancers and is closely associated with clinical outcomes. An inverse correlation between EDIL3 and DNA methylation has been observed in 13 cancers. Enrichment analysis indicated that EDIL3 is correlated with many cellular pathways such as extracellular matrix receptor interactions and focal adhesion. EDIL3 was tightly associated with immune infiltration and immune checkpoints. EDIL3 knockdown can promote GC calls apoptosis while preventing proliferation, migration, and invasion in vitro.

Conclusion

EDIL3 is a promising prognostic, diagnostic, and immunological biomarker in various cancers, which could be applied as a new target for cancer therapy.

Epigenetic silencing of LDHB promotes hepatocellular carcinoma by remodeling the tumor microenvironment

Lactate dehydrogenase B (LDHB) reversibly catalyzes the conversion of pyruvate to lactate or lactate to pyruvate and expressed in various malignancies. However, the role of LDHB in modulating immune responses against hepatocellular carcinoma (HCC) remains largely unknown. Here, we found that down-regulation of lactate dehydrogenase B (LDHB) was coupled with the promoter hypermethylation and knocking down the DNA methyltransferase 3A (DNMT 3A) restored LDHB expression levels in HCC cell lines. Bioinformatics analysis of the HCC cohort from The Cancer Genome Atlas revealed a significant positive correlation between LDHB expression and immune regulatory signaling pathways and immune cell infiltrations. Moreover, immune checkpoint inhibitors (ICIs) have shown considerable promise for HCC treatment and patients with higher LDHB expression responded better to ICIs. Finally, we found that overexpression of LDHB suppressed HCC growth in immunocompetent but not in immunodeficient mice, suggesting that the host immune system was involved in the LDHB-medicated tumor suppression. Our findings indicate that DNMT3A-mediated epigenetic silencing of LDHB may contribute to HCC progression through remodeling the tumor immune microenvironment, and LDHB may become a potential prognostic biomarker and therapeutic target for HCC immunotherapy.

Characterization of zinc finger protein 536, a neuroendocrine regulator, using pan-cancer analysis

BACKGROUND

Previous studies suggested that zinc finger protein 536 (ZNF536) was abundant in the central brain and regulated neuronal differentiation. However, the role of ZNF536 in cancer has remained unclear.

METHODS

ZNF536 mutation, copy number alteration, DNA methylation, and RNA expression were explored using public portals. Data from The Cancer Genome Atlas (TCGA) were utilized to analyze pathways and tumor microenvironment (TME), with a focus on prognosis in both TCGA and immunotherapy pan-cancer cohorts. Methylated ZNF536 from small cell lung cancer (SCLC) cell lines were utilized to train with probes for conducting enrichment analysis. Single-cell RNA profile demonstrated the sublocalization and co-expression of ZNF536, and validated its targets by qPCR.

RESULTS

Genetic alterations in ZNF536 were found to be high-frequency and a single sample could harbor different variations. ZNF536 at chromosome 19q12 exerted a bypass effect on CCNE1, supported by CRISPR data. For lung cancer, ZNF536 mutation was associated with longer survival in primary lung adenocarcinoma (LUAD), but its prognosis was poor in metastatic LUAD and SCLC. Importantly, ZNF536 mutation and amplification had opposite prognoses in Stand Up To Cancer-Mark Foundation (SU2C-MARK) LUAD cohort. ZNF536 mutation altered the patterns of genomic alterations in tumors, and had distinct impacts on the signaling pathways and TME compared to ZNF536 amplification. Additionally, ZNF536 expression was predominantly in endocrine tumors and brain tissues. High-dimensional analysis supported this finding and further revealed regulators of ZNF536. Considering that the methylation of ZNF536 was involved in the synaptic pathway associated with neuroendocrine neoplasms, demonstrating both diagnostic and prognostic value. Moreover, we experimentally verified ZNF536 upregulated neuroendocrine markers.

CONCLUSIONS

Our results showed that ZNF536 alterations in cancer, including variations in copy number, mutation, and methylation. We proved the involvement of ZNF536 in neuroendocrine regulation, and identified highly altered ZNF536 as a potential biomarker for immunotherapy.

Identification of immune-associated biomarker for predicting lung adenocarcinoma: bioinformatics analysis and experiment verification of PTK6

BACKGROUND

Abnormal expression of protein tyrosine kinase 6 (PTK6) has been proven to be involved in the development of gynecological tumors. However, its immune-related carcinogenic mechanism in other tumors remains unclear.

OBJECTIVE

The aim of this study was to identify PTK6 as a novel prognostic biomarker in pan-cancer, especially in lung adenocarcinoma (LUAD), which is correlated with immune infiltration, and to clarify its clinicopathological and prognostic significance.

METHODS

The prognostic value and immune relevance of PTK6 were investigated by using bio-informatics in this study. PTK6 expression was validated in vitro experiments (lung cancer cell lines PC9, NCI-H1975, and HCC827; human normal lung epithelial cells BEAS-2B). Western blot (WB) revealed the PTK6 protein expression in lung cancer cell lines. PTK6 expression was inhibited by Tilfrinib. Colony formation and the Cell Counting Kit-8 (CCK-8) assay were used to detect cell proliferation. The wound healing and trans-well were performed to analyze the cell migration capacity. Then flow cytometry was conducted to evaluate the cell apoptosis. Eventually, the relationship between PTK6 and immune checkpoints was examined. WB was used to estimate the PD-L1 expression at different Tilfrinib doses.

RESULTS

PTK6 was an independent predictive factor for LUAD and was substantially expressed in LUAD. Pathological stage was significantly correlated with increased PTK6 expression. In accordance with survival analysis, poor survival rate in LUAD was associated with a high expression level of PTK6. Functional enrichment of the cell cycle and TGF-β signaling pathway was demonstrated by KEGG and GSEA analysis. Moreover, PTK6 expression considerably associated with immune infiltration in LUAD, as determined by immune analysis. Thus, the result of vitro experiments indicated that cell proliferation and migration were inhibited by the elimination of PTK6. Additionally, PTK6 suppression induced cell apoptosis. Obviously, PD-L1 protein expression level up-regulated while PTK6 was suppressed.

CONCLUSION

PTK6 has predictive value for LUAD prognosis, and could up regulated PD-L1.

Construction of lncRNA/Pseudogene-miRNA Network Based on In Silico Approaches for Glycolysis Pathway to Identify Prostate Adenocarcinoma-Related Potential Biomarkers

LncRNAs, pseudogenes, and miRNAs participate a fundamental function in tumorigenesis, metabolism, and invasion of cancer cells, although their regulation of tumor glycolysis in prostate adenocarcinoma (PRAD) is thoroughly not well studied. In this study, we applied transcriptomic, proteomic, and medical information to identify glycolysis-related key genes and modules associated with PRAD. Then, the glycolysis-related lncRNA/lncRNAs/pseudogenes-miRNA-mRNA network was constructed. Analysis of DNA methylation status and expression data determined a DNA methylation-dysregulated three-DE-mRNAs signature for predicting diagnosis, ANGPTL4, GNE, and HSPA in PRAD patients and healthy control. Several lncRNAs/pseudogenes, significantly correlated with the overall survival PVT1, CA5BP1, MIRLET7BHG, SNHG12, and ZNF37BP and disease-free survival status, MALAT1, GUSBP11, MIRLET7BHG, and SNHG1, of patients with PRAD were determined. The methylation profile of DE-lncRNA/pseudogenes was significantly proper for predicting PRAD prognostic model. The transcription level of 6 DE-mRNA ANGPTL4, QSOX1, BIK, CLDN3, DDIT4, and TFF3 was correlated with cancer-related fibroblast infiltration in PRAD. The mutated form of 7 mRNAs, COL5A1, IDH1, HK2, DDIT4, GNE, and QSOX1, was associated with PRAD. In addition to the glycolysis pathway, DE-RNAs play regulatory roles on several pathways, including DNA damage, RTK, cell cycle, RAS/MAPK, TSC/mTOR and PI3K/AKT, AR hormone, and EMT. Overall, our study improves our knowledge of the relation between lncRNAs/pseudogenes and miRNA related to glycolysis and PRAD pathogenesis.

Potential Role of MAP3K14 in Hepatocellular Carcinoma: A Study Based on Comprehensive Bioinformatical Analysis and Validation

According to reports, MAP3K14 is considered an oncogene and is aberrantly expressed in various types of tumor cells. Its abnormal expression is closely associated with the occurrence and progression of various cancers. MAP3K14 also plays a significant role in the development of non-alcoholic steatohepatitis (NASH)-related hepatocellular carcinoma and its connection to tumor stem cells. The prognostic value of MAP3K14 in HCC, as well as its potential functions and roles, requires further elucidation. We evaluated the potential role of MAP3K14 in HCC based on data mining from a range of public databases. The bioinformatics analysis of TCGA, GEO, TIMER, cBioportal, Kaplan-Meier plotter, MethSurv, ENCORI and CellMiner databases was carried out. The expression of MAP3K14 protein in HCC was detected by immunohistochemical method. The mRNA and protein expression levels of MAP3K14 in tumor tissues were higher than those in normal tissues (p < 0.05). The expression of MAP3K14 was correlated with Pathologic T stage (p=0.026), Pathologic stage (p=0.032), Tumor status (p=0.024) and AFP (p=0.002). HCC patients with high expression of MAP3K14 had poor overall survival (OS), progression free survival (PFS) and recurrence free survival (RFS). Multivariate Cox regression analysis showed that the Pathologic stage (p < 0.001) and MAP3K14 expression levels (p < 0.05) is an independent prognostic factor affecting the survival of patients with liver cancer. GO/KEGG analysis suggested that key biological processes (PI3K-Akt signaling pathway) may be the mechanism promoting HCC development. In addition, MAP3K14 was significantly correlated with the infiltrating levels of B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells (p < 0.05). MAP3K14 is up-regulated in HCC and is closely related to the prognosis of HCC patients. MAP3K14 may serve as a potential biomarker for poor prognosis of HCC.

PARD3 drives tumorigenesis through activating Sonic Hedgehog signalling in tumour-initiating cells in liver cancer

BACKGROUND

Par-3 Family Cell Polarity Regulator (PARD3) is a cellular protein essential for asymmetric cell division and polarized growth. This study aimed to study the role of PARD3 in hepatic tumorigenesis.

METHODS

The essential role of PARD3 in mediating hepatic tumorigenesis was assessed in diet-induced spontaneous liver tumour and syngeneic tumour models. The mechanism of PARD3 was delineated by bulk and single-cell RNA sequencing. The clinical significance of PARD3 was identified by tissue array analysis.

RESULTS

PARD3 was overexpressed in tumour tissues and PARD3 overexpression was positively correlated with high tumour stage as well as the poor prognosis in patients. In models of spontaneous liver cancer induced by choline-deficient, amino acid-defined (CDAA) and methionine-choline-deficient (MCD) diets, upregulation of PARD3 was induced specifically at the tumorigenesis stage rather than other early stages of liver disease progression. Site-directed knockout of PARD3 using an adeno-associated virus 8 (AAV8)-delivered CRISPR/Cas9 single-guide RNA (sgRNA) plasmid blocked hepatic tumorigenesis, while PARD3 overexpression accelerated liver tumour progression. In particular, single-cell sequencing analysis suggested that PARD3 was enriched in primitive tumour cells and its overexpression enhanced tumour-initiating cell (TICs). Overexpression of PARD3 maintained the self-renewal ability of the CD133+ TIC population within hepatocellular carcinoma (HCC) cells and promoted the in vitro and in vivo tumorigenicity of CD133+ TICs. Transcriptome analysis revealed that Sonic Hedgehog (SHH) signalling was activated in PARD3-overexpressing CD133+ TICs. Mechanistically, PARD3 interacted with aPKC to further activate SHH signalling and downstream stemness-related genes. Suppression of SHH signalling and aPKC expression attenuated the in vitro and in vivo tumorigenicity of PARD3-overexpressing CD133+ TICs. Tissue array analysis revealed that PARD3 expression was positively associated with the phosphorylation of aPKC, SOX2 and Gli1 and that the combination of these markers could be used to stratify HCC patients into two clusters with different clinicopathological characteristics and overall survival prognoses. The natural compound berberine was selected as a potent suppressor of PARD3 expression and could be used as a preventive agent for liver cancer that completely blocks diet-induced hepatic tumorigenesis in a PARD3-dependent manner.

CONCLUSION

This study revealed PARD3 as a potential preventive target of liver tumorigenesis via TIC regulation.

mHapBrowser: a comprehensive database for visualization and analysis of DNA methylation haplotypes

DNA methylation acts as a vital epigenetic regulatory mechanism involved in controlling gene expression. Advances in sequencing technologies have enabled characterization of methylation patterns at single-base resolution using bisulfite sequencing approaches. However, existing methylation databases have primarily focused on mean methylation levels, overlooking phased methylation patterns. The methylation status of CpGs on individual sequencing reads represents discrete DNA methylation haplotypes (mHaps). Here, we present mHapBrowser, a comprehensive database for visualizing and analyzing mHaps. We systematically processed data of diverse tissues in human, mouse and rat from public repositories, generating mHap format files for 6366 samples. mHapBrowser enables users to visualize eight mHap metrics across the genome through an integrated WashU Epigenome Browser. It also provides an online server for comparing mHap patterns across samples. Additionally, mHap files for all samples can be downloaded to facilitate local processing using downstream analysis toolkits. The utilities of mHapBrowser were demonstrated through three case studies: (i) mHap patterns are associated with gene expression; (ii) changes in mHap patterns independent of mean methylation correlate with differential expression between lung cancer subtypes; and (iii) the mHap metric MHL outperforms mean methylation for classifying tumor and normal samples from cell-free DNA. The database is freely accessible at http://mhap.sibcb.ac.cn/.

MethMarkerDB: a comprehensive cancer DNA methylation biomarker database

DNA methylation plays a crucial role in tumorigenesis and tumor progression, sparking substantial interest in the clinical applications of cancer DNA methylation biomarkers. Cancer-related whole-genome bisulfite sequencing (WGBS) data offers a promising approach to precisely identify these biomarkers with differentially methylated regions (DMRs). However, currently there is no dedicated resource for cancer DNA methylation biomarkers with WGBS data. Here, we developed a comprehensive cancer DNA methylation biomarker database (MethMarkerDB, https://methmarkerdb.hzau.edu.cn/), which integrated 658 WGBS datasets, incorporating 724 curated DNA methylation biomarker genes from 1425 PubMed published articles. Based on WGBS data, we documented 5.4 million DMRs from 13 common types of cancer as candidate DNA methylation biomarkers. We provided search and annotation functions for these DMRs with different resources, such as enhancers and SNPs, and developed diagnostic and prognostic models for further biomarker evaluation. With the database, we not only identified known DNA methylation biomarkers, but also identified 781 hypermethylated and 5245 hypomethylated pan-cancer DMRs, corresponding to 693 and 2172 genes, respectively. These novel potential pan-cancer DNA methylation biomarkers hold significant clinical translational value. We hope that MethMarkerDB will help identify novel cancer DNA methylation biomarkers and propel the clinical application of these biomarkers.

TET2-mediated tumor cGAS triggers endothelial STING activation to regulate vasculature remodeling and anti-tumor immunity in liver cancer

Induction of tumor vascular normalization is a crucial measure to enhance immunotherapy efficacy. cGAS-STING pathway is vital for anti-tumor immunity, but its role in tumor vasculature is unclear. Herein, using preclinical liver cancer models in Cgas/Sting-deficient male mice, we report that the interdependence between tumor cGAS and host STING mediates vascular normalization and anti-tumor immune response. Mechanistically, TET2 mediated IL-2/STAT5A signaling epigenetically upregulates tumor cGAS expression and produces cGAMP. Subsequently, cGAMP is transported via LRRC8C channels to activate STING in endothelial cells, enhancing recruitment and transendothelial migration of lymphocytes. In vivo studies in male mice also reveal that administration of vitamin C, a promising anti-cancer agent, stimulates TET2 activity, induces tumor vascular normalization and enhances the efficacy of anti-PD-L1 therapy alone or in combination with IL-2. Our findings elucidate a crosstalk between tumor and vascular endothelial cells in the tumor immune microenvironment, providing strategies to enhance the efficacy of combinational immunotherapy for liver cancer.

Carcinogenic roles of MAFG-AS1 in human cancers

The MAF bZIP transcription factor G-antisense RNA 1 (MAFG-AS1) is located on chromosome 17. MAFG-AS1 was upregulated in 15 human cancers. MAFG-AS1 not only suppresses 16 miRNAs but also directly impacts 22 protein-coding genes' expression. Notably, abnormal MAFG-AS1 expression is connected to clinicopathological characteristics and a worse prognosis in a variety of cancers. Moreover, MAFG-AS1 takes its part in the tumorigenesis and progression of various human malignancies by suppressing apoptosis and promoting proliferation, migration, invasion, aerobic glycolysis, ferroptosis, angiogenesis, EMT, and metastasis. Besides, it can predict treatment effectiveness in ER + breast cancer, urothelial bladder carcinoma, and liver cancer by functioning as a trigger of resistance to tamoxifen, sorafenib, and cisplatin. This study systematically presents the functions of MAFG-AS1 in various cancers, as well as the findings of bioinformatics analyses of the MAFG-AS1, which should give clear advice for future research.

Sprouty4 is epigenetically upregulated in human colorectal cancer

Sprouty4 (SPRY4) has been frequently reported as a tumor suppressor and is therefore downregulated in various cancers. For the first time, we report that SPRY4 is epigenetically upregulated in colorectal cancer (CRC). In this study, we explored DNA methylation and hydroxymethylation levels of SPRY4 in CRC cells and patient samples and correlated these findings with mRNA and protein expression levels. Three loci within the promoter region of SPRY4 were evaluated for 5mC levels in CRC using the combined bisulfite restriction analysis. In addition, hydroxymethylation levels within SPRY4 were measured in CRC patients. Lastly, DNA methylation and mRNA expression data were extracted from CRC patients in multiple high-throughput data repositories like Gene Expression Omnibus and The Cancer Genome Atlas. Combined in vitro and in silico analysis of promoter methylation levels of SPRY4 clearly demonstrates that the distal promoter region undergoes hypomethylation in CRC patients and is associated with increased expression. Moreover, a decrease in gene body hydroxymethylation and an increase in gene body methylation within the coding region of SPRY4 were found in CRC patients and correlated with increased expression. SPRY4 is epigenetically upregulated in CRC by promoter hypomethylation and hypermethylation within the gene body that warrants future investigation of atypical roles of this established tumor suppressor.

Integrative Analysis of DNA Methylation and Gene Expression Data Identifies Potential Biomarkers and Functional Epigenetic Modules for SARS-CoV-2

To integrate gene expression and DNA methylation data and find the potential role of DNA methylation in the invasion and replication of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). We first conducted differential expression and methylation analysis between the coronavirus disease of 2019 (COVID-19) and healthy controls. FEM was employed to identify functional epigenetic modules, from which a diagnostic model for COVID-19 was built. SKA1 and WSB1 modules were identified, with SKA1 module enriched in COVID-19 replication and transcription, and WSB1 module related to ubiquitin-protein activity. The differentially expressed or differentially methylated genes in these two modules could be used to distinguish COVID-19 from healthy controls, with AUC reaching 1 and 0.98 for SKA1 and WSB1 modules, respectively. Two epigenetically activated genes (CENPM and KNL1) from the SKA1 module were upregulated in HPV- or HBV-positive tumor samples and were found to be significantly associated with the survival of tumor patients. In conclusion, the identified FEM modules and potential signatures play an essential role in the replication and transcription of coronavirus.

DNA-methylation signature accurately differentiates pancreatic cancer from chronic pancreatitis in tissue and plasma

OBJECTIVE

Pancreatic ductal adenocarcinoma (PDAC) is a lethal malignancy. Differentiation from chronic pancreatitis (CP) is currently inaccurate in about one-third of cases. Misdiagnoses in both directions, however, have severe consequences for patients. We set out to identify molecular markers for a clear distinction between PDAC and CP.

DESIGN

Genome-wide variations of DNA-methylation, messenger RNA and microRNA level as well as combinations thereof were analysed in 345 tissue samples for marker identification. To improve diagnostic performance, we established a random-forest machine-learning approach. Results were validated on another 48 samples and further corroborated in 16 liquid biopsy samples.

RESULTS

Machine-learning succeeded in defining markers to differentiate between patients with PDAC and CP, while low-dimensional embedding and cluster analysis failed to do so. DNA-methylation yielded the best diagnostic accuracy by far, dwarfing the importance of transcript levels. Identified changes were confirmed with data taken from public repositories and validated in independent sample sets. A signature of six DNA-methylation sites in a CpG-island of the protein kinase C beta type gene achieved a validated diagnostic accuracy of 100% in tissue and in circulating free DNA isolated from patient plasma.

CONCLUSION

The success of machine-learning to identify an effective marker signature documents the power of this approach. The high diagnostic accuracy of discriminating PDAC from CP could have tremendous consequences for treatment success, once the result from still a limited number of liquid biopsy samples would be confirmed in a larger cohort of patients with suspected pancreatic cancer.

Identification of TMEM178 as a Potential Prognostic Biomarker and Therapeutic Target for Breast Cancer

Background

The transmembrane protein (TMEM) family plays important roles in cancer. However, the expression pattern and biological roles of TMEM178, a member of TMEM family, remains unclear in breast cancer (BRCA).

Methods

Methylation and RNA-seq data were obtained to explore methylation level. Expression of TMEM178, methylation inhibitor 5-Aza-CdR was used to verify the effect of methylation status on the expression of TMEM178. We comprehensively investigated the prognostic outcomes, biological functions and effects on immune cell infiltration of the TMEM178 in BRCA using multiple bioinformatics methods.

Results

The expression of TMEM178 was downregulated and negatively correlated with the level of DNA methylation and DNA methyltransferase (DNMT1, DNMT3A, and DNMT3B) in BRCA. Consistently, TMEM178 mRNA were confirmed to be downregulated, while upregulated in response to treatment with methylation inhibitor 5-Aza-CdR by RT-qPCR. Patients with high expression of TMEM178 have better prognosis and are more sensitive to targeted drug Pazopanib. Immune infiltration analysis showed that the infiltration levels of CD4+ T cell subsets were reduced in BRAC tissues with high TMEM178 expression, and immunosuppressive molecules of T-cell exhaustion were lower expression level.

Conclusion

Hypermethylation of the TMEM178 promoter region was a contributing factor to the downregulation of its expression, and TMEM178 may reflect a prognostic and immunosuppressive situation in BRCA.

SMC3 epigenetic silencing regulates Rab27a expression and drives pancreatic cancer progression

BACKGROUND

Pancreatic ductal adenocarcinoma (PDAC) is expected to soon surpass colorectal cancer as a leading cause of cancer mortality in both males and females in the US, only lagging behind lung cancer. The lethality of PDAC is driven by late diagnosis and inefficient therapies. The complex biology of PDAC involves various cellular components, including exosomes that carry molecular information between cells. Thus, recipient cells can be reprogrammed, impacting tumorigenesis. Rab27a is a GTPase responsible for the last step of exosomes biogenesis. Hence, dissecting the mechanisms that regulate the expression of Rab27a and that control exosomes biogenesis can provide fundamental insights into the molecular underpinnings regulating PDAC progression.

METHODS

To assess the mechanism that regulates Rab27a expression in PDAC, we used PDAC cell lines. The biological significance of these findings was validated in PDAC genetically engineered mouse models (GEMMs) and human samples.

RESULTS

In this work we demonstrate in human PDAC samples and GEMMs that Rab27a expression decreases throughout the development of the disease, and that Rab27a knockout promotes disease progression. What is more, we demonstrate that Rab27a expression is epigenetically regulated in PDAC. Treatment with demethylating agents increases Rab27a expression specifically in human PDAC cell lines. We found that SMC3, a component of the cohesin complex, regulates Rab27a expression in PDAC. SMC3 methylation is present in human PDAC specimens and treatment with demethylating agents increases SMC3 expression in human PDAC cell lines. Most importantly, high levels of SMC3 methylation are associated with a worse prognosis in PDAC. Mechanistically, we identified an enhancer region within the Rab27a gene that recruits SMC3, and modulates Rab27a expression.

CONCLUSION

Overall, we dissected a mechanism that regulates Rab27a expression during PDAC progression and impacts disease prognosis.

Integrated multiomics analysis of chromosome 19 miRNA cluster in bladder cancer

With 46 microRNAs (miRNAs) embedded tandemly over a distance of ~100 kb, chromosome 19 microRNA cluster (C19MC) is the largest miRNA cluster in the human genome. The C19MC is transcribed from a long noncoding genomic region and is usually expressed simultaneously at a higher level. Hence, we performed an integrative multiomics data analysis to examine C19MC regulation, expression patterns, and their impact on bladder cancer (BCa). We found that 43 members of C19MC were highly expressed in BCa. However, its co-localization with recurrent copy number variation (CNV) gain was not statistically significant to implicate its upregulation. It has been reported that C19MC expression is regulated by a well-established CpG island situated 17.6 kb upstream of the transcription start site, but we found that CpG probes at this island were hypomethylated, which was not statistically significant in the BCa cohort. In addition, the promoter region of C19MC is strongly regulated by a group of seven transcription factors (NR2F6, SREBF1, TBP, GATA3, GABPB1, ETV4, and ZNF444) and five chromatin modifiers (SMC3, KDMA1, EZH2, RAD21, and CHD7). Interestingly, these 12 genes were found to be overexpressed in BCa patients. Further, C19MC targeted 42 tumor suppressor (TS) genes that were downregulated, of which 15 were significantly correlated with patient survival. Our findings suggest that transcription factors and chromatin modifiers at the promoter region may regulate C19MC overexpression. The upregulated C19MC members, transcription regulators, and TS genes can be further exploited as potential diagnostic and prognostic indicators as well as for therapeutic management of BCa.

Identification of KNOP1 as a prognostic marker in hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC) is a malignancy with a poor prognosis. This study aimed to evaluate the role and molecular mechanism of lysine-rich nucleolar protein 1 (KNOP1) in HCC.

Methods

Data from The Cancer Genome Atlas (TCGA), genotype-tissue expression (GTEx), and Gene Expression Omnibus (GEO) databases were used to compare KNOP1 expression in normal and HCC tissues. The Human Protein Atlas (HPA) database was used to verify KNOP1 protein expression. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), gene set enrichment, protein-protein and gene-gene interaction network, DNA methylation, genetic alteration, and immune cell infiltration analyses were used to analyze the function and pathway enrichment of KNOP1. Finally, receiver operating characteristic (ROC) curves, Kaplan-Meier (KM) analysis, univariate/multivariate Cox regression analyses, and nomograms were used to predict the clinical and prognostic significance of KNOP1.

Results

KNOP1 expression was higher in HCC tissue samples than in normal specimens. Additionally, high KNOP1 expression was positively correlated with T helper 2 (Th2) cells and immune checkpoints. KM analysis, Cox regression analysis, and nomogram prognostic model prediction suggested that high KNOP1 expression is a risk factor for poor HCC prognosis.

Conclusions

KNOP1 overexpression is associated with poor HCC prognosis and increased proportions of immune cell infiltration and checkpoints. KNOP1 is a potential biomarker for evaluating HCC prognosis.

Exploration of prognosis and immunometabolism landscapes in ER+ breast cancer based on a novel lipid metabolism-related signature

Introduction

Lipid metabolic reprogramming is gaining attention as a hallmark of cancers. Recent mounting evidence indicates that the malignant behavior of breast cancer (BC) is closely related to lipid metabolism. Here, we focus on the estrogen receptor-positive (ER+) subtype, the most common subgroup of BC, to explore immunometabolism landscapes and prognostic significance according to lipid metabolism-related genes (LMRGs).

Methods

Samples from The Cancer Genome Atlas (TCGA) database were used as training cohort, and samples from the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC), Gene Expression Omnibus (GEO) datasets and our cohort were applied for external validation. The survival-related LMRG molecular pattern and signature were constructed by unsupervised consensus clustering and least absolute shrinkage and selection operator (LASSO) analysis. A lipid metabolism-related clinicopathologic nomogram was established. Gene enrichment and pathway analysis were performed to explore the underlying mechanism. Immune landscapes, immunotherapy and chemotherapy response were further explored. Moreover, the relationship between gene expression and clinicopathological features was assessed by immunohistochemistry.

Results

Two LMRG molecular patterns were identified and associated with distinct prognoses and immune cell infiltration. Next, a prognostic signature based on nine survival-related LMRGs was established and validated. The signature was confirmed to be an independent prognostic factor and an optimal nomogram incorporating age and T stage (AUC of 5-year overall survival: 0.778). Pathway enrichment analysis revealed differences in immune activities, lipid biosynthesis and drug metabolism by comparing groups with low- and high-risk scores. Further exploration verified different immune microenvironment profiles, immune checkpoint expression, and sensitivity to immunotherapy and chemotherapy between the two groups. Finally, arachidonate 15-lipoxygenase (ALOX15) was selected as the most prominent differentially expressed gene between the two groups. Its expression was positively related to larger tumor size, more advanced tumor stage and vascular invasion in our cohort (n = 149).

Discussion

This is the first lipid metabolism-based signature with value for prognosis prediction and immunotherapy or chemotherapy guidance for ER+ BC.

Creation of a Prognostic Model Using Cuproptosis-Associated Long Noncoding RNAs in Hepatocellular Carcinoma

Cuproptosis is an unusual form of cell death caused by copper accumulation in mitochondria. Cuproptosis is associated with hepatocellular carcinoma (HCC). Long noncoding RNAs (LncRNAs) have been shown to be effective prognostic biomarkers, yet the link between lncRNAs and cuproptosis remains unclear. We aimed to build a prognostic model of lncRNA risk and explore potential biomarkers of cuproptosis in HCC. Pearson correlations were used to derive lncRNAs co-expressed in cuproptosis. The model was constructed using Cox, Lasso, and multivariate Cox regressions. Kaplan-Meier survival analysis, principal components analysis, receiver operating characteristic curve, and nomogram analyses were carried out for validation. Seven lncRNAs were identified as prognostic factors. A risk model was an independent prognostic predictor. Among these seven lncRNAs, prostate cancer associated transcript 6 (PCAT6) is highly expressed in different types of cancer, activating Wnt, PI3K/Akt/mTOR, and other pathways; therefore, we performed further functional validation of PCAT6 in HCC. Reverse transcription-polymerase chain reaction results showed that PCAT6 was aberrantly highly expressed in HCC cell lines (HepG2 and Hep3B) compared to LO2 (normal hepatocytes). When its expression was knocked down, cells proliferated and migrated less. PCAT6 might be a potential biomarker for predicting prognosis in HCC.

Identification of MDM2 as a prognostic and immunotherapeutic biomarker in a comprehensive pan-cancer analysis: A promising target for breast cancer, bladder cancer and ovarian cancer immunotherapy

BACKGROUND

The murine double minute 2 (MDM2) gene is a crucial factor in the development and progression of various cancer types. Multiple rigorous scientific studies have consistently shown its involvement in tumorigenesis and cancer progression in a wide range of cancer types. However, a comprehensive analysis of the role of MDM2 in human cancer has yet to be conducted.

METHODS

We used various databases, including TIMER2.0, TCGA, GTEx and STRING, to analyze MDM2 expression and its correlation with clinical outcomes, interacting genes and immune cell infiltration. We also investigated the association of MDM2 with immune checkpoints and performed gene enrichment analysis using DAVID tools.

RESULTS

The pan-cancer MDM2 analysis found that MDM2 expression and mutation status were observably different in 25 types of cancer tissue compared with healthy tissues, and prognosis analysis showed that there was a significant correlation between MDM2 expression and patient prognosis. Furthermore, correlation analysis showed that MDM2 expression was correlated with tumor mutational burden, microsatellite instability and drug sensitivity in certain cancer types. We found that there was an association between MDM2 expression and immune cell infiltration across cancer types, and MDM2 inhibitors might enhance the effect of immunotherapy on breast cancer, bladder cancer and ovarian cancer.

CONCLUSIONS

The first systematic pan-cancer analysis of MDM2 was conducted, and it demonstrated that MDM2 was a reliable prognostic biomarker and was closely related to cancer immunity, providing a potential immunotherapeutic target for breast cancer, bladder cancer and ovarian cancer.

Downregulation of RAB17 have a poor prognosis in kidney renal clear cell carcinoma and its expression correlates with DNA methylation and immune infiltration

BACKGROUND

RAB17 is one of the RAB family members. It has been reported to be closely associated with a variety of tumors and has different roles in various tumors. However, the effect of RAB17 in KIRC remains unclear.

MATERIALS AND METHODS

We analyzed the differential expression of RAB17 in kidney renal clear cell carcinoma (KIRC) tissues and normal tissues using the public databases. The prognostic role of RAB17 in KIRC was analyzed using the Cox regression methods, and a prognostic model was constructed based on the results of the Cox analysis. In addition, further analysis of RAB17 in KIRC was performed in relation to genetic alterations, DNA methylation m6A methylation and immune infiltration. Finally, RAB17 mRNA and protein expression levels were analyzed in tissue samples (KIRC tissues and normal tissues) and cell lines (normal renal tubular cell and KIRC cells), and in vitro functional assays were performed.

RESULTS

RAB17 was low-expressed in KIRC. Downregulation of RAB17 expression is correlated with unfavorable clinicopathological characteristics and a worse prognosis in KIRC. The RAB17 gene alteration in KIRC was primarily characterized by copy number alteration. Six CpG sites of RAB17 DNA methylation levels are higher in KIRC tissues than in normal tissues, and are correlated with RAB17 mRNA expression levels, showing a significant negative correlation. cg01157280 site DNA methylation levels are associated with pathological stage and overall survival, and it may be the only CpG site with independent prognostic significance. Functional mechanism analysis revealed that RAB17 is closely associated with immune infiltration. RAB17 expression was found to be negatively correlated with most immune cell infiltration according to two different methods. Furthermore, most immunomodulators were significantly negatively correlated with RAB17 expression, and significantly positively correlated with RAB17 DNA methylation levels. RAB17 was significantly low expression in KIRC cells and KIRC tissues. In vitro, silencing of RAB17 promoted KIRC cell migration.

CONCLUSION

RAB17 can be used as a potential prognostic biomarker for patients with KIRC and for assessing immunotherapy response.

LEM Domain Containing 1 Acts as a Novel Oncogene and Therapeutic Target for Triple-Negative Breast Cancer

Breast cancer is the most common deadly malignancy in women worldwide. In particular, triple-negative breast cancer (TNBC) exhibits the worst prognosis among four subtypes of breast cancer due to limited treatment options. Exploring novel therapeutic targets holds promise for developing effective treatments for TNBC. Here, we demonstrated for the first time that LEMD1 (LEM domain containing 1) is highly expressed in TNBC and contributes to reduced survival in TNBC patients, through analysis of both bioinformatic databases and collected patient samples. Furthermore, LEMD1 silencing not only inhibited the proliferation and migration of TNBC cells in vitro, but also abolished tumor formation of TNBC cells in vivo. Knockdown of LEMD1 enhanced the sensitivity of TNBC cells to paclitaxel. Mechanistically, LEMD1 promoted the progress of TNBC by activating the ERK signaling pathway. In summary, our study revealed that LEMD1 may act as a novel oncogene in TNBC, and targeting LEMD1 may be exploited as a promising therapeutic approach to enhance the efficacy of chemotherapy against TNBC.

DNA methylation patterns at birth predict health outcomes in young adults born very low birthweight

Background

Individuals born very low birthweight (VLBW) are at increased risk of impaired cardiovascular and respiratory function in adulthood. To identify markers to predict future risk for VLBW individuals, we analyzed DNA methylation at birth and at 28 years in the New Zealand (NZ) VLBW cohort (all infants born < 1500 g in NZ in 1986) compared with age-matched, normal birthweight controls. Associations between neonatal methylation and cardiac structure and function (echocardiography), vascular function and respiratory outcomes at age 28 years were documented.

Results

Genomic DNA from archived newborn heel-prick blood (n = 109 VLBW, 51 controls) and from peripheral blood at ~ 28 years (n = 215 VLBW, 96 controls) was analyzed on Illumina Infinium MethylationEPIC 850 K arrays. Following quality assurance and normalization, methylation levels were compared between VLBW cases and controls at both ages by linear regression, with genome-wide significance set to p < 0.05 adjusted for false discovery rate (FDR, Benjamini-Hochberg). In neonates, methylation at over 16,400 CpG methylation sites differed between VLBW cases and controls and the canonical pathway most enriched for these CpGs was Cardiac Hypertrophy Signaling (p = 3.44E⁻¹¹). The top 20 CpGs that differed most between VLBW cases and controls featured clusters in ARID3A, SPATA33, and PLCH1 and these 3 genes, along with MCF2L, TRBJ2-1 and SRC, led the list of 15,000 differentially methylated regions (DMRs) reaching FDR-adj significance. Fifteen of the 20 top CpGs in the neonate EWAS showed associations between methylation at birth and adult cardiovascular traits (particularly LnRHI). In 28-year-old adults, twelve CpGs differed between VLBW cases and controls at FDR-adjusted significance, including hypermethylation in EBF4 (four CpGs), CFI and UNC119B and hypomethylation at three CpGs in HIF3A and one in KCNQ1. DNA methylation GrimAge scores at 28 years were significantly greater in VLBW cases versus controls and weakly associated with cardiovascular traits. Four CpGs were identified where methylation differed between VLBW cases and controls in both neonates and adults, three reversing directions with age (two CpGs in EBF4, one in SNAI1 were hypomethylated in neonates, hypermethylated in adults). Of these, cg16426670 in EBF4 at birth showed associations with several cardiovascular traits in adults.

Conclusions

These findings suggest that methylation patterns in VLBW neonates may be informative about future adult cardiovascular and respiratory outcomes and have value in guiding early preventative care to improve adult health.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-023-01463-3.

Impact of BSG/CD147 gene expression on diagnostic, prognostic and therapeutic strategies towards malignant cancers and possible susceptibility to SARS-CoV-2

BACKGROUND

BSG (CD147) is a member of the immunoglobulin superfamily that shows roles for potential prognostics and therapeutics for metastatic cancers and SARS-CoV-2 invasion for COVID-19. The susceptibility of malignant cancers to SARS-CoV-2 as well as the correlations between disease outcome and BSG expression in tumor tissues have not been studied in depth.

METHODS

In this study, we explored the BSG expression profile, survival correlation, DNA methylation, mutation, diagnostics, prognostics, and tumor-infiltrating lymphocytes (TILs) from different types of cancer tissues with corresponding healthy tissues. In vitro studies for cordycepin (CD), N6-(2-hydroxyethyl) adenosine (HEA), N6, N6-dimethyladenosine (m⁶₂A) and 5'-uridylic acid (UMP) on BSG expression were also conducted.

RESULTS

We revealed that BSG is conserved among different species, and significantly upregulated in seven tumor types, including ACC, ESCA, KICH, LIHC, PAAD, SKCM and THYM, compared with matched normal tissues, highlighting the susceptibility of these cancer patients to SARS-CoV-2 invasion, COVID-19 severity and progression of malignant cancers. High expression in BSG was significantly correlated with a short OS in LGG, LIHC and OV patients, but a long OS in KIRP patients. Methylation statuses in the BSG promoter were significantly higher in BRCA, HNSC, KIRC, KIRP, LUSC, PAAD, and PRAD tumor tissues, but lower in READ. Four CpGs in the BSG genome were identified as potential DNA methylation biomarkers which could be used to predict malignant cancers from normal individuals. Furthermore, a total of 65 mutation types were found, in which SARC showed the highest mutation frequency (7.84%) and THYM the lowest (0.2%). Surprisingly, both for disease-free and progression-free survival in pan-cancers were significantly reduced after BSG mutations. Additionally, a correlation between BSG expression and immune lymphocytes of CD56bright natural killer cell, CD56dim natural killer cell and monocytes, MHC molecules of HLA-A, HLA-B, HLA-C and TAPBP, immunoinhibitor of PVR, PVRL2, and immunostimulators of TNFRSF14, TNFRSF18, TNFRSF25, and TNFSF9, was revealed in most cancer types. Moreover, BSG expression was downregulated by CD, HEA, m⁶₂A or UMP in cancer cell lines, suggesting therapeutic potentials for interfering entry of SARS-CoV-2.

CONCLUSIONS

Altogether, our study highlights the values of targeting BSG for diagnostic, prognostic and therapeutic strategies to fight malignant cancers and COVID-19. Small molecules CD, HEA, m⁶₂A and UMP imply therapeutic potentials in interfering with entry of SARS-CoV-2 and progression of malignant cancers.

Effect of DPP4/CD26 expression on SARS‑CoV‑2 susceptibility, immune response, adenosine (derivatives m62A and CD) regulations on patients with cancer and healthy individuals

The worldwide COVID‑19 pandemic was brought on by a new coronavirus (SARS Cov‑2). A marker/receptor called Dipeptidyl peptidase 4/CD26(DPP4/CD26) may be crucial in determining susceptibility to tumors and coronaviruses. However, the regulation of DPP4 in COVID‑invaded cancer patients and its role on small molecule compounds remain unclear. The present study used the Human Protein Atlas, Monaco, and Schmiedel databases to analyze the expression of DPP4 in human tissues and immune cells. The association between DPP4 expression and survival in various tumor tissues was compared using GEPIA 2. The DNMIVD database was used to analyze the correlation between DPP4 expression and promoter methylation in various tumors. On the cBioPortal network, the frequency of DPP4 DNA mutations in various cancers was analyzed. The correlation between DPP4 expression and immunomodulators was analyzed by TISIDB database. The inhibitory effects of cordycepin (CD), N6, N6‑dimethyladenosine (m⁶₂A) and adenosine (AD) on DPP4 in cancer cells were evaluated. DPP4 was mainly expressed in endocrine tissue, followed by gastrointestinal tract, female tissue (mainly in placenta), male tissue (mainly in prostate and seminal vesicle), proximal digestive tract, kidney, bladder, liver, gallbladder and respiratory system. In immune cells, DPP4 mRNA was mainly expressed in T cells, and its expression was upregulated in esophageal carcinoma, kidney renal papillary cell carcinoma (KIRP), liver hepatocellular carcinoma (LIHC), lung adenocarcinoma, pancreatic adenocarcinoma, prostate adenocarcinoma, stomach adenocarcinoma, thyroid carcinoma and thymoma. However, it was downregulated in breast invasive carcinoma, kidney chromophobe, lung squamous cell carcinoma and skin cutaneous melanoma. Thus, DPP4 is involved in viral invasion in most types of cancer. The expression of DPP4 could be inhibited by CD, m⁶₂A and AD in different tumor cells. Moreover, CD significantly inhibited the formation of GFP‑positive syncytial cells. In vivo experiments with AD injection further showed that AD significantly inhibited lymphocyte activating factor 3 expression. These drugs may have potential to treat COVID‑19 by targeting DPP4. Thus, DPP4 may be medically significant for SARS‑CoV‑2‑infected cancer patients, providing prospective novel targets and concepts for the creation of drugs against COVID‑19.

Hypermethylated WASF2: tumor suppressive role in head and neck squamous cell carcinoma

Background

WASF2 regulates actin reorganization during cell migration. WASF2 has been identified as a regulator of the development of gastric cancer, breast cancer, and pancreatic cancer. But its regulatory mechanisms remain unknown. Also, its function was absent in head and neck squamous cell carcinoma (HNSCC). Consequently, we examined the effect of DNA methylation on aberrant WASF2 expression in HNSCC.

Methods

TNMplot, TIMER, GSEA pathway analysis, and the Kaplan-Meier Plotter database were used to analyze the expression, function, and prognostic value of WASF2, as well as the correlation between WASF2 and infiltrating immune cells in HNSCC or pan-cancer analysis. WASF2 promoter methylation levels and the correlation between WASF2 expression and WASF2 promoter methylation in HNSCC were evaluated using the DNMIVD database. The effect of DNA methylation inhibitor on WASF2 expression was demonstrated in the GEO database. Finally, the TISIDB database determined the relationships between WASF2 methylation, immune cell infiltration, and immune inhibitors.

Results

WASF2 was significantly downregulated in HNSCC tissues where WASF2 promoter methylation was elevated. According to the GEO database, treatment with a DNA methylation inhibitor notably restored the mRNA expression of WASF2. WASF2 expression was also a favorable indicator of human papilloma virus (HPV)-positive HNSCC. Its level of promoter methylation had detrimental effects on patient survival. In addition, patients with elevated levels of WASF2 demonstrated active G2/M regulation, TGF-β signaling, Kras signaling, fatty acid metabolism, and p53 pathways. WASF2 was positively associated with tumor-killing immune cells, while WASF2 methylation was positively related to immunosuppressive cells and immune-inhibitors.

Conclusions

Hypermethylated WASF2 acted as a tumor suppressor of HNSCC by regulating tumor formation and immune imbalance.

A pan-cancer analysis of the oncogenic role of ERCC6L

BACKGROUND

Excision repair cross-complementation group 6 like (ERCC6L), a polo-like kinase 1 (PLK1)-interacting checkpoint helicase, confers a high risk of cancer and enhances the progression of a variety of cancers. The present investigation aimed to elucidate the pan-cancer expression patterns of ERCC6L and to examine the possibility of using this gene for patient diagnosis and prognosis.

METHODS

The expression patterns of ERCC6L in normal and cancer patients at various clinical stages were explored based on TCGA datasets. Subsequently, Bioinformatics techniques were then used to analyze patient's survival probabilities, Cox multivariate clinical parameters, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) terms related to ERCC6L, the correlation between mRNA expression levels and patient survival, genetic alterations or somatic mutations of ERCC6L, and immune infiltration.

RESULTS

Most cancer types had higher ERCC6L mRNA levels than normal tissue. Higher ERCC6L expression levels were correlated with poor prognosis for cancer patients. Thus, ERCC6L may serve as an effective diagnostic and prognostic marker for multiple cancers. Moreover, ERCC6L expression levels were higher in patients with higher clinical tumor grades and were associated with poor prognoses at these stages. GO and KEGG analyses revealed a correlation between ERCC6L expression levels and chromatin and cell cycle events. We also found that the mRNA expression level of the ERCC6L promoter and a favorable prognosis was negatively correlated with the promoter's methylation but not with copy number variation. A quantitative analysis of immune infiltration suggested a positive correlation between ERCC6L levels and the infiltration of Th2 immune cells in main cancer types. Finally, we examined the ERCC6L somatic mutations, especially single-nucleotide variants, and ERCC6L expression-related drug sensitivity.

CONCLUSIONS

Herein, we reported a comprehensive investigation of the tumor-promoting role of ERCC6L in various cancer types. ERCC6L is a candidate biomarker for diagnosing and unfavorable prognosis of specific cancers.

Antiviral Potential of Small Molecules Cordycepin, Thymoquinone, and N6, N6-Dimethyladenosine Targeting SARS-CoV-2 Entry Protein ADAM17

COVID-19 is an acute respiratory disease caused by SARS-CoV-2 that has spawned a worldwide pandemic. ADAM17 is a sheddase associated with the modulation of the receptor ACE2 of SARS-CoV-2. Studies have revealed that malignant phenotypes of several cancer types are closely relevant to highly expressed ADAM17. However, ADAM17 regulation in SARS-CoV-2 invasion and its role on small molecules are unclear. Here, we evaluated the ADAM17 inhibitory effects of cordycepin (CD), thymoquinone (TQ), and N6, N6-dimethyladenosine (m⁶₂A), on cancer cells and predicted the anti-COVID-19 potential of the three compounds and their underlying signaling pathways by network pharmacology. It was found that CD, TQ, and m⁶₂A repressed the ADAM17 expression upon different cancer cells remarkably. Moreover, CD inhibited GFP-positive syncytia formation significantly, suggesting its potential against SARS-CoV-2. Pharmacological analysis by constructing CD-, TQ-, and m⁶₂A-based drug-target COVID-19 networks further indicated that ADAM17 is a potential target for anti-COVID-19 therapy with these compounds, and the mechanism might be relevant to viral infection and transmembrane receptors-mediated signal transduction. These findings imply that ADAM17 is of potentially medical significance for cancer patients infected with SARS-CoV-2, which provides potential new targets and insights for developing innovative drugs against COVID-19.

Integrated bioinformatic analysis to understand the association between phthalate exposure and breast cancer progression

Phthalates have been extensively used as plasticizers while manufacturing plastic-based consumer products. Estradiol mimicking properties and association studies suggest phthalates may contribute to breast cancer (BC). We performed an in-silico analysis and functional studies to understand the association between phthalate exposure and BC progression. Search for phthalate-responsive genes using the comparative toxicogenomics database identified 20 genes as commonly altered in response to multiple phthalates exposure. Of the 20 genes, 12 were significantly differentially expressed between normal and BC samples. In BC samples, 9 out of 20 genes showed a negative correlation between promoter methylation and its expression. AHR, BAX, BCL2, CAT, ESR2, IL6, and PTGS2 expression differed significantly between metastatic and non-metastatic BC samples. Gene set enrichment analysis identified metabolism, ATP-binding cassette transporters, insulin signaling, and type II diabetes as highly enriched pathways. The diagnostic assessment based on 20 genes expression suggested a sensitivity and a specificity >0.91. The aberrantly expressed phthalate interactive gene influenced the overall survival of BC patients. Drug-gene interaction analysis identified 14 genes and 523 candidate drugs, including 19 BC treatment-approved drugs. Di(2-ethylhexyl) phthlate (DEHP) exposure increased the growth, proliferation, and migration of MCF-7 and MDA-MB-231 cells in-vitro. DEHP exposure induced morphological changes, actin cytoskeletal remodeling, increased ROS content, reduced basal level lipid peroxidation, and induced epithelial to mesenchymal transition (EMT). The present approach can help to explore the potentially damaging effects of environmental agents on cancer risk and understand the underlined pathways and molecular mechanisms.

CanMethdb: a database for genome-wide DNA methylation annotation in cancers

MOTIVATION

DNA methylation within gene body and promoters in cancer cells is well documented. An increasing number of studies showed that cytosine-phosphate-guanine (CpG) sites falling within other regulatory elements could also regulate target gene activation, mainly by affecting transcription factors (TFs) binding in human cancers. This led to the urgent need for comprehensively and effectively collecting distinct cis-regulatory elements and TF-binding sites (TFBS) to annotate DNA methylation regulation.

RESULTS

We developed a database (CanMethdb, http://meth.liclab.net/CanMethdb/) that focused on the upstream and downstream annotations for CpG-genes in cancers. This included upstream cis-regulatory elements, especially those involving distal regions to genes, and TFBS annotations for the CpGs and downstream functional annotations for the target genes, computed through integrating abundant DNA methylation and gene expression profiles in diverse cancers. Users could inquire CpG-target gene pairs for a cancer type through inputting a genomic region, a CpG, a gene name, or select hypo/hypermethylated CpG sets. The current version of CanMethdb documented a total of 38 986 060 CpG-target gene pairs (with 6 769 130 unique pairs), involving 385 217 CpGs and 18 044 target genes, abundant cis-regulatory elements and TFs for 33 TCGA cancer types. CanMethdb might help biologists perform in-depth studies of target gene regulations based on DNA methylations in cancer.

AVAILABILITY AND IMPLEMENTATION

The main program is available at https://github.com/chunquanlipathway/CanMethdb.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Pan-Cancer analysis shows that ACO2 is a potential prognostic and immunotherapeutic biomarker for multiple cancer types including hepatocellular carcinoma

BACKGROUND

Recent evidence increasingly suggests key roles for the tricarboxylic acid cycle and fatty acid metabolism in tumor progression and metastasis. Aconitase 2 (ACO2) is a component of the tricarboxylic acid cycle and represents a key cellular metabolic hub that promotes de novo fatty acid biosynthesis. However, there have been few reports on the role of ACO2 in tumorigenesis and cancer progression.

METHODS

Through the comprehensive use of datasets from The Cancer Genome Atlas, Genotype-Tissue Expression Project, cBioPortal, Human Protein Atlas, UALCAN, Gene Expression Profiling Interactive Analysis, DNA Methylation Interactive Visualization Database, and TIMER2, we adopted bioinformatics methods to uncover the potential carcinogenic roles of ACO2, including by analysing ACO2 expression and correlations between prognosis, genetic mutations, immune cell infiltration, DNA methylation, tumor mutational burden, and microsatellite instability in different tumors. Additionally, the expression level and tumor-promoting effect of ACO2 were verified in hepatocellular carcinoma (HCC) cells. To explore the underlying mechanisms of ACO2 in human cancer, ACO2-related gene enrichment analysis and lipid metabolomics were performed using LM3 cells with or without ACO2 knockdown.

RESULTS

The results indicated that ACO2 was highly expressed in most cancers, showing early diagnostic value in six tumor types, and was positively or negatively associated with prognosis in different tumors. Moreover, ACO2 expression was associated with immune cell infiltration, such as CD8+ T cells and tumor-associated neutrophils, in some cancers. For most cancer types, there was a significant association between immune checkpoint-associated genes and ACO2 expression. Compared with normal hepatocytes, ACO2 was upregulated in HCC cells, which promoted their proliferation and migration. Furthermore, to explore the underlying molecular mechanism, we performed KEGG pathway enrichment analysis of ACO2-associated genes and lipidomics using LM3 cells with or without ACO2 knockdown, which screened 19 significantly altered metabolites, including 17 with reduced levels and 2 with increased levels.

CONCLUSION

Through pan-cancer analysis, we discovered for the first time and verified that ACO2 could be a useful diagnostic biomarker for cancer detection. Additionally, ACO2 could be used as an auxiliary prognostic marker or as a marker for immunotherapy in some tumor types.

MAPK8IP2 is a potential prognostic biomarker and promote tumor progression in prostate cancer

Background

MAPK8IP2 is one of the JNK-interacting proteins (JIPs) family members, and is involved in the regulation of the JNK and P38 MAPK signaling pathways. MAPK8IP2 has been reported to be closely associated with several cancers. However, the biological function of MAPK8IP2 in prostate cancer (PCa) remains unclear.

Methods

MAPK8IP2 expression in PCa and subgroups of PCa was analyzed by public databases. The prognostic role of MAPK8IP2 in prostate cancer was analyzed using the Cox regression method. The potential mechanism by which MAPK8IP2 affects PCa progression was investigated by utilizing public data, including genetic alteration, DNA methylation, m6A methylation, and immune infiltration data. We further performed in vitro assays to validate the effect of MAPK8IP2 on PCa cell proliferation, migration and invasion.

Results

MAPK8IP2 is highly expressed in PCa tissues. Overexpression of MAPK8IP2 is associated with adverse clinicopathological factors and a poor prognosis in PCa. Receiver operating curve analysis showed that MAPK8IP2 can distinguish PCa tissues from non-PCa tissues with a certain accuracy (AUC = 0.814). The MAPK8IP2 genetic alteration rate was 2.6% and MAPK8IP2 alterations correlated with a poor prognosis. We also found that CDK12 and TP53 mutations were associated with MAPK8IP2 expression. The DNA methylation level of MAPK8IP2 was higher in primary tumors than in normal tissues, and the high MAPK8IP2 DNA methylation group of PCa patients had poor survival. Enrichment analysis indicated that MAPK8IP2 was involved in the MAPK signaling pathway. In vitro, knockdown of MAPK8IP2 inhibited PCa cell proliferation, migration and invasion.

Conclusion

MAPK8IP2 is a potential target for PCa treatment and can serve as a novel biomarker for PCa diagnosis and prognosis evaluation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-10259-2.

Prognosis and Immunological Characteristics of PGK1 in Lung Adenocarcinoma: A Systematic Analysis

Background: Aerobic glycolysis plays a key role in tumor metabolic reprogramming to reshape the immune microenvironment. The phosphoglycerate kinase 1 (PGK1) gene codes a glycolytic enzyme that converts 1,3-diphosphoglycerate to 3-phosphoglycerate. However, in lung adenocarcinoma (LUAD), the role of PGK1 in altering the tumor microenvironment (TME) has not yet been determined. Methods: Raw data, including bulk DNA and mRNA-seq data, methylation modification data, single-cell RNA-seq data, proteomics data, clinical case characteristics survival, immunotherapy data, and so on, were obtained from multiple independent public data sets. These data were reanalyzed to uncover the prognosis and immunological characteristics of PGK1 in LUAD. Results: We found that PGK1 mRNA and protein were considerably over-expressed in LUAD compared to normal tissue and that high PGK1 expression is associated with poorer prognostic outcomes in LUAD. The enrichment analysis of PGK1 co-expressed genes in lung adenocarcinoma revealed that PGK1 may be involved in hypoxia, metabolism, DNA synthesis, cell cycle, PI3K/AKT, and various immune and inflammatory signaling pathways. Furthermore, PGK1 is also linked to the recruitment of numerous immune cells, including aDC (dendritic cells), macrophages, and neutrophils. More importantly, PGK1 was highly expressed in immunosuppressive cells, including M2 macrophages, Tregs, and exhausted T cells, among others. Finally, higher PGK1 expression indicated significant correlations to immune checkpoints, TMB (tumor mutation burden), and high response to immunotherapy. Conclusions: The presented findings imply that PGK1, as a glycolysis core gene, may be important for the modification of the immune microenvironment by interacting with the tumor metabolism. The results of this study provide clues for a potential immunometabolic combination therapy strategy in LUAD, for which more experimental and clinical translational research is needed.

Pan-cancer analysis of the prognostic and immunological roles of DEAD-box helicase 5 (DDX5) in human tumors

Background: Recent studies have demonstrated the significance of the DEAD-box helicase 5 (DDX5) gene, which is involved in pathways concerning the modification of RNA structures. DDX5 functions as a coregulator of cellular transcription and splicing, and participates in the processing of small noncoding RNAs. The aberrant regulation of DDX5 expression possibly plays a significant role in the genesis of cancer. However, there are no comprehensive pan-cancer studies on DDX5. This study is the first to conduct a pan-cancer analysis of DDX5 for aiding the diagnosis and treatment of cancer.

Methods: The gene expression, genetic alterations, protein phosphorylation, promoter methylation, immune infiltration, and enrichment analyses of DDX5 were performed using data retrieved from The Cancer Genome Atlas (TCGA), Genotype-tissue Expression (GTEx), Human Protein Atlas (HPA), Tumor Immunological Estimation Resource 2.0 (TIMER2.0), Gene Expression Profiling Interactive Analysis (GEPIA), DNA methylation interactive visualization database (DNMIVD), and Search Tool for the Retrieval of Interaction Genes/Proteins (STRING). Data analyses were performed with the R software and other webtools.

Results: The expression of DDX5 mRNA decreased significantly in 17 cancer types, but increased significantly in eight cancer types. The enhanced expression of DDX5 mRNA in the tumor samples was related to decreased overall survival (OS), progression-free interval (PFI), and disease-specific survival (DSS) in three cancers, but increased OS, PFI, and DSS in other cancers. The DNA promoter methylation level was significantly reduced in eight cancer types, and there were exceptions in the methylation levels of the DDX5 promoter in four cancer types. The expression of DDX5 mRNA was highly correlated with the infiltration of CD8+ T cells, cancer-associated fibroblasts, and B cells in a wide variety of malignancies. The findings revealed a strong association between DDX5 and its co-expressed genes in numerous cancer types. Enrichment analysis suggested that DDX5 was associated with multiple cellular pathways, including RNA splicing, Notch signaling pathway, and viral carcinogenesis, which was consistent with the results of previous studies.

Conclusion: The findings obtained herein provide further information on the oncogenic potential of DDX5 in diverse tumor types. We propose that DDX5 has important roles in tumor immunity and the diagnosis of cancer.

Pan-cancer analysis of cuproptosis regulation patterns and identification of mTOR-target responder in clear cell renal cell carcinoma

BACKGROUND

The mechanism of cuproptosis, a novel copper-induced cell death by regulating tricarboxylic acid cycle (TCA)-related genes, has been reported to regulate oxidative phosphorylation system (OXPHOS) in cancers and can be regarded as potential therapeutic strategies in cancer; however, the characteristics of cuproptosis in pan-cancer have not been elucidated.

METHODS

The multi-omics data of The Cancer Genome Atlas were used to evaluate the cuproptosis-associated characteristics across 32 tumor types. A cuproptosis enrichment score (CEScore) was established using a single sample gene enrichment analysis (ssGSEA) in pan-cancer. Spearman correlation analysis was used to identify pathway most associated with CEScore. Lasso-Cox regression was used to screen prognostic genes associated with OXPHOS and further construct a cuproptosis-related prognostic model in clear cell renal cell carcinoma (ccRCC).

RESULTS

We revealed that most cuproptosis-related genes (CRGs) were differentially expressed between tumors and normal tissues, and somatic copy number alterations contributed to their aberrant expression. We established a CEScore index to indicate cuproptosis status which was associated with prognosis in most cancers. The CEScore was negatively correlated with OXPHOS and significantly featured prognosis in ccRCC. The ccRCC patients with high-risk scores show worse survival outcomes and bad clinical benefits of Everolimus (mTOR inhibitor).

CONCLUSIONS

Our findings indicate the importance of abnormal CRGs expression in cancers. In addition, identified several prognostic CRGs as potential markers for prognostic distinction and drug response in the specific tumor. These results accelerate the understanding of copper-induced death in tumor progression and provide cuproptosis-associated novel therapeutic strategies.

Salivary DNA methylation markers for cancer of oral cavity

PURPOSE: Aberrant DNA methylation plays a crucial role in oral carcinogenesis. Our previous study demonstrated hypermethylation of DAPK1, LRPPRC, RAB6C, and ZNF471 promoters in patients with tongue squamous cell carcinoma compared with normal samples. Methylation profiling using salivary DNA is considered a non-invasive alternative to tissue samples. Hence, the present study tested the DNA methylation status of these four promoters as indicators of oral cancer progression. METHODS: We performed the bisulfite-based targeted next-generation sequencing of four candidate genes in saliva and tissue DNA from normal, premalignant, and squamous cell carcinoma subjects. The clinicopathological association, diagnostic, and prognostic utility of aberrant DNA methylation were evaluated using the TCGA-HNSCC dataset. Using the Xgboost algorithm and logistic regression, CpG sites were prioritized, and Receiver Operating Characteristic was generated. By Log-rank test and Kaplan-Meier (KM) curves, an association between methylation and overall survival (OS), disease-free interval (DFI), and progression-free interval (PFI) were computed. RESULTS: We identified all four genes as significantly hypermethylation in premalignant and malignant samples compared with normal samples. The methylation levels were comparable between saliva and tissue samples with an r-value of 0.6297 to 0.8023 and 0.7823 to 0.9419 between premalignant tissue vs. saliva and OC vs. saliva, respectively. We identified an inverse correlation between DAPK1, LRPPRC, RAB6C, and ZNF471 promoter methylation with its expression. A classifier of 8 differentially methylated CpG sites belonging to DAPK1, RAB6C, and ZNF471 promoters was constructed, showing an AUC of 0.984 to differentiate tumors from normal samples. The differential methylation status of DAPK1, LRPPRC, and ZNF71 promoters was prognostically important. Abnormal expression of all four genes was associated with immune infiltration. CONCLUSIONS: Thus, methylation analysis of these candidate CpG sites from saliva can be helpful as a non-invasive tool for the clinical management of OC.

A critical review of datasets and computational suites for improving cancer theranostics and biomarker discovery

Cancer has been constantly evolving and so is the research pertaining to cancer diagnosis and therapeutic regimens. Early detection and specific therapeutics are the key features of modern cancer therapy. These requirements can only be fulfilled with the integration of diverse high-throughput technologies. Integration of advanced omics methodology involving genomics, epigenomics, proteomics, and transcriptomics provide a clear understanding of multi-faceted cancer. In the past few years, tremendous high-throughput data have been generated from cancer genomics and epigenomic analyses, which on further methodological analyses can yield better biological insights. The major epigenetic alterations reported in cancer are DNA methylation levels, histone post-translational modifications, and epi-miRNA regulating the oncogenes and tumor suppressor genes. While the genomic analyses like gene expression profiling, cancer gene prediction, and genome annotation divulge the genetic alterations in oncogenes or tumor suppressor genes. Also, systems biology approach using biological networks is being extensively used to identify novel cancer biomarkers. Therefore, integration of these multi-dimensional approaches will help to identify potential diagnostic and therapeutic biomarkers. Here, we reviewed the critical databases and tools dedicated to various epigenomic and genomic alterations in cancer. The review further focuses on the multi-omics resources available for further validating the identified cancer biomarkers. We also highlighted the tools for cancer biomarker discovery using a systems biology approach utilizing genomic and epigenomic data. Biomarkers predicted using such integrative approaches are shown to be more clinically relevant.

SLC22A8: An indicator for tumor immune microenvironment and prognosis of ccRCC from a comprehensive analysis of bioinformatics

Clear cell renal cell carcinoma (ccRCC) is one of the most common renal malignancies worldwide. SLC22A8 plays a key role in renal excretion of organic anions. However, its role in ccRCC remains unclear; therefore, this study aimed to elucidate the relationship between SLC22A8 and ccRCC. The The Cancer Genome Atlas-kidney renal clear cell carcinoma cohort was included in this study. The Wilcoxon signed-rank test and logistic regression were used to analyze the relationship between SLC22A8 expression and clinicopathological characteristics. Multifactorial analysis and Kaplan-Meier survival curves were adopted for correlation between SLC22A8 expression and clinicopathological parameters and overall survival. Utilizing the UALCAN database, the correlation of the expression levels of SLC22A8 DNA methylation in ccRCC was explored. Immunological characterization of SLC22A8 regarding the ccRCC tumor microenvironment was carried out by the single sample Gene Set Enrichment Analysis algorithm and the CIBERSORT algorithm. With the CellMiner database, the analysis of the association between SLC22A8 gene expression and drug sensitivity was further performed. Eventually, gene ontology and Kyoto Encyclopedia of Gene and Genome enrichment analyses were applied to identify the functional and signaling pathways involved in SLC22A8. SLC22A8 expression is associated with age, grade, stage, and tumor status. SLC22A8 protein expression levels, phosphorylated protein levels, and DNA methylation expression levels were lower in ccRCC tissues than in normal tissues, and low methylation levels predicted poor overall survival. Comprehensive analysis of tumor immune infiltration and the tumor microenvironment indicated a higher level of overall immunity in the SLC22A8 low expression group. Gene Enrichment Analysis results showed that low expression of SLC22A8 was associated with immune pathways, such as phagocytosis recognition and humoral immune response. SLC22A8 expression was significantly correlated with survival and immune infiltration in ccRCC and can be used as a prognostic biomarker for ccRCC.

The regulation of ISG20 expression on SARS-CoV-2 infection in cancer patients and healthy individuals

ISG20 inhibits viruses such as SARS-CoV-2 invasion; however, details of its expression and regulation with viral susceptibility remain to be elucidated. The present study analyzed ISG20 expression, isoform information, survival rate, methylation patterns, immune cell infiltration, and COVID-19 outcomes in healthy and cancerous individuals. Cordycepin (CD) and N6, N6-dimethyladenosine (m⁶₂A) were used to treat cancer cells for ISG20 expression. We revealed that ISG20 mRNA expression was primarily located in the bone marrow and lymphoid tissues. Interestingly, its expression was significantly increased in 11 different types of cancer, indicating that cancer patients may be less vulnerable to SARS-CoV-2 infection. Among them, higher expression of ISG20 was associated with a long OS in CESC and SKCM, suggesting that ISG20 may be a good marker for both viral prevention and cancer progress. ISG20 promoter methylation was significantly lower in BLCA, READ, and THCA tumor tissues than in the matched normal tissues, while higher in BRCA, LUSC, KIRC, and PAAD. Hypermethylation of ISG20 in KIRC and PAAD tumor tissues was correlated with higher expression of ISG20, suggesting that methylation of ISG20 may not underlie its overexpression. Furthermore, ISG20 expression was significantly correlated with immune infiltration levels, including immune lymphocytes, chemokine, receptors, immunoinhibitors, immunostimulators, and MHC molecules in pan-cancer. STAD exhibited the highest degree of ISG20 mutations; the median progression-free survival time in months for the unaltered group was 61.84, while it was 81.01 in the mutant group. Isoforms ISG20-001 and ISG20−009 showed the same RNase_T domain structure, demonstrating the functional roles in tumorigenesis and SARS-CoV-2 invasion inhibition in cancer patients. Moreover, CD and m⁶₂A increase ISG20 expression in various cancer cell lines, implying the antiviral/anti-SARS-CoV-2 therapeutic potential. Altogether, this study highlighted the value of combating cancer by targeting ISG20 during the COVID-19 pandemic, and small molecules extracted from traditional Chinese medicines, such as CD, may have potential as anti-SARS-CoV-2 and anticancer agents by promoting ISG20 expression.

Hypoxia-Related Signature Is a Prognostic Biomarker of Pancreatic Cancer

Background

Hypoxia plays a significant role in the pathogenesis of pancreatic cancer, but the effect of hypoxia-related genes in pancreatic cancer remains to be elucidated. This study aimed to identify hypoxia-related genes related to pancreatic cancer and construct a prognostic signature.

Methods

Pancreatic cancer datasets were retrieved from TCGA database. Cox regression analyses were used to identify hypoxia-related genes and construct a prognostic signature. Datasets from International Cancer Genome Consortium and GEO databases were used as validated cohorts. The CIBERSORT method was applied to estimate the fractions of immune cell types. DNA methylation and protein levels of the genes in pancreatic cancer were examined.

Results

Three hypoxia-related genes (TES, LDHA, and ANXA2) were identified as associated with patient survival and selected to construct a prognostic signature. Patients were divided into high- and low-risk groups based on the signature. Those in the high-risk group showed worse survival than those in the low-risk group. The signature was shown to be involved in the HIF-1 signaling pathway. The time-dependent ROC analyses of three independent validated cohorts further revealed that this signature had a better prognostic value in the prediction of the survival of pancreatic cancer patients. Immune cells analysis for three datasets demonstrated that high-risk signature was significantly associated with macrophages and T cells. DNA methylation and protein levels of the three genes validated their aberrant expression in pancreatic cancer.

Conclusions

Our research provided a novel and reliable prognostic signature that composes of three hypoxia-related genes to estimate the prognosis of pancreatic cancer.

The Expression and Role Analysis of Methylation-Regulated Differentially Expressed Gene UBE2C in Pan-Cancer, Especially for HGSOC

Simple Summary

DNA methylation has attracted a great deal of scientific interest as an early biomarker and potential therapeutic target. HGSOC result in high mortality due to the absence of reliable biomarkers for early diagnosis and prognosis. In this study, we performed an integrated bioinformatic analysis and found that UBE2C was hypomethylation and overexpression in ovarian cancer, which was associated with advanced cancer stages and poor prognoses. Meantime, this finding was also confirmed in pan-cancer analysis. Furthermore, the experimental validation of the expression and role of UBE2C was performed on HGSOC tissues and cancer cell lines. Importantly, demethylation could upregulate the expression of UBE2C. Taken together, methylation-regulated UBE2C may be a novel biomarker for diagnosis and prognosis, not only for ovarian cancer but a variety of cancers.

Abstract

High-grade serous ovarian cancer (HGSOC) is the most fatal gynecological malignant tumor. DNA methylation is associated with the occurrence and development of a variety of tumor types, including HGSOC. However, the signatures regarding DNA methylation changes for HGSOC diagnosis and prognosis are less explored. Here, we screened differentially methylated genes and differentially expressed genes in HGSOC through the GEO database. We identified that UBE2C was hypomethylation and overexpression in ovarian cancer, which was associated with more advanced cancer stages and poor prognoses. Additionally, the pan-cancer analysis showed that UBE2C was overexpressed and hypomethylation in almost all cancer types and was related to poor prognoses for various cancers. Next, we established a risk or prognosis model related to UBE2C methylation sites and screened out the three sites (cg03969725, cg02838589, and cg00242976). Furthermore, we experimentally validated the overexpression of UBE2C in HGSOC clinical samples and ovarian cell lines using quantitative real-time PCR, Western blot, and immunohistochemistry. Importantly, we discovered that ovarian cancer cell lines had lower DNA methylation levels of UBE2C than IOSE-80 cells (normal ovarian epithelial cell line) by bisulfite sequencing PCR. Consistently, treatment with 5-Azacytidine (a methylation inhibitor) was able to restore the expression of UBE2C. Taken together, our study may help us to understand the underlying molecular mechanism of UBE2C in pan-cancer tumorigenesis; it may be a useful biomarker for diagnosis, treatment, and monitoring, not only of ovarian cancer but a variety of cancers.

TGFβ2 is a Prognostic Biomarker for Gastric Cancer and is Associated With Methylation and Immunotherapy Responses

TGFβ signaling plays a key role in cancer progression and by shaping tumor architecture and inhibiting the anti-tumor activity of immune cells. It was reported that high expression of TGFβ can promote the invasion and metastasis of cancer cells in a variety of tumors. However, there are few studies on TGFβ2 and its methylation in gastric cancer. We analyzed the Harbin Medical University Cancer Hospital (HMUCH) sequencing data and used public data to explore the potential function and prognostic value of TGFβ2 and its methylation in gastric cancer. In this study, we used the ssGSEA algorithm to quantify 23 methylation sites related to TGFβ2. Survival analysis showed that high expression of TGFβ2 and hypomethylation levels of TGFβ2 were negative factors in the prognosis of gastric cancer. Functional enrichment analysis of methylation revealed that methylation of different TGFβ2 methylation scores was mainly involved in energy metabolism, extracellular matrix formation and cell cycle regulation. In the gastric cancer microenvironment TGFβ2 was associated with high levels of multiple immune cell infiltration and cytokine expression, and high TGFβ2 expression was significantly and positively correlated with stemness markers, stromalscore and EMT. Gene set enrichment analysis also revealed an important role of TGFβ2 in promoting EMT. In addition, we discussed the relationship between TGFβ2 and immunotherapy. The expression of PD-1, PD-L1 and CTLA-4 was elevated in the TGFβ2 high expression group. Also when TGFβ2 was highly expressed, the responsiveness of immune checkpoint blockade (ICB) was significantly enhanced. This indicates that TGFβ2 may become an indicator for predicting the efficacy of immunosuppressive agents and a potential target for immunotherapy.

Promoter hypermethylation of GALR1 acts as an early epigenetic susceptibility event in colorectal carcinogenesis

Epigenetics play an essential role in colorectal neoplasia process. There is a need to determine the appropriateness of epigenetic biomarkers for early detection as well as expand our understanding of the carcinogenic process. Therefore, the aim of the study was to assess how DNA methylation pattern of GALR1 gene evolves in a sample set representing colorectal neoplastic progression. The study was designed into three phases. Firstly, Methylation status of GALR1 was assessed with genome-wide DNA methylation beadchip and pyrosequencing assays in colorectal lesions and paired normal tissues. Then, linear mixed-effects modeling analyses were applied to describe the trend of DNA methylation during the progression of colorectal neoplasia. In the third phase, quantitative RT-PCR was used to examine GALR1 expression in patients with precursor lesion and colorectal cancer. We found that significant hypermethylation of GALR1 promoter was a widely existent modification in CRCs (P < 0.001). When further examined methylation pattern of GALR1 during neoplastic progression of CRC, we found that DNA methylation level of GALR1 showed a significant stepwise increase from normal to hyperplastic polyps, to adenomas and to carcinoma samples (P < 0.001). Besides, loss of mRNA expression is a common accompaniment to adenomas and carcinomas. Public omics data analyses showed an inverse correlation between gene expression and DNA methylation (P < 0.001). Our findings indicate that epigenetic alteration of GALR1 promoter is gradually accumulated during the colorectal neoplastic progression. It can potentially be a promising biomarker used for screening and surveillance of colorectal cancer.

Integrated computational analysis reveals HOX genes cluster as oncogenic drivers in head and neck squamous cell carcinoma

Alterations in homeobox (HOX) gene expression are involved in the progression of several cancer types including head and neck squamous cell carcinoma (HNSCC). However, regulation of the entire HOX cluster in the pathophysiology of HNSCC is still elusive. By using different comprehensive databases, we have identified the significance of differentially expressed HOX genes (DEHGs) in stage stratification and HPV status in the cancer genome atlas (TCGA)-HNSCC datasets. The genetic and epigenetic alterations, druggable genes, their associated functional pathways and their possible association with cancer hallmarks were identified. We have performed extensive analysis to identify the target genes of DEHGs driving HNSCC. The differentially expressed HOX cluster-embedded microRNAs (DEHMs) in HNSCC and their association with HOX-target genes were evaluated to construct a regulatory network of the HOX cluster in HNSCC. Our analysis identified sixteen DEHGs in HNSCC and determined their importance in stage stratification and HPV infection. We found a total of 55 HNSCC driver genes that were identified as targets of DEHGs. The involvement of DEHGs and their targets in cancer-associated signaling mechanisms have confirmed their role in pathophysiology. Further, we found that their oncogenic nature could be targeted by using the novel and approved anti-neoplastic drugs in HNSCC. Construction of the regulatory network depicted the interaction between DEHGs, DEHMs and their targets genes in HNSCC. Hence, aberrantly expressed HOX cluster genes function in a coordinated manner to drive HNSCC. It could provide a broad perspective to carry out the experimental investigation, to understand the underlying oncogenic mechanism and allow the discovery of new clinical biomarkers for HNSCC.

A Potential Prognostic Marker PRDM1 in Pancreatic Adenocarcinoma

Pancreatic adenocarcinoma (PAAD) is a major threat to people's health. PRDM1 is a transcription factor with multiple functions, and its functions have been validated in a variety of tumors; however, there are few studies reported on PRDM1 in PAAD. Using the GEPIA2 database, this research found that PRDM1 expression in PAAD was significantly higher than that in normal pancreatic tissue. The Kaplan-Meier Plotter database showed that high expression of PRDM1 in PAAD has a poor prognosis, suggesting that PRDM1 may be a potential prognostic marker in PAAD. The cBioPortal database shows that the expression of PRDM1 in PAAD is significantly correlated with its methylation degree. Further analysis on the coexpressed genes of PRDM1 in PAAD was performed by using LinkedOmics database to explore potential mechanisms. Based on gene enrichment analysis, PRDM1 was implicated in many pathways involved in tumor progression. In the construction of a PPI network of PRDM1 and its coexpressed gene protein via the STRING database, we found that PRDM1 may be involved in the pathogenesis and development of PAAD. TIMER database suggested that a high level of PRDM1 has a significant positive correlation with macrophages, neutrophils, and DCs. Potential methylation sites of PRDM1 were found through DNMIVD database, and MethSurv database explored eight sites which were significantly related with the prognosis of PAAD. In conclusion, PRDM1 may work as a prognostic marker or even provide a potential therapeutic strategy in PAAD.

MrgprF acts as a tumor suppressor in cutaneous melanoma by restraining PI3K/Akt signaling

The incidence of cutaneous melanoma (CM) has been increasing annually worldwide. In this study, we identify that MrgprF, a MAS related GPR family member, is decreased in cutaneous melanoma tissues and cell lines due to hypermethylation of its promoter region, and show that patients with CM expressing high levels of MrgprF exhibit an improved clinical outcome. We demonstrate that MrgprF forced expression inhibits tumor cell proliferation, migration, xenograft tumor growth, and metastasis. On the contrary, MrgprF knockdown promotes tumor cell proliferation and transformation of immortalized human keratinocyte-HaCaT cells, supporting the inhibitory role of MrgprF during tumor progression. Mechanistic studies reveal that MrgprF reduces the phosphoinositol‑3‑kinase (PI3K) complex formation between p101 and p110γ subunits, the critical step for phosphatidylinositol-(3, 4)-P2 (PIP2) conversion to phosphatidylinositol-(3, 4, 5)-P3 (PIP3), and then reduces the activation of PI3K/Akt signaling. This effect can be reversed by Akt specific agonist SC79. In addition, AMG 706, a previously documented inhibitor for endothelial cell proliferation, is identified as a potential agonist for MrgprF, and can impede tumor growth both in vitro and in vivo. Taken together, our findings suggest that MrgprF, a novel tumor suppressor in cutaneous melanoma, may be useful as a therapeutic target in the future.

MIAOME: Human Microbiome Affect The Host Epigenome

Besides the genetic factors having tremendous influences on the regulations of the epigenome, the microenvironmental factors have recently gained extensive attention for their roles in affecting the host epigenome. There are three major types of microenvironmental factors: microbiota-derived metabolites (MDM), microbiota-derived components (MDC) and microbiota-secreted proteins (MSP). These factors can regulate host physiology by modifying host gene expression through the three highly interconnected epigenetic mechanisms (e.g. histone modifications, DNA modifications, and non-coding RNAs). However, no database was available to provide the comprehensive factors of these types. Herein, a database entitled 'Human Microbiome Affect The Host Epigenome (MIAOME)' was constructed. Based on the types of epigenetic modifications confirmed in the literature review, the MIAOME database captures 1068 (63 genus, 281 species, 707 strains, etc.) human microbes, 91 unique microbiota-derived metabolites & components (16 fatty acids, 10 bile acids, 10 phenolic compounds, 10 vitamins, 9 tryptophan metabolites, etc.) derived from 967 microbes; 50 microbes that secreted 40 proteins; 98 microbes that directly influence the host epigenetic modification, and provides 3 classifications of the epigenome, including (1) 4 types of DNA modifications, (2) 20 histone modifications and (3) 490 ncRNAs regulations, involved in 160 human diseases. All in all, MIAOME has compiled the information on the microenvironmental factors influence host epigenome through the scientific literature and biochemical databases, and allows the collective considerations among the different types of factors. It can be freely assessed without login requirement by all users at: http://miaome.idrblab.net/ttd/

Network Approaches for Charting the Transcriptomic and Epigenetic Landscape of the Developmental Origins of Health and Disease

The early developmental phase is of critical importance for human health and disease later in life. To decipher the molecular mechanisms at play, current biomedical research is increasingly relying on large quantities of diverse omics data. The integration and interpretation of the different datasets pose a critical challenge towards the holistic understanding of the complex biological processes that are involved in early development. In this review, we outline the major transcriptomic and epigenetic processes and the respective datasets that are most relevant for studying the periconceptional period. We cover both basic data processing and analysis steps, as well as more advanced data integration methods. A particular focus is given to network-based methods. Finally, we review the medical applications of such integrative analyses.