Epigenetic Inactivation of the Tumor Suppressor IRX1 Occurs Frequently in Lung Adenocarcinoma and Its Silencing Is Associated with Impaired Prognosis

Genome-wide DNA methylation in relation to ARID1A deficiency in ovarian clear cell carcinoma

BACKGROUND

The poor chemo-response and high DNA methylation of ovarian clear cell carcinoma (OCCC) have attracted extensive attentions. Recently, we revealed the mutational landscape of the human kinome and additional cancer-related genes and found deleterious mutations in ARID1A, a component of the SWI/SNF chromatin-remodeling complex, in 46% of OCCC patients. The present study aims to comprehensively investigate whether ARID1A loss and genome-wide DNA methylation are co-regulated in OCCC and identify putative therapeutic targets epigenetically regulated by ARID1A.

METHODS

DNA methylation of ARID1Amt/ko and ARID1Awt OCCC tumors and cell lines were analyzed by Infinium MethylationEPIC BeadChip. The clustering of OCCC tumors in relation to clinical and mutational status of tumors were analyzed by hierarchical clustering analysis of genome-wide methylation. GEO expression profiles were used to identify differentially methylated (DM) genes and their expression level in ARID1Amt/ko vs ARID1Awt OCCCs. Combining three pre-ranked GSEAs, pathways and leading-edge genes epigenetically regulated by ARID1A were revealed. The leading-edge genes that passed the in-silico validation and showed consistent ARID1A-related methylation change in tumors and cell lines were regarded as candidate genes and finally verified by bisulfite sequencing and RT-qPCR.

RESULTS

Hierarchical clustering analysis of genome-wide methylation showed two clusters of OCCC tumors. Tumor stage, ARID1A/PIK3CA mutations and TP53 mutations were significantly different between the two clusters. ARID1A mutations in OCCC did not cause global DNA methylation changes but were related to DM promoter or gene-body CpG islands of 2004 genes. Three pre-ranked GSEAs collectively revealed the significant enrichment of EZH2- and H3K27me3-related gene-sets by the ARID1A-related DM genes. 13 Leading-edge DM genes extracted from the enriched gene-sets passed the expression-based in-silico validation and showed consistent ARID1A-related methylation change in tumors and cell lines. Bisulfite sequencing and RT-qPCR analysis showed promoter hypermethylation and lower expression of IRX1, TMEM101 and TRIP6 in ARID1Amt compared to ARID1Awt OCCC cells, which was reversed by 5-aza-2'-deoxycytidine treatment.

CONCLUSIONS

Our study shows that ARID1A loss is related to the differential methylation of a number of genes in OCCC. ARID1A-dependent DM genes have been identified as key genes of many cancer-related pathways that may provide new candidates for OCCC targeted treatment.

Identification and validation of a prognostic anoikis-related gene signature in papillary thyroid carcinoma by integrated analysis of single-cell and bulk RNA-sequencing

Papillary thyroid carcinoma (PTC) prognosis may be deteriorated due to the metastases, and anoikis palys an essential role in the tumor metastasis. However, the potential effect of anoikis-related genes on the prognosis of PTC was unclear. The mRNA and clinical information were obtained from the cancer genome atlas database. Hub genes were identified and risk model was constructed using Cox regression analysis. Kaplan-Meier (K-M) curve was applied for the survival analysis. Immune infiltration and immune therapy response were calculated using CIBERSORT and TIDE. The identification of cell types and cell interaction was performed by Seurat, SingleR and CellChat packages. GO, KEGG, and GSVA were applied for the enrichment analysis. Protein-protein interaction network was constructed in STRING and Cytoscape. Drug sensitivity was assessed in GSCA. Based on bulk RNA data, we identified 4 anoikis-related risk signatures, which were oncogenes, and constructed a risk model. The enrichment analysis found high risk group was enriched in some immune-related pathways. High risk group had higher infiltration of Tregs, higher TIDE score and lower levels of monocytes and CD8 T cells. Based on scRNA data, we found that 4 hub genes were mainly expressed in monocytes and macrophages, and they interacted with T cells. Hub genes were significantly related to immune escape-related genes. Drug sensitivity analysis suggested that cyclin dependent kinase inhibitor 2A may be a better chemotherapy target. We constructed a risk model which could effectively and steadily predict the prognosis of PTC. We inferred that the immune escape may be involved in the development of PTC.

Biomarker RIPK3 Is Silenced by Hypermethylation in Melanoma and Epigenetic Editing Reestablishes Its Tumor Suppressor Function

For several decades, cancers have demonstrably been one of the most frequent causes of death worldwide. In addition to genetic causes, cancer can also be caused by epigenetic gene modifications. Frequently, tumor suppressor genes are epigenetically inactivated due to hypermethylation of their CpG islands, actively contributing to tumorigenesis. Since CpG islands are usually localized near promoters, hypermethylation of the promoter can have a major impact on gene expression. In this study, the potential tumor suppressor gene Receptor Interacting Serine/Threonine Protein Kinase 3 (RIPK3) was examined for an epigenetic regulation and its gene inactivation in melanomas. A hypermethylation of the RIPK3 CpG island was detected by bisulfite pyrosequencing and was accompanied by a correlated loss of its expression. In addition, an increasing RIPK3 methylation rate was observed with increasing tumor stage of melanomas. For further epigenetic characterization of RIPK3, epigenetic modulation was performed using a modified CRISPR/dCas9 (CRISPRa activation) system targeting its DNA hypermethylation. We observed a reduced fitness of melanoma cells by (re-)expression and demethylation of the RIPK3 gene using the epigenetic editing-based method. The tumor suppressive function of RIPK3 was evident by phenotypic determination using fluorescence microscopy, flow cytometry and wound healing assay. Our data highlight the function of RIPK3 as an epigenetically regulated tumor suppressor in melanoma, allowing it to be classified as a biomarker.

ADRA2A and IRX1 are putative risk genes for Raynaud's phenomenon

Raynaud's phenomenon (RP) is a common vasospastic disorder that causes severe pain and ulcers, but despite its high reported heritability, no causal genes have been robustly identified. We conducted a genome-wide association study including 5,147 RP cases and 439,294 controls, based on diagnoses from electronic health records, and identified three unreported genomic regions associated with the risk of RP (p < 5 × 10-8). We prioritized ADRA2A (rs7090046, odds ratio (OR) per allele: 1.26; 95%-CI: 1.20-1.31; p < 9.6 × 10-27) and IRX1 (rs12653958, OR: 1.17; 95%-CI: 1.12-1.22, p < 4.8 × 10-13) as candidate causal genes through integration of gene expression in disease relevant tissues. We further identified a likely causal detrimental effect of low fasting glucose levels on RP risk (rG = -0.21; p-value = 2.3 × 10-3), and systematically highlighted drug repurposing opportunities, like the antidepressant mirtazapine. Our results provide the first robust evidence for a strong genetic contribution to RP and highlight a so far underrated role of α2A-adrenoreceptor signalling, encoded at ADRA2A, as a possible mechanism for hypersensitivity to catecholamine-induced vasospasms.

Innate differences in the molecular signature of normal inferior & superior human parathyroid glands: potential implications for parathyroid adenoma

Primary hyperparathyroidism is a common endocrine disorder. Interestingly, the majority (75%) of parathyroid tumors are localized to the inferior parathyroid glands. To date, the reason for this natural bias has not been investigated. We assessed the global gene expression profile of superior and inferior glands obtained from forensic autopsies. The genes with significant differential expression between superior and inferior parathyroids were further assessed by RT-PCR in 19 pairs. As an iterative approach, additional genes with an established role in parathyroid disorders, i.e., CASR, MAFB, PAX9, TBCE, TBX1, VDR, MEN1, CCND1, and CDC73 were also evaluated by RT-PCR in all 19 pairs of superior and inferior parathyroid glands. Seven homeobox genes, namely HOXA4, HOXA5, HOXBAS3, HOXB4, HOXB6, HOXB9, IRX1, and one encoding for ALDH1A2 showed a lower expression in the inferior parathyroid glands than in the superior. Conversely, SLC6A1 showed a higher expression in the inferior glands. Of the nine genes with significant differential mRNA expression among superior and inferior glands HOXB9, HOXB4 and IRX1 could be detected by western blotting/mass spectrometry. The study is the first to show the differential expression of nine genes HOXA4, HOXA5, HOXBAS3, HOXB4, HOXB6, HOXB9, IRX1, ALDH1A2, and SLC6A1 in inferior versus the superior parathyroid glands. This could have potential implications for the preferential localization of parathyroid tumors to the inferior parathyroid glands as observed in patients with primary hyperparathyroidism.

Correlation between chromatin epigenetic-related lncRNA signature (CELncSig) and prognosis, immune microenvironment, and immunotherapy in non-small cell lung cancer

Chromatin regulators drive cancer epigenetic changes, and lncRNA can play an important role in epigenetic changes as chromatin regulators. We used univariate Cox, LASSO, and multivariate Cox regression analysis to select epigenetic-associated lncRNA signatures. Twenty-five epigenetic-associated lncRNA signatures (CELncSig) were identified to establish the immune prognostic model. According to Kaplan-Meier analysis, the overall survival of the high-risk group was significantly lower than the low-risk group. Receiver operating characteristic (ROC) curves, C-index, survival curve, nomogram, and principal component analysis (PCA) were performed to validate the risk model. In GO/KEGG analysis, differentially expressed lncRNAs were correlated with the PI3K-Akt pathway, suggesting that they were highly associated with the metastasis of LUAD. Interestingly, in the immune escape analysis, the TIDE score was lower, and the possibility of immune dysfunction is also slighter in the high-risk group, which means they still have the potential to receive immunotherapy. And CELncsig is highly correlated with immune pathways T_cell_co-inhibition and Check-point. Also, the IMvigor210 cohort analysis indicated that our risk-scoring model has significant potential clinical application value in lung cancer immunotherapy. And we also screened out ten potential chemotherapy agents using the 'pRRophetic' package.

The Role of IRX Homeobox Genes in Hematopoietic Progenitors and Leukemia

IRX genes are members of the TALE homeobox gene class and encode six related transcription factors (IRX1-IRX6) controlling development and cell differentiation of several tissues in humans. Classification of TALE homeobox gene expression patterns for the hematopoietic compartment, termed TALE-code, has revealed exclusive IRX1 activity in pro-B-cells and megakaryocyte erythroid progenitors (MEPs), highlighting its specific contribution to developmental processes at these early stages of hematopoietic lineage differentiation. Moreover, aberrant expression of IRX homeobox genes IRX1, IRX2, IRX3 and IRX5 has been detected in hematopoietic malignancies, including B-cell precursor acute lymphoblastic leukemia (BCP-ALL), T-cell ALL, and some subtypes of acute myeloid leukemia (AML). Expression analyses of patient samples and experimental studies using cell lines and mouse models have revealed oncogenic functions in cell differentiation arrest and upstream and downstream genes, thus, revealing normal and aberrant regulatory networks. These studies have shown how IRX genes play key roles in the development of both normal blood and immune cells, and hematopoietic malignancies. Understanding their biology serves to illuminate developmental gene regulation in the hematopoietic compartment, and may improve diagnostic classification of leukemias in the clinic and reveal new therapeutic targets and strategies.

An integrative analysis of DNA methylation and gene expression to predict lung adenocarcinoma prognosis

Background: Abnormal DNA methylation of gene promoters is an important feature in lung adenocarcinoma (LUAD). However, the prognostic value of DNA methylation remains to be further explored. Objectives. We sought to explore DNA methylation characteristics and develop a quantifiable criterion related to DNA methylation to improve survival prediction for LUAD patients.

Methods: Illumina Human Methylation450K array data, level 3 RNA-seq data and corresponding clinical information were obtained from TCGA. Cox regression analysis and the Akaike information criterion were used to construct the best-prognosis methylation signature. Receiver operating characteristic curve analysis was used to validate the prognostic ability of the DNA methylation-related feature score. qPCR was used to measure the transcription levels of the identified genes upon methylation.

Results: We identified a set of DNA methylation features composed of 11 genes (MYEOV, KCNU1, SLC27A6, NEUROD4, HMGB4, TACR3, GABRA5, TRPM8, NLRP13, EDN3 and SLC34A1). The feature score, calculated based on DNA methylation features, was independent of tumor recurrence and TNM stage in predicting overall survival. Of note, the combination of this feature score and TNM stage provided a better overall survival prediction than either of them individually. The transcription levels of all the hypermethylated genes were significantly increased after demethylation, and the expression levels of 3 hypomethylated proteins were significantly higher in tumor tissues than in normal tissues, as indicated by immunohistochemistry data from the Human Protein Atlas. Our results suggested that these identified genes with prognostic features were regulated by DNA methylation of their promoters.

Conclusion: Our studies demonstrated the potential application of DNA methylation markers in the prognosis of LUAD.

The Hematopoietic TALE-Code Shows Normal Activity of IRX1 in Myeloid Progenitors and Reveals Ectopic Expression of IRX3 and IRX5 in Acute Myeloid Leukemia

Homeobox genes encode transcription factors that control basic developmental decisions. Knowledge of their hematopoietic activities casts light on normal and malignant immune cell development. Recently, we constructed the so-called lymphoid TALE-code that codifies expression patterns of all active TALE class homeobox genes in early hematopoiesis and lymphopoiesis. Here, we present the corresponding myeloid TALE-code to extend this gene signature, covering the entire hematopoietic system. The collective data showed expression patterns for eleven TALE homeobox genes and highlighted the exclusive expression of IRX1 in megakaryocyte-erythroid progenitors (MEPs), implicating this TALE class member in a specific myeloid differentiation process. Analysis of public profiling data from acute myeloid leukemia (AML) patients revealed aberrant activity of IRX1 in addition to IRX3 and IRX5, indicating an oncogenic role for these TALE homeobox genes when deregulated. Screening of RNA-seq data from 100 leukemia/lymphoma cell lines showed overexpression of IRX1, IRX3, and IRX5 in megakaryoblastic and myelomonocytic AML cell lines, chosen as suitable models for studying the regulation and function of these homeo-oncogenes. Genomic copy number analysis of IRX-positive cell lines demonstrated chromosomal amplification of the neighboring IRX3 and IRX5 genes at position 16q12 in MEGAL, underlying their overexpression in this cell line model. Comparative gene expression analysis of these cell lines revealed candidate upstream factors and target genes, namely the co-expression of GATA1 and GATA2 together with IRX1, and of BMP2 and HOXA10 with IRX3/IRX5. Subsequent knockdown and stimulation experiments in AML cell lines confirmed their activating impact in the corresponding IRX gene expression. Furthermore, we demonstrated that IRX1 activated KLF1 and TAL1, while IRX3 inhibited GATA1, GATA2, and FST. Accordingly, we propose that these regulatory relationships may represent major physiological and oncogenic activities of IRX factors in normal and malignant myeloid differentiation, respectively. Finally, the established myeloid TALE-code is a useful tool for evaluating TALE homeobox gene activities in AML.

Epigenetically silenced apoptosis-associated tyrosine kinase (AATK) facilitates a decreased expression of Cyclin D1 and WEE1, phosphorylates TP53 and reduces cell proliferation in a kinase-dependent manner

Silencing of the Apoptosis associated Tyrosine Kinase gene (AATK) has been described in cancer. In our study, we specifically investigated the epigenetic inactivation of AATK in pancreatic adenocarcinoma, lower grade glioma, lung, breast, head, and neck cancer. The resulting loss of AATK correlates with impaired patient survival. Inhibition of DNA methyltransferases (DNMTs) reactivated AATK in glioblastoma and pancreatic cancer. In contrast, epigenetic targeting via the CRISPR/dCas9 system with either EZH2 or DNMT3A inhibited the expression of AATK. Via large-scale kinomic profiling and kinase assays, we demonstrate that AATK acts a Ser/Thr kinase that phosphorylates TP53 at Ser366. Furthermore, whole transcriptome analyses and mass spectrometry associate AATK expression with the GO term 'regulation of cell proliferation'. The kinase activity of AATK in comparison to the kinase-dead mutant mediates a decreased expression of the key cell cycle regulators Cyclin D1 and WEE1. Moreover, growth suppression through AATK relies on its kinase activity. In conclusion, the Ser/Thr kinase AATK represses growth and phosphorylates TP53. Furthermore, expression of AATK was correlated with a better patient survival for different cancer entities. This data suggests that AATK acts as an epigenetically inactivated tumor suppressor gene.

In-silico selection of cancer blood plasma proteins by integrating genomic and proteomic databases

Blood protein markers have been studied for the clinical management of cancer. Due to the large number of the proteins existing in blood, it is often necessary to pre-select potential protein markers before experimental studies. However, to date there is a lack of automated method for in-silico selection of cancer blood proteins that integrates the information from both genetic and proteomic studies in a cancer-specific manner. In this work, we synthesized both genomic and proteomic information from several open access databases and established a bioinformatic pipeline for in-silico selection of blood plasma proteins overexpressed in specific type of cancer. We demonstrated the workflow of this pipeline with an example of breast cancer, while the methodology was applicable for other cancer types. With this pipeline we obtained 10 candidate biomarkers for breast cancer. The proposed pipeline provides a useful and convenient tool for in-silico selection of candidate blood protein biomarkers for a variety of cancer research.

Down-regulation of lncRNA MEG3 promotes chronic low dose cadmium exposure-induced cell transformation and cancer stem cell-like property

Cadmium (Cd) is a toxic heavy metal and one of carcinogens that cause lung cancer. However, the exact mechanism of Cd carcinogenesis remains unclear. To investigate the mechanism of Cd carcinogenesis, we exposed human bronchial epithelial cells (BEAS-2B) to a low dose of Cd (2.5 μM, CdCl₂) for 9 months, which caused cell malignant transformation and generated cancer stem cell (CSC)-like cells. The goal of this study is to investigate the underlying mechanism. The long non-coding RNA (lncRNA) microarray analysis showed that the expression level of a tumor suppressive lncRNA maternally expressed 3 (MEG3) is significantly down-regulated in Cd-transformed cells, which is confirmed by further q-PCR analysis. Mechanistically, it was found that chronic Cd exposure up-regulates the levels of DNA methyltransferases (DNMTs), which increases the methylation of the differentially methylated region (DMR) 1.5 kb upstream of MEG3 transcription start site to reduce MEG3 expression. Functional studies showed that stably overexpressing MEG3 in Cd-transformed cells significantly reduces their transformed phenotypes. Moreover, stably overexpressing MEG3 in parental non-transformed human bronchial epithelial cells significantly impaired the capability of chronic Cd exposure to induce cell transformation and CSC-like property. Further mechanistic studies revealed that the cell cycle inhibitor p21 level is reduced and retinoblastoma protein (Rb) phosphorylation is increased in Cd-transformed cells to promote cell cycle progression. In addition, Cd-transformed cells also expressed higher levels of Bcl-xL and displayed apoptosis resistance. In contrast, stably overexpressing MEG3 increased p21 levels and reduced Rb phosphorylation and Bcl-xL levels in Cd-exposed cells and reduced their cell cycle progression and apoptosis resistance. Together, these findings suggest that MEG3 down-regulation may play important roles in Cd-induced cell transformation and CSC-like property by promoting cell cycle progression and apoptosis resistance.

Exploring Differentially Methylated Genes in Vulvar Squamous Cell Carcinoma

Simple Summary

Vulvar squamous cell carcinoma (VSCC) is the most common form of vulvar malignancy, and its incidence has increased in recent years. For better diagnosis and prognostication, and to expand available treatment options, molecular characterization of VSCC is crucial. We sought to identify aberrations in DNA methylation in VSCC, as this has been implicated in the development of several cancers. To this end, we performed genome-wide methylation sequencing on a set of VSCC and normal vulvar tissue using the Infinium MethylationEPIC BeadChip array. We detected 199 genes to be differentially methylated in VSCC compared to normal vulvar tissue. Of these, 194 genes were hyper-methylated, which leads to a loss of function of the genes. As most of these genes are involved in transcription regulator activity, our results suggest that disruption of this process plays an important role in VSCC development.

Abstract

DNA methylation is the most widely studied mechanism of epigenetic modification, which can influence gene expression without alterations in DNA sequences. Aberrations in DNA methylation are known to play a role in carcinogenesis, and methylation profiling has enabled the identification of biomarkers of potential clinical interest for several cancers. For vulvar squamous cell carcinoma (VSCC), however, methylation profiling remains an under-studied area. We sought to identify differentially methylated genes (DMGs) in VSCC, by performing Infinium MethylationEPIC BeadChip (Illumina) array sequencing, on a set of primary VSCC (n = 18), and normal vulvar tissue from women with no history of vulvar (pre)malignancies (n = 6). Using a false-discovery rate of 0.05, beta-difference (Δβ) of ±0.5, and CpG-island probes as cut-offs, 199 DMGs (195 hyper-methylated, 4 hypo-methylated) were identified for VSCC. Most of the hyper-methylated genes were found to be involved in transcription regulator activity, indicating that disruption of this process plays a vital role in VSCC development. The majority of VSCCs harbored amplifications of chromosomes 3, 8, and 9. We identified a set of DMGs in this exploratory, hypothesis-generating study, which we hope will facilitate epigenetic profiling of VSCCs. Prognostic relevance of these DMGs deserves further exploration in larger cohorts of VSCC and its precursor lesions.