Open Access Research article

Genetic lineages of undifferentiated-type gastric carcinomas analysed by unsupervised clustering of genomic DNA microarray data

Ayano Sonoda12, Ken-ichi Mukaisho1, Takahisa Nakayama1, Vo Thi Ngoc Diem1, Takanori Hattori1, Akira Andoh2, Yoshihide Fujiyama2 and Hiroyuki Sugihara1*

Author Affiliations

1 Department of Pathology, Division of Molecular and Diagnostic Pathology, Shiga University of Medical Science, Otsu 520-2192, Japan

2 Department of Internal Medicine, Division of Gastroenterology and Hematology, Shiga University of Medical Science, Otsu 520-2192, Japan

For all author emails, please log on.

BMC Medical Genomics 2013, 6:25  doi:10.1186/1755-8794-6-25

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1755-8794/6/25


Received:20 December 2012
Accepted:11 July 2013
Published:19 July 2013

© 2013 Sonoda et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

It is suspected that early gastric carcinoma (GC) is a dormant variant that rarely progresses to advanced GC. We demonstrated that the dormant and aggressive variants of tubular adenocarcinomas (TUBs) of the stomach are characterized by loss of MYC and gain of TP53 and gain of MYC and/or loss of TP53, respectively. The aim of this study is to determine whether this is also the case in undifferentiated-type GCs (UGCs) of different genetic lineages: one with a layered structure (LS+), derived from early signet ring cell carcinomas (SIGs), and the other, mostly poorly differentiated adenocarcinomas, without LS but with a minor tubular component (TC), dedifferentiated from TUBs (LS−/TC+).

Methods

Using 29 surgically resected stomachs with 9 intramucosal and 20 invasive UGCs (11 LS+ and 9 LS−/TC+), 63 genomic DNA samples of mucosal and invasive parts and corresponding reference DNAs were prepared from formalin-fixed, paraffin-embedded tissues with laser microdissection, and were subjected to array-based comparative genomic hybridization (aCGH), using 60K microarrays, and subsequent unsupervised, hierarchical clustering. Of 979 cancer-related genes assessed, we selected genes with mean copy numbers significantly different between the two major clusters.

Results

Based on similarity in genomic copy-number profile, the 63 samples were classified into two major clusters. Clusters A and B, which were rich in LS+ UGC and LS−/TC+ UGC, respectively, were discriminated on the basis of 40 genes. The aggressive pattern was more frequently detected in LS−/TC+ UGCs, (20/26; 77%), than in LS+UGCs (17/37; 46%; P = 0.0195), whereas no dormant pattern was detected in any of the UGC samples.

Conclusions

In contrast to TUBs, copy number alterations of MYC and TP53 exhibited an aggressive pattern in LS+ SIG at early and advanced stages, indicating that early LS+ UGCs inevitably progress to an advanced GC. Cluster B (enriched in LS−/TC+) exhibited more frequent gain of driver genes and a more frequent aggressive pattern than cluster A, suggesting potentially worse prognosis in UGCs of cluster B.

Background

Gastric carcinoma (GC) have been classified histologically into intestinal, diffuse and unclassified types by Lauren [1] and the unclassified type was further divided into solid and mixed types by Carneiro [2]. The undifferentiated-type gastric carcinoma (UGC) according to the Japanese classification [3] mostly overlaps poorly differentiated GC, which comprises not only the diffuse type including signet ring cell carcinoma (SIG) but also the solid type and the mixed type with minor tubular component (TC).

Recently it has been proposed that advanced diffuse-type GC may derive from either early diffuse-type or intestinal-type GC. Well differentiated tubular adenocarcinoma (TUB) can transform into poorly differentiated adenocarcinoma (POR) after the silencing of cell adhesion-related genes including CDH1[4,5]. Carneiro’s mixed type carcinomas may thus overlap dedifferentiated TUBs. It has been reported that the survival rate of the patients with mixed-type GCs was significantly lower than that of the patients with GCs of other types [2], whereas the survival rate of early GC patients with SIG was higher than that of GC patients without SIG [6]. Thus UGCs may be divided into subgroups with different prognosis. Recently, a mass-screening program for neuroblastomas [7-9] was suspended in Japan because a discontinuous genetic lineage was observed between the early- and the late-presenting neuroblastomas. Negative and late-presenting (≥1 year) neuroblastomas exhibited near-diploidy with terminal 1p deletion, whereas positive neuroblastomas in infants exhibited near-triploidy without 1p deletion [10,11]. To perform such subgrouping, we have classified UGCs based on the continuity of genetic lineages as well as the expression of morphological lineage markers.

Our lineage analysis using chromosomal comparative genomic hybridization (CGH) was based on distinctive morphological lineage markers. A layered structure (LS) represents an incipient phase of SIG development [12] and is commonly retained even at an advanced stage in the human stomach. In tumour regions with LS, the mode of cell proliferation resembles that in the normal gastric mucosa. And it is believed that tumour cells remain confined to the mucosa as far as they grow to form the LS [13]. Our lineage analyses confirmed that POR with LS was derived from intramucosal SIG, whereas POR without LS and with a minor TC (< 30%), was derived from TC [14,15]. However, the TC was not always derived from early TUB but could also be derived from SIG, whereas LS was scarcely derived from TUB [15]. Therefore, as a morphological lineage marker, LS may take priority over TC. In addition, UGCs without LS or TC due to secondary loss of these markers are observed, which prompted us to adopt array CGH (aCGH) and unsupervised cluster analyses of the aCGH data to classify UGCs solely on the basis of similarity in the genomic copy number profile.

In differentiated-type gastric carcinomas (DGCs), our recent aCGH-based lineage analyses revealed two genetic lineages: one with copy-number loss of MYC and copy-number gain of TP53 (MYC− and TP53+), a dormant pattern, and the other with the copy-number gain of MYC and/or copy-number loss of TP53 (MYC+ and/or TP53−), an aggressive pattern. The dormant pattern accounted for 70% of intramucosal carcinoma samples and a half of the intramucosal part samples of invasive carcinomas. The invasive parts of invasive carcinomas mostly exhibited the aggressive copy number alteration (CNA) pattern. When the intramucosal part of an advanced cancer was dormant, the lineage was discontinuous between the mucosal and invasive parts. Therefore, the MYC−/TP53+ and MYC+ and/or TP53− CNA patterns may be signatures of dormant and aggressive TUBs, respectively [16].

In the present study, genomic DNA samples from the mucosal and invasive parts of early and advanced UGCs were prepared and subjected to gene copy-number analyses using aCGH, followed by unsupervised cluster analysis of the aCGH data. Based on these results, we examined relationship between morphological and genetic lineage markers and identified several useful lineage marker genes for UGC.

Methods

The Institutional Review Board on Medical Ethics at Shiga University of Medical Science approved this study on the condition that the UGC samples used were anonymous. Written informed consent was not required because this retrospective study used archival samples.

Tissue samples

This study included 29 surgically resected, buffered formalin-fixed, paraffin embedded UGCs: 20 with LS in at least part of the tumour (LS+, 9 intramucosaltumours and 11 invasive tumours) and 9 without LS but containing a small TC (LS−/TC+, all invasive tumours) (Table 1). TC was defined as a well or moderately differentiated adenocarcinoma component comprising ≤ 30% of the entire tumour [15]. All samples were selected from GC cases diagnosed in our department from 1997 to 2011. Intramucosal LS+ UGC patients averaged 57.6 years of age (range, 48–79) and patients with invasive LS+ UGCs 60.2 years (range, 48–79) and patients with invasive LS−/TC+ UGCs 62.2 years (range; 50–75). The macroscopic classification was determined according to the Japanese Classification of Gastric Cancer with TNM staging [3].

Table 1. Summary of clinicopathological characteristics of 29 UGCs

LS evaluation

LS was defined as in a previous study [17]. In brief, LS+ regions had small carcinoma cells confined to the stroma at the gland-neck level that gradually differentiated to signet ring cells in the superficial (and deep) lamina propria (Figure 1a). The absence of LS in intramucosal regions of the tumour was defined by four patterns: 1) contact of small carcinoma cells to the muscularis mucosae in SIG, 2) mucinous adenocarcinoma, 3) POR and 4) the presence of a TC (Figure 1b-f).

thumbnailFigure 1. Histological appearances of intramucosal parts of undifferentiated-type gastric carcinomas (UGCs). A signet ring cell carcinoma (SIG) component with a layered structure in case A107 (a). Small carcinoma cells are distributed in the deeper part just above or in the muscularis mucosae in a SIG component in case M109 (b). A mucinous adenocarcinoma component in case A107 (c). A poorly differentiated adenocarcinoma component in case SM106 (d). Minor tubular components in cases SM105 and SM201, respectively (e, f).

Laser microdissection and DNA preparation

Tumour tissue samples were obtained from 5-μm-thick tissue sections using a LMD6000 laser microdissection system (Leica Microsystems, Wetzlar, Germany). For invasive cancers, DNA samples were obtained from both the intramucosal and invasive parts. For each sample, cancer tissues were obtained from an area >6 mm2, in which cancer cells accounted for ≥70% of the total cell count. Tissue samples were digested in 200 μg/ml proteinase K solution for approximately 72 hours at 37.0°C and genomic DNA extracted with phenol/chloroform.

Whole genome amplification

Sample DNA was amplified using the GenomePlex Whole Genome Amplification Kit (WGA2 Kit; Sigma, St. Louis, USA) [18]. For some DNA samples that could not be sufficiently amplified, the WGA5 Kit (Sigma) was employed.

Array CGH

An oligo CGH microarray (60K, 60-mer) (Agilent, Santa Clara, USA) was used in this study, according to the manufacturer's instructions. In brief, the amplified tumour and control DNA samples were non-enzymatically labelled with Cy5 and Cy3, respectively, using the Genome DNA ULS Labelling Kit (Agilent) and competitively hybridized to the microarray. The hybridized array images were captured using a DNA microarray scanner (Agilent) and then the fluorescence intensity of the tumour and control at each probe dot was calculated by Feature Extraction Ver.9.5.3 (Agilent). The array data were normalized using Genomic Workbench software Ver.5.0 (Agilent). The positions of oligomers are based on the Human Genome February 2009 assembly (hg19). Copy-number gains and losses were defined as changes in the logarithm to the base 2 of the tumour to reference signal intensity ratio (T/R) greater than 0.3219 and less than −0.3219, respectively.

Cluster analysis

To perform novel subtyping of UGC samples based on genomic profile similarity in this study, an unsupervised hierarchical cluster analysis was applied across 63 samples from 29 UGC cases by using the Cluster 3.0 and TreeView software programs. The clustering algorithm was set to complete linkage clustering using an uncentered correlation. To enable unsupervised cluster analysis, we performed unbiased reduction in probe number from around 60,000 to several thousands of probes. For this purpose, we selected large genes because the greater number of corresponding probes resulted in improved signal-to-noise ratio of the representative gene copy numbers. The unsupervised strategy enabled us to set an internal standard to validate clustering results; the copy number profiles in samples of the same tumour should be more similar than any copy number profiles from another tumour because the gene alterations in the process of carcinogenesis are largely common among the samples from the same tumour.

Statistical analyses

Differences in contingency tables were assessed for statistical significance using Fisher's exact test. A P < 0.05 (2-sided) was considered statistically significant. The Welch’s t test was used to evaluate the difference in mean DNA copy number for each probe between two clusters of samples. The Bonferroni correction was used to correct for multiple comparisons.

Results

Samples analysed with array CGH

Tissue samples were excised from 29 archived GC specimens by laser microdissection. The tissue sample population included 11 regions (from 9 intramucosal SIGs), of which 9 regions were LS+ and the other two LS−, 26 regions (from 11 LS+ invasive UGC), of which 10 were LS+ mucosal regions, 8 were LS− mucosal regions and 8 were invasive regions, and 26 regions (from 9 LS−/TC+ invasive UGCs): 9 intramucosal POR, 9 intramucosal TC and 8 invasive regions.

Genome wide copy number alterations

A plot of the genetic aberration penetrance for all chromosomes is shown for LS+ UGCs and LS−/TC+ UGCs in Figure 2a and Figure 2b, respectively. Copy-number gains and losses were more common in LS−/TC+ UGCs than in LS+ UGCs. The most frequent copy-number gains were detected at 3q26 (7/63 samples), 5p15 (8/63), 8p23 (9/63), 8q24 (7/63) and 12p12 (6/63), while the most frequent copy-number losses were found at 7q36 (5/63) and 12p12 (5/63).

thumbnailFigure 2. Frequency of copy-number alterations at the chromosome level. The percentage of the samples that have CNAs for each chromosome in the LS+ UGCs (a) and LS/TC+ UGCs (b). Gains and losses are indicated with red and green, respectively.

Copy-number alterations (CNAs) common to all the samples from the same tumour were called stemline changes [14] and estimated to occur at the earliest stage of tumourigenesis and to be inherent into tumour lineage. Stemline gains of 3q26 were detected in 2/20 cases of invasive LS+ UGCs and none of invasive LS−/TC+ UGCs. In contrast, stemline gains of 5p15, 8p23 and 12p12 were detected in 2/9 cases of invasive LS−/TC+ UGCs but in no case of invasive LS+ UGCs. No stemline losses were detected in any cases of UGCs.

Previous studies using chromosomal or array CGH analyses [19-27] reported that frequent CNAs in gastric cancers (common to both UGC and DGC) were chromosomal gains at 3q, 5p, 7p, 8q, 13q, 17q, 20p and 20q, and losses at 4q, 5q, 6q, 9p, 17p, 18q and 21q. In the UGCs examined in the present study, all previously reported CNAs were observed except gains at 17q and 20p and losses at 5q and 6q. Gains at 8p and 12p were common in LS−/TC+ UGCs. Copy-number gains at 8q24 were common in both types of UGCs, with 4/20 cases of LS+ UGCs and 3/9 cases of invasive LS−/TC+ UGCs, but these were not stemline changes.

Impartial selection of genes reflecting the whole genome profile

To classify UGC samples based on the overall similarity in the profile of gene copy number changes, we used unsupervised hierarchical cluster analysis. For this purpose, it was necessary to reduce the number of gene probes used in the cluster analysis from 60K to several thousands. The reduced number of genes should still reflect the whole genome profile if impartially selected. To fulfil these conditions, we selected genes based solely on the size of genes (the numbers of corresponding probes). After repeated trials of cluster analyses using genes of various minimum sizes (or probe numbers per gene), we observed that most CNAs from the same tumour were clustered more closely together than any samples from another UGC case when we analysed only genes with 3 or more probes per gene: a total of 5019 genes.

Classification of UGC using hierarchical cluster analysis

We applied an unsupervised two-dimensional hierarchical clustering algorithm, to a total of the 63 DNA samples from 29 UGCs. The samples were classified into two major clusters A and B, based on similarity in the genome profile (Figure 3). Of 63 samples, 30 LS+ UGCs were classified into cluster A and only 7 into cluster B. For LS−/TC+ UGCs, 8 samples were classified into cluster A and 18 into cluster B. All Intramucosal LS+ UGCs were included in cluster A. Clusters A and B had significantly different proportions of morphological subtypes (P = 0.0001).

thumbnailFigure 3. Unsupervised hierarchical cluster analysis of array-based comparative genomic hybridization (aCGH) data. Gene copy-number gains and losses are indicated by red and green, respectively. A total of 63 samples from 29 UGCs were classified into two major clusters: A and B. Most samples of LS+ UGCs were included in cluster A and most LS−/TC+ UGCs samples were in cluster B. All the Intramucosal LS+ UGCs were included in cluster A.

Copy number alterations of MYC and TP53

Gains at 8q24 were common alterations for both LS+ and LS−/TC+ UGCs. The representative genes located at this locus is MYC. Gains of MYC were detected in 2/11 of Intramucosal LS+ UGCs (18.2%), 6/26 of LS+ UGC (23.1%) and 8/26 of invasive LS−/TC+ UGCs (30.1%). The aggressive pattern (MYC+ and/or TP53−) was detected in 6/11 of Intramucosal LS+ UGCs (54.5%), 11/26 of invasive LS+ UGCs (42.3%) and 20/26 of invasive LS−/TC+ UGCs (76.9%; Figure 4). Therefore, the aggressive pattern was more frequently detected in invasive LS−/TC+ UGCs than in LS+ UGCs (P = 0.0195). The dormant pattern (MYC− and TP53+) was not detected in any of the UGC samples, even those from intramucosal GCs (Figure 4).

thumbnailFigure 4. Array CGH data of MYC and TP53 in LS+ UGCs and LS−/TC+ UGCs. LS+ UGCs are divided into intramucosal cancers and invasive cancers. Numerals mean the base 2 logarithm of the test/reference signal intensity ratios of array CGH data. Significant gains and losses are indicated with red and green, respectively. The samples marked with and without grey margin are included in cluster B and cluster A, respectively, in Figure 3.

Copy number alterations of genes other than MYC or TP53

As mentioned above, 5p15 was one of the most frequent gain sites in invasive LS−/TC+ UGCs (8/26; 30.7%), but was not detected in any of the 37 intramucosal and invasive LS+ UGCs (Figure 2). The target genes located at this locus may include the telomerase reverse transcriptase gene (TERT) because a TERT gain was more frequently detected in invasive LS−/TC+ UGCs than intramucosal and invasive LS+ UGCs (16/26 vs. 1/37, P < 0.0001) (Figure 5). In contrast, losses of TERT were detected in 4/37 samples of intramucosal and invasive LS+ UGCs (10.8%) but not in invasive LS−/TC+ UGCs.

thumbnailFigure 5. Array CGH data of genes other than MYC and TP53 with significantly different T/R ratio between clusters A and B. UGCs are divided into clusters A and B that were determined in Figure 3. The heat map indicates the base 2 logarithm of the test/reference signal intensity ratios of array CGH data. Gains and losses are indicated with red and green, respectively.

Welch’s t test was performed to compare the mean T/R ratio between the samples in cluster A and those in cluster B at each 2756 probe loci of 979 cancer-related genes. Fourty-three gene probes, belonging to 40 genes, had significantly different mean T/R ratios between Clusters A and B at a level of P < 0.05 after Bonferroni correction (Table 2). Of the 40 genes, 6 genes (KIT, RAN, RAB39B, RAB9A, RAB37 and TERT), including proto-oncogenes, have been implicated in enhanced tumour growth, and 8 genes (ETS1, SPI1, ETV6, EPHA7, EPHA5, EPHB2, EPHA10 and TRIO) in invasion/metastasis and 3 genes (APC, NF1 and MEN1) in tumour suppression (Table 2). Most of log2 T/R ratios of the 43 distinguishing gene probes were of opposite sign between clusters A and B, with greater in absolute values in cluster B (Figure 5).

Table 2. List of 40 genes that have CNAs significantly different between clusters A and B

Discussion

Based on chromosomal CGH analysis, we have reported that there are two distinct UGC lineages: the LS+ lineage derived from early SIG and LS−/TC+ lineage dedifferentiated from TUB [15]. The former is characterized by LS and the latter by a small TC. However, there are also UGCs without these morphological lineage markers. In the present study, we classified UGC based on similarity in the whole genome copy number profile among samples using unsupervised hierarchical cluster analysis and examined the correlation between this gene-based classification and morphological lineage markers.

Using 5019 large genes and aCGH data from 63 DNA samples from 29 UGCs, we confirmed that most of the samples examined from the same tumour were clustered more closely together than in any other sample, thus fulfilling the criteria for our internal standard. On the basis of this observation, we performed an unsupervised two-dimensional hierarchical cluster analysis. All the samples were classified into two major clusters A and B (Figure 3). Cluster A was rich in LS+ UGCs, whereas cluster B was rich in LS−/TC+ UGCs. This difference was statistically significant (P = 0.0001) and indicates that the classification by the presence or absence of LS and TC is well correlated with the genomic-profile-based classification and validates the LS+ and LS−/TC+ as lineage markers.

All the intramucosal LS+ UGCs were included in cluster A, suggesting that most of UGCs in cluster A were derived from intramucosal SIG, and that the LS−/TC+ UGCs in cluster A may have secondarily lost LS. The LS−/TC+ UGCs in cluster A may also be derived from SIG, as suggested by chromosomal CGH studies [15]. In contrast, LS in advanced LS+ UGCs in cluster B (A107, A108 and A110) was virtually indistinguishable morphologically but showed genomic constitutions different from LS in cluster A. This may be a kind of phenocopy; a fraction of LS+ UGCs were considerably similar in genomic profile to LS−/TC+ UGCs. Although LS exhibits regular cell proliferation and differentiation and a superficially spreading dormant growth [13], it is suggested that LS itself is not a marker of persistent tumour dormancy but has the potential to progress to an advanced stage with the prognosis as poor as that for LS−/TC+ UGCs. This situation may resemble that in chronic myeloid leukaemia, in which blastic transformation occurs after a dormant phase of well retained cellular differentiation.

Most UGCs exhibited the aggressive genomic pattern (TP53− and/or MYC+), even 55% of intramucosal LS+ UGCs, an incidence comparable to that in invasive UGCs. The dormant pattern (MYC− and TP53+) was not detected in any of the UGC samples, even in intramucosal UGCs. These intramucosal UGCs are distinct from early DGCs, in which 70% are of the dormant type [16]. Therefore, TP53 and MYC are not as useful prognostic markers for UGCs.

To explore other genes important for differentiation of genetic lineage and for UGC prognosis, we first compared the profiles of chromosomal copy-number alterations (CNAs) between LS+ and LS−/TC+ UGCs. As shown in Figure 2, CNAs detected in LS+ tumours but not in LS−/TC+ tumours, include 3q26 gain, a locus likely to include SKIL because the average SKIL copy number was greater in LS+ tumours than in LS−/TC+ tumours (P = 0.0060). SKIL encodes SnoN protein that is proto-oncogenic by antagonizing cytostatic responses of TGF-β [28,29] and anti-oncogenic by activating p53 [30]. Those CNAs with the opposite pattern (present in LS−/TC+ tumours but not in LS+ tumours) were gains at 5p15, 8p23 and 12p12. The target genes at 5p15 and 12q12 include TERT, and KRAS, respectively because gains of TERT, and KRAS were more frequently detected in invasive LS−/TC+ UGCs than intramucosal and invasive LS+ UGCs (P < 0.0001 and P = 0.0032, respectively). No target gene was detected at 8q23.

Our second approach to identify lineage-specific CNAs was a screening of genes (from 979 cancer-related genes) that indicated significantly different mean T/R ratios between the samples of clusters A and B. We selected 40 genes that were significantly different between clusters (using t test after Bonferroni correction), of which 6 were related to enhanced tumour growth and 8 to invasion/metastasis (Table 2). As shown in Figure 5, genes that drive tumourigenesis were more common in cluster B and showed larger amplitude CNAs. Thus, UGCs in cluster B may be more dependent on oncogenic genomic alterations and less on environmental and epigenetic alterations than those in cluster A.

The possible drivers of tumour growth screened included KIT, TERT, and RAS family genes. KIT encodes a receptor tyrosine kinase that is activated by stem cell factor binding and initiates numerous signal transduction pathways linked with the process of apoptosis, proliferation and tumorigenesis [31]. RAS family genes encode small GTPase that plays a key role in transduction of signals from receptor kinase to the pathways of various cellular processes [32]. TERT encodes the telomerase catalytic subunit that plays not only an important role in cellular immortalization by telomere elongation [33,34] but also activates cell proliferation [35]. The possible drivers of invasion and metastasis screened include ETS1 and Ephrin receptor genes. ETS1 encodes a transcription factor, Ets1 proto-oncoprotein that promotes invasiveness and is an indicator of poor outcome in epithelial cancers through regulation of MMP1, MMP3, MMP9, uPA, VEGF and VEGF receptor expression [36]. Ephrin receptor genes, EPH39B, A7, A5 and A10 genes encode the ephrin receptor with tyrosine kinase activity that affects tumor growth, invasiveness, angiogenesis, and metastasis [37].

There were no significant differences in the mean copy number of CDH1 and its transcriptional repressor genes (SNAI1, SNAI2, ZEB1, ZEB2, TWIST1, etc.) between the clusters A and B, although these genes were reportedly associated with a poorly differentiated phenotype and poor clinical outcome [38]. However, these genes may still participate in UGC tumourigenesis through epigenetic silencing [39].

We are now extending this study to to validate UGC-associated genes as indicated by aCGH by quantitative PCR and to correlate their genomic copy number to gene expression and prognosis. Thereafter, using quantitative PCR analyses instead of aCGH, similar analyses should be applied to a greater number of tumour cases with known outcomes.

Conclusions

Unsupervised cluster analyses of aCGH data of multiple samples from early and advanced UGCs have demonstrated that early UGCs, including LS+ types in which polarity of cell proliferation and differentiation is well retained, have aggressive potential. Therefore, eradication of UGCs at early stages may thus contribute to better patient survival. In addition, it was observed that the two UGC lineages, one derived from early SIG and the other from TUB, have different genomic copy-number alteration profiles, resulting in different sets of genes contributing to tumourigenesis. The latter lineage from TUB may be more dependent on genomic copy-number alterations and have a poorer outcome than UGCs derived from SIG.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

AS performed most of the experiments, participated in the experimental design performed most of the data analyses and drafted the manuscript. KM participated in sample preparation for array CGH. TN participated partly in cluster analysis. VTND participated in varying selections of genes for cluster analysis, and TH, AA and YF provided critical comments and suggested revisions to the manuscript. HS conceived the study, designed the research route, and guided the experiments and the data analyses. All authors have read and approved the final manuscript.

Acknowledgements

We would like to express cordial thanks to Professor Yataro Daigo, Shiga University of Medical Science, for several suggestions on hierarchical cluster analysis. This work was supported in part by Grants-in-Aid for Scientific Research (C-21590369, C-25460454) from the Japan Society of the Promotion of Science.

References

  1. Lauren P: The two histological main types of gastric carcinoma: diffuse and so-called intestinal-type carcinoma. An attempt at a histoclinical classification.

    Acta Pathol Microbiol Scand 1965, 64:31-49. PubMed Abstract OpenURL

  2. Carneiro F: Classification of gastric carcinomas.

    Curr Diagn Pathol 1997, 4:51-59. Publisher Full Text OpenURL

  3. Japanese Gastric Cancer Association: Japanese classification of gastric carcinoma: 3rd English edition.

    Gastric Cancer 2011, 14:101-112. PubMed Abstract | Publisher Full Text OpenURL

  4. Endoh Y, Tamura G, Watanabe H, Ajioka Y, Motoyama T: The common 18-base pair deletion at codons 418-423 of the E-cadherin gene in differentiated-type adenocarcinomas and intramucosal pre cancerous lesions of the stomach with the features of gastric foveolar epithelium.

    J Pathol 1999, 189:201-206. PubMed Abstract | Publisher Full Text OpenURL

  5. Saito A, Shimoda T, Nakanishi Y, Ochiai A, Toda G: Histologic heterogeneity and mucin phenotypic expression in early gastric cancer.

    Pathol Int 2001, 51:165-171. PubMed Abstract | Publisher Full Text OpenURL

  6. Otsuji E, Yamaguchi T, Sawai K, Takahashi T: Characterization of Signet Ring Cell Carcinoma of the Stomach.

    J Surg Oncol 1998, 67:216-220. PubMed Abstract | Publisher Full Text OpenURL

  7. Sawada T, Nakata T, Takasugi N, Maeda K, Hanawa Y, Shimizu K, Hirayama M, Takeda T, Mori T, Koide R, Tsunoda A, Nagahara N, Yamamoto K: Mass screening for neuroblastoma in infants in Japan. Interim report of a mass screening study group.

    Lancet 1984, 2:271-273. PubMed Abstract OpenURL

  8. Sawada T, Sugimoto T, Tanaka T, Kawakatsu H, Ishii T, Matsumura T, Horii Y: Number and cure rate of neuroblastoma cases detected by the mass screening program in Japan: future aspects.

    Med Pediatr Oncol 1987, 15:14-17. PubMed Abstract | Publisher Full Text OpenURL

  9. Murphy SB, Cohn SL, Craft AW, Woods WG, Sawada T, Castleberry RP, Levy HL, Prorok PC, Hammond GD: Do children benefit from mass screening for neuroblastoma? Consensus Statement from the American Cancer Society Workshop on Neuroblastoma Screening.

    Lancet 1991, 337:344-346. PubMed Abstract | Publisher Full Text OpenURL

  10. Kaneko Y, Kanda N, Maseki N, Sakurai M, Tsuchida Y, Takeda T, Okabe I, Sakurai M: Different karyotypic patterns in early and advanced stage neuroblastomas.

    Cancer Res 1987, 47:311-318. PubMed Abstract | Publisher Full Text OpenURL

  11. Kaneko Y, Kobayashi H, Watanabe N, Tomioka N, Nakagawara A: Biology of neuroblastomas that were found by mass screening at 6 months of age in Japan.

    Pediatr Blood Cancer 2006, 46:285-291. PubMed Abstract | Publisher Full Text OpenURL

  12. Sugihara H, Hattori T, Imamura Y, Noriki S, Fukuda M, Katsura K, Tsuchihashi Y, Fujita S: Morphology and modes of cell proliferation in earliest signet ring cell carcinomas induced in canine stomachs by N-ethyl-N’-nitro-N-nitrosoguanidine.

    J Cancer Res Clin Oncol 1991, 117:197-204. PubMed Abstract | Publisher Full Text OpenURL

  13. Sugihara H, Hattori T, Fukuda M, Fujita S: Cell proliferation and differentiation in intramucosal and advanced signet ring cell carcinomas of the human stomach.

    Virchows Arch A 1987, 411:117-127. Publisher Full Text OpenURL

  14. Peng DF, Sugihara H, Mukaisho K, Tsubosa Y, Hattori T: Alterations of chromosomal copy number during progression of diffuse-type gastric carcinomas: metaphase- and array-based comparative genomic hybridization analyses of multiple samples from individual tumours.

    J Pathol 2003, 201:439-450. PubMed Abstract | Publisher Full Text OpenURL

  15. Peng DF, Sugihara H, Mukaisho K, Ling ZQ, Hattori T: Genetic lineage of poorly differentiated gastric carcinoma with a tubular component analysed by comparative genomic hybridization.

    J Pathol 2004, 203:884-895. PubMed Abstract | Publisher Full Text OpenURL

  16. Nakayama T, Ling ZQ, Mukaisho K, Hattori T, Sugihara H: Lineageeanalysis of early and advanced tubular adenocarcinomas of the stomach: continuous or discontinuous?

    BMC Cancer 2010, 10:311-321. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  17. Natsagdorj L, Sugihara H, Bamba M, Hattori T: Intratumoural heterogeneity of intestinal expression reflects environmental induction and progression-related loss of induction in undifferentiated-type gastric carcinomas.

    Histopathology 2008, 53:685-697. PubMed Abstract | Publisher Full Text OpenURL

  18. Little SE, Vuononvirta R, Reis-Filho JS, Natrajan R, Iravani M, Fenwick K, Mackay A, Ashworth A, Pritchard-Jones K, Jones C: Array CGH using whole genome amplification of fresh-frozen and formalin-fixed, paraffin-embedded tumor DNA.

    Genomics 2006, 87:298-306. PubMed Abstract | Publisher Full Text OpenURL

  19. Koizumi Y, Tanaka S, Mou R, Koganei H, Kokawa A, Kitamura R, Yamauchi H, Ookubo K, Saito T, Tominaga S, Matsumura K, Shimada H, Tsuchida N, Sekihara H: Changes in DNA copy number in primary gastric carcinomas by comparative genomic hybridization.

    Clin Cancer Res 1997, 3:1067-1076. PubMed Abstract | Publisher Full Text OpenURL

  20. Sakakura C, Mori T, Sakabe T, Ariyama Y, Shinomiya T, Date K, Hagiwara A, Yamaguchi T, Takahashi T, Nakamura Y, Abe T, Inazawa J: Gains, losses, and amplifications of genomic materials in primary gastric cancers analyzed by comparative genomic hybridization.

    Genes Chromosomes Cancer 1999, 24:299-305. PubMed Abstract | Publisher Full Text OpenURL

  21. van Grieken NC, Weiss MM, Meijer GA, Hermsen MA, Scholte GH, Lindeman J, Craanen ME, Bloemena E, Meuwissen SG, Baak JP, Kuipers EJ: Helicobacter pylori-related and -non-related gastric cancers do not differ with respect to chromosomal aberrations.

    J Pathol 2000, 192:301-306. PubMed Abstract | Publisher Full Text OpenURL

  22. Kimura Y, Noguchi T, Kawahara K, Kashima K, Daa T, Yokoyama S: Genetic alterations in 102 primary gastric cancers by comparative genomic hybridization: gain of 20q and loss of 18q are associated with tumor progression.

    Mod Pathol 2004, 17:1328-1337. PubMed Abstract | Publisher Full Text OpenURL

  23. Hidaka S, Yasutake T, Kondo M, Takeshita H, Yano H, Haseba M, Tsuji T, Sawai T, Nakagoe T, Tagawa Y: Frequent gains of 20q and losses of 18q are associated with lymph node metastasis in intestinal-type gastric cancer.

    Anticancer Res 2003, 23:3353-3357. PubMed Abstract OpenURL

  24. Takada H, Imoto I, Tsuda H, Sonoda I, Ichikura T, Mochizuki H, Okanoue T, Inazawa J: Screening of DNA copy-number aberrations in gastric cancer cell lines by array-based comparative genomic hybridization.

    Cancer Sci 2005, 96:100-110. PubMed Abstract | Publisher Full Text OpenURL

  25. Yang S, Jeung HC, Jeong HJ, Choi YH, Kim JE, Jung JJ, Rha SY, Yang WI, Chung HC: Identification of genes with correlated patterns of variations in DNA copy number and gene expression level in gastric cancer.

    Genomics 2007, 89:451-459. PubMed Abstract | Publisher Full Text OpenURL

  26. Fan B, Dachrut S, Coral H, Yuen ST, Chu KM, Law S, Zhang L, Ji J, Leung SY, Chen X: Integration of DNA copy number alterations and transcriptional expression analysis in human gastric cancer.

    PLoS One 2012, 7:e29824. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Rossi E, Klersy C, Manca R, Zuffardi O, Solcia E: Correlation between genomic alterations assessed by array comparative genomic hybridization, prognostically informative histologic subtype, stage, and patient survival in gastric cancer.

    Hum Pathol 2011, 42:1937-1945. PubMed Abstract | Publisher Full Text OpenURL

  28. Deheuninck J, Luo K: Ski and SnoN, potent negative regulators of TGF-beta signaling.

    Cell Res 2009, 19:47-57. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Jahchan NS, Luo K: SnoN in mammalian development, function and diseases.

    Curr Opin Pharmacol 2010, 10(6):670-675. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Jahchan NS, Ouyang G, Luo K: Expression Profiles of SnoN in Normal and Cancerous Human Tissues Support Its Tumor Suppressor Role in Human Cancer.

    PLoS One 2013, 8(2):e55794. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Heinrich MC, Rubin BP, Longley BJ, Fletcher JA: Biology and genetic aspects of gastrointestinal stromal tumors: KIT activation and cytogenetic alterations.

    Hum Pathol 2002, 33:484-495. PubMed Abstract | Publisher Full Text OpenURL

  32. Agarwal R, Jurisica I, Mills GB, Cheng KW: The Emerging Role of the RAB25 Small GTPase in Cancer.

    Traffic 2009, 10(11):1561-1568. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  33. Cohen SB, Graham ME, Lovrecz GO, Bache N, Robinson PJ, Reddel RR: Protein composition of catalytically active human telomerase from immortal cells.

    Science 2007, 315:1850-1853. PubMed Abstract | Publisher Full Text OpenURL

  34. Colgin LM, Reddel RR: Telomere maintenance mechanisms and cellular immortalization.

    Curr Opin Genet Dev 1999, 9:97-103. PubMed Abstract | Publisher Full Text OpenURL

  35. Sarin KY, Cheung P, Gilison D, Lee E, Tennen RI, Wang E, Artandi MK, Oro AE, Artandi SE: Conditional telomerase induction causes proliferation of hair follicle stem cells.

    Nature 2005, 436:1048-1052. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Jürgen D: The Biology of the Ets1 Proto-Oncogene.

    Mol Cancer 2003, 2:29. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  37. Pasquale EB: Eph receptors and ephrins in cancer: bidirectional signaling and Beyond.

    Nat Rev Cancer 2010, 10(3):165-180. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  38. Humar B, Guilford P: Hereditary diffuse gastric cancer: a manifestation of lost cell polarity.

    Cancer Sci 2009, 100:1151-1157. PubMed Abstract | Publisher Full Text OpenURL

  39. Dumont N, Wilson MB, Crawford YG, Reynolds PA, Sigaroudinia M, Tlsty TD: Sustained induction of epithelial to mesenchymal transition activates DNA methylation of genes silenced in basal-like breast cancers.

    Proc Natl Acad Sci USA 2008, 105:14867-14872. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

Pre-publication history

The pre-publication history for this paper can be accessed here:

http://www.biomedcentral.com/1755-8794/6/25/prepub