|
|
||||||||
Proteomics and Protein Markers |
1 British Columbia Cancer Agency, Vancouver, British Columbia, Canada.
2 Ciphergen Biosystems, Inc., Fremont, CA.
aAddress correspondence to this author at: Genomics Sciences Centre, British Columbia Cancer Agency, 600 West 10th Ave., Vancouver, British Columbia, V5Z 4E6 Canada. Fax 604-877-6011; e-mail msadar{at}bccrc.ca.
| Abstract |
|---|
|
|
|---|
Methods: We used SELDI-TOF MS for protein profiling of sera from prostate cancer patients (n = 38) with and without bone metastases in our effort to identify individual or multiple serum markers that may be of added benefit to those in current use. Serum was applied to ProteinChip® surfaces (H4 and IMAC) to quickly screen samples and detect peaks predominating in the samples obtained from patients with bone metastases. Unique proteins in the bone metastasis cohort observed by SELDI-TOF MS were identified by two-dimensional gel electrophoresis, in-gel trypsin digestion, and tandem MS. The identities of the proteins were confirmed by ELISA and immunodepletion assays.
Results: The cluster of unique proteins in the sera of patients with bone metastases was identified as isoforms of serum amyloid A. Machine-learning algorithms were also used to identify patients with bone metastases with a sensitivity and specificity of 89.5%.
Conclusions: SELDI-TOF MS protein profiling in combination with other proteomic approaches may provide diagnostic tools with potential clinical applications and serve as tools to aid in the discovery of biomarkers associated with various diseases.
| Introduction |
|---|
|
|
|---|
20% of cases exhibit a combination of blastic and lytic (bone destroying) lesions. Lytic lesions occur in 1530% of patients (5). Bone metastases are associated with a poor prognosis and with the development of hormone-refractory disease, for which there is currently no effective treatment (6)(7)(8)(9). Once the disease is hormone refractory, the typical survival time is <2 years. Bone tissue is dynamically maintained by two tightly coupled processes of resorption and formation, which are controlled by osteoclasts and osteoblasts, respectively. Bone remodeling liberates potential biomarkers. Examples of biomarkers for resorption are associated with collagen degradation, such as the carboxy-terminal telopeptides of type I collagen, pyridinoline, and deoxypyridinoline. Markers for bone formation, or osteogenesis, include the carboxy-terminal propeptide of type I collagen, which is cleaved from procollagen molecules during the formation of type I collagen; alkaline phosphatase; and osteocalcin. Bone formation involves three steps: The first is synthesis of collagen, which is associated with increased proliferation and increased concentrations of carboxy-terminal propeptide of type I collagen. The second phase involves matrix maturation, at which time concentrations of alkaline phosphatase increase. The third phase involves mineralization, which is marked by increased concentrations of osteocalcin. Studies examining the potential usefulness of biochemical markers for the detection of bone metastasis have generally focused on markers for bone turnover in prostate cancer patients and have shown a lack of sensitivity compared with bone scintigraphy scans.
Serum concentrations of prostate-specific antigen (PSA)1 are useful in combination with a physical examination for the initial diagnosis of prostate cancer, but when used alone, PSA is associated with a high degree of false positives. Increased PSA concentrations in the serum of men with prostate cancer have been suggested to be directly correlated to the tumor volume (10)(11). One of the best applications of serum PSA measurements is for the monitoring of treatment response and relapse in prostate cancer patients. However, serum PSA concentrations cannot predict the presence of bone or soft tissue metastases (12). At the molecular level, PSA gene expression is regulated by the androgen receptor via androgen response elements (13). Immunohistochemical studies confirmed the presence of an androgen receptor in metastatic prostate cancer biopsies (14)(15).
Bone-derived factors such as interleukin-6 (IL-6) (16)(17)(18)(19)(20) and insulin-like growth factor-I (21) have been shown to stimulate androgen-responsive reporters or PSA gene expression in the absence of androgen [for a review, see Ref. (22)]. IL-6 is also increased in the sera of patients with metastatic prostate cancer (23)(24)(25) and hormone-refractory disease(26), and may be a surrogate marker for androgen-independent disease. Other potential markers of prostate cancer currently being investigated include cadherin (27), human kallikrein 2(28), and chromogranin A (29).
Despite a growing list of potential markers of prostate cancer, sensitive assays to detect, and ultimately predict susceptibility to, metastases by use of a blood test are still urgently required. Serum concentrations of markers for bone turnover, such as alkaline phosphatase, become increased in response to systemic therapy (30), reducing the clinical application of this biomarker for the monitoring of prostate cancer patients, particularly in the case of those patients potentially receiving intermittent androgen suppression. Detection of bone lesions is generally made by clinical symptoms and bone scans, but the latter may not reflect response to therapy (31). Identification of a biomarker for bone or other metastatic lesions that could be measured concurrently with PSA in the blood of patients would enable more accurate clinical assessment of the stage of the disease and may also be of use in the stratification of patients with respect to anticipated or actual response to therapy where current methods lack reliability. We applied surface-enhanced laser desorption/ionization time-of-flight mass spectrometry (SELDI-TOF MS) profiling to sera from prostate cancer patients with and without positive bone scans in our effort to identify individual or multiple serum markers that may be complementary to those in current use as well as to provide further insights into the mechanisms governing the apparent prostate/bone symbiosis. Here we show that rapid SELDI-TOF MS profiling, combined with MS identification of proteins that were isolated by two-dimensional polyacrylamide gel electrophoresis designed to resolve the observed cluster seen with SELDI-TOF MS, correlated with the bone metastases cohort. The cluster of proteins were identified as isoforms of serum amyloid A (SAA). Immunologic quantification validated a correlation of SAA within the bone metastasis group. Machine-learning algorithms were used to identify patients with bone metastases. To our knowledge, this is the first study that compares SELDI-TOF MS spectra, ELISA, and machine-learning algorithms to differentiate metastatic disease.
| Materials and Methods |
|---|
|
|
|---|
|
seldi-tof ms analysis of serum proteins
Aliquots of serum samples were thawed on ice, and protein concentrations were measured using the bicinchoninic acid assay (Pierce). All samples for MS analysis were diluted to a final concentration of 1 g/L in HEPES (pH 7.4), 30 mmol/L NaCl. Diluted serum samples were analyzed on both reversed-phase (hydrophobic surface) H4 and Immobilized Metal Affinity Capture (IMAC) ProteinChip Arrays (Ciphergen Biosystems, Inc.).
Reversed-phase (acyl chains of 16 methylene groups) H4 ProteinChip Arrays were preactivated with 2 µL of acetonitrile (500 mL/L) in a humidity chamber for 5 min at room temperature. The acetonitrile was removed, and spots were allowed to briefly air dry before application of diluted serum samples (2 µL). Samples were allowed to dry on the spot, thereby forcing an interaction with the reversed-phase surface. Spots were then washed twice with 2 µL of acetonitrile (300 mL/L). A saturated solution of sinapinic acid in acetonitrile (500 mL/L)trifluoroacetic acid (TFA; 5 mL/L) was then applied twice (0.5 µL each time) as the chemical matrix and allowed to air dry between applications.
IMAC3 (nitrilotriacetic acid functional group) ProteinChip Arrays were assembled into a bioprocessor (Ciphergen Biosystems), a device that holds 12 ProteinChip Arrays and allows application of sample volumes up to 350 µL, before spots were charged with 100 µL of CuSO4 (50 mmol/L) for 5 min. The ProteinChip Arrays were then washed four times with 100 µL of water, fixed with 100 µL of 100 mmol/L sodium acetate, 200 mmol/L NaCl, and equilibrated with three 100-µL volumes of phosphate-buffered saline, pH 7.2. Diluted serum samples (20 µL) were denatured with 30 µL of 8 mol/L urea, 10 g/L CHAPS and vortex-mixed for 15 min at 4 °C. To lower the urea concentration, we added 100 µL of 1 mol/L urea, 1.25 g/L CHAPS and vortex-mixed the sample. Finally, we added 600 µL of binding buffer (phosphate-buffered saline) to the denatured serum samples and loaded 50 µL per spot. The sample was allowed to bind to the ProteinChip Array for 30 min with shaking at room temperature. A 30-min binding time has previously been found to be optimal with no statistically significant gain in information obtained when binding times are increased up to 2 h in length (data not shown).
The spots were washed with three 50-µL volumes of phosphate-buffered saline to remove proteins interacting nonspecifically with the IMAC-Cu++ surface, followed by five 50-µL washes with water to remove salts. Arrays were then air-dried, and two 0.5-µL volumes of sinapinic acid (prepared as described previously) were added to each spot and allowed to air dry before mass analysis using a ProteinChip Reader (Ciphergen Biosystems).
The ProteinChip reader, or linear SELDI-TOF MS, has delayed extraction (time lag focusing) and a mass analysis range from <0.1 to >500 kDa. The sensitivity is better than 700:1 mean signal-to-noise ratio on analysis of 150 fmol of bovine IgG when digitally smoothed with a (0.2 times peak width) sliding window average. Mass accuracy with internal calibration is better than ± 200 ppm (0.02%) at m/z 100010 000 and better than ± 500 ppm (0.05%) at m/z 10 000300 000. The mass accuracy with external calibration is better than ±1000 ppm (0.1%) at m/z 100010 000 and better than ±2000 ppm (0.2%) at m/z 10 000300 000. Spectra were generated by use, on average, of 90105 laser shots using laser intensities of 215 with a sensitivity of 8 (IMAC-Cu++) or laser of 250 with a sensitivity of 9 (H4). The instrument was externally calibrated by use of bovine ubiquitin, cytochrome C, superoxide dismutase, ß-lactoglobulin, and horseradish peroxidase (Ciphergen Biosystems) with matrix application identical to that for the serum samples. Before data analysis, spectra were baseline subtracted by use of a modified piecewise convex-hull algorithm, which attempts find the bottom of the spectrum and correct the peak height and area, and normalized by use of total ion current.
two-dimensional gel electrophoresis
Each serum sample (500 µg of protein) was solubilized in 8 mol/L urea, 40 mL/L CHAPS, 40 mmol/L Tris, 2 mL/L Bio-Lyte 3/10 ampholytes (Sequential Extraction Buffer 2; Bio-Rad) with 5 mmol/L tris(2-carboxyethyl)phosphine reducing agent and incubated for 1 h at room temperature. The samples were loaded on 7-cm immobilized pH gradient strips (ReadyStrip IPG strips; linear gradient, pH 310; Bio-Rad). Isoelectric focusing of samples (Bio-Rad Protean Cell) was carried out as follows: 250 V for 15 min, followed by a slow voltage ramp over 8000 V-h to 4000 V, with a maximum of 50 µA per strip. Voltage was held at 4000 V, with a maximum of 50 µA per strip, for 20 000 V-h, followed by a 500-V hold. Before the second dimension, the strips were reduced and alkylated for 20 min in 62.5 mmol/L Tris (pH 8.45), 5 mmol/L tris(2-carboxyethyl)phosphine, 25 g/L acrylamide, 200 mL/L glycerol, 20 g/L sodium dodecyl sulfate (SDS), 6 mol/L urea. The second dimension consisted of a 16.5% Separation Gel (4.5 cm), 10% spacer gel (1 cm), and a 4% stacking gel (32). The IPG strips were held in place with 1% agarose in 62.5 mmol/L Tris (pH 8.45), 2 g/L SDS. The second dimension was run at 150 V for
60 min with 100 mmol/L Tris, 100 mmol/L tricine, 10 g/L SDS as cathode buffer (pH 8.3) and 200 mmol/L Tris (pH 8.9) as the anode buffer. The gels were fixed in 500 mL/L ethanol, 30 mL/L phosphoric acid for 2 h, followed by washes with MilliQ water (320 min). The gels were equilibrated for 1 h in Neuhoff solution (33) (250 mL/L methanol, 50 mL/L phosphoric acid, 160 g/L ammonium sulfate) before addition of 0.1 g of Serva Blue G per 100 mL and stained for 24 h. Gels were stored in 200 g/L ammonium sulfate at 4 °C until bands were excised.
identification of proteins by ms
Protein bands plus a blank region from the edge of the gel were excised from two-dimensional gels. Individual spots were washed several times for at least 2 h with 500 mL/L methanol, 50 mL/L acetic acid, reduced with 10 mmol/L dithiothreitol, and alkylated with 100 mmol/L iodoacetamide before in-gel trypsin digestion with 10 µL of 20 µg/L trypsin (Boehringer) overnight at 37 °C. Supernatants were then collected and gel pieces extracted with additional ammonium bicarbonate (10 µL of a 25 mmol/L solution) for 15 min, which was pooled with the original supernatant (sample). The samples were zip-tipped (Millipore uC18); washed with 1 mL/L TFA in water; eluted with 2 µL of 500 mL/L acetonitrile, 2,5-dihydroxybenzoic acid, 1 mL/L TFA onto a 100-well stainless steel plate (Applied Biosystems); and analyzed via matrix-assisted laser desorption/ionization (MALDI) on a QStar Pulsar I (Applied Biosystems/MDS SCIEX) equipped with a MALDI I ion source. The QStar was calibrated by use of two tandem MS (MS/MS) fragments of Glu-fibrinogen peptide (Sigma) with [M+H]+ masses of 175.11195 and 1570.6774. After manual inspection of the survey scan to prevent the selection of keratin and trypsin autolysis peaks, the most intense ions were selected for MS/MS fragmentation with an operator-controlled collision energy. All spectra were acquired in positive-ion mode. All product ion scans were acquired from 100 m/z to an upper range that included the peptide mass selected for MS/MS fragmentation. The most intense peptide peaks were selected for fragmentation, and analysis of the resulting data was carried out with the Mascot search engine.
psa, alkaline phosphatase, il-6, and saa
PSA was detected with the enzymatic immunoassay (Abbott IMx). Alkaline phosphatase was measured by an absorbance assay using the Dimension® clinical chemistry system (Dade Behring). SAA concentrations in diluted serum samples were measured in triplicate by ELISA (Biosource International), according to the manufacturers instructions. The detection limit of the assay is 5 µg/L, and all SAA isoforms were detected. IL-6 concentrations in serum were measured in triplicate by ELISA according to the manufacturers instructions (Biosource International). The detection limit of the assay is <2 ng/L. The assay detects both free and receptor-associated IL-6.
western blot analysis and immunodepletion assay
A representative sample from each group was separated by 20% Tris-tricine-SDSpolyacrylamide electrophoresis and transferred to an Immobilon-P membrane (Millipore) in 25 mmol/L Tris, 192 mmol/L glycine, 100 mL/L methanol. Western blotting was performed on 31 µg/lane. The immunoblot was blocked with 50 g/L milk in 20 mmol/L Tris, 500 mmol/L NaCl, 3 mL/L Tween 20 (Tris-buffered salineTween) for 1 h and then incubated overnight at 4 °C with 0.83 µg/L antibody to SAA (Biosource International). The immunoblot was washed and incubated for 1 h with a 1:5000 dilution of secondary antibody (SC-2006; Santa Cruz Biotechnology). The antibodies were diluted in 50 g/L milk in Tris-buffered salineTween. All incubations were performed at room temperature unless otherwise stated. The SAA protein was detected by enhanced chemiluminescence (Amersham).
A sample from each patient group was depleted of SAA by immunoprecipitation using the Seize Primary Immunoprecipitation Kit (Pierce Endogen). Each sample was immunodepleted with an 8 molar excess of anti-SAA (Biosource International) coupled to AminoLink Plus Coupling Gel (Pierce Endogen), according to manufacturer instructions. The eluate (depleted serum) and washes were collected and subsequently analyzed by SELDI-TOF MS using a WCX2 (weak cation-exchange surface) ProteinChip Array.
support vector machine analysis
Spectra obtained from IMAC-Cu++ ProteinChip Arrays were analyzed by mapping the raw (nonfiltered and non-baseline-subtracted) TOF spectra to mass spectra consisting of 16 384 channels with a mass calibration given by m/z = aC2, where C is the channel number and a = 0.0001 Da. Masses above 26.8 kDa were not considered in the analysis. All mass spectra were then normalized to the same total area. An automated procedure to find and fit the significant peaks in the mass spectra was developed and implemented in the C language. The resulting open-source X Windows application, named Nonprofit (for nonlinear protein fit), should compile and run on any UNIX operating system. It is freely available and can be obtained by contacting one of the authors (sflibotte{at}bcgsc.ca). In brief, a spectrum representing the average of all samples was first divided in several sections, the width and position of each section being automatically selected by the software, by use of a heuristic approach to obtain the best possible fits in the following step. Each section was then fitted with the appropriate number of gaussian peaks superimposed on a quadratic background, the number of peaks producing the best fit being found with an iterative procedure. The positions, widths, and heights of the peaks were fitted at the same time as the background parameters. The relative peak width was fixed, which means that all peaks in a given section had the same width but that this width was allowed to vary in the fitting procedure. For each section, the result of the fit of the average spectrum was then used as a template to set the initial parameters for fitting all spectra in the data set, the number peaks and the relative positions of the peaks being fixed. This procedure assumed that all peaks to be fitted appeared in the average spectrum, which is a reasonable assumption because it is unlikely that a peak appearing in only a few spectra would be useful in a classification scheme. This procedure corrected for an imperfect gain matching between spectra because the absolute positions of the peaks were free to vary during the fit; only their relative positions were fixed. For each spectrum, the position and area of each gaussian peak were saved for subsequent analysis and data mining. Such a fitting procedure has the advantage of being able to extract meaningful abundance even for proteins associated with mass peaks that were not completely resolved. Examples of fits obtained with the Nonprofit software are shown in Fig. s1 of the Data Supplement that accompanies this article at http://www.clinchem.org/content/vol51/issue4/.
Support vector machines (SVMs) were used to determine whether machine-learning algorithms could differentiate prostate cancer patients with bone lesions from those without bone lesions. Peaks under m/z 975 were not fitted because many of them showed obvious signs of saturation. The automated procedure applied divided the spectra into 39 sections for fitting, and a total of 270 peaks were fitted for each spectrum, to give >10 000 gaussian peaks. With the help of the fitted peak areas, a classifier was built to distinguish between the spectra associated with the patients who tested positive for bone metastasis and those who did not. SVM was chosen because this learning scheme performed very well as a classifier in association with our peak-fitting procedure on larger SELDI-TOF MS data sets publicly available (data not shown). Furthermore, SVMs are known to perform well in the challenging situation in which the number of samples in the data set is not large compared with the number of attributes. For example, SVMs have been successfully used as a classification scheme in microarray experiments (34). In our case, the attributes were the peak areas. Suppose that n attributes are being considered to build a classifier; each sample (or instance in the machine-learning terminology) can then be represented by a point in an n-dimensional space. SVMs find the maximum margin hyperplane, which is the hyperplane separating the two classes of samples in this n-dimensional space while maximizing the distance between the hyperplane and the closest training point. With our fitting procedure on SELDI-TOF MS data, we found that it is not usually necessary to select the attributes to reduce the dimensionality before building a good classifier with SVM. We therefore do not have to select which type of attribute filter or wrapper would be the most efficient. In practice, we have used the sequential minimal optimization algorithm as implemented in the Weka machine-learning software (35). To maximize the use of this relatively small data set, we have assessed the performance of the classifier with a leave-one-out cross-validation scheme. In other words, 37 of the samples were used to build a classifier and predict the class of the remaining sample, and this procedure was repeated 38 times, leaving out each sample one at a time.
statistical analysis
For our purposes, specificity is defined as the ratio of the patients who do not have the protein peak (or above-normal concentrations) and do not have bone metastases to the total number of patients that do not have bone metastases. Sensitivity is defined as the ratio of the patients with bone metastases who had the biomarker (or above-normal concentrations) to the total number of patients with bone metastases. Comparisons between patients with and without bone metastases were performed with the Welch T-test and the nonparametric Wilcoxon test. The latter test does not assume that the variables follow gaussian distributions. Statistical significance was assumed when P was <0.05.
Univariate statistical analysis of SELDI-TOF MS peak masses and relative intensity values was conducted with the nonparametric MannWhitney U-test.
| Results |
|---|
|
|
|---|
m/z 11 700 that are predominantly detected in samples from patients with bone metastases (lanes 2738). This cluster was also detected in two samples obtained from patients not considered to have bone metastases (lanes 18 and 19). Enhanced analysis of this region of the spectra for lane 4 (without bone metastasis) and lane 37 (with bone metastasis) is shown in shown in Fig. 1B
|
Statistical analysis of peak intensities of the H4 ProteinChip array spectra revealed that the differences observed of
11.7 kDa in sera of patients with bone metastases were statistically significantly increased [m/z 11 488 (P = 0.0020); m/z 11 537 (P = 0.0004); m/z 11 639 (P = 0.0022); m/z 11 680 (P = 0.0051); Fig. 2
]. A plot of the peak intensities of the neighboring statistically insignificant (P = 0.3812) m/z 13 885 peak is included for reference. The dominant peak of the four peaks is
11.7 kDa. The patients with bone metastases who had this cluster had variable Gleason scores (510), current PSA concentrations (402600 µg/L), and alkaline phosphatase (152726 U/L). The Gleason score is the sum of the numbers assigned to the two glandular patterns that predominate in a biopsy of the prostate malignancy and represents the microscopic appearance of the tumor (38). Tumors with Gleason scores of 24 are considered less aggressive than tumors assigned a Gleason score of 810. The two patients not considered to have bone metastases who had this cluster at
11.7 kDa (lanes 18 and 19) had Gleason scores of 7 and 9, respectively.
|
identification of the cluster of proteins by two-dimensional gel electrophoresis and ms
Preparative two-dimensional electrophoresis was used to isolate the cluster of proteins observed at 11.7 kDa in the SELDI-TOF MS spectra. Serum samples chosen for this experiment were from one patient not considered to have bone metastases and from one who did have bone metastasis. The samples correspond to lanes 8 and 37 in Fig. 1
. The gels are shown in Fig. s2A of the online Data Supplement and exhibited four differences in the mass region
11.7 kDa when we compared the serum from the patient with bone metastases with the serum obtained from the patient who did not have bone metastases. These protein spots, plus a blank region at the edge of the gel, were excised and subjected to an in-gel trypsin digestion protocol. The resulting trypsin fragments from each spot were analyzed in a QStar Pulsar I with the MALDI I ion source to identify the unknown proteins, which were identified to be variants of SAA. The amino acid sequence of unprocessed SAA is shown in Fig. s2B of the online Data Supplement, with the variable amino acids boxed. A MALDI-TOF mass spectrum of peptides derived from spot 4 of the two-dimensional gel is shown in Fig. s2C of the online Data Supplement. A Mascot search using the peptide mass fingerprinting (PMF) data indicated that four peptides matched with peptides for SAA, giving a sequence coverage of 40%. Significant confidence (P <0.05) in the identity was indicated by a Mascot score as high as 80. Further confirmation of protein identification was provided by MS/MS fragmentation analysis, which provided a sequence tag for database searching. Identification of the m/z 2178 ion confirmed its identity as a SAA peptide with a Mascot score of 131 (Fig. s2D in the online Data Supplement). MS/MS analysis of m/z 2178 detected most b and y ions (Fig. s2D in the online Data Supplement). Combined with PMF data, all four spots from the two-dimensional gel yielded very high confidence hits with SAA variants.
western blot analysis and immunodepletion assays
A representative sample from each group was selected for Western blot analysis based on whether the SELDI spectrum was negative or positive for what we putatively identified as SAA proteins. For this analysis, the samples corresponded to lanes 8 and 37 of Fig. 1
(negative and positive for SAA, respectively). Western blot analysis with an antibody against SAA detected a single band at a molecular mass of
12 kDa in sample 37, whereas nothing was detected in serum from a patient not considered to have bone metastases (Fig. 3
). Only one band could be detected in serum from a patient with confirmed bone metastases by Western blot when using a 20% tricine gel. Although the antibody used recognizes all four isoforms of SAA, these isoforms could not be resolved by Western blot analysis.
|
To confirm that the peaks detected by SELDI-TOF MS (Fig. 1
) were indeed isoforms of SAA, we next analyzed serum from a representative sample of each cohort (lanes 7 and 35 in Fig. 1
) both before and after immunodepletion of SAA. Selection of samples was based on the same rationale as described above. The spectrum obtained for serum from the patient with bone metastases had the cluster of peaks at 11.7 kDa, with intensities ranging from
5 to 13 (Fig. 4
). After immunodepletion, the intensities for these peaks were all <1 (compare panels A and B in Fig. 4
). The spectra obtained with serum from the patient without bone metastases did not have the cluster of peaks at 11.7 kDa, either before or after immunodepletion of SAA (panels C and D in Fig. 4
). A protein cluster at
13.8 kDa was only slightly diminished after immunodepletion, indicating the specificity of the assay.
|
confirmation of increased saa concentrations in sera from patients with bone metastases
Taken together, the results of the two-dimensional gel electrophoresis, Western blotting, and immunodepletion experiments described above demonstrate that the 11.7-kDa cluster of peaks in the SELDI-TOF MS mass spectra correspond to SAA. To validate that the serum SAA protein concentrations correlated with peak intensities in the SELDI-TOF MS spectra, we performed an ELISA for SAA in all serum samples and plotted the values vs peak intensity obtained by SELDI-TOF MS for each patient. This plot showed that peak intensity correlated to SAA values (Fig. 5
).
|
correlation of saa and il-6 in serum of prostate cancer patients
SAA has been reported to be increased by IL-6 (39). IL-6 has also been reported to be increased in patients with recurrent prostate cancer (26). We therefore measured IL-6 in these same serum samples and compared the results to the SAA concentrations detected (Table 2
). The results (Fig. 6A
) indicated that the two of the patients with bone metastases (lanes 37 and 38) and the highest SAA concentrations also have the highest IL-6 concentrations in this cohort (Fig. 6B
). One patient considered not to have bone metastases (Fig. 6
, lane 19) also had a very high IL-6 concentration when compared with others in this cohort. This particular patient also had the highest SAA concentrations and was the only patient to have PSA >100 µg/L at the time of diagnosis in the cohort with no bone metastases. The probability of metastatic disease increases with increasing PSA concentrations. No specific trend could be readily observed between SAA and IL-6; some patients with very high SAA concentrations did not have high IL-6 concentrations (e.g., Fig. 6
, lane 35).
|
|
Serum alkaline phosphatase is a marker for osteoblastic lesions (40); however, in patients with confirmed bone metastases, it was increased to above-normal concentrations (>200 U/L) in only 44% (8 of 18) of the samples. Serum concentrations of alkaline phosphatase were generally higher in patients with positive bone scans (Fig. 6C
, lanes 2038), but there was no obvious trend relating these concentrations to SAA (Fig. 6A
), IL-6 (Fig. 6B
), or PSA (Fig. 6
, D and E; see Tables 1
and 2
).
svm
Using linear class boundaries and a leave-one-out cross-validation scheme with spectra obtained from IMAC surfaces, SVMs correctly classified 34 of 38 samples, for a classification accuracy of 89.5%. Two samples associated with a positive test for bone metastasis and two samples associated with a negative test were misclassified. This corresponds to a sensitivity and specificity of 89.5%. The analysis pipeline described was applied to data obtained from other ProteinChip surfaces, but the best classification accuracy was obtained with the IMAC data. The simultaneous use of all data from various ProteinChip surfaces did not improve classification accuracy. No obvious unique peaks for bone metastases samples could be manually detected with IMAC-Cu++ ProteinChip Arrays.
| Discussion |
|---|
|
|
|---|
SELDI-TOF MS is a high-throughput proteomic approach that has been successfully applied for proteomic profiling in a variety of clinical samples, including serum (45)(46)(47) and laser-capture-microdissected tissues (48)(49)(50). These profiling studies have shown the ability of SELDI-TOF MS to produce spectra that can be analyzed by various machine-learning algorithms to differentiate cancer cases from noncancer cases. Only one report exists on the application of SELDI-TOF MS profiling of clinical samples for metastatic disease (49). The authors of that report used serum from patients with hepatocellular carcinoma but did not go on to identify potential biomarkers associated with the protein profile for metastasis. In the current study, we first applied SELDI-TOF MS technology, using reversed-phase H4 ProteinChip Arrays to rapidly screen unfractionated serum samples obtained from prostate cancer patients both with and with out bone metastases. On identification of a cluster of proteins that predominated in the spectra obtained from patients with bone metastases we used preparative two-dimensional gel electrophoresis followed by in-gel trypsin digestion, peptide mapping, and MS/MS to identify this cluster of proteins as isoforms of SAA. Evidence that these proteins detected in the SELDI-TOF MS spectra were indeed correctly identified as SAA was obtained by immunodepletion assay. Validation that SAA concentrations were increased in the sera obtained from patients with bone metastases was confirmed by ELISA. Spectra obtained from a different ProteinChip Array (IMAC-Cu++) with the same samples were analyzed by SVM, and this analysis was able to distinguish prostate cancer patients with bone metastases from those without metastases with a sensitivity of 89.5%.
SAA is an acute-phase protein [for a review, see Ref. (51)]. Circulating concentrations of SAA are transiently increased as much as 1000-fold in response to inflammation (52). SAA is produced predominantly in the liver as well as in other epithelial tissue, including prostate tissue (53). It is unlikely that the increased SAA concentrations detected here were attributable to enhanced production in prostatic tissue. There is an inherent lack of specificity for SAA; it is generally increased in patients with a broad spectrum of primary and disseminated diseases (54)(55). SAA has received considerable attention for being a potential marker for neoplastic activity, however, because of its inverse correlation with patient survival in spite of the lack of elucidation of its role in metastatic processes [for a review see Ref. (51)]. SAA is induced by IL-6 (39), and although IL-6 is reported to be increased in patients with metastatic prostate cancer, we were unable to confirm this finding in our samples (see Fig. 6, A and B
), thereby negating a possible mechanistic role correlating IL-6 to increased serum concentrations of SAA in prostate cancer patients with bone metastases. In 1984, SAA was proposed for use in the monitoring of active prostate cancer that was not in remission (56). Since that initial publication, SAA has not been reported for monitoring of prostate cancer, probably because of its lack of specificity and the successful emergence of PSA as a biomarker for this disease.
Nevertheless, to our knowledge, the work presented here is the first to identify a biomarker for metastasis based on a combination of SELDI-TOF MS, two-dimensional gel electrophoresis, in-gel trypsin digestion, PMF, and MS/MS followed by validation by ELISA and immunodepletion assay, with the results obtained with the various approaches compared with results obtained with machine-learning algorithms. Results from these comparison are shown in Table 3
and indicate that SVM was the most sensitive approach for differentiating prostate cancer patients with bone metastases, with a sensitivity of 89.5% compared with 84% for ELISA and 58% for SELDI-TOF MS alone. Application of SVM involves the analysis of
270 peaks representing peptides or proteins per spectrum. Thus, SVM analysis uses multiple proteins/peptides rather than just the single SAA cluster of isoforms, as are used with manual inspection of SELDI-TOF MS spectra or ELISA. These data support that measurement of multiple proteins or a panel of proteins can enhance the detection rate [for a review, see Ref. (57)]. Thus, SELDI-TOF MS protein profiling may provide a diagnostic tool with potential clinical applications as well as aid in the discovery of biomarkers associated with various diseases.
|
| Acknowledgments |
|---|
| Footnotes |
|---|
| References |
|---|
|
|
|---|
The following articles in journals at HighWire Press have cited this article:
![]() |
J. S. Vermaat, I. van der Tweel, N. Mehra, S. Sleijfer, J. B. Haanen, J. M. Roodhart, J. Y. Engwegen, C. M. Korse, M. H. Langenberg, W. Kruit, et al. Two-protein signature of novel serological markers apolipoprotein-A2 and serum amyloid alpha predicts prognosis in patients with metastatic renal cell cancer and improves the currently used prognostic survival models Ann. Onc., December 18, 2009; (2009) mdp559v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. I. van Leuven, R. Hezemans, J. H. Levels, S. Snoek, P. C. Stokkers, G. K. Hovingh, J. J. P. Kastelein, E. S. Stroes, E. de Groot, and D. W. Hommes Enhanced atherogenesis and altered high density lipoprotein in patients with Crohn's disease J. Lipid Res., December 1, 2007; 48(12): 2640 - 2646. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Nedelkov, U. A. Kiernan, E. E. Niederkofler, K. A. Tubbs, and R. W. Nelson Population Proteomics: The Concept, Attributes, and Potential for Cancer Biomarker Research Mol. Cell. Proteomics, October 1, 2006; 5(10): 1811 - 1818. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. L. Hortin The MALDI-TOF Mass Spectrometric View of the Plasma Proteome and Peptidome Clin. Chem., July 1, 2006; 52(7): 1223 - 1237. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |