Abstract
Purpose and Experimental Design: To identify cancer-related genes, the expression profiles of colorectal cancer cells and normal epithelial cells were examined and compared using laser microdissection and cDNA microarray analysis. From these combined techniques, several cancer-related genes, including TROP2, were identified. TROP2 is known as a calcium signal transducer and is highly expressed in several types of tumors. However, no studies have investigated the significance of TROP2 expression in colorectal cancer. Thus, the expression status of TROP2 was investigated in 74 colorectal cancer samples by quantitative real-time reverse transcription-PCR and immunohistochemical studies.
Results: Laser microdissection and cDNA microarray analysis showed that there were 84 overexpressed genes in cancer cells. One of the highly overexpressed genes was TROP2. Quantitative real-time reverse transcription-PCR showed that TROP2 expression in cancer samples was significantly higher than in normal samples (P < 0.001). The samples were divided into high (n = 26) and low (n = 48) TROP2 expression groups. The cases with high TROP2 expression showed a higher frequency of liver metastasis (P = 0.005) and more cancer-related death (P = 0.046). Those cases also had an inclination of deeper depth of invasion (P = 0.064) and more lymph node metastasis (P = 0.125). Interestingly, the patients with high TROP2 expression tumors had poorer prognosis (P = 0.0036). Multivariate analysis showed that TROP2 expression status was an independent prognostic factor (relative risk, 2.38; 95% confidence interval, 1.29-4.74; P < 0.01).
Conclusion:TROP2 is one of the cancer-related genes that correlates with biological aggressiveness and poor prognosis of colorectal cancer. Thus, TROP2 is a possible candidate gene for diagnosis and molecular target therapy of colorectal cancer.
Colorectal cancer is one of the most prevalent cancers in the world. In Japan, the disease rate of colorectal cancer patients has doubled over the past 20 years, with ∼75 of 100,000 people suffering from the disease today. Additionally, colorectal cancer has been the second cause of death in neoplastic disease (1). Recently, molecular target therapy and cancer immunotherapy for solid cancers have been introduced to the clinic (2–5). However, indications for these therapies are limited due to the low frequency of target gene expression, unstable effectiveness, and/or severe side effects (5, 6). Thus, there is a pressing need to explore novel cancer-specific genes to serve as molecular targets for therapy and cancer specific immunotherapy.
In this study, the gene expression profiles were compared between cancer cells and normal epithelial cells using the combined techniques of laser microdissection and cDNA microarray analysis to explore cancer-related or cancer-specific genes. From this analysis, several cancer-related genes were identified. Among them, TROP2 showed markedly different expression between cancer cells and normal epithelial cells.
Fornaro et al. cloned the TROP2 gene, which encodes a 35,709-Da type 1 transmembrane protein with a single transmembrane domain and is homologous to TROP1/KSA/GA733-2. Moreover, TROP2, a cell surface receptor, has been shown to play a role in regulating the growth of carcinoma cells (7). TROP2 is also responsible for gelatinous drop-like corneal dystrophy (8). TROP2 is expressed at high levels in human trophoblast cells and it has been reported as one of the highly specific genes to the tumors, such as ovarian and bladder cancers (9–12). On the other hand, our institute developed tumor-specific immunotherapy using MAGE peptide-pulsed dendritic cells (13). We are planning to expand the therapy using other tumor-specific genes. Depending on the outcome of this analysis, TROP2 may become a candidate gene. In addition, the clinical significance of the gene has not been studied in colorectal cancer. Thus, we investigated TROP2 expression in clinical samples of colorectal cancer to determine its clinicopathologic significance.
Materials and Methods
Tissue sampling, laser microdissection, and cDNA microarray analysis. The samples of cancer tissues and noncancerous tissues were obtained from 16 patients with colorectal cancer who underwent surgical resection in Kyushu University Hospital (Beppu, Japan). Written informed consent was obtained from all patients. Tumors and adjacent normal tissues were immediately embedded in Tissue-Tek OCT compound medium (Sakura, Tokyo, Japan) and were kept frozen at −80°C until laser microdissection was done. Serial 8-μm frozen sections were generated by a cryostat. Sections were mounted onto a foil-coated glass slide, 90 FOIL-SL25 (Leica Microsystems, Wetzlar, Germany) for laser microdissection. Slides were stained with H&E at room temperature and dehydrated with ethanol. The Application Solutions Laser Microdissection System (Leica Microsystems) was introduced for laser microdissection to obtain the cancer cells and normal epithelial cells and to discard the mesenchymal tissues. Laser microdissection was done for several sequential sections and total RNA was extracted from each section (14). As the extracted total RNA was insufficient for hybridization to the cDNA microarray, the RNA was subjected to T7-based RNA amplification (15). The purity and concentration of the amplified RNA were determined by an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA) as described previously (14). In brief, high-quality amplified RNA run on a bioanalyzer typically has the shape of a hump peak and one marker peak, indicating no contamination of rRNA. Of 16 cancer and 16 normal microdissected and T7-based amplified RNA samples, 8 samples from the cancer sections and 10 samples from the normal sections were determined to be of sufficient quality. Each 8 RNA samples from cancer sections and the mixture of 10 RNA samples from normal epithelium were hybridized competitively to a cDNA microarray containing 12,814 genes. A list of genes on this cDNA microarray is available from http://www.agilent.com/chem/genelists.
Microarray analysis. After subtracting the local and global background signals, the expression values were calculated as the intensity of the dye-normalized red (Cy5) and green (Cy3) channel signals. The data flagged as poor quality according to the Agilent data extraction software were removed from the analysis. All data calculated by data extraction software were imported to the Rosetta Luminator system version 2.0 (Rosetta Biosoftware, Kirkland, WA). Candidate genes were selected that fulfilled the following criteria: the control (Cy3) intensity was <700, the fold changes were >2.5, and the P was <0.01 (14). Moreover, within the selected genes that met these criteria, genes that were up-regulated in three or more RNA samples (eight total samples) were further analyzed.
Semiquantitative real-time PCR. The semiquantitative real-time reverse transcription-PCR assay used 74 operatively resected paired cancer and normal samples that were not used for microarray analysis. Total RNA was extracted from each bulk sample and cDNA was synthesized from 8.0 μg total RNA as described previously (16). The purity and concentration of total RNA were determined using an Agilent 2100 Bioanalyzer. The following primers were used to amplify the TROP2 gene: sense primer 5′-GCCTACTACTTCGAGAGGGACA-3′ and antisense primer 5′-CAGTTCCTTGATCTCCACCTTC-3′. The glyceraldehyde-3-phosphate dehydrogenase (sense primer 5′-TTGGTATCGTGGAAGGACTCA-3′ and antisense primer 5′-TGTCATCATATTTGGCAGGTTT-3′) gene was used as an internal control. The reaction was done in a LightCycler system (Roche Applied Science, Indianapolis, IN) using the LightCycler FastStart DNA Master SYBR Green I kit (Roche Diagnostics, Mannheim, Germany). Details of each reaction are described elsewhere (17). Briefly, thermal cycling for all genes was initiated with a denaturation step of 95°C for 10 minutes and then consisted of 40 cycles at 95°C for 10 seconds, 65°C (60°C for glyceraldehyde-3-phosphate dehydrogenase) for 10 seconds, and 72°C for each optimal length (1 second/25 bp). All calculated concentrations of target genes were divided by the amount of the endogenous reference (glyceraldehyde-3-phosphate dehydrogenase) to obtain the normalized TROP2 expression values. Each assay was done in triplicate.
Immunohistochemistry. Of 74 cases, immunohistochemical staining was done in 34 selected cases. Immunohistochemical studies of TROP2 were done on formalin-fixed, paraffin-embedded surgical sections. After deparaffinization and blocking, the antigen-antibody reaction was incubated overnight at 4°C. The LSAB2 kit (DAKO, Kyoto, Japan) was applied to detect the signal of the TROP2 antigen-antibody reaction. All sections were counterstained with hematoxylin. The purified goat polyclonal antibody against the purified recombinant human TROP2 extracellular domain (R&D Systems, Inc., Minneapolis, MN) was used at 5 μg/mL.
Statistical analysis. Quantitative real-time reverse transcription-PCR data were calculated with JMP 5 for Windows software (SAS Institute, Inc., Cary, NC). Differences between groups were estimated using the Student's t test and the χ2 test. The survival curves were estimated by the Kaplan-Meier method and the comparison between the curves was made by the log-rank test. The relative risk (RR) was calculated using the Cox proportional hazards model. A probability level of 0.05 was chosen for statistical significance.
Results
Microarray analysis. In the microarray analysis, 84 genes were identified that had a higher expression level in cancer cells than in normal epithelial cells. Among these genes, some are involved with signal transduction (27.4%), transcription (16.7%), transport (13.1%), metabolism (4.8%), cell adhesion (4.8%), secretion (3.6%), apoptosis (3.6%), and an unknown function (16.7%; Table 1). Of these genes, some were already reported as associated with colorectal cancer aggressiveness, such as MMP7 (18, 19), MMP11 (20, 21), CD44 (22–24), CD61 (25, 26), CRIPTO (27, 28), ENC1 (29), SLC7A5 (30), DPEP1 (31), and activin A (32, 33).
Gene name . | Symbols . | Systematic code . | Locus . | Fold overexpression . |
---|---|---|---|---|
Dipeptidase 1 (renal) | DPEP1 | J05257 | 16q24,3 | 17.59 |
Teratocarcinoma-derived growth factor 1 | CRIPTO | AU124747 | 3p21,31 | 8.35 |
Human megakaryocyte-enhanced gene transcript 1 protein | MEGT1 | AF195764 | 6p21,33 | 8.22 |
Cell surface glycoprotein Trop-2 | TROP2 | BG259957 | 1p32,1 | 7.84 |
Integrin β3 (platelet glycoprotein IIIa, antigen CD61) | CD61 | NM_000212 | 17q21,32 | 7.12 |
Matrix metalloproteinase-7 (matrilysin, uterine) | MMP7 | NM_002423 | 11q21-q22 | 6.71 |
Phosphoprotein C8FW | TRIB3 | BG387820 | 20p13 | 6.56 |
Human cell surface glycoprotein CD44 (CD44) gene | CD44 | L05411 | 11p13 | 6.08 |
Epiregulin | EREG | BG235918 | 4q13,3 | 5.71 |
Matrix metalloproteinase-11 (stromelysin 3) | MMP11 | AW301093 | 22q11,23 | 5.70 |
Human clone PP1195 unknown mRNA | — | AF217970 | 8p23,1 | 4.80 |
Human cDNA FLJ20093fis, clone COL04263 | ANKRD10 | AK000100 | 13q34 | 4.73 |
Cadherin 3, type 1, P-cadherin (placental) | CDH3 | NM_001793 | 1q22,1 | 4.63 |
Fatty acid–binding protein 6, ileal (gastrotropin) | FABP6 | X90908 | 5q23,35 | 4.48 |
KIAA0546 protein | — | AB011118 | 12q15 | 4.36 |
Ataxia telangiectasia group D–associated protein | TRIM29 | AF230388 | 11q22-q23 | 4.24 |
TAR (HIV) RNA-binding protein 1 | TARBP1 | U38847 | 1q42,3 | 4.06 |
Expressed sequence tags | PPAT | AA442070 | 4q12 | 3.82 |
Expressed sequence tags | AK126318 | BE672109 | 17q23,2 | 3.68 |
Human mRNA; cDNA DKFZp564C053 | — | AL049246 | 3q23 | 3.60 |
Solute carrier family 12, member 2 | SLC12A2 | U30246 | 5q23,2 | 3.58 |
Human Na,K-ATPase α-1 subunit gene | ATP1A1 | M30309 | 1p13,1 | 3.56 |
Ectodermal-neural cortex (with BTB-like domain) | ENC1 | BC000418 | 5q12-13,3 | 3.47 |
Bone morphogenetic protein 7 (osteogenic protein 1) | BMP7 | BE395650 | 20q13,31 | 3.42 |
Solute carrier family 7, member 5 | SLC7A5 | M80244 | 16q24,2 | 3.42 |
Human DNA for apoER2 | LRP8 | D86407 | 1p32,3 | 3.41 |
Human mRNA for KIAA0761 protein | MCLC | AB018304 | 1p13,3 | 3.40 |
Human activin β-A subunit | INHBA | X57579 | 7p14,1 | 3.39 |
Deleted in lymphocytic leukemia, 2 | DLEU2 | AW978447 | 13q14,2 | 3.38 |
Human mRNA for nuclear pore complex protein | NUP107 | AJ295745 | 12q15 | 3.33 |
Growth factor receptor-bound protein 7 | GRB7 | AU148656 | 17q12 | 3.32 |
Human mRNA for KIAA0619 protein | ROCK2 | AB014519 | 2p25,1 | 3.32 |
Meiotic recombination (Saccharomyces cerevisiae) 11 homologue A | MRE11A | AF073362 | 11q21 | 3.29 |
Chromosome 20 open reading frame 119 | TOMM34 | AK026760 | 20q13,12 | 3.26 |
ZFM1 protein alternatively spliced product | SF1 | D26121 | 11q13 | 3.24 |
Mevalonate kinase (mevalonic aciduria) | MVLK | BG474232 | 12q24,11 | 3.22 |
Oviductal glycoprotein 1, 120 kDa (mucin 9, oviductin) | OVGP1 | NM_002557 | 1p13,2 | 3.16 |
Zinc finger protein 195 | BRDT | AW025438 | 1p22,1 | 3.14 |
Centromere protein F (350/400 kDa, mitosin) | CENPF | U19769 | 1q41 | 3.12 |
Cadherin 6, type 2, K-cadherin (fetal kidney) | CDH6 | AU149929 | 5p15,1-p14 | 3.09 |
Solute carrier family 11, member 2 | SLC11A2 | AI888673 | 12q13,12 | 3.08 |
Zinc finger, X-linked, duplicated A | ZXDA | AL031115 | Xp11,1 | 3.07 |
Runt-related transcription factor 1 | RUNX1 | D43969 | 21q22,12 | 3.06 |
Expressed sequence tags | TDRD9 | AA844124 | 14q32 | 3.03 |
Spectrin, α, erythrocytic 1 (elliptocytosis 2) | SPTA1 | AA703344 | 1q23,1 | 3.02 |
Human transcription factor SL1 (by similarity) | TAF1A | AK001054 | 1q41 | 3.01 |
Glycoprotein glucosyltransferase precursor (by similarity) | UGCGL2 | AK001735 | 13q32,1 | 2.98 |
Formyl peptide receptor-like 1 | FPRL1 | BG541691 | 19q13,41 | 2.93 |
Guanine nucleotide binding protein 4 | GNG4 | AW593228 | 1q42,3 | 2.93 |
Prostaglandin E receptor 4 (subtype EP4) | PTGER4 | NM_000958 | 5p13,1 | 2.93 |
Human cDNA FLJ10517fis, clone NT2RP2000812 | ASPM | AK001379 | 1q31,3 | 2.92 |
Human G-protein-coupled receptor gene | GPR19 | U55312 | 12p132 | 2.88 |
Human mRNA; cDNA DKFZp566P1124 | — | AL110236 | 11q14,1 | 2.87 |
Human thiazide-sensitive NaCl cotransporter | SLC12A3 | U44128 | 16q13 | 2.87 |
KIAA0410 gene product | NUPL1 | BI599177 | 13q12,13 | 2.87 |
High-mobility group protein isoform I-C | HMGA2 | BG250825 | 12q14,3 | 2.83 |
Human, paired box gene 9 | PAX9 | BC001159 | 14q13,3 | 2.82 |
Guanidinoacetate N-methyltransferase | GAMT | NM_000156 | 19p13,3 | 2.81 |
Karyopherin α3 | KPNA3 | D89618 | 13q14,2 | 2.80 |
KIAA0008 gene product | DLG7 | BI087140 | 14q22,3 | 2.79 |
Nucleophosmin/nucleoplasmin 3 | NPM3 | AI631542 | 10q24,31 | 2.77 |
Zinc finger protein 200 | ZNF200 | NM_003454 | 16p13,3 | 2.75 |
Ras-related associated with diabetes | RRAD | AI186786 | 16q22,1 | 2.73 |
Human metabotropic glutamate receptor 8 | GRM8 | U92459 | 7q31,33 | 2.70 |
Zinc finger protein homologous to Zfp91 in mouse | ZFP91 | AB057443 | 11q12,1 | 2.70 |
Human G-protein-coupled receptor GPR86 | GPR86 | AF295368 | 3q25,1 | 2.67 |
Nuclear receptor subfamily 1, group I, member 3 | NR1I3 | Z30425 | 1q23,3 | 2.66 |
Human homeobox gene | HOXB | AF287967 | 17q21,32 | 2.62 |
Human interleukin-17 | IL17 | U32659 | 6p12 | 2.62 |
Zinc finger protein 184 (Kruppel-like) | ZNF184 | BG254958 | 6p22,1 | 2.61 |
Macaque somatostatin I | SST | M19318 | 3q27,3 | 2.60 |
Expressed sequence tags | PTPN22? | AA401425 | 1p13,2 | 2.58 |
Human ZNF43 | ZNF43 | X59244 | 19p12 | 2.57 |
Homer, neuronal immediate early gene, 1B | HOMER1 | BI858644 | 5q14,1 | 2.55 |
Homo sapiens clone 23903 mRNA sequence | — | AF035281 | 7q36,1 | 2.55 |
Natural killer cell group 7 sequence | NKG7 | S69115 | 19q13,41 | 2.55 |
Expressed sequence tags | — | AI025099 | 19p13,2 | 2.54 |
Human mRNA for CSR1 | SCARA3 | AB007829 | 8p21,1 | 2.53 |
Neuronal pentraxin II | NPTX | BC009924 | 7q22,1 | 2.52 |
KIAA0322 protein | NEDL1 | AB002320 | 7p14,1-p13 | 2.51 |
Cyclin E1 | CCNE1 | BG761079 | 19q12 | 2.50 |
Death-associated transcription factor 1 | DATF1 | AB002331 | 20q13,33 | 2.50 |
Human DNA for single-minded gene 2 | — | D85922 | 21q22,13 | 2.50 |
RasGAP-related protein | IQGAP2 | AAB37765 | 5q13,3 | 2.50 |
Gene name . | Symbols . | Systematic code . | Locus . | Fold overexpression . |
---|---|---|---|---|
Dipeptidase 1 (renal) | DPEP1 | J05257 | 16q24,3 | 17.59 |
Teratocarcinoma-derived growth factor 1 | CRIPTO | AU124747 | 3p21,31 | 8.35 |
Human megakaryocyte-enhanced gene transcript 1 protein | MEGT1 | AF195764 | 6p21,33 | 8.22 |
Cell surface glycoprotein Trop-2 | TROP2 | BG259957 | 1p32,1 | 7.84 |
Integrin β3 (platelet glycoprotein IIIa, antigen CD61) | CD61 | NM_000212 | 17q21,32 | 7.12 |
Matrix metalloproteinase-7 (matrilysin, uterine) | MMP7 | NM_002423 | 11q21-q22 | 6.71 |
Phosphoprotein C8FW | TRIB3 | BG387820 | 20p13 | 6.56 |
Human cell surface glycoprotein CD44 (CD44) gene | CD44 | L05411 | 11p13 | 6.08 |
Epiregulin | EREG | BG235918 | 4q13,3 | 5.71 |
Matrix metalloproteinase-11 (stromelysin 3) | MMP11 | AW301093 | 22q11,23 | 5.70 |
Human clone PP1195 unknown mRNA | — | AF217970 | 8p23,1 | 4.80 |
Human cDNA FLJ20093fis, clone COL04263 | ANKRD10 | AK000100 | 13q34 | 4.73 |
Cadherin 3, type 1, P-cadherin (placental) | CDH3 | NM_001793 | 1q22,1 | 4.63 |
Fatty acid–binding protein 6, ileal (gastrotropin) | FABP6 | X90908 | 5q23,35 | 4.48 |
KIAA0546 protein | — | AB011118 | 12q15 | 4.36 |
Ataxia telangiectasia group D–associated protein | TRIM29 | AF230388 | 11q22-q23 | 4.24 |
TAR (HIV) RNA-binding protein 1 | TARBP1 | U38847 | 1q42,3 | 4.06 |
Expressed sequence tags | PPAT | AA442070 | 4q12 | 3.82 |
Expressed sequence tags | AK126318 | BE672109 | 17q23,2 | 3.68 |
Human mRNA; cDNA DKFZp564C053 | — | AL049246 | 3q23 | 3.60 |
Solute carrier family 12, member 2 | SLC12A2 | U30246 | 5q23,2 | 3.58 |
Human Na,K-ATPase α-1 subunit gene | ATP1A1 | M30309 | 1p13,1 | 3.56 |
Ectodermal-neural cortex (with BTB-like domain) | ENC1 | BC000418 | 5q12-13,3 | 3.47 |
Bone morphogenetic protein 7 (osteogenic protein 1) | BMP7 | BE395650 | 20q13,31 | 3.42 |
Solute carrier family 7, member 5 | SLC7A5 | M80244 | 16q24,2 | 3.42 |
Human DNA for apoER2 | LRP8 | D86407 | 1p32,3 | 3.41 |
Human mRNA for KIAA0761 protein | MCLC | AB018304 | 1p13,3 | 3.40 |
Human activin β-A subunit | INHBA | X57579 | 7p14,1 | 3.39 |
Deleted in lymphocytic leukemia, 2 | DLEU2 | AW978447 | 13q14,2 | 3.38 |
Human mRNA for nuclear pore complex protein | NUP107 | AJ295745 | 12q15 | 3.33 |
Growth factor receptor-bound protein 7 | GRB7 | AU148656 | 17q12 | 3.32 |
Human mRNA for KIAA0619 protein | ROCK2 | AB014519 | 2p25,1 | 3.32 |
Meiotic recombination (Saccharomyces cerevisiae) 11 homologue A | MRE11A | AF073362 | 11q21 | 3.29 |
Chromosome 20 open reading frame 119 | TOMM34 | AK026760 | 20q13,12 | 3.26 |
ZFM1 protein alternatively spliced product | SF1 | D26121 | 11q13 | 3.24 |
Mevalonate kinase (mevalonic aciduria) | MVLK | BG474232 | 12q24,11 | 3.22 |
Oviductal glycoprotein 1, 120 kDa (mucin 9, oviductin) | OVGP1 | NM_002557 | 1p13,2 | 3.16 |
Zinc finger protein 195 | BRDT | AW025438 | 1p22,1 | 3.14 |
Centromere protein F (350/400 kDa, mitosin) | CENPF | U19769 | 1q41 | 3.12 |
Cadherin 6, type 2, K-cadherin (fetal kidney) | CDH6 | AU149929 | 5p15,1-p14 | 3.09 |
Solute carrier family 11, member 2 | SLC11A2 | AI888673 | 12q13,12 | 3.08 |
Zinc finger, X-linked, duplicated A | ZXDA | AL031115 | Xp11,1 | 3.07 |
Runt-related transcription factor 1 | RUNX1 | D43969 | 21q22,12 | 3.06 |
Expressed sequence tags | TDRD9 | AA844124 | 14q32 | 3.03 |
Spectrin, α, erythrocytic 1 (elliptocytosis 2) | SPTA1 | AA703344 | 1q23,1 | 3.02 |
Human transcription factor SL1 (by similarity) | TAF1A | AK001054 | 1q41 | 3.01 |
Glycoprotein glucosyltransferase precursor (by similarity) | UGCGL2 | AK001735 | 13q32,1 | 2.98 |
Formyl peptide receptor-like 1 | FPRL1 | BG541691 | 19q13,41 | 2.93 |
Guanine nucleotide binding protein 4 | GNG4 | AW593228 | 1q42,3 | 2.93 |
Prostaglandin E receptor 4 (subtype EP4) | PTGER4 | NM_000958 | 5p13,1 | 2.93 |
Human cDNA FLJ10517fis, clone NT2RP2000812 | ASPM | AK001379 | 1q31,3 | 2.92 |
Human G-protein-coupled receptor gene | GPR19 | U55312 | 12p132 | 2.88 |
Human mRNA; cDNA DKFZp566P1124 | — | AL110236 | 11q14,1 | 2.87 |
Human thiazide-sensitive NaCl cotransporter | SLC12A3 | U44128 | 16q13 | 2.87 |
KIAA0410 gene product | NUPL1 | BI599177 | 13q12,13 | 2.87 |
High-mobility group protein isoform I-C | HMGA2 | BG250825 | 12q14,3 | 2.83 |
Human, paired box gene 9 | PAX9 | BC001159 | 14q13,3 | 2.82 |
Guanidinoacetate N-methyltransferase | GAMT | NM_000156 | 19p13,3 | 2.81 |
Karyopherin α3 | KPNA3 | D89618 | 13q14,2 | 2.80 |
KIAA0008 gene product | DLG7 | BI087140 | 14q22,3 | 2.79 |
Nucleophosmin/nucleoplasmin 3 | NPM3 | AI631542 | 10q24,31 | 2.77 |
Zinc finger protein 200 | ZNF200 | NM_003454 | 16p13,3 | 2.75 |
Ras-related associated with diabetes | RRAD | AI186786 | 16q22,1 | 2.73 |
Human metabotropic glutamate receptor 8 | GRM8 | U92459 | 7q31,33 | 2.70 |
Zinc finger protein homologous to Zfp91 in mouse | ZFP91 | AB057443 | 11q12,1 | 2.70 |
Human G-protein-coupled receptor GPR86 | GPR86 | AF295368 | 3q25,1 | 2.67 |
Nuclear receptor subfamily 1, group I, member 3 | NR1I3 | Z30425 | 1q23,3 | 2.66 |
Human homeobox gene | HOXB | AF287967 | 17q21,32 | 2.62 |
Human interleukin-17 | IL17 | U32659 | 6p12 | 2.62 |
Zinc finger protein 184 (Kruppel-like) | ZNF184 | BG254958 | 6p22,1 | 2.61 |
Macaque somatostatin I | SST | M19318 | 3q27,3 | 2.60 |
Expressed sequence tags | PTPN22? | AA401425 | 1p13,2 | 2.58 |
Human ZNF43 | ZNF43 | X59244 | 19p12 | 2.57 |
Homer, neuronal immediate early gene, 1B | HOMER1 | BI858644 | 5q14,1 | 2.55 |
Homo sapiens clone 23903 mRNA sequence | — | AF035281 | 7q36,1 | 2.55 |
Natural killer cell group 7 sequence | NKG7 | S69115 | 19q13,41 | 2.55 |
Expressed sequence tags | — | AI025099 | 19p13,2 | 2.54 |
Human mRNA for CSR1 | SCARA3 | AB007829 | 8p21,1 | 2.53 |
Neuronal pentraxin II | NPTX | BC009924 | 7q22,1 | 2.52 |
KIAA0322 protein | NEDL1 | AB002320 | 7p14,1-p13 | 2.51 |
Cyclin E1 | CCNE1 | BG761079 | 19q12 | 2.50 |
Death-associated transcription factor 1 | DATF1 | AB002331 | 20q13,33 | 2.50 |
Human DNA for single-minded gene 2 | — | D85922 | 21q22,13 | 2.50 |
RasGAP-related protein | IQGAP2 | AAB37765 | 5q13,3 | 2.50 |
TROP2 expression in colorectal cancer samples and corresponding normal tissues. Quantitative real-time reverse transcription-PCR was done on 74 paired samples to show TROP2 mRNA expression in clinical samples. Quantitative real-time reverse transcription-PCR showed that TROP2 expression in cancer samples was significantly higher (averaged expression values of cancer were 8.34-fold higher; P < 0.0001) than those in normal samples (Fig. 1). Figure 2 shows the results of immunohistochemical studies of TROP2 expression in representative clinical samples of well-differentiated adenocarcinoma (Fig. 2A), moderately differentiated adenocarcinoma (Fig. 2B), poorly differentiated adenocarcinoma (Fig. 2C), and mucinous adenocarcinoma (Fig. 2D). The majority of the TROP2 expression was observed in the cancer cells, the minority in stromal cells, and none in the normal colonic epithelium. Immunohistochemical studies revealed that the staining was strong (n = 7), moderate (n = 9), or weak (n = 18) in the tumor cells, whereas very weak or none in the normal cells in all 34 cases. There was a significant difference in immuohistochemical staining between the tumor and the normal samples (P < 0.01), and the data are similar to those obtained from mRNA expression analysis. All 16 tumors with strong or moderate immunohistochemical expression showed higher mRNA expression values (>0.2). The expression of TROP2 mRNA relatively associated with protein expression. Among the undifferentiated cell types, several regions of the cancer had strong immunohistochemical staining for TROP2 (Fig. 2C and D).
High TROP2 expression correlates with clinicopathologic variables. The experimental samples were divided into two groups [the high expression group with TROP2 expression values (>0.2; n = 26) and the remaining samples in the low expression group (n = 48)] to investigate TROP2 expression in association with clinicopathologic variables (Table 2). The border of the two groups was defined by an upper limit, including 95% of the expression values of the normal samples. The incidence in liver metastasis was significantly higher (P = 0.005) in the high expression group (9 of 26, 34.6%) than in the low expression group (4 of 48, 8.3%), and the incidence of cancer death was significantly higher (P = 0.046) in the high expression group (9 of 26, 34.6%) than in the low expression group (7 of 48, 14.6%). The high expression group also had inclinations of deeper invasion (P = 0.064) and more lymph node metastasis (P = 0.125) than the low expression group.
Variables . | Expression . | . | P . | |||
---|---|---|---|---|---|---|
. | High (n = 26) . | Low (n = 48) . | . | |||
Age | 66.6 ± 3.8 | 67.5 ± 2.8 | 0.701 | |||
Sex | ||||||
Male | 14 | 30 | 0.470 | |||
Female | 12 | 18 | ||||
Histologic cell type | ||||||
Well | 6 | 18 | 0.420 | |||
Moderately | 19 | 29 | ||||
Poorly and others | 1 | 1 | ||||
Tumor site* | ||||||
Right colon | 9 | 13 | 0.498 | |||
Left colon | 17 | 35 | ||||
Serosal invasion | ||||||
Absent | 14 | 36 | 0.064 | |||
Present | 12 | 12 | ||||
Lymph node metastasis | ||||||
Absent | 12 | 31 | 0.125 | |||
Present | 14 | 17 | ||||
Lymphatic permeation | ||||||
Absent | 13 | 29 | 0.388 | |||
Present | 13 | 19 | ||||
Venous permeation | ||||||
Absent | 18 | 39 | 0.241 | |||
Present | 8 | 9 | ||||
Liver metastasis | ||||||
Absent | 17 | 44 | 0.005† | |||
Present | 9 | 4 | ||||
Cancer-related death | ||||||
Alive | 17 | 41 | 0.046† | |||
Death | 9 | 7 |
Variables . | Expression . | . | P . | |||
---|---|---|---|---|---|---|
. | High (n = 26) . | Low (n = 48) . | . | |||
Age | 66.6 ± 3.8 | 67.5 ± 2.8 | 0.701 | |||
Sex | ||||||
Male | 14 | 30 | 0.470 | |||
Female | 12 | 18 | ||||
Histologic cell type | ||||||
Well | 6 | 18 | 0.420 | |||
Moderately | 19 | 29 | ||||
Poorly and others | 1 | 1 | ||||
Tumor site* | ||||||
Right colon | 9 | 13 | 0.498 | |||
Left colon | 17 | 35 | ||||
Serosal invasion | ||||||
Absent | 14 | 36 | 0.064 | |||
Present | 12 | 12 | ||||
Lymph node metastasis | ||||||
Absent | 12 | 31 | 0.125 | |||
Present | 14 | 17 | ||||
Lymphatic permeation | ||||||
Absent | 13 | 29 | 0.388 | |||
Present | 13 | 19 | ||||
Venous permeation | ||||||
Absent | 18 | 39 | 0.241 | |||
Present | 8 | 9 | ||||
Liver metastasis | ||||||
Absent | 17 | 44 | 0.005† | |||
Present | 9 | 4 | ||||
Cancer-related death | ||||||
Alive | 17 | 41 | 0.046† | |||
Death | 9 | 7 |
Relative to splenic flexure.
Cox proportional model.
Survival analysis. The 5-year actuarial overall survival rates in patients with high TROP2 mRNA expression levels and those with low levels were 48.5% and 83.3%, respectively. Moreover, the high expression group also showed poorer prognosis in the Kaplan-Meier survival curve (P = 0.0036; Fig. 3). Table 3 revealed the result of RR and the 95% confidence interval (95% CI) by the Cox proportional hazards model. Table 3A shows the univariate analysis for all clinicopathologic variables. Regarding these variables, lymph node metastasis (RR, 2.711; 95% CI, 1.526-5.681; P = 0.0004), liver metastasis (RR, 2.409; 95% CI, 1.381-4.059; P = 0.003), depth of invasion (RR, 2.127; 95% CI, 1.276-3.675; P = 0.0041), histologic grade (RR, 2.097; 95% CI, 1.098-5.313; P = 0.0225), and high TROP2 expression (RR, 2.026; 95% CI, 1.216-3.501; P = 0.0071) were statistically significant. Table 3B shows the result of multivariate analysis in the final model, which included lymph node metastasis, depth of invasion, histologic grade, tumor site, and lymphatic permeation. In this model, the variable of high TROP2 expression was an independent prognostic predictor for the patients with colorectal cancer (RR, 2.38; 95% CI, 1.29-4.74; P = 0.005; Table 3B). Of the variables that were entered in the multivariate analysis, the variable of liver metastasis was excluded because there was a significant correlation with the variable of high TROP2 expression.
Variable . | RR (95% CI) . | P* . | ||
---|---|---|---|---|
(A) Univariate analysis | ||||
Lymph node metastasis (present) | 2.71 (1.53-5.68) | 0.000 | ||
Liver metastasis (present) | 2.41 (1.38-4.06) | 0.003 | ||
Serosal invasion (present) | 2.13 (1.28-3.68) | 0.004 | ||
Histologic cell type† | 2.10 (1.10-5.31) | 0.023 | ||
High TROP2 | 2.03 (1.22-3.50) | 0.007 | ||
Tumor site† (right colon) | 1.55 (0.92-2.59) | 0.098 | ||
Lymphatic permeation (present) | 1.54 (0.93-2.67) | 0.095 | ||
Gender (female) | 1.30 (0.78-2.20) | 0.310 | ||
Age at surgery (>65 y) | 1.27 (0.74-2.41) | 0.402 | ||
Venous permeation (present) | 1.13 (0.59-1.93) | 0.686 | ||
(B) Multivariate analysis | ||||
High TROP2 | 2.38 (1.29-4.74) | 0.005 | ||
Lymph node metastasis (present) | 2.04 (1.08-4.48) | 0.027 | ||
Histologic cell type† | 2.00 (0.97-5.29) | 0.060 | ||
Tumor site (right colon)‡ | 1.99 (1.06-3.89) | 0.031 | ||
Serosal invasion (present) | 1.70 (0.98-3.08) | 0.059 | ||
Lymphatic permeation (present) | 0.92 (0.50-1.70) | 0.788 |
Variable . | RR (95% CI) . | P* . | ||
---|---|---|---|---|
(A) Univariate analysis | ||||
Lymph node metastasis (present) | 2.71 (1.53-5.68) | 0.000 | ||
Liver metastasis (present) | 2.41 (1.38-4.06) | 0.003 | ||
Serosal invasion (present) | 2.13 (1.28-3.68) | 0.004 | ||
Histologic cell type† | 2.10 (1.10-5.31) | 0.023 | ||
High TROP2 | 2.03 (1.22-3.50) | 0.007 | ||
Tumor site† (right colon) | 1.55 (0.92-2.59) | 0.098 | ||
Lymphatic permeation (present) | 1.54 (0.93-2.67) | 0.095 | ||
Gender (female) | 1.30 (0.78-2.20) | 0.310 | ||
Age at surgery (>65 y) | 1.27 (0.74-2.41) | 0.402 | ||
Venous permeation (present) | 1.13 (0.59-1.93) | 0.686 | ||
(B) Multivariate analysis | ||||
High TROP2 | 2.38 (1.29-4.74) | 0.005 | ||
Lymph node metastasis (present) | 2.04 (1.08-4.48) | 0.027 | ||
Histologic cell type† | 2.00 (0.97-5.29) | 0.060 | ||
Tumor site (right colon)‡ | 1.99 (1.06-3.89) | 0.031 | ||
Serosal invasion (present) | 1.70 (0.98-3.08) | 0.059 | ||
Lymphatic permeation (present) | 0.92 (0.50-1.70) | 0.788 |
Cox proportional model.
Indicates moderately and poorly differentiated adenocarcinoma.
Relative to splenic flexure.
Discussion
For further understanding of cancer biology, it is important to identify cancer-related genes. To find such genes, a combination of techniques, laser microdissection and cDNA microarray analysis, was applied to clinical samples. Laser microdissection can isolate cancer cells from tumor tissue that consists of mixed populations of carcinoma cells and stromal cells, such as fibroblasts, macrophages, and lymphocytes (16). RNA extracted from these tissues was subjected to rigorous procedures for quality to ensure its utility in cDNA microarray analysis. Almost 50% to 60% of the frozen surgical samples passed the RNA quality check according to our experience. Through using both laser microdissection and microarray analysis, 84 cancer-related genes were identified that were overexpressed in cancer cells compared with normal cells in the colon.
This study showed that TROP2 was one of the most highly expressed genes in colorectal cancer cells compared with normal colon tissue. TROP2 is reported to be highly expressed in several types of cancers. For example, Nakajima et al. revealed that TROP2 is overexpressed in esophageal cancer, and the titer of serum TROP2 antibody correlated with tumor size (34). With respect to colon cancer, Kanai et al. studied the expression of mRNA of genes altered by 5-azacytidine treatment in cancer cell lines and found that TROP2 was overexpressed in clinical colorectal cancers (35). However, they did not investigate the correlation between the TROP2 expression status and the clinicopathologic factors. Our study of a comparison of both groups showed a significant difference in liver metastasis, and this was associated with cancer-related death. Although there were no significant differences in the other pathologic factors examined, there was a tendency in serosal invasion between the groups. Therefore, these results suggest that high TROP2 expression has an association with biological aggressiveness of colorectal cancer.
The reason why the cancers with high TROP2 expression show aggressive behavior is unclear. There are some studies suggesting the biological mechanism of high TROP2 expression associated with the aggressiveness of cancer. Fornaro et al. suggested that TROP2 not only acts as a calcium signal transducer but also has a receptor function activity. TROP2 has homology to serum insulin-like growth factor-II–binding proteins and has a conservative structure for potential cytoplasmic phosphorylation sites (7). Basu et al. revealed that the phosphorylation occurred on Ser303 by protein kinase C (36).
We are investigating the predicted function of TROP2 gene and taking notice of its intracellular domain of phosphatidylinositol 4,5-bisphosphate–binding sequence. Phosphatidylinositol 4,5-bisphosphate is the most important lipid in the cytoplasmic leaflet of the plasma membrane and is responsible for a wide range of membrane-related phenomena that consist of the functions of endocytosis, exocytosis, cytoskeletal attachment, enzyme activation, actin-binding protein, production of three second messengers, ion-channel activation, and binding site for PH and other domains (37). EI Sewedy et al. investigated that the phosphatidylinositol 4,5-bisphosphate–binding site of TROP2 has an ability of being phosphorylated by protein kinase C (38); thus, the signal transduction through TROP2 has a possibility of causing the aggressiveness of cancer cells. Consequently, we assume that the TROP2 gene has a receptor activity and influences the aggressive behavior of colorectal cancer cells during tumor progression. To investigate its biological function, we are planning the following experiments: (a) Establishment of the transfectant of the TROP2 gene to clarify the dominant proliferation and progression of invasiveness. (b) Interference of the TROP2 gene expression by small interfering RNA (or short hairpin RNA) to examine the alteration of its proliferation and invasiveness. (c) Verification of the binding activity between TROP2 and phosphatidylinositol 4,5-bisphosphate and exploration of the alteration of the binding activity by protein kinase C. (d) Utilization of the transfectant with mutant TROP2 gene, which lacks the phosphatidylinositol 4,5-bisphosphate–binding site to compare the result of in vitro analysis with that of experiment 1. Thus, high TROP2 expression would contribute to biological aggressiveness of cancer cells through paracrine and autocrine signaling pathway. Nonetheless, more precise understanding of the mechanism that is responsible for the aggressiveness is necessary.
Because TROP2 is overexpressed in colon cancer cells and not expressed in normal cells, it is a novel target for treatment. In fact, Mangino G et al. showed that TROP2 is one of the target molecules recognized by human CTLs (39). We investigated the effectiveness of cancer-specific immunotherapy using MAGE antigen as a dendritic cell therapy (13, 40–42). However, the patients that can be treated with this therapy are limited because the patients must meet conditions, such as adequate HLA typing and positive MAGE gene expression in the tumor. The percentage of such patients is ∼10% to 40% in several types of cancers. If TROP2 can become an immunotherapeutic molecular target, the combined specific immunotherapy with MAGE and TROP2 genes would be a powerful new treatment modality for colorectal cancer patients.
In conclusion, our results indicated that TROP2 mRNA and protein were overexpressed in colorectal cancer cells and high TROP2 expression status correlated with liver metastasis and poor prognosis. These findings suggest that TROP2 is one of the cancer-related genes that are associated with cancer aggressiveness and is one of the predictors of survival in patients with colorectal cancer. The overexpression of this gene has a possible role for assessment of the postoperative adjuvant therapy of colorectal cancer patients. Furthermore, this gene will be useful not only for diagnostics but also molecular targeted therapy.
Grant support: Core Research for Evolutional Science and Technology, Japan Science and Technology Agency; Uehara Memorial Foundation; Japan Society for the Promotion of Science Grant-in-Aid for Scientific Research (grants 17015032, 17109013, 17591411, 17591413, 17015032, and 16390381).
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.
Acknowledgments
We thank Dr. S. Tanaka for helpful discussion and T. Shimooka, K. Ogata, M. Oda, M. Kasagi, and Y. Nakagawa for excellent technical assistance.