We studied patients with Familial Adenomatous Polyposis (FAP) because they are virtually certain to develop colon cancer, and because much is known about the causative APC gene. We hypothesized that the inherited heterozygous mutation itself leads to changes in the proteome of morphologically normal crypts and the proteins that changed may represent targets for preventive and therapeutic agents. We determined the differential protein expression of morphologically normal colon crypts of FAP patients versus those of individuals without the mutation, using two-dimensional gel electrophoresis, mass spectrometry, and validation by two-dimensional gel Western blotting. Approximately 13% of 1,695 identified proteins were abnormally expressed in the morphologically normal crypts of APC mutation carriers, indicating that a colon crypt cell under the one-hit state is already abnormal. Many of the expression changes affect pathways consistent with the function of the APC protein, including apoptosis, cell adhesion, cell motility, cytoskeletal organization and biogenesis, mitosis, transcription, and oxidative stress response. Thus, heterozygosity for a mutant APC tumor suppressor gene alters the proteome of normal-appearing crypt cells in a gene-specific manner, consistent with a detectable one-hit event. These changes may represent the earliest biomarkers of colorectal cancer development, potentially leading to the identification of molecular targets for cancer prevention. [Cancer Res 2008;68(18):7579–86]

Colorectal cancer represents the second leading cause of cancer-related death in the United States (1). Its development is due to the accumulation of mutations of a number of tumor suppressor genes and oncogenes, as well as DNA repair genes controlling genomic stability (2). Patients with Familial Adenomatous Polyposis (FAP), a hereditary disease characterized by the presence of multiple polyps in the colon, carry germline mutations in one allele of the adenomatous polyposis coli (APC) tumor suppressor gene. Colon cancer in FAP patients requires at least one additional somatic genetic event at the other APC allele in such individuals, consistent with the “two-hit” hypothesis for cancer arising from defective tumor suppressor genes that are recessive at the cellular level (3). The APC gene is also mutated in most cases of sporadic colorectal cancer.

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 (47). 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.

Sample acquisition. Subjects were recruited for this research following approval by the respective Institutional Review Boards of Fox Chase Cancer Center (protocol number 00–852) and Thomas Jefferson University (protocol number 00–0023). Morphologically normal mucosa was obtained from eligible men and women who were previously diagnosed with FAP and whose corresponding APC gene mutations had been documented by DNA sequencing of circulating white cells (Table 1), including 4 male and 3 female FAP patients with a mean age of 30 y (range, 17–48 y). No correlation was observed between the location of the APC mutation and severity of the FAP phenotype. In view of the logistical difficulty in obtaining colon specimens from unaffected normal controls, morphologically normal mucosa was obtained from patients undergoing surgery for colorectal cancer who, however, had no family or personal history of colon cancer, or any other cancer, and were therefore classified as patients with sporadic colorectal cancer. The mucosal samples were obtained 10 to 20 cm away from the primary lesion in the descending colon and referred to, in this report, as control crypts (9). The morphologically normal colonic mucosa from sporadic colorectal cancer patients was used as a surrogate for the normal colonic mucosa of cancer-free individuals. The control patients included 5 male and 2 female patients ages 41 to 58 y, with a mean age of 51.1 y (Table 1). Because a FAP patient often undergoes colectomy at an early age, it was not entirely possible to age match the FAP cases with sporadic cancer patients (see below).

Table 1.

Patient samples used in the comparison of the colon crypt proteomes of FAP patients and controls (distal mucosa from sporadic cancer cases)

Patient codePathology diagnosisGenderAgeAPC Mutation
CP129 Sporadic, rectal cancer 47 Not determined 
CP203 Sporadic, rectal cancer 49 Not determined 
CP70 Sporadic, rectal cancer 54 Not determined 
CP144 Sporadic, rectal cancer 58 Not determined 
CP206 Sporadic, rectal cancer 41 Not determined 
CP73 Sporadic, sigmoid cancer 56 Not determined 
SID514 FAP, no cancer 17 Codon 1148 APC3443delCT 
SID516 FAP, no cancer 39 Codon 178 (genetically attenuated) 
SID601 FAP 48 IVS4+1>A 
SID602 FAP 42 APC3714delT 
SID608 FAP 24 APC3927del5 
SID618 FAP, no cancer 35 Codon 1148 APC3443delCT 
SID622 FAP, no cancer 26 APC3183del5 
Patient codePathology diagnosisGenderAgeAPC Mutation
CP129 Sporadic, rectal cancer 47 Not determined 
CP203 Sporadic, rectal cancer 49 Not determined 
CP70 Sporadic, rectal cancer 54 Not determined 
CP144 Sporadic, rectal cancer 58 Not determined 
CP206 Sporadic, rectal cancer 41 Not determined 
CP73 Sporadic, sigmoid cancer 56 Not determined 
SID514 FAP, no cancer 17 Codon 1148 APC3443delCT 
SID516 FAP, no cancer 39 Codon 178 (genetically attenuated) 
SID601 FAP 48 IVS4+1>A 
SID602 FAP 42 APC3714delT 
SID608 FAP 24 APC3927del5 
SID618 FAP, no cancer 35 Codon 1148 APC3443delCT 
SID622 FAP, no cancer 26 APC3183del5 

Two-dimensional gel electrophoresis. We adhered strictly to our established Standard Operating Procedures for two-dimensional gel electrophoresis and their analyses (9, 10). Patient samples were randomized during the gel runs to avoid systematic artifacts. The FAP and control samples were analyzed concurrently. Our procedure for colonic crypts preparation from morphologically normal colon was previously described in detail (9, 10). The first-dimension isoelectric focusing (IEF) was performed using in-gel rehydration method for linear pH 4 to 7 and pH 5 to 8 immobilized pH gradient (IPG) strips (ReadyStrip; 0.5 × 3 × 170 mm; Bio-Rad), whereas cuploading method was used for linear pH 6 to 11 IPG strips (Immobiline DryStrip; 0.5 x 3 x 180 mm; Amersham Pharmacia Biotech). For both IEF methods, 100 μg of protein sample diluted in two-dimensional sample buffer [7 mol/L urea, 2 mol/L thiourea, 4% (w/v) CHAPS, 1.6% (w/v) DTT, 2% (v/v) biolytes, and 1.2% (v/v) DeStreak reagent] was included for each IPG strip. After IEF, protein reduction was performed by equilibration of IPG strip in 0.4% (w/v) DTT for 15 min followed by protein alkylation in 4% (w/v) iodoacetamide for 15 min. The second dimension was performed using 12% polyacrylamide gel (acrylamide/bisacrylamide 37.5:1, 2.6% cross-linker; gel dimension, 20 cm × 20 cm × 1 mm). Two-dimensional gels were fixed for 1 h (25% ethyl alcohol and 7% acetic acid) and stained overnight with the Sypro Ruby protein gel stain (Bio-Rad).

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.

Clinical Specimens

The crypts prepared in this study (Fig. 1A) were epithelial cell type specific that avoided the high contamination of blood and stromal cells found in ordinary colon mucosa preparations (16). The heterozygous nature of the APC mutations in FAP colonic crypt cells were confirmed directly by resequencing the mutation sites of the genomic DNA from crypt cells of four FAP individuals studied in this report (Supplementary Data S1–3; Supplementary Data S4 Table S1). The purity of the crypts enabled a meaningful proteomic analysis, particularly of the lower abundance proteins, a subset that may include many potential biomarkers. Both control and FAP samples were initially used to construct a protein and isoform database for the colon crypts (Fig. 1B; ref. 9). About 1,695 unique protein spots, representing ∼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.

Figure 1.

FAP colonic crypt and its two-dimensional gel proteome. A, a crypt isolated from the morphologically normal colonic mucosa of a FAP patient. B, 3 overlapping pH range two-dimensional gels across pI 4-11 resolve >4,000 unique protein spots of FAP crypt proteome. Two hundred ninety-seven protein spots that changed with a P value of <0.01 are identified by their unique gel spot number in the Supplementary Data S5 to S7 and at our Web site. These maps are searchable and hyperlink-enabled to gene ontology databases.

Figure 1.

FAP colonic crypt and its two-dimensional gel proteome. A, a crypt isolated from the morphologically normal colonic mucosa of a FAP patient. B, 3 overlapping pH range two-dimensional gels across pI 4-11 resolve >4,000 unique protein spots of FAP crypt proteome. Two hundred ninety-seven protein spots that changed with a P value of <0.01 are identified by their unique gel spot number in the Supplementary Data S5 to S7 and at our Web site. These maps are searchable and hyperlink-enabled to gene ontology databases.

Close modal

The files for 3 searchable point-and-click proteome maps of the human colon crypts representing the 214 proteins that separate FAP from control with a P value of <0.01, complete with hyperlinks to gene ontology Web sites, are provided in the Supplementary Data S5 to S7 and at our Web site. This proteomic report complies with the recommendations for the analysis and documentation of peptide and protein identification (17).

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.

Figure 2.

Hierarchical Clustering analysis of the proteomes of FAP crypts versus control (cont) crypts. Data from two-dimensional gels of the following pH ranges: A, pH 4 to 7; B, pH 5 to 8; C, pH 6 to 11; D, pH 5 to 8 with the protein intensities in each column randomized. Red, FAP > control; green, FAP < control; black, no change. There are 5 shades of red and 5 shades of green. The sample identification numbers are listed with each sample.

Figure 2.

Hierarchical Clustering analysis of the proteomes of FAP crypts versus control (cont) crypts. Data from two-dimensional gels of the following pH ranges: A, pH 4 to 7; B, pH 5 to 8; C, pH 6 to 11; D, pH 5 to 8 with the protein intensities in each column randomized. Red, FAP > control; green, FAP < control; black, no change. There are 5 shades of red and 5 shades of green. The sample identification numbers are listed with each sample.

Close modal

The validity of this “unsupervised hierarchical clustering” was tested in Fig. 2D. Here, the spots of each patient of the pH 5 to 8 range in Fig. 2B were scrambled by randomizing the data matrix one column at a time, unlinking the relationship of each protein name to its spot intensity in each gel. The result of the unsupervised hierarchical clustering of this randomized matrix was a disordered display of spot intensities, with consequent lack of separation between FAP and control crypt, thus validating that the order observed in Fig. 2B could not have arisen from statistical chance and bias. Accordingly, the first three components of a PCA cleanly separated FAP and control samples (Fig. 3).

Figure 3.

Three-dimensional display of PCA of the first 3 components of the pH 4 to 11 proteomes of FAP crypts (right) versus control crypts (left). The larger white squares denote the samples, and the smaller squares are individual proteins in the Principal Component space; gender is indicated as M and F. Red squares, 335 proteins that distinguish FAP from control with a P value of <0.01 from two-tailed t test.

Figure 3.

Three-dimensional display of PCA of the first 3 components of the pH 4 to 11 proteomes of FAP crypts (right) versus control crypts (left). The larger white squares denote the samples, and the smaller squares are individual proteins in the Principal Component space; gender is indicated as M and F. Red squares, 335 proteins that distinguish FAP from control with a P value of <0.01 from two-tailed t test.

Close modal

The Genesis one-way ANOVA tool was used to calculate the two-tailed t test P value for each protein in terms of its significance in separating FAP samples from control samples (Supplementary Data S8 Table S2). With this approach, 122, 207, and 65 proteins, total 394, with P values of <0.01 were obtained for pH 4 to 7, 5 to 8, and 6 to 11 range, respectively. Two hundred fourteen of these protein spots have been identified by mass spectrometry.

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.

Figure 4.

Two-dimensional gel Western blotting validation of ECHS1 differential expression in FAP versus control crypts, using antibodies to ECHS1. The immunodetection was obtained by chemiluminescence on film and scanned for densitometry integration of the spots by the Progenesis image analysis software. Spot 1 in control and spot 2 in FAP correspond to the protein spot found by proteomics to be >2-fold elevated in FAP patients.

Figure 4.

Two-dimensional gel Western blotting validation of ECHS1 differential expression in FAP versus control crypts, using antibodies to ECHS1. The immunodetection was obtained by chemiluminescence on film and scanned for densitometry integration of the spots by the Progenesis image analysis software. Spot 1 in control and spot 2 in FAP correspond to the protein spot found by proteomics to be >2-fold elevated in FAP patients.

Close modal
Figure 5.

Two-dimensional gel Western blotting validation of SOD2 differential expression in FAP versus control crypts, using antibodies to SOD2. The immunodetection was obtained by chemiluminescence on film and scanned for densitometry integration of the spots by the Progenesis image analysis software. Spot 2 in control and spot 6 in FAP correspond to a protein spot found by proteomics to be 2-fold elevated in FAP patients.

Figure 5.

Two-dimensional gel Western blotting validation of SOD2 differential expression in FAP versus control crypts, using antibodies to SOD2. The immunodetection was obtained by chemiluminescence on film and scanned for densitometry integration of the spots by the Progenesis image analysis software. Spot 2 in control and spot 6 in FAP correspond to a protein spot found by proteomics to be 2-fold elevated in FAP patients.

Close modal

Effect of Gender or Age

The influence of patient gender on the separation of two patient groups was ascertained by performing a two-tailed t test of male samples versus female samples for the 2,556 protein spots as performed above for FAP versus control, using the Genesis software. A P value of <0.01 was used to obtain a list of protein names that significantly separated male from female. This list was compared with the above list of proteins that separated FAP from control with a P value of <0.01. With the division of the patients according to gender, only 24 identified proteins separated the 2 groups. Of these 24 proteins, none were in common with the 214 identified proteins that separated FAP from control group (Supplementary Data S9).

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:

  • Superoxide dismutase [Mn]: This mitochondria-specific SOD2 (18), was observed as having four isoforms. Spots 2612 (pH 5–8 gel) and 2380 (pH 4–7 gel) were unchanged. Spots 2623 and 2628 were 1-fold and 5-fold higher, respectively, in FAP than in controls (Supplementary Data S8 Table S2).

  • Catalase: Two isoforms of catalase were detected (Supplementary Data S8 Table S2). Spots 1033 and 1024 showed opposing changes in expression in the same pH range gel in FAP versus control crypts, suggesting perhaps that regulation may have occurred by posttranslational modification.

  • Protein disulfide-isomerase A3: Three isoforms of protein disulfide isomerase A3, found in the endoplasmic reticulum, were observed to be elevated in the FAP patient samples (Supplementary Data S8 Table S2).

  • Oncogene DJ-1 (PARK7): Oncogene DJ-1 isoform spot 2319 increased in FAP samples (Supplementary Data S8 Table S2).

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.

Characterization of FAP versus controls.APC plays a dominant role in the etiology of colorectal cancer. Although colorectal adenocarcinomas represent a primary clinical outcome in FAP patients, stromal involvement and stromalepithelial cell interactions are also affected (47). The latter is clearly manifested clinically in the Gardner syndrome variant of FAP where, in addition to colon cancer, desmoid tumors are prevalent (19). It is of interest, however, that in the present study, APC mutations that occur within discrete functional domains of the gene product do not necessarily conform to a clinical phenotype that has been ascribed to these functional domains (2022). APC mutations are also present in ∼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.

No potential conflicts of interest were disclosed.

Note: Supplementary data for this article are available at Cancer Research Online (http://cancerres.aacrjournals.org/).

Grant support: National Cancer Institute, Contracts No. NO1-CN-15103 (A. Knudson) and N01-CN-43309, Grant DK063014 (B. Boman), Institutional Core Grant P30CA006927, CA119242 and the Ewing Trust (A. Yeung), the Fannie E. Rippel Foundation, The Shöller Foundation, Tobacco Settlement Funds from the Commonwealth of Pennsylvania, the Pew Charitable Trust, and the Kresge Foundation.

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.

1
Leading Sites of New Cancer Cases and Deaths 2006 Estimates. 2007 [cited March 29, 2007]. Available from: http://www.cancer.org/downloads/STT/CAFF2006PWSecured.pdf.
2
Fearon ER, Vogelstein B. A genetic model for colorectal tumorigenesis.
Cell
1990
;
61
:
759
–67.
3
Knudson AG. Two genetic hits (more or less) to cancer.
Nat Rev Cancer
2001
;
1
:
157
–62.
4
Kopelovich L. Phenotypic markers in human skin fibroblasts as possible diagnostic indices of hereditary adenomatosis of the colon and rectum.
Cancer
1977
;
40
:
2534
–41.
5
Kopelovich L. Heritable colorectal cancer and cancer genes: systemic expressions.
Mol Carcinog
1993
;
8
:
3
–6.
6
Danes BS. Increased in vitro tetraploidy: tissue specific within the heritable colorectal cancer syndromes with polyposis coli.
Cancer
1978
;
41
:
2330
–4.
7
Stoyanova R, Clapper ML, Bellacosa A, et al. Altered Gene Expression in Phenotypically Normal Renal Cells from Carriers of Tumor Suppressor Gene Mutations. Cancer Biol Ther 
2004
;
3
:
1313
–21.
8
Carrette O, Burkhard PR, Sanchez JC, Hochstrasser DF. State-of-the-art twodimensional gel electrophoresis: a key tool of proteomics research.
Nat Protocol
2006
;
1
:
812
–23.
9
Patel BB, Li XM, Dixon MP, et al. Searchable high-resolution 2D gel proteome of the human colon crypt.
J Proteome Res
2007
;
6
:
2232
–8.
10
Li XM, Patel BB, Blagoi EL, et al. Analyzing alkaline proteins in human colon crypt proteome.
J Proteome Res
2004
;
3
:
821
–33.
11
Eisen MB, Spellman PT, Brown PO, Botstein D. Cluster analysis and display of genome-wide expression patterns.
Proc Natl Acad Sci U S A
1998
;
95
:
14863
–8.
12
Sturn A, Quackenbush J, Trajanoski Z. Genesis: cluster analysis of microarray data.
Bioinformatics
2002
;
18
:
207
–8.
13
Misra J, Schmitt W, Hwang D, Hsiao LL, Gullans S, Stephanopoulos G. Interactive exploration of microarray gene expression patterns in a reduced dimensional space.
Genome Res
2002
;
12
:
1112
–20.
14
Peterson LE. Partitioning large-sample microarray-based gene expression profiles using principal components analysis.
Comput Methods Programs Biomed
2003
;
70
:
107
–19.
15
Al-Shahrour F, Diaz-Uriarte R, Dopazo J. FatiGO: a web tool for finding significant associations of Gene Ontology terms with groups of genes.
Bioinformatics
2004
;
20
:
578
–80.
16
Wildrick DM, Lointier P, Nichols DH, Roll R, Quintanilla B, Boman BM. Isolation of normal human colonic mucosa: comparison of methods.
In Vitro Cell Dev Biol Anim
1997
;
33
:
18
–27.
17
Carr S, Aebersold R, Baldwin M, Burlingame A, Clauser K, Nesvizhskii A. The need for guidelines in publication of peptide and protein identification data: Working Group on Publication Guidelines for Peptide and Protein Identification Data.
Mol Cell Proteomics
2004
;
3
:
531
–3.
18
Fridovich I. Superoxide radical and superoxide dismutases.
Annu Rev Biochem
1995
;
64
:
97
–112.
19
Simpson RD, Harrison EG, Jr., Mayo CW. Mesenteric Fibromatosis in Familial Polyposis. A Variant of Gardner's Syndrome.
Cancer
1964
;
17
:
526
–34.
20
Qian J, Steigerwald K, Combs KA, Barton MC, Groden J. Caspase cleavage of the APC tumor suppressor and release of an amino-terminal domain is required for the transcription-independent function of APC in apoptosis.
Oncogene
2007
;
26
:
4872
–6.
21
Moseley JB, Bartolini F, Okada K, Wen Y, Gundersen GG, Goode BL. Regulated binding of adenomatous polyposis coli protein to actin.
J Biol Chem
2007
;
282
:
12661
–8.
22
Mimori-Kiyosue Y, Matsui C, Sasaki H, Tsukita S. Adenomatous polyposis coli (APC) protein regulates epithelial cell migration and morphogenesis via PDZ domainbased interactions with plasma membranes.
Genes Cells
2007
;
12
:
219
–33.
23
Perucho M. Cancer of the microsatellite mutator phenotype.
Biol Chem
1996
;
377
:
675
–84.
24
Kopelovich L, Conlon S, Pollack R. Defective organization of actin in cultured skin fibroblasts from patients with inherited adenocarcinoma.
Proc Natl Acad Sci U S A
1977
;
74
:
3019
–22.
25
Kopelovich L, Lipkin M, Blattner WA, Fraumeni JF, Jr. Organization of actincontaining cables in cultured skin fibroblasts from individuals at high risk of colon cancer.
Int J Cancer
1980
;
26
:
301
–7.
26
Gonzalez-Quintela A, Garcia J, Campos J, et al. Serum cytokeratins in alcoholic liver disease: contrasting levels of cytokeratin-18 and cytokeratin-19.
Alcohol
2006
;
38
:
45
–9.
27
Barak V, Goike H, Panaretakis KW, Einarsson R. Clinical utility of cytokeratins as tumor markers.
Clin Biochem
2004
;
37
:
529
–40.
28
Prasannan L, Misek DE, Hinderer R, Michon J, Geiger JD, Hanash SM. Identification of β-tubulin isoforms as tumor antigens in neuroblastoma.
Clin Cancer Res
2000
;
6
:
3949
–56.
29
Rock CO, Karim MA, Zhang YM, Jackowski S. The murine pantothenate kinase (Pank1) gene encodes two differentially regulated pantothenate kinase isozymes.
Gene
2002
;
291
:
35
–43.
30
Higashitsuji H, Higashitsuji H, Nagao T, et al. A novel protein overexpressed in hepatoma accelerates export of NF-κ B from the nucleus and inhibits p53-dependent apoptosis.
Cancer Cell
2002
;
2
:
335
–46.
31
Higashitsuji H, Higashitsuji H, Masuda T, Liu Y, Itoh K, Fujita J. Enhanced deacetylation of p53 by the anti-apoptotic protein HSCO in association with histone deacetylase 1.
J Biol Chem
2007
;
282
:
13716
–25.
32
Kuninaka S, Ichinose Y, Koja K, Toh Y. Suppression of manganese superoxide dismutase augments sensitivity to radiation, hyperthermia and doxorubicin in colon cancer cell lines by inducing apoptosis.
Br J Cancer
2000
;
83
:
928
–34.
33
Cao C, Leng Y, Kufe D. Catalase activity is regulated by c-Abl and Arg in the oxidative stress response.
J Biol Chem
2003
;
278
:
29667
–75.
34
Tu BP, Weissman JS. Oxidative protein folding in eukaryotes: mechanisms and consequences.
J Cell Biol
2004
;
164
:
341
–6.
35
Zhou W, Zhu M, Wilson MA, Petsko GA, Fink AL. The oxidation state of DJ-1 regulates its chaperone activity toward α-synuclein.
J Mol Biol
2006
;
356
:
1036
–48.
36
Martinat C, Shendelman S, Jonason A, et al. Sensitivity to oxidative stress in DJ-1-deficient dopamine neurons: an ES- derived cell model of primary Parkinsonism.
PLoS Biol
2004
;
2
:
e327
.
37
Takahashi K, Taira T, Niki T, Seino C, Iguchi-Ariga SM, Ariga H. DJ-1 positively regulates the androgen receptor by impairing the binding of PIASx α to the receptor.
J Biol Chem
2001
;
276
:
37556
–63.

Supplementary data