Genome-wide screen of promoter methylation identifies novel markers in melanoma

Development and validation of a prognostic model for colon cancer based on mitotic gene signatures and immune microenvironment analysis

BACKGROUND

Mitotic processes play a pivotal role in tumor progression and immune responses. However, the correlation between mitosis-related genes, clinical outcomes, and the tumor microenvironment (TME) in colon cancer remains unclear. This study aims to develop a prognostic and therapeutic significance model for colon cancer based on mitosis-related genes.

METHODS

RNA expression profiles and clinical data of 453 colon cancer patients were downloaded from The Cancer Genome Atlas (TCGA). Mitosis-related genes were selected from the MsigDb database. The gene model was constructed using differential analysis, univariate and multivariate Cox regression, and Lasso regression analyses. The predictive model was validated using data from the GSE17536, GSE17537, and GSE39582 datasets. Predictive accuracy was evaluated via Receiver Operating Characteristic (ROC) curves, while nomograms were developed by integrating clinical and pathological features. Gene set enrichment analysis explored biological processes and pathways linked to the model. TME was assessed using ESTIMATE, and the proportion and function of immune cells were analyzed through CIBERSORT. Drug sensitivity analysis was conducted using the CTRP database.

RESULTS

A predictive model based on 17 mitosis-related genes (KIFC1, CCNF, EME1, CDC25C, ORC1, CCNJL, ANKRD53, MEIS2, FZD3, TPD52L1, MAPK3, CDKN2A, EDN3, NPM2, PSRC1, INHBA, BIRC5) was created. The model exhibited robust predictive performance across both training and validation cohorts. Nomograms for predicting 3-, 5-, and 7-year survival rates in colon cancer (COAD) patients were generated. The model's correlation with immune cell infiltration and function was highlighted.

CONCLUSION

The mitosis-related gene model serves as a valuable indicator for predicting survival outcomes in colon cancer patients.

Uncovering Minimal Pathways in Melanoma Initiation

Cutaneous melanomas are clinically and histologically heterogeneous. Most display activating mutations in Braf or Nras and complete loss of function of one or more tumor suppressor genes. Mouse models that replicate such mutations produce fast-growing, pigmented tumors. However, mice that combine Braf activation with only heterozygous loss of Pten also produce tumors and, as we show here, in an Albino background this occurs even with Braf activation alone. Such tumors arise rarely, grow slowly, and express low levels of pigmentation genes. The timing of their appearance was consistent with a single step stochastic event, but no evidence could be found that it required de novo mutation, suggesting instead the involvement of an epigenetic transition. Single-cell transcriptomic analysis revealed such tumors to be heterogeneous, including a minor cell type we term LNM ( L ow-pigment, N eural- and extracellular M atrix-signature) that displays gene expression resembling "neural crest"-like cell subsets detected in the fast-growing tumors of more heavily-mutated mice, as well as in human biopsy and xenograft samples. We provide evidence that LNM cells pre-exist in normal skin, are expanded by Braf activation, can transition into malignant cells, and persist with malignant cells through multiple rounds of transplantation. We discuss the possibility that LNM cells not only serve as a pre-malignant state in the production of some melanomas, but also as an important intermediate in the development of drug resistance.

Cellular functions of heat shock protein 20 (HSPB6) in cancer: A review

Heat shock proteins (HSP) are a large family of peptide proteins that are widely found in cells. Studies have shown that the expression and function of HSPs in cells are very complex, and they can participate in cellular physiological and pathological processes through multiple pathways. Multiple heat shock proteins are associated with cancer cell growth, proliferation, metastasis, and resistance to anticancer drugs, and they play a key role in cancer development by ensuring the correct folding or degradation of proteins in cancer cells. As research hotspots, HSP90, HSP70 and HSP27 have been extensively studied in cancer so far. However, HSP20, also referred to as HSPB6, as a member of the small heat shock protein family, has been shown to play an important role in the cardiovascular system, but little research has been conducted on HSP20 in cancer. This review summarizes the current cellular functions of HSP20 in different cancer types, as well as its effects on cancer proliferation, progression, prognosis, and its other functions in cancer, to illustrate the close association between HSP20 and cancer. We show that, unlike most HSPs, HSP20 mainly plays an active anticancer role in cancer development, which is expected to provide new ideas and help for cancer diagnosis and treatment and research.

Functional analysis of recurrent CDC20 promoter variants in human melanoma

Small nucleotide variants in non-coding regions of the genome can alter transcriptional regulation, leading to changes in gene expression which can activate oncogenic gene regulatory networks. Melanoma is heavily burdened by non-coding variants, representing over 99% of total genetic variation, including the well-characterized TERT promoter mutation. However, the compendium of regulatory non-coding variants is likely still functionally under-characterized. We developed a pipeline to identify hotspots, i.e. recurrently mutated regions, in melanoma containing putatively functional non-coding somatic variants that are located within predicted melanoma-specific regulatory regions. We identified hundreds of statistically significant hotspots, including the hotspot containing the TERT promoter variants, and focused on a hotspot in the promoter of CDC20. We found that variants in the promoter of CDC20, which putatively disrupt an ETS motif, lead to lower transcriptional activity in reporter assays. Using CRISPR/Cas9, we generated an indel in the CDC20 promoter in human A375 melanoma cell lines and observed decreased expression of CDC20, changes in migration capabilities, increased growth of xenografts, and an altered transcriptional state previously associated with a more proliferative and less migratory state. Overall, our analysis prioritized several recurrent functional non-coding variants that, through downregulation of CDC20, led to perturbation of key melanoma phenotypes.

Novel CAF-identifiers via transcriptomic and protein level analysis in HNSC patients

Cancer-associated fibroblasts (CAFs), a prominent component of the tumor microenvironment, play an important role in tumor development, invasion, and drug resistance. The expression of distinct "CAF-markers" which separates CAFs from normal fibroblasts and epithelial cells, have traditionally been used to identify them. These commonly used CAF-markers have been reported to differ greatly across different CAF subpopulations, even within a cancer type. Using an unbiased -omic approach from public data and in-house RNAseq data from patient derived novel CAF cells, TIMP-1, SPARC, COL1A2, COL3A1 and COL1A1 were identified as potential CAF-markers by differential gene expression analysis using publicly available single cell sequencing data and in-house RNAseq data to distinguish CAF populations from tumor epithelia and normal oral fibroblasts. Experimental validation using qPCR and immunofluorescence revealed CAF-specific higher expression of TIMP-1 and COL1A2 as compared to other markers in 5 novel CAF cells, derived from patients of diverse gender, habits and different locations of head and neck squamous cell carcinoma (HNSC). Upon immunohistochemical (IHC) analysis of FFPE blocks however, COL1A2 showed better differential staining between tumor epithelia and tumor stroma. Similar data science driven approach utilizing single cell sequencing and RNAseq data from stabilized CAFs can be employed to identify CAF-markers in various cancers.

HSPB6 Is Depleted in Colon Cancer Patients and Its Expression Is Induced by 5-aza-2'-Deoxycytidine In Vitro

Background and Objectives: Colon cancer (CC) is the second most common cancer in Saudi Arabia, and the number of new cases is expected to increase by 40% by 2040. Sixty percent of patients with CC are diagnosed in the late stage, causing a reduced survival rate. Thus, identifying a new biomarker could contribute to diagnosing CC in the early stages, leading to delivering better therapy and increasing the survival rate. Materials and Methods: HSPB6 expression was investigated in extracted RNA taken from 10 patients with CC and their adjacent normal tissues, as well as in DMH-induced CC and a colon treated with saline taken from a male Wistar rat. Additionally, the DNA of the LoVo and Caco-2 cell lines was collected, and bisulfite was converted to measure the DNA methylation level. This was followed by applying 5-aza-2'-deoxycytidine (AZA) to the LoVo and Caco-2 cell lines for 72 h to see the effect of DNA methylation on HSPB6 expression. Finally, the GeneMANIA database was used to find the interacted genes at transcriptional and translational levels with HSPB6. Results: We found that the expression of HSPB6 was downregulated in 10 CC tissues compared to their adjacent normal colon tissues, as well as in the in vivo study, where its expression was lower in the colon treated with the DMH agent compared to the colon treated with saline. This suggests the possible role of HSPB6 in tumor progression. Moreover, HSPB6 was methylated in two CC cell lines (LoVo and Caco-2), and demethylation with AZA elevated its expression, implying a mechanistic association between DNA methylation and HSPB6 expression. Conclusions: Our findings indicate that HSPB6 is adversely expressed with tumor progression, and its expression may be controlled by DNA methylation. Thus, HSPB6 could be a good biomarker employed in the CC diagnostic process.

Epigenetically regulated PCDHB15 impairs aggressiveness of metastatic melanoma cells

The protocadherin proteins are cell adhesion molecules at the crossroad of signaling pathways playing a major role in neuronal development. It is now understood that their role as signaling hubs is not only important for the normal physiology of cells but also for the regulation of hallmarks of cancerogenesis. Importantly, protocadherins form a cluster of genes that are regulated by DNA methylation. We have identified for the first time that PCDHB15 gene is DNA-hypermethylated on its unique exon in the metastatic melanoma-derived cell lines and patients' metastases compared to primary tumors. This DNA hypermethylation silences the gene, and treatment with the DNA demethylating agent 5-aza-2'-deoxycytidine reinduces its expression. We explored the role of PCDHB15 in melanoma aggressiveness and showed that overexpression impairs invasiveness and aggregation of metastatic melanoma cells in vitro and formation of lung metastasis in vivo. These findings highlight important modifications of the methylation of the PCDHβ genes in melanoma and support a functional role of PCDHB15 silencing in melanoma aggressiveness.

Current understanding of epigenetics role in melanoma treatment and resistance

Melanoma is the most aggressive form of skin cancer resulting from genetic mutations in melanocytes. Several factors have been considered to be involved in melanoma progression, including genetic alteration, processes of damaged DNA repair, and changes in mechanisms of cell growth and proliferation. Epigenetics is the other factor with a crucial role in melanoma development. Epigenetic changes have become novel targets for treating patients suffering from melanoma. These changes can alter the expression of microRNAs and their interaction with target genes, which involves cell growth, differentiation, or even death. Given these circumstances, we conducted the present review to discuss the melanoma risk factors and represent the current knowledge about the factors related to its etiopathogenesis. Moreover, various epigenetic pathways, which are involved in melanoma progression, treatment, and chemo-resistance, as well as employed epigenetic factors as a solution to the problems, will be discussed in detail.

Characterization of the CpG Island Hypermethylated Phenotype Subclass in Primary Melanomas

Cutaneous melanoma can be lethal even if detected at an early stage. Epigenetic profiling may facilitate the identification of aggressive primary melanomas with unfavorable outcomes. We performed clustering of whole-genome methylation data to identify subclasses that were then assessed for survival, clinical features, methylation patterns, and biological pathways. Among 89 cutaneous primary invasive melanomas, we identified three methylation subclasses exhibiting low methylation, intermediate methylation, or hypermethylation of CpG islands, known as the CpG island methylator phenotype (CIMP). CIMP melanomas occurred as early as tumor stage 1b and, compared with low-methylation melanomas, were associated with age at diagnosis ≥65 years, lentigo maligna melanoma histologic subtype, presence of ulceration, higher American Joint Committee on Cancer stage and tumor stage, and lower tumor-infiltrating lymphocyte grade (all P < 0.05). Patients with CIMP melanomas had worse melanoma-specific survival (hazard ratio = 11.84; confidence interval = 4.65‒30.20) than those with low-methylation melanomas, adjusted for age, sex, American Joint Committee on Cancer stage, and tumor-infiltrating lymphocyte grade. Genes hypermethylated in CIMP compared with those in low-methylation melanomas included PTEN, VDR, PD-L1, TET2, and gene sets related to development/differentiation, the extracellular matrix, and immunity. CIMP melanomas exhibited hypermethylation of genes important in melanoma progression and tumor immunity, and although present in some early melanomas, CIMP was associated with worse survival independent of known prognostic factors.

A Novel Multi-Omics Analysis Model for Diagnosis and Survival Prediction of Lower-Grade Glioma Patients

Background

Lower-grade gliomas (LGGs) are characterized by remarkable genetic heterogeneity and different clinical outcomes. Classification of LGGs is improved by the development of molecular stratification markers including IDH mutation and 1p/19q chromosomal integrity, which are used as a hallmark of survival and therapy sensitivity of LGG patients. However, the reproducibility and sensitivity of the current classification remain ambiguous. This study aimed to construct more accurate risk-stratification approaches.

Methods

According to bioinformatics, the sequencing profiles of methylation and transcription and imaging data derived from LGG patients were analyzed and developed predictable risk score and radiomics score. Moreover, the performance of predictable models was further validated.

Results

In this study, we determined a cluster of 6 genes that were correlated with IDH mutation/1p19q co-deletion status. Risk score model was calculated based on 6 genes and showed gratifying sensitivity and specificity for survival prediction and therapy response of LGG patients. Furthermore, a radiomics risk score model was established to noninvasively assist judgment of risk score in pre-surgery. Taken together, a predictable nomogram that combined transcriptional signatures and clinical characteristics was established and validated to be preferable to the histopathological classification. Our novel multi-omics nomograms showed a satisfying performance. To establish a user-friendly application, the nomogram was further developed into a web-based platform: https://drw576223193.shinyapps.io/Nomo/, which could be used as a supporting method in addition to the current histopathological-based classification of gliomas.

Conclusions

Our novel multi-omics nomograms showed the satisfying performance of LGG patients and assisted clinicians to draw up individualized clinical management.

NPM2 in malignant peritoneal mesothelioma: from basic tumor biology to clinical medicine

BACKGROUND

This review systematically summarizes gene biology features and protein structure of nucleoplasmin2 (NPM2) and the relationship between NPM2 and malignant peritoneal mesothelioma (MPM), in order to explore the molecular pathological mechanism of MPM and explore new therapeutic targets.

METHODS

NCBI PubMed database was used for the literature search. NCBI Gene and Protein databases, Ensembl Genome Browser, UniProt, and RCSB PDB database were used for gene and protein review. Three online tools (Consurf, DoGSiteScorer, and ZdockServer), the GEPIA database, and the Cancer Genome Atlas were used to analyze bioinformatics characteristics for NPM2 protein.

RESULTS

The main structural domains of NPM2 protein include the N-terminal core region, acidic region, and motif and disordered region. The N-terminal core region, involved in histone binding, is the most conserved domain in the nucleoplasmin (NPM) family. NPM2 with a large acidic tract in its C-terminal tail (NPM2-A2) is able to bind histones and form large complexes. Bioinformatics results indicated that NPM2 expression was correlated with the pathology of multiple tumors. Among mesothelioma patients, 5-year survival of patients with low-NPM2-expression was significantly higher than that of the high-NPM2-expression patients. NPM2 can facilitate the formation of histone deacetylation. NPM2 may promote histone deacetylation and inhibit the related-gene transcription, thus leading to abnormal proliferation, invasion, and metastasis of MPM.

CONCLUSION

NPM2 may play a key role in the development and progression of MPM.

The Role of Extracellular Matrix Remodeling in Skin Tumor Progression and Therapeutic Resistance

The extracellular matrix remodeling in the skin results from a delicate balance of synthesis and degradation of matrix components, ensuring tissue homeostasis. These processes are altered during tumor invasion and growth, generating a microenvironment that supports growth, invasion, and metastasis. Apart from the cellular component, the tumor microenvironment is rich in extracellular matrix components and bound factors that provide structure and signals to the tumor and stromal cells. The continuous remodeling in the tissue compartment sustains the developing tumor during the various phases providing matrices and proteolytic enzymes. These are produced by cancer cells and stromal fibroblasts. In addition to fostering tumor growth, the expression of specific extracellular matrix proteins and proteinases supports tumor invasion after the initial therapeutic response. Lately, the expression and structural modification of matrices were also associated with therapeutic resistance. This review will focus on the significant alterations in the extracellular matrix components and the function of metalloproteinases that influence skin cancer progression and support the acquisition of therapeutic resistance.

Epigenetics in the Diagnosis and Therapy of Malignant Melanoma

Epigenetic mechanisms are fundamentally important for cancer initiation and development. However, a survey of the literature reveals that, to date, they appear less comprehensively investigated in melanoma than in many other cancers, e.g., prostate, breast, and colon carcinoma. The aim of this review is to provide a short summary of epigenetic aspects of functional relevance for melanoma pathogenesis. In addition, some new perspectives from epigenetic research in other cancers with potential for melanoma diagnosis and therapy are introduced. For example, the PrimeEpiHit hypothesis in urothelial carcinoma, which, similarly to malignant melanoma, can also be triggered by a single exogenous noxa, states that one of the first steps for cancer initiation could be epigenetic changes in key genes of one-carbon metabolism. The application of such insights may contribute to further progress in the diagnosis and therapy of melanoma, a deadly type of cancer.

Signal pathways of melanoma and targeted therapy

Melanoma is the most lethal skin cancer that originates from the malignant transformation of melanocytes. Although melanoma has long been regarded as a cancerous malignancy with few therapeutic options, increased biological understanding and unprecedented innovations in therapies targeting mutated driver genes and immune checkpoints have substantially improved the prognosis of patients. However, the low response rate and inevitable occurrence of resistance to currently available targeted therapies have posed the obstacle in the path of melanoma management to obtain further amelioration. Therefore, it is necessary to understand the mechanisms underlying melanoma pathogenesis more comprehensively, which might lead to more substantial progress in therapeutic approaches and expand clinical options for melanoma therapy. In this review, we firstly make a brief introduction to melanoma epidemiology, clinical subtypes, risk factors, and current therapies. Then, the signal pathways orchestrating melanoma pathogenesis, including genetic mutations, key transcriptional regulators, epigenetic dysregulations, metabolic reprogramming, crucial metastasis-related signals, tumor-promoting inflammatory pathways, and pro-angiogenic factors, have been systemically reviewed and discussed. Subsequently, we outline current progresses in therapies targeting mutated driver genes and immune checkpoints, as well as the mechanisms underlying the treatment resistance. Finally, the prospects and challenges in the development of melanoma therapy, especially immunotherapy and related ongoing clinical trials, are summarized and discussed.

Predicting the possible effect of miR-203a-3p and miR-29a-3p on DNMT3B and GAS7 genes expression

Aberrant expression of genes involved in methylation, including DNA methyltransferase 3 Beta (DNMT3B), can cause hypermethylation of various tumor suppressor genes. In this regard, various molecular factors such as microRNAs can play a critical role in regulating these methyltransferase enzymes and eventually downstream genes such as growth arrest specific 7 (GAS7). Accordingly, in the present study we aimed to predict regulatory effect of miRNAs on DNMT3B and GAS7 genes expression in melanoma cell line. hsa-miR-203a-3p and hsa-miR-29a-3p were predicted and selected using bioinformatics software. The Real-time PCR technique was performed to investigate the regulatory effect of these molecules on the DNMT3B and GAS7 genes expression. Expression analysis of DNMT3B gene in A375 cell line showed that there was a significant increase compared to control (p value = 0.0015). Analysis of hsa-miR-203a-3p and hsa-miR-29a-3p indicated the insignificant decreased expression in melanoma cell line compared to control (p value < 0.05). Compared to control, the expression of GAS7 gene in melanoma cells showed a significant decrease (p value = 0.0323). Finally, our findings showed that the decreased expression of hsa-miR-203a-3p and hsa-miR-29a-3p can hypothesize that their aberrant expression caused DNMT3B dysfunction, possible methylation of the GAS7 gene, and ultimately decreased its expression. However, complementary studies are necessary to definite comment.

Epigenetic Regulation in Melanoma: Facts and Hopes

Cutaneous melanoma is a lethal disease, even when diagnosed in advanced stages. Although recent progress in biology and treatment has dramatically improved survival rates, new therapeutic approaches are still needed. Deregulation of epigenetics, which mainly controls DNA methylation status and chromatin remodeling, is implied not only in cancer initiation and progression, but also in resistance to antitumor drugs. Epigenetics in melanoma has been studied recently in both melanoma preclinical models and patient samples, highlighting its potential role in different phases of melanomagenesis, as well as in resistance to approved drugs such as immune checkpoint inhibitors and MAPK inhibitors. This review summarizes what is currently known about epigenetics in melanoma and dwells on the recognized and potential new targets for testing epigenetic drugs, alone or together with other agents, in advanced melanoma patients.

An integrated epigenomic-transcriptomic landscape of lung cancer reveals novel methylation driver genes of diagnostic and therapeutic relevance

Background: Aberrant DNA methylation occurs commonly during carcinogenesis and is of clinical value in human cancers. However, knowledge of the impact of DNA methylation changes on lung carcinogenesis and progression remains limited.

Methods: Genome-wide DNA methylation profiles were surveyed in 18 pairs of tumors and adjacent normal tissues from non-small cell lung cancer (NSCLC) patients using Reduced Representation Bisulfite Sequencing (RRBS). An integrated epigenomic-transcriptomic landscape of lung cancer was depicted using the multi-omics data integration method.

Results: We discovered a large number of hypermethylation events pre-marked by poised promoter in embryonic stem cells, being a hallmark of lung cancer. These hypermethylation events showed a high conservation across cancer types. Eight novel driver genes with aberrant methylation (e.g., PCDH17 and IRX1) were identified by integrated analysis of DNA methylome and transcriptome data. Methylation level of the eight genes measured by pyrosequencing can distinguish NSCLC patients from lung tissues with high sensitivity and specificity in an independent cohort. Their tumor-suppressive roles were further experimentally validated in lung cancer cells, which depend on promoter hypermethylation. Similarly, 13 methylation-driven ncRNAs (including 8 lncRNAs and 5 miRNAs) were identified, some of which were co-regulated with their host genes by the same promoter hypermethylation. Finally, by analyzing the transcription factor (TF) binding motifs, we uncovered sets of TFs driving the expression of epigenetically regulated genes and highlighted the epigenetic regulation of gene expression of TCF21 through DNA methylation of EGR1 binding motifs.

Conclusions: We discovered several novel methylation driver genes of diagnostic and therapeutic relevance in lung cancer. Our findings revealed that DNA methylation in TF binding motifs regulates target gene expression by affecting the binding ability of TFs. Our study also provides a valuable epigenetic resource for identifying DNA methylation-based diagnostic biomarkers, developing cancer drugs for epigenetic therapy and studying cancer pathogenesis.

Global DNA methylation pattern involved in the modulation of differentiation potential of adipogenic and myogenic precursors in skeletal muscle of pigs

Background

Skeletal muscle is a complex and heterogeneous tissue accounting for approximately 40% of body weight. Excessive ectopic lipid accumulation in the muscle fascicle would undermine the integrity of skeletal muscle in humans but endow muscle with marbling-related characteristics in farm animals. Therefore, the balance of myogenesis and adipogenesis is of great significance for skeletal muscle homeostasis. Significant DNA methylation occurs during myogenesis and adipogenesis; however, DNA methylation pattern of myogenic and adipogenic precursors derived from skeletal muscle remains unknown yet.

Methods

In this study, reduced representation bisulfite sequencing was performed to analyze genome-wide DNA methylation of adipogenic and myogenic precursors derived from the skeletal muscle of neonatal pigs. Integrated analysis of DNA methylation and transcription profiles was further conducted. Based on the results of pathway enrichment analysis, myogenic precursors were transfected with CACNA2D2-overexpression plasmids to explore the function of CACNA2D2 in myogenic differentiation.

Results

As a result, 11,361 differentially methylated regions mainly located in intergenic region and introns were identified. Furthermore, 153 genes with different DNA methylation and gene expression level between adipogenic and myogenic precursors were characterized. Subsequently, pathway enrichment analysis revealed that DNA methylation programing was involved in the regulation of adipogenic and myogenic differentiation potential through mediating the crosstalk among pathways including focal adhesion, regulation of actin cytoskeleton, MAPK signaling pathway, and calcium signaling pathway. In particular, we characterized a new role of CACNA2D2 in inhibiting myogenic differentiation by suppressing JNK/MAPK signaling pathway.

Conclusions

This study depicted a comprehensive landmark of DNA methylome of skeletal muscle-derived myogenic and adipogenic precursors, highlighted the critical role of CACNA2D2 in regulating myogenic differentiation, and illustrated the possible regulatory ways of DNA methylation on cell fate commitment and skeletal muscle homeostasis.

Supplementary information

The online version contains supplementary material available at 10.1186/s13287-020-02053-3.

Combined Analysis of DNA Methylome and Transcriptome Reveal Novel Candidate Genes Related to Porcine Escherichia coli F4ab/ac-Induced Diarrhea

Enterotoxigenic Escherichia coli (ETEC) that express F4 (K88) fimbriae are the principal microorganisms responsible for bacterial diarrhea in neonatal and pre-weaning piglets. To better understand the molecular effects of ETEC F4ab/ac infection, we performed a genome-wide comparison of the changes in DNA methylation and gene expression in ETEC F4ab/ac infected porcine intestinal epithelial cells. We characterized the pattern of changes in methylation and found 3297 and 1593 differentially methylated regions in cells infected with F4ab and F4ac, respectively. Moreover, 606 and 780 differentially expressed genes (DEGs) in ETEC F4ab and F4ac infected cells were detected and these genes were highly enriched in immune/defense response related pathways. Integrative analysis identified 27 and 10 genes showing inverse correlations between promoter methylation and expression with ETEC F4ab/ac infection. Altered DNA methylation and expression of various genes suggested their roles and potential functional interactions upon ETEC F4ab/ac infection. Further functional analyses revealed that three DEGs (S100A9, SGO1, and ESPL1) in F4ab infected cells and three DEGs (MAP3K21, PAK6, and MPZL1) in F4ac infected cells are likely involved in the host cells response to ETEC infection. Our data provides further insight into the epigenetic and transcriptomic alterations of ETEC F4ab/ac infected porcine intestinal epithelial cells, and may advance the identification of biomarkers and drug targets for predicting susceptibility to and controlling ETEC F4ab/ac induced diarrhea.

Higher methylation subtype of malignant melanoma and its correlation with thicker progression and worse prognosis

Malignant melanoma (MM) is the most life‐threatening disease among all skin malignancies, and recent genome‐wide studies reported BRAF, RAS, and NF1 as the most frequently mutated driver genes. While epigenetic aberrations are known to contribute to the oncogenic activity seen in various cancers, their role in MM has not been fully investigated. To investigate the role of epigenetic aberrations in MM, we performed genome‐wide DNA methylation analysis of 51 clinical MM samples using Infinium 450k beadarray. Hierarchical clustering analysis stratified MM into two DNA methylation epigenotypes: high‐ and low‐methylation subgroups. Tumor thickness was significantly greater in case of high‐methylation tumors than low‐methylation tumors (8.3 ± 5.3 mm vs 4.5 ± 2.9 mm, P = .003). Moreover, prognosis was significantly worse in high‐methylation cases (P = .03). Twenty‐seven genes were found to undergo significant and frequent hypermethylation in high‐methylation subgroup, where TFPI2 was identified as the most frequently hypermethylated gene. MM cases with lower expression levels of TFPI2 showed significantly worse prognosis (P = .001). Knockdown of TFPI2 in two MM cell lines, CHL‐1 and G361, resulted in significant increases of cell proliferation and invasion. These indicate that MM can be stratified into at least two different epigenetic subgroups, that the MM subgroup with higher DNA methylation shows a more progressive phenotype, and that methylation of TFPI2 may contribute to the tumor progression of MM.

Clinical significance of COL1A1 and COL1A2 expression levels in hypopharyngeal squamous cell carcinoma

Alterations in collagen type I α1 (COL1A1) and collagen type I α 2 (COL1A2) expression levels have been reported to predict prognosis in various types of cancer. However, the effect of these biomarkers on hypopharyngeal squamous cell carcinoma (HPSCC) is yet to be fully elucidated. The present study aimed to explore the prognostic significance of COL1A1 and COL1A2 expression levels in HPSCC. The expression levels of COL1A1 and COL1A2 in 67 patients with HPSCC were examined using an immunohistochemical assay in a tissue microarray. The associations between COL1A1/COL1A2 expression levels and patient clinicopathological features were analyzed using ANOVA, Pearson's χ² or Fisher's exact test. The Cox proportional hazard models and Kaplan-Meier survival analysis with log-rank tests were used to analyze the significance of COL1A1/COL1A2 as prognostic markers for patients with HPSCC. As a result, immunohistochemical staining revealed that COL1A1 was positively expressed in all cases, among which 40.3% were strong positive, while COL1A2 was positively expressed in 76.1% of the HPSCC cases with 6.0% of the samples exhibiting strong staining. Further analysis revealed no significant association between the expression levels of COL1A1/COL1A2 and other clinicopathological features. Cox regression analysis revealed that a high COL1A2 expression level predicted a high locoregional recurrence and a less favorable disease-free survival rate (P=0.042 and 0.020, respectively). Overall, the present study indicated that COL1A2 expression levels may have value as a prognostic indicator in HPSCC.

DNA hypermethylation is associated with invasive phenotype of malignant melanoma

Tumor cell invasion is one of the key processes during cancer progression, leading to life-threatening metastatic lesions in melanoma. As methylation of cancer-related genes plays a fundamental role during tumorigenesis and may lead to cellular plasticity which promotes invasion, our aim was to identify novel epigenetic markers on selected invasive melanoma cells. Using Illumina BeadChip assays and Affymetrix Human Gene 1.0 microarrays, we explored the DNA methylation landscape of selected invasive melanoma cells and examined the impact of DNA methylation on gene expression patterns. Our data revealed predominantly hypermethylated genes in the invasive cells affecting the neural crest differentiation pathway and regulation of the actin cytoskeleton. Integrative analysis of the methylation and gene expression profiles resulted in a cohort of hypermethylated genes (IL12RB2, LYPD6B, CHL1, SLC9A3, BAALC, FAM213A, SORCS1, GPR158, FBN1 and ADORA2B) with decreased expression. On the other hand, hypermethylation in the gene body of the EGFR and RBP4 genes was positively correlated with overexpression of the genes. We identified several methylation changes that can have role during melanoma progression, including hypermethylation of the promoter regions of the ARHGAP22 and NAV2 genes that are commonly altered in locally invasive primary melanomas as well as during metastasis. Interestingly, the down-regulation of the methylcytosine dioxygenase TET2 gene, which regulates DNA methylation, was associated with hypermethylated promoter region of the gene. This can probably lead to the observed global hypermethylation pattern of invasive cells and might be one of the key changes during the development of malignant melanoma cells.

Evaluation of MT Family Isoforms as Potential Biomarker for Predicting Progression and Prognosis in Gastric Cancer

Background

Metallothioneins (MTs) family comprises many isoforms, most of which are frequently dysregulated in a wide range of cancers. However, the expression pattern and exact role of each distinct MT family isoform which contributes to tumorigenesis, progression, and drug resistance of gastric cancer (GC) are still unclear.

Methods

Publicly available databases including Oncomine, Gene Expression Profiling Interactive Analysis (GEPIA), Kaplan-Meier plotter, SurvExpress, MethHC, cBioportal, and GeneMANIA were accessed to perform an integrated bioinformatic analysis and try to detect fundamental relationships between each MT family member and GC.

Results

Bioinformatic data indicated that the mRNA expression of all MT family members was almost lowly expressed in GC compared with normal gastric tissue (P<0.05), and patients with reduced mRNA expression of each individual MT member had inconsistent prognostic value (OS, FP, PPS), which depended on the individual isoform of MT. A negative correlation between the methylation in promoter region of majority of MT members and their mRNA expression was detected from MethHC database (p<0.001). Data downloaded from TCGA revealed that MTs were rarely mutated in GC patients and MT2A was frequently regulated by other three genes (FOS, JUN, SP1) in GC patients.

Conclusion

MTs were nearly downregulated, and distinct type of MT harbored different prognostic role in GC patients. Methylation in gene promoter region of MTs partially contributed to their reduced expression in GC. Our comprehensive analyses from multiple independent databases may further lead researches to explore MT-targeting reagents or potential diagnostic and prognostic markers for GC patients.

Clinical protein science in translational medicine targeting malignant melanoma

Melanoma of the skin is the sixth most common type of cancer in Europe and accounts for 3.4% of all diagnosed cancers. More alarming is the degree of recurrence that occurs with approximately 20% of patients lethally relapsing following treatment. Malignant melanoma is a highly aggressive skin cancer and metastases rapidly extend to the regional lymph nodes (stage 3) and to distal organs (stage 4). Targeted oncotherapy is one of the standard treatment for progressive stage 4 melanoma, and BRAF inhibitors (e.g. vemurafenib, dabrafenib) combined with MEK inhibitor (e.g. trametinib) can effectively counter BRAFV600E-mutated melanomas. Compared to conventional chemotherapy, targeted BRAFV600E inhibition achieves a significantly higher response rate. After a period of cancer control, however, most responsive patients develop resistance to the therapy and lethal progression. The many underlying factors potentially causing resistance to BRAF inhibitors have been extensively studied. Nevertheless, the remaining unsolved clinical questions necessitate alternative research approaches to address the molecular mechanisms underlying metastatic and treatment-resistant melanoma. In broader terms, proteomics can address clinical questions far beyond the reach of genomics, by measuring, i.e. the relative abundance of protein products, post-translational modifications (PTMs), protein localisation, turnover, protein interactions and protein function. More specifically, proteomic analysis of body fluids and tissues in a given medical and clinical setting can aid in the identification of cancer biomarkers and novel therapeutic targets. Achieving this goal requires the development of a robust and reproducible clinical proteomic platform that encompasses automated biobanking of patient samples, tissue sectioning and histological examination, efficient protein extraction, enzymatic digestion, mass spectrometry-based quantitative protein analysis by label-free or labelling technologies and/or enrichment of peptides with specific PTMs. By combining data from, e.g. phosphoproteomics and acetylomics, the protein expression profiles of different melanoma stages can provide a solid framework for understanding the biology and progression of the disease. When complemented by proteogenomics, customised protein sequence databases generated from patient-specific genomic and transcriptomic data aid in interpreting clinical proteomic biomarker data to provide a deeper and more comprehensive molecular characterisation of cellular functions underlying disease progression. In parallel to a streamlined, patient-centric, clinical proteomic pipeline, mass spectrometry-based imaging can aid in interrogating the spatial distribution of drugs and drug metabolites within tissues at single-cell resolution. These developments are an important advancement in studying drug action and efficacy in vivo and will aid in the development of more effective and safer strategies for the treatment of melanoma. A collaborative effort of gargantuan proportions between academia and healthcare professionals has led to the initiation, establishment and development of a cutting-edge cancer research centre with a specialisation in melanoma and lung cancer. The primary research focus of the European Cancer Moonshot Lund Center is to understand the impact that drugs have on cancer at an individualised and personalised level. Simultaneously, the centre increases awareness of the relentless battle against cancer and attracts global interest in the exceptional research performed at the centre.

Effect of Sulforaphane and 5-Aza-2'-Deoxycytidine on Melanoma Cell Growth

Background: UV exposure-induced oxidative stress is implicated as a driving mechanism for melanoma. Increased oxidative stress results in DNA damage and epigenetic dysregulation. Accordingly, we explored whether a low dose of the antioxidant sulforaphane (SFN) in combination with the epigenetic drug 5-aza-2’-deoxycytidine (DAC) could slow melanoma cell growth. SFN is a natural bioactivated product of the cruciferous family, while DAC is a DNA methyltransferase inhibitor. Methods: Melanoma cell growth characteristics, gene transcription profiles, and histone epigenetic modifications were measured after single and combination treatments with SFN and DAC. Results: We detected melanoma cell growth inhibition and specific changes in gene expression profiles upon combinational treatments with SFN and DAC, while no significant alterations in histone epigenetic modifications were observed. Dysregulated gene transcription of a key immunoregulator cytokine—C-C motif ligand 5 (CCL-5)—was validated. Conclusions: These results indicate a potential combinatorial effect of a dietary antioxidant and an FDA-approved epigenetic drug in controlling melanoma cell growth.

Key candidate genes of STAT1 and CXCL10 in melanoma identified by integrated bioinformatical analysis

The underlying mechanisms and gene signatures of melanoma are unknown. In this study, three expression profile data sets (GSE65568, GSE100050, GSE114445) were integrated to identify candidate genes explaining the pathways and functions of melanoma. Expression data sets including 24 melanoma tumours and 13 normal skin samples were merged and analysed in detail. The three GSE profiles shared 431 differentially expressed genes (DEGs), including 227 upregulated genes, 200 downregulated genes and 4 differentially regulated genes. Moreover, the functions and signalling pathways of the shared DEGs with significant p-values were identified. The two most significant modules were filtered from the DEGs protein-protein interaction (PPI) network, which consisted of 284 nodes. We also plotted the prognostic value of hub genes from an online database. In summary, using integrated bioinformatic analysis, we have identified candidate DEGs and pathways in melanoma that could improve our understanding of the causes and underlying molecular events of melanoma, and these candidate genes and pathways could be therapeutic targets for melanoma.

Recently Evolved Tumor Suppressor Transcript TP73-AS1 Functions as Sponge of Human-Specific miR-941

MicroRNA (miRNA) sponges are vital components of posttranscriptional gene regulation. Yet, only a limited number of miRNA sponges have been identified. Here, we show that the recently evolved noncoding tumor suppressor transcript, antisense RNA to TP73 gene (TP73-AS1), functions as a natural sponge of human-specific miRNA miR-941. We find unusually nine high-affinity miR-941 binding sites clustering within 1 kb region on TP73-AS1, which forms miR-941 sponge region. This sponge region displays increased sequence constraint only in humans, and its formation can be traced to the tandem expansion of a 71-nt-long sequence containing a single miR-941 binding site in old world monkeys. We further confirm TP73-AS1 functions as an efficient miR-941 sponge based on massive transcriptome data analyses, wound-healing assay, and Argonaute protein immunoprecipitation experiments conducted in cell lines. The expression of miR-941 and its sponge correlate inversely across multiple healthy and cancerous tissues, with miR-941 being highly expressed in tumors and preferentially repressing tumor suppressors. Thus, the TP73-AS1 and miR-941 duo represents an unusual case of the extremely rapid evolution of noncoding regulators controlling cell migration, proliferation, and tumorigenesis.

Gene Expression and Methylation Analysis in Melanomas and Melanocytes From the Same Patient: Loss of NPM2 Expression Is a Potential Immunohistochemical Marker for Melanoma

DNA methylation is considered the primary epigenetic mechanism underlying the development of malignant melanoma. Since DNA methylation can be influenced by environmental factors, it is preferable to compare cancer and normal cells from the same patient. In order to compare the methylation status in melanoma tissues and melanocytes from the same individuals, we employed a novel epidermal sheet cultivation technique to isolate normal melanocytes from unaffected sites of melanoma patients. We also analyzed primary and metastatic melanoma samples, three commercially available melanocytes, and four melanoma cell lines. Cluster analysis of DNA methylation data classified freshly isolated melanomas and melanocytes into the same group, whereas the four melanoma cell lines were clustered together in a distant clade. Moreover, our analysis discovered methylation at several novel loci (KRTCAP3, AGAP2, ZNF490), in addition to those identified in previous studies (COL1A2, GPX3); however, the latter two were not observed in fresh melanoma samples. Subsequent studies revealed that NPM2 was hypermethylated and downregulated in melanomas, which was consistent with previous reports. In many normal melanocytes, NPM2 showed distinct immunohistochemical staining, while its expression was lost in malignant melanoma cells. In particular, intraepithelial lesions of malignant melanoma, an important challenge in clinical practice, could be distinguished from benign nevi. The present findings indicate the importance of using fresh melanoma samples, not melanoma cell lines and melanocytes in epigenetic studies. In addition, NPM2 immunoreactivity could be used to differentiate melanomas from normal melanocytes or benign disease.

Demethylation by low-dose 5-aza-2'-deoxycytidine impairs 3D melanoma invasion partially through miR-199a-3p expression revealing the role of this miR in melanoma

BACKGROUND

Efficient treatments against metastatic melanoma dissemination are still lacking. Here, we report that low-cytotoxic concentrations of 5-aza-2'-deoxycytidine, a DNA demethylating agent, prevent in vitro 3D invasiveness of metastatic melanoma cells and reduce lung metastasis formation in vivo.

RESULTS

We unravelled that this beneficial effect is in part due to MIR-199A2 re-expression by promoter demethylation. Alone, this miR showed an anti-invasive and anti-metastatic effect. Throughout integration of micro-RNA target prediction databases with transcriptomic analysis after 5-aza-2'-deoxycytidine treatments, we found that miR-199a-3p downregulates set of genes significantly involved in invasion/migration processes. In addition, analysis of data from melanoma patients showed a stage- and tissue type-dependent modulation of MIR-199A2 expression by DNA methylation.

CONCLUSIONS

Thus, our data suggest that epigenetic- and/or miR-based therapeutic strategies can be relevant to limit metastatic dissemination of melanoma.

Mass spectrometric imaging of cysteine rich proteins in human skin

Looking insight pathological processes, metallothioneins (MTs) are considered to be potential biomarkers for monitoring of a development of various types of diseases, such as cancer. The early identification of the MTs in biological tissues could be important tool for the estimation of appropriate clinical therapy. Therefore, here we investigated the application of matrix assisted laser desorption/ionization mass spectrometry imaging (MALDI MSI) together with immunohistochemical analyses (IHC) using MT-1/2 antibody for MT detection in formalin-fixed paraffin-embedded (FFPE) biopsy specimens of human skin. Principal component analyses revealed differences in the peptide/protein profiles separating healthy skin from the carcinoma specimens. Statistically significant ion peaks at m/z 6038, 6300, 6676, and 7026 were more frequently detected in squamous cell carcinoma (SCC), basal cell carcinoma (BCC) and melanoma. Using IHC, we found that MT-1/2 was significantly higher in SCC and melanoma compared to healthy skin. Surprisingly, significantly low levels of MT-1/2 were found in BCC. On one side, the results indicate important role of MTs in melanoma occurrence and progression, as on the second side, there are hidden processes associated with MTs based on differences of the occurrence of the MS peaks, which could be associated with cycling of MTs isoforms.

Promoter methylation as biomarkers for diagnosis of melanoma: A systematic review and meta-analysis

Melanoma is one of the most common skin cancer that is characterized by rapid growth, early metastasis, high malignant, and mortality. Accumulating evidence demonstrated that promoter methylation of tumor-suppressor genes is implicated in the pathogenesis of melanoma. In the current study, we performed a meta-analysis to identify promising methylation biomarkers in the diagnosis of melanoma. We carried out a systematic literature search using Pubmed, Embase, and ISI web knowledge database and found that gene promoter methylation of 50 genes was reported to be associated with the risk of melanoma. Meta-analysis revealed that hypermethylation of claudin 11 (CLDN11; odds ratio [OR], 16.82; 95% confidence interval [CI], 1.97-143.29; p = 0.010), O-6-methylguanine-DNA methyltransferase (MGMT; OR, 5.59; 95% CI, 2.51-12.47; p < 0.0001), cyclin-dependent kinase inhibitor 2A (p16; OR, 6.57; 95% CI, 2.19-19.75; p = 0.0008), retinoic acid receptor β (RAR-β2; OR, 24.31; 95% CI, 4.58-129.01; p = 0.0002), and Ras association domain family member (RASSF1A; OR, 9.35; 95% CI, 4.73-18.45; p < 0.00001) was significantly higher in melanoma patients compared with controls. CLDN11 (OR, 14.52; 95% CI, 1.84-114.55; p = 0.01), MGMT (OR, 8.08; 95% CI, 1.84-35.46; p = 0.006), p16 (OR, 9.44; 95% CI, 2.68-33.29; p = 0.0005), and RASSF1A (OR, 7.72; 95% CI, 1.05-56.50; p = 0.04) hypermethylation was significantly increased in primary melanoma compared with controls. Methylation frequency of CLDN11 (OR, 25.56; 95% CI, 2.32-281.66; p = 0.008), MGMT (OR, 4.64; 95% CI, 1.98-10.90; p = 0.0004), p16 (OR, 4.31; 95% CI, 1.33-13.96; p = 0.01), and RASSF1A (OR, 10.10; 95% CI, 2.87-35.54; p = 0.0003) was significantly higher in metastasis melanoma compared with controls. These findings indicated that CLDN11, MGMT, p16, RAR-β2, and RASSF1A hypermethylation is a risk factor and a potential biomarker for melanoma. CLDN11, MGMT, p16, and RASSF1A promoter methylation may take part in the development of melanoma and become useful biomarkers in the early diagnosis of the disease.

Identification of DNA methylation regulated novel host genes relevant to inhibition of virus replication in porcine PK15 cell using double stranded RNA mimics and DNA methyltransferase inhibitor

During RNA viruses's replication, double-stranded RNA (dsRNA) is normally produced and induce host innate immune response. Most of gene activation due cytokine mediated but which are due to methylation mediated is still unknown. In the study, DNA methylome was integrated with our previous transcriptome data to investigate the differentially methylated regions and genes using MeDIP-chip technology. We found that the transcriptional expressions of 15, 37 and 18 genes were negatively related with their promoter DNA methylation levels in the cells treated by PolyI:C, Aza-CdR, as well as PolyI:C plus Aza-CdR, respectively, compared with the untreated cells. GO analysis revealed hypo-methylated genes (BNIP3L and CDK9) and a hyper-methylated gene (ZC3HAV1) involved in the host response to viral replication. Our results suggest that these novel genes targeted by DNA methylation can be potential markers relevant to virus replication and host innate immune response to set up a medical model of infectious diseases.

Integrated analysis of DNA methylation profiling and gene expression profiling identifies novel markers in lung cancer in Xuanwei, China

BACKGROUND

Aberrant DNA methylation occurs frequently in cancer. The aim of this study was to identify novel methylation markers in lung cancer in Xuanwei, China, through integrated genome-wide DNA methylation and gene expression studies.

METHODS

Differentially methylated regions (DMRs) and differentially expressed genes (DEGs) were detected on 10 paired lung cancer tissues and noncancerous lung tissues by methylated DNA immunoprecipitation combined with microarray (MeDIP-chip) and gene expression microarray analyses, respectively. Integrated analysis of DMRs and DEGs was performed to screen out candidate methylation-related genes. Both methylation and expression changes of the candidate genes were further validated and analyzed.

RESULTS

Compared with normal lung tissues, lung cancer tissues expressed a total of 6,899 DMRs, including 5,788 hypermethylated regions and 1,111 hypomethylated regions. Integrated analysis of DMRs and DEGs identified 45 tumor-specific candidate genes: 38 genes whose DMRs were hypermethylated and expression was downregulated, and 7 genes whose DMRs were hypomethylated and expression was upregulated. The methylation and expression validation results identified 4 candidate genes (STXBP6, BCL6B, FZD10, and HSPB6) that were significantly hypermethylated and downregulated in most of the tumor tissues compared with the noncancerous lung tissues.

CONCLUSIONS

This integrated analysis of genome-wide DNA methylation and gene expression in lung cancer in Xuanwei revealed several genes regulated by promoter methylation that have not been described in lung cancer before. These results provide new insight into the carcinogenesis of lung cancer in Xuanwei and represent promising new diagnostic and therapeutic targets.

Comprehensive transcriptomic analysis of heat shock proteins in the molecular subtypes of human breast cancer

BACKGROUND

Heat Shock Proteins (HSPs), a family of genes with key roles in proteostasis, have been extensively associated with cancer behaviour. However, the HSP family is quite large and many of its members have not been investigated in breast cancer (BRCA), particularly in relation with the current molecular BRCA classification. In this work, we performed a comprehensive transcriptomic study of the HSP gene family in BRCA patients from both The Cancer Genome Atlas (TCGA) and the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) cohorts discriminating the BRCA intrinsic molecular subtypes.

METHODS

We examined gene expression levels of 1097 BRCA tissue samples retrieved from TCGA and 1981 samples of METABRIC, focusing mainly on the HSP family (95 genes). Data were stratified according to the PAM50 gene expression (Luminal A, Luminal B, HER2, Basal, and Normal-like). Transcriptomic analyses include several statistical approaches: differential gene expression, hierarchical clustering and survival analysis.

RESULTS

Of the 20,531 analysed genes we found that in BRCA almost 30% presented deregulated expression (19% upregulated and 10% downregulated), while of the HSP family 25% appeared deregulated (14% upregulated and 11% downregulated) (|fold change| > 2 comparing BRCA with normal breast tissues). The study revealed the existence of shared HSP genes deregulated in all subtypes of BRCA while other HSPs were deregulated in specific subtypes. Many members of the Chaperonin subfamily were found upregulated while three members (BBS10, BBS12 and CCTB6) were found downregulated. HSPC subfamily had moderate increments of transcripts levels. Various genes of the HSP70 subfamily were upregulated; meanwhile, HSPA12A and HSPA12B appeared strongly downregulated. The strongest downregulation was observed in several HSPB members except for HSPB1. DNAJ members showed heterogeneous expression pattern. We found that 23 HSP genes correlated with overall survival and three HSP-based transcriptional profiles with impact on disease outcome were recognized.

CONCLUSIONS

We identified shared and specific HSP genes deregulated in BRCA subtypes. This study allowed the recognition of HSP genes not previously associated with BRCA and/or any cancer type, and the identification of three clinically relevant clusters based on HSPs expression patterns with influence on overall survival.

A framework for analyzing DNA methylation data from Illumina Infinium HumanMethylation450 BeadChip

Background

DNA methylation has been identified to be widely associated to complex diseases. Among biological platforms to profile DNA methylation in human, the Illumina Infinium HumanMethylation450 BeadChip (450K) has been accepted as one of the most efficient technologies. However, challenges exist in analysis of DNA methylation data generated by this technology due to widespread biases.

Results

Here we proposed a generalized framework for evaluating data analysis methods for Illumina 450K array. This framework considers the following steps towards a successful analysis: importing data, quality control, within-array normalization, correcting type bias, detecting differentially methylated probes or regions and biological interpretation.

Conclusions

We evaluated five methods using three real datasets, and proposed outperform methods for the Illumina 450K array data analysis. Minfi and methylumi are optimal choice when analyzing small dataset. BMIQ and RCP are proper to correcting type bias and the normalized result of them can be used to discover DMPs. R package missMethyl is suitable for GO term enrichment analysis and biological interpretation.

Attenuation of genome-wide 5-methylcytosine level is an epigenetic feature of cutaneous malignant melanomas

Epigenetic modification of DNA, namely covalent changes of cytosine residues, plays a key role in the maintenance of inactive chromatin regions, both in health and in disease. In the vast majority of malignant melanomas, the most notable known epigenetic abnormality is the attenuation of 5-hydroxymethylcytosine (5-hmC) residues. However, it remains unknown whether a decrease in 5-hmC represents a primary defect of melanoma cancer epigenome or whether it is secondary to the loss of 5-methylcytosine (5-mC), a chemical substrate for 5-hmC. Here, we evaluated 5-mC levels in a spectrum of melanocytic proliferations. To study the epigenetic features of melanocytic nuclei, we began by measuring 5-mC levels in histologic specimens semiquantitatively by immunohistochemistry. We next treated established melanoma cell lines with S-adenosyl methionine (SAM), a universal methyl group donor, in an effort to cause changes in 5-mC levels. We detected a marked reduction in 5-mC levels in both primary and metastatic melanomas compared with 5-mC levels in benign melanocytic nevi. We also empirically induced changes in 5-mC in melanoma cell lines by incubation with SAM. To our surprise, we observed a significant cytoreductive effect of SAM on all melanoma cell lines examined. At subcytotoxic levels, SAM treatment is accompanied by a genome-wide increase in 5-mC. Moreover, we recorded a dose-dependent increase in genome-wide 5-mC levels in melanoma cell lines following SAM treatment. Taken together, we report that genome-wide attenuation of 5-mC is a hallmark of malignant melanomas. We propose that genome-wide attenuation of 5-mC is not merely an epiphenomenon as it is required for melanoma cell growth, albeit by an as of yet undetermined mechanism. Given its potential benefit in slowing down the growth of melanoma cells, SAM should be studied further to determine its role in epigenome modulation.

Prognostic value of aberrant promoter hypermethylation of tumor-related genes in early-stage head and neck cancer

Staging and pathological grading are useful, but imperfect predictors of recurrence in head and neck squamous cell carcinoma (HNSCC). Accordingly, molecular biomarkers that predict the risk of recurrence are necessary to improve clinical outcomes. The methylation statuses of the promoters of 11 tumor-related genes (p16, RASSF1A, E-cadherin, H-cadherin, MGMT, DAPK, DCC, COL1A2, TAC1, SST, and GALR1) were analyzed in 133 HNSCC cases using quantitative methylation-specific PCR. We detected frequent methylation of p16 (44%), RASSF1A (18%), E-cadherin (53%), H-cadherin (35%), MGMT (35%), DAPK (53%), DCC (42%), COL1A2 (44%), TAC1 (61%), SST (64%), and GALR1 (44%) in HNSCC. Disease-free survival was lower in patients with 6–11 methylated genes than in those with 0–5 methylated genes (log-rank test, P = 0.001). In a multivariate Cox proportional hazards analysis, the methylation of E-cadherin, COL1A2, TAC1, and GALR1 was associated with poor survival, with hazard ratios of 4.474 (95% CI, 1.241–16.124). In a joint analysis of these four genes, patients with 2–4 methylated genes had a significantly lower survival rate than those with 0–1 methylated genes in early-stage HNSCC. Importantly, the methylation of some genes was closely related to poor prognosis in early-stage HNSCC, providing strong evidence that these hypermethylated genes are valuable biomarkers for prognostic evaluation.

Integrative analysis of the epigenetic basis of muscle-invasive urothelial carcinoma

Background

Elucidation of epigenetic alterations in bladder cancer will lead to further understanding of the biology of the disease and hopefully improved therapies. Our aim was to perform an integrative epigenetic analysis of invasive urothelial carcinoma of the bladder to identify the epigenetic abnormalities involved in the development and progression of this cancer.

Methods

Pre-processed methylation data and RNA-seq data were downloaded from The Cancer Genome Atlas (TCGA) and processed using the R package TCGA-Assembler. An R package MethylMix was used to perform an analysis incorporating both methylation and gene expression data on all samples, as well as a subset analysis comparing patients surviving less than 2 years and patients surviving more than 2 years. Genes associated with poor prognosis were individually queried. Pathway analysis was performed on statistically significant genes identified by MethylMix criteria using ConsensusPathDB. Validation was performed using flow cytometry on bladder cancer cell lines.

Results

A total of 408 patients met all inclusion criteria. There were a total of 240 genes differentially methylated by MethylMix criteria. Review of individual genes specific to poor-prognosis patients revealed the majority to be candidate tumor suppressors in other cancer types. Pathway analysis showed increase in methylation of genes involved in antioxidant pathways including glutathione and NRF2. Genes involved in estrogen metabolism were also hypermethylated while genes involved in the EGFR pathway were found to be hypomethylated. EGFR expression was confirmed to be elevated in six bladder cancer cell lines.

Conclusions

In patients with invasive urothelial carcinoma, we found differential methylation in patients with better and worse prognosis after cystectomy. Differentially methylated genes are involved in many relevant oncologic pathways, including EGFR and antioxidant pathways, that may be a target for therapy or chemoprevention.

Homocysteine inhibits endothelial progenitor cells proliferation via DNMT1-mediated hypomethylation of Cyclin A

Endothelial progenitor cells (EPCs) contribute to neovasculogenesis and reendothelialization of damaged blood vessels to maintain the endothelium. Dysfunction of EPCs is implicated in the pathogenesis of vascular injury induced by homocysteine (Hcy). We aimed to investigate the role of Cyclin A in Hcy-induced EPCs dysfunction and explore its molecular mechanism. In this study, by treatment of EPCs with Hcy, we found that the expression of Cyclin A mRNA and protein were significantly downregulated in a dose-dependent manner. Knockdown of Cyclin A prominently reduced proliferation of EPCs, while over-expression of Cyclin A significantly promoted the cell proliferation, suggesting that Hcy inhibits EPCs proliferation through downregulation of Cyclin A expression. In addition, epigenetic study also demonstrated that Hcy induces DNA hypomethylation of the Cyclin A promoter in EPCs through downregulated expression of DNMT1. Moreover, we found that Hcy treatment of EPCs leads to increased SAM, SAH and MeCP2, while the ratio of SAM/SAH and MBD expression decrease. In summary, our results indicate that Hcy inhibits Cyclin A expression through hypomethylation of Cyclin A and thereby suppress EPCs proliferation. These findings demonstrate a novel mechanism of DNA methylation mediated by DNMT1 in prevention of Hcy associated cardiovascular disease.

An insight into the complex roles of metallothioneins in malignant diseases with emphasis on (sub)isoforms/isoforms and epigenetics phenomena

Metallothioneins (MTs) belong to a group of small cysteine-rich proteins that are ubiquitous throughout all kingdoms. The main function of MTs is scavenging of free radicals and detoxification and homeostating of heavy metals. In humans, 16 genes localized on chromosome 16 have been identified to encode four MT isoforms labelled by numbers (MT-1-MT-4). MT-2, MT-3 and MT-4 proteins are encoded by a single gene. MT-1 comprises many (sub)isoforms. The known active MT-1 genes are MT-1A, -1B, -1E, -1F, -1G, -1H, -1M and -1X. The rest of the MT-1 genes (MT-1C, -1D, -1I, -1J and -1L) are pseudogenes. The expression and localization of individual MT (sub)isoforms and pseudogenes vary at intra-cellular level and in individual tissues. Changes in MT expression are associated with the process of carcinogenesis of various types of human malignancies, or with a more aggressive phenotype and therapeutic resistance. Hence, MT (sub)isoform profiling status could be utilized for diagnostics and therapy of tumour diseases. This review aims on a comprehensive summary of methods for analysis of MTs at (sub)isoforms levels, their expression in single tumour diseases and strategies how this knowledge can be utilized in anticancer therapy.

Bisulfite-independent analysis of CpG island methylation enables genome-scale stratification of single cells

Conventional DNA bisulfite sequencing has been extended to single cell level, but the coverage consistency is insufficient for parallel comparison. Here we report a novel method for genome-wide CpG island (CGI) methylation sequencing for single cells (scCGI-seq), combining methylation-sensitive restriction enzyme digestion and multiple displacement amplification for selective detection of methylated CGIs. We applied this method to analyzing single cells from two types of hematopoietic cells, K562 and GM12878 and small populations of fibroblasts and induced pluripotent stem cells. The method detected 21 798 CGIs (76% of all CGIs) per cell, and the number of CGIs consistently detected from all 16 profiled single cells was 20 864 (72.7%), with 12 961 promoters covered. This coverage represents a substantial improvement over results obtained using single cell reduced representation bisulfite sequencing, with a 66-fold increase in the fraction of consistently profiled CGIs across individual cells. Single cells of the same type were more similar to each other than to other types, but also displayed epigenetic heterogeneity. The method was further validated by comparing the CpG methylation pattern, methylation profile of CGIs/promoters and repeat regions and 41 classes of known regulatory markers to the ENCODE data. Although not every minor methylation differences between cells are detectable, scCGI-seq provides a solid tool for unsupervised stratification of a heterogeneous cell population.

Genome-wide DNA methylation and transcriptome analyses reveal genes involved in immune responses of pig peripheral blood mononuclear cells to poly I:C

DNA methylation changes play essential roles in regulating the activities of genes involved in immune responses. Understanding of variable DNA methylation linked to immune responses may contribute to identifying biologically promising epigenetic markers for pathogenesis of diseases. Here, we generated genome-wide DNA methylation and transcriptomic profiles of six pairs of polyinosinic-polycytidylic acid-treated pig peripheral blood mononuclear cell (PBMC) samples and corresponding controls using methylated DNA immunoprecipitation sequencing and RNA sequencing. Comparative methylome analyses identified 5,827 differentially methylated regions and 615 genes showing differential expression between the two groups. Integrative analyses revealed inverse associations between DNA methylation around transcriptional start site and gene expression levels. Furthermore, 70 differentially methylated and expressed genes were identified such as TNFRSF9, IDO1 and EBI3. Functional annotation revealed the enriched categories including positive regulation of immune system process and regulation of leukocyte activation. These findings demonstrated DNA methylation changes occurring in immune responses of PBMCs to poly I:C stimulation and a subset of genes potentially regulated by DNA methylation in the immune responses. The PBMC DNA methylome provides an epigenetic overview of this physiological system in response to viral infection, and we expect it to constitute a valuable resource for future epigenetic epidemiology studies in pigs.

Current trends in electrochemical sensing and biosensing of DNA methylation

DNA methylation plays an important role in physiological and pathological processes. Several genetic diseases and most malignancies tend to be associated with aberrant DNA methylation. Among other analytical methods, electrochemical approaches have been successfully employed for characterisation of DNA methylation patterns that are essential for the diagnosis and treatment of particular diseases. This article discusses current trends in the electrochemical sensing and biosensing of DNA methylation. Particularly, it provides an overview of applied electrode materials, electrode modifications and biorecognition elements applications with an emphasis on strategies that form the core DNA methylation detection approaches. The three main strategies as (i) bisulfite treatment, (ii) cleavage by restriction endonucleases, and (iii) immuno/affinity reaction were described in greater detail. Additionally, the availability of the reviewed platforms for early cancer diagnosis and the approval of methylation inhibitors for anticancer therapy were discussed.

Comprehensive DNA methylation study identifies novel progression-related and prognostic markers for cutaneous melanoma

BACKGROUND

Cutaneous melanoma is the deadliest skin cancer, with an increasing incidence and mortality rate. Currently, staging of patients with primary melanoma is performed using histological biomarkers such as tumor thickness and ulceration. As disruption of the epigenomic landscape is recognized as a widespread feature inherent in tumor development and progression, we aimed to identify novel biomarkers providing additional clinical information over current factors using unbiased genome-wide DNA methylation analyses.

METHODS

We performed a comprehensive DNA methylation analysis during all progression stages of melanoma using Infinium HumanMethylation450 BeadChips on a discovery cohort of benign nevi (n = 14) and malignant melanoma from both primary (n = 33) and metastatic (n = 28) sites, integrating the DNA methylome with gene expression data. We validated the discovered biomarkers in three independent validation cohorts by pyrosequencing and immunohistochemistry.

RESULTS

We identified and validated biomarkers for, and pathways involved in, melanoma development (e.g., HOXA9 DNA methylation) and tumor progression (e.g., TBC1D16 DNA methylation). In addition, we determined a prognostic signature with potential clinical applicability and validated PON3 DNA methylation and OVOL1 protein expression as biomarkers with prognostic information independent of tumor thickness and ulceration.

CONCLUSIONS

Our data underscores the importance of epigenomic regulation in triggering metastatic dissemination through the inactivation of central cancer-related pathways. Inactivation of cell-adhesion and differentiation unleashes dissemination, and subsequent activation of inflammatory and immune system programs impairs anti-tumoral defense pathways. Moreover, we identify several markers of tumor development and progression previously unrelated to melanoma, and determined a prognostic signature with potential clinical utility.

DNA methylation/hydroxymethylation in melanoma

Melanoma is a malignant tumor of melanocytes and is considered to be the most aggressive cancer among all skin diseases. The pathogenesis of melanoma has not been well documented, which may restrict the research and development of biomarkers and therapies. To date, several genetic and epigenetic factors have been identified as contributing to the development and progression of melanoma. Besides the findings on genetic susceptibilities, the recent progress in epigenetic studies has revealed that loss of the DNA hydroxymethylation mark, 5-hydroxymethylcytosine (5-hmC), along with high levels of DNA methylation at promoter regions of several tumor suppressor genes in melanoma, may serve as biomarkers for melanoma. Moreover, 5-Aza-2′-deoxycytidine, an epigenetic modifier causing DNA demethylation, and ten-eleven translocation family dioxygenase (TET), which catalyzes the generation of 5-hmC, demonstrate therapeutic potential in melanoma treatment. In this review, we will summarize the latest progress in research on DNA methylation/hydroxymethylation in melanoma, and we will discuss and provide insight for epigenetic biomarkers and therapies for melanoma. Particularly, we will discuss the role of DNA hydroxymethylation in melanoma infiltrating immune cells, which may also serve as a potential target for melanoma treatment.

Genome-wide DNA methylation analysis reveals hypomethylation in the low-CpG promoter regions in lymphoblastoid cell lines

BACKGROUND

Epidemiological studies of DNA methylation profiles may uncover the molecular mechanisms through which genetic and environmental factors contribute to the risk of multifactorial diseases. There are two types of commonly used DNA bioresources, peripheral blood cells (PBCs) and EBV-transformed lymphoblastoid cell lines (LCLs), which are available for genetic epidemiological studies. Therefore, to extend our knowledge of the difference in DNA methylation status between LCLs and PBCs is important in human population studies that use these DNA sources to elucidate the epigenetic risks for multifactorial diseases. We analyzed the methylation status of the autosomes for 192 and 92 DNA samples that were obtained from PBCs and LCLs, respectively, using a human methylation 450 K array. After excluding SNP-associated methylation sites and low-call sites, 400,240 sites were subjected to analysis using a generalized linear model with cell type, sex, and age as the independent variables.

RESULTS

We found that the large proportion of sites showed lower methylation levels in LCLs compared with PBCs, which is consistent with previous reports. We also found that significantly different methylation sites tend to be located on the outside of the CpG island and in a region relatively far from the transcription start site. Additionally, we observed that the methylation change of the sites in the low-CpG promoter region was remarkable. Finally, it was shown that the correlation between the chronological age and ageing-associated methylation sites in ELOVL2 and FHL2 in the LCLs was weaker than that in the PBCs.

CONCLUSIONS

The methylation levels of highly methylated sites of the low-CpG-density promoters in PBCs decreased in the LCLs, suggesting that the methylation sites located in low-CpG-density promoters could be sensitive to demethylation in LCLs. Despite being generated from a single cell type, LCLs may not always be a proxy for DNA from PBCs in studies of epigenome-wide analysis attempting to elucidate the role of epigenetic change in disease risks.