Genome-wide profiling of DNA methylation and gene expression in esophageal squamous cell carcinoma

Investigation of the Mechanism of Periploca forrestii against Rheumatoid Arthritis with Network Pharmacology-Based Analysis

Periploca forrestii Schltr (P. forrestii) is an edible medicinal herb with various health benefits such as treating antirheumatoid arthritis (RA), reducing inflammation, and preventing tumor growth. The active ingredients in P. forrestii responsible for its protective effect against RA, however, remain unknown. In this study, the active ingredient of P. forrestii and its potential mechanism of action against RA were investigated by network pharmacology and enrichment analysis. The methods included predicting target genes of P. forrestii, constructing a protein interaction network, and performing gene-ontology (GO) and Kyoto-encyclopedia of genes and genomes (KEGG) enrichment analysis. We discovered targets of RA through retrieval of OMIM and GeneCards public databases. Cardiac glycosides (CGs) are considered the primarily active ingredients of P. forrestii, and the target genes of GCs were discovered to be overlapped with relevant targets of RA using the Venn diagram. After that, prediction of relevant targets of P. forrestii was accomplished with a network pharmacology-based approach. Through the Venn diagram, we discovered 99 genes shared in the target genes of P. forrestii and RA. Gene enrichment analysis showed that the mechanisms of CGs against RA are associated with 55 signaling pathways, including endocrine resistance, Epstein-Barr virus infection, bladder cancer, prostate cancer, and coronavirus disease (COVID-19) signaling pathways. Coexpression analysis indicated ADSL, ATIC, AR, CCND1, MDM2, and HSP90AA1 as the hub genes between putative targets of P. forrestii-derived CGs and known therapeutic targets of RA. In conclusion, we clarified the mechanism of action of P. forrestii against RA, which would provide a basis for further understanding the clinical application of P. forrestii.

Multi-omic characterization of genome-wide abnormal DNA methylation reveals diagnostic and prognostic markers for esophageal squamous-cell carcinoma

This study investigates aberrant DNA methylations as potential diagnosis and prognosis markers for esophageal squamous-cell carcinoma (ESCC), which if diagnosed at advanced stages has <30% five-year survival rate. Comparing genome-wide methylation sites of 91 ESCC and matched adjacent normal tissues, we identified 35,577 differentially methylated CpG sites (DMCs) and characterized their distribution patterns. Integrating whole-genome DNA and RNA-sequencing data of the same samples, we found multiple dysregulated transcription factors and ESCC-specific genomic correlates of identified DMCs. Using featured DMCs, we developed a 12-marker diagnostic panel with high accuracy in our dataset and the TCGA ESCC dataset, and a 4-marker prognostic panel distinguishing high-risk patients. In-vitro experiments validated the functions of 4 marker host genes. Together these results provide additional evidence for the important roles of aberrant DNA methylations in ESCC development and progression. Our DMC-based diagnostic and prognostic panels have potential values for clinical care of ESCC, laying foundations for developing targeted methylation assays for future non-invasive cancer detection methods.

Association of dietary intakes of vitamin B12, vitamin B6, folate, and methionine with the risk of esophageal cancer: the Japan Public Health Center-based (JPHC) prospective study

Background

B vitamins and methionine are essential substrates in the one-carbon metabolism pathway involved in DNA synthesis and methylation. They may have essential roles in cancer development. We aimed to evaluate the associations of dietary intakes of vitamin B12, vitamin B6, folate, and methionine with the risk of esophageal cancer (EC) using data from the Japan Public Health Center-based Prospective Study.

Methods

We included 87,053 Japanese individuals who completed a food frequency questionnaire and were followed up from 1995–1998 to 2013 and 2015. Hazard ratios (HRs) and 95% confidence intervals (CIs) were calculated by Cox proportional-hazard regression across quintiles of dietary intakes of B vitamins and methionine.

Results

After 1,456,678 person-years of follow-up, 427 EC cases were documented. The multivariable HR (95% CI) of incident EC in the highest versus lowest quintile of dietary intake of vitamin B12 was 1.75 (1.13–2.71; p-trend=0.01). Stratification analysis based on alcohol consumption showed that higher dietary intakes of vitamin B12 and methionine were associated with an increased risk of EC among never-drinkers; HRs (95% CIs) were 2.82 (1.18–6.74; p-trend=0.009; p-interaction=0.18) and 3.45 (1.32–9.06; p-trend=0.003; p-interaction 0.02) for vitamin B12 and methionine, respectively. Meanwhile, there was no association between vitamin B12 and methionine intake with the risk of EC among drinkers. There were no associations between dietary intake of folate or vitamin B6 and the risk of EC.

Conclusion

Dietary intake of vitamin B12 was positively associated with the risk of EC in the Japanese population.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08721-8.

Genome-Wide DNA Methylation Profiling of Esophageal Squamous Cell Carcinoma from Global High-Incidence Regions Identifies Crucial Genes and Potential Cancer Markers

Epigenetic mechanisms such as aberrant DNA methylation (DNAme) are known to drive esophageal squamous cell carcinoma (ESCC), yet they remain poorly understood. Here, we studied tumor-specific DNAme in ESCC cases from nine high-incidence countries of Africa, Asia, and South America. Infinium MethylationEPIC array was performed on 108 tumors and 51 normal tissues adjacent to the tumors (NAT) in the discovery phase, and targeted pyrosequencing was performed on 132 tumors and 36 NAT in the replication phase. Top genes for replication were prioritized by weighting methylation results using RNA-sequencing data from The Cancer Genome Atlas and GTEx and validated by qPCR. Methylome analysis comparing tumor and NAT identified 6,796 differentially methylated positions (DMP) and 866 differential methylated regions (DMR), with a 30% methylation (Δβ) difference. The majority of identified DMPs and DMRs were hypermethylated in tumors, particularly in promoters and gene-body regions of genes involved in transcription activation. The top three prioritized genes for replication, PAX9, SIM2, and THSD4, had similar methylation differences in the discovery and replication sets. These genes were exclusively expressed in normal esophageal tissues in GTEx and downregulated in tumors. The specificity and sensitivity of these DNAme events in discriminating tumors from NAT were assessed. Our study identified novel, robust, and crucial tumor-specific DNAme events in ESCC tumors across several high-incidence populations of the world. Methylome changes identified in this study may serve as potential targets for biomarker discovery and warrant further functional characterization. SIGNIFICANCE: This largest genome-wide DNA methylation study on ESCC from high-incidence populations of the world identifies functionally relevant and robust DNAme events that could serve as potential tumor-specific markers. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/81/10/2612/F1.large.jpg.

Identification and Validation of Potential Pathogenic Genes and Prognostic Markers in ESCC by Integrated Bioinformatics Analysis

Esophageal squamous cell carcinoma (ESCC) is one of the most fatal malignancies of the digestive tract, but its underlying molecular mechanisms are not known. We aim to identify the genes involved in ESCC carcinogenesis and discover potential prognostic markers using integrated bioinformatics analysis. Three pairs of ESCC tissues and paired normal tissues were sequenced by high-throughput RNA sequencing (RNA-seq). Integrated bioinformatics analysis was used to identify differentially expressed coding genes (DECGs) and differentially expressed long non-coding RNA (lncRNA) genes (DELGs). A protein-protein interaction (PPI) network of DECGs was established using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) website and visualized with Cytoscape. Survival analysis was conducted by log-rank tests to identify "hub" genes with potential prognostic value, and real-time reverse transcription-quantitative polymerase chain reaction (RT-qPCR) was conducted to assess expression of these genes in ESCC tissues. Transwell^TM assays were employed to examine the migration ability of cells after knockdown of LINC01614 expression, followed by investigation of epithelial-mesenchymal transition (EMT) by western blotting (WB). A total of 106 upregulated genes and 42 downregulated genes were screened out from the ESCC data sets. Survival analysis showed two hub protein-coding genes with higher expression in module 1 of the PPI network (SPP1 and BGN) and another three upregulated lncRNAs (LINC01614, LINC01415, NKILA) that were associated with a poor prognosis. High expression of SPP1, BGN, LINC01614, and LINC01415 in tumor samples was validated further by RT-qPCR. In vitro experiments show that knockdown of LINC01614 expression could significantly inhibit the migration of ESCC cells by regulating EMT, which was confirmed by WB. These results indicate that BGN, SPP1, LINC01614, and LINC01415 might be critical genes in ESCC and potential prognostic biomarkers.

DNA Methylome Distinguishes Head and Neck Cancer from Potentially Malignant Oral Lesions and Healthy Oral Mucosa

There is a strong need to find new, good biomarkers of head and neck squamous cell carcinoma (HNSCC) because of the bad prognoses and high mortality rates. The aim of this study was to identify the potential biomarkers in HNSCC that have differences in their DNA methylome and potentially premalignant oral lesions, in comparison to healthy oral mucosa. In this study, 32 oral samples were tested: nine healthy oral mucosae, 13 HNSCC, and 10 oral lesions for DNA methylation by the Infinium MethylationEPIC BeadChip. Our findings showed that a panel of genes significantly hypermethylated in their promoters or specific sites in HNSCC samples in comparison to healthy oral samples, which are mainly oncogenes, receptor, and transcription factor genes, or genes included in cell cycle, transformation, apoptosis, and autophagy. A group of hypomethylated genes in HNSCC, in comparison to healthy oral mucosa, are mainly involved in the host immune response and transcriptional regulation. The results also showed significant differences in gene methylation between HNSCC and potentially premalignant oral lesions, as well as differently methylated genes that discriminate between oral lesions and healthy mucosa. The given methylation panels point to novel potential biomarkers for early diagnostics of HNSCC, as well as potentially premalignant oral lesions.

Methylation silencing of TGF-β receptor type II is involved in malignant transformation of esophageal squamous cell carcinoma

Background

Although massive studies have been conducted to investigate the mechanisms of esophageal squamous cell carcinoma (ESCC) carcinogenesis, the understanding of molecular alterations during the malignant transformation of epithelial dysplasia is still lacking, especially regarding epigenetic changes.

Results

To better characterize the methylation changes during the malignant transformation of epithelial dysplasia, a whole-genome bisulfite sequencing analysis was performed on a series of tumor, dysplastic, and non-neoplastic epithelial tissue samples from esophageal squamous cell carcinoma (ESCC) patients. Promoter hypermethylation in TGF-β receptor type II (TGFBR2), an important mediator of TGF-β signaling, was identified. Further, we evaluated the methylation and expression of TGFBR2 in tumor samples through The Cancer Genome Atlas multiplatform data as well as immunohistochemistry. Moreover, treatment of ESCC cell lines with5-Aza-2′-deoxycytidine, a DNA methyltransferase inhibitor, reactivated the expression of TGFBR2. The lentiviral mediating the overexpression of TGFBR2 inhibited the proliferation of ESCC cell line by inducing cell cycle G2/M arrest. Furthermore, the overexpression of TGFBR2 inhibited the tumor growth obviously in vivo.

Conclusions

The characterization of methylation silencing of TGFBR2 in ESCC will enable us to further explore whether this epigenetic change could be considered as a predictor of malignant transformation in esophageal epithelial dysplasia and whether use of a TGFBR2 agonist may lead to a new therapeutic strategy in patients with ESCC.

Methylation silencing of TGF-β receptor type II is involved in malignant transformation of esophageal squamous cell carcinoma

BACKGROUND

Although massive studies have been conducted to investigate the mechanisms of esophageal squamous cell carcinoma (ESCC) carcinogenesis, the understanding of molecular alterations during the malignant transformation of epithelial dysplasia is still lacking, especially regarding epigenetic changes.

RESULTS

To better characterize the methylation changes during the malignant transformation of epithelial dysplasia, a whole-genome bisulfite sequencing analysis was performed on a series of tumor, dysplastic, and non-neoplastic epithelial tissue samples from esophageal squamous cell carcinoma (ESCC) patients. Promoter hypermethylation in TGF-β receptor type II (TGFBR2), an important mediator of TGF-β signaling, was identified. Further, we evaluated the methylation and expression of TGFBR2 in tumor samples through The Cancer Genome Atlas multiplatform data as well as immunohistochemistry. Moreover, treatment of ESCC cell lines with5-Aza-2'-deoxycytidine, a DNA methyltransferase inhibitor, reactivated the expression of TGFBR2. The lentiviral mediating the overexpression of TGFBR2 inhibited the proliferation of ESCC cell line by inducing cell cycle G2/M arrest. Furthermore, the overexpression of TGFBR2 inhibited the tumor growth obviously in vivo.

CONCLUSIONS

The characterization of methylation silencing of TGFBR2 in ESCC will enable us to further explore whether this epigenetic change could be considered as a predictor of malignant transformation in esophageal epithelial dysplasia and whether use of a TGFBR2 agonist may lead to a new therapeutic strategy in patients with ESCC.

Cancer mutational burden is shaped by G4 DNA, replication stress and mitochondrial dysfunction

A hallmark of cancer is genomic instability, which can enable cancer cells to evade therapeutic strategies. Here we employed a computational approach to uncover mechanisms underlying cancer mutational burden by focusing upon relationships between 1) translocation breakpoints and the thousands of G4 DNA-forming sequences within retrotransposons impacting transcription and exemplifying probable non-B DNA structures and 2) transcriptome profiling and cancer mutations. We determined the location and number of G4 DNA-forming sequences in the Genome Reference Consortium Human Build 38 and found a total of 358,605 covering ∼13.4 million bases. By analyzing >97,000 unique translocation breakpoints from the Catalogue Of Somatic Mutations In Cancer (COSMIC), we found that breakpoints are overrepresented at G4 DNA-forming sequences within hominid-specific SVA retrotransposons, and generally occur in tumors with mutations in tumor suppressor genes, such as TP53. Furthermore, correlation analyses between mRNA levels and exome mutational loads from The Cancer Genome Atlas (TCGA) encompassing >450,000 gene-mutation regressions revealed strong positive and negative associations, which depended upon tissue of origin. The strongest positive correlations originated from genes not listed as cancer genes in COSMIC; yet, these show strong predictive power for survival in most tumor types by Kaplan-Meier estimation. Thus, correlation analyses of DNA structure and gene expression with mutation loads complement and extend more traditional approaches to elucidate processes shaping genomic instability in cancer. The combined results point to G4 DNA, activation of cell cycle/DNA repair pathways, and mitochondrial dysfunction as three major factors driving the accumulation of somatic mutations in cancer cells.

Overexpression of long noncoding RNA LINC01419 in esophageal squamous cell carcinoma and its relation to the sensitivity to 5-fluorouracil by mediating GSTP1 methylation

BACKGROUND

Genome-wide sequencing investigations have identified numerous long noncoding RNAs (lncRNAs) among mammals, many of which exhibit aberrant expression in cancers, including esophageal squamous cell carcinoma (ESCC). Herein, this study elucidates the role and mechanism by which LINC01419 regulates the DNA methylation of glutathione S-transferase pi 1 (GSTP1) in relation to ESCC progression and the sensitivity of ESCC cells to 5-fluorouracil (5-FU).

METHODS

LINC01419 and GSTP1 levels were quantified among 38 paired ESCC and adjacent tissue samples collected from patients with ESCC. To ascertain the contributory role of LINC01419 in the progression of ESCC and identify the interaction between LINC01419 and GSTP1 promoter methylation, LINC01419 was overexpressed or silenced, and the DNA methyltransferase inhibitor 5-Aza-CdR was treated.

RESULTS

Data from the GEO database (GSE21362) and the Cancer Genome Atlas displayed elevated levels of LINC01419 and downregulated levels of GSTP1 in the ESCC tissues and cells. The silencing of LINC01419 led to decreased proliferation, increased apoptosis, and enhanced sensitivity to 5-FU in ESCC cells. Notably, LINC01419 could bind to the promoter region of the GSTP1 gene, resulting in elevated GSTP1 methylation and reduced GSTP1 levels via the recruitment of DNA methyltransferase among ESCC cells, whereby ESCC progression was stimulated accompanied by reduced ESCC cell sensitivity to 5-FU. GSTP1 demethylation by 5-Aza-CdR was observed to reverse the effects of LINC01419 overexpression in ESCC cells and the response to 5-FU.

CONCLUSION

Highly expressed LINC01419 in ESCC promotes GSTP1 methylation, which ultimately acts to promote the event of ESCC and diminish the sensitivity of ESCC cells to 5-FU, highlighting a novel potential strategy to improve 5-FU-based chemotherapy in ESCC.

Molecular landscape of esophageal cancer: implications for early detection and personalized therapy

Esophageal cancer (EC) is one of the most lethal cancers and a public health concern worldwide, owing to late diagnosis and lack of efficient treatment. Esophageal squamous cell carcinoma (ESCC) and esophageal adenocarcinoma (EAC) are main histopathological subtypes of EC that show striking differences in geographical distribution, possibly due to differences in exposure to risk factors and lifestyles. ESCC and EAC are distinct diseases in terms of cell of origin, epidemiology, and molecular architecture of tumor cells. Past efforts aimed at translating potential molecular candidates into clinical practice proved to be challenging, underscoring the need for identifying novel candidates for early diagnosis and therapy of EC. Several major international efforts have brought about important advances in identifying molecular landscapes of ESCC and EAC toward understanding molecular mechanisms and critical molecular events driving the progression and pathological features of the disease. In our review, we summarize recent advances in the areas of genomics and epigenomics of ESCC and EAC, their mutational signatures and immunotherapy. We also discuss implications of recent advances in characterizing the genome and epigenome of EC for the discovery of diagnostic/prognostic biomarkers and development of new targets for personalized treatment and prevention.

Prognostic and diagnostic potential of isocitrate dehydrogenase 1 in esophageal squamous cell carcinoma

We aimed to investigate the pattern of expression and clinical significance of isocitrate dehydrogenase 1(IDH1) in esophageal squamous cell carcinoma (ESCC). The IDH1 expression was determined by quantitative real-time polymerase chain reaction, immunohistochemistry, and Western blot analysis using 38 pairs of frozen tissues. Enzyme-linked immunosorbent assay was employed to measure 67 pairs of serum samples from patients and their controls to evaluate its diagnostic value. Immunohistochemistry analysis of 111 formalin-fixed paraffin embedded tissue samples was conducted for explaining its prognostic value. After shRNA transfection, CCK8 and clonal efficiency assays were carried on for verifying the function of IDH1 in vitro. Increased expression at mRNA (P < 0.001) and protein levels (immunohistochemistry: P < 0.001, Western blot analysis: P < 0.001) were observed. Similarly, the IDH1 expression in serum from patients with ESCC was significantly upregulated relative to that from healthy controls (P < 0.001). Kaplan–Meier curve indicated that IDH1 upregulation predicted worse overall survival (OS) and progression-free survival (PFS). Univariate and multivariate analyses identified IDH1 expression as an independent prognostic factor for OS and PFS. Furthermore, OD450 values and colony numbers were decreased in sh-IDH1 groups (all P < 0.05). In conclusion, IDH1 is upregulated in patients with ESCC and can be used as a good potential biomarker for diagnosis and prognosis.

International cancer seminars: a focus on esophageal squamous cell carcinoma

The International Agency for Research on Cancer (IARC) and the US National Cancer Institute (NCI) have initiated a series of cancer-focused seminars [Scelo G, Hofmann JN, Banks RE et al. International cancer seminars: a focus on kidney cancer. Ann Oncol 2016; 27(8): 1382-1385]. In this, the second seminar, IARC and NCI convened a workshop in order to examine the state of the current science on esophageal squamous cell carcinoma etiology, genetics, early detection, treatment, and palliation, was reviewed to identify the most critical open research questions. The results of these discussions were summarized by formulating a series of 'difficult questions', which should inform and prioritize future research efforts.

Non-blood circulating tumor DNA detection in cancer

Tumor DNA contains specific somatic alterations that are crucial for the diagnosis and treatment of cancer. Due to the spatial and temporal intra-tumor heterogeneity, multi-sampling is needed to adequately characterize the somatic alterations. Tissue biopsy, however, is limited by the restricted access to sample and the challenges to recapitulate the tumor clonal diversity. Non-blood circulating tumor DNA are tumor DNA fragments presents in non-blood body fluids, such as urine, saliva, sputum, stool, pleural fluid, and cerebrospinal fluid (CSF). Recent studies have demonstrated the presence of tumor DNA in these non-blood body fluids and their application to the diagnosis, screening, and monitoring of cancers. Non-blood circulating tumor DNA has an enormous potential for large-scale screening of local neoplasms because of its non-invasive nature, close proximity to the tumors, easiness and it is an economically viable option. It permits longitudinal assessments and allows sequential monitoring of response and progression. Enrichment of tumor DNA of local cancers in non-blood body fluids may help to archive a higher sensitivity than in plasma ctDNA. The direct contact of cancerous cells and body fluid may facilitate the detection of tumor DNA. Furthermore, normal DNA always dilutes the plasma ctDNA, which may be aggravated by inflammation and injury when very high amounts of normal DNA are released into the circulation. Altogether, our review indicate that non-blood circulating tumor DNA presents an option where the disease can be tracked in a simple and less-invasive manner, allowing for serial sampling informing of the tumor heterogeneity and response to treatment.

The Genomic Impact of DNA CpG Methylation on Gene Expression; Relationships in Prostate Cancer

The process of DNA CpG methylation has been extensively investigated for over 50 years and revealed associations between changing methylation status of CpG islands and gene expression. As a result, DNA CpG methylation is implicated in the control of gene expression in developmental and homeostasis processes, as well as being a cancer-driver mechanism. The development of genome-wide technologies and sophisticated statistical analytical approaches has ushered in an era of widespread analyses, for example in the cancer arena, of the relationships between altered DNA CpG methylation, gene expression, and tumor status. The remarkable increase in the volume of such genomic data, for example, through investigators from the Cancer Genome Atlas (TCGA), has allowed dissection of the relationships between DNA CpG methylation density and distribution, gene expression, and tumor outcome. In this manner, it is now possible to test that the genome-wide correlations are measurable between changes in DNA CpG methylation and gene expression. Perhaps surprisingly is that these associations can only be detected for hundreds, but not thousands, of genes, and the direction of the correlations are both positive and negative. This, perhaps, suggests that CpG methylation events in cancer systems can act as disease drivers but the effects are possibly more restricted than suspected. Additionally, the positive and negative correlations suggest direct and indirect events and an incomplete understanding. Within the prostate cancer TCGA cohort, we examined the relationships between expression of genes that control DNA methylation, known targets of DNA methylation and tumor status. This revealed that genes that control the synthesis of S-adenosyl-l-methionine (SAM) associate with altered expression of DNA methylation targets in a subset of aggressive tumors.