Comparison of different normalization assumptions for analyses of DNA methylation data from the cancer genome

Nowadays, some researchers normalized DNA methylation arrays data in order to remove the technical artifacts introduced by experimental differences in sample preparation, array processing and other factors. However, other researchers analyzed DNA methylation arrays without performing data normalization considering that current normalizations for methylation data may distort real differences between normal and cancer samples because cancer genomes may be extensively subject to hypomethylation and the total amount of CpG methylation might differ substantially among samples. In this study, using eight datasets by Infinium HumanMethylation27 assay, we systemically analyzed the global distribution of DNA methylation changes in cancer compared to normal control and its effect on data normalization for selecting differentially methylated (DM) genes. We showed more differentially methylated (DM) genes could be found in the Quantile/Lowess-normalized data than in the non-normalized data. We found the DM genes additionally selected in the Quantile/Lowess-normalized data showed significantly consistent methylation states in another independent dataset for the same cancer, indicating these extra DM genes were effective biological signals related to the disease. These results suggested normalization can increase the power of detecting DM genes in the context of diagnostic markers which were usually characterized by relatively large effect sizes. Besides, we evaluated the reproducibility of DM discoveries for a particular cancer type, and we found most of the DM genes additionally detected in one dataset showed the same methylation directions in the other dataset for the same cancer type, indicating that these DM genes were effective biological signals in the other dataset. Furthermore, we showed that some DM genes detected from different studies for a particular cancer type were significantly reproducible at the functional level.     

High-throughput method for detecting DNA methylation

Aberrant DNA methylation of CpG site is among the earliest and most frequent alterations in cancer. Detection of promoter hypermethylation of cancer-related gene may be useful for cancer diagnosis or the detection of recurrence. However, most of the studies have focused on a single gene only and gave little information about the concurrent methylation status of multiple genes. In this study, we attempted to develop a microarray method coupled with linker-PCR for detecting methylation status of multiple genes in the tumor tissue. A series of synthesized oligonucleotides were synthesised and purified to completely match with 16 investigated targets. Then they were immobilized on the aldehyde-coated glass slide to fabricate a DNA microarray for detecting methylation status of these genes. The results indicated that these genes were all methylated in the positive control. However, no methylated was found in these genes for the negative control. Only p16 and p15 genes were methylated in investigated genes for the gastric tumor tissue, whereas others were not methylated. The above results were validated by bisulfite DNA sequencing. Our experiments successfully demonstrated that the DNA microarray could be applied as a high-throughput tool to determine methylation status of the investigated genes.     

A DNA methylation microarray-based study identifies ERG as a gene commonly methylated in prostate cancer

DNA methylation of promoter regions is a common event in prostate cancer, one of the most common cancers in men worldwide. Because prior reports demonstrating that DNA methylation is important in prostate cancer studied a limited number of genes, we systematically quantified the DNA methylation status of 1505 CpG dinucleotides for 807 genes in 78 paraffin-embedded prostate cancer samples and three normal prostate samples. The ERG gene, commonly repressed in prostate cells in the absence of an oncogenic fusion to the TMPRSS2 gene, was one of the most commonly methylated genes, occurring in 74% of prostate cancer specimens. In an independent group of patient samples, we confirmed that ERG DNA methylation was common, occurring in 57% of specimens, and cancer-specific. The ERG promoter is marked by repressive chromatin marks mediated by polycomb proteins in both normal prostate cells and prostate cancer cells, which may explain ERG's predisposition to DNA methylation and the fact that tumors with ERG DNA methylation were more methylated, in general. These results demonstrate that bead arrays offer a high-throughput method to discover novel genes with promoter DNA methylation such as ERG, whose measurement may improve our ability to more accurately detect prostate cancer.     

Functional impact bias reveals cancer drivers

Identifying cancer driver genes and pathways among all somatic mutations detected in a cohort of tumors is a key challenge in cancer genomics. Traditionally, this is done by prioritizing genes according to the recurrence of alterations that they bear. However, this approach has some known limitations, such as the difficulty to correctly estimate the background mutation rate, and the fact that it cannot identify lowly recurrently mutated driver genes. Here we present a novel approach, Oncodrive-fm, to detect candidate cancer drivers which does not rely on recurrence. First, we hypothesized that any bias toward the accumulation of variants with high functional impact observed in a gene or group of genes may be an indication of positive selection and can thus be used to detect candidate driver genes or gene modules. Next, we developed a method to measure this bias (FM bias) and applied it to three datasets of tumor somatic variants. As a proof of concept of our hypothesis we show that most of the highly recurrent and well-known cancer genes exhibit a clear FM bias. Moreover, this novel approach avoids some known limitations of recurrence-based approaches, and can successfully identify lowly recurrent candidate cancer drivers.     

A DNA methylation microarray-based study identifies ERG as a gene commonly methylated in prostate cancer

DNA methylation of promoter regions is a common event in prostate cancer, one of the most common cancers in men worldwide. Because prior reports demonstrating that DNA methylation is important in prostate cancer studied a limited number of genes, we systematically quantified the DNA methylation status of 1505 CpG dinucleotides for 807 genes in 78 paraffin-embedded prostate cancer samples and three normal prostate samples. The ERG gene, commonly repressed in prostate cells in the absence of an oncogenic fusion to the TMPRSS2 gene, was one of the most commonly methylated genes, occurring in 74% of prostate cancer specimens. In an independent group of patient samples, we confirmed that ERG DNA methylation was common, occurring in 57% of specimens, and cancer-specific. The ERG promoter is marked by repressive chromatin marks mediated by polycomb proteins in both normal prostate cells and prostate cancer cells, which may explain ERG’s predisposition to DNA methylation and the fact that tumors with ERG DNA methylation were more methylated, in general. These results demonstrate that bead arrays offer a high-throughput method to discover novel genes with promoter DNA methylation such as ERG, whose measurement may improve our ability to more accurately detect prostate cancer.     

Epigenetic analysis of childhood acute lymphoblastic leukemia

We used a chromosome 3 wide NotI microarray for identification of epigenetically inactivated genes in childhood acute lymphoblastic leukemia (ALL). Three novel genes demonstrated frequent methylation in childhood ALL. PPP2R3A (protein phosphatase 2, regulatory subunit B'', alpha) was frequently methylated in T (69%) and B (82%)-ALL. Whilst FBLN2 (fibulin 2) and THRB (thyroid hormone receptor, beta) showed frequent methylation in B-ALL (58%; 56% respectively), but were less frequently methylated in T-ALL (17% for both genes). Recently it was demonstrated that BNC1 (Basonuclin 1) and MXS1 (msh homeobox 1) were frequently methylated across common epithelial cancers. In our series of childhood ALL BNC1 was frequently methylated in both T (77%) and B-ALL (79%), whilst MSX1 showed T-ALL (25%) specific methylation. The methylation of the above 5 genes was cancer specific and expression of the genes could be restored in methylated leukemia cell lines treated with 5-aza-2’-deoxycytidine. This is the first report demonstrating frequent epigenetic inactivation of PPP2R3A, FBLN2, THRB, BNC1 and MSX1 in leukemia. The identification of frequently methylated genes showing cancer specific methylation will be useful in developing early cancer detection screens and for targeted epigenetic therapies.     

Whole genome DNA methylation profiling of oral cancer in ethnic population of Meghalaya,North East India reveals novel genes

Oral Squamous Cell Carcinoma (OSCC) is a serious and one of the most common and highly aggressive malignancies. Epigenetic factors such as DNA methylation have been known to be implicated in a number of cancer etiologies. The main objective of this study was to investigate physiognomies of Promoter DNA methylation patterns associated with oral cancer epigenome with special reference to the ethnic population of Meghalaya, North East India. The present study identifies 27,205 CpG sites and 3811 regions that are differentially methylated in oral cancer when compared to matched normal. 45 genes were found to be differentially methylated within the promoter region, of which 38 were hypermethylated and 7 hypomethylated. 14 of the hypermethylated genes were found to be similar to that of the TCGA-HNSCC study some of which are TSGs and few novel genes which may serve as candidate methylation biomarkers for OSCC in this poorly characterized ethnic group.     

Oligonucleotide-based microarray for DNA methylation analysis: principles and applications

Gene silencing via promoter CpG island hypermethylation offers tumor cells growth advantages. This epigenetic event is pharmacologically reversible, and uncovering a unique set of methylation-silenced genes in tumor cells can bring a new avenue to cancer treatment. However, high-throughput tools capable of surveying the methylation status of multiple gene promoters are needed for this discovery process. Herein we describe an oligonucleotide-based microarray technique that is both versatile and sensitive in revealing hypermethylation in defined regions of the genome. DNA samples are bisulfite-treated and PCR-amplified to distinguish CpG dinucleotides that are methylated from those that are not. Fluorescently labeled PCR products are hybridized to arrayed oligonucleotides that can discriminate between methylated and unmethylated alleles in regions of interest. Using this technique, two clinical subtypes of non-Hodgkin's lymphomas, mantle cell lymphoma, and grades I/II follicular lymphoma, were further separated based on the differential methylation profiles of several gene promoters. Work is underway in our laboratory to extend the interrogation power of this microarray system in multiple candidate genes. This novel tool, therefore, holds promise to monitor the outcome of various epigenetic therapies on cancer patients.     

Frequent epigenetic inactivation of RASSF10 in thyroid cancer

The Ras association domain family (RASSF) encodes for distinct tumor suppressors and several members are frequently silenced in human cancer. In our study, we analyzed the role of a novel RASSF member termed RASSF10 in thyroid carcinogenesis. The RASSF10 CpG island promoter was intensively methylated in nine thyroid cancer cell lines and in 66% of primary thyroid carcinomas. RASSF10 methylation was significantly increased in primary thyroid carcinoma compared to normal thyroid and follicular adenoma (0% and 10%, respectively; p     

Formamide as a denaturant for bisulfite conversion of genomic DNA: Bisulfite sequencing of the GSTPi and RARβ2 genes of 43 formalin-fixed paraffin-embedded prostate cancer specimens

Analysis of methylated DNA, which refers to 5-methycytosine (5mC) versus cytosine (C) at specific loci in genomic DNA (gDNA), has received increased attention in epigenomics, particularly in the area of cancer biomarkers. Many different methods for analysis of methylated DNA rely on initial reaction of gDNA with concentrated acidic sodium bisulfite to quantitatively convert C to uracil (U) via sulfonation of denatured, single-stranded gDNA under conditions where 5mC is resistant to analogous sulfonation leading to thymine (T). These methods typically employ polymerase chain reaction (PCR) amplification after bisulfite conversion, thereby leading to readily detectable amounts of amplicons where T and C are measured as surrogates for C and 5mC in the original unconverted gDNA. However, incomplete bisulfite conversion of C in gDNA has been reported to be a common source of error in analysis of methylated DNA. Incomplete conversion can be revealed during the course of bisulfite sequencing, which is the generally accepted “gold standard” for analysis of methylated DNA. Previous bisulfite sequencing investigations of conventional predenaturation of gDNA with NaOH followed by the use of bisulfite containing added urea to maintain denaturation and thus mitigate incomplete conversion of C have been reported to give conflicting results. The current study describes a new approach where conventional predenaturation of gDNA with NaOH is instead achieved with formamide and maintains denaturation during subsequent sample handling and sulfonation. This formamide-based method was applied to 46 formalin-fixed/paraffin-embedded (FFPE) biopsy tissue specimens from well-characterized patients with primary prostate cancer. These specimens were representative of difficult-to-analyze samples due to the chemically compromised nature of the gDNA, which was recovered by modifying the protocol for a commercially available total RNA/DNA extraction kit (RecoverALL). An additional novel aspect of this study was analysis of CpG-rich promoter regions of two prostate cancer-related genes: glutathione S-transferase pi (GSTPi) and retinoic acid receptor beta2 (RARβ2). High-quality bisulfite sequencing results were obtained for both genes in 43 of 46 (93%) specimens. Detection of methylated GSTPi and RARβ2 genes was significantly associated with primary prostate cancer as compared with the benign prostate (Fisher’s exact test, P < 0.001). The sensitivity and specificity of detection of methylated GSTPi and RARβ2 genes were 86% and 100% and 91% and 100%, respectively. Moreover, the presence of either methylated gene was detected in primary prostate cancer with sensitivity and specificity of 100% and 100%, respectively. The results demonstrated a high degree of reliability of formamide-based denaturation and bisulfite conversion that should extend, generally, to FFPE and other types of samples intended for any analytical method predicated on bisulfite conversion. This pilot study also demonstrated the efficacy of determining methylation of these two genes with high sensitivity and specificity in FFPE biopsy tissue specimens. Moreover, the results showed a highly significant association of methylated GSTPi and RARβ2 genes with primary prostate cancer. Finally, this improved procedure for determining these two methylated genes may allow the detection of prostate cancer cells in core biopsy specimens with insufficient numbers of cells and poor morphology.     

Identification of New Candidate Genes and Chemicals Related to Esophageal Cancer Using a Hybrid Interaction Network of Chemicals and Proteins

Cancer is a serious disease responsible for many deaths every year in both developed and developing countries. One reason is that the mechanisms underlying most types of cancer are still mysterious, creating a great block for the design of effective treatments. In this study, we attempted to clarify the mechanism underlying esophageal cancer by searching for novel genes and chemicals. To this end, we constructed a hybrid network containing both proteins and chemicals, and generalized an existing computational method previously used to identify disease genes to identify new candidate genes and chemicals simultaneously. Based on jackknife test, our generalized method outperforms or at least performs at the same level as those obtained by a widely used method - the Random Walk with Restart (RWR). The analysis results of the final obtained genes and chemicals demonstrated that they highly shared gene ontology (GO) terms and KEGG pathways with direct and indirect associations with esophageal cancer. In addition, we also discussed the likelihood of selected candidate genes and chemicals being novel genes and chemicals related to esophageal cancer.     

Constitutive activation of peroxisome proliferator-activated receptor-gamma suppresses pro-inflammatory adhesion molecules in human vascular endothelial cells

Peroxisome proliferator-activated receptor-gamma (PPAR-gamma) is a ligand-activated nuclear receptor that has an essential role in adipogenesis and glucose homeostasis. PPAR-gamma is expressed in vascular tissues including endothelial cells (ECs). PPAR-gamma activity can be regulated by many pathophysiological and pharmacological agonists. However, the role of PPAR-gamma activation in ECs remains unclear. In this study, we examined the effect of the constitutive activation of PPAR-gamma on the phenotypic modulation of ECs. Adenovirus-mediated expression of a constitutively active mutant of PPAR-gamma resulted in significant ligand-independent activation of PPAR-gamma and specific induction of the PPAR-gamma target genes. However, PPAR-gamma activation significantly suppressed the expression of vascular adhesion molecules in ECs and the ensuing leukocyte recruitment. Furthermore, constitutive activation of PPAR-gamma resulted in simultaneous repression of AP-1 and NF-kappaB activity, which suggests that PPAR-gamma may reduce pro-inflammatory phenotypes via, at least in part, suppression of the AP-1 and NF-kappaB pathways. Therefore, using a gain-of-function approach, our study provides novel evidence showing that constitutive activation of PPAR-gamma is sufficient to prevent ECs from converting into a pro-inflammatory phenotype. These results also suggest that, in addition to pharmacological agonists, the genetic modification of the PPAR-gamma activity in ECs may be a potential approach for therapeutic intervention in various inflammatory disorders.     

Genome-wide DNA methylation profiling of stomach cancer in the ethnic population of Mizoram,North East India

Stomach cancer is the fifth most common cancer in terms of prevalence and incidence and the fourth leading cause of mortality in men and women worldwide. It is well-established that aberrant DNA methylation in cells can lead to carcinogenesis. The primary objective of our study was to investigate the aberrant DNA methylation status of genes associated with stomach cancer with a particular reference to the ethnic population of Mizoram, North East India. The site-level analysis identified 2883 CpG sites differentially methylated, representing ～922 genes. Out of which 476 Differentially Methylated Positions (DMPs) were promoter-associated, 452 DMPs were hypermethylated, and 24 were hypomethylated. The region-level analysis identified 462 Differentially Methylated Regions (DMRs) corresponding to ～320 genes, of which ～281 genes were hypermethylated and ～40 genes were hypomethylated. TCGA analysis showed that some of the genes had been previously implicated in other cancers including stomach cancer. Five hypermethylated genes were selected as candidate genes for further investigations and they have shown to be novel and could serve as candidate hypermethylation biomarkers for stomach cancer in this particular ethnic group.     

A novel method to assess the full genome methylation profile using monoclonal antibody combined with the high throughput based microarray approach

In this work we demonstrate a novel method of methylation detection that utilises immunoaffinity to detect the presence of methylated DNA hybridised to a cDNA microarray. We use a monoclonal antibody specific to 5 methyl cytidine to detect the presence of 5 methyl cytosine in genomic DNA from human fibroblasts bearing the karyotype 45 XO. We report that over 2900 genes show the presence of methylation in this condition. We also report that 165 genes are consistently methylated in all replicates of these experiments. The methylated genes show a uniform distribution over all the chromosomes. The gene ontology of these also indicates no functional correlation between the genes that are methylated. We detect the presence of methylation in IGF2, an imprinted gene and thus known to harbour DNA methylation. The method is extremely specific and offers a quick and efficient way to analyse the methylation landscape on a high throughput scale. This method uses existing technology to assess methylation and thus can integrate very efficiently into any platform used.     

Methylomic analysis identifies frequent DNA methylation of zinc finger protein 582 (ZNF582) in cervical neoplasms

Background

Despite of the trend that the application of DNA methylation as a biomarker for cancer detection is promising, clinically applicable genes are few. Therefore, we looked for novel hypermethylated genes for cervical cancer screening.

Methods and Findings

At the discovery phase, we analyzed the methylation profiles of human cervical carcinomas and normal cervixes by methylated DNA immunoprecipitation coupled to promoter tiling arrays (MeDIP-on-chip). Methylation-specific PCR (MSP), quantitative MSP and bisulfite sequencing were used to verify the methylation status in cancer tissues and cervical scrapings from patients with different severities. Immunohistochemical staining of a cervical tissue microarray was used to confirm protein expression. We narrowed to three candidate genes: DBC1, PDE8B, and ZNF582; their methylation frequencies in tumors were 93%, 29%, and 100%, respectively. At the pre-validation phase, the methylation frequency of DBC1 and ZNF582 in cervical scraping correlated significantly with disease severity in an independent cohort (n = 330, both P<0.001). For the detection of cervical intraepithelial neoplasia 3 (CIN3) and worse, the area under the receiver operating characteristic curve (AUC) of ZNF582 was 0.82 (95% confidence interval  = 0.76–0.87).

Conclusions

Our study shows ZNF582 is frequently methylated in CIN3 and worse lesions, and it is demonstrated as a potential biomarker for the molecular screening of cervical cancer.     

Expression of HYOU1 via Reciprocal Crosstalk between NSCLC Cells and HUVECs Control Cancer Progression and Chemoresistance in Tumor Spheroids

Among all cancer types, lung cancer ranks highest worldwide in terms of both incidence and mortality. The crosstalk between lung cancer cells and their tumor microenvironment (TME) has begun to emerge as the “Achilles heel” of the disease and thus constitutes an attractive target for anticancer therapy. We previously revealed that crosstalk between lung cancer cells and endothelial cells (ECs) induces chemoresistance in multicellular tumor spheroids (MCTSs). In this study, we demonstrated that factors secreted in response to crosstalk between ECs and lung cancer cells play pivotal roles in the development of chemoresistance in lung cancer spheroids. We subsequently determined that the expression of hypoxia up-regulated protein 1 (HYOU1) in lung cancer spheroids was increased by factors secreted in response to crosstalk between ECs and lung cancer cells. Direct interaction between lung cancer cells and ECs also caused an elevation in the expression of HYOU1 in MCTSs. Inhibition of HYOU1 expression not only suppressed stemness and malignancy, but also facilitated apoptosis and chemosensitivity in lung cancer MCTSs. Inhibition of HYOU1 expression also significantly increased the expression of interferon signaling components in lung cancer cells. Moreover, the activation of the PI3K/AKT/mTOR pathway was involved in the HYOU1-induced aggression of lung cancer cells. Taken together, our results identify HYOU1, which is induced in response to crosstalk between ECs and lung cancer cells within the TME, as a potential therapeutic target for combating the aggressive behavior of cancer cells.     

Attenuated NER Expressions of XPF and XPC Associated with Smoking Are Involved in the Recurrence of Bladder Cancer

The varied NER genes and smoking are two important risk factors of bladder cancer, but the mechanism of the NER protein and smoking in cancer progression, however, remains unclear. In this report, we compared the expressions of NER genes in 79 bladder cancer tissues with or without any recurrence by real-time PCR and then analyzed the varied NER genes by immunochemistry in 219 bladder cancer tissue samples. Based on the clinical data, we analyzed the clinical value of varied NER genes and smoking in 219 bladder cancers by the Kaplan-Meier method and Cox proportional hazards regression. We found the expressions of the NER gene XPF and XPC were significantly lower in bladder cancer tissues with a recurrence compared with those without a recurrence at mRNA level. Also, the patients with the XPF and XPC defect had a statistically significant lower median recurrence-free survival time than those without the XPF and XPC defect, and smoking can make this difference more remarkable. Our results suggest that XPF and XPC expression may be a potential predictive factor for bladder cancer, and smoking can not only influence the recurrence of bladder cancer as a single factor but also aggravate the results of the XPF defect and XPC defect.