Increased expression of Six1 correlates with progression and prognosis of prostate cancer
© Zeng et al. 2015
Received: 11 May 2015
Accepted: 8 June 2015
Published: 20 June 2015
Sineoculis homeobox homolog 1 (Six1), normally a developmentally restricted transcriptional regulator, is frequently dysregulated in mutiple cancers. Increasing evidences show that overexpression of Six1 plays a key role in tumorigenesis. However, the Six1 expression status and its relationship with the clinicopathological characteristics in prostate cancer were unclear. In this study, the mRNA and protein levels of Six1 in prostate cancer tissues and normal prostate tissues were evaluated. The clinicopathological significance of Six1 was investigated by immunohistochemistry (IHC) on a prostate cancer tissue microarray. The cut-off score for high expression of Six1 was determined by the receiver-operating characteristic (ROC) analysis. The correlation between Six1 protein expression and clinicopathological characteristics of prostate cancer was analyzed by Chi-square test. Increased expression of Six1 protein was observed in the majority of prostate cancer, compared with their paired adjacent normal prostate tissues. When Six1 high expression percentage was determined to be above 55 % (area under ROC curve = 0.881, P = 0.000), high expression of Six1 was observed in 55.6 % (80/144) of prostate cancer tissues and low expression of Six1 was observed in all normal prostate tissues by IHC. Increased expression of Six1 in patients was correlated with high histological grade (χ2 = 58.651, P = 0.00), advanced clinical stage (χ2 = 57.330, P = 0.000), high Gleason score (χ2 = 63.480, P = 0.000), high primary tumor grade (χ2 = 57.330, P = 0.000) and positive regional lymph node metastasis (χ2 = 19.294, P = 0.000). Furthermore, univariate and multivariate survival analysis suggested that Six1 was an independent prognostic indicator for overall survival (P < 0.05). This study suggests that Six1 could be served as an additional biomarker in identifying prostate cancer patients at risk of tumor progression, might potentially be used for predicting survival outcome of patients with prostate cancer.
Current diagnosis for prostate cancer includes digital rectal examination (DRE), prostate-specific antigen (PSA) and needle biopsy . The appearance of prostate-specific antigen (PSA) testing has revolutionarily improved early prostate cancer detection. If rising PSA levels are detected, a needle biopsy of the prostate is recommended to inspect for histologic evidence of prostate cancer. If cancer is inspected, the patient can choose either active surveillance or one of several given treatment options, such as surgery, radiotherapy, or brachytherapy. However, current reports have raised concern over the efficacy of PSA testing. The U.S. prostate cancer screening trial report found that the mortality caused by prostate cancer was not reduced by PSA screening . These studies suggested that PSA testing cannot distinguish between aggressive and indolent prostate cancer; therefore, many patients accepted overtreatment. As a result of these kind reports, the U.S. Preventive Services Task Force recommended to stop routine screening by PSA testing on all males . Regardless of how much controversy about the PSA screening, there is an agreement that a definitive test is needed to distinguish patients with aggressive prostate cancer from the patients that have latent or indolent prostate cancer. Therefore, a critical question in the clinical activity of prostate cancer is how to distinguish the patients with indolent prostate cancer from patients with highly aggressive prostate cancer who would benefit from definitive treatment.
Homeobox genes encode transcription factors that are essential for the development of numerous organs and control processes such as proliferation, apoptosis, migration, and invasion [4–7]. The Six1 homeoprotein, a member of the Six family of homeodomain transcription factors, has been found overexpressed in mutiple cancers, including breast [8–10], rhabdomyosarcomas [11–13], hepatocellular carcinomas , ovarian  and Wilms tumors . Recent evidence demonstrates that Six1 plays a role in cellular migration and invasion during embryogenesis [4, 17–20] and in breast cancer [21, 22] through a mechanism that may involve an EMT, which was associated with increased TGF-β receptor type I (TβRI) expression and Smad-dependent transforming growth factor-beta (TGF-β) signaling. In addition, Six1 overexpression enhances lymph node metastasis by increasing VEGF-C depending on TGF-β signaling in cervical cancer [23–26]. Forementioned studies suggested that Six1 enhances cancer cell proliferation and shortens survival and its overexpression in immortalized mammary epithelial cells induces EMT, leading to highly aggressive and invasive tumors when transplanted into nude mice .
Noteworthy, HD Jin et al. reported that Six1 overexpression positively correlated with clinical stage, lymph node metastasis, overall survival (OS) and disease-free survival (DFS) rates of patients with breast cancer. Moreover, patients with high Six1 expression had poorer prognosis than those with low Six1 expression in late stage breast cancer cases . However, the relationship between Six1 and the clinicopathological characteristics of prostate cancer has not previously been examined.
Materials and methods
Patient specimens and prostate cancer cohorts
A total of 144 paraffin-embedded tissues diagnosed from 2002 to 2012 were retrieved for tissue microarray (TMA) construction and immunohistochemistry analysis. All the samples were collected from the Department of Urology, Nanfang Hospital (Guangzhou, China). The samples selected were pathologically diagnosed with prostate cancer, having received no prior chemotherapy or radiotherapy before surgery. Ages of the 144 patients with prostate cancer varied from 44 to 77 years (median, 65 years), clinicopathological features of patients including age at diagnosis, histologic grade, clinical stage, gleason score and pTNM stages. Written informed consent was obtained from all patients for use of the tissue samples and clinical records. The study protocol was performed under the approval by the Ethic Committee of the Nanfang Hospital. All cases were evaluated by experienced pathologists for histopathological grading.
Quantitative real-time PCR
qRT-PCR was performed to determine the expression of Six1 messenger RNA (mRNA). Briefly, we isolated total RNAs from frozen tissues using TRIzol Reagent according to the manufacturer’s protocol (Life Technologies) and reverse transcribed to generate cDNA (PrimeScript RT-PCR Kit, Takara Bio). β-Actin was used as an internal control. The levels of mRNA encoding were quantified by real-time PCR with the Applied Biosystems 7500 Fast Real-Time PCR System using SYBR Premix Ex Taq (Applied Takara Bio). The sequences of the primers were as follows: Six1 forward 5’- TTCTCGCCTCACAACCACCCCA-3’ and reverse 5’-TACCACTCCCGCAGGACACCCC-3’ and β-actin forward 5’-TGGCACCCAGCACAATGAA-3’ and reverse 5’-CTAAGTCATAGTCCGCCTAGAAGCA-3’. The PCR conditions were as following procedures: initial step 95 °C for 2 min, followed by 35 cycles of 95 °C for 30 s, 56 °C for 30 s, and 72 °C for 2 min and a final elongation step of 72 °C for 5 min. All qRT-PCRs were repeated three times. Relative quantification of Six1 mRNA expression was calculated using the 2−ΔΔCT method.
Western blotting analysis
Total proteins from 8 pairs of fresh prostate cancer and adjacent normal tissues were extracted by radio-immunoprecipitation assay (RIPA) buffer containing 1 mM phenylmethanesulfonyl fluoride (PMSF) (Beyotime, Haimen, China). After centrifugation, supernatant was collected and treated with BCA Protein Assay Kit (Beyotime, Haimen, China). The protein concentration was then measured at 562 nm -by the TECAN Infinite 200 microplate reader (TECAN, Switzerland). Tissue homogenates (35 μg protein for each sample) were separated by 10 % sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE), and the resolved proteins were transferred onto a polyvinylidene difluoride (PVDF) membrane (Millipore, USA) by a Trans-Blot SD Semi-dry transfer cell machine (Bio-Rad, USA). After the blots were washed with 1 × TBST buffer (10 mM Tris-HCl [pH 7.6], 150 mM NaCl, and 0.05 % Tween-20), the membranes were blocked overnight with 5 % skim milk and incubated with the appropriate primary antibody at room temperature for 2 h. Polyclonal rabbit anti-human antibody against Six1 (Atlas antibody, Sweden, 1:500), and monoclonal rabbit anti-human antibody against β-actin (Cell Signaling Technology, Danvers, MA, USA, 1:3000) were used for detecting the protein level of Six1 and β-actin in each sample. The membranes were then washed by 1 × TBST, primary antibodies were detected with goat anti-rabbit IgG conjugated to horseradish peroxidase (HRP) (Santa Cruz Biotechnology, Cal, USA, 1:5000), and the bands were detected by BeyoECL Plus Kit (Beyotime, Haimen, China). Finally, the results of western blotting were visualized by the Image Station 4000R PRO scanner (CareStream Health, Rochester, NY, U.S.A.).
Tissue microarray (TMA) construction and immunohistochemistry (IHC)
Association of Six1 expression with patients’ clinicopathologic features in prostate cancer
G1 + G2
I + II
III + IV
T1 + T2
T3 + T4
N1 + N2
Evaluation of IHC
Immunoreactivity for the Six1 protein was scored by a semi-quantitative method by recording the proportion of positive tumor cells over the total number of tumor cells. Scores were assigned by using 5 % increments (0, 5…100 %). The reproducibility of the scoring method between pathologists has been described previously for TMAs [28–31]. The scores were accepted if two of the three investigators (M.W, R.S and J.Z) agree with the values. Otherwise, the values were re-estimated until a consensus was reached. Our conclusions were in complete agreement in 80 % of the cases, which indicated that the scoring method was highly reproducible.
Selection of Cut-off scores using receiver- operating characteristic (ROC) curves
Receiver operating characteristic curve analysis was utilized to determine the cut-off score by using the 0, 1-criterion . At the Six1 score, the sensitivity and specificity for each outcome under study was plotted, thus generating various ROC curves. The score was selected as the cut-off score, which was closest to the point with both maximum sensitivity and specificity. Tumors designated as “low” for Six1 were those with scores below and equal to the cut-off value, whereas “high” tumors were those with scores above the value . In order to perform ROC curve analysis, the clinicopathological features were dichotomized: age (≤65 or >65), histological grade (low [G1 + G2] or high G3), clinical stage ([I + II] or above [III + IV]), Gleason stage (1–7 or 8–10), T stage (early [T1 + T2] or late [T3 + T4]), N stage (N0 [no lymph node involvement] or N1 + N2 [any lymph node involvement]), M stage (M0 [no distant metastasis] or M1 [distant metastasis]).
Statistical analysis was performed by using the SPSS statistical software program (standard version 18.0; SPSS, Chicago, USA). The relationship between Six1 protein expression and clinicopathological data of Prostate cancer patients was estimated using the χ 2 test. ROC curve analysis was applied to determine the cut-off score for Six1 positivity, and the areas under curves (AUCs) were then calculated. The association between survival and each variable was determined with the Kaplan-Meier method (the log-rank test). Multiple Cox proportional hazards regression was carried out to identify the independent prognostic factors. P < 0.05 (two-tailed) denotes the presence of a statistically significant difference.
The expression level of Six1 in prostate cancer and adjacent normal prostate tissues detected by western blotting and IHC
Selection of Six1 cut-off scores
Area under the ROC curve values for each clinicopathological feature
AUC (95 % CI)
Correlations between Six1 protein expression and clinicopathological parameters of prostate cancer
Six1 protein overexpression was significantly correlated with histological grade, clinical stage, Gleason score, primary tumor statues and regional lymph node metastasis of prostate cancer (P < 0.001). However, Six1 protein expression was not related with patient age, distant metastasis in prostate cancer (P > 0.05) (Table 1). For histological grades, the high expression rate of Six1 protein was significantly higher in late grade (G3) prostate cancers (89.6 %, 60/67) than in early histological grade cases (G1 + G2) (26.0 %, 20/77) (P < 0.001). Similarly, the high expression rate of Six1 protein was higher in breast cancers with high Clinical stage (75.0 %, 78/104) compared with those with low Clinical stage (0.5 %, 2/40) (P < 0.001). The high expression rate of Six1 protein was significantly higher in high Gleason score (8–10) prostate cancers (81.1 %, 73/90) than in low Gleason score cases (1–7) (13.0 %, 7/54) (P < 0.001). The high expression rate of Six1 protein was 75.0 % (78/104) in high pT grade, which was significantly higher than in low pT grade (0.5 %, 2/40) (P < 0.001). Meanwhile, the high expression rate of Six1 protein was significantly higher in high pN grade (N1 + N2) prostate cancers (81.1 %, 73/90) than in low pN grade cases (N0) (13.0 %, 7/54) (P < 0.001).
Relationship between clinicopathologic variables, Six1 protein expression and prostate patients’ survival
Univariate and Multivariate Analysis of Different Prognostic Features in 144 Patients with prostate cancer
HR (95 % CI)
G1 + G2
I + II
III + IV
T1 + T2
T3 + T4
N1 + N2
Independent prognostic factors of prostate cancer: multivariate survival analysis
A multivariate analysis based on the Cox proportional hazard model was applied to test the independent value of each parameter predicting overall survival (Table 3). Expressions of Six1 as well as other clinicopathological features that were significant by univariate analysis (histological grade, clinical stage, Gleason score, pT stage and pN stage) were included in multivariate analysis. The expression of Six1 was found to be an independent prognostic factor for adverse overall survival (relative risk, 2.322; 95 % CI, 1.163-4.636; P = 0.017). Of the other parameters, Clinical stage (P = 0.044), pT stage (P = 0.044) were demonstrated as well an independent prognostic factor for overall survival.
Prostate cancer is one of the most frequently diagnosed malignancies and second only to lung and bronchus cancer in men worldwide. Approximately 258,400 men dead from this disease in 2008 . At current rates of diagnosis, one-sixth of men will be diagnosed with prostate cancer during his lifetime . From the clinical perspective, the tumors in different patients must be varying in the molecular level and the aim of personalized medicine is to generate individual risk profiles from the primary prostate cancer that could distinguish high-risk individuals who should accept aggressive therapeutic treatment and clinical follow-up from those who have indolent prostate cancer to avoid these individuals from undergoing overtreatment. Recently, several studies have investigated molecular and genetic characteristics of prostate cancer to develop both predictive and prognostic biomarkers [35–37].
It was reported that in multiple types of human cancer, such as breast, rhabdomyosarcomas, hepatocellular carcinomas, ovarian, and Wilms tumors, overexpression of Six1 was frequently identified [8–16]. In the study, we investigated the expression patterns of Six1, by Western blotting using fresh prostate tissues, and by IHC using TMA containing a large cohort of prostate tumor samples (144 cases) with complete clinicopathological and follow-up data. Western blotting revealed that up-regulation expression of Six1 was detected in prostate cancer, when compared with adjacent normal prostate tissues. Moreover, the IHC results demonstrated that an increasing expression of Six1 was observed from low malignant tumors to high malignant tumors. Six1 immunoreactivity was assessed using a scoring method described in previous study . The reliability of this scoring method for Six1 was assessed by three investigators and was again found to be highly reproducible. To select the IHC cut-off scores for Six1 positivity, ROC analysis was carried out for each of the clinicopathological parameters, including age, histological grade, clinical stage, Gleason score, pTNM stage. ROC analysis for different clinicopathological features and Six1 expression was also used to evaluate the survival status. Six1 expression demonstrated significant statistical associations with the survival status of patients with prostate cancer. These data provided evidence that the unregulated expression of Six1 played an important role in tumorigenic process of multiple human cancers, including prostate cancer.
In summary, the study suggests that increased expression of Six1 may represent a more aggressive status of prostate cancer. In addition, these results show that increased expression of Six1 predicts poorer outcome of the cancer for the patients compared to those patients with decreased expression of Six1. Further studies are needed to clarify the precise mechanism of Six1 in prostate cancer and to develop potential therapies targeting Six1 in prostate cancer.
- Velonas VM, Woo HH, Remedios CG, Assinder SJ. Current status of biomarkers for prostate cancer. Int J Mol Sci. 2013;14(6):11034–60.PubMed CentralPubMedView ArticleGoogle Scholar
- Andriole GL, Crawford ED, Grubb III RL, Buys SS, Chia D, Church TR, et al. Mortality results from a randomized prostate-cancer screening trial. N Engl J Med. 2009;360(13):1310–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Moyer VA. Screening for prostate cancer: U.S. Preventive Services Task Force recommendation statement. Ann Intern Med. 2012;157(2):120–34.PubMedView ArticleGoogle Scholar
- Zheng W, Huang L, Wei ZB, Silvius D, Tang B, Xu PX. The role of Six1 in mammalian auditory system development. Development. 2003;130(17):3989–4000.PubMed CentralPubMedView ArticleGoogle Scholar
- Ikeda K, Kageyama R, Suzuki Y, Kawakami K. Six1 is indispensable for production of functional progenitor cells during olfactory epithelial development. Int J Dev Biol. 2010;54(10):1453–64.PubMedView ArticleGoogle Scholar
- Xu PX, Zheng W, Huang L, Maire P, Laclef C, Silvius D. Six1 is required for the early organogenesis of mammalian kidney. Development. 2003;130(14):3085–94.PubMed CentralPubMedView ArticleGoogle Scholar
- Coletta RD, Christensen KL, Micalizzi DS, Jedlicka P, Varella-Garcia M, Ford HL. Six1 overexpression in mammary cells induces genomic instability and is sufficient for malignant transformation. Cancer Res. 2008;68(7):2204–13.PubMedView ArticleGoogle Scholar
- Ford HL, Kabingu EN, Bump EA, Mutter GL, Pardee AB. Abrogation of the G2 cell cycle checkpoint associated with overexpression of HSIX1: a possible mechanism of breast carcinogenesis. Proc Natl Acad Sci U S A. 1998;95(21):12608–13.PubMed CentralPubMedView ArticleGoogle Scholar
- Coletta RD, Christensen K, Reichenberger KJ, Lamb J, Micomonaco D, Huang L, et al. The Six1 homeoprotein stimulates tumorigenesis by reactivation of cyclin A1. Proc Natl Acad Sci U S A. 2004;101(17):6478–83.PubMed CentralPubMedView ArticleGoogle Scholar
- Reichenberger KJ, Coletta RD, Schulte AP, Varella-Garcia M, Ford HL. Gene amplification is a mechanism of Six1 overexpression in breast cancer. Cancer Res. 2005;65(7):2668–75.PubMedView ArticleGoogle Scholar
- Khan J, Bittner ML, Saal LH, Teichmann U, Azorsa DO, Gooden GC, et al. cDNA microarrays detect activation of a myogenic transcription program by the PAX3-FKHR fusion oncogene. Proc Natl Acad Sci U S A. 1999;96(23):13264–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Yu Y, Khan J, Khanna C, Helman L, Meltzer PS, Merlino G. Expression profiling identifies the cytoskeletal organizer ezrin and the developmental homeoprotein Six-1 as key metastatic regulators. Nat Med. 2004;10(2):175–81.PubMedView ArticleGoogle Scholar
- Yu Y, Davicioni E, Triche TJ, Merlino G. The homeoprotein six1 transcriptionally activates multiple protumorigenic genes but requires ezrin to promote metastasis. Cancer Res. 2006;66(4):1982–9.PubMedView ArticleGoogle Scholar
- Ng KT, Man K, Sun CK, Lee TK, Poon RT, Lo CM, et al. Clinicopathological significance of homeoprotein Six1 in hepatocellular carcinoma. Br J Cancer. 2006;95(8):1050–5.PubMed CentralPubMedView ArticleGoogle Scholar
- Behbakht K, Qamar L, Aldridge CS, Coletta RD, Davidson SA, Thorburn A, et al. Six1 overexpression in ovarian carcinoma causes resistance to TRAIL-mediated apoptosis and is associated with poor survival. Cancer Res. 2007;67(7):3036–42.PubMedView ArticleGoogle Scholar
- Li CM, Guo M, Borczuk A, Powell CA, Wei M, Thaker HM, et al. Gene expression in Wilms’ tumor mimics the earliest committed stage in the metanephric mesenchymal-epithelial transition. Am J Pathol. 2002;160(6):2181–90.PubMed CentralPubMedView ArticleGoogle Scholar
- Ozaki H, Nakamura K, Funahashi J, Ikeda K, Yamada G, Tokano H, et al. Six1 controls patterning of the mouse otic vesicle. Development. 2004;131(3):551–62.PubMedView ArticleGoogle Scholar
- Li X, Oghi KA, Zhang J, Krones A, Bush KT, Glass CK, et al. Eya protein phosphatase activity regulates Six1-Dach-Eya transcriptional effects in mammalian organogenesis. Nature. 2003;426(6964):247–54.PubMedView ArticleGoogle Scholar
- Grifone R, Demignon J, Houbron C, Souil E, Niro C, Seller MJ, et al. Six1 and Six4 homeoproteins are required for Pax3 and Mrf expression during myogenesis in the mouse embryo. Development. 2005;132(9):2235–49.PubMedView ArticleGoogle Scholar
- Ikeda K, Ookawara S, Sato S, Ando Z, Kageyama R, Kawakami K. Six1 is essential for early neurogenesis in the development of olfactory epithelium. Dev Biol. 2007;311(1):53–68.PubMedView ArticleGoogle Scholar
- Micalizzi DS, Wang CA, Farabaugh SM, Schiemann WP, Ford HL. Homeoprotein Six1 increases TGF-beta type I receptor and converts TGF-beta signaling from suppressive to supportive for tumor growth. Cancer Res. 2010;70(24):10371–80.PubMed CentralPubMedView ArticleGoogle Scholar
- Micalizzi DS, Christensen KL, Jedlicka P, Coletta RD, Baron AE, Harrell JC, et al. The Six1 homeoprotein induces human mammary carcinoma cells to undergo epithelial-mesenchymal transition and metastasis in mice through increasing TGF-beta signaling. J Clin Invest. 2009;119(9):2678–90.PubMed CentralPubMedView ArticleGoogle Scholar
- Zheng XH, Liang PH, Guo JX, Zheng YR, Han J, Yu LL, et al. Expression and clinical implications of homeobox gene Six1 in cervical cancer cell lines and cervical epithelial tissues. Int J Gynecol Cancer. 2010;20(9):1587–92.PubMedGoogle Scholar
- Tan J, Zhang C, Qian J. Expression and significance of Six1 and Ezrin in cervical cancer tissue. Tumour Biol. 2011;32(6):1241–7.PubMedView ArticleGoogle Scholar
- Liu D, Li L, Zhang XX, Wan DY, Xi BX, Hu Z, et al. SIX1 Promotes Tumor Lymphangiogenesis by Coordinating TGFbeta Signals That Increase Expression of VEGF-C. Cancer Res. 2014;74(19):5597–607.PubMedView ArticleGoogle Scholar
- Wang CA, Jedlicka P, Patrick AN, Micalizzi DS, Lemmer KC, Deitsch E, et al. SIX1 induces lymphangiogenesis and metastasis via upregulation of VEGF-C in mouse models of breast cancer. J Clin Invest. 2012;122(5):1895–906.PubMed CentralPubMedView ArticleGoogle Scholar
- Jin H, Cui M, Kong J, Cui X, Lin Z, Wu Q, et al. Sineoculis homeobox homolog 1 protein is associated with breast cancer progression and survival outcome. Exp Mol Pathol. 2014;97(2):247–52.PubMedView ArticleGoogle Scholar
- Zlobec I, Minoo P, Baker K, Haegert D, Khetani K, Tornillo L, et al. Loss of APAF-1 expression is associated with tumour progression and adverse prognosis in colorectal cancer. Eur J Cancer. 2007;43(6):1101–7.PubMedView ArticleGoogle Scholar
- Zlobec I, Steele R, Terracciano L, Jass JR, Lugli A. Selecting immunohistochemical cut-off scores for novel biomarkers of progression and survival in colorectal cancer. J Clin Pathol. 2007;60(10):1112–6.PubMed CentralPubMedView ArticleGoogle Scholar
- Zlobec I, Vuong T, Hayashi S, Haegert D, Tornillo L, Terracciano L, et al. A simple and reproducible scoring system for EGFR in colorectal cancer: application to prognosis and prediction of response to preoperative brachytherapy. Br J Cancer. 2007;96(5):793–800.PubMed CentralPubMedView ArticleGoogle Scholar
- Zlobec I, Steele R, Michel RP, Compton CC, Lugli A, Jass JR. Scoring of p53, VEGF, Bcl-2 and APAF-1 immunohistochemistry and interobserver reliability in colorectal cancer. Mod Pathol. 2006;19(9):1236–42.PubMedView ArticleGoogle Scholar
- Cai MY, Zhang B, He WP, Yang GF, Rao HL, Rao ZY, et al. Decreased expression of PinX1 protein is correlated with tumor development and is a new independent poor prognostic factor in ovarian carcinoma. Cancer Sci. 2010;101(6):1543–9.PubMedView ArticleGoogle Scholar
- Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global cancer statistics. CA Cancer J Clin. 2011;61(2):69–90.PubMedView ArticleGoogle Scholar
- Jemal A, Murray T, Samuels A, Ghafoor A, Ward E, Thun MJ. Cancer statistics, 2003. CA Cancer J Clin. 2003;53(1):5–26.PubMedView ArticleGoogle Scholar
- Glinsky GV, Glinskii AB, Stephenson AJ, Hoffman RM, Gerald WL. Gene expression profiling predicts clinical outcome of prostate cancer. J Clin Invest. 2004;113(6):913–23.PubMed CentralPubMedView ArticleGoogle Scholar
- Markert EK, Mizuno H, Vazquez A, Levine AJ. Molecular classification of prostate cancer using curated expression signatures. Proc Natl Acad Sci U S A. 2011;108(52):21276–81.PubMed CentralPubMedView ArticleGoogle Scholar
- Ding Z, Wu CJ, Chu GC, Xiao Y, Ho D, Zhang J, et al. SMAD4-dependent barrier constrains prostate cancer growth and metastatic progression. Nature. 2011;470(7333):269–73.PubMed CentralPubMedView ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.