|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Molecular & Cellular Proteomics 7:163-169, 2008.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
We have developed a different protein microarray platform, dissociable antibody microarray (DAMA)1 staining (26). This technology combines the power of immunohistochemical staining and the parallel analysis of antibody microarrays. DAMA staining provides a new approach in the global analysis of protein expression and subcellular localization. In DAMA staining, targeted cells are grown on a coverslip in a culture dish or mounted on a coverslip and fixed and permeabilized with the standard protocol for traditional immunostaining. In the next step, instead of adding a primary antibody, an array of antibodies immobilized on a membrane is placed on top of the specimen. Pressure is applied to maintain close contact between the antibody array and cells. During incubation, antibodies dissociate from the array support and bind to their respective antigens in the cells without significant lateral diffusion. In this way, hundreds of antibodies are delivered to the targeted cells or tissues in a position-dependent manner. The bound antibodies are then detected either with enzyme- or fluorophore-conjugated secondary antibodies. The expression profiles of hundreds of proteins can thus be determined simultaneously. Furthermore when stained with fluorophore-conjugated secondary antibodies, subcellular localization profiles of hundreds of proteins can be obtained from a single staining by a fluorescence microscope equipped with a computer-controlled motorized stage.
Here we report the development of DAMA staining technology and its application in identifying potential biomarkers for breast cancer. Breast cancer is the most common form of cancer in women with
210,000 new cases diagnosed annually in the United States alone (27). This disease causes significant mortality, accounting for
40,000 deaths in the United States per year (28) and many more fatalities worldwide. Earlier detection and better treatment will improve prognosis and survival of the disease (27). Detection of new molecular biomarkers will not only be very useful for breast cancer diagnosis but also may identify novel and key targets for therapy. In this report, we compared the expression profiles of 312 proteins among 10 different normal breast and breast cancer cell lines using DAMA staining. We also developed a data analysis program to identify differentially expressed proteins and validated the prediction by Western analysis in the same cell lines.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
500 µm in diameter. The antibody list for Array-320 is shown in the supplemental materials. Antibodies were selected from Hypromatrix's collection. All antibodies have been characterized and demonstrated to bind their targets in various assays.
Cell Culture—
Ten different breast cell lines, MCF10A, MCF12A, Hs578Bst, MCF7, T-47D, ZR-75-1, MDA-MB-231, BT549, Hs578T, and MDA-MB-435S, were purchased from American Type Culture Collection (ATCC). Cells were maintained and propagated as recommended by the ATCC and were grown on 10-cm cell culture dishes until 90–95% confluence for DAMA staining.
The DAMA Staining Procedure—
Cells were fixed with formaldehyde solution, permeabilized with Triton X-100, and blocked with goat serum. A 320-antibody microarray (Array-320) was placed over the cells and was incubated at room temperature under the optimal pressure of
100g/cm2 as measured by a pressure detection sensor in an Economical Load and Force (ELF) single handle system (Tekscan, Inc.). Bound antibodies were detected by alkaline phosphatase-conjugated secondary antibodies (both goat anti-rabbit and goat anti-mouse) with 1-stepTM NBT/BCIP (nitro blue tetrazolium/5-bromo-4-chloro-3-indolyl phosphate) substrate (Pierce). Images of protein expression profiles were scanned by an HP Scanjet 4890 and processed by Photoshop for intensity integration.
Intensity Integration and Initial Data Analysis—
The ScanAlyze program developed for DNA microarray analysis was used for intensity integration of the DAMA staining images. One image (e.g. the reference cell line, such as MCF10A) was used as the channel 1 data, and the other image (e.g. other sample cell lines, such as MCF7) was used as the channel 2 data. The scanned DAMA staining images were first compressed into eight-bit grayscale images with a resolution of 600 d.p.i. and color-inverted by Photoshop. The preprocessed images were then combined into an RGB overlay image with one in red channel 1 and another in green channel 2. A grid was created with a circular mask defining the boundary for each spot. For the purpose of intensity integration, the location of the mask was refined for every spot. MRAT(i, j) values, the median ratio of intensity between channel 1 and channel 2 within the mask, were determined and exported to an Excel file. The subsequent analysis was done with a customized program (DAMAPEP). Briefly MRAT(i, j) was converted to log2(MRAT(i, j)), normalized for all 312 proteins by subtracting the median log2(MRAT(i, j)) value, and scaled by dividing the root mean square of all normalized log2(MRAT(i, j)) values. The median intensity of every spot and their corresponding median and mean background values, CH1I, CH1B, and CH1AB for the reference cell line R and CH2I, CH2B, and CH2AB for the sample cell line S, were also exported for the analysis by the DAMAPEP program.
Western Blotting Analysis—
Ten breast cell lines were grown on 10-cm cell culture dishes in corresponding medium to around 90–95% confluence. The whole cell lysates were extracted, concentrated, then quantified with the Bradford assay, and adjusted to the same values. Equal amounts of total protein were loaded in each lane (29).
Gel Quantification and Statistical Analysis—
For every protein, relative band densities of 10 cell lines were scanned and digitized by UN-SCAN-IT software (Silk Scientific, Orem, UT). These values were normalized against the highest value and averaged over all experiments. Relative expression levels of five proteins between normal and cancer cells were analyzed by using the Student's t test. For every protein, expression levels in three normal cells and in seven cancer cells were analyzed by a two-sample t test assuming equal variances. The degree of freedom is 9 (n1 = 3, n2 = 7). The one-tailed p values were used for statistical inference.
| RESULTS |
|---|
|
|
|---|
For this purpose, protein expression profiles of 10 different cell lines from human mammary glands were obtained from DAMA staining (Fig. 1). The 10 breast cell lines include three normal cell lines (MCF12A, MCF10A, and Hs578Bst), three estrogen receptor-positive carcinomas (T-47D, MCF7, and ZR-75-1), and four estrogen receptor-negative carcinomas (MDA-MB-231, BT549, Hs578T, and MDA-MB-435S). The cells were grown on a 10-cm cell culture dish to 90–95% confluence and fixed and permeabilized by a standard protocol. Array-320 was then used to deliver an array of 312 primary antibodies to those fixed cells. Bound antibodies were detected by using alkaline phosphatase-conjugated secondary antibodies. The resulting array of gray dots was imaged and subjected to intensity integration and data analysis to generate expression profiles for the 312 proteins. The experiments were repeated at least twice for each cell line. The representative images and the Pearson's correlation coefficients (Rr) of the replicated images for all 10 breast cell lines are summarized in Fig. 1. Most cell lines had reproducible DAMA staining images with the Pearson's correlation coefficients higher than 0.5.
|
Initial Data Analysis for the DAMA Staining Images with Program DAMAPEP—
A program, DAMAPEP (DAMA protein expression profiling), was developed to retrieve, normalize, and scale the data from the exported ScanAlyze files. The method for this data analysis is shown in Fig. 2. As the protein expression profiles for each cell line were determined at least twice, there are two independent DAMA staining images for each spot: two for the sample cell line (S1 and S2) and two for the reference cell line (R1 and R2). Therefore, four different sets of MRAT values, MRAT(1, 1), MRAT(1, 2), MRAT(2, 1) and MRAT(2, 2), corresponding to the intensity ratios of S1 to R1, S1 to R2, S2 to R1, and S2 to R2, respectively, were obtained for every spot(i, j) (Fig. 2a). The logarithms of those MRAT values were normalized and scaled. The average of those scaled log2(MRAT(i, j)) values, defined as SigS(i, j), represents the intensity change between the sample cell line (S) and the reference cell line (R) for the protein at spot(i, j) (Fig. 2c, Equation 1).
|
Comparison of the Protein Expression Profiles with Program DAMAPEP—
After calculating signal (SigS(i, j)) and background (BKS(i, j)) values for all spots, the program DAMAPEP identifies those proteins differentially expressed between normal and cancer cell lines. For this purpose, three different criteria were utilized to identify those candidate proteins from the DAMA staining images (Fig. 2b).
First a decision was made according to the average intensity change of the protein at spot(i, j) between all cancer cells and all normal cells. SigCan(i, j), defined as the average SigS(i, j) for all seven cancer cells (T47D, MCF7, Zr-75-1, MDA-MB-231, BT549, HS578T, and MDA-MB-435s), was calculated for every spot(i, j) (Fig. 2, b and c, Equation 3). This value represents the average intensity change of the protein at spot(i, j) in all cancer cell lines relative to the reference cell line MCF10A. Similarly SigNor(i, j), the average SigS(i, j) for two normal cells (MCF12A and HS578Bst), was also calculated for every spot(i, j) (Fig. 2c, Equation 4). SigNor(i, j) corresponds to the average intensity change of the protein at spot(i, j) in other normal cell lines relative to the intensity in the reference cell MCF10A. The Ratio(i, j) calculated (Fig. 2c, Equation 5) thus represents the average intensity change of the protein at spot(i, j) between cancer cells and normal cells relative to that of the normal cells (Fig. 2b). The higher the absolute value of Ratio(i, j), the larger the intensity difference. Therefore, DAMAPEP can predict a list of proteins with different expression levels between normal and cancer cells based on the selected cutoff value of Ratio(i, j).
When 2.0 was used as a cutoff value of Ratio(i, j), 54 proteins were predicted to be differentially expressed proteins between normal and cancer cells. However, 22 of the 54 proteins had lower spot intensity in at least one of the 10 cell lines. The contribution of those weak intensity spots to the corresponding value of Ratio(i, j) could be enlarged. To test the possibility, the expression levels of eight proteins, randomly selected from the 22 proteins, were compared among 10 cell lines by Western blotting analysis. All eight proteins showed low or undetectable expression levels (data not shown). Therefore, to decrease false positives caused by lowered spot intensity, DAMAPEP automatically compares the intensity of each spot with their background from the ScanAlyze exported files (CH1I versus CH1B and CH1AB; CH2I versus CH2B and CH2AB). The spots with intensities lower than their corresponding backgrounds in any one of the cell lines were excluded from the final prediction list (Fig. 2b).
The third criteria, Sigcan(i, j) (the averaged SigS(i, j) values for all seven cancer cell lines), could easily be biased by large numbers in one or two cell lines. To eliminate this bias, the SigS(i, j) value was compared with the corresponding background BKS(i, j) for all seven cancer cell lines. Only those proteins at spot(i, j) having more than four cell lines with higher SigS(i, j) value than the corresponding background value BKS(i, j) were included in the final prediction list (Fig. 2b).
Identifying the Differentially Expressed Proteins between Normal and Cancer Cell Lines by DAMAPEP—
As discussed above, the expression profiles of 312 proteins obtained by the DAMA staining for 10 breast cell lines were repeated at least twice and were analyzed by the programs ScanAlyze and DAMAPEP. When using a minimum cutoff value of 2.0 for the absolute value of Ratio(i, j), 10 of 312 proteins (leukemia inhibitory factor receptor alpha, Tyk2, Rb p130, SRF, c-Kit, Rb p107, RAIDD, Mos proto-oncogene, I
B-β, and IL-2Rβ) were predicted by DAMAPEP to have different expression levels in cancer cells versus normal cells (Fig. 3a). The first seven proteins were predicted to have higher expression levels in cancer cells, and the latter three were predicted to have lower expression levels. Their corresponding dots in the DAMA images for all 10 cell lines demonstrated the consistency between the DAMAPEP prediction and the original DAMA staining data especially for those confirmed proteins (Fig. 3b and supplemental Fig. S1). When the cutoff value for Ratio(i, j) was decreased from 2.0 to 1.5 (representing a 50% intensity change between normal and cancer cells) six extra proteins were added to the prediction list. Among them, ErbB3, c-Raf, and Rad52 were predicted to have higher expression in cancer cells, and Fas (TNFRSF6)-associated death domain protein, FLICE inhibitory protein (short/long), and IL-1R were predicted to have lower expression. However, none of these six proteins was confirmed to have the predicted difference by Western blotting analysis (data not shown). Therefore, 2.0 was used as the cutoff value for Ratio(i, j) in the DAMAPEP prediction.
|
B-β were similar in both normal and cancer cell lines (data not shown). The expression of RAIDD was almost undetectable in the three normal cell lines. The expression of Rb p107, SRF, and Tyk2 in normal cells was between 15 and 40% of their expression in cancer cells. The expression of Rb p130 in normal cells was about 20% of its expression in cancer cells (Fig. 4b). Mean values of the normalized expression of RAIDD, Rb p107, Rb p130, SRF, and Tyk2 in three normal cell lines are 5 versus 61% (p = 0.0032), 26 versus 81% (p = 0.0041), 18 versus 69% (p = 0.0055), 25 versus 73% (p = 0.0033), and 30% (p = 0.0012), respectively (Fig. 4c). Those data confirm that RAIDD, Rb p107, Rb p130, SRF, and Tyk2 are indeed overexpressed in the tested cancer cells.
|
| DISCUSSION |
|---|
|
|
|---|
The observed higher expression levels of two proteins in breast cancer cells are consistent with published reports. The increased expression of Rb2/p130 gene was observed in the majority of breast cancers (30). Up-regulated and highly active SRF was observed in squamous epithelial tumor cells (31). However, we are not aware of any studies addressing the expression changes of RAIDD, Rb p107, and Tyk2 in breast cancer.
Among the seven proteins predicted with higher expression level in breast cancer cells, five were confirmed by Western blot analysis. On the other hand, none of the three proteins predicted to have lower expression levels were confirmed by Western blotting. These errors could be caused by the quality of the DAMA staining images. For example, images of MCF7, ZR-75-1, and BT549 have a higher density in the lower part of the array (Fig. 1). This uneven intensity distribution could be caused by uneven pressure applied during primary antibody staining. In addition, the lower part of the array corresponds to the center of the culture dish, which could have higher cell density and therefore might result in the observed uneven distribution.
Two cell lines with lower Pearson's correlation coefficients, Hs578Bst and ZR-75-1, could also affect the prediction accuracy. The overlap coefficients for these two cell lines, a modified Pearson's correlation coefficient where the average gray value of each image was not subtracted from the original gray value of each pixel, were also lower. Those data suggested that superimposition of those images was not perfect possibly due to their lower spot intensities. The situation was improved by increasing the spot diameters for intensity integration by ScanAlyze. Therefore, the prediction accuracy could be improved by better DAMA staining images after further optimizing the experimental conditions.
Another possible reason for the inaccurate prediction of decreased expression levels is systematic bias in data analysis especially in the normalization procedure. As there are only 312 proteins in the dataset and most proteins are highly regulated, the assumption that those 312 proteins should have the same expression levels among different samples could generate systematic bias during data normalization. This systematic bias could be decreased by including additional antibodies on the array. In addition, deploying multiple replicate spots for each antibody on the array could also increase data accuracy, decreasing experimental error and increasing prediction accuracy.
The alkaline phosphatase-based detection with its narrow detection range of signal intensity could be another reason for the inaccuracy of the method. The intensity difference between the strongest spot and the weakest spot of the staining images as determined by the protein expression level and antibody affinity could be hundreds of thousands-fold. The weak spot may obscure the difference. For example, both ER-positive and ER-negative cell lines were utilized. However, the DAMA staining for antibody against ER at spot F5 in all cell lines showed weak signal. This could be improved by using a fluorescent conjugated secondary antibody and a fluorescence scanner. We tried several commercial microarray scanners but failed to obtain adequate signals. A more sensitive scanner may be needed for our purpose. In addition, the homogeneity of cell samples and the percentage of cell confluence could also affect the accuracy of protein expression profiling. For this reason, cells were grown to 90–95% confluence for consistency. Under current experimental conditions, 10 of 320 proteins were identified by the DAMA staining as differentially expressed proteins, and their expression patterns were consistent with those obtained by Western blot analysis, confirming the accuracy and broad application of the DAMA staining technology.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
Published, MCP Papers in Press, October 13, 2007, DOI 10.1074/mcp.M700115-MCP200
1 The abbreviations used are: DAMA, dissociable antibody microarray; DAMAPEP, DAMA protein expression profiling; ER, estrogen receptor; RIP, receptor-interacting protein; RAIDD, RIP-associated ICH-1/CED-3-homologus protein with a death domain; SRF, serum response factor. ![]()
2 X. C. Song, G. Fu, X. Yang, Z. Jiang, Y. Wang, G. W. Zhou, unpublished results. ![]()
* This work was supported by National Institutes of Health Grant CA110047 and a grant from the Susan G. Komon Breast Cancer 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. ![]()
S The on-line version of this article (available at http://www.mcponline.org) contains supplemental material. ![]()
Both authors contributed equally this work. ![]()
** To whom correspondence should be addressed. Tel.: 225-578-6733; Fax: 225-578-8011; E-mail: zhouw{at}lsu.edu
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
O. Simma, E. Zebedin, N. Neugebauer, C. Schellack, A. Pilz, S. Chang-Rodriguez, K. Lingnau, E. Weisz, E. M. Putz, W. F. Pickl, et al. Identification of an Indispensable Role for Tyrosine Kinase 2 in CTL-Mediated Tumor Surveillance Cancer Res., January 1, 2009; 69(1): 203 - 211. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Gakovic, J. Ragimbeau, V. Francois, S. N. Constantinescu, and S. Pellegrini The Stat3-activating Tyk2 V678F Mutant Does Not Up-regulate Signaling through the Type I Interferon Receptor but Confers Ligand Hypersensitivity to a Homodimeric Receptor J. Biol. Chem., July 4, 2008; 283(27): 18522 - 18529. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| All ASBMB Journals | Journal of Biological Chemistry |
| Journal of Lipid Research | ASBMB Today |