Abstract
Aberrant methylations in human gastric cancers were searched for by a genome scanning technique, methylation-sensitive representational difference analysis. Nine CpG islands (CGIs) in the 5′ regions of nine genes, LOX, HRASLS, bA305P22.2.3, FLNc (γ-filamin/ABPL), HAND1, a homologue of RIKEN 2210016F16, FLJ32130, PGAR (HFARP/ANGPTL4/ARP4), and thrombomodulin, were found to be methylated in two gastric cancer cell lines, MKN28 and MKN74, but unmethylated in the normal sample. Their expressions were lost in the two cell lines, and the causal roles of the methylations in gene silencing were shown by treatment with 5-aza-2′-deoxycytidine. In 41 primary gastric cancers, methylations of these CGIs in association with reduced expressions were observed at high incidences (29–41%) for five genes, including LOX and HRASLS, which have been reported to have tumor-suppressive activities. The other four genes were methylated at low incidences (0–7%). By analysis of the numbers of aberrant methylations in individual cancers, a subset of cancers with high prevalence of aberrant methylations was identified, and all of the 11 cancers in the subset were diffuse type. To analyze the possible involvement of decreased fidelity of maintenance methylation in this subset of cancers, aberrant hypomethylations of three normally methylated CGIs were examined. Cancers with high prevalence of hypomethylations of normally methylated CGIs, however, constituted a different subset. It is expected that these nine genes may include important genes in gastric cancer development and would be useful to identify a distinct subset of gastric cancers.
INTRODUCTION
Epigenetic alterations, as well as genetic alterations, are involved in cancer development and progression (1, 2). Methylation of CGIs4 in the 5′ regions of tumor suppressor genes is known to inhibit transcriptional initiation and cause permanent silencing of the genes. Gene silencing is now regarded as one of the major mechanisms to inactivate tumor suppressor genes, along with gene mutations and deletions. On the other hand, methylation of CGIs in downstream exons and introns does not block elongation of mRNA in mammalian cells (1, 3). To identify genes involved in human cancers, such as novel tumor suppressor genes, it seems reasonable to search for aberrant methylations in the 5′ regions of genes as their markers. For this purpose, we previously established a genome-scanning method, MS-RDA (4). By this method, DNA fragments that are differentially methylated in two genomes can be isolated after amplification of DNA fragments obtained by digestion with a methylation-sensitive restriction enzyme and a series of subtractive hybridizations and selective amplifications.
Gastric cancer is one of the major current causes of cancer death in Asian countries (5). As for its molecular basis, silencing of the p16 gene has been reported in 42.0–42.2% of cancers (6, 7). Microsatellite instability is observed at incidences ranging between 31 and 67% (8, 9), and the majority of the microsatellite instability-positive gastric cancers show methylation of the hMLH1 gene (10, 11). RUNX3 was shown recently to be involved in human gastric cancers, and the major mechanism for its inactivation is also promoter methylation (12). The most prevalent genetic alteration in human gastric cancers is the p53 mutation, which is observed in 37–43% of intestinal-type gastric cancers (13). E-cadherin is inactivated by mutations in 17–56% of diffuse gastric cancers (14, 15) and by promoter methylation in a similar fraction (15).
In the present study, we applied the MS-RDA method to human gastric cancer cell lines using three restriction enzymes, HpaII, SacII, and NarI, and identified silencing of nine genes.
MATERIALS AND METHODS
Tissue Samples and Cell Lines.
Forty-one primary gastric cancer samples were obtained from 40 patients undergoing gastrectomy, with informed consent. Two independent gastric cancer samples (4C-1 and 4C-2) were obtained from case 4. For 28 of the 41 cancers, noncancerous gastric epithelial tissues (normal samples) were also obtained by scraping off the mucosae. These samples were frozen and stored at −80°C until extraction of DNA or RNA. Six gastric cancer cell lines, KATOIII, MKN28, MKN45, MKN74, NUGC3, and AGS, were obtained from the Japanese Collection of Research Bioresources (Tokyo, Japan) and American Type Culture Collection (Manassas, VA). Three gastric cancer cell lines, HSC39, HSC44, and HSC57, were established by K. Y. DNA was extracted by standard phenol/chloroform procedures. RNA was isolated with ISOGEN (Nippon Gene, Tokyo, Japan).
MS-RDA.
Three series of MS-RDA were performed using HpaII, with some modifications on the original method (4, 16), and two other methylation-sensitive restriction enzymes, SacII and NarI, respectively.
Genomic DNAs of a normal sample (sample 21N), MKN28 (established from an intestinal-type gastric cancer, 70-year-old female), and MKN74 (established from an intestinal-type gastric cancer, 37-year-old male) were digested with a methylation-sensitive restriction enzyme. An R adaptor specific to the restriction enzyme was ligated to the digestion product of the normal sample and to the mixture of digestion products of MKN28 and MKN74. An amplicon was prepared by PCR using an R primer in the presence of 1 m betaine (Sigma Chemical Co., St. Louis, MO), which facilitated amplification of G + C-rich sequences (17, 18). The R adaptor was then removed by digestion and gel filtration. For only the normal sample, a new J adaptor was ligated. The ligation product was mixed with an excess amount of the amplicon of the cancer cell lines and underwent denaturation and reannealing. Using the product, PCR was performed with a J primer, and DNA fragments present in only the amplicon of the normal sample were selectively amplified. After switching the J adaptor to an N adaptor, the second cycle of competitive hybridization and selective amplification was performed in the same manner. After the third competitive hybridization and selective amplification using the J adaptor again, the PCR product was cloned into the pGEM-T Easy vector (Promega Corp., Madison, WI). Ninety-six clones were sequenced and tested for independence. Adaptors for SacII and NarI were as follows: RSac26/9, AGC ACT CTC CAG CCT CTC ACG ACC GC/GGT CGT GAG; JSac26/9, ACC GAC GTC GAC TAT CCA TGA ACC GC/GGT TCA TGG; NSac26/9, AGG CAA CTG TGC TAT CCG AGG ACC GC/GGT CCT CGG; RNar24/11, AGC ACT CTC CAG GCA CTC ACC AGG/CGC CTG GTG AG; JNar24/11, ACC GAC GTC GAC TAT CCA TGA AGG/CGC CTT CAT GG; and NNar24/11, AGG CAA CTG TGC TAT CCG AGG AGG/CGC CTC CTC GG. Primers for SacII and NarI were as follows: RSac24, AGC ACT CTC CAG CCT CTC ACG ACC; JSac24, ACC GAC GTC GAC TAT CCA TGA ACC; NSac24, AGG CAA CTG TGC TAT CCG AGG ACC; RNar24, AGC ACT CTC CAG GCA CTC ACC AGG; JNar24, ACC GAC GTC GAC TAT CCA TGA AGG; and NNar24, AGG CAA CTG TGC TAT CCG AGG AGG.
DNA Sequencing and Database Search.
Cycle-sequencing was performed with a BigDye Terminator kit (PE Biosystems, Foster City, CA) and an ABI automated DNA sequencer (PE Biosystems). Repetitive sequences were screened by the Repeat Masker software (University of Washington Genome Center5), and sequence homology was searched using the BLASTN algorithm.
Sodium Bisulfite Modification, Bisulfite Sequencing, and MSP.
One μg of DNA, digested with BamHI, underwent sodium bisulfite modification as described previously (19) and was suspended in 20 μl of H2O. For bisulfite sequencing, 1 μl of the solution containing bisulfite-modified DNA was used for PCR with the primers common for methylated and unmethylated DNA sequences. The primers and PCR conditions are available as supplementary data. PCR products were cloned into pGEM-T Easy vector, and 10 clones from each sample were cycle-sequenced using T7 and Sp6 primers.
For MSP, 1 μl of the solution was used for PCR with primers specific for methylated (M) or unmethylated (U) sequences. Using the DNA from the normal sample and DNA that was methylated with SssI methylase, annealing temperatures specific for M and U primers were determined. Minimal PCR cycles to yield bands for these positive controls were determined, and four cycles were added for test samples. These primers and PCR conditions are available as supplementary data. MSP for p16 and hMLH1 were performed as described previously (11, 20). For normally unmethylated CGIs, a sample was regarded as methylation (−) when PCR product was obtained only with the U set and was regarded as methylation (+) when PCR product was obtained with the M set or with both the U set and M set. For NM-CGIs, a sample was regarded as demethylation (−) when PCR product was obtained only with the M set and was regarded as demethylation (+) when PCR product was obtained with the U set or with both the M set and U set.
Quantitative RT-PCR.
Total RNA was treated with DNase I (Ambion, Austin, TX), and cDNA was synthesized from 3 μg of total RNA using a Superscript II kit (Life Technologies, Inc., Rockville, MD). Real-time PCR was performed using SYBR Green PCR Core Reagents (PE Biosystems, Warrington, United Kingdom) and an iCycler Thermal Cycler (Bio-Rad Laboratories, Hercules, CA). The number of molecules of a specific gene in a sample was measured by comparing its amplification with the amplification of standard samples that contained 101 to 106 copies of the gene. The quantity of mRNA of each gene was normalized with that of PCNA. The primers and PCR conditions are available as supplementary data.
5-aza-dC Treatment.
Cells were seeded at a density of 3 × 105 cells/10 cm dish on day 0 and treated with freshly prepared 1 μm 5-aza-dC (Sigma) for 24 h on days 1, 3, and 5. After each treatment, the cells were placed in fresh medium and harvested on day 6.
RESULTS
MS-RDA.
DNA fragments aberrantly methylated in the gastric cancer cell lines MKN28 and MKN74 were searched for by three series of MS-RDA using three methylation-sensitive restriction enzymes, HpaII, SacII, and NarI, respectively. Of the 96 clones sequenced in each series, 36, 17, and 18 clones were found to be independent, respectively. A database search revealed that 17, 17, and 17, respectively, of them were derived from CGIs, and that 3, 7, and 6, respectively, were derived from CGIs in the 5′ regions of genes (Table 1).
Methylation Status of the 5′ Regions of 16 CGIs.
As for the 16 CGIs in the 5′ regions, their methylation status in the normal sample (21N) and two gastric cancer cell lines (MKN28 and MKN74) were analyzed by bisulfite sequencing. Nine CGIs, which were in the 5′ regions of LOX, HRASLS, bA305P22.2.3, FLNc (γ-filamin/ABPL), HAND1, a homologue of RIKEN 2210016F16, FLJ32130, PGAR (HFARP/ANGPTL4/ARP4), and TM genes, were methylated in the two cancer cell lines but not in the normal sample (Fig. 1). For all of these nine genes, all of the DNA molecules analyzed were methylated in the cancer cell lines, and it was indicated that no unmethylated allele was retained. Especially for the HRASLS gene, a cytosine at −292 had a C/G polymorphism (Fig. 1,B), and it was confirmed that both the C allele and G allele were methylated. MSP primers were developed for these nine CGIs to detect their methylation (Fig. 2).
As for the remaining seven CGIs, five CGIs at the 5′ regions of XPCT, pcHDP7, a homologue of Chimerin 1, DKFZp566K023, and KIAA0247 were unmethylated both in the normal sample and in the cancer cell lines, and two CGIs of RX and FHF-4 were methylated in both (data not shown).
Confirmation of Silencing for the Nine Genes.
Expressions of the nine genes were analyzed in the normal sample (21N) and nine gastric cancer cell lines, including MKN28 and MKN74, by quantitative RT-PCR (Fig. 3). Their expressions in MKN28 and MKN74 were reduced to less than 1/10, generally less than 1/100, of 21N, whereas their expressions were retained in other gastric cancer cell lines without methylation of the CGIs in the 5′ regions.
To examine the role of methylation of the CGIs in gene silencing, MKN28 and MKN74 were treated with 5-aza-dC, a demethylating agent. Demethylation of the 5′ regions of the nine genes was observed in both cell lines (Fig. 2), and reexpression was induced (Fig. 3). These data showed that the methylation of the 5′ regions of the nine genes caused their silencing.
Presence of Methylation and Reduced Expression in Primary Gastric Cancers.
Using 41 primary gastric cancers and nine gastric cancer cell lines, MSP was performed to examine the methylation status of the nine CGIs identified here and two CGIs in the promoter regions of p16 and hMLH1 (Figs. 2 and 4). The CGIs in the 5′ regions of FLNc, TM, HRASLS, HAND1, and LOX were methylated at high incidences (29–41% in primary cancers), whereas those of the other four genes were methylated at low incidences (0–7%). As for the normal samples, methylation was detected in 0–5 of the 28 samples (Fig. 4).
Expression levels of the nine genes, as well as p16 and hMLH1, were examined in the 41 primary gastric cancers and eight normal samples by quantitative RT-PCR (Fig. 5). Expression levels in cancer samples with methylation were generally reduced, but each gene had a unique expression profile. As for LOX, HRASLS, HAND1, PGAR, and TM, both overexpression and reduced expression were observed in the cancers without methylation, and reduced expression was observed in the majority of cancers with methylation. This pattern was similar to p16. As for bA305P.22.2.3 and FLJ32130, reduced expression was observed in the cancers without methylation, and the reduction was clearer in the cancers with methylation. This pattern was similar to hMLH1. A homologue of RIKEN 220016F16 showed a similar tendency, although primary cancers with methylation were not available. As for FLNc, remarkable overexpression was observed in some cancers, both without methylation and with methylation.
Methylation and Hypomethylation Profile in the Primary Gastric Cancers.
When the number of aberrant methylations in individual cancers was analyzed using the nine CGIs identified in this study, 11 cancers had aberrant methylations in four or more CGIs, whereas 16 cancers had no aberrant methylation at all (Fig. 4). Aberrant methylations of p16 and/or hMLH1 were detected in 10 of the 11 cancers with high prevalence of aberrant methylations but in 2 of the 16 cancers with no aberrant methylation (P < 0.0001, χ2 test). All of the 11 cancers with high prevalence of aberrant methylations were diagnosed as diffuse-type gastric cancers, whereas 6 of the 16 cancers with no aberrant methylation were diffuse type (P < 0.001, χ2 test).
To examine the possible presence of aberrant hypomethylation of the NM-CGIs in the 5′ regions, MAGE-A1, MAGE-A3, and MAGE-B2 (21, 22, 23) were selected. A high prevalence of hypomethylation of the NM-CGIs was observed in a different subset of gastric cancers (Figs. 2 and 4).
DISCUSSION
Nine CGIs in the 5′ regions of nine genes were identified by MS-RDA as being completely methylated in two gastric cancer cell lines, MKN28 and MKN74 (Fig. 1). The nine genes were silenced in both cancer cell lines by the CGI methylation, but the normal sample and other unmethylated cancer cell lines retained expression (Figs. 2 and 3). The gene expressions in methylated cancer cell lines were not detected at all, or at very low levels less than 1/100 of cell lines without methylation. The residual expression could be derived from a tiny fraction of cells whose methylation was removed spontaneously, or leaky expression could take place even if the promoter region was methylated.
Five of the nine genes were methylated at high incidences (29–41%) in 41 primary cancers, and two of the five, LOX and HRASLS, were previously reported to have tumor-suppressive activities. LOX is known to be expressed in NIH3T3 cells, down-regulated after the transformation of the cells, but reexpressed in revertants (24). A selective inactivation of LOX caused transformation of rat fibroblasts (25). HRASLS is a human homologue of mouse A-C1, which has been reported to inhibit growth of ras-transformed NIH3T3 cells (26, 27). Although the roles of LOX and HRASLS in human gastric carcinogenesis still need to be studied, their methylations in cancers were first demonstrated here.
FLNc is a member of the filamin family, which is known to organize actin polymerization in response to various signals (28). In our previous study (29), a subunit of Arp2/3 complex, p41-Arc, was found to be reduced in human gastric cancers partly because of promoter methylation, and this complex is known to join the filamin family to regulate actin polymerization (28). These findings suggest that abrogation of normal actin polymerization may be important in gastric carcinogenesis. HAND1 encodes a basic helix-loop-helix transcriptional factor, which is essential for placental development and cardiac morphogenesis (30). Its expression in normal gastric mucosa and silencing in gastric cancers suggested that HAND1 might have a role in the maintenance of differentiated status of gastric epithelium. TM encodes an endothelial cell receptor that binds thrombin to activate protein C and works as a member of anticoagulant pathway (31). The significance of its expression in gastric epithelium and its silencing in cancers was unclear.
Analysis of primary cancers is difficult because contamination of normal cells in cancer tissues is inevitable, and “cancers with methylation” still have a chance to retain a normal allele. Therefore, the expressions of LOX, HRASLS, bA305P22.2.3, HAND1, FLJ32130, PGAR, and TM can be regarded as reduced in association with their methylation, as in the case of p16 and hMLH1 (Fig. 5). This showed that the silencing found in the cell lines were also present in the primary cancers. LOX, HRASLS, HAND1, PGAR, and TM showed both overexpression and reduced expression in cancers without methylation, similar to p16. An increase in the number of cells in cell cycle causes overexpression of p16, whereas loss of p16 contributes to tumorigenesis (32, 33). It was speculated that some genes with this type of expression profile, such as p16, might have roles in the negative regulation of cell growth. In contrast, bA305P22.2.3, a homologue of RIKEN 2210016F16 and FLJ32130, similar to hMLH, showed only reduced expression in cancers without methylation. It was speculated that genes in this group may not be directly involved in the regulation of cell growth but that some of them might be involved in carcinogenesis, e.g., by maintaining genome stability.
Eleven of the 41 primary cancers showed aberrant methylations frequently (4 CGIs or more), and a significantly high incidence of p16 and hMLH1 methylations was observed in this subset of cancers (Fig. 4). This was in good agreement with a previous report that the “CGI methylator phenotype (CIMP)” was present in 41% of human gastric cancers (34). CIMP is known to be associated with p16 and hMLH1 methylations and with loss of imprinting of IGF2 (35). It is noteworthy, however, that the 11 cancers with high prevalence of aberrant methylations here were significantly associated with diffuse-type gastric cancers, although the original MS-RDA was performed using cell lines from intestinal-type cancers. It could be interpreted as: (a) genetic or epigenetic instability underlying CIMP led to the poorly differentiated histology; and/or (b) some of the methylations that had been considered as “aberrant” reflected the methylations that were present in the precursor cells of diffuse-type gastric cancers.
The molecular basis underlying CIMP is unknown. It was hypothesized that if decreased fidelity of maintenance methylation was involved in CIMP, aberrant hypomethylation would be also observed in CIMP(+) cancers. We examined aberrant hypomethylation of three NM-CGIs at the 5′ regions of three MAGE genes (21, 22, 23). Five gastric cancers showed hypomethylation of all of the three but constituted a subset different from CIMP(+) cancers. Unusual clustering of hypomethylations was observed, and it was significant (P < 0.05) by a method reported previously (34). This indicated that a new phenotype, which could be designated as “CGI hypomethylator phenotype (CHOP),” is present and that CIMP and CHOP are caused by mechanisms other than decreased fidelity of maintenance methylation. Interestingly, four of five CHOP(+) gastric cancers were intestinal type, whereas 5 of 17 cancers with no aberrant hypomethylation were intestinal type (P < 0.05, χ2 test). MAGE products are known as carcinotesticular antigens recognized by autologous cytotoxic T cells, which are used for cancer immunotherapy (36, 37). Therefore, there is a possibility that CHOP(+) cancers constitute a group of cancers where immunotherapy is effective.
MS-RDA was performed using a mixture of DNAs of two cancer cell lines as the driver. This strategy might have neglected infrequent methylations that occurred in one cell line only, because MS-RDA recognizes differentially methylated regions by selectively amplifying unmethylated fragments, and thus MS-RDA preferentially isolates DNA fragments methylated almost completely in the driver.6 But this strategy was effective to isolate aberrant methylations that commonly occurred in the two cell lines, and thus possibly in other cell lines and primary cancers as well. However, confirmation of aberrant methylation in primary cancers is still necessary to exclude methylation restricted to cell lines.
In the present study, we found that silencing of nine genes is present in human gastric cancers, that gastric cancers with high prevalence of aberrant methylations are associated with diffuse-type histology, and that a new phenotype, CHOP, exists.
Genomic structures of the nine CGIs identified, isolated MS-RDA clones, and their methylation status. Densities of GpC, CpG, and recognition sites for the restriction enzyme used for MS-RDA are shown. For each gene, methylation status of the 5′ CGI, from which the MS-RDA clone was derived, was determined by bisulfite sequencing in the normal sample and two gastric cancer cell lines. For each sample, 10 clones were analyzed for methylation statuses. Open box, clone isolated by MS-RDA; closed circle, methylated CpG site; open circle, unmethylated CpG site. As for HRASLS, the cytosine at the 17th CpG site had a C/G polymorphism in the two cancer cell lines, and biallelic methylation of the 5′ region was clearly demonstrated.
Genomic structures of the nine CGIs identified, isolated MS-RDA clones, and their methylation status. Densities of GpC, CpG, and recognition sites for the restriction enzyme used for MS-RDA are shown. For each gene, methylation status of the 5′ CGI, from which the MS-RDA clone was derived, was determined by bisulfite sequencing in the normal sample and two gastric cancer cell lines. For each sample, 10 clones were analyzed for methylation statuses. Open box, clone isolated by MS-RDA; closed circle, methylated CpG site; open circle, unmethylated CpG site. As for HRASLS, the cytosine at the 17th CpG site had a C/G polymorphism in the two cancer cell lines, and biallelic methylation of the 5′ region was clearly demonstrated.
MSP analysis of the 5′ CGIs of the nine genes identified (A–I) and three MAGE genes (J–L). 1N, 5N, 6N, 17N, 18N, 21N, and 28N, normal samples. 1C, 5C, 6C, 17C, 18C, and 28C, cancer samples. Aza, treated with 5-aza-dC. SssI, genomic DNA of 21N treated with SssI methylase as a positive control for methylated DNA. A–I, specific amplification of unmethylated DNA by U primers and that of methylated DNA by M primers were confirmed by observing a band with U primers in 21N and a band with M primers in SssI-treated DNA. The 5′ CGIs of the nine genes were methylated in both MKN28 and MKN74, and their demethylation was induced by treatment with 5-aza-dC. Representative analysis of a cell line without methylation and primary samples with or without methylation are also shown. J–L, the 5′ CGIs of the three MAGE genes were methylated in normal samples (17N and 18N). Demethylation (−) in 17C and demethylation (+) in 18C were observed for all of the three genes.
MSP analysis of the 5′ CGIs of the nine genes identified (A–I) and three MAGE genes (J–L). 1N, 5N, 6N, 17N, 18N, 21N, and 28N, normal samples. 1C, 5C, 6C, 17C, 18C, and 28C, cancer samples. Aza, treated with 5-aza-dC. SssI, genomic DNA of 21N treated with SssI methylase as a positive control for methylated DNA. A–I, specific amplification of unmethylated DNA by U primers and that of methylated DNA by M primers were confirmed by observing a band with U primers in 21N and a band with M primers in SssI-treated DNA. The 5′ CGIs of the nine genes were methylated in both MKN28 and MKN74, and their demethylation was induced by treatment with 5-aza-dC. Representative analysis of a cell line without methylation and primary samples with or without methylation are also shown. J–L, the 5′ CGIs of the three MAGE genes were methylated in normal samples (17N and 18N). Demethylation (−) in 17C and demethylation (+) in 18C were observed for all of the three genes.
Expression of the nine genes in the normal sample, an unmethylated gastric cancer cell line, and MKN28 and MKN74 before (−) and after (+) 5-aza-dC treatment. Expression levels (shown in the logarithmic scale) were analyzed by quantitative RT-PCR and normalized with that of PCNA. The methylation status of each sample is shown in Fig. 2. For each of the nine genes, the expression in MKN28 and MKN74 was reduced to less than 1/10, generally less than 1/100, of 21N, whereas the expression was retained in a gastric cancer cell line without methylation. Treatment of MKN28 and MKN74 with 5-aza-dC resulted in reexpression of a gene to a level comparable with the expression level of a cell line without methylation.
Expression of the nine genes in the normal sample, an unmethylated gastric cancer cell line, and MKN28 and MKN74 before (−) and after (+) 5-aza-dC treatment. Expression levels (shown in the logarithmic scale) were analyzed by quantitative RT-PCR and normalized with that of PCNA. The methylation status of each sample is shown in Fig. 2. For each of the nine genes, the expression in MKN28 and MKN74 was reduced to less than 1/10, generally less than 1/100, of 21N, whereas the expression was retained in a gastric cancer cell line without methylation. Treatment of MKN28 and MKN74 with 5-aza-dC resulted in reexpression of a gene to a level comparable with the expression level of a cell line without methylation.
Distribution of the aberrant methylations and hypomethylations in primary gastric cancers (A) and gastric cancer cell lines (B). Methylation status was analyzed for the nine genes identified, p16, hMLH1, and three MAGE genes. bA305P, bA305P22.2.3; RIKEN, a homologue of RIKEN 2210016F16; dif, diffuse-type gastric cancer; int, intestinal-type gastric cancer. Closed boxes for the left 11 genes, which are normally unmethylated, denote aberrant methylations in cancers. Shaded boxes in the right three MAGE genes, which are normally methylated, denote aberrant hypomethylations in cancers. The incidences of aberrant methylation or hypomethylation of each gene in the 28 normal samples are shown at the bottom of A. The 5′ regions of FLNc, TM, HRASLS, HAND1, and LOX were methylated at high incidences, ranging from 12 to 17 (29–41%) of 41 primary cancers. The 5′ regions of bA305P22.2.3, PGAR, FLJ32130, and a homologue of RIKEN 2210016F16 were methylated at lower incidences, ranging from 0 to 3 (0–7%) of 41 primary cancers. A subset of 11 cancers had aberrant methylations in four of the nine genes or more. The subset showed significantly high incidences of methylation of p16 and hMLH1 and was associated with diffuse-type histology. A high prevalence of hypomethylation was observed in a different subset of gastric cancers.
Distribution of the aberrant methylations and hypomethylations in primary gastric cancers (A) and gastric cancer cell lines (B). Methylation status was analyzed for the nine genes identified, p16, hMLH1, and three MAGE genes. bA305P, bA305P22.2.3; RIKEN, a homologue of RIKEN 2210016F16; dif, diffuse-type gastric cancer; int, intestinal-type gastric cancer. Closed boxes for the left 11 genes, which are normally unmethylated, denote aberrant methylations in cancers. Shaded boxes in the right three MAGE genes, which are normally methylated, denote aberrant hypomethylations in cancers. The incidences of aberrant methylation or hypomethylation of each gene in the 28 normal samples are shown at the bottom of A. The 5′ regions of FLNc, TM, HRASLS, HAND1, and LOX were methylated at high incidences, ranging from 12 to 17 (29–41%) of 41 primary cancers. The 5′ regions of bA305P22.2.3, PGAR, FLJ32130, and a homologue of RIKEN 2210016F16 were methylated at lower incidences, ranging from 0 to 3 (0–7%) of 41 primary cancers. A subset of 11 cancers had aberrant methylations in four of the nine genes or more. The subset showed significantly high incidences of methylation of p16 and hMLH1 and was associated with diffuse-type histology. A high prevalence of hypomethylation was observed in a different subset of gastric cancers.
Expression analysis of the nine genes, p16 and hMLH1. Quantitative RT-PCR was performed using 41 primary gastric cancers and eight normal samples. The expression levels were normalized with that of PCNA and shown by dots with their means; bars, SE. N, normal samples. C(U), cancer samples without methylation. C(M), cancer samples with methylation. Numbers of samples were shown at the bottom of each graph. In C(M), expression levels were generally reduced for the nine of the 11 genes. In C(U), expression levels of LOX, HRASLS, HAND1, PGAR, and TM, similar to p16, ranged from overexpression to reduced expression. Expression levels of bA305P.22.2.3, a homologue of RIKEN 220016F16, and FLJ32130 were consistently reduced, similar to hMLH1. Overexpression of FLNc was observed in some cancer samples.
Expression analysis of the nine genes, p16 and hMLH1. Quantitative RT-PCR was performed using 41 primary gastric cancers and eight normal samples. The expression levels were normalized with that of PCNA and shown by dots with their means; bars, SE. N, normal samples. C(U), cancer samples without methylation. C(M), cancer samples with methylation. Numbers of samples were shown at the bottom of each graph. In C(M), expression levels were generally reduced for the nine of the 11 genes. In C(U), expression levels of LOX, HRASLS, HAND1, PGAR, and TM, similar to p16, ranged from overexpression to reduced expression. Expression levels of bA305P.22.2.3, a homologue of RIKEN 220016F16, and FLJ32130 were consistently reduced, similar to hMLH1. Overexpression of FLNc was observed in some cancer samples.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
This study was supported by a Grant-in-Aid from the Ministry of Education, Science, Culture and Sport; by a Grant-in-Aid for Human Genome and Tissue Regeneration from the Ministry of Health, Labor and Welfare; and by Public Trust Haraguchi Memorial Cancer Research Fund.
Supplementary data for this article is available at Cancer Research Online (http://cancerres.aacrjournals.org).
The abbreviations used are: CGI, CpG island; NM-CGI, normally methylated CGI; MS-RDA, methylation-sensitive representational difference analysis; MSP, methylation-specific PCR; RT-PCR, reverse transcription-PCR; PCNA, proliferating cell nuclear antigen; LOX, lysyl oxidase; HRASLS, HRAS-like suppressor; FLNc, filamin C; TM, thrombomodulin; PGAR, peroxisome proliferator-activated receptor-γ angiopoietin related; 5-aza-dC, 5-aza-2′-deoxycytidine; CIMP, CpG island methylator phenotype; CHOP, CpG island hypomethylator phenotype.
Internet address: http://ftp.genome.washington.edu/cgi-bin/RepeatMasker.
A. Kaneda, D. Takai, M. Kaminishi, E. Okochi, and T. Ushijima. Methylation-sensitive-representational difference analysis and its application to cancer research, submitted for publication.
MS-RDA clones derived from CpG island in the 5′ region of a gene
Clone no. . | Length (bp) . | G + C content . | CpG scorea . | Accession no. . | Chromosome . | Gene . | ||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
MS-RDA with HpaII | ||||||||||||
C2 | 388 | 65.7% | 0.29 | AL121673 | 20 | bA305P22.2.3 | ||||||
G8 | 410 | 57.3% | 0.64 | AC068162 | 3 | HRAS-like suppressor (HRASLS) | ||||||
G11 | 338 | 62.4% | 0.82 | AC090062 | 5 | Lysyl oxidase (LOX) | ||||||
MS-RDA with SacII | ||||||||||||
A1 | 373 | 67.8% | 0.81 | AC026688 | 5 | Basic helix-loop-helix protein HAND1 (HAND1) | ||||||
A2 | 328 | 70.1% | 0.71 | AC067859 | 18 | Retinal homeobox protein (RX) | ||||||
A9 | 367 | 63.5% | 0.83 | AL354733 | 9 | A homologue of RIKEN 2210016F16 | ||||||
A10 | 903 | 58.4% | 0.65 | AC004073 | X | Membrane transporter gene XPCT | ||||||
B3 | 456 | 63.6% | 0.77 | AC008063 | Liver dipeptidyl peptidase (pcHDP7) | |||||||
F2 | 360 | 53.1% | 0.79 | AL136965 | 13 | Fibroblast growth factor homologue (FHF-4) | ||||||
H3 | 827 | 66.9% | 0.75 | AF252549 | 7 | Filamin C (FLNc/γ-filamin/ABPL) | ||||||
MS-RDA with NarI | ||||||||||||
C12 | 353 | 68.3% | 0.88 | AL049651 | 20 | Thrombomodulin (TM) | ||||||
E12 | 864 | 60.0% | 0.57 | AC002310 | 16 | FLJ32130 | ||||||
F4 | 485 | 68.3% | 0.98 | AC074291 | 2 | A homologue of Chimerin 1 | ||||||
F11 | 893 | 66.1% | 0.81 | AC010323 | 19 | PPARγ angiopoietin related (PGAR/HFARP/ANGPTL4/ARP4) | ||||||
G6 | 612 | 62.1% | 0.70 | AL356152 | 1 | DFKZp566K023 | ||||||
H8 | 783 | 54.3% | 0.71 | AL133445 | 14 | KIAA0247 |
Clone no. . | Length (bp) . | G + C content . | CpG scorea . | Accession no. . | Chromosome . | Gene . | ||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
MS-RDA with HpaII | ||||||||||||
C2 | 388 | 65.7% | 0.29 | AL121673 | 20 | bA305P22.2.3 | ||||||
G8 | 410 | 57.3% | 0.64 | AC068162 | 3 | HRAS-like suppressor (HRASLS) | ||||||
G11 | 338 | 62.4% | 0.82 | AC090062 | 5 | Lysyl oxidase (LOX) | ||||||
MS-RDA with SacII | ||||||||||||
A1 | 373 | 67.8% | 0.81 | AC026688 | 5 | Basic helix-loop-helix protein HAND1 (HAND1) | ||||||
A2 | 328 | 70.1% | 0.71 | AC067859 | 18 | Retinal homeobox protein (RX) | ||||||
A9 | 367 | 63.5% | 0.83 | AL354733 | 9 | A homologue of RIKEN 2210016F16 | ||||||
A10 | 903 | 58.4% | 0.65 | AC004073 | X | Membrane transporter gene XPCT | ||||||
B3 | 456 | 63.6% | 0.77 | AC008063 | Liver dipeptidyl peptidase (pcHDP7) | |||||||
F2 | 360 | 53.1% | 0.79 | AL136965 | 13 | Fibroblast growth factor homologue (FHF-4) | ||||||
H3 | 827 | 66.9% | 0.75 | AF252549 | 7 | Filamin C (FLNc/γ-filamin/ABPL) | ||||||
MS-RDA with NarI | ||||||||||||
C12 | 353 | 68.3% | 0.88 | AL049651 | 20 | Thrombomodulin (TM) | ||||||
E12 | 864 | 60.0% | 0.57 | AC002310 | 16 | FLJ32130 | ||||||
F4 | 485 | 68.3% | 0.98 | AC074291 | 2 | A homologue of Chimerin 1 | ||||||
F11 | 893 | 66.1% | 0.81 | AC010323 | 19 | PPARγ angiopoietin related (PGAR/HFARP/ANGPTL4/ARP4) | ||||||
G6 | 612 | 62.1% | 0.70 | AL356152 | 1 | DFKZp566K023 | ||||||
H8 | 783 | 54.3% | 0.71 | AL133445 | 14 | KIAA0247 |
CpG score = (number of CpG) × (length)/(number of C) × (number of G).
Acknowledgments
We thank Dr. Y. Nakanishi for help in collecting tumor samples and histological examination.