About this Journal Submit a Manuscript Table of Contents
Journal of Oncology
Volume 2012 (2012), Article ID 413256, 9 pages
http://dx.doi.org/10.1155/2012/413256
Clinical Study

Salivary Protein Profiles among HER2/neu-Receptor-Positive and -Negative Breast Cancer Patients: Support for Using Salivary Protein Profiles for Modeling Breast Cancer Progression

Department of Diagnostic Sciences, UTHSC, Dental Branch, Behavioral and Biomedical Sciences Building, Rm. 5322, Houston, TX 77054, USA

Received 25 October 2011; Accepted 17 January 2012

Academic Editor: Esra Gunduz

Copyright © 2012 Charles F. Streckfus et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Purpose. The objective of this study was to compare the salivary protein profiles from individuals diagnosed with breast cancer that were either HER2/neu receptor positive or negative. Methods. Two pooled saliva specimens underwent proteomic analysis. One pooled specimen was from women diagnosed with stage IIa HER2/neu-receptor-positive breast cancer patients ( 𝑛 = 1 0 ) and the other was from women diagnosed with stage IIa HER2/neu-receptor-negative cancer patients ( 𝑛 = 1 0 ). The pooled samples were trypsinized and the peptides labeled with iTRAQ reagent. Specimens were analyzed using an LC-MS/MS mass spectrometer. Results. The results yielded approximately 71 differentially expressed proteins in the saliva specimens. There were 34 upregulated proteins and 37 downregulated proteins.

1. Introduction

Clinicopathologic factors such as histological type, tumor size, tumor grade, hormone receptor status, lymph node involvement, and HER-2/neu overexpression are recognized as having prognostic use in breast cancer management. HER-2/neu (HER2), also known as c-erbB-2, is a biomarker assayed in tissue biopsies from women diagnosed with malignant breast tumors [1, 2]. Used primarily as a prognostic indicator, HER2/neu protein is overexpressed in approximately 20%–30% of malignant breast tumors and has been used in postoperative followup evaluation as an indicator of patient relapse [36].

The evolution of HER2 testing, first as a prognostic marker assay and later as a diagnostic test to determine eligibility for trastuzumab-targeted therapy, has expanded the role of traditional diagnostic pathology. Unlike most testing performed by anatomic pathologists, which serves as an adjunct to establishing a diagnosis, the results of HER2 testing stand alone in determining which patients are likely to respond to trastuzumab therapy. HER2 status may also predict sensitivity to certain cytotoxic drugs and antiestrogens [36].

Currently, two testing methods are approved by the Food and Drug Administration (FDA) for HER2 assessment in the laboratory [1, 2]. They are immunohistochemical analysis (IHC) and fluorescence in situ hybridization (FISH). Commercially available, FDA-approved HER2 assays are available for both methods. Immunohistochemical analysis and FISH have the advantage over other assay methods (i.e., those requiring homogenization) because they are morphologically driven. This allows for the direct evaluation of tumor cells, correlation with other morphologic features, and the ability to assay smaller patient samples such as needle core biopsy specimens [1, 2].

While clinical treatment choices are critical, the actual tests used to determine HER2 status have demonstrated a number of pitfalls. One such pitfall is the number of false negatives and false positives associated with the tests. This creates a treatment dilemma as it can result in situations where patients requiring trastuzumab-targeted therapy may not receive it, while those receiving trastuzumab-targeted therapy should not receive it [68].

Not everyone is convinced that the problem is as simple as false negatives and false positives. One problem may be with the cutoff points that scientists have established to delineate between HER2 negative and positive in IHC testing. Data from several preclinical and clinical studies suggest that trastuzumab activity does not strictly require HER2 overexpression or gene amplification, as is currently thought. Instead, even tumor cells that express a lower level of the protein might respond to a trastuzumab-chemotherapy combination. Currently, tumors that have moderate amounts of HER2 protein would be scored as 1 + or 2 + on IHC tests and would be called negative because they are below the predetermined cutoff value for the test [18].

Even FISH testing, which is considered the “gold standard” for detecting gene amplification, has problems. The test measures the ratio between the area surrounding the HER2 gene on chromosome 17 and other parts of chromosome 17. That means a cell that has extra copies of chromosome 17, called polysomy 17, appears negative by FISH but actually has extra copies of HER2, and probably expresses more protein than a wild-type cell would [1, 2].

The inconsistency in the test may stem, in part, from heterogeneity of HER2 expression in the tumor sample. If so, then multiple tests on the same tumor could yield different results, even though the tests are working as designed. The testing problem has been a matter of discussion, and how it will be resolved remains unclear. However, clinicians and researchers agree that current technology to assess HER2 status needs improvement [18].

Therefore, the purpose of this study was twofold. The first objective was to compare salivary protein profiles among Her2-receptor-positive and -negative breast cancer patients and secondly, to support the theory for using salivary protein profile expression as a method for modeling breast cancer progression [9, 10].

2. Methods

2.1. Design

The investigators protein profiled three pooled, stimulated whole saliva specimens. One specimen consisted of pooled saliva from 10 stage IIa (T2N0M0) HER2-receptor-status-positive invasive ductal carcinoma patients (IDC), and the second pooled specimen was from 10 subjects diagnosed with stage IIb (T2N0M0), HER2-receptor-status-positive invasive ductal carcinoma. The cancer cohorts were estrogen, progesterone receptor status negative as determined by the pathology report. Histological grade was not available for this study. The subjects were matched for age and race and were nontobacco users.

The participating subjects were given an explanation about their participation rights, and they signed an IRB consent form. The saliva specimens and related patient data are nonlinked and bar-coded in order to protect patient confidentiality. This study was performed under the UTHSC IRB-approved protocol number HSC-DB-05-0394. All procedures were in accordance with the ethical standards of the UTHSC IRB and with the Helsinki Declaration of 1975, as revised in 1983.

2.2. Saliva Collection and Sample Preparation

Stimulated whole salivary gland secretion is based on the reflex response occurring during the mastication of a bolus of food. Usually, a standardized bolus (1 gram) of paraffin or a gum base (generously provided by the Wrigley Co., Peoria, IL) is given to the subject to chew at a regular rate. The individual, upon sufficient accumulation of saliva in the oral cavity, expectorates periodically into a preweighed disposable plastic cup. This procedure is continued for a period of five minutes. The volume and flow rate is then recorded along with a brief description of the specimen’s physical appearance [9, 10]. The cup with the saliva specimen is reweighed and the flow rate determined gravimetrically. A protease inhibitor from Sigma Co. (St. Louis, MI, USA) is added along with enough orthovanadate from a 100 mM stock solution to bring its concentration to 1 mM. The treated samples were centrifuged for 10 minutes at top speed in a tabletop centrifuge. The supernatant was divided into 1 mL aliquots and frozen at −80.

2.3. Two-Dimensional Gel Analysis (2D DIGE)
2.3.1. Sample Preparation

2D DIGE and protein ID was performed by Applied Biomics, Inc. (Hayward, CA). Proteins from the saliva were precipitated by methanol then resuspended in a 2D cell lysis buffer (30 mM Tris-HCl, pH 8.8, containing 7 M urea, 2 M thiourea, and 4% CHAPS). Protein concentration was measured using Bio-Rad protein assay method.

2.3.2. CyDye Labeling

For each sample, 30 ug of protein was mixed with 1.0 μL of diluted CyDye and kept in dark on ice for 30 min. Samples from each pair were labeled with Cy2 and Cy5, respectively. The labeling reaction was stopped by adding 1.0 μL of 10 mM Lysine to each sample, and incubating in dark on ice for additional 15 min. The labeled samples were then mixed together. The 2X 2D sample buffer (8 M urea, 4% CHAPS, 20 mg/mL DTT, 2% pharmalytes, and trace amount of bromophenol blue), 100 μL destreak solution, and Rehydration buffer (7 M urea, 2 M thiourea, 4% CHAPS, 20 mg/mL DTT, 1% pharmalytes, and trace amount of bromophenol blue) were added to the labeling mix to make the total volume of 250 μL. They were Mixed well and spined before loading the labeled samples into strip holder.

2.3.3. IEF and SDS-PAGE

After loading the labeled samples, IEF (pH3-10 Linear) was run following the protocol provided by Amersham BioSciences. Upon finishing the IEF, the IPG strips were incubated in the freshly made equilibration buffer 1 (50 mM Tris-HCl, pH 8.8, containing 6 M urea, 30% glycerol, 2% SDS, trace amount of bromophenol blue, and 10 mg/mL DTT) for 15 minutes with gentle shaking. Then the strips were rinsed in the freshly made equilibration buffer 2 (50 mM Tris-HCl, pH 8.8, containing 6 M urea, 30% glycerol, 2% SDS, trace amount of bromophenol blue, and 45 mg/mL DTT) for 10 minutes with gentle shaking. Next, the IPG strips were rinsed in the SDS-gel running buffer before transferring into 13.5% SDS gels. The SDS gels were run at 15°C until the dye front exuded out of the gels.

2.3.4. Image Scan and Data Analysis

Gel images were scanned immediately following the SDS-PAGE using Typhoon TRIO (Amersham BioSciences). The scanned images were then analyzed by Image Quant software (version 6.0, Amersham BioSciences), followed by in-gel analysis using DeCyder software version 6.0 (Amersham BioSciences). The fold change of the protein expression levels was obtained from in-gel DeCyder analysis.

2.4. Top-Down Mass Spectrometry Using iTRAQ Labeling

A thorough explanation for the top-down mass spectrometry using iTRAQ labeling can be found in detail in previous publications [9, 10]. Briefly, the saliva samples were thawed and immediately centrifuged to remove insoluble materials. The supernatant was assayed for protein using the Bio-Rad protein assay (Hercules, CA, USA), and an aliquot containing 100 μg of each specimen was precipitated with six volumes of −20°C acetone. The precipitate was resuspended and treated according to the manufacturer’s instructions. Protein digestion and reaction with iTRAQ labels was carried out as previously described and according to the manufacturer’s instructions (Applied Biosystems, Foster City, CA). Briefly, the acetone precipitable protein was centrifuged in a tabletop centrifuge at 15,000 ×g for 20 minutes. The acetone supernatant was removed and the pellet resuspended in 20 ŲL dissolution buffer. The soluble fraction was denatured, and disulfides were reduced by incubation in the presence of 0.1% SDS and 5 mM TCEP (tris-(2-carboxyethyl) phosphine) at 60°C for one hour. Cysteine residues were blocked by incubation at room temperature for 10 minutes with MMTS (methyl methane-thiosulfonate). Trypsin was added to the mixture to a protein: trypsin ratio of 10 : 1. The mixture was incubated overnight at 37°C. The protein digests were labeled by mixing with the appropriate iTRAQ reagent and incubating at room temperature for one hour. On completion of the labeling reaction, the four separate iTRAQ reaction mixtures were combined. Since there are a number of components that can interfere with the LCMSMS analysis, the labeled peptides are partially purified by a combination of strong cation exchange followed by reverse-phase chromatography on preparative columns. The combined peptide mixture is diluted 10-fold with loading buffer (10 mM, KH2PO4 in 25% acetonitrile at pH 3.0) and applied by syringe to an ICAT Cartridge-Cation Exchange column (Applied Biosystems, Foster City, CA) column that has been equilibrated with the same buffer. The column is washed with 1 mL loading buffer to remove contaminants. To improve the resolution of peptides during LCMSMS analysis, the peptide mixture is partially purified by elution from the cation exchange column in three fractions. Stepwise elution from the column is achieved with sequential 0.5 mL aliquots of 10 mM KH2PO4 at pH 3.0 in 25% acetonitrile containing 116 mM, 233 mM, and 350 mM KCl, respectively. The fractions are evaporated by Speed Vac to about 30% of their volume to remove the acetonitrile and then slowly applied to an Opti-Lynx Trap C18 100 μL reverse-phase column (Alltech, Deerfield, IL) with a syringe. The column was washed with 1 mL of 2% acetonitrile in 0.1% formic acid and eluted in one fraction with 0.3 mL of 30% acetonitrile in 0.1% formic acid. The fractions were dried by lyophilization and resuspended in 10 μL 0.1% formic acid in 20% acetonitrile solution. Each of the three fractions was analyzed by reverse-phase nano-LCMS/MS on an API QSTAR XL mass spectrometer (ABS Sciex Instruments).

2.5. Bioinformatics

The Swiss-Prot database was employed for protein identification, while the PathwayStudio bioinformatics software package was used to determine Venn diagrams that were also constructed using the NIH software program (http://ncrr.pnl.gov/). Graphic comparisons with log conversions and error bars for protein expression were produced using the ProQuant software. Candidates with either protein score C.I. percentage or Ion C.I. percentage greater than 95 were considered significant.

3. Results

3.1. 2D Gel Results

Figures 1 and 2 illustrate the results of the 2D gel analyses. Figure 1 demonstrates the protein comparisons between the pooled HER2/neu-receptor-positive and the HER2/neu-receptor-negative pooled specimens. Figure 2 represents the spots of interest that were selected up by Ettan Spot Picker (Amersham BioSciences) based on the in-gel analysis and spot picking design by DeCyder software. As shown, there are 96 spots of interest illustrated on the 2D gel analysis. This visually indicates the differing salivary protein patterns between HER2/neu-positive and HER2/neu-negative patients.

fig1
Figure 1: This figure illustrates the protein profiles for HER2/neu-receptor-positive and HER2/neu-receptor-negative samples. As shown in the far right red and green dyed gel comparisons, there are numerous differences between the two profiles.
413256.fig.002
Figure 2: This figure demonstrates the differences in salivary protein profiles between HER2/neu -positive and HER2/neu -negative samples. Please change.
3.2. LC-MS/MS Mass Spectrometry Results

The results yielded 188 comparative salivary proteins among the HER2/neu-positive and HER2/neu-negative samples. Among the total number of proteins, 71 were significantly differentially expressed between the two specimens. There were 34 upregulated proteins and 37 downregulated proteins. Listed in Table 1 is the 34 upregulated proteins, and in Table 2, the 37 downregulated proteins.

tab1
Table 1: Upregulated salivary proteins.
tab2
Table 2: Downregulated salivary proteins.

Of the 34 proteins listed in Table 1, the mean percent peptide coverage for the complete panel of proteins was 50.3% ( ± 1 9 . 6 ) with a range of 15.3% to 88.7% coverage. The mean protein ratio was 1 . 6 4 ( ± 5 4 . 2 ) and ranged in value with a maximum of 3.13 to a maximum minimum of 1.12. Likewise, the mean alpha level was 𝑃 < 0 . 0 0 7 ( ± 0 . 0 1 3 ) and ranged in value with a maximum of 𝑃 < 0 . 0 4 to a maximum minimum of 𝑃 < 0 . 0 0 0 0 0 1 .

Of the 37 downregulated salivary proteins listed in Table 2, the mean percent peptide coverage for the complete panel of proteins was 40.1% ( ± 2 0 . 7 ) with a range of 9.7% to 100% coverage. The mean protein ratio was 1.64 ( ± 0 . 5 6 1 ) and ranged in value with a maximum of 0.89 to a maximum minimum of 0.10. Likewise, the mean alpha level was 𝑃 < 0 . 0 0 9 ( ± 0 . 0 1 5 ) and ranged in value with a maximum of 𝑃 < 0 . 0 5 to a maximum minimum of 𝑃 < 0 . 0 0 0 0 0 1 . Figure 3 illustrates the comparison of the log ratio of the relative intensity (HER2+/HER2−) of the total 71 up- and downregulated salivary proteins. Additionally, Figure 4 illustrates the protein function of the panel of 71 proteins. Proteins related to cellular metabolism and immune response comprised nearly 42% of the protein panel. Cellular structure constituted 17% of the protein panel, which is consistent with SKBR3 and other cell line protein analyses [15]. Additionally, there were a considerable number of histones that were present in the functional protein profile.

413256.fig.003
Figure 3: The figure represents the natural logarithm differential expression of salivary proteins. To the right and left of the figure listed in rank order of expression are the up- and downregulated proteins, respectively.
413256.fig.004
Figure 4: This figure represents the percentage of expressed proteins according to their cellular function.

4. Discussion

To the best of our knowledge, this is the first attempt to determine salivary protein profile alterations related to HER2/neu receptor status. Therefore, we have only a few references by which to compare our data. Additionally, a complete proteomic catalog of the SKBR3 cell line is also not available to compare to the salivary protein profiles, which are altered secondary to HER2 receptor status. Of the articles that were identified through the PubMed Central (United States National Library of Medicine) search engine [1117], twenty (28%) of the 71 of the salivary proteins were found cited among the manuscripts [1117] and reported to be altered in the SKBR3 cell lines. These salivary proteins are listed in Table 3 with the references that cited the SKBR3 protein phenotype alterations.

tab3
Table 3: Altered protein in saliva and in SKBR3 cell lines.

The numerous proteins in the panel need validation; however, the authors selected Profilin-1 as a test case to determine if the results of the proteomic analyses are feasible. As illustrated in Figure 5, Profilin-1 is present in both the SKBR3 and the human salivary gland (HSG) epithelial cell lines. Additionally, Profilin-1 is present in saliva and is downregulated to HER2-receptor-positive status. The pilot evidence supports the proteomic prediction.

fig5
Figure 5: It illustrates the presence of profilin-1 in the SKBR3 and HSG cell lysates and saliva sampled from healthy, benign, and malignant tumor patients. Profilin is a downregulated protein in the presence of malignancy, and it is visualized by the lighter bands associated with malignancy. It is also worth noting that the Her2/neu-receptor-negative band is darker than the Her2/neu-receptor-positive counterpart suggesting further downregulation of the profiling-1 protein.

It is worth noting that among the list (Table 3) of SKBR3 cell lysates proteins, there is a distinct absence of any proteins related to immune response to inflammatory cancer activity as compared to the salivary proteins listed in both Tables 1 and 2. The findings at this point suggest the strength of this in vivo model which could be indicative of response to therapy in the event these proteins are diminished in activity during treatment, thereby, indicating a response to therapy.

Further research is required to support the theories presented in this paper. For example, proteomic analyses of low-abundance proteins in the SKRB3 cell line lysates and saliva are required in order to address gaps in varying molecular pathways. Studies validating the panel of markers are also necessary, and an assessment of salivary protein modulation during trastuzumab therapy is essential.

5. Conclusions

The results of the study suggest salivary protein alterations secondary to HER2 receptor status. This is not surprising considering that the ductal cells of the salivary glands contain HER2/neu receptors. More importantly, the study raises the notion that salivary gland protein secretions may be used as a “real-time”, in vivo model for studying breast cancer progression [9, 10]. Currently, there are three major methods for creating models for studying breast cancer progression [18]. The three methods utilize either breast cancer tumor cell lines xenografts of cell lines, and the third method uses animals—in this case genetically engineered mice [19] for creating various models for studying breast cancer [18]. All three models have generated useful insight into cancer progression; however, despite their utility, no individual model recapitulates all aspects of cancer progression [18, 19].

Hence, an adjunct, in vivo model system is needed for breast cancer tumorigenesis and predictive modeling for treatment response [18, 19]. The authors have demonstrated in previous studies that the salivary protein profiles are altered in the presence of ductal carcinoma in situ and are further altered in the presence of lymph node involvement [9, 10]. The preliminary findings of this paper coupled with previous studies do imply that this in vivo experimental model system, which utilizes one of the most easily obtained body fluids for marker analysis, may fill in the current gaps in our understanding of breast cancer pathogenesis, signaling pathways, the efficacy of varying chemotherapeutics, and identifying novel therapies. Most importantly, this new approach may shed new light on metastatic progression that is the principle cause of patient mortality.

Acknowledgments

The Avon Breast Cancer Foundation (no. 07-2007-071), Komen Foundation (KG080928), Gillson-Longenbaugh Foundation, and the Texas Ignition Fund supported the research presented in this paper. The authors would also like to thank Dr. William Dubinsky for the LC-MS/MS salivary mass spectrometry analyses.

References

  1. A. C. Wolff, M. E. H. Hammond, J. N. Schwartz et al., “American Society of Clinical Oncology/College of American Pathologists guideline recommendations for human epidermal growth factor receptor 2 testing in breast cancer,” Journal of Clinical Oncology, vol. 25, no. 1, pp. 118–145, 2007. View at Publisher · View at Google Scholar · View at Scopus
  2. M. E. Hammond, D. F. Hayes, and A. C. Wolff, “Clinical Notice for American Society of Clinical Oncology-College of American Pathologists guideline recommendations on ER/PgR and HER2 testing in breast cancer,” Journal of Clinical Oncology, vol. 29, no. 15, p. e458, 2011. View at Publisher · View at Google Scholar
  3. J. S. Ross and J. A. Fletcher, “The HER-2/neu oncogene in breast cancer: prognostic factor, predictive factor, and target for therapy,” Stem Cells, vol. 16, no. 6, pp. 413–428, 1998. View at Publisher · View at Google Scholar · View at Scopus
  4. S. Koscielny, P. Terrier, M. Spielmann, and J. C. Delarue, “Prognostic importance of low c-erbB2 expression in breast tumors,” Journal of the National Cancer Institute, vol. 90, no. 9, p. 712, 1998. View at Publisher · View at Google Scholar · View at Scopus
  5. W. Y. Huang, B. Newman, R. C. Millikan et al., “Risk of breast cancer according to the status of HER-2/neu oncogene amplification,” Cancer Epidemiology Biomarkers & Prevention, vol. 9, no. 1, pp. 65–71, 2000. View at Scopus
  6. L. R. Howe and P. H. Brown, “Targeting the HER/EGFR/ErbB family to prevent breast cancer,” Cancer Prevention Research, vol. 4, no. 8, pp. 1149–1157, 2011. View at Publisher · View at Google Scholar
  7. J. Baselga, “Phase I and II clinical trials of trastuzumab,” Annals of Oncology, vol. 12, supplement 1, pp. S49–S55, 2001. View at Scopus
  8. S. Ménard, P. Valagussa, S. Pilotti et al., “Response to cyclophosphamide, methotrexate, and fluorouracil in lymph node-positive breast cancer according to HER2 overexpression and other tumor biologic variables,” Journal of Clinical Oncology, vol. 19, no. 2, pp. 329–335, 2001. View at Scopus
  9. C. F. Streckfus, O. Mayorga-Wark, D. Arreola, C. Edwards, L. Bigler, and W. P. Dubinsky, “Breast cancer related proteins are present in saliva and are modulated secondary to ductal carcinoma in situ of the breast,” Cancer Investigation, vol. 26, no. 2, pp. 159–167, 2008. View at Publisher · View at Google Scholar · View at Scopus
  10. C. F. Streckfus, K. A. Storthz, L. Bigler, and W. P. Dubinsky, “A comparison of the proteomic expression in pooled saliva specimens from individuals diagnosed with ductal carcinoma of the breast with and without lymph node involvement,” Journal of Oncology, vol. 2009, Article ID 737619, 11 pages, 2009. View at Publisher · View at Google Scholar
  11. Y. M. Kulkarni, V. Suarez, and D. J. Klinke, “Inferring predominant pathways in cellular models of breast cancer using limited sample proteomic profiling,” BMC Cancer, vol. 10, article 291, 2010. View at Publisher · View at Google Scholar · View at Scopus
  12. S. J. Pitteri, K. S. Kelly-Spratt, K. E. Gurley et al., “Tumor microenvironment-derived proteins dominate the plasma proteome response during breast cancer induction and progression,” Cancer Research, vol. 71, no. 15, pp. 5090–5100, 2011. View at Publisher · View at Google Scholar
  13. S. L. Wu, W. S. Hancock, G. G. Goodrich, and S. T. Kunitake, “An approach to the proteomic analysis of a breast cancer cell line (SKBR-3),” Proteomics, vol. 3, no. 6, pp. 1037–1046, 2003. View at Publisher · View at Google Scholar · View at Scopus
  14. E. Charafe-Jauffret, C. Ginestier, F. Monville et al., “Gene expression profiling of breast cell lines identifies potential new basal markers,” Oncogene, vol. 25, no. 15, pp. 2273–2284, 2006. View at Publisher · View at Google Scholar · View at Scopus
  15. R. A. Harris, A. Yang, R. C. Stein et al., “Cluster analysis of an extensive human breast cancer cell line protein expression map database,” Proteomics, vol. 2, no. 2, pp. 212–223, 2002. View at Publisher · View at Google Scholar · View at Scopus
  16. D. Zhang, L. K. Tai, L. L. Wong, L. L. Chiu, S. K. Sethi, and E. S. C. Koay, “Proteomic study reveals that proteins involved in metabolic and detoxification pathways are highly expressed in HER-2/neu-positive breast cancer,” Molecular & Cellular Proteomics, vol. 4, no. 11, pp. 1686–1696, 2005. View at Publisher · View at Google Scholar · View at Scopus
  17. J. Kao, K. Salari, M. Bocanegra et al., “Molecular profiling of breast cancer cell lines defines relevant tumor models and provides a resource for cancer gene discovery,” PLoS ONE, vol. 4, no. 7, Article ID e6146, 2009. View at Publisher · View at Google Scholar · View at Scopus
  18. T. Vargo-Gogola and J. M. Rosen, “Modelling breast cancer: one size does not fit all,” Nature Reviews Cancer, vol. 7, no. 9, pp. 659–672, 2007. View at Publisher · View at Google Scholar · View at Scopus
  19. L. Hennighausen, “Mouse models for breast cancer,” Breast Cancer Research, vol. 2, no. 1, pp. 2–7, 2000. View at Publisher · View at Google Scholar · View at Scopus