| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Epidemiology and Prevention |
1 Laboratory of Proteomics and Analytical Technologies and 2 Advanced Biomedical Computing Center, SAIC-Frederick Inc., National Cancer Institute at Frederick, Frederick, Maryland; 3 Division of Cancer Prevention, National Cancer Institute, Bethesda, Maryland; 4 University of Texas M. D. Anderson Cancer Center, Houston, Texas; and 5 Imperial Cancer Research Fund, St. Marks Hospital, London, United Kingdom
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The cancer prevention trial using celecoxib in FAP patients in a double-blinded, placebo-controlled study (14) demonstrated a significant reduction both in the number and burden of polyps. Even so, the individual response of patients in this trial was remarkably heterogeneous. An extensive profiling of the molecular circuitry affected by celecoxib would be helpful in understanding the heterogeneity of patients responses in this trial. Currently, high throughput proteomic technologies are being used to identify diagnostic markers for early cancer detection and plausible therapeutic targets responsive to specific treatments. In particular, surface-enhanced laser desorption/ionization time of flight mass spectrometry (SELDI-TOF MS) is being explored extensively for its use in diagnostic medicine. This technique circumvents sample purification/fractionation procedures in a direct one-step affinity capture of a subset of proteins from highly complex biological fluids, such as serum and plasma (15 , 16) , which are believed to reflect the physiological as well as pathological molecular events. Recently, SELDI-TOF MS, together with bioinformatics tools, has been successfully used for differential serum proteomic profiling of various cancers. In each case, a panel of markers was identified with the ability to discriminate between individuals with no cancer, cancer, and among various stages of cancer development (17, 18, 19, 20, 21, 22) . These proteomic profiling approaches can be instrumental in identifying molecular signatures of cancers; furthermore, these can also be used in monitoring molecular alterations that ensue after any intervention with the potential of impacting drug development.
We used SELDI-TOF MS to obtain proteomic profiles of sera from patients of the FAP/celecoxib clinical trial with the assumption that celecoxib-modulated changes may be reflected in the alteration of serum protein components. The objectives of this study were to delineate serum proteomic changes in the FAP patients after receiving celecoxib and to identify serum markers to predict response to treatment with celecoxib.
| MATERIALS AND METHODS |
|---|
|
|
|---|
SELDI Proteomic Profiling.
Weak cationic exchange (WCX) and strong anionic exchange (SAX) protein chips (Ciphergen Biosystems, Inc., Fremont, CA) were used to obtain protein profiles from serum samples analyzed in duplicate. SAX chips were first pre-equilibrated twice with 200 µl of binding buffer [20 mM Tris (pH 8.0), 0.1% Triton X-100] and agitated for 5 min each time. Equilibration of WCX chips consisted of treating with 100 µl of 10 mM HCl for 5 min on a shaker, rinsing twice with deionized H2O for 5 min, and coating with 150 µl of binding buffer [50 mM sodium acetate (pH 4.5), 0.1% Triton X-100] for 10 min. Serum samples (5 µl of each) were thawed on ice, diluted to 100 µl in respective binding buffers, applied on the chip arrays, and incubated on a shaker for 1 h. The chips were then washed with 200 µl of the respective binding buffers three times for 5 min each. After the final wash, chips were rinsed twice with deionized H2O for 30 s each and air-dried.
Before SELDI-TOF MS analysis, 0.5 µl of saturated solution of sinapinic acid (Fluka, Milwaukee, WI) in 50% (v/v) acetonitrile, and 0.5% trifluoroacetic acid was applied onto each spot twice and air-dried between applications. Chips were then placed in the PBS-II SELDI-TOF mass spectrometer (Ciphergen Biosystems, Inc.) operated in the positive ion mode. Time-of-flight spectra were generated by averaging 70 laser shots collected on each spot with laser intensity setting of 220, detector sensitivity setting of 9, and a lag time focusing of 900 ns. The spectra were calibrated using the All-in-1 protein molecular mass standard (Ciphergen Biosystems, Inc.). To compensate for slight spot-to-spot variations, if any, the spectra were also normalized using the total ion current method in the m/z range of 1,500 to 50,000 (SAX surface) or 1,500 to 20,000 (WCX surface) with subtracted baseline.
Peaks were identified after mass calibration, background subtraction, and normalization using the clustering and alignment function of ProteinChip software 3.0 (Ciphergen). The settings used were as follows: signal/noise ratio in the first pass, 3; minimum peak threshold, 5%; cluster mass window, 0.3%; signal/noise in the second pass, 1.5. The peaks detected within ±0.15% mass to charge (m/z) units of each other across the spectra were considered as one cluster, and a particular cluster was represented by its average m/z value. The peak information including m/z and intensity values was exported into tables in Microsoft Excel for statistical analysis.
Statistical Analysis.
For pre- versus post-treatment analysis, four spectra, duplicates of pre- and post-treatment samples, from each of the 31 patients enrolled at M. D. Anderson Cancer Center were used. The logarithm of the intensities for each peak cluster in a given treatment arm was used in a linear mixed effects model (23)
to construct P in t test. The celecoxib-modulated changes were identified as the peaks with P < 0.05 in both the 100- and 400-mg treatment arms and >0.7 in the placebo arm.
To find response predictors, 55 patients whose pretreatment sera were available for analyses, were classified into four groups (please see text and Table 3
for details). The random forest (RF) classification algorithm (24, 25, 26)
was used to search for serum markers (peak clusters) that distinguish between responders and nonresponders. This method uses an ensemble of 500 decision trees that "vote" on the classification of a sample. The method also prevents over-fitting without relying on tuning parameters, which are common problems associated with other standard decision tree algorithms. Each tree was grown to full-length without pruning, but only a random subset of the samples was used. For example, the group 1 set contains 110 samples (two spectra from all 55 patients) from which a random subset was chosen to construct each decision tree. This set was chosen from the full set with replacement until it contained 110 samples. This procedure selects approximately two thirds of the samples to build each tree. The excluded, or "out-of-bag," samples became the test set for that particular tree. The algorithm also used random feature selection to prevent false inflation of the sensitivity and specificity estimates and to give more independent weight to each trees vote. Specifically, 15% of the features were randomly chosen each time a node was to be split, and the selection of the best feature was made within this random subset and not from among all features. Because a small subset of the features was examined each time a node in a decision tree needed to be split, the number of splits needed to grow the tree to full-length, and therefore the number of features used in the construction, varied from tree to tree. This meant that approximately all of the features had a chance of being used at least once in the collection of 500 decision trees, and the "best overall feature" had a reasonable chance of not being used until a later split or not being used at all.
|
| RESULTS |
|---|
|
|
|---|
Because the serum from each patient was analyzed in duplicate, a linear mixed effects model (23) was used to remove the random variation in each patients duplicate spectra from the total variation between pre- and post-treatment spectra thereby obtaining a better estimate of the effect of the treatment. The peak clusters with P < 0.05 were identified from both chip surfaces. The WCX spectra contained seven significant peaks from the placebo arm and 10 and 26 peaks from the 100-mg bid and 400-mg bid arms, respectively. The spectra from the SAX chip produced three, nine, and seven significant peaks for the respective treatments. The differentially expressed markers in the placebo arm probably represented time-related changes between pre- and post-treatment (data not shown).
Six differentially expressed peak clusters from the WCX chip and two from the SAX chip were common to both 100- and 400-mg arms and absent from the placebo arm. The abundance level of each of the eight clusters followed the same direction in pre- versus post-treatment change with both doses of celecoxib and on both chip surfaces (Tables 1
and 2
). Fig. 1, A and B
, show examples of two patients each on the WCX and SAX chips, respectively, with the gel views and trace views of increased expression of the marker at m/z 3,158.7/3,158.8 and reduced expression of the marker at m/z 5,089.5/5,093.6 before and after treatment with celecoxib.
|
|
|
|
The classification yielded 94 and 97% sensitivity but only a 52 and 50% specificity for group 1 and group 2, respectively. The low specificity could be attributable to the potentially erroneous classification of subjects in the placebo arm and nonresponders in the 100-mg treatment arm. In group 3, the algorithm classified all 30 responders correctly, but three of the nonresponders were misclassified. This reduced the specificity to 40%, but the sensitivity and negative-predictive values increased to 100%, and the positive-predictive value increased to 91%. The low specificity of 40% may be attributable to the highly imbalanced numbers of responders (n = 30) versus nonresponders (n = 5) producing a 6:1 ratio of responders to nonresponders in this group. This means that the randomly chosen samples used to build each decision tree were very heavily weighted toward responders. Any nonresponders in the out-of-bag samples were then tested in their corresponding responder-biased decision trees, and this is the most likely reason why two of the nonresponders were misclassified in the group 3 run. When the responders and nonresponders of the 400-mg arm alone were examined (group 4), one subject from each set was misclassified. This produced a sensitivity of 94%, a specificity of 80%, and positive- and negative-predictive values of 94 and 80%, respectively.
The peaks with sufficiently high importance scores are also listed in Table 3
. It is noticeable that the classifiers for groups 1 (m/z 7,567.48 and 7,928.24) and 2 (m/z 7,928.24 and 3,975.82) and groups 3 (m/z 16,961.4) and 4 (m/z 16,961.4 and 4,647.23) are different with one shared peak each between groups 1 and 2 (m/z 7,928.24) and groups 3 and 4 (m/z 16,961.4). The fact that the important spectral features were different for groups 1 and 2 than groups 3 and 4 supports the contention that the algorithm was simply trying to fit less precise data. The degradation in the assigned class of the samples resulted in new features being chosen as the most important ones.
The peak cluster m/z 16,961.4 with a very high importance score is an important response discriminator for both groups 3 and 4. Fig. 3A
displays the region of the SAX spectra between m/z 16,500 and 17,500 from two strong responders and nonresponders. The spectra exhibit the presence of m/z 16,961.4 in nonresponders and its absence in the strong responders. Although there is a pattern of several differentially expressed markers of much higher intensities in this region, Fig. 3B
demonstrates that the peak cluster at m/z 16,961.4 is by far the best discriminator of celecoxib response in group 4 patients treated with a 400-mg bid dose. A separate classification algorithm using a distance-dependent K-nearest neighbor algorithm also demonstrated that the m/z 16,961.4 SAX peak was an important discriminator between responders and nonresponders (data not shown).
|
| DISCUSSION |
|---|
|
|
|---|
In this study, we, for the first time, describe serum proteomic profiling to analyze drug response in a clinical trial of cancer prevention. We examined serum proteomic changes in FAP patients in response to celecoxib and attempted to identify predictive molecular patterns for the response or nonresponse of patients. Comparative profiling of sera from pre- and post-treatment samples identified a panel of markers with specific m/z ratios, the expression of which was significantly modulated by celecoxib. More importantly, these proteomic markers were over- or underexpressed in both treatment arms with low and high doses of celecoxib suggesting that at least some of them are relevant to the celecoxib-mediated regression of adenomas. Information on the identities of these celecoxib-modulated proteomic markers is expected to provide an understanding of the complex cellular networks affected by celecoxib, which may act in multiple ways to restrain carcinogenic process via COX-2-dependent as well as COX-2-independent mechanisms.
The second objective of our study was to use pretreatment sera for the identification of markers, which may serve as classifiers for response to celecoxib. The RF algorithm identified two proteomic markers, m/z 16,961.4 from the SAX chip and m/z 4,647.23 from the WCX chip, especially m/z 16,961.4 with a very high importance score, that could distinguish between responders and nonresponders of the treatment arm receiving 400 mg of celecoxib with misclassification of only one subject each of the responders and nonresponders (Table 3)
. Fig. 4
shows that the features m/z 16,961.4 and 4,647.23 can effectively separate the responders from the nonresponders in groups 3 and 4. The fact that m/z 16,961.4 and 4,647.23 markers from the pretreatment sera are recognized by the RF algorithm as important classifiers for response to celecoxib and are not among the markers in Tables 1
and 2
suggests that they are relevant for the chemopreventive effect of celecoxib without being significantly modulated by it. Identification and characterization of the proteomic markers that are modulated by celecoxib and/or aid in response prediction are clearly warranted.
|
Our study is a significant first attempt in using proteomic tools to predict response outcome in a cancer prevention trial. It also addresses the issue of inter-patient variability in response resulting in variable benefit of the same drug to individual patients. Identification of molecular predictors would be beneficial in maximizing potential response to specific prevention/treatment agents and simultaneously minimizing undesired side effects. The efficacy of celecoxib is currently being evaluated in a much larger clinical trial of sporadic colorectal adenomas. The use of celecoxib for the prevention and treatment of a variety of other cancers is also being investigated. Serum expression profiling to identify drug-modulated proteomic markers and clinical response predictors may revolutionize treatment strategies based on our ability to predict response of individual patients to various treatments. There is certainly a need for predictive molecular factors if personalized medicine is to become a reality.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
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.
Note: By acceptance of this article, the publisher or recipient acknowledges the right of the United States Government to retain a nonexclusive, royalty-free license and any copyright covering the article. The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does mention of trade names, commercial products, or organization imply endorsement by the United States Government.
Requests for reprints: Iqbal U. Ali, Division of Cancer Prevention, National Cancer Institute, Bethesda, MD 20893. E-mail: alii{at}mail.nih.gov
Received 12/ 2/03. Revised 2/23/04. Accepted 2/23/04.
| REFERENCES |
|---|
|
|
|---|
716 knockout mice by inhibition of cyclooxygenase 2 (COX-2). Cell, 87: 803-9, 1996.[CrossRef][Medline]
This article has been cited by other articles:
![]() |
P. J. O'Dwyer, S. G. Eckhardt, D. G. Haller, J. Tepper, D. Ahnen, S. Hamilton, A. B. Benson III, M. Rothenberg, N. Petrelli, H.-J. Lenz, et al. Priorities in Colorectal Cancer Research: Recommendations From the Gastrointestinal Scientific Leadership Council of the Coalition of Cancer Cooperative Groups J. Clin. Oncol., June 1, 2007; 25(16): 2313 - 2321. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. S. Tan, A. Ploner, A. Quandt, J. Lehtio, and Y. Pawitan Finding regions of significance in SELDI measurements for identifying protein biomarkers Bioinformatics, June 15, 2006; 22(12): 1515 - 1523. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Fosslien Cardiovascular Complications of Non-Steroidal Anti-Inflammatory Drugs Ann. Clin. Lab. Sci., October 1, 2005; 35(4): 347 - 385. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. J. Robbins, J. Villanueva, and P. Tempst Distilling Cancer Biomarkers From the Serum Peptidome: High Technology Reading of Tea Leaves or an Insight to Clinical Systems Biology? J. Clin. Oncol., August 1, 2005; 23(22): 4835 - 4837. [Full Text] [PDF] |
||||
![]() |
S. M. Lippman and B. Levin Cancer Prevention: Strong Science and Real Medicine J. Clin. Oncol., January 10, 2005; 23(2): 249 - 253. [Full Text] [PDF] |
||||
![]() |
E. T. Hawk and B. Levin Colorectal Cancer Prevention J. Clin. Oncol., January 10, 2005; 23(2): 378 - 391. [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 |