This article has Open Peer Review reports available.
Inferring predominant pathways in cellular models of breast cancer using limited sample proteomic profiling
© Kulkarni et al; licensee BioMed Central Ltd. 2010
Received: 10 December 2009
Accepted: 15 June 2010
Published: 15 June 2010
Molecularly targeted drugs inhibit aberrant signaling within oncogenic pathways. Identifying the predominant pathways at work within a tumor is a key step towards tailoring therapies to the patient. Clinical samples pose significant challenges for proteomic profiling, an attractive approach for identifying predominant pathways. The objective of this study was to determine if information obtained from a limited sample (i.e., a single gel replicate) can provide insight into the predominant pathways in two well-characterized breast cancer models.
A comparative proteomic analysis of total cell lysates was obtained from two cellular models of breast cancer, BT474 (HER2+/ER+) and SKBR3 (HER2+/ER-), using two-dimensional electrophoresis and MALDI-TOF mass spectrometry. Protein interaction networks and canonical pathways were extracted from the Ingenuity Pathway Knowledgebase (IPK) based on association with the observed pattern of differentially expressed proteins.
Of the 304 spots that were picked, 167 protein spots were identified. A threshold of 1.5-fold was used to select 62 proteins used in the analysis. IPK analysis suggested that metabolic pathways were highly associated with protein expression in SKBR3 cells while cell motility pathways were highly associated with BT474 cells. Inferred protein networks were confirmed by observing an up-regulation of IGF-1R and profilin in BT474 and up-regulation of Ras and enolase in SKBR3 using western blot.
When interpreted in the context of prior information, our results suggest that the overall patterns of differential protein expression obtained from limited samples can still aid in clinical decision making by providing an estimate of the predominant pathways that underpin cellular phenotype.
Breast cancers are clinically heterogeneous . Particular molecules have been identified that are associated with clinical prognosis. For instance, 20% to 25% of breast cancers are associated with the overexpression of HER2, and its presence is associated with poor prognosis [2, 3]. In addition, increased ER/PR expression has been identified in 70% of breast cancer patients. These biomarkers have motivated a shift away from "one size fits all" approach of treating breast cancer to developing therapies that target specific molecules. In particular, tamoxifen, a selective ER modulator (SERM), improves survival for patients with ER/PR positive tumors [4, 5]. Trastuzumab, a human monoclonal antibody was developed to bind the HER2 receptor and block its activity . However, de novo or acquired resistance to tamoxifen [7–9] and trastuzumab  has been an emerging problem.
Molecularly targeted drugs, like trastuzumab and tamoxifen, are designed to block aberrant signaling within oncogenic pathways. Cell signaling pathways direct the flow of information (i.e. flux) from an extracellular stimulus to the corresponding cellular response (e.g., cellular proliferation, contact inhibition, or cellular death). By analogy with metabolic control analysis, control of flow of information is distributed among all the steps in a network [11–13]. This implies that an increase in one protein does not necessarily correspond to an increased pathway flux. Conversely, the decrease in expression of one protein via therapeutic modification does not necessarily lead to a decrease in pathway flux. Conceptually, this leads to the hypothesis that combining gene expression measurements over group of genes that fall within common pathways will be more effective means of marker identification. In fact, recently it was shown that breast cancer genes that do not exhibit a change in their expression profile still play a central role interconnecting deregulated genes in a protein network . The observation that onset and progression of many diseases arises from the interactions of a number of interconnected genes has shifted the drug discovery perspective from a molecule-centric to a network/pathway-centric approach . Proteomics provides an attractive platform for interrogating pathway flux as measuring actual protein levels instead of measuring proxy mRNA levels maybe more informative in spite of added experimental complexity .
One of the most commonly used techniques for proteomic profiling is 2DE based protein separation in combination with mass spectrometry based identification. Using this approach, in addition to analyzing proteins in the blood, tumor tissues are being examined to yield insights about molecular pathways that are altered in cancer progression. While 2-DE based high-throughput proteomic data reveal proteins that are differentially regulated, different sources of biological and analytical variations affect the statistical importance of these results . These can be addressed using an experimental design that incorporates several technical and biological replicates to account for variations at two levels, within gels and within samples, respectively [18–20]. This puts a demand on the sample amount and composition as proteomic analysis of breast cancer biopsies is complicated due to heterogeneity of cellular phenotypes contained in the sample . While laser capture microdissection (LCM) can provide a relatively homogeneous sample by concentrating on the cell type of interest, generating enough sample for a conventional proteomic study is laborious with a minimum of 100,000 cells and a dissection time in tens of hours required for one 2D-PAGE . Given the desire to aid in clinical decision-making, the ability to obtain sufficient clinical sample presents a significant challenge. In the case of sample that is insufficient to carry out a proteomic study with multiple replicates, like in the case of an early stage breast tumor, does the information obtained from a single gel replicate still provide insight into the predominant cell signaling pathway at work in a cell?
Thus the objective of this study was to identify predominant pathways and protein interaction networks in two breast cancer phenotypes using prior information. Given the conservation of genetic information and marker expression between tumors and their corresponding cell lines [23–25], we have used well characterized model systems in our study. The central hubs of protein interaction networks obtained using prior information were validated to establish confidence in the protein expression patterns.
Cell culture and reagents
The human breast cancer cell lines (BT-474 and SK-BR-3) were kindly provided by Dr. Jia Luo (Health Sciences Center; West Virginia University, WV). Cells were grown in 75-cm2 plastic tissue culture flasks (Costar Corning; Corning, NY) in a humidified incubator at 37°C and 5% (v/v) CO2. The BT-474 cells were routinely maintained in Rosewell Park Memorial Institute (RPMI) 1640 medium (Mediatech, Inc., Herndon, VA) supplemented with 10% (v/v) heat inactivated fetal bovine serum (FBS) (Hyclone, Inc., Logan, UT), 0.3% (w/v) L-glutamine, 1% (v/v) penicillin/streptomycin (BioWhittaker, Walkersville, MD) and 10 ng/mL insulin (Sigma, St Louis, MO). SK-BR-3 cells were maintained in Improved Modified Eagle Medium (IMEM) Zn2+ option (Invitrogen) containing 4 mM L-glutamine, 2 ml/L L-proline, 50 μg/mL gentamicin sulfate supplemented with 10% FBS (Hyclone) and 1% penicillin/streptomycin (BioWhittaker). Cells were passaged at 1:5 dilution with fresh medium every 5 days.
Preparation of cell lines for 2-DE
Cells were grown to approximately 80% confluence. Growth medium was removed from dishes and cells were washed twice with 10 mL Phosphate Buffered Saline (PBS) to remove dead cells as many extracellular proteins as possible. Cells were made non-adherent by incubating the flasks at 37°C for 10 min in the presence of trypsin (BioWhittaker). Trypsin was neutralized by the addition of FBS. Cells were then washed twice with PBS and harvested at 1,200 rpm at 4°C for 10 min. Sufficient precaution was taken to get rid of PBS to eliminate salts that could possibly interfere with the 2DE. Cells were incubated in lysis buffer (7M Urea, 2M thiourea, 2% (w/v) CHAPS) for 30 min on ice and sonicated five times in an ultrasonic water bath, where each sonication was performed for 10 s followed by 10 s cooling interval on ice. Cell debris were pelleted by centrifugation at 14,000 rpm for 40 min at 4°C. The supernatant was aliquoted in fresh tubes and stored at -80°C. The protein concentration was determined using BCA protein assay kit (Pierce).
For each cell line, 500 μg of cell lysate was mixed with rehydration buffer (7M urea, 2M thiourea, 2% CHAPS, 1% DTT, 2% IPG buffer, 0.002% bromophenol blue) and incubated for 1 h at room temperature prior to rehydration on Immobilized pH Gradient (IPG) strips pH 3-10 NL, 24 cm, (GE Healthcare, Uppsala, Sweden) for 12 h at 25°C. Isoelectric focusing was done using Ettan IPGphor apparatus (Amersham Biosciences) for a total of 90 kVh at 50 μA per strip at 20°C. Thereafter, IPG strips were equilibrated in 75 mM Tris-HCl pH 8.8, 6M urea, 30% (v/v) glycerol, 2% (w/v) SDS, 0.002% (w/v) bromophenol blue and 1% (w/v) DTT for 30 min. A second equilibration step was done for another 30 min by replacing the DTT with 2.5% iodoacetamide. Equilibrated strips were transferred onto 24 cm 12% uniform precast SDS-polyacrylamide gels (Jule, Inc., Milford, CT) poured between non-fluorescent glass plates. IPG strips were sealed with 0.5% (w/v) low melting point agarose in SDS running buffer containing bromophenol blue. Gels were run in Ettan DALTsix Larger Vertical System (Amersham Biosciences) at 30mA per gel at room temperature, until the dye front had run off the bottom of the gels.
Gels were fixed in 10% (v/v) methanol, 7% (v/v) acetic acid overnight, washed in 18 MΩ water, and stained overnight with SYPRORuby dye (Bio-Rad). Excess dye was removed by washing twice with 18 MΩ water in a dark room. Gels were imaged using the Typhoon 9400 scanner (Amersham Biosciences) at 200 μm resolution with a 488nm laser with 610nm band pass filter at normal sensitivity under fluorescence acquisition mode. Data were saved in .gel format using ImageQuant software (Amersham Biosciences). The 2-DE results are representative of three biological replicates.
The images were analyzed using SameSpots software from Nonlinear Dynamics. Saturated and damaged areas of the gels were ignored in the analysis by selecting a region of interest (ROI). The images were warped using automatic and manual vectors to a reference image that was automatically selected based on the gel containing the most spots. Normalized spot volumes were generated from the optical densities for each individual spot to the ratio of the total spot volume in each gel. 304 differentially expressed protein spots were chosen for further analysis.
The gel spots of interest were excised using an Ettan Spot Picker (Amersham Biosciences) fitted with a 1.5-mm spot picker head. Briefly, specified excised spots were reduced in DTT (10 mM, 60°C, 10 min) and alkylated with iodoacetamide (100 mM, room temperature, 45 min) in a dark room. The gel pieces were dehydrated in acetonitrile for 10 min. Then the gel pieces were vacuum dried and rehydrated with 10 μL of digestion buffer (10 ng/μL of trypsin (Promega; Madison, WI) in 25 mM NH4HCO3) and covered with 10 μL of NH4HCO3. The samples were incubated for 16 h at 37°C to allow for complete digestion. Peptides were extracted from gel plugs by sonication in 2.5 μL 5% formic acid.
MALDI-TOF MS analysis
MALDI-TOF-MS system model Micromass MALDI-R (Waters®) was used to obtain the peptide mass fragment spectra as recommended by the manufacturer. Protein digest solutions were mixed at a 1:1 ratio with the MALDI matrix α-cyano-4-hydroxycinnamic acid (CHCA) (Sigma-Aldrich Fluka; St. Louis, MO). 1 μL of tryptic peptide sample was applied to the MALDI plate and allowed to dry. The MALDI-TOF MS was operated in the positive ion delayed extraction reflector mode for highest resolution and mass accuracy. Peptides were ionized/desorbed with a 337-nm laser and spectra were acquired at 15 kV accelerating potential with optimized parameters. The close external calibration method employing a mixture of standard peptides (Applied Biosystems) provided mass accuracy of 25-50 ppm. Internal calibration was performed with the monoisotopic peak of adrenocorticotropic hormone (ACTH) (18-39) peptide (m/z: 2465.1989). Mass spectral analysis for each sample was based on the average of 300 laser shots. Peptide masses were measured from m/z: 800 to 3,000. The peak lists containing the m/z ratio and corresponding intensity values were exported to Microsoft Excel for further processing.
Protein identification using peptide mass fingerprinting (PMF)
Peptide mass fingerprints for each of the 304 proteins were entered in an Excel spreadsheet along side each other. To optimize the database searching, the list of peptide mass peaks from the spectrum of each sample was processed and background peaks that were observed in greater than 10% of the PMF's were eliminated to improve the efficiency of database searching . MASCOT http://www.matrixscience.com, Aldente (ExPASy) and MS-Fit (Protein Prospector; University of California, San Francisco) were each used to query the UniProtKB/Swiss-Prot human database with the corresponding monoisotopic peptide mass fingerprints with the following settings: peptide mass tolerance of 50 ppm, one missed cleavage site, one fixed modification of carboxymethyl cysteine, one variable modification of methionine oxidation, and no restrictions on protein molecular mass or isoelectric point. The protein identities reported were ranked high in at least two of the three algorithms used.
Ingenuity pathway Analysis
Differentially regulated proteins identified by 2DE and PMF were further analyzed using Ingenuity Pathway Analysis (IPA; Ingenuity Systems, Mountain View, CA; http://www.ingenuity.com). IPA was used to interpret the differentially expressed proteins in terms of an interaction network and predominant canonical pathways. The Ingenuity Pathways Knowledge Base (IKB) is a regularly updated curated database that consists of interactions between different proteins culled from scientific literature. IPA uses this database to construct protein interaction clusters that involve direct and indirect interactions, physical binding interactions, enzyme-substrate relationships, and cis-trans relationships in transcriptional control. The networks are displayed graphically as nodes (proteins) and edges (the biological relationship between the proteins).
A protein interaction network was generated as follows. A dataset containing the upregulated proteins, called the focus proteins, for a particular cell line was uploaded into the IPA. These focus proteins were overlaid onto a global molecular network developed from the information in the IKB. Networks of these focus proteins were then algorithmically generated by including as many focus proteins as possible and other non-focus proteins from the IKB that are needed to generate the network based on connectivity.
2 × 2 contingency table used for testing the significance of gene/protein enrichment in all IKB pathway annotations.
Genes associated with pathway
n - k
Genes not associated with pathway
K - k
(N - n) - (K - k)
N - n
N - K
For western blot analysis, 10-30 μg of total cell lysate was separated by SDS-PAGE using a 12% Tris polyacrylamide gel with a 4% stacking gel at 75 V for 4 h. Proteins were transferred onto Bio Trace PVDF membrane (PALL Life Sciences; Pensacola, FL) at 42 V for 1.5 h. Blots were washed in Tris Buffered Saline (TBS) for 5 min at room temperature, blocked for 1 h in TBS + 0.1% Tween 20 (TBS/T) plus 5% dry milk at room temperature and then washed three times in TBS/T. Blots were incubated overnight at 4°C with primary antibodies specific for IGF-1R (sc-9038), α-Enolase (sc-100812), GAPDH (sc-25778) (all from Santa Cruz Santa Cruz, CA), Ras (BD Biosciences, 610001) and Profilin (Millipore, AB3891) in TBS/T plus 5% dry milk. The next day, blots were washed three times in TBS/T, incubated for 1 h at room temperature with anti-biotin (Cell Signaling Technology, Inc., Danvers, MA, 7727) and either a goat anti-mouse IgG-horseradish peroxidase (HRP) (BD BioSciences, 554002) or a goat anti-rabbit IgG-HRP (Sigma-Aldrich, A0545). Finally, the blots were washed three times in TBS/T, developed using LumiGLO reagent (Cell Signaling Technology, Inc., Danvers, MA, 7003) and bands were visualized on KODAK Biomax light film (Fisher Scientific). Densitometric analysis was performed using ImageJ software (National Institute of Health) and protein levels were normalized to GAPDH protein levels for each sample. Given the uncertainty in estimating the level of expression in both cell lines, an empirical Bayesian approach was used to establish the level of confidence associated with differential expression between BT474 and SKBR-3 given the available data . Levels of expression were log-transformed to minimize potential bias in estimating the expression ratio, R, as follows: log10 (XBT474) = R + log10 (XSKBR3). A Markov chain Monte Carlo algorithm was used to estimate the posterior distribution in the differential expression coefficient. An initial unbiased gaussian prior distribution was used to propose new steps in the Markov chain. The prior distribution was scaled to achieve an acceptance fraction of 0.4. The Gelman-Rubin potential scale reduction factor was used to estimate convergence of three independent Markov chains to the posterior distribution . Posterior estimates of the expression ratio were obtained from the tails of the three independent chains following convergence.
Identification of differentially expressed proteins
2DE and Image Analysis
Peptide mass fingerprinting
Of the 866 spots identified to be common between both the gels using Progenesis SameSpots, spots that were not differentially expressed or detected in streaks were excluded from further analysis. 304 well resolved protein spots that were differentially expressed by at least 1.1-fold were excised from the gel, subjected to trypsin digestion and MALDI-TOF analysis. Using the resulting peptide mass fingerprints, 167 proteins were identified using Mascot as a primary database search algorithm, out of which 135 were unique. A subset of 62 identified proteins that were differentially regulated by at least 1.5-fold is shown in Additional file 1, Table II. More than 95% of the spots had sequence coverage exceeding 25%. In a majority of the cases, the identified proteins were the only candidate with significant score leading to their unambiguous identification. However, in some cases two different algorithms, MS-Fit and Aldente, were used to provide support for the protein identity from the PMF in order to minimize false positive identifications. Agreement between the apparent Mr and pI observed from 2-DE gels and the theoretical values of the identified proteins provided additional support for positive identification. A small subset of the identified proteins exists in multiple forms, as they were associated with multiple protein spots. This may be due to either protein degradation by proteases and formation of protein fragments or post-translational modifications  such as phosphorylation and carbamylation which is common with proteins in urea buffer and/or different isotypes. A list of all the identified proteins with their molecular weight, isoelectric point, rank, score and sequence coverage using the three algorithms is provided (Additional file 2, Table S1).
Protein distribution according to function, subcellular localization and fold-change
Pathway analysis and network generation using IPA
Validation of results by Western Blotting Analysis
Findings from the western blot were in accordance with the results from proteomic analysis. IGF1R and profilin were up-regulated in BT474 and Ras and alpha-enolase were up-regulated in SKBR3. The change in density-fold from 2DE to western blots can be attributed to the semi-quantitative nature of immunoblotting as well as the differences in the dynamic ranges of photographic quantification of chemiluminiscent and fluorescent methods of detection for immunoblotting and 2-DE respectively.
Discussion and Conclusion
It is the specific biological question that shapes the design of experiments . For example, a common approach is to detect changes in expression in each spot, considered individually, that is consistently above a threshold determined by the system's experimental noise. Implicitly this work reflected a desire to strike an optimal balance between the amount of data required and the ability to infer, with predictive potential, differentially activated pathways in two cancer models. Duncan and Hunsucker have used an engineering term - fitness-for-purpose - to characterize experimental design in proteomics where different constraints; such as limits on biological samples, effective use of resources, and how the information will be used; shape the experimental design . To strike a meaningful balance between these constraints, we hypothesized that a single gel replicate was sufficient to infer, with predictive potential, cell signaling pathways and protein networks that are differentially regulated between two breast cancer models. We demonstrated the predictive potential by validating the inferred differentially regulated cell signaling pathways using previously reported gene expression data  and by validating the inferred protein networks by western blot. Using previously reported gene expression data, we have attempted to minimize potential bias in our results introduced by our protocols (i.e., subtle differences in tissue culture or proteomics workflow). In cases where there is less prior information, such as a proteomic analysis of primary tissue, independent analysis of other gel replicates could help establish the confidence in the inferred differentially regulated cell signaling pathways. Weitzel et al.  recently described such an approach, where they used an additional gel replicate to confirm the upregulated pathways.
The cellular origins and predominant signaling pathways within these two cell lines are quite different. BT474 is derived from a solid invasive ductal carcinoma in the breast  while the SKBR3 cell line is derived from pleural effusion adenocarcinoma . The aggressiveness of these cell lines are different as BT474 is ER+/PR+ with a high in vitro invasion capability whereas SKBR3 is ER-/PR- with a low in vitro invasion capability. The aggressiveness of BT474 is supported by finding that the regulation of actin based motility by Rho and actin cytoskeleton signaling pathways were enhanced. In contrast, metabolic pathways like amino acid biosynthesis and glycolysis/gluconeogenesis were more pronounced in SKBR3 cell line.
To compare our results against gene expression data for the BT474 and SKBR3 cell lines, we used the mRNA expression data from the study of a collection of breast cancer cell lines. The gene expression data was reported as a matrix of probe sets by cell lines in which value is the calculated log abundance of each probe set gene for each cell line. Gene expression values were centered by subtracting the mean value of each probe set across the cell line from each measured value. To calculate the fold-difference between the two cell lines, these log abundance gene expression values were subtracted from each other. All the up-regulated genes for BT474 (8663 genes) and SKBR3 (8237 genes) were uploaded into the IPA. As shown in Additional File 4, Table S2, the top 5 associated network functions for BT474 had an identical significance score of 27 and were as diverse as cancer, skeletal disorder and dermatological diseases. Similar analysis for SKBR3 gene expression showed the top 5 functions had an identical significance score of 25 and were as varied as embryonic development, hematological system development, and lipid metabolism. The canonical pathway analysis applied to all of the differentially expressed genes revealed that none of the canonical pathways for either cell lines were significant (Additional file 5, Figure S1). This is because the multiple-testing criteria raise the threshold for significance such that none of the embedded gene expression patterns within the dataset provide a sufficiently strong signal to surpass this increased threshold. The threshold for gene expression data was increased by setting the cut-off value to be the same as protein expression data which was 1.5-fold. The resulting analysis for BT474 dataset which consisted of 506 genes showed 4 out of top 5 associated network functions to be related to cancer with scores better than 13 as shown in Additional File 6, Table S3. The most significant pathway for this analysis was IGF-1 signaling pathway with p-value < 1.2 × 10-4 (Figure 3B). For 304 genes which were up-regulated by a factor of 1.5-fold in SKBR3, 4 out of 5 associated network functions were related to cell death and cancer and had scores better than 11. The most significant pathway was urea cycle and metabolism of amino groups with a p-value < 6.45 × 10-4 (Figure 3B).
The group of pathways in SKBR3 cell line associated with our protein expression data was in agreement with the group of pathways associated with gene expression data. In BT474 cell line, though the top pathways for protein expression data were associated with cell motility and the pathway for gene expression data was associated with IGF-1 signaling, they are both associated with proliferation  and resistance to apoptosis [43, 44] in a broader sense. The number of focus molecules for the protein datasets involved in the top two networks for BT474 and SKBR3 were 25 and 29 which was 89.2% and 85.2% respectively of the total dataset; similar number for the gene datasets were 8.3% and 10.8% for BT474 and SKBR3 respectively, suggests that proteomics provides greater information per observation relative to gene expression.
In summary, experimental designs that consider each protein individually place a high burden on clinical samples. In cases where a sample is limited, a single replicate is unable to establish whether a single protein can be used as a biomarker. Our results do suggest that in a non-ideal case scenario, the overall pattern of differential protein expression can still be used, in conjunction with prior information, to infer pathways that underpin differences in cell phenotype. This information may prove helpful in tailoring therapies to the patient.
This work was supported by the PhRMA Foundation, National Cancer Institute (NCI) R15CA123123, and the National Institute of Allergy and Infectious Disease (NIAID) R56AI076221. The content is solely the responsibility of the authors and does not represent the official views of the NCI, the NIAID, or the National Institutes of Health. We thank Dr. Jia Luo for kindly providing the breast cancer cell lines.
- Henderson IC, Patek AJ: The relationship between prognostic and predictive factors in the management of breast cancer. Breast Cancer Res Treat. 1998, 52: 261-288. 10.1023/A:1006141703224.View ArticlePubMedGoogle Scholar
- Slamon DJ, Clark GM, Wong SG, Levin WJ, Ullrich A, McGuire WL: Human breast cancer: Correlation of relapse and survival with amplification of the HER-2/neu oncogene. Science. 1987, 235 (4785): 177-182. 10.1126/science.3798106.View ArticlePubMedGoogle Scholar
- Winston JS, Ramanaryanan J, Levine E: HER-2/neu evaluation in breast cancer are we there yet?. Am J Clin Pathol. 2004, 121 (suppl 3): 33-49.Google Scholar
- Osborne CK: Tamoxifen in the treatment of breast cancer. N Engl J Med. 1998, 339: 1609-1618. 10.1056/NEJM199811263392207.View ArticlePubMedGoogle Scholar
- Osborne CK, Zhao H, Fuqua SA: Selective estrogen receptor modulators: Structure, function, and clinical use. J Clin Oncol. 2000, 18: 3172-3186.PubMedGoogle Scholar
- Carter P, Presta L, Gorman CM, Ridgway JB, Henner D, Wong WL, Rowland AM, Kotts C, Carver ME, Shepard HM: Humanization of an anti-p185HER2 antibody for human cancer therapy. Proc Natl Acad Sci USA. 1992, 89 (10): 4285-4289. 10.1073/pnas.89.10.4285.View ArticlePubMedPubMed CentralGoogle Scholar
- Clarke R, Liu MC, Bouker KB, Zhiping G, Lee RY, Zhu Y, Skaar TC, Gomez B, O'Brien K, Wang Y, Hilakivi-Clarke LA: Antiestrogen resistance in breast cancer and the role of estrogen receptor signaling. Oncogene. 2003, 22: 7316-7339. 10.1038/sj.onc.1206937.View ArticlePubMedGoogle Scholar
- Osborne CK, Schiff R: Growth factor receptor cross-talk with estrogen receptor as a mechanism for tamoxifen resistance in breast cancer. Breast. 2003, 12: 362-367. 10.1016/S0960-9776(03)00137-1.View ArticlePubMedGoogle Scholar
- Nicholson RI, McClelland RA, Gee JM, Manning DL, Cannon P, Robertson JFR, Ellis IO, Blmaey RW: Epidermal growth factor receptor expression in breast cancer: association with response to endocrine therapy. Breast Cancer Res Treat. 1994, 29: 117-125. 10.1007/BF00666187.View ArticlePubMedGoogle Scholar
- Kute T, Lack CM, Willingham M, Bishwokama B, Williams H, Barrett K, Mitchell T, Vaughn JP: Development of herceptin resistance in breast cancer cells. Cytometry A. 2004, 57 (2): 86-93. 10.1002/cyto.a.10095.View ArticlePubMedGoogle Scholar
- Bruggeman FJ, Snoep JL, Westerhoff HV: Control, responses and modularity of cellular regulatory networks: acontrol analysis perspective. IET Syst Biol. 2008, 2 (6): 397-410. 10.1049/iet-syb:20070065.View ArticlePubMedGoogle Scholar
- Kholodenko BN, Hoek JB, Westerhoff HV, Brown GC: Quantification of information transfer via cellular signal transduction pathways. FEBS Lett. 1997, 414 (2): 430-434. 10.1016/S0014-5793(97)01018-1.View ArticlePubMedGoogle Scholar
- Westerhoff HV, Kell DB: What BioTechnologists knew all along...?. J Theoret Biol. 1996, 182: 411-420. 10.1006/jtbi.1996.0181.View ArticleGoogle Scholar
- Chuang HY, Lee E, Liu YT, Lee D, Ideker T: Network-based classification of breast cancer metastasis. Mol Syst Biol. 2007, 3: 140-10.1038/msb4100180.View ArticlePubMedPubMed CentralGoogle Scholar
- Russell RB, Aloy P: Targeting and tinkering with interaction networks. Nat Chem Biol. 2008, 4 (11): 666-673. 10.1038/nchembio.119.View ArticlePubMedGoogle Scholar
- Cox J, Mann M: Is proteomics the new genomics?. Cell. 2007, 130: 395-398. 10.1016/j.cell.2007.07.032.View ArticlePubMedGoogle Scholar
- Karp NA, Lilley KS: Design and Analysis Issues in Quantitative Proteomics Studies. Proteomics. 2007, 7 (1): 42-50. 10.1002/pmic.200700683.View ArticlePubMedGoogle Scholar
- Chich JF, David O, Villers F, Schaeffer B, Lutomski D, Huet S: Statistics for proteomics: experimental design and 2-DE differential analysis. J Chromatogr B. 2007, 849: 261-272. 10.1016/j.jchromb.2006.09.033.View ArticleGoogle Scholar
- Hunt SM, Thomas MR, Sebastian LT, Pedersen SK, Harcourt RL, Sloane AJ, Wilkins MR: Optimal replication and the importance of experimental design for gel-based quantitative proteomics. J Proteome Res. 2005, 4 (3): 809-819. 10.1021/pr049758y.View ArticlePubMedGoogle Scholar
- Karp NA, Spencer M, Lindsay H, O'Dell K, Lilley KS: Impact of replicate types on proteomic expression analysis. J Proteome Res. 2005, 4 (5): 1867-1871. 10.1021/pr050084g.View ArticlePubMedGoogle Scholar
- Hondermarck H: Breast cancer: when proteomics challenges biological complexity. Mol Cell Proteomics. 2003, 2: 281-291.PubMedGoogle Scholar
- Craven R, Banks R: Laser campture microdissection and proteomics: Possibilities and limitation. Proteomics. 2001, 1: 1200-1204. 10.1002/1615-9861(200110)1:10<1200::AID-PROT1200>3.0.CO;2-Q.View ArticlePubMedGoogle Scholar
- Shin KH, Ku JL, Kim WH, Lee SE, Lee C, Kim SW, Park JG: Establishment and characterization of seven human renal cell carcinoma cell lines. Br J Urol Int. 2000, 85 (1): 130-138.View ArticleGoogle Scholar
- Wistuba II, Behrens C, Milchgrub S, Syed S, Ahmadian M, Virmani AK, Kurvari V, Cunningham TH, Ashfaq R, Minna JD, Gazdar AF: Comparison of features of human breast cancer cell lines and their corresponding tumors. Clin Cancer Res. 1998, 4: 2931-2938.PubMedGoogle Scholar
- Wistuba II, Bryant D, Behrens C, Milchgrub S, Virmani AK, Ashfaq R, Minna JD, Gazdar AF: Comparison of features of human lung cancer cell lines and their corresponding tumors. Clin Cancer Res. 1999, 5: 991-1000.PubMedGoogle Scholar
- Henkin JA, Jennings ME, Matthews DE, Vigoreaux JO: Mass Processing - An Improved Technique for Protein Identification with Mass Spectrometry Data. J Biomol Tech. 2004, 4: 230-237.Google Scholar
- Klinke DJ: An empirical Bayesian approach for model-based inference of cellular signaling networks. BMC Bioinformatics. 2009, 10 (1): 371-10.1186/1471-2105-10-371.View ArticlePubMedPubMed CentralGoogle Scholar
- Gelman A, Rubin DB: Inference from iterative simulation using multiple sequences. Stat Sci. 1992, 7: 457-511. 10.1214/ss/1177011136.View ArticleGoogle Scholar
- Harry JL, Wilkins MR, Herbert BR, Packer NH, Gooley AA, Williams KL: Proteomics: Capacity versus utility. Electrophoresis. 2000, 21: 1071-1081. 10.1002/(SICI)1522-2683(20000401)21:6<1071::AID-ELPS1071>3.0.CO;2-M.View ArticlePubMedGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol. 1995, 57 (1): 125-133.Google Scholar
- Pezzino V, Papa V, Milazzo G, Gliozzo B, Russo P, Scalia PL: Insulin-like growth factor-I (IGF-I) receptors in breast cancer. Ann N Y Acad Sci. 1996, 784: 189-201. 10.1111/j.1749-6632.1996.tb16236.x.View ArticlePubMedGoogle Scholar
- Papa V, Gliozzo B, Clark GM, McGuire WL, Moore D, Fujita-Yamaguchi Y, Vigneri R, Goldfine ID, Pezzino V: Insulin-like growth factor-I receptors are overexpressed and predict a low risk in human breast cancer. Cancer Res. 1993, 53 (16): 3736-3740.PubMedGoogle Scholar
- Peyrat JP, Bonneterre J: Type 1 IGF receptor in human breast diseases. Breast Cancer Res Treat. 1992, 22: 59-67. 10.1007/BF01833334.View ArticlePubMedGoogle Scholar
- Resnik JL, Reichart DB, Huey K, Webster N, Seely BL: Elevated Insulin-like Growth Factor I Receptor Autophosphorylation and Kinase Activity in Human Breast Cancer. Cancer Res. 1998, 58: 1159-1164.PubMedGoogle Scholar
- von Lintig FC, Dreilinger AD, Varki NM, Wallace AM, Casteel DE, Boss GR: Ras activation in human breast cancer. Breast Cancer Res Treat. 2000, 62 (1): 51-62. 10.1023/A:1006491619920.View ArticlePubMedGoogle Scholar
- Molloy MP, Herbert BR, Williams KL, Gooley AA: Extraction of Escherichia coli proteins with organic solvents prior to two-dimensional electrophoresis. Electrophoresis. 1999, 20: 701-704. 10.1002/(SICI)1522-2683(19990101)20:4/5<701::AID-ELPS701>3.0.CO;2-5.View ArticlePubMedGoogle Scholar
- Duncan MW, Hunscuker SW: Comments on standards in proteomics and the concept of fitness-for-purpose. Proteomics. 2006, 6 (Suppl.1-2 (Practical Proteomics)): 45-47. 10.1002/pmic.200600559.View ArticlePubMedGoogle Scholar
- Neve RM, Chin K, Fridlyand J, Yeh J, Baehner FL, Fevr T, Clark L, Bayani N, Coppe JP, Tong F, Speed T, Spellman PT, DeVries S, Lapuk A, Wang NJ, Kuo WL, Stilwell JL, Pinkel D, Albertson DG, Waldman FM, McCormick F, Dickson RB, Johnson MD, Lippman M, Ethier S, Gazdar A, Gray JW: A collection of breast cancer cell lines for the study of functionally distinct cancer subtypes. Cancer Cell. 2006, 10 (6): 515-527. 10.1016/j.ccr.2006.10.008.View ArticlePubMedPubMed CentralGoogle Scholar
- Weitzel LB, Byers T, Allen J, Finlayson C, Helmke SM, Hokanson JE, Hunsucker SW, Murphy JR, Newell K, Queensland KM, Singh M, Wischmeyer PE, Duncan MW, Elias A: Discovery and verification of protein differences between Er positive/Her2/neu negative breast tumor tissue and matched adjacent normal breast tissue. Breast Cancer Res Treat. 2010Google Scholar
- Lasfargues EY, Coutinho WG, Redfield ES: Isolation of two human tumor epithelial cell lines from solid breast carcinomas. J Natl Cancer Inst. 1978, 61 (4): 961-978.Google Scholar
- Trempe GL: Human breast cancer in culture. Recent Results Cancer Res. 1976, 57: 33-41.PubMedGoogle Scholar
- Karey KP, Sirbasku DA: Differential responsiveness of human breast cancer cell lines MCF-7 and T47D to growth factors and 17 beta-estradiol. Cancer Res. 1998, 48: 4083-4092.Google Scholar
- Ciampolillo A, De Tullio C, Giorgino F: The IGF-I/IGF-I receptor pathway: implications in the pathophysiology of thyroid cancer. Curr Med Chem. 2005, 12: 2881-2891. 10.2174/092986705774454715.View ArticlePubMedGoogle Scholar
- Tanno S, Mitsuuchi Y, Altomare DA, Xiao GH, Testa JR: AKT activation up-regulates insulin-like growth factor I receptor expression and promotes invasiveness of human pancreatic cancer cells. Cancer Res. 2001, 61: 589-593.PubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://0-www.biomedcentral.com.brum.beds.ac.uk/1471-2407/10/291/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.