| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Molecular Biology, Pathobiology, and Genetics |
Divisions of 1 Basic Science, 2 Medical Science, and 3 Population Science, Fox Chase Cancer Center, and 4 Division of Genetic and Preventive Medicine, Thomas Jefferson University, Philadelphia, Pennsylvania; 5 University of California at San Francisco, San Francisco, California; and 6 Division of Cancer Prevention, National Cancer Institute, Bethesda, Maryland
Requests for reprints: Anthony T. Yeung, Room R404, Fox Chase Cancer Center, 333 Cottman Avenue, Philadelphia, PA 19111-2497. Phone: 215-728-2488; Fax: 215-728-3647; E-mail: AT_Yeung{at}fccc.edu.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
Alterations in the gene expression pattern associated with single-hit mutations of the APC tumor suppressor gene are likely to represent the earliest molecular changes during colon carcinogenesis. The idea that there might be "one-hit" effects of APC mutation in heterozygous cells was supported years ago by Kopelovich and colleagues (4, 5) who observed profound genetic alterations in morphologically normal skin fibroblasts that were derived from FAP patients, including increased sensitivity to transformation by the KiMSV and SV40 viruses, altered cytoskeletal actin distribution, and increased expression of p53. Danes and colleagues (6) reported increased tetraploidy in morphologically normal colonic epithelial cells derived from FAP patients of the Gardner syndrome variant. Both of these reports indicated one-hit effects. These observations have recently been extended using RNA expression profiles in morphologically normal renal epithelial cells that were derived from patients affected with tuberous sclerosis complex and von Hippel-Lindau syndrome (7). Some or all of these early changes may have a direct bearing on subsequent tumor induction because similar expression profiles have been identified in the corresponding cancer cells (4–7). Thus, even a small growth advantage could increase the number of one-hit cells available for conversion to two-hit tumor cells that represent the paradigm for adenoma formation providing some selective advantage during colorectal cancer development.
Proteomics provides a direct approach to study thousands of proteins in a given tissue at the same time, facilitating the determination of critical pathways during cancer development. We have chosen to study the colonic crypt proteome by two-dimensional gel electrophoresis because this approach also provides global quantitative protein isoform information (8). We showed that, compared with the proteome profile of morphologically normal crypts sampled at
10 to 20 cm from the tumor of sporadic colorectal cancer patients (9), the morphologically normal FAP colonic crypts exhibit an altered proteome in a gene-specific manner, consistent with detectable single-hit effect associated with heterozygosity for a mutant APC tumor suppressor gene. These protein aberrations seem to be directly related to cancer and may represent some of the earliest known biomarkers during colorectal cancer development, potentially leading to the identification of molecular targets and agents that might inhibit or delay adenoma-carcinoma transition.
| Materials and Methods |
|---|
|
|
|---|
|
Data analysis. Two-dimensional gel Images (16 bits) were recorded by ProXPRESS Proteomic Imaging System and analyzed using the Progenesis Discovery Workstation Software v2004 (Nonlinear Dynamics, Inc.). For each colonic crypt sample, we obtained at least three two-dimensional gels per pH range, of which two gels were used for further image analysis. Total of 12 images in control group and 12 images in FAP group were used for each pH range image analysis. Five FAP patients each provided enough samples for all three pH ranges. However, the sixth patient for each pH range came from one of two additional FAP patients, giving a total of seven FAP patients in this report.
The Progenesis output files for spot intensities of gel replicates were averaged for corresponding protein spots so that each patient was represented by one file. We rejected any spot that was missing in >4 of the 24 gels in each pH range unless all but two of the spots absent were found in the same class of patients. The data were transformed into log base 2 space, median centered, and normalized for each patient, and then median centered and normalized for each protein. The results were subjected to analysis by Hierarchical Clustering (Cluster 3.0; visualized with Java TreeView software; ref. 11) and Genesis 1.6.0 software (Institute for Genomics and Bioinformatics, Graz University of Technology; ref. 12). The results are provided at our web site7 as Genesis project files accessible through the Genesis software that can be downloaded.8 Principal Component Analysis (PCA; refs. 13, 14) was performed using Genesis 1.6.0 software.
Protein identification. Robotic cutting of protein spots, trypsin digestion, automated MALDI-TOF mass spectrometry, and protein identification from Swiss-Prot database were previously described in detail (9, 10). The identity and function of the genes separating FAP from control were explored at the gene ontology databases FatiGO+ (15) and Database for Annotation, Visualization and Integrated Discovery.9
Western blots. One hundred micrograms of protein of each sample were resolved on a one-dimensional or a two-dimensional gel. Blotting was performed at 200 mA for 2 h in 5% Nonfat Dry milk/TBS. Primary antibodies shown were Anti-ECHS1 (#11305-1-AP rabbit polyclonal antibody; 1:1,000 dilution; Proteintech Group, Inc.) and anti-SOD2 (#ab13533 rabbit polyclonal antibody; 1:5,000 dilution; Abcam, Inc.). Secondary antibody was Goat Anti rabbit IgG (#62-6120; Zymed, Inc.), 1:2,000 dilution. Imaging was obtained on Kodak BioMax film, which was developed with Bio-Rad ECL reagent. The films were scanned by the ProXPRESS scanner and quantified by Progenesis image analysis software.
| Results |
|---|
|
|
|---|
808 nonredundant protein entries and
900 of their isoforms were identified with <2% false discovery rate (FDR) as determined by a database search to which was appended a decoy copy of the same database wherein the amino acid sequences of each protein entry had been randomized (9). Size and isoelectric point (pI) information on each spot allowed the detection and classification of posttranslational modifications.
|
Statistical Analyses
Although some interindividual variability in the colonic proteome of patients was apparent, the group differences among FAP crypts and the nonaffected crypts of sporadic colon cancer patients were highly significant. Initially, "unsupervised hierarchical clustering" was performed with each gel as an independent sample. Invariably, gel replicates clustered together correctly (data not shown). The conclusions in Fig. 2
are supported by the comparison of 7 FAP and 7 control patients for each of which we used the average of multiple gel replicates for statistical analysis. The samples from FAP cases, with no exception, were properly separated from controls for each of the three pH ranges (Fig. 2A–C), indicating that the phenotypic differences of FAP versus control crypts are the major factors distinguishing these two classes of patients.
|
|
For the purpose of demonstrating that proteins over the complete range of pIs can separate FAP from non-FAP patients, we selected 5 FAP samples and 5 control samples for which two-dimensional gels of all 3 pH ranges were available and mapped the protein spots of a P value of <0.01 (214 proteins with a P value of <0.01 and fold change of >0.5.xls: Supplementary Data S9) back onto the PCA analysis to test their agreement. The raw data and analyses are provided in Supplementary Data S10 and S11. With the provision that some of the protein spots in the overlapping regions of two pH ranges may have been counted more than once in this analysis, of a total of 2,556 protein spots included in the analysis, 335 proteins (13%) provided a P value of <0.01 as being differentially expressed in FAP versus control. These 335 spots were mapped into the PCA analysis in Fig. 3 as small squares in red color. The clear separation of FAP versus control and the lack of outlier red spots indicate that the two statistical methods are congruent.
FDR
The foregoing lists of P < 0.01 significant proteins were again used to perform PCA (Supplementary Data S12) and hierarchical cluster analysis (Supplementary Data S13) to illustrate the successful grouping of important proteins that distinguish these phenotypes. The PCA analysis in Supplementary Data S12 is the result of coupling two statistical evaluations, t test of individual proteins, and PCA of all the significant proteins. The spot distribution in PCA confirms that the t test results were correct. This congruency would not be automatically true had the t test P < 0.01 list been short relative to the list that can be produced from random chance. For example, when the randomized data file used in Fig. 2D was used in the t test, the number of proteins showing a P value of <0.01 was only 7 (data not shown), which equals a FDR of 3.4%. Thus, the low FDR of our data set was due to the high percentage of significant protein expression shifts detectable in our proteomic analysis.
The list of proteins that were up- or down-regulated as a result of the heterozygous APC mutation obtained by Progenesis image analysis and the fold changes in expression in combination with the P values of the t test performed by Genesis are shown in Supplementary Data S8 Table S2. Congruent conclusion of 214 proteomic changes, including 148 nonredundant protein entries and 66 protein isoforms by these two software analyses plus PCA, provides confidence that the assignments were statistically significant.
Validation of Two-Dimensional Gel Fluorescence Staining by Two-Dimensional Gel Western Blots
It is likely that the differential protein expression seen for 172 of the 214 biomarker proteins that distinguish FAP crypt from control crypt with a P value of <0.01 involved not the total protein of each biomarker but one or more isoforms of each protein, a situation that can arise from shifting in posttranslational modifications. We therefore chose to further validate several biomarker proteins by two-dimensional gel Western blotting experiments. Examples shown are for ECHS1 (Fig. 4
) and SOD2 (Fig. 5
). Although two-dimensional gel Western blotting experiments are less quantitative than two-dimensional gels stained with the fluorescent dye Sypro Ruby, they do validate that the protein identification for each isoform was correct and proportional to the changes observed with fluorescence staining presented in Supplementary Data S8 Table S2.
|
|
The difference between the mean ages of the FAP group versus the control group may be a confounding factor. A test of the effects of age was performed by exchanging two samples of the same gender and similar age: CP206 (female, age 41 years) from the control group with SID516 (female, age 39 years) from the FAP group. The resultant table of 2,556 protein spots for pH 4 to 11 used in the above FAP versus control analysis was used to calculate the P values separating group 1 from group 2 using the two-tailed t test of each protein. With the exchange of these two patients, only 36 identified proteins separated the two 5 groups; of these 36 proteins, only 9 proteins were in common with the 214 identified proteins that separated FAP from control (Supplementary Data S9).
Characterization of FAP Versus Controls
Apoptosis-related proteins. One of the proteomic changes with the greatest magnitude was the major isoform of ETHE1 protein (spot 2470; pH 5–8), elevated
27-fold in FAP crypts from nongene carriers (Supplementary Data S8 Table S2). We also observed a 4.3-fold increase of Baculoviral IAP repeat–containing protein 8 [Inhibitor of apoptosis-like protein 2: spot 367 (pH 5–8), 1.3-fold decrease and 1.8-fold increase for isoform spot 2399 (pH 4–7) and spot 1847 (pH 5–8), respectively] in FAP crypts, a 5.2-fold decrease in the Programmed cell death 6–interacting protein (isoform spot 633, pH 5–8), a 0.5-fold decrease in the 14-3-3 protein
(Stratifin; spot 2110, pH 4–7), a 1.7-fold decrease of nucleophosmin (nuclear phosphoprotein B23; spot 1823, pH 5–8), and a 4.5-fold decrease of an isoform of Annexin A11 (Calcyclin-associated annexin 50; spot 1161, pH 5–8; Supplementary Data S8 Table S2).
Proteins relevant to oxidative stress. Another significant FAP-specific phenotype is several-fold elevation in the isoforms of 18 mitochondrial enzymes in FAP crypts, including many enzymes of the tricarboxylic acid cycle, concomitant with a decrease in a number of glycolytic enzymes (Supplementary Data S8 Table S2).
Four additional proteins that are presumably associated with oxidative stress were observed in this study to be elevated in FAP crypts, as follows:
Changes in apoptosis, cell proliferation, and migration may require significant alteration of the cytoskeleton. For example, several cytokeratins known to be expressed in colonic mucosa are down-regulated in FAP crypts compared with control crypts (Supplementary Data S8 Table S2).
Other significant protein changes, notably FK506-binding protein 4, and Carbonic Anhydrase 1 (up to 10-fold increase in FAP), among others are shown in Supplementary Data S8 Table S2.
| Discussion |
|---|
|
|
|---|
20% of hereditary nonpolyposis colorectal cancer (HNPCC) and they are highly prevalent in sporadic forms of colorectal cancers, suggesting the involvement of similar and divergent pathways in the etiology of colorectal cancer (23). For example, although fibroblasts from FAP patients show significant deformities in cytoskeletal actin (24), fibroblasts from HNPCC patients seemed normal for cytoskeletal actin (4, 25), pointing to a gene-specific effect in single-hit cells by APC, which is not seen in mismatch repair–deficient cells. This result further supports the value of biomarkers, whether by phenotype, RNA expression arrays, or proteomics, in elucidating mechanisms of cancer progression using morphologically normal tissue from patients at varying genetic risk for cancer. Accordingly, a significant observation in this study is that the morphologically normal FAP and the non-FAP control colon crypt proteomes are substantially different, providing biomarkers that can potentially be used to reveal the pathogenesis and pathophysiology underlying the FAP condition. Importantly, this report has shown that the colonic crypt cells from the morphologically normal mucosa in adult FAP patients are one hit. With regard to the observed changes in the isoforms of cytokeratins, there has been substantial interest in using TPS (determining fragments of cytokeratin 18) and CYFRA 21.1 (determining fragments of cytokeratin 19) as serum markers of cancer, although recent studies seem inconclusive (26, 27). Based on the current studies of the cytokeratins and other examples, vide infra, it would seem that it is important to consider both quantitative and posttranslational modification changes to accurately evaluate the FAP phenotype. For example, neuroblastoma patients form autoantibodies against specific truncated isoforms I and III of β-tubulin that do not react with the β-tubulin isoforms I and III from noncancer patients (28). One of the two alternate spliced isoform of Pantothenate kinase, a rate-determining enzyme in CoA biosynthesis, confers feedback inhibition on the enzyme, leading to regulation of the enzyme activity and the level of acetyl-coA in the cell (29).
Two phenotypes of FAP crypts suggested by this study are the attenuation of apoptosis and the elevation of oxidative stress. ETHE1 (HSCO) increase was recently shown to down-regulate p53-dependent apoptosis (30), and its elevation may lead to attenuation of apoptosis during the aging of the crypt cells. ETHE1 specifically associates with histone deacetylase 1 (HDAC1) independently of Mdm2 and facilitates deacetylation of p53 at Lys373/382 by HDAC1 (31). By increasing ubiquitination and degradation of p53, ETHE1 reduces p53 protein levels. Consistent with this hypothesis of attenuation of apoptosis in FAP crypts were the observed changes in Baculoviral IAP repeat–containing protein 8, the Programmed cell death 6–interacting protein, the 14-3-3
protein, and Annexin A11.
The oxidative stress damage caused by a high level of oxidative phosphorylation in FAP crypts may be attenuated by SOD2, catalase, protein disulfide isomerase, and DJ-1. SOD2, elevated in FAP crypts, catalyzes the dismutation of O2-superoxide into oxygen and hydrogen peroxide. Elevated SOD2 expression in colon cancer cell lines is also correlated with enhanced survival after cytotoxic exposures (32), consistent with the antiapoptosis effect vide ante from elevated ETHE1. Catalase catalyzes the decomposition of hydrogen peroxide to water and oxygen. Y231 and Y386 residues in human catalase can be phosphorylated by c-Abl and Arg nonreceptor tyrosine kinases to regulate the catalase activity (33). Oxidative stress causes oxidation of disulfide bonds in proteins, which must be restored by protein disulfide isomerases (34). Another effect of oxidative damage is the oxidation of oncogene DJ-1 at cys106, which results in pI change to more acidic range, consistent with the observed DJ-1 isoform spot 2319 increase in FAP samples. Only this oxidized form of DJ-1 seems to have significant antiaggregation properties against
-synuclein (35) and protection against oxidative damage in neurons (36). This potential role for DJ-1 in the FAP crypt may be distinct from another role for DJ-1 as a regulator of the androgen receptor (37).
In the current report, we found the FAP crypt proteome to be significantly different from control crypts isolated from morphologically normal mucosa of patients with nonhereditary cancer. The proteomic alterations observed in the FAP syndrome encompassed many proteins in diverse pathways that are also seen in the corresponding tumors themselves, suggesting that, in general, the phenotypic effects of a heterozygous tumor suppressor gene mutation on the path to cancer involves similar proteome-wide changes; i.e., there is a one-hit phenotype effect of germline mutations in this tumor suppressor gene. Collectively, these alterations in the gene expression repertoire associated with single-hit mutations of tumor suppressor genes are likely to represent the earliest molecular changes during tumorigenesis. Furthermore, there may well be a further clarification of optimal targets at the two-hit stage of tumorigenesis in that those revealed in one-hit lesions would not include confounding secondary tumor effects. Thus, >60% of the most significant protein differences in the morphologically normal FAP crypts compared with control crypts are consistent with known or proposed physiology of colorectal cancer development. Moreover, because 172 of 214 significant FAP-related proteome changes involved proteins that had at least one more isoform identified in these gels, the FAP phenotype may involve critical change in one or more posttranslational modification systems that produced these pleotrophic proteome changes. In this regard, a notable, albeit a meaningful "exception," is the presumed lack of a "Warburg effect" and the dominant role of mitochondrial proteins in FAP crypts. The elevation of mitochondrial marker proteins and oxidative stress response proteins suggests that oxidative stress response may be an important early event that is associated with a one-hit mutation of the APC gene during the very early stages of carcinogenesis. Also, p53 in the morphologically normal fibroblasts from FAP patients, although abnormally stable, was not mutated in these cells (5), again pointing to subtle but significant occurrences very early in cancer development. Thus, there are many opportunities for the discovery of biomarkers for cancer intervention and prevention by studying the morphologically normal crypts of patients with varying risk for colon cancer. These one-hit cells might, therefore, be an important experimental system for chemoprevention studies.
| Disclosure of Potential Conflicts of Interest |
|---|
|
|
|---|
| Acknowledgments |
|---|
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.
We thank Drs. David Weinberg, Lisa Henske, Avram Hershko, and Joseph Testa for critical readings of this article; the physicians who enrolled the majority of FAP and control cases: Dr. Jonathan Terdiman, University of California at San Francisco, Dr. Daniel Chung, Massachusetts General Hospital, and Dr. James Church, Cleveland Clinic; and Dr. Margie Clapper of Fox Chase Cancer Center for her assistance in the recruitment of FAP research participants.
| Footnotes |
|---|
9 http://david.abcc.ncifcrf.gov/ ![]()
Received 3/14/08. Revised 6/ 3/08. Accepted 6/ 7/08.
| References |
|---|
|
|
|---|
B from the nucleus and inhibits p53-dependent apoptosis. Cancer Cell 2002;2:335–46.[CrossRef][Medline]
-synuclein. J Mol Biol 2006;356:1036–48.[CrossRef][Medline]
to the receptor. J Biol Chem 2001;276:37556–63.This article has been cited by other articles:
![]() |
E. L. Huttlin, X. Chen, G. A. Barrett-Wilt, A. D. Hegeman, R. B. Halberg, A. C. Harms, M. A. Newton, W. F. Dove, and M. R. Sussman Discovery and validation of colonic tumor-associated proteins via metabolic labeling and stable isotopic dilution PNAS, October 6, 2009; 106(40): 17235 - 17240. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |