- Research article
- Open Access
- Open Peer Review
DNA methylation subgroups and the CpG island methylator phenotype in gastric cancer: a comprehensive profiling approach
BMC Gastroenterologyvolume 14, Article number: 55 (2014)
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.
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.
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.
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.
Gastric cancer (GC) is a complex disease that involves risk factors such as Helicobacter pylori (H. pylori) infection, family history of cancer, environment, diet and genetic susceptibility variants. GC typically has poor prognosis due to late clinical presentation at an advanced stage of disease . Improvements in early detection via screening and the reduction of known risk factors such as chronic H. pylori infection and consumption of preserved/salted food [2–6] has resulted in significantly lower incidence rates in most parts of the world . However, GC remains a major public health issue and is the fourth most common cancer type and the second leading cause of cancer death worldwide [8, 9].
Transcriptional inactivation by cytosine methylation at promoter CpG islands of tumor suppressor genes is an important mechanism contributing to the development of human cancer. In several cancer types, subgroups defined by distinctive methylation patterns have been linked to features such as tumor size in breast cancer , tumor type in lung  and tumor histology in glioma . The most well studied methylation-defined subgroup is the CpG Island Methylator Phenotype (CIMP) in colorectal cancer (CRC) first proposed in 1999 by Toyota et al. . CIMP + CRC exhibit widespread CpG island methylation in gene promoter regions and are characterized by distinct clinical, pathological and molecular features. These include a higher incidence in females and in the proximal colon, poor histological differentiation and frequent association with microsatellite instability (MSI) and BRAF mutations [14, 15]. A panel of five methylation markers has been proposed to standardize the evaluation of CIMP in CRC .
The existence of GC subgroups that are characterized by distinct methylation patterns and/or CIMP-like properties has been explored in several studies [17–26]. However, a standard panel of methylation markers has yet to be proposed for GC and technical issues remain concerning the use of non-quantitative analytical methods and the limited number of genes investigated for methylation. To consolidate knowledge on DNA methylation in GC, we recently performed a meta-analysis of 106 case–control studies that reported on the methylation of 122 candidate genes . A total of 77 genes were found to be differentially methylated between tumor and normal tissue, including genes involved in apoptosis (APAF2, BCL2), cell cycle regulation (p15, p16) and DNA repair (XRCC1). Some studies alluded to the existence of CIMP by referring to a distinct subset of GC that exhibited a high frequency of concurrent gene promoter CpG island hypermethylation. However, the existence and phenotypic properties of CIMP in GC remain controversial, with major confounding factors likely to be the number and identity of CpG sites interrogated for methylation and the GC sample size and quality.
We previously demonstrated that the level of tumor cell content in GC has a major impact on the hierarchical clustering of methylation data . We established that a tumor cell content of 70% was the minimum level required for the reliable analysis of methylation. In the current study, 60 GC samples underwent prior review by pathologists to ensure this minimum tumor cell content was met prior to methylation analysis using the GoldenGate Methylation BeadArray (Illumina) platform. This system allows simultaneous quantification of the methylation level at 1,421 autosomal CpG sites located within 768 cancer-related genes. The aim of our study was therefore to use a comprehensive genome-wide approach to investigate in an unbiased fashion whether methylation subgroups including CIMP occur in GC.
Formalin-fixed and paraffin-embedded (FFPE) primary tumor and matched tumor-adjacent gastric tissue samples from 60 patients with GC were obtained from the Department of Pathology at the National University Hospital System, Singapore, under an institutionally approved protocol. The tumor samples were reviewed for their tumor content and scored in deciles independently by two experienced gastrointestinal pathologists (CLC, BY). All 60 cases included in the study had a tumor cell content of >70% . This cohort has well-annotated clinicopathological information including age, gender, ethnicity, stage, location, tumor size, adenocarcinoma subtype, differentiation, Lauren classification, lymphoid invasion, perineural invasion, H. pylori status, history of chronic gastritis/atrophic gastritis/intestinal metaplasia/dysplasia, overall survival (OS), disease-specific survival (DSS), disease-free survival (DFS) and molecular features such as BRAF V600E, KRAS (codons 12 and 13) mutation and microsatellite instability (MSI).
DNA was extracted from 20 μm sections and verified for DNA quantity and quality as described earlier . The sections were incubated for 3 days at 55°C in 200 μl of digestion buffer (10 mM Tris-hydrochloric acid, pH8.3; 1 mM EDTA; 0.5% Tween 20) and 45 μl of Proteinase K (20 mg/ml, Promega, Madison, WI) without prior dewaxing. The enzyme was inactivated by heating for 10 minutes at 94°C and then samples were centrifuged at 12,000 g for 10 minutes and stored at 4°C without further DNA purification. DNA quantity and quality were determined spectrophotometrically using the NanoDrop ND-1000 (Wilmington, DE). Five hundred nanograms of DNA was bisulfite-converted using the EZ DNA Methylation kit (Zymo Research, Orange, CA) as per the manufacturer’s instructions.
For the validation of candidates, frozen tumour and matched tumor-adjacent tissue from an independent sample series of gastric cancers were obtained from the National University Health System under an institutionally approved protocol. DNA was extracted using the DNeasy Blood and Tissue Kit (Qiagen, Hilden, Germany), quantified and bisulfite converted as described above.
BRAF mutation, KRASmutations and microsatellite instability (MSI)
Hotspot mutations in BRAF (V600E) and KRAS (codons 12 and 13) were detected using direct sequencing as described previously [30, 31]. MSI was determined by analysis of 5 mononucleotide repeats, including BAT-25, BAT-26, NR21, NR22, NR24 and NR27, as reported by Buhard et al. , with tumors being defined as MSI when ≥3 markers showed instability.
Illumina GoldenGate® methylation technology
Comprehensive DNA methylation profiling at 1,505 individual CpG loci contained within 807 genes using the Illumina GoldenGate Methylation Cancer Panel I (Illumina, San Diego, CA) was carried out as described by Bibikova and Fan . Human sperm DNA and Universal methylated DNA (Chemicon, Temcula, CA) were included in each run as unmethylated and methylated controls, respectively. Hybridized arrays were scanned using the BeadArray Reader (Illumina). Normalization of background intensity was estimated from a set of built-in negative controls and subtracted from each methylation data point. To assess sample quality, only those samples having >75% loci with a detection p-value of less than 0.05 were included for analysis. The methylation level at each CpG site or the β-value was defined as the ratio of the methylated allele to the sum of the methylated and unmethylated alleles, and ranged from 0 (completely unmethylated) to 1 (completely methylated).
PCR was performed using 2 μl bisulfite DNA, 1 x reaction buffer with 1 mM MgCl2, 0.8 mM deoxynucleotide triphosphates, 1 unit of FastStart Taq DNA polymerase (Roche Diagnostics, Mannheim, Germany) and 400 mM each of PCR forward primer and a 1:9 mixture of PCR reverse primer and universal biotinylated PCR primer. PCR cycling consisted of incubation at 95°C for 4 min, 50 cycles of 95°C for 30 sec, 54°C for 30 sec and 72°C for 30 sec, followed by a final extension at 72°C for 1 min. Pyrosequencing was performed using the PyroMark annealing buffer (Qiagen) and PyroMark binding buffer (Qiagen), 3 μL Streptavidin Sepharose High Performance beads (GE Healthcare, Stockholm, Sweden) and 350 mM pyrosequencing primer on the PyroMark Q24 (Qiagen) according to manufacturer’s instructions. The PCR forward, reverse and sequencing primers, and pyrosequencing dispensation order were 5′-TTT GGA AGT TAG GAT TTT GG-3′, 5′-GGG ACA CCG CTG ATC GTT TAT CAA TAA AAA AAA AAC AAC CTC AA-3′, 5′- GTT TAT TTA GGG TTG TAA TGT TTT A-3′ and CTA CGA TCT GTC AGT CGT AG respectively for HOXA5, and 5′-GGA GTA AAA TAG GTG AAA GT-3′, 5′-GCC CTT CCC CAA CCT C-3′, 5′-GGT TTT TTT TTT TTA TTA CGT ATT-3′ and GTC AGT TGG TGA respectively for WNT5A. The sequence of the universal biotinylated PCR primer was 5′-GGG ACA CCG CTG ATC GTT TA-3′.
Data from a total of 84 CpG sites contained within 39 X-chromosome genes on the array were removed from the analysis to eliminate gender-specific bias. Thus, 1,421 probes across 768 genes were included for the analyses. All statistical analyses were done in R version 2.14.2 at 5% significance level unless otherwise stated (The R Foundation for Statistical Computing).
The rpmm function in the RPMM library was used for the identification of methylation subgroups. Recursively partitioned mixture model (RPMM) is a model-based unsupervised clustering approach developed for beta-distributed DNA methylation measurements that lie between 0 and 1 . A fanny algorithm was used for initialization and level-weighted version of Bayesian information criterion (BIC) as a split criterion for an existing cluster as implemented in the R-based RPMM package . For the purpose of comparison, classification of tumor samples was also performed with the optimal number of clusters determined using the Calinski-Harabasz pseudo F-statistic , and the robustness evaluated by bootstrap resampling analysis (n = 1000). Graphical representations of the β-values were achieved by the heatmap.plus function with the gplots and heatmap.plus libraries.
Identification of CpG sites that were differentially methylated between tumors and matched tumor-adjacent gastric tissues was performed using the paired sample t-test, while that between methylation subgroups was done using the ANOVA-test. A Benjamini and Hochberg false discovery rate (FDR) cut-off of 0.001 was used, with a Supplementary filter of a minimum difference of 0.15 in the average β-value between the two groups, as described earlier [37, 38]. The associations of methylation subgroups with clinicopathological and molecular factors were compared with the likelihood ratio or Fisher’s Exact test where appropriate. Average methylation level and frequency of methylation (on binarized data) across methylation subgroups was compared using ANOVA and the likelihood ratio test respectively.
All statistical analyses were carried out using the β-value as a continuous variable unless specified otherwise. When β-values were binarized, a methylated threshold of 0.297 was used . Pathway analysis was performed for KEGG pathway mapping using DAVID with a threshold of EASE score ≤ 0.05 [40–42].
DNA methylation patterns in gastric tumor and tumor-adjacenttissues
Unsupervised hierarchical clustering of methylation levels from all 1,421 autosomal CpG sites in 60 tumor samples revealed five distinct subgroups [Additional file 1]. No distinct subgroups were observed for the corresponding tumor-adjacent gastric tissues, with RPMM analysis computing the number of distinct subgroups was equal to the number of cases. These findings support the occurrence of non-random methylation events in tumorigenesis.
A total of 219 CpG sites (185 hypermethylated and 34 hypomethylated) in 147 unique genes were significantly differentially methylated between tumor and tumor-adjacent gastric tissue (FDR = 0.001), with the top three CpG sites located in HOXA5, SFRP1 and CCNA1 [Additional file 2]. Analysis by DAVID revealed that “Pathways in cancer”, the “Hedgehog signalling pathway” and “Cytokine-cytokine receptor interaction” were the top three significant pathways revealed by genes with tumor-specific methylation in GC (Table 1).
GC subgroups revealed by tumor-specific CpG methylation
Unsupervised clustering of the 219 tumor-specific CpG methylation sites by RPMM analysis revealed six GC clusters, denoted A-F (Figure 1). Bootstrap resampling analysis (n = 1,000) indicated two groups were the mode (91%) for the optimum number, hence the six clusters were consolidated into two major subtypes denoted as L (low methylation; clusters A-C; 28% of cases) and H (high methylation; clusters D-F; 72% of cases). In support of this classification, the mean methylation value (β-value) in subtype H tumors was twice that observed in subtype L (0.603 vs. 0.305, respectively; P < 0.001). Using a β-value threshold of ≥0.297 to binarize DNA methylation levels as described previously , H tumors also showed twice as many methylated CpG sites (180/219, 82%) compared to L tumors (89/219, 40%; P < 0.001). When RPMM was performed on the matched tumor-adjacent gastric mucosa using the tumor-specific CpG sites, the samples did not cluster according to the subtypes (L or H) of their corresponding paired tumors [Additional file 3]. This suggests the methylation patterns observed in tumors did not pre-exist in tumor-adjacent gastric tissue and were likely to result from somatic events.
Of the 219 CpG sites showing differential methylation between tumor-adjacent gastric and tumor tissue, 114 were also significantly different between the H and L tumor subgroups (FDR = 0.001) and all of these were hypermethylated in H but not L tumors [Additional file 4]. The three CpG sites showing the most significant difference in methylation level between H and L tumors were located in SEZ6L, FLT4 and ALK. Interestingly, the four most significant KEGG pathways identified by the differentially methylated genes between tumor and tumor-adjacent gastric tissues were also identified by the differentially methylated genes between H and L tumors (Table 1).
Associations between clinicopathological and molecular features and the H and L GC subtypes defined by methylation are shown in Figure 1 and Additional file 5. GC from female patients were almost all subtype H (16/17, 94%) and this was significantly higher than for male patients (27/43, 63%; P = 0.024). No significant associations were observed between the H and L subgroups and any of the other clinicopathological or molecular features of GC.
Methylation status and genome location, polycomb occupancy and histone modification
CpG sites that were differentially methylated between tumor and tumor-adjacent gastric tissue as well as between tumor subtypes H and L were referred to as Group HG (114 loci in 78 genes). CpG sites that were differentially methylated between tumor and tumor-adjacent gastric tissue, but not between tumor subtypes H and L, were referred to as Group LG (105 loci in 69 genes), while CpG sites not differentially methylated between tumor and tumor-adjacent gastric tissue were classified as Group NG (1,202 loci in 626 genes).
The distribution of CpG loci within CpG islands and promoters according to their annotation in the GoldenGate information sheet is displayed in Figure 2A for the NG, LG and HG groups. The proportion of CpG loci located within CpG islands increased progressively from NG (66%) to LG (71%) and HG (99%), with the differences between HG and LG (P < 0.001) and HG and NG (P < 0.001) being significant. This result shows that almost all differentially methylated CpG sites in the H group of GC were located within CpG islands.
Polycomb receptor complex-2 (PRC2) occupancy of gene promoters by the components Suz12 and Eed has been linked to susceptibility to methylation in human embryonic stem cells . In agreement with this, the proportion of genes with occupancy of both factors steadily increased from Groups NG (10%) to LG (19%) and HG (28%) (Figure 2B). The differences between groups NG and LG (P = 0.034) and between NG and HG (P <0.001) were significant.
The H3K4 and K3K27 trimethylation status in human embryonic stem cells has been shown to reflect modes of gene regulation in differentiated cells, namely constitutive expression (H3K4+ K3K27-), constitutive repression (H3K4-K3K27+), and bivalent, or “primed” for repression and expression (H3K4+ K3K27+) . A significantly higher proportion of Group HG genes showed H3K4+ K3K27+ bivalent marks (62%) compared to Group LG (39%, P = 0.008) and Group NG genes (19%, P < 0.001) (Figure 2C).
CpG sites methylated according to H. pyloristatus
None of the CpG sites investigated showed significantly different methylation between tumor samples of HP+ and HP- patients. However, comparison of the tumor-adjacent gastric tissue between HP+ and HP- patients revealed 8 differentially methylated CpG sites located within 7 genes (CCNA1, CSPG2, DAB2IP, DIO3, FLT1, STAT5A and TWIST1) [Additional file 6]. All 8 sites were hypomethylated in HP+ compared to HP- cases.
Verification of differential methylation
To verify the robustness of the observed differential methylation, pyroseqeuencing was performed to quantify methylation at the same CpG site of the top differentially methylated gene (HOXA5) and a gene from the top differentially methylated pathway (WNT5A, hedgehog signalling) in an independent series of tumor and matched tumor-adjacent gastric tissues from 60 subjects with GC. The higher level of methylation in tumors compared to tumor-adjacent gastric tissues that was observed in GoldenGate analysis was observed again for both HOXA5 (mean difference = 16.4%, P < 0.001 by paired t-test) and WNT5A (20.0%, P < 0.001) in this independent series (Additional file 7).
Candidate gene studies have so far identified 77 genes that are differentially methylated between normal and malignant gastric tissue . In the present work, the methylation of 1,421 autosomal CpG sites located within 768 cancer-related genes was evaluated in 60 pairs of GC and matched tumor-adjacent gastric tissue. A total of 219 CpG sites within 147 genes were found to be differentially methylated. Only 27/77 (35%) of the genes previously identified as being differentially methylated in the candidate gene studies were included in the GoldenGate methylation arrays used here. Hence, with the exception of 6 genes (CHFR, DAB2IP, DLC1, SFRP1, TCF4 and TFPI2), almost all of the 147 genes identified in the present study are novel methylation markers for GC that could be investigated further for potential roles in gastric tumorigenesis and for early screening.
GoldenGate methylation arrays have previously been used for comprehensive methylation studies of several cancer types including colorectal, head and neck, renal, breast and non-small cell lung cancer [10, 39, 45–47]. They have also been used to study non-cancerous gastric mucosa with respect to H. pylori infection and the presence of malignant tissue . However, the current study is the first to apply GoldenGate methylation arrays to investigate differential methylation between GC and matched tumor-adjacent gastric tissues. We believe that the high-quality tissues used in our study allows us to add value to the available scientific knowledge. Specifically, this unbiased, genome-wide approach revealed the existence of six methylation subgroups contained within two distinct clusters that comprised 28% (L) and 72% (H) of GC tumors (Figure 1). The mean methylation level of CpG loci in H tumors was twice that of L tumors. Compared to L tumors, H tumors were significantly over-represented with female patients (37% vs. 6%). Despite the small number of MSI cases (n = 4), a trend was also observed for more frequent association of the MSI phenotype with H tumors (9% vs. 0%).
The associations of female gender and MSI with the high methylation (H) subgroup of GC are in line with CIMP+ CRC, where these associations are reported consistently. The observation that hypermethylated CpG loci in the H subgroup are almost exclusively located in CpG islands (Figure 2A) is also consistent with the definition of CIMP and with previous findings in CIMP+ CRC [13–16, 49]. In general, however, the evidence in support of a distinct CIMP+ GC subgroup has so far been unconvincing. Highly methylated GC subgroups have shown contradictory associations with the stage of disease [18, 19, 25, 26, 50], histological type [17–19, 22, 50] and patient survival [17, 19, 20, 50]. Reported associations with older age , proximal tumor location  and poor differentiation  have not been confirmed by others, although more consistent associations have been reported with Epstein Barr virus infection [17, 18, 50], lymph node metastasis [25, 51] and MSI [20, 52].
A likely explanation for the inconsistent results to date on CIMP in GC is that studies have been limited to a small number of genes used different methylation assays and thresholds . Array-based studies have enabled a more comprehensive analysis of the DNA methylome. In addition to the present study using GoldenGate methylation arrays, Kim et al. recently published results using Infinium HumanMethylation450 BeadChip arrays that evaluate almost 500,000 CpG sites . They reported a CIMP+ subgroup in 11 of 30 (37%) GC samples analysed. These tumors displayed a higher frequency of oncogene mutations including KRAS and PIK3CA. Zouridis et al. also recently published their results on 203 GC obtained using Infinium HumanMethylation27 BeadChip arrays that evaluate 27,578 CpG sites . These workers reported a CIMP+ subgroup that comprised a similar proportion (35%) to that reported by Kim et al. but was considerably smaller than the H subgroup found in the current study (72%). The CIMP+ subgroup described by Zouridis et al. was also characterized by younger patient age and worse survival. A possible source of bias in our study was that all CpG sites evaluated by the Illumina GoldenGate Methylation Cancer Panel I array were within cancer-related genes. Moreover, only two or less CpG sites were evaluated for most (86%) of these genes. Together with the larger Illumina Infinium BeadChip arrays, next generation sequencing will enable even more comprehensive profiling of the methylome in GC. However, this does not necessarily imply that a distinctive CIMP+ GC subgroup will be identified in a consistent fashion.
HOXA5 was found here to be the most differentially methylated gene between gastric tumor and tumor-adjacent gastric tissue (Additional file 2), with the observation of a higher methylation level in tumor replicated in an independent series using a different experimental technique (P < 0.001; pyrosequencing). The HOX gene family, of which HOXA5 is a member, is known to play important roles in embryonic development and adult cell differentiation . HOXA5 is temporally expressed in the mesenchymal component of the developing gut  and a loss of HOXA5 function can perturb intestinal maturation in mice . Hypermethylation of HOXA5 has been reported in several cancer types [57–61] and is associated with decreased expression [59, 62]. The present study is the first to our knowledge to report HOXA5 methylation in GC. This finding warrants further functional studies to determine whether methylation-induced silencing of HOXA5 is a driver event for gastric tumorigenesis. The second most differentially methylated gene observed here, SFRP1, has previously been reported to be methylated in over 90% of primary GC .
Recursive partitioning identified a single CpG site within SEZ6L whose methylation status could differentiate the L and H GC subgroups. SEZ6L methylation has previously been reported in GC  and CRC . SEZ6L methylation in the gastric mucosa of non-GC subjects has also been associated with H. pylori infection . The role of this gene in gastric tumorigenesis is currently unclear, although it has been implicated as a risk factor for lung cancer [65, 66].
The associations between methylated genes, polycomb occupancy and H3K4/H3K27 modifications observed here for GC (Figure 2B and C) and elsewhere for other cancer types suggests that aberrations in chromatin regulation could underlie the hypermethylation phenotypes observed in cancer. The recent introduction of standardized methylation assay platforms with genome-wide coverage, such as the Illumina Infinium BeadChip arrays, should allow this area to be investigated in much more detail in future studies.
Aberrant methylation of gastric mucosa has been implicated in the elevated risk of GC in HP-infected individuals [67, 68]. Other studies have reported differentially methylated genes between HP- and HP+ GC tissue [48, 69]. The current study did not find any genes that were differentially methylated between HP- and HP+ GC tissue, but did find 8 genes that were hypomethylated in the tumor-adjacent gastric tissue of HP+ GC patients [Additional file 5]. Interestingly, amongst these 8 genes were DAB2IP and TWIST1, both of which have been implicated in gastric tumorigenesis [70, 71]. STAT5A was also previously observed to be hypomethylated in HP+ compared to HP- tissue from non-GC subjects , thus mirroring the present results in GC patients. CDH1, FLNC and HAND1 were previously reported to be methylated in HP+ GC tissues [67–69] using the GoldenGate array, but were not differentially methylated in the current study. This may due to the use of continuous versus binary values for methylation and to the thresholds used for statistical testing. Differential gene methylation in normal gastric mucosa between HP+ and HP- individuals may reflect the fact this pathogen is an initiating factor in the neoplastic transformation of gastric mucosa.
Methylation analysis of 1,421 CpG loci led to the identification of 147 genes that were differentially methylated in GC. Two major subgroups of GC were identified according to unbiased analysis of methylation levels. Tumors with high levels of methylation (subgroup H) shared some features consistent with CIMP in CRC. The methylation status of a single CpG site in SEZ6L was sufficient to allow absolute discrimination of the L and H subgroups of GC. Methylated genes in subtype H were characterized by higher frequencies of polycomb occupancy and H3K4+/H3K27+ bivalent marks, thus providing evidence of links between the hypermethylated phenotype and chromatin dysregulation. Further investigations of large population-based series are required to validate these findings and to assess the clinical utility of subgroups defined by methylation status.
Bayesian information criterion
CpG island methylator phenotype
False discovery rate
Formalin-fixed and paraffin-embedded
- H. pylori:
Recursively partitioned mixture model.
Yeoh KG: How do we improve outcomes for gastric cancer?. J Gastroenterol Hepatol. 2007, 22 (7): 970-972. 10.1111/j.1440-1746.2007.04956.x.
Chie WC, Chang YH, Chen HH: A novel method for evaluation of improved survival trend for common cancer: early detection or improvement of medical care. J Eval Clin Pract. 2007, 13 (1): 79-85. 10.1111/j.1365-2753.2006.00656.x.
Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D: Global cancer statistics. CA Cancer J Clin. 2011, 61 (2): 69-90. 10.3322/caac.20107.
Chen J, Bu XL, Wang QY, Hu PJ, Chen MH: Decreasing seroprevalence of Helicobacter pylori infection during 1993–2003 in Guangzhou, southern China. Helicobacter. 2007, 12 (2): 164-169. 10.1111/j.1523-5378.2007.00487.x.
Kawakami E, Machado RS, Ogata SK, Langner M: Decrease in prevalence of Helicobacter pylori infection during a 10-year period in Brazilian children. Arq Gastroenterol. 2008, 45 (2): 147-151.
Tkachenko MA, Zhannat NZ, Erman LV, Blashenkova EL, Isachenko SV, Isachenko OB, Graham DY, Malaty HM: Dramatic changes in the prevalence of Helicobacter pylori infection during childhood: a 10-year follow-up study in Russia. J Pediatr Gastroenterol Nutr. 2007, 45 (4): 428-432. 10.1097/MPG.0b013e318064589f.
Bertuccio P, Chatenoud L, Levi F, Praud D, Ferlay J, Negri E, Malvezzi M, La Vecchia C: Recent patterns in gastric cancer: a global overview. Int J Cancer. 2009, 125 (3): 666-673. 10.1002/ijc.24290.
Parkin DM: International variation. Oncogene. 2004, 23 (38): 6329-6340. 10.1038/sj.onc.1207726.
Parkin DM, Bray FI, Devesa SS: Cancer burden in the year 2000. The global picture. Eur J Cancer. 2001, 37 (Suppl 8): S4-S66.
Christensen BC, Kelsey KT, Zheng S, Houseman EA, Marsit CJ, Wrensch MR, Wiemels JL, Nelson HH, Karagas MR, Kushi LH, Kwan ML, Wiencke JK: Breast cancer DNA methylation profiles are associated with tumor size and alcohol and folate intake. PLoS Genet. 2010, 6 (7): e1001043-10.1371/journal.pgen.1001043.
Christensen BC, Marsit CJ, Houseman EA, Godleski JJ, Longacker JL, Zheng S, Yeh RF, Wrensch MR, Wiemels JL, Karagas MR, Bueno R, Sugarbaker DJ, Nelson HH, Wiencke JK, Kelsey KT: Differentiation of lung adenocarcinoma, pleural mesothelioma, and nonmalignant pulmonary tissues using DNA methylation profiles. Cancer Res. 2009, 69 (15): 6315-6321. 10.1158/0008-5472.CAN-09-1073.
Christensen BC, Smith AA, Zheng S, Koestler DC, Houseman EA, Marsit CJ, Wiemels JL, Nelson HH, Karagas MR, Wrensch MR, Kelsey KT, Wiencke JK: DNA methylation, isocitrate dehydrogenase mutation, and survival in glioma. J Natl Cancer Inst. 2010, 103 (2): 143-153.
Toyota M, Ahuja N, Ohe-Toyota M, Herman JG, Baylin SB, Issa JP: CpG island methylator phenotype in colorectal cancer. Proc Natl Acad Sci U S A. 1999, 96 (15): 8681-8686. 10.1073/pnas.96.15.8681.
Issa JP: CpG island methylator phenotype in cancer. Nat Rev Cancer. 2004, 4 (12): 988-993. 10.1038/nrc1507.
Ogino S, Goel A: Molecular classification and correlates in colorectal cancer. J Mol Diagn. 2008, 10 (1): 13-27. 10.2353/jmoldx.2008.070082.
Weisenberger DJ, Siegmund KD, Campan M, Young J, Long TI, Faasse MA, Kang GH, Widschwendter M, Weener D, Buchanan D, Koh H, Simms L, Barker M, Leggett B, Levine J, Kim M, French AJ, Thibodeau SN, Jass J, Haile R, Laird PW: CpG island methylator phenotype underlies sporadic microsatellite instability and is tightly associated with BRAF mutation in colorectal cancer. Nat Genet. 2006, 38 (7): 787-793. 10.1038/ng1834.
Enomoto S, Maekita T, Tsukamoto T, Nakajima T, Nakazawa K, Tatematsu M, Ichinose M, Ushijima T: Lack of association between CpG island methylator phenotype in human gastric cancers and methylation in their background non-cancerous gastric mucosae. Cancer Sci. 2007, 98 (12): 1853-1861. 10.1111/j.1349-7006.2007.00625.x.
Kusano M, Toyota M, Suzuki H, Akino K, Aoki F, Fujita M, Hosokawa M, Shinomura Y, Imai K, Tokino T: Genetic, epigenetic, and clinicopathologic features of gastric carcinomas with the CpG island methylator phenotype and an association with Epstein-Barr virus. Cancer. 2006, 106 (7): 1467-1479. 10.1002/cncr.21789.
Park SY, Kook MC, Kim YW, Cho NY, Jung N, Kwon HJ, Kim TY, Kang GH: CpG island hypermethylator phenotype in gastric carcinoma and its clinicopathological features. Virchows Arch. 2010, 457 (4): 415-422. 10.1007/s00428-010-0962-0.
An C, Choi IS, Yao JC, Worah S, Xie K, Mansfield PF, Ajani JA, Rashid A, Hamilton SR, Wu TT: Prognostic significance of CpG island methylator phenotype and microsatellite instability in gastric carcinoma. Clin Cancer Res. 2005, 11 (2 Pt 1): 656-663.
Oue N, Mitani Y, Motoshita J, Matsumura S, Yoshida K, Kuniyasu H, Nakayama H, Yasui W: Accumulation of DNA methylation is associated with tumor stage in gastric cancer. Cancer. 2006, 106 (6): 1250-1259. 10.1002/cncr.21754.
Oue N, Oshimo Y, Nakayama H, Ito R, Yoshida K, Matsusaki K, Yasui W: DNA methylation of multiple genes in gastric carcinoma: association with histological type and CpG island methylator phenotype. Cancer Sci. 2003, 94 (10): 901-905. 10.1111/j.1349-7006.2003.tb01373.x.
Kang GH, Lee S, Cho NY, Gandamihardja T, Long TI, Weisenberger DJ, Campan M, Laird PW: DNA methylation profiles of gastric carcinoma characterized by quantitative DNA methylation analysis. Lab Invest. 2008, 88 (2): 161-170. 10.1038/labinvest.3700707.
Kang GH, Lee S, Kim WH, Lee HW, Kim JC, Rhyu MG, Ro JY: Epstein-barr virus-positive gastric carcinoma demonstrates frequent aberrant methylation of multiple genes and constitutes CpG island methylator phenotype-positive gastric carcinoma. Am J Pathol. 2002, 160 (3): 787-794. 10.1016/S0002-9440(10)64901-2.
Chen HY, Zhu BH, Zhang CH, Yang DJ, Peng JJ, Chen JH, Liu FK, He YL: High CpG island methylator phenotype is associated with lymph node metastasis and prognosis in gastric cancer. Cancer Sci. 2011, 103 (1): 73-79.
Toyota M, Ahuja N, Suzuki H, Itoh F, Ohe-Toyota M, Imai K, Baylin SB, Issa JP: Aberrant methylation in gastric cancer associated with the CpG island methylator phenotype. Cancer Res. 1999, 59 (21): 5438-5442.
Sapari NS, Loh M, Vaithilingam A, Soong R: Clinical potential of DNA methylation in gastric cancer: a meta-analysis. PLoS One. 2012, 7 (4): e36275-10.1371/journal.pone.0036275.
Loh M, Liem N, Lim PL, Vaithilingam A, Cheng CL, Salto-Tellez M, Yong WP, Soong R: Impact of sample heterogeneity on methylation analysis. Diagn Mol Pathol. 2010, 19 (4): 243-247. 10.1097/PDM.0b013e3181de4396.
Soong R, Iacopetta BJ: A rapid and nonisotopic method for the screening and sequencing of p53 gene mutations in formalin-fixed, paraffin-embedded tumors. Mod Pathol. 1997, 10 (3): 252-258.
Davies H, Bignell GR, Cox C, Stephens P, Edkins S, Clegg S, Teague J, Woffendin H, Garnett MJ, Bottomley W, Davis N, Dicks E, Ewing R, Floyd Y, Gray K, Hall S, Hawes R, Hughes J, Kosmidou V, Menzies A, Mould C, Parker A, Stevens C, Watt S, Hooper S, Wilson R, Jayatilake H, Gusterson BA, Cooper C, Shipley J, et al: Mutations of the BRAF gene in human cancer. Nature. 2002, 417 (6892): 949-954. 10.1038/nature00766.
Krypuy M, Newnham GM, Thomas DM, Conron M, Dobrovic A: High resolution melting analysis for the rapid and sensitive detection of mutations in clinical samples: KRAS codon 12 and 13 mutations in non-small cell lung cancer. BMC Cancer. 2006, 6: 295-10.1186/1471-2407-6-295.
Buhard O, Cattaneo F, Wong YF, Yim SF, Friedman E, Flejou JF, Duval A, Hamelin R: Multipopulation analysis of polymorphisms in five mononucleotide repeats used to determine the microsatellite instability status of human tumors. J Clin Oncol. 2006, 24 (2): 241-251.
Bibikova M, Fan JB: GoldenGate assay for DNA methylation profiling. Methods Mol Biol. 2009, 507: 149-163. 10.1007/978-1-59745-522-0_12.
Houseman EA, Christensen BC, Yeh RF, Marsit CJ, Karagas MR, Wrensch M, Nelson HH, Wiemels J, Zheng S, Wiencke JK, Kelsey KT: Model-based clustering of DNA methylation array data: a recursive-partitioning algorithm for high-dimensional data arising as a mixture of beta distributions. BMC Bioinformatics. 2008, 9: 365-10.1186/1471-2105-9-365.
Hinoue T, Weisenberger DJ, Lange CP, Shen H, Byun HM, Van Den Berg D, Malik S, Pan F, Noushmehr H, van Dijk CM, Tollenaar RA, Laird PW: Genome-scale analysis of aberrant DNA methylation in colorectal cancer. Genome Res. 2011, 22 (2): 271-282.
Calinski RB, Harabasz J: A dendrite method for cluster analysis. Comm Stat. 1974, 3: 1-27.
Bibikova M, Lin Z, Zhou L, Chudin E, Garcia EW, Wu B, Doucet D, Thomas NJ, Wang Y, Vollmer E, Goldmann T, Seifart C, Jiang W, Barker DL, Chee MS, Floros J, Fan JB: High-throughput DNA methylation profiling using universal bead arrays. Genome Res. 2006, 16 (3): 383-393. 10.1101/gr.4410706.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995, 57: 289-300.
Ang PW, Loh M, Liem N, Lim PL, Grieu F, Vaithilingam A, Platell C, Yong WP, Iacopetta B, Soong R: Comprehensive profiling of DNA methylation in colorectal cancer reveals subgroups with distinct clinicopathological and molecular features. BMC Cancer. 2010, 10: 227-10.1186/1471-2407-10-227.
Hosack DA, Dennis G, Sherman BT, Lane HC, Lempicki RA: Identifying biological themes within lists of genes with EASE. Genome Biol. 2003, 4 (10): R70-10.1186/gb-2003-4-10-r70.
Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4 (5): 3-10.1186/gb-2003-4-5-p3.
Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M: KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 1999, 27 (1): 29-34. 10.1093/nar/27.1.29.
Lee TI, Jenner RG, Boyer LA, Guenther MG, Levine SS, Kumar RM, Chevalier B, Johnstone SE, Cole MF, Isono K, Koseki H, Fuchikami T, Abe K, Murray HL, Zucker JP, Yuan B, Bell GW, Herbolsheimer E, Hannett NM, Sun K, Odom DT, Otte AP, Volkert TL, Bartel DP, Melton DA, Gifford DK, Jaenisch R, Young RA: Control of developmental regulators by Polycomb in human embryonic stem cells. Cell. 2006, 125 (2): 301-313. 10.1016/j.cell.2006.02.043.
Kim TY, Jong HS, Jung Y, Kim TY, Kang GH, Bang YJ: DNA hypermethylation in gastric cancer. Aliment Pharmacol Ther. 2004, 20 (Suppl 1): 131-142.
Marsit CJ, Christensen BC, Houseman EA, Karagas MR, Wrensch MR, Yeh RF, Nelson HH, Wiemels JL, Zheng S, Posner MR, McClean MD, Wiencke JK, Kelsey KT: Epigenetic profiling reveals etiologically distinct patterns of DNA methylation in head and neck squamous cell carcinoma. Carcinogenesis. 2009, 30 (3): 416-422. 10.1093/carcin/bgp006.
McRonald FE, Morris MR, Gentle D, Winchester L, Baban D, Ragoussis J, Clarke NW, Brown MD, Kishida T, Yao M, Latif F, Maher ER: CpG methylation profiling in VHL related and VHL unrelated renal cell carcinoma. Mol Cancer. 2009, 8: 31-10.1186/1476-4598-8-31.
Moran A, Fernandez-Marcelo T, Carro J, De Juan C, Pascua I, Head J, Gomez A, Hernando F, Torres AJ, Benito M, Iniesta P: Methylation profiling in non-small cell lung cancer: clinical implications. Int J Oncol. 2011, 40 (3): 739-746.
Shin CM, Kim N, Jung Y, Park JH, Kang GH, Park WY, Kim JS, Jung HC, Song IS: Genome-wide DNA methylation profiles in noncancerous gastric mucosae with regard to Helicobacter pylori infection and the presence of gastric cancer. Helicobacter. 2011, 16 (3): 179-188. 10.1111/j.1523-5378.2011.00838.x.
Iacopetta B, Kawakami K, Watanabe T: Predicting clinical outcome of 5-fluorouracil-based chemotherapy for colon cancer patients: is the CpG island methylator phenotype the 5-fluorouracil-responsive subgroup?. Int J Clin Oncol. 2008, 13 (6): 498-503. 10.1007/s10147-008-0854-3.
Chang MS, Uozaki H, Chong JM, Ushiku T, Sakuma K, Ishikawa S, Hino R, Barua RR, Iwasaki Y, Arai K, Fujii H, Nagai H, Fukayama M: CpG island methylation status in gastric carcinoma with and without infection of Epstein-Barr virus. Clin Cancer Res. 2006, 12 (10): 2995-3002. 10.1158/1078-0432.CCR-05-1601.
Zhang KL, Sun Y, Li Y, Liu M, Qu B, Cui SH, Kong QY, Chen XY, Li H, Liu J: Increased frequency of CpG island methylator phenotype and CDH1 methylation in a gastric cancer high-risk region of china. Transl Oncol. 2008, 1 (1): 28-35. 10.1593/tlo.07103.
Kim H, Kim YH, Kim SE, Kim NG, Noh SH: Concerted promoter hypermethylation of hMLH1, p16INK4A, and E-cadherin in gastric carcinomas with microsatellite instability. J Pathol. 2003, 200 (1): 23-31. 10.1002/path.1325.
Kim JG, Takeshima H, Niwa T, Rehnberg E, Shigematsu Y, Yoda Y, Yamashita S, Kushima R, Maekita T, Ichinose M, Katai H, Park WS, Hong YS, Park CH, Ushijima T: Comprehensive DNA methylation and extensive mutation analyses reveal an association between the CpG island methylator phenotype and oncogenic mutations in gastric cancers. Cancer Lett. 2013, 330 (1): 33-40. 10.1016/j.canlet.2012.11.022.
Zouridis H, Deng N, Ivanova T, Zhu Y, Wong B, Huang D, Wu YH, Wu Y, Tan IB, Liem N, Gopalakrishnan V, Luo Q, Wu J, Lee M, Yong WP, Goh LK, Teh BT, Rozen S, Tan P: Methylation subtypes and large-scale epigenetic alterations in gastric cancer. Sci Transl Med. 2012, 4 (156)): 156ra140-
Byun HM, Siegmund KD, Pan F, Weisenberger DJ, Kanel G, Laird PW, Yang AS: Epigenetic profiling of somatic tissues from human autopsy specimens identifies tissue- and individual-specific DNA methylation patterns. Hum Mol Genet. 2009, 18 (24): 4808-4817. 10.1093/hmg/ddp445.
Berman DM, Karhadkar SS, Maitra A, Montes De Oca R, Gerstenblith MR, Briggs K, Parker AR, Shimada Y, Eshleman JR, Watkins DN, Beachy PA: Widespread requirement for Hedgehog ligand stimulation in growth of digestive tract tumours. Nature. 2003, 425 (6960): 846-851. 10.1038/nature01972.
Taniguchi H, Yamamoto H, Akutsu N, Nosho K, Adachi Y, Imai K, Shinomura Y: Transcriptional silencing of hedgehog-interacting protein by CpG hypermethylation and chromatic structure in human gastrointestinal cancer. J Pathol. 2007, 213 (2): 131-139. 10.1002/path.2216.
Vismodegib granted FDA approval for treatment of basal cell carcinoma. Oncology. 2012, 26 (2): 174-213.
Strathdee G, Sim A, Soutar R, Holyoake TL, Brown R: HOXA5 is targeted by cell-type-specific CpG island methylation in normal cells and during the development of acute myeloid leukaemia. Carcinogenesis. 2007, 28 (2): 299-309.
Ma X, Chen K, Huang S, Zhang X, Adegboyega PA, Evers BM, Zhang H, Xie J: Frequent activation of the hedgehog pathway in advanced gastric adenocarcinomas. Carcinogenesis. 2005, 26 (10): 1698-1705. 10.1093/carcin/bgi130.
Martin J, Donnelly JM, Houghton J, Zavros Y: The role of sonic hedgehog reemergence during gastric cancer. Dig Dis Sci. 2010, 55 (6): 1516-1524. 10.1007/s10620-010-1252-z.
Houghton J, Stoicov C, Nomura S, Rogers AB, Carlson J, Li H, Cai X, Fox JG, Goldenring JR, Wang TC: Gastric cancer originating from bone marrow-derived cells. Science. 2004, 306 (5701): 1568-1571. 10.1126/science.1099513.
Nojima M, Suzuki H, Toyota M, Watanabe Y, Maruyama R, Sasaki S, Sasaki Y, Mita H, Nishikawa N, Yamaguchi K, Hirata K, Itoh F, Tokino T, Mori M, Imai K, Shinomura Y: Frequent epigenetic inactivation of SFRP genes and constitutive activation of Wnt signaling in gastric cancer. Oncogene. 2007, 26 (32): 4699-4713. 10.1038/sj.onc.1210259.
Suzuki H, Gabrielson E, Chen W, Anbazhagan R, van Engeland M, Weijenberg MP, Herman JG, Baylin SB: A genomic screen for genes upregulated by demethylation and histone deacetylase inhibition in human colorectal cancer. Nat Genet. 2002, 31 (2): 141-149. 10.1038/ng892.
Nishioka M, Kohno T, Takahashi M, Niki T, Yamada T, Sone S, Yokota J: Identification of a 428-kb homozygously deleted region disrupting the SEZ6L gene at 22q12.1 in a lung cancer cell line. Oncogene. 2000, 19 (54): 6251-6260. 10.1038/sj.onc.1204031.
Raji OY, Agbaje OF, Duffy SW, Cassidy A, Field JK: Incorporation of a genetic factor into an epidemiologic model for prediction of individual risk of lung cancer: the Liverpool Lung Project. Cancer Prev Res (Phila). 2010, 3 (5): 664-669. 10.1158/1940-6207.CAPR-09-0141.
Maekita T, Nakazawa K, Mihara M, Nakajima T, Yanaoka K, Iguchi M, Arii K, Kaneda A, Tsukamoto T, Tatematsu M, Tamura G, Saito D, Sugimura T, Ichinose M, Ushijima T: High levels of aberrant DNA methylation in Helicobacter pylori-infected gastric mucosae and its possible association with gastric cancer risk. Clin Cancer Res. 2006, 12 (3 Pt 1): 989-995.
Nakajima T, Maekita T, Oda I, Gotoda T, Yamamoto S, Umemura S, Ichinose M, Sugimura T, Ushijima T, Saito D: Higher methylation levels in gastric mucosae significantly correlate with higher risk of gastric cancers. Canc Epidemiol Biomarkers Prev: a publication of the American Association for Cancer Research, cosponsored by the American Society of Preventive Oncology. 2006, 15 (11): 2317-2321. 10.1158/1055-9965.EPI-06-0436.
Chan AO, Lam SK, Wong BC, Wong WM, Yuen MF, Yeung YH, Hui WM, Rashid A, Kwong YL: Promoter methylation of E-cadherin gene in gastric mucosa associated with Helicobacter pylori infection and in gastric cancer. Gut. 2003, 52 (4): 502-506. 10.1136/gut.52.4.502.
Dote H, Toyooka S, Tsukuda K, Yano M, Ota T, Murakami M, Naito M, Toyota M, Gazdar AF, Shimizu N: Aberrant promoter methylation in human DAB2 interactive protein (hDAB2IP) gene in gastrointestinal tumour. Br J Cancer. 2005, 92 (6): 1117-1125. 10.1038/sj.bjc.6602458.
Schneider BG, Peng DF, Camargo MC, Piazuelo MB, Sicinschi LA, Mera R, Romero-Gallo J, Delgado AG, Bravo LE, Wilson KT, Peek RM, Correa P, El-Rifai W: Promoter DNA hypermethylation in gastric biopsies from subjects at high and low risk for gastric cancer. Int J Cancer. 2010, 127 (11): 2588-2597. 10.1002/ijc.25274.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-230X/14/55/prepub
This work is supported by a grant from the National Medical Research Council of Singapore (NMRC/TCR/001-NUH/2007).
The authors declare that they have no competing interests.
ML performed all statistical analysis and drafted the manuscript. NL, AV, PLL, NSS, EE, and ZUM carried out the experimental work. CLC, BY, BP, MST, and WPY contributed critical clinical resources, assessment and perspective. BI and RS co-ordinated the study and compiled the manuscript. All authors read and approved the final version of the manuscript.