DNA methylation subgroups and the CpG island methylator phenotype in gastric cancer: a comprehensive profiling approach

BACKGROUND: Methylation-induced silencing of promoter CpG islands in tumor suppressor genes plays an important role in human carcinogenesis. In colorectal cancer, the CpG island methylator phenotype (CIMP) is defined as widespread and elevated levels of DNA methylation and CIMP+ tumors have distinctive clinicopathological and molecular features. In contrast, the existence of a comparable CIMP subtype in gastric cancer (GC) has not been clearly established. To further investigate this issue, in the present study we performed comprehensive DNA methylation profiling of a well-characterised series of primary GC.

METHODS: The methylation status of 1,421 autosomal CpG sites located within 768 cancer-related genes was investigated using the Illumina GoldenGate Methylation Panel I assay on DNA extracted from 60 gastric tumors and matched tumor-adjacent gastric tissue pairs. Methylation data was analysed using a recursively partitioned mixture model and investigated for associations with clinicopathological and molecular features including age, Helicobacter pylori status, tumor site, patient survival, microsatellite instability and BRAF and KRAS mutations.

RESULTS: A total of 147 genes were differentially methylated between tumor and matched tumor-adjacent gastric tissue, with HOXA5 and hedgehog signalling being the top-ranked gene and signalling pathway, respectively. Unsupervised clustering of methylation data revealed the existence of 6 subgroups under two main clusters, referred to as L (low methylation; 28% of cases) and H (high methylation; 72%). Female patients were over-represented in the H tumor group compared to L group (36% vs 6%; P = 0.024), however no other significant differences in clinicopathological or molecular features were apparent. CpG sites that were hypermethylated in group H were more frequently located in CpG islands and marked for polycomb occupancy.

CONCLUSIONS: High-throughput methylation analysis implicates genes involved in embryonic development and hedgehog signaling in gastric tumorigenesis. GC is comprised of two major methylation subtypes, with the highly methylated group showing some features consistent with a CpG island methylator phenotype.

Comprehensive analysis of CpG island methylator phenotype (CIMP)-high, -low, and -negative colorectal cancers based on protein marker expression and molecular features

CpG island methylator phenotype (CIMP) is being investigated for its role in the molecular and prognostic classification of colorectal cancer patients but is also emerging as a factor with the potential to influence clinical decision-making. We report a comprehensive analysis of clinico-pathological and molecular features (KRAS, BRAF and microsatellite instability, MSI) as well as of selected tumour- and host-related protein markers characterizing CIMP-high (CIMP-H), -low, and -negative colorectal cancers. Immunohistochemical analysis for 48 protein markers and molecular analysis of CIMP (CIMP-H: ? 4/5 methylated genes), MSI (MSI-H: ? 2 instable genes), KRAS, and BRAF were performed on 337 colorectal cancers. Simple and multiple regression analysis and receiver operating characteristic (ROC) curve analysis were performed. CIMP-H was found in 24 cases (7.1%) and linked (p < 0.0001) to more proximal tumour location, BRAF mutation, MSI-H, MGMT methylation (p = 0.022), advanced pT classification (p = 0.03), mucinous histology (p = 0.069), and less frequent KRAS mutation (p = 0.067) compared to CIMP-low or -negative cases. Of the 48 protein markers, decreased levels of RKIP (p = 0.0056), EphB2 (p = 0.0045), CK20 (p = 0.002), and Cdx2 (p < 0.0001) and increased numbers of CD8+ intra-epithelial lymphocytes (p < 0.0001) were related to CIMP-H, independently of MSI status. In addition to the expected clinico-pathological and molecular associations, CIMP-H colorectal cancers are characterized by a loss of protein markers associated with differentiation, and metastasis suppression, and have increased CD8+ T-lymphocytes regardless of MSI status. In particular, Cdx2 loss seems to strongly predict CIMP-H in both microsatellite-stable (MSS) and MSI-H colorectal cancers. Cdx2 is proposed as a surrogate marker for CIMP-H.

The impact of CpG island methylator phenotype and microsatellite instability on tumour budding in colorectal cancer

In colorectal cancer, tumour budding, a process likened to epithelial mesenchymal transition, is an adverse prognostic factor which is rarely found in tumours with high-level microsatellite instability (MSI-H). Cases with MSI-H or high-level CpG island methylator phenotype (CIMP-H) have similar histomorphological features, yet seemingly opposite prognosis. We hypothesized that tumour budding is related to CIMP, thus partially explaining this prognostic difference.

Possible role of Cdx2 in the serrated pathway of colorectal cancer characterized by BRAF mutation, high-level CpG Island methylator phenotype and mismatch repair-deficiency

Colorectal cancer is a heterogeneous disease at the histomorphological, clinical and molecular level. Approximately 20% of cases may progress through the "serrated" pathway characterized by BRAF mutation and high-level CpG Island Methylator Phenotype (CIMP). A large subgroup are additionally microsatellite instable (MSI) and demonstrate significant loss of tumor suppressor Cdx2. The aim of this study is to determine the specificity of Cdx2 protein expression and CpG promoter hypermethylation for BRAF(V600E) and high-level CIMP in colorectal cancer. Cdx2, Mlh1, Msh2, Msh6, and Pms2 were analyzed by immunohistochemistry using a multi-punch tissue microarray (TMA; n = 220 patients). KRAS and BRAF(V600E) mutation analysis, CDX2 methylation and CIMP were investigated. Loss of Cdx2 was correlated with larger tumor size (P = 0.0154), right-sided location (P = 0.0014), higher tumor grade (P < 0.0001), more advanced pT (P = 0.0234) and lymphatic invasion (P = 0.0351). Specificity was 100% for mismatch repair (MMR)-deficiency (P < 0.0001), 92.2% (P < 0.0001) for BRAF(V600E) and 91.8% for CIMP-high. Combined analysis of BRAF(V600E) /CIMP identified Cdx2 loss as sensitive (80%) and specific (91.5%) for mutation/high status. These results were validated on eight well-established colorectal cancer cell lines. CDX2 methylation correlated with BRAF(V600E) (P = 0.0184) and with Cdx2 protein loss (P = 0.0028). These results seem to indicate that Cdx2 may play a role in the serrated pathway to colorectal cancer as underlined by strong relationships with BRAF(V600E) , CIMP-high and MMR-deficiency. Whether this protein can only be used as a "surrogate" marker, or is functionally involved in the progression of these tumors remains to be elucidated.

RNA-SEQUENCING APPLICATIONS: GENE EXPRESSION QUANTIFICATION AND METHYLATOR PHENOTYPE IDENTIFICATION

My dissertation focuses on two aspects of RNA sequencing technology. The first is the methodology for modeling the overdispersion inherent in RNA-seq data for differential expression analysis. This aspect is addressed in three sections. The second aspect is the application of RNA-seq data to identify the CpG island methylator phenotype (CIMP) by integrating datasets of mRNA expression level and DNA methylation status. Section 1: The cost of DNA sequencing has reduced dramatically in the past decade. Consequently, genomic research increasingly depends on sequencing technology. However it remains elusive how the sequencing capacity influences the accuracy of mRNA expression measurement. We observe that accuracy improves along with the increasing sequencing depth. To model the overdispersion, we use the beta-binomial distribution with a new parameter indicating the dependency between overdispersion and sequencing depth. Our modified beta-binomial model performs better than the binomial or the pure beta-binomial model with a lower false discovery rate. Section 2: Although a number of methods have been proposed in order to accurately analyze differential RNA expression on the gene level, modeling on the base pair level is required. Here, we find that the overdispersion rate decreases as the sequencing depth increases on the base pair level. Also, we propose four models and compare them with each other. As expected, our beta binomial model with a dynamic overdispersion rate is shown to be superior. Section 3: We investigate biases in RNA-seq by exploring the measurement of the external control, spike-in RNA. This study is based on two datasets with spike-in controls obtained from a recent study. We observe an undiscovered bias in the measurement of the spike-in transcripts that arises from the influence of the sample transcripts in RNA-seq. Also, we find that this influence is related to the local sequence of the random hexamer that is used in priming. We suggest a model of the inequality between samples and to correct this type of bias. Section 4: The expression of a gene can be turned off when its promoter is highly methylated. Several studies have reported that a clear threshold effect exists in gene silencing that is mediated by DNA methylation. It is reasonable to assume the thresholds are specific for each gene. It is also intriguing to investigate genes that are largely controlled by DNA methylation. These genes are called “L-shaped” genes. We develop a method to determine the DNA methylation threshold and identify a new CIMP of BRCA. In conclusion, we provide a detailed understanding of the relationship between the overdispersion rate and sequencing depth. And we reveal a new bias in RNA-seq and provide a detailed understanding of the relationship between this new bias and the local sequence. Also we develop a powerful method to dichotomize methylation status and consequently we identify a new CIMP of breast cancer with a distinct classification of molecular characteristics and clinical features.

Genome-wide methylation profiling identifies an essential role of reactive oxygen species in pediatric glioblastoma multiforme and validates a methylome specific for H3 histone family 3A with absence of G-CIMP/isocitrate dehydrogenase 1 mutation

Background. Pediatric glioblastoma multiforme (GBM) is rare, and there is a single study, a seminal discovery showing association of histone H3.3 and isocitrate dehydrogenase (IDH) 1 mutation with a DNA methylation signature. The present study aims to validate these findings in an independent cohort of pediatric GBM, compare it with adult GBM, and evaluate the involvement of important functionally altered pathways. Methods. Genome-wide methylation profiling of 21 pediatric GBM cases was done and compared with adult GBM data (GSE22867). We performed gene mutation analysis of IDH1 and H3 histone family 3A (H3F3A), status evaluation of glioma cytosine-phosphate-guanine island methylator phenotype (G-CIMP), and Gene Ontology analysis. Experimental evaluation of reactive oxygen species (ROS) association was also done. Results. Distinct differences were noted between methylomes of pediatric and adult GBM. Pediatric GBM was characterized by 94 hypermethylated and 1206 hypomethylated cytosine-phosphate-guanine (CpG) islands, with 3 distinct clusters, having a trend to prognostic correlation. Interestingly, none of the pediatric GBM cases showed G-CIMP/IDH1 mutation. Gene Ontology analysis identified ROS association in pediatric GBM, which was experimentally validated. H3F3A mutants (36.4%; all K27M) harbored distinct methylomes and showed enrichment of processes related to neuronal development, differentiation, and cell-fate commitment. Conclusions. Our study confirms that pediatric GBM has a distinct methylome compared with that of adults. Presence of distinct clusters and an H3F3A mutation-specific methylome indicate existence of epigenetic subgroups within pediatric GBM. Absence of IDH1/G-CIMP status further indicates that findings in adult GBM cannot be simply extrapolated to pediatric GBM and that there is a strong need for identification of separate prognostic markers. A possible role of ROS in pediatric GBM pathogenesis is demonstrated for the first time and needs further evaluation.

A High Degree of LINE-1 Hypomethylation Is a Unique Feature of Early-Onset Colorectal Cancer

12 p.

The topography of DNA methylation in the non-neoplastic colonic mucosa surrounding colorectal cancers

The degree of gene hypermethylation in non-neoplastic colonic mucosa (NNCM) is a potentially important event in the development of colorectal cancer (CRC), particularly for the subgroup with a CpG island methylator phenotype (CIMP). In this study, we aimed to use an unbiased and high-throughput approach to evaluate the topography of DNA methylation in the non-neoplastic colonic mucosa (NNCM) surrounding colorectal cancer (CRC). A total of 61 tissue samples comprising 53 NNCM and 8 tumor samples were obtained from hemicolectomy specimens of two CRC patients (Cases 1 and 2). NNCM was stripped from the underlying colonic wall and samples taken at varying distances from the tumor. The level of DNA methylation in NNCM and tumor tissues was assessed at 1,505 CpG sites in 807 cancer-related genes using Illumina GoldenGate® methylation arrays. Case 1 tumor showed significantly higher levels of methylation compared to surrounding NNCM samples (P?

Stratification and Prognostic Relevance of Jass's Molecular Classification of Colorectal Cancer

Background: The current proposed model of colorectal tumorigenesis is based primarily on CpG island methylator phenotype (CIMP), microsatellite instability (MSI), KRAS, BRAF, and methylation status of 0-6-Methylguanine DNA Methyltransferase (MGMT) and classifies tumors into five subgroups. The aim of this study is to validate this molecular classification and test its prognostic relevance. Methods: Three hundred two patients were included in this study. Molecular analysis was performed for five CIMP-related promoters (CRABP1, MLH1, p16INK4a, CACNA1G, NEUROG1), MGMT, MSI, KRAS, and BRAF. Methylation in at least 4 promoters or in one to three promoters was considered CIMP-high and CIMP-low (CIMP-H/L), respectively. Results: CIMP-H, CIMP-L, and CIMP-negative were found in 7.1, 43, and 49.9% cases, respectively. One hundred twenty-three tumors (41%) could not be classified into any one of the proposed molecular subgroups, including 107 CIMP-L, 14 CIMP-H, and two CIMP-negative cases. The 10 year survival rate for CIMP-high patients [22.6% (95%CI: 7-43)] was significantly lower than for CIMP-L or CIMP-negative (p = 0.0295). Only the combined analysis of BRAF and CIMP (negative versus L/H) led to distinct prognostic subgroups. Conclusion: Although CIMP status has an effect on outcome, our results underline the need for standardized definitions of low- and high-level CIMP, which clearly hinders an effective prognostic and molecular classification of colorectal cancer.

Loss of Cdx2 Expression in Primary Tumors and Lymph Node Metastases is Specific for Mismatch Repair-Deficiency in Colorectal Cancer

Background: Approximately 20% of all colorectal cancers are hypothesized to arise from the "serrated pathway" characterized by mutation in BRAF, high-level CpG Island Methylator Phenotype, and microsatellite instability/mismatch repair (MMR)-deficiency. MMR-deficient cancers show frequent losses of Cdx2, a homeodomain transcription factor. Here, we determine the predictive value of Cdx2 expression for MMR-deficiency and investigate changes in expression between primary cancers and matched lymph node metastases. Methods: Immunohistochemistry for Cdx2, Mlh1, Msh2, Msh6, and Pms2 was performed on whole tissue sections from 201 patients with primary colorectal cancer and 59 cases of matched lymph node metastases. Receiver operating characteristic curve analysis and Area under the Curve (AUC) were investigated; association of Cdx2 with clinicopathological features and patient survival was carried out. Results: Loss of Cdx2 expression was associated with higher tumor grade (p = 0.0002), advanced pT (p = 0.0166), and perineural invasion (p = 0.0228). Cdx2 loss was an unfavorable prognostic factor in univariate (p = 0.0145) and multivariate [p = 0.0427; HR (95% CI): 0.58 (0.34-0.98)] analysis. The accuracy (AUC) for discriminating MMR-proficient and - deficient cancers was 87% [OR (95% CI): 0.96 (0.95-0.98); p < 0.0001]. Specificity and negative predictive value for MMR-deficiency was 99.1 and 96.3%. One hundred and seventy-four patients had MMR-proficient cancers, of which 60 (34.5%) showed Cdx2 loss. Cdx2 loss in metastases was related to MMR-deficiency (p < 0.0001). There was no difference in expression between primary tumors and matched metastases. Conclusion: Loss of Cdx2 is a sensitive and specific predictor of MMR-deficiency, but is not limited to these tumors, suggesting that events "upstream" of the development of microsatellite instability may impact Cdx2 expression.

Tyrosine kinase receptor B (TrkB) expression in colorectal cancers highlights anoikis resistance as a survival mechanism of tumour budding cells.

AIMS Tumour buds in colorectal cancer represent an aggressive subgroup of non-proliferating and non-apoptotic tumour cells. We hypothesize that the survival of tumour buds is dependent upon anoikis resistance. The role of tyrosine kinase receptor B (TrkB), a promoter of epithelial-mesenchymal transition and anoikis resistance, in facilitating budding was investigated. METHODS AND RESULTS Tyrosine kinase receptor B immunohistochemistry was performed on a multiple-punch tissue microarray of 211 colorectal cancer resections. Membranous/cytoplasmic and nuclear expression was evaluated in tumour and buds. Tumour budding was assessed on corresponding whole tissue slides. Relationship to Ki-67 and caspase-3 was investigated. Analysis of Kirsten Ras (KRAS), proto-oncogene B-RAF (BRAF) and cytosine-phosphate-guanosine island methylator phenotype (CIMP) was performed. Membranous/cytoplasmic and nuclear TrkB were strongly, inversely correlated (P < 0.0001; r = -0.41). Membranous/cytoplasmic TrkB was overexpressed in buds compared to the main tumour body (P < 0.0001), associated with larger tumours (P = 0.0236), high-grade budding (P = 0.0011) and KRAS mutation (P = 0.0008). Nuclear TrkB was absent in buds (P <0.0001) and in high-grade budding cancers (P =0.0073). Among patients with membranous/cytoplasmic TrkB-positive buds, high tumour membranous/cytoplasmic TrkB expression was a significant, independent adverse prognostic factor [P = 0.033; 1.79, 95% confidence interval (CI) 1.05-3.05]. Inverse correlations between membranous/cytoplasmic TrkB and Ki-67 (r = -0.41; P < 0.0001) and caspase-3 (r =-0.19; P < 0.05) were observed. CONCLUSIONS Membranous/cytoplasmic TrkB may promote an epithelial-mesenchymal transition (EMT)-like phenotype with high-grade budding and maintain viability of buds themselves.